pyoomph.equations.darcy module

class pyoomph.equations.darcy.DarcyEquation(fluid_props, permeability=1.0000000000000000777E-15 * meter**2, porosity=0.3, solve_also_velocity=False)

Bases: Equations

Represents the Darcy equation for porous media flow. The governing equations are given by:

d(rho_f*eta)/dt + div(rho_f*eta*u) = 0 eta*u = -kappa/mu * grad(p)

where rho_f is the fluid density, eta is the porosity, u is the velocity, kappa is the permeability, mu is the dynamic viscosity and p is the porous pressure. Optionally, the velocity field can be not solved for, in which case it is computed from the porous pressure, such that the equation for pressure is:

div(rho_f*kappa/mu*grad(p)) = d(rho_f*eta)/dt

Parameters:

• fluid_props (Union[PureLiquidProperties, PureGasProperties, MixtureLiquidProperties, MixtureGasProperties]) – The fluid properties of the fluid in the porous media.

• permeability (Union[Expression, int, float]) – The permeability of the porous media. Default is 1e-15 m^2.

• porosity (Union[Expression, int, float]) – The porosity of the porous media. Default is 0.3.

• solve_also_velocity (bool) – If True, the velocity field is also solved for. Default is False.

class pyoomph.equations.darcy.PorousFront(*args: Any, **kwargs: Any)

Bases: InterfaceEquations

Move the mesh along with the porous velocity.

This class requires the parent equations to be of type DarcyEquation, meaning that if DarcyEquation (or subclasses) are not defined in the parent domain, an error will be raised.

required_parent_type

alias of DarcyEquation

class pyoomph.equations.darcy.PorousNavierStokesConnection(*args: Any, **kwargs: Any)

Bases: InterfaceEquations

Imposes the droplet pressure to the porous domain and enforces the z-velocity of the droplet to be the porous velocity.

This class requires the parent equations to be of type DarcyEquation, meaning that if DarcyEquation (or subclasses) are not defined in the parent domain, an error will be raised.

required_parent_type

alias of DarcyEquation