solvers_steadyGroundwater2DFoam

steadyGroundwater2DFoam solver

Description

Steady solver for free-surface flow in porous media (steady version of groundwater2DFoam)

Solver uses the OpenFOAM "SIMPLE" structure

Transient solver for free-surface flow though porous media (2D modeling)

Solver requires several geographical information (see below).

A seepage option can be activated with Digital Elevation Model to constrain maximal water height (if potential > potentialDEM) and cell flux < 0 : solver sets potential=potentialDEM at this cell and computes exfiltration outflow)

The file constant/potentialDEM gives the maximal potential values used for seepage. This file can be generated automatically by bilinear interpolation using a DEM file at the XYZ format ASCII.

Explicit fixed potentials (locations and values) can be given as a list in constant/transportProperties. The user-defined values in list can be overriden by bilinear interpolation using the field constant/potentialDEM (requires activation of useDEMtoFixPotential option, see below).

Configuration files

constant/transportProperties

// fluid properties
rho rho [1 -3 0 0 0 0 0] 1000;
mu mu [1 -1 -1 0 0 0 0] 1e-03;

infiltration infiltration [0 1 -1 0 0 0 0] -5e-9; // constant uniform infiltration velocity

seepage yes; // Seepage activation (requires constant/potentialDEM field)

fileDEM dem.xyz; //(optional) mnt file read by solver to construct constant/potentialDEM if not present

fixedPotentialList // list of the fixed potential positions (and user-defined potential values)
(
    ( (584256.50 2413589.00 1) 145.2)
    ( (584556.50 2413889.00 1) 149.3)
);

useDEMtoFixPotential yes; // (use potentialDEM values to impose potential (instead of user-defined potential values)

system/fvSolution

solvers
{
    potential
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-12;
        relTol          0;
    }

}

SIMPLE
{
    nNonOrthogonalCorrectors 0;

    residualControl
    {
        potential          1e-10; //residual control for solving potential
    }
}

relaxationFactors
{
    fields
    {
        potential          0.5; // relaxation factor for updating solution
    }
}

Required fields

• 0/potential : The potential field
• 0/U : The velocity field
• constant/K : permeability field
• constant/z0 : height of the impermeable layer

Optional fields

• 0/infiltration : infiltration velocity field (replace infiltration in transport properties)
• constant/potentialDEM : field used by forcing potential modes: fixedPotential and seepage

Output computed fields

• 0/hwater : water height in the watershed

Return to porousMultiphaseFoam main page