Dear list, I'm trying to solve the density driven natural convection due to the dissolution of CO2 in water in an aquifer. Please see the attached script. The script works fine for a Cartesian grid. However, it does not work on a cylindrical grid. I installed fipy from the 'develop' branch, or at least I think I did, and the problem still exist. My other issue is that all my boundary conditions for the continuity (pressure) equation are Neumann, which makes the Matrix of coefficient ill conditioned. Is there a way to fix the value of pressure (variable p) at a point? How can I check the condition number of the matrix of coefficients?
With kind regards, Ehsan P.S. The physics of the problem can be seen as a water saturated porous medium which is connected to a CO2 layer on top. CO2 is dissolved in water, increase its density and makes the system unstable due to the presence of gravity. We use the Boussinesq approximation, i.e., consider the effect of density only in the gravity term.
case05_v3.py
Description: Binary data
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
