CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Identification of the porous region in powerLawLopesdaCosta porosity model

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 7, 2020, 02:19
Default Identification of the porous region in powerLawLopesdaCosta porosity model
Senior Member
Join Date: Aug 2010
Location: State College, PA
Posts: 111
Rep Power: 11
Tarak is on a distinguished road

The standard porosity models implemented in OpenFOAM, e.g DarcyForchheimer and powerLaw, identify the regions in the computational domain where the body force term has to be added using the cellZone supplied by the user. This offers a very convenient way to prescribe the porous regions in the domain.

But, in the newly added "powerLawLopesdaCosta porosity model", instead of utilizing user-prescribed cellZones, the model explicitly identifies the cells in the porous region using a convoluted process (using information from surface files representing the top and bottom regions of the porous region, and generating cellZone while running the model). Maybe, there's an option to not create a new cellZone if there's already an existing one (created by the user manually, corresponding to the else branch of if(zoneID == -1), but it still requires surface files corresponding to the top and bottom of the porous regions to run the model.

I'm not sure why this length process has been pursued for the powerLawLopesdaCosta porosity model, is it only to allow the porosity surface area per unit volume parameter (Sigma) to be a function of height?
Tarak is offline   Reply With Quote

Old   February 11, 2020, 17:32
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Posts: 298
Rep Power: 3
HPE is on a distinguished road
i think the term 'porosity' for that model is very unfortunate. that term should be 'canopy' instead, since that particular model is for modelling vegetation canopy of atmospheric boundary layers. I dont think therefore this model relates to what other porosity models do. i might be wrong as well, but quite unlikely

hope it helps.
Kind regards,

The OpenFOAM community is the biggest contributor to OpenFOAM; not person A or B. Let's report bugs or issue feature requests in GitLab-OpenFOAM.

Please `design` your questions well. Such effort always pays off well for inquirers and `helpmate`s:
How to create a Minimal, Reproducible Example,
How to give enough info to get help.
HPE is online now   Reply With Quote


porosity model, powerlawlopesdacosta

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
Thermal non-equilibrium porous media model with conjugate heat transfer Hexahedron FLUENT 4 January 26, 2020 11:46
Question about porous term and porosity value in OpenFOAM sazabi2001 OpenFOAM Running, Solving & CFD 11 May 3, 2019 02:49
chtMultiRegionSimpleFoam with porous solid region. srikanthallu OpenFOAM 14 September 14, 2018 12:23
Can I use fvOptions to couple a solid region and a fluid region? titanchao OpenFOAM Running, Solving & CFD 1 March 7, 2018 08:46
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 08:00

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