Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Causality and the inhomogeneous wave equation with a moving source

  1. Jan 18, 2012 #1

    I am confused about my solutions to the following governing equation:

    [tex] u_{tt}-c^2u_{xx}=F_{xx} [/tex]

    For [tex]F=A(x)sech^2\left(\frac{x-c_gt}{B}\right) [/tex]

    Where [itex]c,c_g,B \in \mathbb{R} [/itex] and [itex] A(x)[/itex] is a linear function. Also, we have [itex]c_g<c[/itex]. Substituting physical values for the parameters, I can solve the equation using a 2-d integrator to implement the method of Duhamel. My solutions are unclear to me though.

    I expect to see the intensity of the response increase with position, as it does along the ray traveling at [itex]c_g[/itex]. I do not, however, see the same response out of the ray traveling at [itex]c[/itex]. I also do not see a ray traveling at [itex]-c[/itex], as I naively would expect. In the plot attached, the ray traveling at speed c is the lower one, since this is an x-t diagram.

    I also would think that the ray traveling at speed c would constantly be generated at the location of the forcing, otherwise it is unclear to me how the information from the forcing would get to the faster wave traveling at speed c.

    I think my numerical implementation is fine. Perhaps the fact that I'm working on the half domain, [itex] x \ge 0 [/itex] is messing things up?

    Any help is appreciated,


    Attached Files:

  2. jcsd
  3. Jan 18, 2012 #2
    Unfortunately, I can not grasp the physical sense of your problem out of hand from your picture, but it seems to me that you may analyse the analytical solution of your problem.

    For linear [itex]A(x)=ax+b[/itex] the general solution to your PDE is relatively easy:

    [itex]u(t,x) = C_1(x+ct)+C_2(x-ct)+4\{[-a(c^2-c_g^2)x-bc^2+c_g^2(aB+b)]\exp[-2(x-c_gt)/B]+aBc_g^2\}\{exp[-2(x-c_gt)/B]+1\}^{-2}(c^2-c_g^2)^{-2},[/itex]

    where [itex]C_1,C_2[/itex] are arbitrary functions.

    I hope that the like solution can be found for more complicated [itex]A(x)[/itex].

    Then you have to impose the initial conditions of your physical problem to obtain particular solution which then can be analysed under your goals.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook