2D magnet simulation - unphysical results
Posted: 23 Sep 2019, 19:10
I have been trying to simulate a simple 2D magnet and I referred horseshoe magnet tutorial for this, however while the horseshoe simulation shows sensible result, I can't get 2D magnet to work.
The potential seems to be "caged" inside the top and bottom magnetization bodies and doesn't spread through inner iron body or air.
sif file is as below:
What am I missing here?
Mesh is shown below
The potential seems to be "caged" inside the top and bottom magnetization bodies and doesn't spread through inner iron body or air.
sif file is as below:
Code: Select all
Header
CHECK KEYWORDS Warn
Mesh DB "." "."
Include Path ""
Results Directory ""
End
Simulation
Max Output Level = 5
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
Simulation Type = Steady state
Steady State Max Iterations = 1
Output Intervals = 1
Timestepping Method = BDF
BDF Order = 1
Solver Input File = case.sif
Post File = case.ep
End
Constants
Gravity(4) = 0 -1 0 9.82
Stefan Boltzmann = 5.67e-08
Permittivity of Vacuum = 8.8542e-12
Boltzmann Constant = 1.3807e-23
Unit Charge = 1.602e-19
End
Body 1
Target Bodies(1) = 1
Name = "Body Property 1"
Equation = 1
Material = 1
End
Body 2
Target Bodies(1) = 2
Name = "Body Property 2"
Equation = 1
Material = 2
End
Body 3
Target Bodies(1) = 3
Name = "Body Property 3"
Equation = 1
Material = 3
End
Body 4
Target Bodies(1) = 4
Name = "Body Property 4"
Equation = 1
Material = 4
End
Solver 2
Equation = MgDyn2DPost
Procedure = "MagnetoDynamics2D" "BSolver"
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 20
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-10
BiCGstabl polynomial degree = 2
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 10
Linear System Precondition Recompute = 1
End
Solver 1
Equation = MgDyn2D
Procedure = "MagnetoDynamics2D" "MagnetoDynamics2D"
Variable = Potential
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 20
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-10
BiCGstabl polynomial degree = 2
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 10
Linear System Precondition Recompute = 1
End
Equation 1
Name = "MgDyn2D"
Active Solvers(2) = 2 1
End
Equation 2
Name = "MgDyn2DPost"
Active Solvers(2) = 2 1
End
Material 1
Name = "Air"
Relative Permeability = 1
End
Material 2
Name = "Iron"
Relative Permeability = 5000
End
Material 3
Name = "IronPlus"
Relative Permeability = 5000
Magnetization 2 = Real 8000
End
Material 4
Name = "IronMinus"
Relative Permeability = 5000
Magnetization 2 = Real -8000
End
Boundary Condition 1
Target Boundaries(4) = 1 2 3 4
Name = "FarField"
Infinity BC = True
End
Mesh is shown below