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

Help with linear PDE problem

  1. Sep 1, 2009 #1
    Here is the problem:
    ∂v(s,n)/∂n + ∂u(s,n)/∂s + ∂ξ(s,n)/∂s + An dc(s)/ds = 0 (1)
    A1 ∂ξ(s,n)/∂n + ∂v(s,n)/∂s -c(s)+A2 v(s,n) + A3 c(s) = 0 (2)
    ∂u(s,n)/∂s + 2A2 u(s,n)=A2(ξ(s,n) + Anc(s)) -A1 ∂ξ(s,n)/∂s-A2nc(s) (3)

    Unknowns: u(s,n),v(s,n),ξ(s,n)
    Boundary conditions: v(s,1)=v(s,-1)= 0

    I am trying to solve this set of PDEs by iteration and I am not sure if I am going about it the correct way. I have attached my attempt at a solution and it seems i am going in circles. Does anyone have a better idea? Someone suggested using Fourier analysis to solve the problem. I am reading that up now but i would really appreciate any ideas on how to start.

    Thank you.

    Attached Files:

  2. jcsd
  3. Sep 10, 2009 #2
    Dear tsetty2000,

    Your PDE system (if I do not misprint it) can be solved exactly by Maple 13 in the following way

    sys:=[diff(v(s,n),n) + diff(u(s,n),s) +diff(xi(s,n),s) + A*n*diff(c(s),s) = 0,
    A1*diff(xi(s,n),n) + diff(v(s,n),s) -c(s)+A2*v(s,n) + A3* c(s) = 0,
    diff(u(s,n),s) + 2*A2* u(s,n)=A2*(xi(s,n) + A*n*c(s)) -A1*diff(xi(s,n),s)-A2*n*c(s)];

    > ans:=pdsolve(sys,[xi,{u,v}]);

    > pdetest(ans,sys);

    The answer is bulky enough to be reproduced here (see att in Maple).

    Attached Files:

Share this great discussion with others via Reddit, Google+, Twitter, or Facebook