ShirgaonkarIB

Syntax

Defined in couplingProperties dictionary.

forceModels
(
    ShirgaonkarIB
);

ShirgaonkarIBProps
{
    velFieldName         "U";
    pressureFieldName    "p";
    voidfractionFieldName "voidfraction";
    twoDimensional;
    depth                scalar;
    useTotalViscosity    Switch;

    // (optional) forceSubModel switches
    verbose              switch;
    treatForceExplicit   switch;
    implForceDEM         switch;
};

  • U = name of the finite volume fluid velocity field (default: “U”)

  • p = name of the finite volume pressure field (default: “p”)

  • voidfraction = name of the finite volume voidfraction field (default: “voidfraction”)

  • twoDimensional = optional keyword for conducting a two dimensional calculation. Keyword only.

  • depth = depth of domain for two dimensional cases. Only read if twoDimensional is active.

  • useTotalViscosity = use the total viscosity, e.g. when using a non-Newtonian viscosity (default: false)

This forceModel reads the following forceSubModel switches and overwrites the defaults as indicated in parentheses:

  • scaleDrag

  • treatForceExplicit

  • implForceDEM

  • verbose

  • useTorque

Examples

forceModels
(
    ShirgaonkarIB
);

ShirgaonkarIBProps
{
    velFieldName "U";
    pressureFieldName "p";
    voidfractionFieldName "voidfraction";
}

Description

Calculates the drag force acting on each particle in a resolved manner by summing viscous and drag forces over all cells that the particle covers, see [1].

Note

This model is only suited for resolved CFD-DEM simulations where the particle is represented by the immersed boundary method. If the hybrid voidfracion model is used the voidfractionFieldName here should match the one set for the IB voidfracion model.

Literature

[1] A. A. Shirgaonkar, M. A. MacIver and N. A. Patankar. “A new mathematical formulation and fast algorithm for fully resolved simulation of self-propulsion”. Journal of Compuation Physics 228 (2009): 2366-2390

Restrictions

Only for immersed boundary solvers like cfdemSolverIB or cfdemSolverPimple.