Boundary conditions for 3D WhitneyAVHarmonicSolver

Numerical methods and mathematical models of Elmer
Post Reply
vencels
Posts: 24
Joined: 20 Sep 2016, 17:05
Antispam: Yes
Location: Latvia
Contact:

Boundary conditions for 3D WhitneyAVHarmonicSolver

Post by vencels » 13 Jun 2017, 10:47

Hi,

I am trying to solve inductor in the air, see picture attached.
Inductor is made from highly conductive material and frequencies are about 1kHz so skin effect should be taken in account. Lets say I have 1kA current flowing in the inductor.

I see two options to model this:
1) Set appriximate current distribution as a analytic function.
2) Set guessed potential to inductor's terminals. Then in postprocessing integrate and compute resulting current, and scale potential to match that specified 1kA current.

Q1: Is there a better way to do this?
Q2: What boundary condition is the most appropriate on YZ plane where terminals are located and the outer airbox?

I triedBC's from mgdyn_harmonic_wire (see code below), but I lack physical interpretation whether these BC's are correct.

Thanks!

Code: Select all

Boundary Condition 1
  Name = "terminal1"
  Target Boundaries(1) = 1 

  Coil Start = Logical True

  AV re {e} 2 = Real 0.0
  AV re {e} 3 = Real 0.0
  AV re = Real 0.0

  AV im {e} 3 = Real 0.0
  AV im {e} 4 = Real 0.0
  AV im = Real 0.0
End

Boundary Condition 2
  Name = "terminal2"
  Target Boundaries(1) = 2

  Coil End = Logical True

  AV re {e} 2 = Real 0.0
  AV re {e} 3 = Real 0.0
  AV re = Real 1

  AV im {e} 2 = Real 0.0
  AV im {e} 3 = Real 0.0
  AV im = Real 0.0
End

Boundary Condition 3
  Name = "terminalPlane"
  Target Boundaries(1) = 3

  AV re {e} 2 = Real 0.0
  AV re {e} 3 = Real 0.0

  AV im {e} 2 = Real 0.0
  AV im {e} 3 = Real 0.0

Boundary Condition 4
  Name = "remainingAirBox"
  Target Boundaries(1) = 4

  AV re {e} = Real 0.0
  AV im {e} = Real 0.0
End
Attachments
inductorInAir.png
inductorInAir_picture
(26.55 KiB) Not downloaded yet

Post Reply