Explicitly Solving a PDE: Is There a Solution?

  • Thread starter Thread starter Clever-Name
  • Start date Start date
  • Tags Tags
    Pde
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 3K views
Clever-Name
Messages
378
Reaction score
1
I just had a 'quiz' in my PDE class today and there was a problem my friends and I are convinced has no explicit solution. I want to know if maybe we are doing something wrong?

Homework Statement


[tex](x+y)u_{x} + yu_{y} = 0 [/itex]<br /> <br /> [tex]u(1,y) = \frac{1}{y} + ln(y) [/itex]<h2>Homework Equations</h2><br /> ...<h2>The Attempt at a Solution</h2><br /> Here was my approach via method of characteristics:<br /> <br /> [tex]\frac{dy}{dt} = y \to y(t) = Se^{t}[/tex]<br /> [tex]y(t=0) = S[/tex]<br /> <br /> [tex]\frac{dx}{dt} = x+y[/tex]<br /> [tex]x'-x = Se^{t}[/tex]<br /> Via integrating factor we arrive at:<br /> [tex]x(t) = Ste^{t} + g(s)e^{t}[/tex]<br /> [tex]x(t=0) = 1 \to 1 = g(s)[/tex]<br /> <br /> Thus we are left with:<br /> <br /> [tex]x(t) = Ste^{t} + e^{t}[/tex]<br /> and<br /> [tex]y(t) = Se^{t}[/tex]<br /> <br /> You can check for yourself that the following is true:<br /> <br /> [tex]\frac{dU(x(t),y(t))}{dt} = (x+y)u_{x} + yu_{y} = 0[/tex]<br /> <br /> Thus our general solution is:<br /> <br /> [tex]U(t) = f(s)[/tex]<br /> <br /> A general function of s.<br /> <br /> However, we have found that it is impossible to explicitly solve for s. Am I missing something?[/tex][/tex]
 
Physics news on Phys.org
I suspect the problem cannot be solved. Along the initial line [itex]y=1[/itex] you can differentiate the initial conditions and get

[itex]u_y(1,y)=-1/y^2+1/y=(y-1)/y^2=0[/itex]

that you can substitute back into the PDE in order to obtain [itex]u_x[/itex]

[itex](x+1)u_x=0 → u_x=0[/itex]

which means [itex]u(x,y)=u(y)[/itex]. Hence, [itex]u(x,y)[/itex] should be equal to the initial conditions everywhere. But this is not true, as one can check substituting
[itex]u_x=0[/itex] into the PDE:

[itex]yu_y=y(1-y)/y^2=0[/itex]

true only when [itex]y=0[/itex] or [itex]y=1[/itex].

I am not an expert of the method of characteristics, but I believe the curve on which the initial conditions are given is "noncharacteristic".