\phi \cdot c_t \cdot \mu \cdot \partial_t \Psi -   
 \nabla \cdot \left( k \cdot  \Big( \vec \nabla \Psi  - \frac{\rho^2}{\mu}  \, {\bf g} \Big)  \right) 
  = \sum_k q({\bf r}) \cdot \delta({\bf r}-{\bf r}_k)

where

	reservoir pressure		time
	fluid density		position vector
	effective porosity		position vector of the -th source
	total compressibility		Dirac delta function
	sandface flowrates of the -th source		gradient operator
	phase mobility		gravity vector
	formation permeability to a given fluid		fluid velocity under Darcy flow
	dynamic viscosity of a given fluid
	fluid compressibility factor
	Pseudo-Pressure

Derivation of Single-phase pseudo-pressure diffusion @model

See also