How can I implement this PDE system ?
14 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Clément Douarre
le 15 Mar 2016
Modifié(e) : Clément Douarre
le 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 commentaires
Réponse acceptée
Torsten
le 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 commentaire
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Eigenvalue Problems dans Help Center et File Exchange
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!