VxB fof WhitneyAVHarmonic

Discussion about coding and new developments

VxB fof WhitneyAVHarmonic

Postby vencels » 06 Oct 2017, 12:11


I am using WhitneyAVSolver as a reference and trying to implement VxB term into WhitneyAVHarmonicSolver and CalcFields.
Next step is sending conductive fluid velocity from OpenFOAM to Elmer. I think that from Elmer side I need to feed this velovity to the variable "lorentz_velo(3,27)".

1) Isn't it coded in such way that lorentz_velo is homogenious? In WhitneyAVSolver I cannot see any spatial dependence for lorentz_velo.
2) If lorentz_velo is implemented in such way that it is inhomogenious, will the solver converge?

Posts: 22
Joined: 20 Sep 2016, 17:05
Location: Latvia

Re: VxB fof WhitneyAVHarmonic

Postby kataja » 10 Oct 2017, 09:11


Yeah, the "lorentz_velo" refers to nodal interpolant to nodes of the current element of the real valued vector variable "Lorentz velocity" appearing in the sif file.

So in the sif file you could have
Code: Select all
Lorentz velocity 1 = variable coordinate 2
real matc "-tx"
Lorentz velocity 2 = variable coordinate 1
real matc "tx"

and this will correspond to field rotating around z-axis.

For rotating fields there's also the "Angular velocity" where you give w in $v = w \times r$ (where r is the position vector).

So the Lorentz velocity can be pretty much anything because it is treated as a normal elmer variable.

Posts: 64
Joined: 09 May 2014, 16:06

Return to Software development

Who is online

Users browsing this forum: No registered users and 1 guest