Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.
Hydrodynamic dispersion in porous media
Posted May 19, 2019, 4:03 AM EDT Fluid, Computational Fluid Dynamics (CFD), Subsurface Flow 0 Replies
Please login with a confirmed email address before reporting spam
I'm currently working on a mass induced natural convection 2D problem in porous media (single phase miscible displacement).
At the initial condition, component 1 fills a rectangular domain. On top, a Dirichlet BC is set with mass fraction(component 2)=1. A non-monotonic density-mass fraction function is implemented, such that when component 1 is pure, it is lighter than component 2. However, when they mix due to diffusion, the mixture becomes heavier forming instabilities which lead to the classical fingering patterns.
In order to do so i'm solving the two components mass conservation. The conservation of momentum in done via the Darcy law for porous media. Hydrodynamic dispersion is represented by the scheidegger model which states linear dependency between the mechanical dispersion and velocity (multiplied by the dispersivity coefficient). Finally, a linear viscosity-mass fraction function is implemented.
So far so good. I currently use Comsol (finite element and implicit scheme) to do so, and manage to solve the system to get reasonable results (compared to experimental results).
Problem starts when i increase the dispersivity coefficient. For some reason, which is counter intuitive to me, when i increase dispersivity, for the same element size, convergence becomes worse up to a point where it doesn't converge at all. I assumed that increasing dispersivity would play a stabilizing factor (decreasing the elemental Peclet number), but for some reason it is not.
Hello Avihai Tsinober
Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.
If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.