Skip to Main Content U.S. Department of Energy


Short Course Example Problem 6

Download problem description | Download input file (txt)

NAPL infiltration and redistribution in a 2D aquifer

Abstract:The first objective of this two-part example is to investigate the effects of fluid density and viscosity on the movement of NAPLs after a spill in a partly saturated, hypothetical, aquifer. The user will investigate the effects of NAPL fluid properties on infiltration and redistribution. The objective of the second part is to compute steady-state initial conditions through a separate simulation. The Restart file generated by this simulation is used to define the initial conditions for the actual infiltration event. The Water-Oil mode (STOMP-WO) is used for this application.

Problem Description

Part 1:
A known amount of LNAPL, with a density of 800 kg/m3 and a viscosity of 0.002 Pa s, is injected into a hypothetical two-dimensional aquifer. The simulation domain is 50-m long and 10-m high. A total of 25 and 20 uniform nodes are used in the x- and z-direction, respectively. The water table is at approximately 6 m from the surface. The NAPL is injected with a rate of 0.25 m/day at the top of node (13,1,20) for a period of 10 days. After the infiltration period, the NAPL is allowed to redistribute for 190 days. The total simulation time is 200 days.

Part 2:
This problem differs from the problem in Part 1 in that the steady-state initial conditions prior to the NAPL infiltration cannot be specified in the Initial Conditions Card because precipitation (100 ml/year) is included. When more complex initial conditions are required, a separate run has to be completed that yields a restart file that contains the proper initial conditions.

Back to Problem Table