CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Porous Modeling of Energy equation in OpenFOAM

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 28, 2015, 08:32
Default Porous Modeling of Energy equation in OpenFOAM
  #1
New Member
 
mohammad amin mohammadi
Join Date: Mar 2012
Location: Iran
Posts: 2
Rep Power: 0
mohammad_kordo is on a distinguished road
Hi to all In momentum equation effect of porous with considering volume averaged equation, is a source term named Darcy-Forchheimer term. It is considered as source term that added to momentum in UEqn.H file that can be seen in porousSimpleFoam or rhoPorousSimpleFoam solvers. But in energy equation, porous modeling is based on tow model. One is local thermal equilibrium and another non-thermal equilibrium. In local thermal equilibrium the temperature of solid (porous) and fluid is same and effect of porous properties (density, specific heat coefficient and thermal conductivity) is considered in fluid energy equation. Term of “rhof*cpf” (multiply of density in specific heat coefficient of fluid) in energy equation is modified to
rhof*cpf+rhos*cps*(1-epsilon)/epsilon : f:fluid, s:solid.
Also conductivity (k) is modified from kf to kf+ks*(1-epsilon)/epsilon.
In rhoPorousSimpleFoam solver and in EEqn.H file, these effects don’t exist. EEqn.H or energy equation for this solver is same to other solvers that don’t consider porous media (like rhoSimpleFoam). My question is about that.
How I can change the condition to that one? How I can consider the effect of porous properties in energy equation? How I can say to OepnFoam that in porous zone, conductivity and rho*cp is modified?
mohammad_kordo is offline   Reply With Quote

Old   September 6, 2016, 02:53
Default
  #2
New Member
 
Bibin K.S
Join Date: Oct 2014
Posts: 19
Rep Power: 12
bibin.sme is on a distinguished road
Hi

Could you find some solution to the above problem?
please if so please give me some idea...

Regards
Bibin
bibin.sme is offline   Reply With Quote

Old   February 19, 2018, 10:24
Default
  #3
New Member
 
Adri
Join Date: Sep 2017
Posts: 24
Rep Power: 9
Adri_12 is on a distinguished road
Hi,

I am already wondering if that is not more interesting to start from chtMultiRegionFoam solvers implementing a porousity model in fvoption and applying it only in the porous domain.

Any ideas about this issue ?

adrià
Adri_12 is offline   Reply With Quote

Old   February 21, 2018, 09:09
Default
  #4
Member
 
Anurag
Join Date: Aug 2014
Location: Germany
Posts: 57
Rep Power: 12
anuragm is on a distinguished road
Quote:
Originally Posted by Adri_12 View Post
Hi,

I am already wondering if that is not more interesting to start from chtMultiRegionFoam solvers implementing a porousity model in fvoption and applying it only in the porous domain.

Any ideas about this issue ?

adrià
But why start with chtMultiRegionFoam if there is only a need for porosity through fvOptions? The porosity can be defined only within a certain cellZone, for example.
anuragm is offline   Reply With Quote

Old   February 22, 2018, 09:16
Default
  #5
New Member
 
Adri
Join Date: Sep 2017
Posts: 24
Rep Power: 9
Adri_12 is on a distinguished road
I want to modelize heat transfers between air and porous zone considering two thermal equations (1 for the air and 1 for the solid porous zone), like mentionned in the post of mohammad amin.

Quote:
Originally Posted by mohammad_kordo View Post
Hi to all In momentum equation effect of porous with considering volume averaged equation, is a source term named Darcy-Forchheimer term. It is considered as source term that added to momentum in UEqn.H file that can be seen in porousSimpleFoam or rhoPorousSimpleFoam solvers. But in energy equation, porous modeling is based on tow model. One is local thermal equilibrium and another non-thermal equilibrium. In local thermal equilibrium the temperature of solid (porous) and fluid is same and effect of porous properties (density, specific heat coefficient and thermal conductivity) is considered in fluid energy equation. Term of “rhof*cpf” (multiply of density in specific heat coefficient of fluid) in energy equation is modified to
rhof*cpf+rhos*cps*(1-epsilon)/epsilon : f:fluid, s:solid.
Also conductivity (k) is modified from kf to kf+ks*(1-epsilon)/epsilon.
In rhoPorousSimpleFoam solver and in EEqn.H file, these effects don’t exist. EEqn.H or energy equation for this solver is same to other solvers that don’t consider porous media (like rhoSimpleFoam). My question is about that.
How I can change the condition to that one? How I can consider the effect of porous properties in energy equation? How I can say to OepnFoam that in porous zone, conductivity and rho*cp is modified?

that is why I mentionned the possibility to start with chtMultiRegionFoam. But I am wondering if the fvOptions defining a porous zone would be enough to consider the heat transfer + the Darcy behaviour of the fluid.

I hope it is more clear ?

Thanks for your help !

Adrià
Adri_12 is offline   Reply With Quote

Old   February 23, 2018, 08:09
Post
  #6
Member
 
Anurag
Join Date: Aug 2014
Location: Germany
Posts: 57
Rep Power: 12
anuragm is on a distinguished road
Quote:
Originally Posted by Adri_12 View Post
that is why I mentionned the possibility to start with chtMultiRegionFoam. But I am wondering if the fvOptions defining a porous zone would be enough to consider the heat transfer + the Darcy behaviour of the fluid.
Starting with the basics:

  1. The porous medium effects the momentum and temperature profiles in the flow
  2. If you have a perfect description of the porous medium in your geometry, you do not need Darcy or any other modifications
  3. However, in many cases porous media are non-uniform or just not worth the effort to describe accurately with all the details
  4. So it makes sense to approximate their effect indirectly, by modifying the momentum and temperature equations
  5. This modification can be achieved by either
    a. directly editing the UEqns.H and EEqns.H or equivalent files (requires writing equations & recompiling the solver)
    b. through the use of fvOptions (if solver includes support for it, then no recompilation needed)
  6. In case of 5b. you will still need to look at the way the equations are implemented in the solver to add a source term that is consistent with the rest of the equation. You will need source terms for both energy and momentum in this case.
