The finite-difference time-domain (FDTD) method requires a minimum of ten cells per wavelength to assume a linear field transition between adjacent cells in order to consistently achieve accuracy. Maintaining this relationship for good conductors is not practical because the fields are decaying exponentially over the skin depth, which is a short distance when compared to the free space wavelength. The surface conductivity correction in XF can be enabled for good conductors to more accurately model the losses and reflections at the boundary of a metal.

Background

FDTD's discrete representation of Maxwell's equations creates a need for XF's surface connectivity correction. In order to accurately model field interactions at metal boundaries, the conductivity needs to be modified so the impedance used in FDTD matches the theoretical value.

In this case, assume a time-harmonic transverse electromagnetic wave propagating in lossy media. Further, simplify the wave to traveling in the z direction and containing only the x component of the electric field.

For the theoretical case, the equation for the electric field is written as

\begin{equation} {\bf E}(z) = \hat a_xE_x(z) = \hat a_x\left(E_0^+e^{-\gamma z}+E_0^-e^{+\gamma z}\right) \end{equation}

where the two terms represent the positive and negative traveling waves. $\gamma$ is the propagation constant and takes the form

\begin{equation} \gamma=\sqrt{-\omega^2\mu\epsilon+j\omega\mu\sigma_t} \end{equation}

$\sigma_t$ is the theoretical conductivity of the material. The corresponding magnetic field is determined from Maxwell's equations as

\begin{eqnarray} {\bf H}&=&-\frac{1}{j\omega\mu}\nabla\times{\bf E} \\ &=&-\hat a_y\frac{1}{j\omega\mu}\frac{\partial E_x}{\partial z} \\ &=&\hat a_y\frac{\gamma}{j\omega\mu}\left(E_0^+e^{-\gamma z}-E_0^-e^{+\gamma z}\right) \\ &=&\hat a_y\sqrt{\frac{\sigma_t+j\omega\epsilon}{j\omega\mu}}\left(E_0^+e^{-\gamma z}-E_0^-e^{+\gamma z}\right) \end{eqnarray}

From the above equation, the theoretical wave impedance, or ratio of electric to magnetic field, is represented as

\begin{equation} Z_t=\sqrt{\frac{j\omega\mu}{\sigma_t+j\omega\epsilon}}=\eta_t \end{equation}

This is equivalent to the intrinsic impedance, $\eta_t$ of the lossy medium.

FDTD approximates the partial derivatives in Maxwell's equation by discretizing them, causing the representation of the material's intrinsic impedance to be different. Starting with Maxwell's equation, the z directed wave is represented

\begin{equation} \frac{\partial H_y}{\partial z} = \epsilon \frac{\partial E_x}{\partial t} + \sigma_f E_x \end{equation}

where $\sigma_f$ is the conductivity of the material in FDTD. In partial FDTD notation, Maxwell's equation takes the form

\begin{equation} \frac{\left(H_y\big|_{i,j,k+1/2}-H_y\big|_{i,j,k-1/2}\right)}{\Delta z} = \epsilon\frac{E_x\big|^{n-1/2}_{i,j,k} - E_x\big|^{n+1/2}_{i,j,k}}{\Delta t} + \sigma_f E_x\big|_{i,j,k} \end{equation}

where $i,j,k$ represent the spatial index of the fields. Notice the use of $z-1/2$ because the magnetic fields are separated from the electric fields in space by a half cell.

Rearranging the terms and simplifying the equation with assumptions for a good conductor allows the relation between the electric and magnetic fields to be written as

\begin{equation} E_x\big|_{i,j,k}\simeq\frac{\Delta t}{(\epsilon+\sigma_f\Delta t)}\left[\frac{H_y\big|_{i,j,k+1/2}}{\Delta z}\right] \end{equation}

Since $\sigma\Delta t\gg\epsilon$ for good conductors, the intrinsic impedance in FDTD is written

\begin{equation} \eta_f\simeq\frac{1}{\sigma_f\Delta z} \end{equation}

In order to accurately model reflections and losses at good conductor boundaries, XF's surface conductivity correction setting is necessary for matching the impedance used in FDTD to that from theory, $\eta_f = \eta_t$.

\begin{equation} \frac{1}{\sigma_f\Delta z} = \sqrt{\frac{j\omega\mu}{\sigma_t+j\omega\epsilon}} \end{equation}

Since FDTD is a time-domain formulation, the impedance can be maintained by solving for the real part of $\sigma_f$.

\begin{equation} \sigma_f=\frac{1}{\Delta z}\sqrt{\frac{2\sigma_t}{\omega\mu}} \label{eq:sigmaf} \end{equation}

The surface conductivity correction setting adjusts the material conductivity in the calculation engine to $\sigma_f$ from the theoretical value, $\sigma_t$, specified by the user in the user interface.

Material Editor

Surface Conductivity Connection enabled.

The Surface Conductivity Correction setting is available in the Material Editor and should be checked when defining the material. Regardless of the Entry Method, the theoretical values for the material properties should be entered. The Evaluation Frequency determines $\omega$ in the previous equation for maintaining the impedance. The calculation engine determines and uses $\sigma_f$.