Re: WhitneyAV 3D coupled current emag of a coil
Posted: 13 Sep 2018, 12:00
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
**********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