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

[waves2Foam] Understanding porousWaveFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 22, 2023, 16:15
Default Understanding porousWaveFoam
  #1
New Member
 
saeed barzegar
Join Date: Feb 2012
Posts: 19
Rep Power: 14
saeed.barzegar.v is on a distinguished road
Hi All,

I've spent a lot of time, and I'm at a point where I have no clue. I was wondering if someone could help me understand porousWaveFoam. In the code (which I'm guessing is built based on porousSimpleFoam), I see the term "pm->addResistance(UEqn)", which represents the resistance forces and the porosity term, n, is only included in the convention term. Based on waves2Foam's manual (eq. 2.4), Jensen et al. (2014), Higuera et al. (2014), Jacobsen et al., 2015, and many other references, shouldn't there be a porosity value behind the transient and viscous terms in UEqn.H file?


What we have in porousWaveFoam is:
fvVectorMatrix UEqn
(
pm->ddt(rho, U)
+ 1.0/porosity*fvm::div(rhoPhi/(porosityFace), U)
+ MRF.DDt(rho, U)
+ turbulence->divDevRhoReff(rho, U)
==
fvOptions(rho, U)
);


Shouldn't it be:
fvVectorMatrix UEqn
(
pm->(1.0 + cPorField) / porosity *ddt(rho, U)
+ 1.0/porosity*fvm::div(rhoPhi/(porosityFace), U)
+ (1.0 + cPorField) / porosity * MRF.DDt(rho, U)
+ 1.0/porosity * turbulence->divDevRhoReff(rho, U)
==
fvOptions(rho, U)
);


Where cPorField is Cm in those references. When comparing it to olaFlow (link below), I can't understand why the porosity representative parameters behind transient and viscous terms are missing.

https://github.com/phicau/olaFlow/bl..._latest/UEqn.H

Thank you,
Saeed
saeed.barzegar.v is offline   Reply With Quote

Old   November 23, 2023, 14:15
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Saeed,


Try to follow the "pm->" pointer. It will be a deep dive, where everything happens in the class jensenJacobsenChristensen2014. It is unfortunately a more messy corner of 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

Old   November 23, 2023, 16:26
Default
  #3
New Member
 
saeed barzegar
Join Date: Feb 2012
Posts: 19
Rep Power: 14
saeed.barzegar.v is on a distinguished road
Quote:
Originally Posted by ngj View Post
Hi Saeed,

Try to follow the "pm->" pointer. It will be a deep dive, where everything happens in the class jensenJacobsenChristensen2014. It is unfortunately a more messy corner of waves2Foam.

Kind regards
Niels
Thanks Niels! Appreciate your quick response. I'll dig more into it and try to understand the jensenJacobsenChristensen2014 class you mentioned. I'll keep you posted on my findings.

Cheers,
Saeed
saeed.barzegar.v 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
Understanding hardware & parallel processing requirements yk2359 Hardware 3 May 24, 2018 13:34
Simple parabolic velocity UDF with symmetry - what am I not understanding!?! NPGrad FLUENT 1 September 30, 2016 11:42
Static or Total Pressure Inlet Boundary, Understanding problem Overdue CFX 13 August 4, 2016 07:25
Understanding pimpleFoam convergence criterion Nucleophobe OpenFOAM Running, Solving & CFD 0 March 13, 2013 18:46
y+Factor: Help understanding? eRzBeNgEl STAR-CCM+ 2 January 26, 2012 16:23


All times are GMT -4. The time now is 17:30.