T3P is a 3D parallel finite-element time-domain solver to calculate transient field response of a electromagnetic structure to imposed fields, and dipole or beam excitations.
In our approach, Ampere's and Faraday's laws are combined and integrated over time to yield the inhomogeneous vector wave equation for the time integral of the electric field E:
Unknown macro: {latex} \begin
Unknown macro: {eqnarray*} \left( \epsilon \frac
Unknown macro: {partial^2}
Unknown macro: {partial t^2}
+ \sigma_
Unknown macro: {eff}
\frac
Unknown macro: {partial}
Unknown macro: {partial t}
+ \nabla \times \mu^{-1}\nabla\times \right) \int^t \emph
Unknown macro: {E}
d \tau = -\emph
Unknown macro: {J}
\end
with permittivity
Unknown macro: {latex} $\epsilon = \epsilon_0 \epsilon_r $
and permeability
Unknown macro: {latex}
$\mu = \mu_0 \mu_r$
. In the current implementation, a constant value of the effective conductivity
Unknown macro: {latex}
$\sigma_
Unknown macro: {eff}
= \tan \delta \cdot 2 \pi f \cdot \epsilon $
is assumed by fixing a frequency
Unknown macro: {latex}
$f$
, and the losses are specified by the loss tangent
Unknown macro: {latex} $\tan \delta $
. As is common for Wakefield computations of rigid beams, the electric current source density J is given by a one-dimensional Gaussian particle distribution, moving at the speed of light along the beam line.
The computational domain is discretized into curved tetrahedral elements and
Unknown macro: {latex} $\int^t \emph
Unknown macro: {E}
d \tau $
in Eq. (1) is represented as an expansion in hierarchical Whitney vector basis functions
Unknown macro: {latex} $ \emph
Unknown macro: {N}
_i $
Unknown macro: {latex}
\begin
Unknown macro: {eqnarray*} \int^t \emph
Unknown macro: {E}
(\emph
Unknown macro: {x}
, \tau) d \tau = \sum^
Unknown macro: {N_p}
_
Unknown macro: {i=1}
e_i (t) \cdot \emph
Unknown macro: {N}
_i \end
up to order
Unknown macro: {latex} $p$
within each element. For illustration, the numbers of basis functions for first, second and sixth order are
Unknown macro: {latex} $N_1 = 6 $
,
Unknown macro: {latex} $N_2 = 20 $
and
Unknown macro: {latex} $N_6 = 216 $
, respectively. Higher order elements (both curved and with higher-order basis functions) not only significantly improve field accuracy and dispersive properties, but also generally lead to higher-order accurate particle-field coupling equivalent to, but much less laborious than, complicated higher-order interpolation schemes commonly found in finite-difference methods.