25.01.2013 Views

Hydro-Mechanical Properties of an Unsaturated Frictional Material

Hydro-Mechanical Properties of an Unsaturated Frictional Material

Hydro-Mechanical Properties of an Unsaturated Frictional Material

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

9.3. SETUP OF NUMERICAL SIMULATIONS 187<br />

9.3 Setup <strong>of</strong> Numerical Simulations<br />

For the numerical simulation, a regular rect<strong>an</strong>gular FE mesh is used. The spacing <strong>of</strong> the<br />

grid is 5 mm along the z-axis (direction <strong>of</strong> gravity). The width <strong>of</strong> the finite elements on the<br />

x-axis is chosen such that the bottom area <strong>of</strong> the discretized model domain corresponds to the<br />

bottom area <strong>of</strong> the cylindrical s<strong>an</strong>d column described in the previous section. In the numerical<br />

simulation, although the problem is one-dimensional, the model domain has a depth <strong>of</strong> 1.0 m.<br />

Thus, the width <strong>of</strong> the domain in the numerical model is chosen to be 0.073062 m (see Fig.<br />

9.2). A timestep <strong>of</strong> 4 s is chosen.<br />

The initial <strong>an</strong>d boundary conditions in the numerical simulation correspond to those <strong>of</strong><br />

the experiment. The primary variables are the water pressure uw <strong>an</strong>d the air saturation Sn.<br />

The specimen is initially saturated <strong>an</strong>d thus the effective saturation <strong>of</strong> water in the model<br />

domain is equal to 1.0. A hydrostatic distribution is given for the water pressure. At the top<br />

boundary, a Dirichlet boundary condition is applied for the water <strong>an</strong>d the air phase. This<br />

boundary condition defines directly a water pressure value equal to −3400 Pa at the top <strong>of</strong><br />

the column. This water pressure is equal to the maximum value <strong>of</strong> matric suction measured<br />

during the s<strong>an</strong>d column test I experiment. In combination with the Dirichlet boundary value<br />

1.0 for the air saturation, the definition <strong>of</strong> phase pressure difference (ψ = ua − uw) delivers<br />

that the pressure <strong>of</strong> air at the top boundary is equal to zero (atmospheric condition). The<br />

applied Dirichlet boundary condition at the top <strong>of</strong> the numerical model does not avoid water<br />

to flow through this boundary <strong>an</strong>d therefore might introduce <strong>an</strong> error in the mass bal<strong>an</strong>ce.<br />

To prevent <strong>an</strong>y mass <strong>of</strong> water from leaving the domain through the top boundary, the wetting<br />

phase relative permeability at the finite elements <strong>of</strong> the upper layer is explicitly set to zero.<br />

Outflow mass (kg)<br />

12<br />

10<br />

8<br />

6<br />

4<br />

2<br />

0<br />

0 20000 40000 60000 80000 100000<br />

Time (s)<br />

Figure 9.2: Domain <strong>an</strong>d boundary conditions (left) <strong>an</strong>d water outflow <strong>an</strong>d inflow (right) in<br />

the numerical simulation

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!