31.08.2018 Views

Eduardo Kausel-Fundamental solutions in elastodynamics_ a compendium-Cambridge University Press (2006)

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

102 Solution to the Helmholtz and wave equations<br />

and substitut<strong>in</strong>g <strong>in</strong>to the Helmholtz equation, we f<strong>in</strong>d that it is satisfied if<br />

k 2 x + k2 y + k2 z = k2 P<br />

(8.47)<br />

The general solution to the scalar Helmholtz equation is then<br />

(x, y, z,ω) = ( a 1 e ikxx + a 2 e −ikxx)( a 3 e iky y + a 4 e −iky y)( a 5 e ikzz + a 6 e −ikzz) (8.48)<br />

<strong>in</strong> which the a i are arbitrary constants. In comb<strong>in</strong>ation with the implied factor e iωt , and for<br />

real wavenumbers k x , k y , k z , this solution for represents compressional (P) waves that<br />

propagate <strong>in</strong> the (locally homogeneous) medium with velocity α.<br />

For two-dimensional problems <strong>in</strong> which the wave velocity α does not change <strong>in</strong> horizontal<br />

planes (i.e., for a laterally homogeneous medium), and consider<strong>in</strong>g only waves that<br />

propagate and/or decay <strong>in</strong> the positive x direction, the solution simplifies to<br />

(x, z,ω) = ( a 1 e ikzz + a 2 e −ikzz) √<br />

e −ikxx , k z = kP 2 − k2 x (8.49)<br />

8.3 Vector Helmholtz equation <strong>in</strong> Cartesian coord<strong>in</strong>ates<br />

Consider the vector Helmholtz equation<br />

∇ 2 Ψ + k 2 S Ψ = 0 (8.50)<br />

with parameter k s = ω/β. More precisely,<br />

∇ · ∇Ψ + k 2 S Ψ = 0 (8.51)<br />

which <strong>in</strong>volves the vector Ψ = ψ x î + ψ y ĵ + ψ z ˆk. We shall assume Ψ to be solenoidal, which<br />

means that it satisfies the gaugecondition ∇ · Ψ = 0. Hence, only two of the three components<br />

of Ψ are <strong>in</strong>dependent, so the solution can only conta<strong>in</strong> two <strong>in</strong>dependent functions.<br />

However,<br />

∇ · ∇Ψ =∇(∇ · Ψ) −∇×∇×Ψ<br />

Hence<br />

=−∇×∇×Ψ (8.52)<br />

∇×∇×Ψ = k 2 S Ψ (8.53)<br />

This equation admits <strong>solutions</strong> of the form<br />

Ψ = Ψ 1 + Ψ 2<br />

= ψ ˆk + 1 k 2 S<br />

∇<br />

( ) ∂ψ<br />

+ 1 ( )<br />

∇× χ ˆk<br />

∂z k S<br />

(8.54)<br />

<strong>in</strong> which ψ, χ are <strong>solutions</strong> to the scalar Helmholtz equations<br />

∇ 2 ψ + k 2 S ψ = 0 and ∇2 χ + k 2 S χ = 0 (8.55)<br />

To prove that this is <strong>in</strong>deed the solution, we make use of some well-known vector operation<br />

identities, namely ∇×∇f = 0, ∇ · ∇×f = 0, which are true for any scalar and vector

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!