SimPEG.electromagnetics.natural_source.Simulation1DPrimarySecondary

class SimPEG.electromagnetics.natural_source.Simulation1DPrimarySecondary(mesh, survey=None, sigmaPrimary=None, **kwargs)

Bases: SimPEG.electromagnetics.natural_source.simulation.Simulation1DElectricField

A NSEM problem soving a e formulation and primary/secondary fields decomposion.

By eliminating the magnetic flux density using

$$\mathbf{b}=\frac{1}{i\omega }(-\mathbf{C}\mathbf{e})$$

we can write Maxwell’s equations as a second order system in $\mathbf{e}$ only:

$$[{\mathbf{C}}^{\mathrm{\top }}{\mathbf{M}}_{{\mu }^{\mathbf{-}\mathbf{1}}}^{\mathbf{e}}\mathbf{C}+i\omega {\mathbf{M}}_{\sigma }^{\mathbf{f}}]{\mathbf{e}}_s=i\omega {\mathbf{M}}_{{\sigma }_{\mathbf{s}}}^{\mathbf{f}}{\mathbf{e}}_p$$

which we solve for ${\mathbf{e}}_{\mathbf{s}}$. The total field $\mathbf{e}={\mathbf{e}}_{\mathbf{p}}+{\mathbf{e}}_{\mathbf{s}}$.

The primary field is estimated from a background model (commonly half space ).

Attributes

sigmaPrimary

A background model, use for the calculation of the primary fields.

Methods

fieldsPair	alias of SimPEG.electromagnetics.natural_source.fields.Fields1DPrimarySecondary
getADeriv(freq, u, v[, adjoint])	The derivative of A wrt sigma
getRHS(freq)	Function to return the right hand side for the system.
getRHSDeriv(freq, src, v[, adjoint])	The derivative of the RHS wrt sigma