SimPEG.potential_fields.gravity.Simulation3DDifferential

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

Bases: SimPEG.base.pde_simulation.BasePDESimulation

Finite volume simulation class for gravity.

Notes

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

$${\mathrm{\nabla }}^2\varphi =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:

$$[\mathbf{D}{\mathbf{M}}_{\mathbf{f}}{\mathbf{D}}^{\mathbf{T}}]\mathbf{u}=-{\mathbf{M}}_{\mathbf{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