Free Books

Localized Velocity Excitations

Initial velocity excitations are straightforward in the DW paradigm, but can be less intuitive in the FDTD domain. It is well known that velocity in a displacement-wave DW simulation is determined by the difference of the right- and left-going waves [437]. Specifically, initial velocity waves $ v^{\pm}$ can be computed from from initial displacement waves $ y^\pm$ by spatially differentiating $ y^\pm$ to obtain traveling slope waves $ y'^\pm$, multiplying by minus the tension $ K$ to obtain force waves, and finally dividing by the wave impedance $ R=\sqrt{K\epsilon }$ to obtain velocity waves:

$\displaystyle v^{+}$ $\displaystyle =$ $\displaystyle -cy'^{+}= \frac{f^{{+}}}{R}$  
$\displaystyle v^{-}$ $\displaystyle =$ $\displaystyle \;cy'^{-}= -\frac{f^{{-}}}{R},
\protect$ (E.13)

where $ c=\sqrt{K/\epsilon }$ denotes sound speed. The initial string velocity at each point is then $ v(nT,mX)=v^{+}(n-m)+v^{-}(n+m)$. (A more direct derivation can be based on differentiating Eq.$ \,$(E.4) with respect to $ x$ and solving for velocity traveling-wave components, considering left- and right-going cases separately at first, and arguing the general case by superposition.)

We can see from Eq.$ \,$(E.11) that such asymmetry can be caused by unequal weighting of $ y_{n,m}$ and $ y_{n,m\pm1}$. For example, the initialization

y_{n-1,m+1} &=& +1\\
y_{n-1,m} &=& -1

corresponds to an impulse velocity excitation at position $ m+1/2$. In this case, both interleaved grids are excited.

Next Section:
More General Velocity Excitations
Previous Section:
Localized Displacement Excitations