SimPEG.electromagnetics.static.spontaneous_potential.HydraulicHeadMap#

class SimPEG.electromagnetics.static.spontaneous_potential.HydraulicHeadMap(mesh, L)[source]#

Bases: LinearMap

Maps hydraulic head to charge density accumulation rate.

Parameters:
meshdiscretize.base.BaseMesh
Lfloat or (mesh.n_cells,) array_like

Cross coupling property model (A/m^2).

Notes

This map relates hydraulic head, \(h\), in m to the volumetric current density, \(q\), (A/m^3) using a cross coupling term to convert the negative gradient of hydrualic head into secondary current densities. This cross coupling parameter \(L\) is a multiplication of the hydraulic conductivity (m/s) of the material and the ionic density (C/m^3) of the fluid.

The mapping is:

\[q = \nabla \cdot L \nabla h\]