How can I implement this PDE system ?
15 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Clément Douarre
il 15 Mar 2016
Modificato: Clément Douarre
il 22 Mar 2016
Hello all,
I am trying to implement a system of Partial Differential Equations. i am using the PDE Toolbox which seems quite complete.
The 2-D system is the following :
dn/dt + nabla * (nv) = 0 [1]
dv/dt + v * nabla v = nabla c [2]
dc/dt = D nabla^2 c + alpha n - c/tau [3]
Where the variables are n, v and c, and depend on x,y and t. The derivatives are of course partial ones. alpha, tau and D are parameters. Variable v is a vector.
This system comes from a scientific article, and the following indications are given :
[1] is a mass conservation equation
[2] is a Burgers equation
[3] is a reaction-diffusion equation
The main problem I have is that [3] is the only equation I managed to get in shape for the PDE toolbox to be able to solve.
Indeed, the form that the toolbox takes is : d * du/dt nabla·(nabla c u)+au=f, where d, c, a and f are coefficients which may depend on u, du/dx, du/dy, x and y.
I cannot manage to see how to write [1] and [2] in this form, especially for the second term.
This makes me think that perhaps I am considering the problem all wrong. These PDEs seem quite different to me, and maybe I'm not taking the right approach by trying to jumble them all together.
Any help on how to get started on implementing this would be a great help.
Thank you in advance, Clément
0 Commenti
Risposta accettata
Torsten
il 15 Mar 2016
I think a good starting point are the Euler equations of fluid dynamics.
In 1d, a well-written code to solve them is CLAWPACK available under
https://depts.washington.edu/clawpack/
I don't know of any MATLAB implementations to solve the Euler equations.
Good luck !
Best wishes
Torsten.
1 Commento
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Eigenvalue Problems in Help Center e File Exchange
Prodotti
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!