Re: Solving Parabolic non-linear PDE using MATLAB


In article <1174479720.669669.118750@xxxxxxxxxxxxxxxxxxxxxxxxxxxx>,
"Grazia" <dalgranati@xxxxxxxxxxx> writes:
Hi guys
I need to solve a one-dimensional non linear parabolic PDE, using
MATLAB
PDE looks like this:
d(U^2)/dt = f(t)*d2(U^6)/dx2

I tried using pdepe.m function. However, this function requires a
certain format. Specifically:

C(x,t,u,du/dx)d(u)/dt = d[F(x,t,u,du/dx)]/dx+S(x,t,u,du/dx). Does
anyone have any idea how I translate my function to fit this format?

i.e. What are the function C,S and F that I need to input pdepe.m?

Thanks,

Grazia


d(U^2)/dt = 2*U*(dU/dt) => C(....)=2*U

d[F(x,t,u,du/dx)]/dx =
F_x + F_u*(du/dx) + F_{u_x} u_xx
and
f(t) d2(U^6)/dx2 = f(t)* d(6*U^5*U_x)/dx
=> F(t,x,u,u_x) = f(t)*6*U^5*U_x
F(t,x,u,v) = f(t)*6*u^5*v
hth
peter
.

Relevant Pages