choclo.prism.magnetic_en

Contents

choclo.prism.magnetic_en#

choclo.prism.magnetic_en(easting, northing, upward, prism_west, prism_east, prism_south, prism_north, prism_bottom, prism_top, magnetization_east, magnetization_north, magnetization_up)[source]#

Northing derivative of the easting component of the magnetic field.

Returns the northing derivative of the easting component of the magnetic field due to a single rectangular prism on a single computation point.

Parameters:
easting, northing, upwardfloat

Easting, northing and upward coordinates of the observation point. Must be in meters.

prism_west, prism_east, prism_south, prism_north, prism_bottom, prism_topfloat

The boundaries of the prism. Must be in meters.

magnetization_eastfloat

The East component of the magnetization vector of the prism. Must be in \(A m^{-1}\).

magnetization_northfloat

The North component of the magnetization vector of the prism. Must be in \(A m^{-1}\).

magnetization_upfloat

The upward component of the magnetization vector of the prism. Must be in \(A m^{-1}\).

Returns:
b_enfloat

Northing derivative of the easting component of the magnetic field generated by the prism on the observation point in \(\text{T}\). Return numpy.nan if the observation point falls in a singular point: prism vertices, prism edges or interior points.

Notes

Computes the northing derivative of the easting component of the magnetic field \(\mathbf{B}(\mathbf{p})\) generated by a rectangular prism \(R\) with a magnetization vector \(M\) on the observation point \(\mathbf{p}\) as follows:

\[B_{xy}(\mathbf{p}) = \frac{\mu_0}{4\pi} \left( M_x u_{xxy} + M_y u_{xyy} + M_z u_{xyz} \right)\]

Where \(u_{ijk}\) are:

\[u_{ijk} = \frac{\partial^3}{\partial i \partial j \partial k} \int\limits_R \frac{1}{\lVert \mathbf{p} - \mathbf{q} \rVert} dv\]

with \(i,j,k \in \{x, y, z\}\).

References