SimPEG.potential_fields.gravity.Simulation3DDifferential#

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

Bases: SimPEG.base.pde_simulation.BasePDESimulation

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

GetA creates and returns the A matrix for the Gravity nodal problem

getRHS()

Return right-hand side for the linear system