I am working with porous mediums as well, but I have only implemented an additional momentum source using the fvOptions at the moment. I will add the term for the energy equation in the near future. If you want to do it as well, the following resource will help.

https://web.stanford.edu/~csoulain/O...G_PART_3v5.pdf
anuragm is offline   Reply With Quote

Old   March 1, 2018, 12:04
Default
  #7
New Member
 
Adri
Join Date: Sep 2017
Posts: 24
Rep Power: 9
Adri_12 is on a distinguished road
Quote:
Originally Posted by anuragm View Post
  1. The porous medium effects the momentum and temperature profiles in the flow
  2. This modification can be achieved by either
    a. directly editing the UEqns.H and EEqns.H or equivalent files (requires writing equations & recompiling the solver)
    b. through the use of fvOptions (if solver includes support for it, then no recompilation needed)
  3. In case of 5b. you will still need to look at the way the equations are implemented in the solver to add a source term that is consistent with the rest of the equation. You will need source terms for both energy and momentum in this case.
Thanks for your time and you answer and sorry for my late response. Your comments match with both possibilities I will try. I will give you more details when I get more results.

Adrià
Adri_12 is offline   Reply With Quote

Old   September 18, 2020, 12:56
Default
  #8
New Member
 
Yumeko
Join Date: Sep 2020
Posts: 2
Rep Power: 0
CriCri is on a distinguished road
Hello,
Did you find a solution to the above problem? I'm working on a similar case using the ChTmultiRegion solver but I'm not sure how to implement an additional source term on the energy equation (as far as I understood the darcy forcheimer source term acts only on the momentum equation)
thanks for your help!
CriCri is offline   Reply With Quote

Old   November 21, 2020, 05:36
Default
  #9
New Member
 
wanghongjie
Join Date: Apr 2020
Posts: 28
Rep Power: 6
wanghongjie is on a distinguished road
Quote:
Originally Posted by anuragm View Post
Starting with the basics:

  1. The porous medium effects the momentum and temperature profiles in the flow
  2. If you have a perfect description of the porous medium in your geometry, you do not need Darcy or any other modifications
  3. However, in many cases porous media are non-uniform or just not worth the effort to describe accurately with all the details
  4. So it makes sense to approximate their effect indirectly, by modifying the momentum and temperature equations
  5. This modification can be achieved by either
    a. directly editing the UEqns.H and EEqns.H or equivalent files (requires writing equations & recompiling the solver)
    b. through the use of fvOptions (if solver includes support for it, then no recompilation needed)
  6. In case of 5b. you will still need to look at the way the equations are implemented in the solver to add a source term that is consistent with the rest of the equation. You will need source terms for both energy and momentum in this case.
I am working with porous mediums as well, but I have only implemented an additional momentum source using the fvOptions at the moment. I will add the term for the energy equation in the near future. If you want to do it as well, the following resource will help.

https://web.stanford.edu/~csoulain/O...G_PART_3v5.pdf
Did you have add the energy equation to the porous..Foam? I have a similar question, and I only have a solution to compile cht-like Foam which includes fluid and solid. Maybe you have a better solution, help me, please.
wanghongjie is offline   Reply With Quote

Old   November 22, 2020, 08:18
Default
  #10
New Member
 
Utkan Erdem
Join Date: Sep 2020
Posts: 19
Rep Power: 6
Utkan is on a distinguished road
Quote:
Originally Posted by anuragm View Post
Starting with the basics:

  1. The porous medium effects the momentum and temperature profiles in the flow
  2. If you have a perfect description of the porous medium in your geometry, you do not need Darcy or any other modifications
  3. However, in many cases porous media are non-uniform or just not worth the effort to describe accurately with all the details
  4. So it makes sense to approximate their effect indirectly, by modifying the momentum and temperature equations
  5. This modification can be achieved by either
    a. directly editing the UEqns.H and EEqns.H or equivalent files (requires writing equations & recompiling the solver)
    b. through the use of fvOptions (if solver includes support for it, then no recompilation needed)
  6. In case of 5b. you will still need to look at the way the equations are implemented in the solver to add a source term that is consistent with the rest of the equation. You will need source terms for both energy and momentum in this case.
I am working with porous mediums as well, but I have only implemented an additional momentum source using the fvOptions at the moment. I will add the term for the energy equation in the near future. If you want to do it as well, the following resource will help.

https://web.stanford.edu/~csoulain/O...G_PART_3v5.pdf
I am also curious whether you were able to find a way to implement the energy equation.
Utkan is offline   Reply With Quote

Reply


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
Map of the OpenFOAM Forum - Understanding where to post your questions! wyldckat OpenFOAM 10 September 2, 2021 06:29
energy equation in porous media? Tony FLUENT 1 June 9, 2015 14:41
Definition of energy equation using UDF aestas Fluent UDF and Scheme Programming 15 April 28, 2015 04:33
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
Solving the total energy equation dohnie OpenFOAM 21 June 13, 2014 01:44


All times are GMT -4. The time now is 09:38.