System of PDEs which is tricky for PDEPE
14 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
I have a system of PDEs, mainly diffusion equations of the form:
T_t-(k(x)(T_x)_x=a*E^2
c_t-(D(x)*c_x)_x=d*(E_x+c_x-T_x)
(epsilon*E)_x=-b*c
Where a and d are constants and _t,_x represent partial differentiation w.r.t. t and x respectively. In terms of pdepe, I would have: f=0 for the third equation. Would this cause errors in the code?
I'm also struggling to see how I can input the boundary conditions.
Any suggestions?
5 comentarios
Bill Greene
el 24 de Mayo de 2018
Do you mean T=0 at the boundaries? What about the initial conditions on the three variables? At t=0, the boundary and initial conditions must be consistent and also satisfy your differential equations.
Respuestas (5)
Torsten
el 23 de Mayo de 2018
Editada: Torsten
el 23 de Mayo de 2018
"pdepe" is designed to solve systems of parabolic-elliptic pdes. Your third pde is hyperbolic in nature. Thus "pdepe" is not suited to solve your system.
You will have to discretize your equations in space and solve the resulting system of ordinary differential equations using ODE15S.
Look up "method-of-lines" for more details.
Best wishes
Torsten.
2 comentarios
Torsten
el 23 de Mayo de 2018
Since "pdefun" is called only for a single value of x, you don't have the complete vector for c available. Thus no chance to calculate E, I guess.
Torsten
el 24 de Mayo de 2018
You could try to use "pdepe" with the third equation differentiated:
(epsilon*E)_xx=-b*c_x
with Dirichlet boundary condition
E = E0
at one end of the interval and
(epsilon*E)_x + b*c=0
at the opposite end.
Best wishes
Torsten.
3 comentarios
Torsten
el 24 de Mayo de 2018
But the E_x and T_x terms in the diffusion equation for c will be put in the source term s for this equation, not in the f-term.
Precise Simulation
el 9 de Jun. de 2018
Systems of PDEs might also be easier to solve with the FEATool FEM Toolbox which features a GUI and easy syntax for defining custom PDEs and equations. This example of heat transport and diffusion might be a good start, otherwise you could perhaps use the built-in convection and diffusion physics mode and modify it according to your equations.
0 comentarios
Ver también
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!