!!------ Skeleton for body section -----
!Body 1 = air
!Body 2 = trace 1
!Body 3 = trace 2
!------ Skeleton for boundary section -----
!BC1 = Current in coil 1
!BC2 = Current out coil 1
!BC3 = Current in coil 2
!BC4 = Current out coil 2
!BC5 = Air_BC (external air, where vector potential is assumed to be 0)
! All coils excited
Header
Check Keywords "Warn"
Mesh DB "." "adjacent_open_coil_V01"
Include Path ""
Results Directory ""
End
Simulation
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
! Coordinate Scaling=0.001
Simulation Type = Steady
Steady State Max Iterations = 1
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 = "Air"
Equation = 1
Material = 1
End
Body 2
Target Bodies(1) = 2
Name = "Coil 1"
Equation = 2
Material = 2
Body Force = 1
End
Body 3
Target Bodies(1) = 3
Name = "Coil 2"
Equation = 2
Material = 2
Body Force = 1
End
Material 1
Name = "Air (room temperature)"
Viscosity = 1.983e-5
Heat expansion Coefficient = 3.43e-3
Heat Conductivity = 0.0257
Sound speed = 343.0
Density = 1.205
Heat Capacity = 1005.0
Relative Permittivity = Real 1
Relative Permeability = Real 1
Electric Conductivity = Real 0
End
Material 2
Name = "Copper (generic)"
Heat expansion Coefficient = 16.5e-6
Heat Conductivity = 401.0
Sound speed = 3810.0
Heat Capacity = 385.0
Density = 8960.0
Mesh Poisson ratio = 0.34
Poisson ratio = 0.34
Youngs modulus = 115.0e9
Relative Permittivity = Real 1
Relative Permeability = Real 1
Electric Conductivity = Real 5.96e7
End
Equation 1
Name = "Equation for air"
Active Solvers(2) = 2 3
End
Equation 2
Name = "Equation for coils 1 and 2"
Active Solvers(3) = 1 2 3
End
Solver 1
Equation = "solving currents"
Calculate Volume Current = True
Calculate Elemental Fields = True
Procedure = "StatCurrentSolveVec" "StatCurrentSolver"
Variable = Potential
Nonlinear System Max Iterations = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 1000
Linear System Convergence Tolerance = 1.0e-10
Linear System Preconditioning = ILU1
Linear System Abort Not Converged = False
Linear System Residual Output = 10
End
!===============================================
!=========== MAGNETIC FIELD SOLVER =============
Solver 2
! Exec Solver = Never
Equation = "MGDynamics"
Variable = P
Procedure = "MagnetoDynamics" "WhitneyAVSolver"
Fix Input Current Density = Logical True
Nonlinear System Max Iterations = 1
Linear System Solver = "Iterative"
Linear System Iterative Method = BiCGStabl
BicGStabL Polynomial Degree = 4
Linear System Max Iterations = 1000
Linear System Convergence Tolerance = 1.0e-8
Linear System Preconditioning = none
Linear System Abort Not Converged = False
Linear System Residual Output = 10
End
Solver 3
! Exec Solver = Never
Equation = "MGDynamicsCalc"
Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields"
Potential Variable = P
! Calculate Current Density = Logical True
! Calculate Electric Field = Logical True
! Calculate Magnetic Field Strength = Logical True
Calculate Magnetic Field Density = Logical True
Steady State Convergence Tolerance = 1e-09
Linear System Solver = "Iterative"
Linear System Iterative Method = CG
Linear System Preconditioning = ILU0
Linear System Max Iterations = 1000
Linear System Convergence Tolerance = 1.0e-8
End
Solver 4
! Equation = Result Output
! Exec Solver = Never
Procedure = "ResultOutputSolve" "ResultOutputSolver"
Binary Output = True
Output Format = vtu
Output File Name = results
Save Geometry Ids = Logical True
End
Solver 5
! Equation = "scalars"
! Exec Solver = Never
procedure = "SaveData" "SaveScalars"
Filename = "data_scalars.dat"
End
Boundary Condition 1
Target Boundaries(1) = 5
Name = "Outer space"
P {e}= real 0
End
Boundary Condition 2
Target Boundaries(1) = 1
Name = "Coil 1 current in"
Current Density = 40.0
Current Density Normalize by Area = Logical True
End
Boundary Condition 3
Target Boundaries(1) = 2
Name = "Coil 1 current out"
Potential = real 0
End
Boundary Condition 4
Target Boundaries(1) = 3
Name = "Coil 2 current in"
Current Density = 10.0
Current Density Normalize by Area = Logical True
End
Boundary Condition 5
Target Boundaries(1) = 4
Potential = real 0
End
Body Force 1
Name = "Coils Current Density"
! Elemental current works better
Current Density 1 = Equals "Volume current e 1"
Current Density 2 = Equals "Volume current e 2"
Current Density 3 = Equals "Volume current e 3"
End