LorentzForce

Syntax

Defined in couplingProperties dictionary.

meshMotionModel LorentzForce;
LorentzForceProps
{
    sigma        scalar;
}

sigma : electric conductivity of the fluid

Examples

meshMotionModel LorentzForce;
LorentzForceProps
{
    sigma         1;
}

Description

The LorentzForceModel realizes an electric potential method solver allowing for accounting for the effect of a static magnetic field on flow (i.e. 1-way MHD coupling). An exlicit Lorentz force term, calculated from a Poisson equation for the electric potential, is accounted for in the momentum balance equation, following Ferroni.

The user must provide the follwing fields (initial solutian and boundary conditions):

B: the magnetic field (more specifically, the flux density) (units: kg / (A s^2))
PotE : the calculated electric potential (units: J / (A s))
Lorentz: the calculated volumetric Lorentz force (units: N / m^3)

The user must, moreover, provide the following solver settings in fvSolution:

Settings on how to solve the equation for the electirf field PotE, e.g.:
```
"(PotE|PotEFinal)"
{
    $p;
};
```

PIMPLE-like solver settings:

PotE
{
    nCorrectors     0;
    nNonOrthogonalCorrectors 0;
    PotERefCell        0;
    PotERefValue       0;
}

Restrictions

Requires usage of cfdemSolverPimple, cfdemSolverRhoPimple, or cfdemSolverChem.

Literature

[1] Ferroni, F. (2012). Magneto-hydrodynamic simulations of liquid metal flows in fusion reactors (Doctoral dissertation, Master Thesis, Imperial College London).