WhitneyAV 3D coupled current emag of a coil

Elmer cases by the users for the users
jtapas3761
Posts: 8
Joined: 16 Jul 2018, 10:17
Antispam: Yes

Re: WhitneyAV 3D coupled current emag of a coil

Post by jtapas3761 »

If this will help to find the error !!!!!

**********Sif file
Header
CHECK KEYWORDS "Warn"
Mesh DB "." "exportmesh"
End

Simulation
Max Output Level = 9
Coordinate System = "Cartesian 3D"
Coordinate Mapping(3) = 1 2 3
Simulation Type = Steady state
Steady State Max Iterations = 1
Output Intervals = 1
Post File= "xyz.ep"
End

Constants
Permittivity of Vacuum = 8.8542e-12
End

Solver 1
Exec Solver = Before Simulation

Procedure = "StatCurrentSolve" "StatCurrentSolver"
Equation = "Stat Current Solver"
Variable = Potential
Variable DOFs = 1
Calculate Volume Current = True
Calculate Joule Heating = False
Current Control = Real 12.5 ! 1/4th of 50 A
Linear System Solver = Iterative
Linear System Iterative Method = CG
Linear System Max Iterations = 1000
Linear System Convergence Tolerance = 1.0e-6
Linear System Preconditioning = ILU0
Linear System Abort Not Converged = True
Linear System Residual Output = 1
End

Solver 2
Equation = "MGDynamics"

Variable = "P"

Procedure = "MagnetoDynamics" "WhitneyAVSolver"
Fix Input Current Density = Logical True

Newton-Raphson Iteration = Logical True
Nonlinear System Max Iterations = 30
Nonlinear System Convergence Tolerance = 1e-6

Linear System Symmetric = Logical True
Linear System Solver = "Iterative"
Linear System Preconditioning = None
Linear System Convergence Tolerance = 1e-8
Linear System Residual Output = 100
Linear System Max Iterations = 5000
Linear System Iterative Method = CG
Steady State Convergence Tolerance = 1e-6
!Automated Source Projection BCs= Logical False

End

Solver 3
Equation = "MGDynamicsCalc"
Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields"
Linear System Symmetric = True
Potential Variable = String "P"

Calculate Magnetic Field Strength = Logical True
Steady State Convergence Tolerance = 0
Linear System Solver = "Iterative"
Linear System Preconditioning = None
Linear System Residual Output = 0
Linear System Max Iterations = 5000
Linear System Iterative Method = CG
Steady State Convergence Tolerance = 1e-6
Linear System Convergence Tolerance = 1.0e-8
End

Solver 4
Exec Solver = String "after all"
exec interval = 1
Equation = String "ResultOutput"
Procedure = File "ResultOutputSolve" "ResultOutputSolver"
Output File Name = File "emag."
Output Format = String "vtk"
Vtk Format = Logical True
End


Equation 1
Name = "Coupled Equations"
Active Solvers(2) = 2 3
End

Equation 2
Name = "Coupled Equations"
Active Solvers(3) = 1 2 3
End

Material 1
Name = "Air"
Electric Conductivity = 0.0
Relative Permeability = 1.0
End

Material 2
Name = "Copper"
Electric Conductivity = 58.14e6
Relative Permeability = 1.0
End


Body 1
Name = "Air"
Target Bodies(1) = 1
Equation = 1
Material = 1
End

Body 2
Name = "Coil"
Target Bodies(1) = 2
Equation = 2
Material = 2
Body Force = 1
End


Boundary Condition 1
Name = "BC Potential"
Target Boundaries(1) =7
P {e}= real 0
P = Real 0
Potential = 0.001
End

Boundary Condition 2
Name = "BC Voltage Ground"
Target Boundaries(1) = 8
P {e} = Real 0
P = Real 0
Potential = 0.0
End

Boundary Condition 3
Name = "BCn Flux Parallel"
Target Boundaries(6) = 1 2 3 4 5 6 7 8
P {e} = Real 0
P = Real 0.0
End

Body Force 1
Name = "Current Density"
Current Density 1 = Equals Volume current 1
Current Density 2 = Equals Volume current 2
Current Density 3 = Equals Volume current 3
End
Phlix
Posts: 11
Joined: 31 Aug 2018, 13:51
Antispam: Yes

Re: WhitneyAV 3D coupled current emag of a coil

Post by Phlix »

Hi,

I am have the same problem.

Is there a working solution so far?

What are other possibilities to get convergence in the current solver?

I appreciate every hint to solve the problem!

Best regards

Phlix
lazyTom
Posts: 56
Joined: 05 Jun 2013, 01:32
Antispam: Yes

Re: WhitneyAV 3D coupled current emag of a coil

Post by lazyTom »

Hi all,

out of curiosity, did anybody check how the mesh looks like before running the simulation?

I mean, have a look at the pictures below, it doesn't look good... I suspect this is the reason why the solver is not converging.
brokenMesh1.png
brokenMesh1.png (103.52 KiB) Viewed 7202 times
brokenMesh2.png
brokenMesh2.png (201.84 KiB) Viewed 7202 times
brokenMesh3.png
brokenMesh3.png (169.48 KiB) Viewed 7202 times
If you are still interested in this model I have prepared something similar here: viewtopic.php?f=3&t=6427&p=20338#p20262

There you will find a lo of information about current solvers and how do they work coupled with the magnetic solver.

Regards,

Marco
Last edited by lazyTom on 10 Jun 2019, 13:16, edited 1 time in total.
jtapas3761
Posts: 8
Joined: 16 Jul 2018, 10:17
Antispam: Yes

Re: WhitneyAV 3D coupled current emag of a coil

Post by jtapas3761 »

Hello Everyone,

I am using this module to generate magnetic field in 3D around a current carrying coil. In this the first solver being used is a STAT CURRENT SOLVER, where we apply some Voltage as boundary conditions. Is it possible to set a constant current in the coil instead of applying a voltage BC? Is there a different solver for this or do we need to make some modifications in the STAT CURRENT Solver? Could anyone provide any guidelines or send a link to a sample project for reference? Thanks!

Tapas
mzenker
Posts: 1999
Joined: 07 Dec 2009, 11:49
Location: Germany

Re: WhitneyAV 3D coupled current emag of a coil

Post by mzenker »

Hi,

you can set a constant current as boundary condition, see models manual.

HTH,
Matthias
Post Reply