Skip to content
Get started

1. Continuum and statics

A tank is motionless, yet the pressure on its wall grows with depth. Before simulating moving water, predict that pressure and the force on a vertical gate. This establishes the sign and units of the stresses we will later compute.

A continuum model assigns density ρ\rho, velocity u\boldsymbol u and pressure pp to points. A point represents an average over a region containing many molecules, small compared with the length on which the flow varies. The model loses this interpretation if molecular scales become comparable to the geometry.

For the initial route, assume a single Newtonian fluid with constant density and viscosity, no free-surface motion and no temperature dependence. These assumptions are suitable for studying the equations; they must be reconsidered for a boiling liquid, a rarefied gas or a strongly compressible flow.

Quantity Meaning SI unit
ρ\rho Mass per volume kg/m3\mathrm{kg/m^3}
pp Isotropic compressive stress Pa=kg/(ms2)\mathrm{Pa=kg/(m\,s^2)}
μ\mu Dynamic viscosity Pas\mathrm{Pa\,s}
ν=μ/ρ\nu=\mu/\rho Momentum diffusivity m2/s\mathrm{m^2/s}

Pressure and viscosity are different physical quantities. Pressure can exist at rest. A Newtonian viscous stress requires a velocity gradient, as we derive in Chapter 2.

Choose zz upward and gravitational acceleration g=gez\boldsymbol g=-g\boldsymbol e_z. Take a stationary slice of area AA and height dz\mathrm dz. Pressure below pushes upward, pressure above pushes downward, and weight acts downward:

p(z)Ap(z+dz)AρgAdz=0.p(z)A-p(z+\mathrm dz)A-\rho gA\,\mathrm dz=0.

Dividing by AdzA\,\mathrm dz and taking the limit gives

dpdz=ρg.\frac{\mathrm dp}{\mathrm dz}=-\rho g.

If the free surface is at z=Hz=H and has pressure patmp_{\rm atm}, then

p(z)=patm+ρg(Hz).p(z)=p_{\rm atm}+\rho g(H-z).
surface z = H: p = p_atm
     │
     │ depth h = H − z       gravity ↓
     │
point  z:     p = p_atm + ρgh
     │
bottom z = 0: p = p_atm + ρgH
Depth increases downward while the chosen coordinate increases upward.

For an illustrative constant density 1000kg/m31000\,\mathrm{kg/m^3} and g=9.81m/s2g=9.81\,\mathrm{m/s^2}, the gauge pressure one metre below the surface is 9810Pa9810\,\mathrm{Pa}. This value follows from the balance and chosen inputs; it is not a numerical-solver result.

Consider a vertical rectangular gate of width bb, extending from the surface to depth HH, with atmospheric pressure on its other side. At depth hh, the net force on a strip of height dh\mathrm dh is dF=ρghbdh\mathrm dF=\rho ghb\,\mathrm dh. Therefore

F=0Hρghbdh=ρgbH22.F=\int_0^H\rho ghb\,\mathrm dh=\frac{\rho gbH^2}{2}.

The depth of the resultant follows by equating moments about the surface:

hcp=0HhdFF=2H3.h_{\rm cp}=\frac{\int_0^H h\,\mathrm dF}{F}=\frac{2H}{3}.

For b=H=1mb=H=1\,\mathrm m, the load is 4905N4905\,\mathrm N, acting two thirds of a metre below the surface. The resultant lies below the midpoint because the lower strips carry larger loads.

Gauge pressure subtracts a reference. When both sides of the gate see the same atmosphere, that contribution cancels in the net force. It would not cancel if the back side were evacuated. A constant shift in a mathematical pressure unknown is harmless only when the problem’s boundary data permit the same shift. Absolute pressure also matters for physical questions such as cavitation, outside this incompressible teaching model.

The cylinder workflow later uses a prescribed outlet traction. Do not add an arbitrary pressure shift while holding that traction fixed: you would change the boundary-value problem.

  1. Double the liquid depth. Predict the bottom gauge pressure and gate force.
  2. Integrate the pressure over a gate extending from depths h0h_0 to h1h_1.
  3. Replace the upward coordinate with downward depth. Derive the pressure derivative again before writing a sign.
  4. Someone multiplies bottom pressure by the whole gate area. Why is the result wrong, and by what factor for a gate starting at the surface?

Check your reasoning. In exercise 1, pressure doubles and force quadruples. In exercise 4, the result is twice the correct force because the mean gauge pressure is half the bottom value. These are checks of the model, independent of mesh or software.

For the continuum framework and equilibrium reasoning, consult the instructor’s fluid mechanics notes.

Book map · Next: Mass, momentum and stress