Model induction between two coils

[mrbearb]

GoodViewInduction.png (884.78 KiB) Viewed 223 times

Induction? About to add magnetic cores and re-run

[mrbearb]

Ok, here's the updated run with iron cores:

MagFieldsAndCurrentDensity.png (656.32 KiB) Viewed 212 times

Any idea why I am unable to see Electric Fields? I have the option enabled in the sif, but the electric fields do not show up in Paraview

My sif:

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) = 1
  Solver Input File = case.sif
  Post File = case.vtu
  Calculate Mesh Pieces = True
  Desired Mesh Pieces = 1
End

Constants
  Gravity(4) = 0 -1 0 9.82
  Stefan Boltzmann = 5.670374419e-08
  Permittivity of Vacuum = 8.85418781e-12
  Permeability of Vacuum = 1.25663706e-6
  Boltzmann Constant = 1.380649e-23
  Unit Charge = 1.6021766e-19
End

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

Body 2
  Target Bodies(1) = 2
  Name = "CoilBody2"
  Equation = 2
  Material = 1
  Initial condition = 1
End

Body 3
  Target Bodies(1) = 3
  Name = "Body Property 3"
  Equation = 2
  Material = 2
End

Solver 1
  Equation = CoilSolver
  Fix Input Current Density = True
  Procedure = "CoilSolver" "CoilSolver"
  Coil Center(3) = 0 0 0
  Coil Normal(3) = 1 0 0
  Desired Coil Current = 1
  Calculate Nodal Fields = True
  Normalize Coil Current = True
  Coil Closed = True
  Narrow Interface = True
  Desired Current Density = 100000
  Exec Solver = Always
  Stabilize = True
  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 = 2000
  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 2
  Equation = MgDyn
  Use Elemental CoilCurrent = True
  Fix Input Current Density = True
  Procedure = "MagnetoDynamics" "WhitneyAVSolver"
  Exec Solver = Always
  Stabilize = True
  Optimize Bandwidth = True
  Steady State Convergence Tolerance = 1.0e-6
  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 = 10000
  Linear System Convergence Tolerance = 1.0e-10
  BiCGstabl polynomial degree = 6
  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 3
  Equation = MgDynPost
  Calculate Current Density = True
  Calculate Electric Field = True
  Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields"
  Calculate Magnetic Field Strength = True
  Skip Nodal Fields = False
  Exec Solver = Always
  Stabilize = True
  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 = "ActiveEquation"
  Active Solvers(3) = 1 2 3
End

Equation 2
  Name = "PassiveEquation"
  Active Solvers(2) = 3 2
End

Material 1
  Name = "Copper (generic)"
  Density = 8960.0
  Electric Conductivity = 59.59e6
  Heat Conductivity = 401.0
  Sound speed = 3810.0
  Youngs modulus = 115.0e9
  Heat expansion Coefficient = 16.5e-6
  Heat Capacity = 385.0
  Poisson ratio = 0.34
  Relative Permeability = 0.999994
End

Material 2
  Name = "Air (room temperature)"
  Heat Conductivity = 0.0257
  Heat Capacity = 1005.0
  Viscosity = 1.983e-5
  Heat expansion Coefficient = 3.43e-3
  Relative Permittivity = 1.00059
  Relative Permeability = 1.00000037
  Density = 1.205
  Sound speed = 343.0
End

Material 3
  Name = "Iron (generic)"
  Poisson ratio = 0.29
  Youngs modulus = 193.053e9
  Heat Conductivity = 80.2
  Density = 7870.0
  Electric Conductivity = 10.30e6
  Sound speed = 5000.0
  Relative Permeability = 6000
  Heat expansion Coefficient = 11.8e-6
  Heat Capacity = 449.0
End

Body Force 1
  Name = "Source"
  Current Density 2 = Equals "CoilCurrent e 2"
  Current Density 1 = Equals "CoilCurrent e 1"
  Current Density 3 = Equals "CoilCurrent e 3"
End

Initial Condition 1
  Name = "InitialCondition 2"
  AV = 0
  AV {e} = 0
End

Boundary Condition 1
  Target Boundaries(6) = 13 14 15 16 17 18 
  Name = "BoxCondition"
  AV = 0
  AV {e} = 0
End