Simple Flow Problem ExampleΒΆ

This input file will run a vertical 1D column with variable saturated flow using RICHARDS Mode.


#Description: 1D infiltration

SIMULATION
SIMULATION_TYPE SUBSURFACE
PROCESS_MODELS
SUBSURFACE_FLOW flow
MODE RICHARDS
/
/
END

SUBSURFACE

#=========================== discretization ===================================
TYPE structured
ORIGIN 0.d0 0.d0 0.d0
NXYZ 1 1 10
DXYZ
1.d0
1.d0
1.d0
/
END

#=========================== material properties ==============================
ID 1
POROSITY 0.33
TORTUOSITY 0.5d0
PERMEABILITY
PERM_X 1.d-10 ! ~100 m/d
PERM_Y 1.d-10
PERM_Z 1.d-10
/
/

#=========================== saturation functions =============================
SATURATION_FUNCTION VAN_GENUCHTEN
LIQUID_RESIDUAL_SATURATION 0.1d0
M 0.8d0
ALPHA 1.d-4
/
PERMEABILITY_FUNCTION MUALEM_VG_LIQ
LIQUID_RESIDUAL_SATURATION 0.1d0
M 0.8d0
/
/

#=========================== output options ===================================
FORMAT TECPLOT POINT
/

#=========================== times ============================================
FINAL_TIME 10.d0 y
MAXIMUM_TIMESTEP_SIZE 0.1d0 y
/

#=========================== regions ==========================================
REGION all
COORDINATES
0.d0 0.d0 0.d0
1.d0 1.d0 10.d0
/
END

REGION top
COORDINATES
0.d0 0.d0 10.d0
1.d0 1.d0 10.d0
/
FACE TOP
END

REGION bottom
COORDINATES
0.d0 0.d0 0.d0
1.d0 1.d0 0.d0
/
FACE BOTTOM
END

#=========================== flow conditions ==================================
TYPE
FLUX neumann
/
FLUX 3.171d-10 ! 1 cm/yr
END

TYPE
PRESSURE hydrostatic
/
PRESSURE 101325.d0
END

#=========================== condition couplers ===============================
# initial condition
REGION all
END

# top boundary condition
END

# bottom boundary condition
REGION bottom
END

#=========================== stratigraphy couplers ============================
REGION all
MATERIAL soil
END

END_SUBSURFACE