CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

darcy and continuity equations coupling

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree3Likes
  • 1 Post By mprinkey
  • 1 Post By mprinkey
  • 1 Post By mprinkey

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 6, 2015, 18:37
Default darcy and continuity equations coupling
  #1
New Member
 
Morocco
Join Date: Mar 2015
Posts: 11
Rep Power: 11
zaynab is on a distinguished road
Hi everyone,

I am working on a finite volume modeling, the system is a porous media, composed of a solid phase and a fluid phase.
unstead of the momentum equation, the velocity of the fluid is described by the darcy equation, where the velocity is directly related to the pressure gradient.
my question is:
should i apply the same methods (simple, simpler,piso..) algorithms for the coupling of darcy's law and the continuity equation?
zaynab is offline   Reply With Quote

Old   April 7, 2015, 09:05
Default
  #2
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
I think everything is much easier than general pressure/velocity coupling with Navier-Stokes. The face mass/volume fluxes arise from the face gradient of the pressure--that makes it a function of the two cells on either side of face (with, perhaps skew terms added in explicitly). Then substitute those face fluxes in the continuity equation. That gives you a Poisson(-like) equation for pressure. Solve it, and explicitly update the velocities at the faces using the grad(p). There is your Darcy flow.

The SIMPLEx/PISO techniques are algorithms attempting to converge the effects of inertia and convective transport manifested in the momentum equations with the incompressibility constraint in the continuity equation. In Darcy flow, the momentum equation reduces so significantly that it no long includes the competing effects. Pressure dictates velocity directly.
zaynab likes this.
mprinkey is offline   Reply With Quote

Old   April 7, 2015, 18:54
Default
  #3
New Member
 
Morocco
Join Date: Mar 2015
Posts: 11
Rep Power: 11
zaynab is on a distinguished road
Thank you, dear Michael, for your detailed answer

So, there will be no checkboard pressure problem if I use the pressure values of the discretised poisson's equation (1D)(d²p/dx²= 0) to calculate the velocity values using the discretised form of the darcy equation?
zaynab is offline   Reply With Quote

Old   April 7, 2015, 19:10
Default
  #4
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
No, because you are not solving a linear system for velocity as you do with SIMPLEx/PISO applied to NS. Making the substitution I outlined leaves you with only one independent variable--the pressure. There is nothing left to get out of sync. The velocity field (at the faces) uniquely follows from the pressure field (on the cells).
zaynab likes this.
mprinkey is offline   Reply With Quote

Old   April 7, 2015, 19:15
Default
  #5
New Member
 
Morocco
Join Date: Mar 2015
Posts: 11
Rep Power: 11
zaynab is on a distinguished road
Ok, Thank you so much.
zaynab is offline   Reply With Quote

Old   April 7, 2015, 19:31
Default
  #6
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
I should point out a few things. There may be a gravity term in the Darcy law. That could give you something more complicated than a Laplace equation for pressure. The other issue is with the wildly varying permeability from cell to cell that can occur in real simulations. You will need to do something like harmonic averaging to get that formulated correctly. Good luck.
zaynab likes this.
mprinkey is offline   Reply With Quote

Reply

Tags
continuity, darcy, porous media, simple algorithm

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 06:18.