Steady state analysis¶

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


		(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:


		(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:


		(206)

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


		(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:


		(208)