Is It Possible to Solve This Diffusion Equation via Separation of Variables?

Remixex
Messages
57
Reaction score
4

Homework Statement


$$\frac{\partial U}{\partial t}=\nu \frac{\partial^{2} U}{\partial y^{2}}$$
$$U(0,t)=U_0 \quad for \quad t>0$$
$$U(y,0)=0 \quad for \quad y>0$$
$$U(y,t) \rightarrow {0} \quad \forall t \quad and \quad y \rightarrow \infty$$

Homework Equations


This is a diffusion problem on fluid mechanics, but it's more of a math problem so i posted it here.

The Attempt at a Solution


I'm trying to solve this via separation of variables (the textbook uses a "similarity" method I've never seen before, and concludes the function U must be erf) is it even possible to reach an analytic result via SV?
The first boundary condition is what gets me, I tried
$$U_{0} e^{{k^{2}t}} e^{{-\frac{k}{\sqrt{\nu}}y}}$$
But it clearly doesn't work for any boundary condition except the last.
I don't think sinusoidal is the answer here either because it must eventually converge to zero, for every t.
Is there really no analytic answer?
 
Physics news on Phys.org
Please show your separation of variables and the resulting differential equations.
 
Remixex said:

Homework Statement


$$\frac{\partial U}{\partial t}=\nu \frac{\partial^{2} U}{\partial y^{2}}$$
$$U(0,t)=U_0 \quad for \quad t>0$$
$$U(y,0)=0 \quad for \quad y>0$$
$$U(y,t) \rightarrow {0} \quad \forall t \quad and \quad y \rightarrow \infty$$

Homework Equations


This is a diffusion problem on fluid mechanics, but it's more of a math problem so i posted it here.

The Attempt at a Solution


I'm trying to solve this via separation of variables (the textbook uses a "similarity" method I've never seen before, and concludes the function U must be erf) is it even possible to reach an analytic result via SV?
The first boundary condition is what gets me, I tried
$$U_{0} e^{{k^{2}t}} e^{{-\frac{k}{\sqrt{\nu}}y}}$$
But it clearly doesn't work for any boundary condition except the last.
I don't think sinusoidal is the answer here either because it must eventually converge to zero, for every t.
Is there really no analytic answer?

Use Laplace transforms with respect to ##t##. Let
$$W(y,s) = \int_0^{\infty} e^{-st} U(y,t) \, dt $$
be the Laplace transform. Then, using standard properties of Laplace transforms, we get the DE
$$\nu W_{yy}(y,s) = s W(y,s) - U(y,0) = s W(y,s),$$
where ##W_{yy} = \partial^2 W / \partial y^2##.
Also: ##U(0,t) = U_0## implies that
$$W(0,s) = \frac{U_0}{s} $$
Finally, the initial value theorem requires that ##\lim_{s \to \infty} s W(y,s) = 0## for ##y > 0##.

These are enough to determine ##W(y,s)##. Then it is just a matter of taking the inverse Laplace transform of ##W(y,s)## to get ##U(y,t)##.

Separation of variables will never work in this example, simply because it leads to the wrong kind of function.
 
Last edited:
Ray Vickson said:
Separation of variables will never work in this example, simply because it leads to the wrong kind of function.
THANK YOU :D
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top