Mass balance issue when solving an elliptic diffusion problem using pde toolbox
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
Mohammad Monfared
il 21 Gen 2014
Modificato: Mohammad Monfared
il 30 Gen 2014
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!
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/155389/image.jpeg)
Does anyone have an idea where I'm going wrong?
Thanks a lot.
0 Commenti
Risposta accettata
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Eigenvalue Problems in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!