choclo.prism.magnetic_eu#
- choclo.prism.magnetic_eu(easting, northing, upward, prism_west, prism_east, prism_south, prism_north, prism_bottom, prism_top, magnetization_east, magnetization_north, magnetization_up)[source]#
Upward derivative of the easting component of the magnetic field.
Returns the upward derivative of the easting component of the magnetic field due to a single rectangular prism on a single computation point.
- Parameters:
- easting, northing, upward
float
Easting, northing and upward coordinates of the observation point. Must be in meters.
- prism_west, prism_east, prism_south, prism_north, prism_bottom, prism_top
float
The boundaries of the prism. Must be in meters.
- magnetization_east
float
The East component of the magnetization vector of the prism. Must be in \(A m^{-1}\).
- magnetization_north
float
The North component of the magnetization vector of the prism. Must be in \(A m^{-1}\).
- magnetization_up
float
The upward component of the magnetization vector of the prism. Must be in \(A m^{-1}\).
- easting, northing, upward
- Returns:
- b_eu
float
Upward 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.
- b_eu
See also
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_{xz}(\mathbf{p}) = \frac{\mu_0}{4\pi} \left( M_x u_{xxz} + M_y u_{xyz} + M_z u_{xzz} \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