BC problem: modelling ground as a porous material
 

Hi,
As shown in the attached figure, ground and well is modeled (axisymmetric) as a porous material using laminar flow, energy equation and some UDF for the porous resistances. There is only one phase (water).
Because of ground water flow from the outside boundary, I choose the pressure inlet BC for the outside and bottom boundary (upper boundary is modeled as a wall).
: 1 velocity inlet, 1 outlet at well surface. porosity = 1 in the well, porosity < 1 in the ground.

Q1: To model the hydrostatic pressure due to the gravity, I checked the gravity box and specified operating density (=0). In that case, I can see a pressure distribution with depth direction. However, if I use 998.2 kg/m3 as the operating density, can I get pressure head through the analysis? Because, in this case, pressure = 0 at the whole domain from the result.
[gradient of groundwater flow is calculated in the total head difference in the aspect of groundwater. total head = potential head (due to the gravity, location) + pressure head + dynamic head (in most cases, = 0) ]


Q2: For the outside and bottom boundary, I used a gauge pressure = 0 at boundaries. Is that reasonable?:confused:
I can see many reverse flow faces in this case at the pressure inlet boundary with operating density = 998.2 kg/m3.
Also, I've tried to give a reverse pressure to the boundary for the equal pressure in the ground using simple UDF ( (height - location) * 998.2 * 9.81 ). This is for the equal pressure at the boundary to model pressure head.
How can I model the pressure which is only due to the pressure head not static pressure? Any suggestions?

Thanks :)

Problem solved using UDF, Thanks

excuse me sir, I'm also working in porous media, but I can't definite it. May be can you help me? my email is abdulhalim060189@yahoo.co.id