boundary conditions

Slicing of the 3D Mandelbrot set, and analysis.

February 8, 2021 math and physics play , , , , , , ,

Some slices.

As followup to:

here is a bit more experimentation with Paraview slice filtering. This time I saved all the point data with the escape time counts, and rendered it with a few different contours

The default slice filter places the plane in the x-y orientation:

but we can also tilt it in the Paraview render UI

and suppress the contour view to see just the slice

As a very GUI challenged user, I don’t find the interface particularly intuitive, but have at least figured out this one particular slicing task, which is kind of cool.  It’s impressive that the UI can drive interesting computational tasks without having to regenerate or reload any of the raw data itself.  This time I was using the MacOSX Paraview client, which is nicer looking than the Windows version, but has some weird glitches in the file dialogues.

Analysis.

The graphing play above shows some apparent rotational symmetry our vector equivalent to the Mandelbrot equation
\begin{equation}
\Bx \rightarrow \Bx \Be_1 \Bx + \Bc.
\end{equation}
It was not clear to me if this symmetry existed, as there were artifacts in the plots that made it appear that there was irregularity. However, some thought shows that this irregularity is strictly due to sampling error, and perhaps also due to limitations in the plotting software, as such an uneven surface is probably tricky to deal with.

To see this, here are the first few iterations of the Mandlebrot sequence for an arbitary starting vector \( \Bc \).
\begin{equation}
\begin{aligned}
\Bx_0 &= \Bc \\
\Bx_1 &= \Bc \Be_1 \Bc + \Bc \\
\Bx_2 &= \lr{ \Bc \Be_1 \Bc + \Bc } \Be_1 \lr{ \Bc \Be_1 \Bc + \Bc } + \Bc \Be_1 \Bc + \Bc.
\end{aligned}
\end{equation}

Now, what happens when we rotate the starting vector \( \Bc \) in the \( y-z \) plane. The rotor for such a rotation is
\begin{equation}
R = \exp\lr{ e_{23} \theta/2 },
\end{equation}
where
\begin{equation}
\Bc \rightarrow R \Bc \tilde{R}.
\end{equation}
Observe that if \( \Bc \) is parallel to the x-axis, then this rotation leaves the starting point invariant, as \( \Be_1 \) commutes with \( R \). That is
\begin{equation}
R \Be_1 \tilde{R} =
\Be_1 R \tilde{R} = \Be_1.
\end{equation}
Let \( \Bc’ = R \Bc \tilde{R} \), so that
\begin{equation}
\Bx_0′ = R \Bc \tilde{R} = R \Bx_0 \tilde{R} .
\end{equation}
\begin{equation}
\begin{aligned}
\Bx_1′
&= R \Bc \tilde{R} \Be_1 R \Bc \tilde{R} + R \Bc \tilde{R} \\
&= R \Bc \Be_1 \Bc \tilde{R} + R \Bc \tilde{R} \\
&= R \lr{ \Bc \Be_1 \Bc R + \Bc } \tilde{R} \\
&= R \Bx_1 \tilde{R}.
\end{aligned}
\end{equation}
\begin{equation}\label{eqn:m2:n}
\begin{aligned}
\Bx_2′
&= \Bx_1′ \Be_1 \Bx_1′ + \Bc’ \\
&= R \Bx_1 \tilde{R} \Be_1 R \Bx_1 \tilde{R} + R \Bc \tilde{R} \\
&= R \Bx_1 \Be_1 \Bx_1 \tilde{R} + R \Bc \tilde{R} \\
&= R \lr{ \Bx_1 \Be_1 \Bx_1 + \Bc } \tilde{R} \\
&= R \Bx_2 \tilde{R}.
\end{aligned}
\end{equation}

The pattern is clear. If we rotate the starting point in the y-z plane, iterating the Mandelbrot sequence results in precisely the same rotation of the x-y plane Mandelbrot sequence. So the apparent rotational symmetry in the 3D iteration of the Mandelbrot vector equation is exactly that. This is an unfortunately boring 3D fractal. All of the interesting fractal nature occurs in the 2D plane, and the rest is just a consequence of rotating that image around the x-axis. We get some interesting fractal artifacts if we slice the rotated Mandelbrot image.

Maxwell equation boundary conditions

September 6, 2016 math and physics play , , , , , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

Motivation

