porous_sat#

EQ = porous_sat {Galerkin_wt} P_liq {Interpol_fnc} <floatlist>

Description / Usage#

This card provides information for solving a differential equation for saturated porous flow. Definitions of the input parameters are defined below. Note that <floatlist> has five parameters to define the constant multipliers in front of each type of term in the equation.The Galerkin weight and the interpolation function must be the same for the code to work properly. If upwinding is desired for advection dominated problems, we can set this through a Petrov-Galerkin weight function in the material file.

porous_sat

Name of the equation to be solved.

{Galerkin_wt}

Two-character value that defines the type of weighting function for this equation, where:

  • Q1-Linear

  • Q2-Quadratic

P_SAT

Name of the variable associated with this equation.

{Interpol_fnc}

Two-character value that defines the interpolation function used to represent the variable P_SAT, where:

  • Q1-Linear Continuous

  • Q2-Quadratic Continuous

<float1>

Multiplier on mass matrix term ( d ⁄dt ).

<float2>

Multiplier on advective term.

<float3>

Multiplier on boundary term ( \(\underline{n}\) • flux ).

<float4>

Multiplier on diffusion term.

<float5>

Multiplier on source term.

Note: These multipliers are intended to provide a means of activating or deactivating terms of an equation, and hence should be set to zero or one. If a multiplier is zero, the section of code that evaluates the corresponding term will be skipped.

Examples#

The following is a sample card that uses a linear continuous interpolation and weight function for the saturated porous equation and has all the term multipliers on except for the mass matrix for time derivatives:

EQ = porous_sat Q1 P_liq Q1 0. 1. 1. 1. 1.

Technical Discussion#

This card is not yet in use.

References#

No References.