Steady state analysis

<< Click to Display Table of Contents >>

Navigation:  »No topics above this level«

Steady state analysis

 

The gouverning partial differential equation for a confined aquifier with flow in the horizontal (x,y) plane is:

 

 

p_70_image004

(204)

                               

Where ϕ is the fluid potential or total head measured from the bottom of the aquifier, kx and ky are the coefficients of the permeability in the x and y directions and Q is the recharge. Pumping is a negative Q. The finite element discretization process reduces the differential equation to a set of equilibrium type equations of the form:

 

 

p_70_image006

(205)

                                             

Where kc is the symmetrical conductivity matrix, Φ is a vector of nodal potential (total head) values, and q is a vector of nodal inflows/outflows. The element matrices are calculated using:

 

 

p_71_image002

(206)

                             

Assuming that the principal axes of the permeability tensor coincide with x and y, the property matrix K is:

 

 

p_71_image004

(207)

                                               

and the T matrix is similar with the strain-displacement matrix B in the stress analysis. For instance for 4 noded quadrilateral element the matrix T is given by:

 

 

p_71_image006

 

 

(208)

 

 

 

 

 

 

 

 

 

 

© 2020 GeoStru Software