1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Solution of unsteady linearized potential flow PDE

  1. Oct 1, 2015 #1

    I have a problem following the solution of a linearized potential flow equation in a publication by Fung.

    The problem describes potential flow over an oscillating plate. A boundary layer is approximated by defining a subsonic layer over the panel and supersonic flow above the subsonic flow. From the equation of motion (1) and (2) in combination with a standing wave condition of the wall (8) and traveling wave of the perturbations (9) and (10) it seems to be easy to get the solutions (12) and (13).

    https://dl.dropboxusercontent.com/u/20358584/fung1.png [Broken]
    https://dl.dropboxusercontent.com/u/20358584/fung2.png [Broken]

    Can anyobody give me hint of how to get to this solution? The paper is the following:
    [/PLAIN] [Broken]

    Many thanks in advance!!!

    https://dl.dropboxusercontent.com/u/20358584/fung3.png [Broken]
    Last edited by a moderator: May 7, 2017
  2. jcsd
  3. Oct 1, 2015 #2

    Andy Resnick

    User Avatar
    Science Advisor
    Education Advisor

    When I am confronted with "it's easy to see that..." I usually first try substituting the answer into the expression and seeing what happens- sometimes there's an oddball change of variables or trig identity involved.
  4. Oct 29, 2015 #3
    I don't know. This doesn't help. What could he have done? I've heard that Duhamel's principle could be an approach of solving non-homogenious PDEs like the wave equation. Could the solution have something to do with this approach?

    Substituting (13) in (2) gives:
    [itex] \begin{equation} [-\frac{1}{a_{\delta}^2} \omega^2 - \frac{2M_{\delta}}{a_{\delta}} \alpha_{\nu}\omega + \beta_{\delta}^2\alpha_{\nu}^2 + \zeta_{\delta}^2] \cdot e^{i(\omega t + \alpha_{\nu} x)} \cdot [C_{\nu} sin(\zeta_{\nu}y ) + D_{\nu} cos(\zeta_{\nu}y) ] = 0 \end{equation} [/itex]
    Last edited: Oct 29, 2015
  5. Nov 2, 2015 #4
    Oh... the substitution was absolutely wrong. I need to find the correlation between the potential and z..
  6. Nov 4, 2015 #5
    I made some progress to get equation (14) and (15). Not sure if I can already corellate the fourer constant alpha with the wave number.

    The equation of motion:
    \frac{1}{a^2}\frac{\partial^2\phi}{\partial t^2} + \frac{2M}{a} \frac{\partial^2 \phi}{\partial x \partial t} + \overline{\beta}^2 \frac{\partial^2 \phi}{\partial x^2} = \frac{\partial^2 \phi}{\partial y^2}
    The potential must oscillate harmonically:
    \phi = \Psi(x,y)e^{i\omega t}
    This yields:
    \Big(\frac{i\omega}{a}\Big)^2\Psi + 2i\frac{\omega M}{a} \frac{\partial \Psi}{\partial x} + \overline{\beta}^2 \frac{\partial^2 \Psi}{\partial x^2} = \frac{\partial^2 \Psi}{\partial y^2}
    A double Fourier transformation in x and y:
    \Psi^* = \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} e^{-i(\gamma y + \alpha x)} \Psi(x,y) dx dy
    If this Fourier transformation is applied to all terms of Eq.~(\ref{eq:03}) then $\Phi^*$ cancels out and (\ref{eq:03}) can be written as:
    \frac{\omega^2}{a^2} + 2 \frac{\omega M}{a} \alpha + \overline{\beta}^2 \alpha^2 = \gamma^2
    This is equation (14) in the publication from Fung. The same approach for Fungs equation (2) yields (15). Can I already assume that $\gamma$ is $\gamma_{\nu}$ and $\alpha$ is $\alpha_{\nu}$?
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook