Mass balance issue when solving an elliptic diffusion problem using pde toolbox
1 回表示 (過去 30 日間)
古いコメントを表示
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.
0 件のコメント
採用された回答
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Eigenvalue Problems についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!