How to use ode23s to solve heat equation
5 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Xiaoyan
am 13 Nov. 2014
Kommentiert: Torsten
am 24 Nov. 2014
assume ut=D*uxx D is constant, x = [-5:0.1:5], with delta_x=0.1 initial condition u(x0,0) = derivative of x boundary condition if x=-5 or 5 u(x,t)=0
hence can i write ODE that satisfy u1 to u99? then use ODE23s solve heat equation with initial data u50=1 and uj=0 where j=/50?
thank you!!
0 Kommentare
Akzeptierte Antwort
Torsten
am 21 Nov. 2014
If you are still interested in a solution for your problem:
Choose
u(x)=1/eps*max(1-abs(x/eps),0)
for eps=0.001, e.g. as initial condition for your Problem (this is an approximation to the Dirac Delta function).
Choose a very fine grid around x=0 such that the initial guess function is sufficiently resolved.
Then write the usual ODEs that result from
ut=D*uxx
and solve them using ODE15s.
If needed, I could also give you the analytical solution for comparison.
Best wishes
Torsten.
2 Kommentare
Torsten
am 24 Nov. 2014
The analytical solution for the problem
ut=D*uxx
u(x,0)=delta(x)
u(-5,t)=u(5,t)=0
is given by
u(x,t)=sum_{n=-oo}^{n=+oo}1/sqrt(Pi*D*t)*(exp(-(x-20*n)^2/(4*D*t))-exp(-(x-20*n+10)^2/(4*D*t)))
Best wishes
Torsten.
Weitere Antworten (4)
Orion
am 13 Nov. 2014
Xiaoyan,
You're problem is a partial differential equation (PDE), so ode23 won't work, it's meant to solves ordinary differential equation (ODE).
If you don't have it, then you're gonna have to code your own finite element method (such as finite difference). you could probably some opensource code on the web for those methods.
1 Kommentar
Torsten
am 13 Nov. 2014
Xiaoyan wants to apply the method of lines to the heat equation and solve the resulting system of ODEs using ODE23.
But I don't understand which Initial and boundary conditions he/she wants to impose.
Best wishes
Torsten.
Torsten
am 13 Nov. 2014
I guess numerical methods are not adequate if Dirac-Delta function is used as initial condition.
You should look up the analytical solution for this problem.
Best wishes
Torsten.
0 Kommentare
Zoltán Csáti
am 13 Nov. 2014
I would suggest an other way. A time-dependent differential equation can be solved not only by the method of lines, but also using the so-called Rothe's approach, ie. you discretize first in time and then is space. So use for example the forward Euler-method to prepare the semi-discrete equation that can be discretized by FDM/FEM or a spectral method.
0 Kommentare
Siehe auch
Kategorien
Mehr zu Ordinary Differential Equations finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!