CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[waves2Foam] VRANS vs OpenFOAM Sink Term

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 3, 2023, 13:19
Default VRANS vs OpenFOAM Sink Term
  #1
Member
 
Felix S.
Join Date: Feb 2021
Location: Germany, Braunschweig
Posts: 86
Rep Power: 6
Felix_Sp is on a distinguished road
Hello everyone,

I guess this is not a specific thread only in regard to waves2Foam, but also in regard to OlaFLOW or any other solver using VRANS.
One thing that's been bugging me for the past few weeks and I can't understand is the difference between VRANS and specifying a local sink term for the momentum equation, such as OpenFOAM's explicitPorositySource. Since most models use the Darcy-Forchheimer equation, I will assume the same for the sink terms for my questions below. I think it shouldn't matter how you quantify the pressure drop, as long as you quantify it .

Now when there is a pressure drop, e.g. a liquid flowing through a porous medium in a pipe, then one can often quantify the pressure drop by a polynomial of some order. Often second order is enough, fulfilling a version of the Darcy-Forchheimer equation. Comparing the coefficients of the polynomial and the Darcy-Forchheimer equation, one can calculate the laminar and turbulent pressure drop and can then apply the sink term to model the porous medium.

Now, why do we need to include the porosity in the momentum equation leading to the VRANS equations. Is the sink term a problem if one wants to model the actual flow inside the porous medium, as the acceleration is wrong otherwise? What problem arises by modelling the pressure drop inside a porous medium by just a sink term. Because now I cannot simply use the polynomial equation as basis for the sink term, as the porosity in the momentum equation is changing the acceleration of the fluid inside the porous zone.

If anyone has a good read comparing both methods, I'd be happy to look around the literature, as I haven't found anything along those lines.

Thanks for your help in advance
Felix
Felix_Sp is offline   Reply With Quote

Old   August 3, 2023, 13:52
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,901
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Felix,


It is simply the way the momentum equation comes out, when you do the correct averaging over the pore space (see e.g. https://www.sciencedirect.com/scienc...78383913001816). Furthermore, for the free surface models, you need to account for the solids in the computational cells, as your free surface does not move fast enough otherwise (see again Jensen et al. above and https://journal.openfoam.com/index.p...rticle/view/72).


Kind regards


Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   August 3, 2023, 14:07
Default
  #3
Member
 
Felix S.
Join Date: Feb 2021
Location: Germany, Braunschweig
Posts: 86
Rep Power: 6
Felix_Sp is on a distinguished road
Hey Niels,

thank you for your reply. I guess I understand, that it is the correct mathematical derivation. However, if this is the case, when can you apply the simple source term integrated into OpenFOAM (https://www.openfoam.com/documentati...orosity.html)? Or would such an application always be non-physical, as one only models the pressure drop and does not account for the blocking solids inside the porous medium?

Thanks again for your help!
Felix
Felix_Sp is offline   Reply With Quote

Old   August 5, 2023, 05:43
Default
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,901
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Felix,


I would say that it is formally non-physical, but it is for you to assess, whether the bias is acceptable for your application. Personally, I would always apply the porousWaveFoam in waves2Foam.


Kind regards


Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Reply

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Question about porous term and porosity value in OpenFOAM sazabi2001 OpenFOAM Running, Solving & CFD 12 May 17, 2020 20:57
problems creating a sink term with scalarCodedSource Stephen Waite OpenFOAM Programming & Development 4 June 7, 2018 06:26
OpenFOAM Training Jan-Apr 2017, Virtual, London, Houston, Berlin cfd.direct OpenFOAM Announcements from Other Sources 0 September 21, 2016 11:50
OpenFOAM Training, London, Chicago, Munich, Houston 2016-2017 cfd.direct OpenFOAM Announcements from Other Sources 0 September 14, 2016 03:19
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 3 September 8, 2010 06:25


All times are GMT -4. The time now is 13:08.