:uid: PSfdem :title: Primary-Secondary Approaches for Frequency Domain Maxwell :description: 1+1=2, but sometimes it is more useful to think of 2=a+b where a=... :tooltip:/img/twoab :tag: em,fdem,finite volume,primary-secondary :group: simpeg :license: CC-BY-4.0 :source: https://api.github.com/gists/b10a48dfea9af70aeb52
When solving any math or physics problem, there are usually a few tricks that can take you a long ways. Some of the favorites are: multiplying by 1 (or the identity), adding zero, and breaking things up in to sums (i.e. 2 = a+b, where a=1). I am sure these will be recurring themes among our articles. Here, I am going to spend some time on last point in the context of the frequency domain Maxwell's equations. In this case, 2 = a+b goes by the name of Primary-Secondary. I will develop this for the frequency domain electromagnetic (FDEM) problem, but keep in mind that it is quite general.
The FDEM is governed by Maxwell's equations, which, under the the quasi-static approximation are given by:
where
There are a large number of options for primary-secondary: 2 = a + b leaves the door wide open... there are actually an infinite number of possibilities, but not all of them will be useful. I will give you a couple of examples, and hopefully that will give you a feel for where this can be handy.
All we are going to do here is define
$\sigma = \sigma_p + \sigma_s$ $\vec{E} = \vec{E}_p + \vec{E}_s$ $\vec{B} = \vec{B}_p + \vec{B}_s$
where the superscript
We haven't actually done anything yet. To make this a useful approach, we need to define our primary. We are free to choose anything, the idea here is that we gain something (ie. easier or faster to compute, or the ability to take advantage of an analytical solution), so a useful choice turns out to be choosing a primary which satisfies:
where
we expand it out:
and use our definition of the primary:
which is a system of equations in the secondary fields with the source
So what we have done is replace our source term with this new source:
Often when using a magnetic dipole source, it is useful to use a zero-frequency primary, then we can define the source analytically as a curl of a vector potential. This ensures that the divergence of the magnetic field is zero (ie. we are not numerically creating magnetic monopoles) if we are using a mimetic discretization. We start with:
Here, we choose a primary such that $$ \nabla \times \mu^{-1} \vec{B}_p = \vec{S}_e $$
and
Since
For the case of a static magnetic dipole in a space with homogeneous permeability
See for example: Griffiths, 1999 section 5.4
Using this, we define
Now to set-up the secondary problem, we play the same game as before and split up the fields in Maxwell's equations into primary and secondary:
using the definition of the primary
That is, we have a frequency-dependent magnetic source term for the secondary problem.
A primary-secondary approach is often used in Magnetotellurics (MT) to simplify the boundary conditions. If performing a full 3D simulation for MT, the boundary conditions drive the system, and must be accounted for in the definition of the discrete differential operators. However, primary-secondary approach can be used to allow natural boundary conditions (ie. the fields have sufficiently decayed before encountering the boundary) to be employed. A 1D background is often assumed for the primary, and the primary fields can be computed analytically using a propagator matrix approach (see Ward and Hohmann, 1988 page 194), or with a 1D forward simulation. The source term is found using the Primary-Secondary approach for conductivity shown above.
Recently, I have been looking at using Primary-Secondary in settings with steel-cased wells. If the source is inside of the well, and the geologic setting is a layered 1D earth, then the problem is cylindrically symmetric and a 2D problem can be solved. However, if there are 3D structures, a primary-secondary approach can be used: solve the cylindrically symmetric problem first, and use these as a primary for the 3D problem. If you want to know more... shameless self-plug: I will be presenting at SEG in New Orleans on Tuesday, October 20 at 9:45 in the Borehole Geophysics 1 session.
We have broken the problem up into two: the primary and the secondary. The primary fields serve as a source for the secondary, thus connecting the two problems. Breaking this problem up in to two steps affords us some freedom as to how we choose to solve each of these problems. It is note necessarily the approach for every problem, but where symmetry can be exploited or analytics employed, it can prove to be a useful approach.