How can I implement this PDE system ?

조회 수: 14 (최근 30일)
Clément Douarre
Clément Douarre 2016년 3월 15일
편집: Clément Douarre 2016년 3월 22일
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

채택된 답변

Torsten
Torsten 2016년 3월 15일
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
Clément Douarre
Clément Douarre 2016년 3월 22일
편집: Clément Douarre 2016년 3월 22일
Thank you for your quick answer (and sorry for my late acceptance). I have accordingly looked away from Matlab, at Clawpack, but have found that FiPy seems more complete. Anyway, thanks for confirming that it would be very difficult -if possible at all - to implement in Matlab.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Eigenvalue Problems에 대해 자세히 알아보기

제품

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by