Magnetic Induction and Electric Potential Solvers For Incompressible MD Flows
Magnetic Induction and Electric Potential Solvers For Incompressible MD Flows
Magnetic Induction and Electric Potential Solvers For Incompressible MD Flows
Alessandro Tassone
05/12/2016
Fusion reactors
u=0 (1)
Du
= (p/) + 2 u + (J B)/ (2)
Dt
B 1 2
= (u B) + B (3)
t
1
J= B (4)
This set is called the B-formulation of the MHD governing equations
Inductionless approximation
Rm 1
A Poisson equation for the electric potential and the Ohms law substitute
(3) and (4)
u=0 (1)
Du
= (p/) + 2 u + (J B)/ (2)
Dt
2 = (u B) (5)
J = ( + u B) (6)
The new set is called the -formulation of the MHD equations
(B0 + b)t+1 = 0
Application
PISO(u, p) BPISO(B, pB)
Figure: Solenoid magnetic field
cd $FOAM_SOLVERS/electromagnetics/mhdFoam
ls
gedit mhdFoam.C
if (bpiso.finalNonOrthogonalIter())
{
phiB -= pBEqn.flux();
}
mag(fvc::div(phiB))().weightedAverage(mesh.V()).value()
Case parameters
M = 20
Re = 2
c=0
B|wall = B0
Available in
$FOAM_TUTORIALS
/electromagnetics
/mhdFoam/hartmann
Figure: Hartmann problem
Drawbacks
Requirement on J interpolation
Nonconservative treatment of Lorentz force
Worse convergence behavior
Constrain on charge conservation J = 0
Source code
jfinal.correctBoundaryConditions();
Compile epotFoam
mkdir -p $WM_PROJECT_USER_DIR/applications/solvers/MHD
mv epotFoam $WM_PROJECT_USER_DIR/applications/solvers/MHD
cd $WM_PROJECT_USER_DIR/applications/solvers/MHD/epotFoam
wmake
Case parameters
M = 20
Re = 2
cs = 0
ch =
blockMesh
epotFoam >& log &
blockMeshDict
Boundary Conditions
Insulating wall
=0 type zeroGradient;
n
For other patches the setting is zeroGradient, except for the symmetry
plane. Velocity and pressure fields employ usual BCs (i.e. noSlip,
inlet-velocity/outlet-pressure)
constant
Magnetic induction intensity is defined in electromagneticProperties
All transportProperties are assumed to be constant and equal to 1
system
t = 2.5 e-3 s
Temporal discretization scheme: Crank-Nicholson (explicit, second order)
Spacial discretization scheme: Gauss liner (central differencing)
Residual tolerance: 1e-06 (p), 1e-05 (u, PotE)
Post Process
The sampling data utility produces data for every written time step. The
plot.sh is a script for the generation of the comparison picture employing
gnuplot
Results
Results
Figure: Velocity contour Hunt flow Figure: Potential contour Hunt flow
Follow-up activities
Tutorials for elemental 2D MHD flows (Hunt and Shercliff cases) were
developed for both mhdFoam and epotFoam