# SimPEG.potential_fields.gravity.Simulation3DDifferential#

class SimPEG.potential_fields.gravity.Simulation3DDifferential(mesh, rho=1.0, rhoMap=None, **kwargs)[source]#

Finite volume simulation class for gravity.

Notes

From Blakely (1996), the scalar potential $$\phi$$ outside the source region is obtained by solving a Poisson’s equation:

$\nabla^2 \phi = 4 \pi \gamma \rho$

where $$\gamma$$ is the gravitational constant and $$\rho$$ defines the distribution of density within the source region.

Applying the finite volumn method, we can solve the Poisson’s equation on a 3D voxel grid according to:

$\big [ \mathbf{D M_f D^T} \big ] \mathbf{u} = - \mathbf{M_c \, \rho}$

Attributes

 rho Specific density (g/cc) physical property model. rhoDeriv Derivative of Specific density (g/cc) wrt the model. rhoMap Mapping of the inversion model to Specific density (g/cc).

Methods

 fields([m]) Compute fields GetA creates and returns the A matrix for the Gravity nodal problem Return right-hand side for the linear system