Mass balance issue when solving an elliptic diffusion problem using pde toolbox

1 次查看（过去 30 天）
Mohammad Monfared
Mohammad Monfared 2014-1-21
编辑： Mohammad Monfared 2014-1-30
Hi
I'm trying to solve a diffusion like, elliptic equation using pde toolbox. I used 'pednonlin' function but the answer I get, gets me a bit confused regarding mass conservation issue since the solution should be steady-state.
Here are more about problem:
The equation is simply -div(K(u)*del(u)) = 0 on a rectangular domain with Dirichlet B.C. on top and bottom and no flux B.C. on side boundaries. But when I plot the solution using:
[cux,cuz] = pdecgrad(p,t,c,u) ;
pdeplot(p,e,t,'flowdata',[-cux;-cuz])
the flux has a pattern which induces mass is not conserved!
Does anyone have an idea where I'm going wrong?
Thanks a lot.

请先登录，再回答此问题。

采纳的回答

Mohammad Monfared
Mohammad Monfared 2014-1-30
编辑：Mohammad Monfared 2014-1-30
The issue was that the Jacobian has to be calculated using 'full' parameter, not 'fixed' or 'lump'.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Geometry and Mesh 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by