simpeg.electromagnetics.time_domain.Simulation3DCurrentDensity.getRHS#

Simulation3DCurrentDensity.getRHS(tInd)[source]#

Right-hand sides for the given time index.

This method returns the right-hand sides for the time index provided. The right-hand side for each source is constructed according to:

\[\mathbf{q}_k = - \frac{1}{\Delta t_k} \big [ \mathbf{s}_{\mathbf{e},k} + \mathbf{s}_{\mathbf{e},k-1} \big ] \; - \; \frac{1}{\Delta t_k} \mathbf{C M_{e\mu}^{-1}} \big [ \mathbf{s}_{\mathbf{m},k} + \mathbf{s}_{\mathbf{m},k-1} \big ]\]

where

  • \(\Delta t_k\) is the step length

  • \(\mathbf{C}\) is the discrete curl operator

  • \(\mathbf{s_m}\) and \(\mathbf{s_e}\) are the integrated magnetic and electric source terms, respectively

  • \(\mathbf{M_{e\mu}}\) is the inner-product matrix for permeabilities projected to edges

See the Notes section of the doc strings for Simulation3DCurrentDensity for a full description of the formulation.

Parameters:
tIndint

The time index; between [0, n_steps].

Returns:
(n_faces, n_sources) numpy.ndarray

The right-hand sides.