0
$\begingroup$

I'm trying to use the Fourier transform method to solve the following PDE:

Original Question

This is a an infinite string with a pulse for it's initial condition. (At $t=0$, the string is stricken sharply so that the impulse $\Delta J$ is transferred instantaneously at position $x=0$).

I wrote the PDE in terms of the Fourier Transform:

Work 1

And then did algebra to get it here:

Fourier Transform U

My issue comes with doing the inverse Fourier transform, as I am running into 2 integrals I don't know how to solve.

First Integral

Piecing together my professors notes, somehow this is the solution to the first integral?

Solution(1) to first integral

But even then, I'm left with a second integral that has no solution that makes sense:

Second Integral

If I run this through Mathematica, I get a Sign(c) Sign(t) function, which has no $x$ dependence. Not only is this not a function of $U(x, t)$, neither second derivative is also not equal to the delta function, so I don't believe this is the correct answer.

If anyone knows how to approach these integrals, I would greatly appreciate it.

Edit: I think something was up with my assumptions in Mathematica, because now I am getting the following:

Solution?

Which looks more like a wave equation? Issue is, I'm not sure what the second* derivative for $\text{Sign}(c \cdot t \pm x)$ looks like with respect to either $t$ or $x$, so I am not able to confirm if this is correct.

$\endgroup$
3
  • $\begingroup$ The derivative of sign function is a Dirac-delta function: $\frac{d \mathrm{sgn}(x)}{d x} = 2 \delta(x)$ $\endgroup$
    – Galilean
    Commented Jan 17 at 6:35
  • $\begingroup$ Yes, but I need the second derivative of that. And as far as I can find, the derivative derivative of the delta function / second derivative of the Sign() function isn't very well defined. $\endgroup$
    – MoreDust
    Commented Jan 17 at 7:04
  • 2
    $\begingroup$ Use MathJax, not images. $\endgroup$ Commented Jan 17 at 8:57

2 Answers 2

2
$\begingroup$

Hints:

  1. Instead of Fourier transformation, it is easier to introduce light-cone coordinates $$x^{\pm}~=~x\pm ct,\tag{A}$$ which factorize OP's PDE: $$ \frac{\partial^2u}{\partial x^+\partial x^-}~=~-\frac{\Delta J}{2c\rho} \delta(x^+)\delta(x^-). \tag{B}$$

  2. The complete solution can now be easily obtained: $$ u~=~-\frac{\Delta J}{8c\rho}{\rm sgn}(x^+){\rm sgn}(x^-)+f_L(x^+) + f_R(x^-), \tag{C}$$ where $f_L,f_R$ are arbitrary functions.

  3. Next invoke various boundary conditions. Presumably OP is interested in the causal/retarded solution $$ u~=~\frac{\Delta J}{2c\rho}\theta(x^+)\theta(-x^-), \tag{D}$$ which vanishes for $t<0$.

$\endgroup$
2
$\begingroup$

Note that any function that is a linear combination of terms of the form $f(x \pm ct)$ satisfies the homogeneous wave equation (with no source term). Your proposed solution meets this criterion, so it cannot be a solution to the inhomogeneous equation with the impulse source term.

The solution Mathematica proposes is almost correct, we just have to "turn it on" at $t = 0$, so that $U = 0$ for $t < 0$. The correct solution can be written as $$ U = \frac1{4} \frac{\Delta J}{\rho c} u(t)\left[\mathrm{sgn}(ct - x) + \mathrm{sgn}(ct + x)\right]$$ where $u(t) = \frac12[1+\mathrm{sgn}(t)]$ is the unit (or Heaviside) step function.

To verify this solution, plugging it into the left-hand side of the wave equation gives $$(\partial_t^2-c^2\partial_x^2)U=\frac{\Delta J}{\rho}\left\{2\delta(t)\delta(x)+\frac1{4c}\delta'(t)[\mathrm{sgn}(ct-x)+\mathrm{sgn}(ct+x)\right\}.$$ Using the identity $\delta'(y)f(y) = f(0)\delta'(y) - f'(0)\delta(y)$, the second term becomes $-\delta(t)\delta(x)$, so we end up with the correct source term $\frac{\Delta J}{\rho}\delta(t)\delta(x).$

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.