boundaryConditionsTwoSurfacesFig1

fig 1. Two surfaces normal to the interface.

Most electrodynamics textbooks either start with or contain a treatment of boundary value conditions. These typically involve evaluating Maxwell’s equations over areas or volumes of decreasing height, such as those illustrated in fig. 1, and fig. 2. These represent surfaces and volumes where the height is allowed to decrease to infinitesimal levels, and are traditionally used to find the boundary value constraints of the normal and tangential components of the electric and magnetic fields.

boundaryConditionsPillBoxFig2

fig 2. A pillbox volume encompassing the interface.

More advanced topics, such as evaluation of the Fresnel reflection and transmission equations, also rely on similar consideration of boundary value constraints. I’ve wondered for a long time how the Fresnel equations could be attacked by looking at the boundary conditions for the combined field \( F = \BE + I c \BB \), instead of the considering them separately.

A unified approach.

The Geometric Algebra (and relativistic tensor) formulations of Maxwell’s equations put the electric and magnetic fields on equal footings. It is in fact possible to specify the boundary value constraints on the fields without first separating Maxwell’s equations into their traditional forms. The starting point in Geometric Algebra is Maxwell’s equation, premultiplied by a stationary observer’s timelike basis vector

\begin{equation}\label{eqn:maxwellBoundaryConditions:20}
\gamma_0 \grad F = \inv{\epsilon_0 c} \gamma_0 J,
\end{equation}

or

\begin{equation}\label{eqn:maxwellBoundaryConditions:40}
\lr{ \partial_0 + \spacegrad} F = \frac{\rho}{\epsilon_0} – \frac{\BJ}{\epsilon_0}.
\end{equation}

The electrodynamic field \(F = \BE + I c \BB\) is a multivector in this spatial domain (whereas it is a bivector in the spacetime algebra domain), and has vector and bivector components. The product of the spatial gradient and the field can still be split into dot and curl components \(\spacegrad M = \spacegrad \cdot M + \spacegrad \wedge M \). If \(M = \sum M_i \), where \(M_i\) is an grade \(i\) blade, then we give this the Hestenes’ [1] definitions

\begin{equation}\label{eqn:maxwellBoundaryConditions:60}
\begin{aligned}
\spacegrad \cdot M &= \sum_i \gpgrade{\spacegrad M_i}{i-1} \\
\spacegrad \wedge M &= \sum_i \gpgrade{\spacegrad M_i}{i+1}.
\end{aligned}
\end{equation}

With that said, Maxwell’s equation can be rearranged into a pair of multivector equations

\begin{equation}\label{eqn:maxwellBoundaryConditions:80}
\begin{aligned}
\spacegrad \cdot F &= \gpgrade{-\partial_0 F + \frac{\rho}{\epsilon_0} – \frac{\BJ}{\epsilon_0 c}}{0,1} \\
\spacegrad \wedge F &= \gpgrade{-\partial_0 F + \frac{\rho}{\epsilon_0} – \frac{\BJ}{\epsilon_0 c}}{2,3},
\end{aligned}
\end{equation}

The latter equation can be integrated with Stokes theorem, but we need to apply a duality transformation to the latter in order to apply Stokes to it

\begin{equation}\label{eqn:maxwellBoundaryConditions:120}
\begin{aligned}
\spacegrad \cdot F
&=
-I^2 \spacegrad \cdot F \\
&=
-I^2 \gpgrade{\spacegrad F}{0,1} \\
&=
-I \gpgrade{I \spacegrad F}{2,3} \\
&=
-I \spacegrad \wedge (IF),
\end{aligned}
\end{equation}

so

\begin{equation}\label{eqn:maxwellBoundaryConditions:100}
\begin{aligned}
\spacegrad \wedge (I F) &= I \lr{ -\inv{c} \partial_t \BE + \frac{\rho}{\epsilon_0} – \frac{\BJ}{\epsilon_0 c} } \\
\spacegrad \wedge F &= -I \partial_t \BB.
\end{aligned}
\end{equation}

Integrating each of these over the pillbox volume gives

\begin{equation}\label{eqn:maxwellBoundaryConditions:140}
\begin{aligned}
\oint_{\partial V} d^2 \Bx \cdot (I F)
&=
\int_{V} d^3 \Bx \cdot \lr{ I \lr{ -\inv{c} \partial_t \BE + \frac{\rho}{\epsilon_0} – \frac{\BJ}{\epsilon_0 c} } } \\
\oint_{\partial V} d^2 \Bx \cdot F
&=
– \partial_t \int_{V} d^3 \Bx \cdot \lr{ I \BB }.
\end{aligned}
\end{equation}

In the absence of charges and currents on the surface, and if the height of the volume is reduced to zero, the volume integrals vanish, and only the upper surfaces of the pillbox contribute to the surface integrals.

\begin{equation}\label{eqn:maxwellBoundaryConditions:200}
\begin{aligned}
\oint_{\partial V} d^2 \Bx \cdot (I F) &= 0 \\
\oint_{\partial V} d^2 \Bx \cdot F &= 0.
\end{aligned}
\end{equation}

With a multivector \(F\) in the mix, the geometric meaning of these integrals is not terribly clear. They do describe the boundary conditions, but to see exactly what those are, we can now resort to the split of \(F\) into its electric and magnetic fields. Let’s look at the non-dual integral to start with

\begin{equation}\label{eqn:maxwellBoundaryConditions:160}
\begin{aligned}
\oint_{\partial V} d^2 \Bx \cdot F
&=
\oint_{\partial V} d^2 \Bx \cdot \lr{ \BE + I c \BB } \\
&=
\oint_{\partial V} d^2 \Bx \cdot \BE + I c d^2 \Bx \wedge \BB \\
&=
0.
\end{aligned}
\end{equation}

No component of \(\BE\) that is normal to the surface contributes to \(d^2 \Bx \cdot \BE \), whereas only components of \(\BB\) that are normal contribute to \(d^2 \Bx \wedge \BB \). That means that we must have tangential components of \(\BE\) and the normal components of \(\BB\) matching on the surfaces

\begin{equation}\label{eqn:maxwellBoundaryConditions:180}
\begin{aligned}
\lr{\BE_2 \wedge \ncap} \ncap – \lr{\BE_1 \wedge (-\ncap)} (-\ncap) &= 0 \\
\lr{\BB_2 \cdot \ncap} \ncap – \lr{\BB_1 \cdot (-\ncap)} (-\ncap) &= 0 .
\end{aligned}
\end{equation}

Similarly, for the dot product of the dual field, this is

\begin{equation}\label{eqn:maxwellBoundaryConditions:220}
\begin{aligned}
\oint_{\partial V} d^2 \Bx \cdot (I F)
&=
\oint_{\partial V} d^2 \Bx \cdot (I \BE – c \BB) \\
&=
\oint_{\partial V} I d^2 \Bx \wedge \BE – c d^2 \Bx \cdot \BB.
\end{aligned}
\end{equation}

For this integral, only the normal components of \(\BE\) contribute, and only the tangential components of \(\BB\) contribute. This means that

\begin{equation}\label{eqn:maxwellBoundaryConditions:240}
\begin{aligned}
\lr{\BE_2 \cdot \ncap} \ncap – \lr{\BE_1 \cdot (-\ncap)} (-\ncap) &= 0 \\
\lr{\BB_2 \wedge \ncap} \ncap – \lr{\BB_1 \wedge (-\ncap)} (-\ncap) &= 0.
\end{aligned}
\end{equation}

This is why we end up with a seemingly strange mix of tangential and normal components of the electric and magnetic fields. These constraints can be summarized as

\begin{equation}\label{eqn:maxwellBoundaryConditions:260}
\begin{aligned}
( \BE_2 – \BE_1 ) \cdot \ncap &= 0 \\
( \BE_2 – \BE_1 ) \wedge \ncap &= 0 \\
( \BB_2 – \BB_1 ) \cdot \ncap &= 0 \\
( \BB_2 – \BB_1 ) \wedge \ncap &= 0
\end{aligned}
\end{equation}

These relationships are usually expressed in terms of all of \(\BE, \BD, \BB\) and \(\BH \). Because I’d started with Maxwell’s equations for free space, I don’t have the \( \epsilon \) and \( \mu \) factors that produce those more general relationships. Those more general boundary value relationships are usually the starting point for the Fresnel interface analysis. It is also possible to further generalize these relationships to include charges and currents on the surface.

References

[1] D. Hestenes. New Foundations for Classical Mechanics. Kluwer Academic Publishers, 1999.