CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Buoyancy of the buoyantFOAM case (

liugx212 June 12, 2006 23:58

Hi friends, I want to use t
Hi friends,

I want to use the buoyantFOAM model to simulate the CO2 sequestration. But I have questions on it.

Firstly, the buoyantFOAM case is driven by buoyancy. I want to control the effect of this buoyancy such as decreasing or increasing. But I really do not know where and how to do this. Also, in this model, I found maybe state equation was solved. But I can not find the variable T in source code.

Secondely, I want to modify this model as the density-driven flow in porous media. So the Darcy`s law will be considered. Since the P (pressure) is solved by PISO algorithm. The problem is whether any state or energy equation will be added to buoyancyFOAM case for P?

Could any friend give me advice?

Thanks in advance.

david_h June 13, 2006 08:29

The energy equation, "hEqn.H",
The energy equation, "hEqn.H", is solved in terms of the internal enthalpy, "h". The temperature and several other derived thermodynamic quanties are calculated(updated) using "thermo->correct()". You can get access to the temperature using "thermo->T()". There are also two densities. The first, "rho" which is initialized in "createFields.H", and a second, calculated by thermodynamic class, "thermo->rho()". "rho" is initially calculated from mass conservation ("rhoEqn.H"). The weighted difference between the two is the continuity error.("compressibleContinuityErrors.H")

The gravity vector can be adjusted in:

The "reactingFoam" application provides an example of using FOAM for multispecies flow calculations.

good luck

liugx212 June 13, 2006 11:48

Hi Dave: Thank you very muc
Hi Dave:

Thank you very much.

I will follow you to try.

Also, if I want to change the propersties of the room, which is air, As you said, maybe I should change this air to coal. How should I do? Also, Could you please give help to explain the fields of "<case>/constant/thermophysicalProperties" ?

Thanks again.

gschaider June 13, 2006 15:56

Hi Guoxiang! When I think o
Hi Guoxiang!

When I think of coal, I think of black solid stuff that I buy in shops for barbecues or for ovens. I don't think that the assumptions made by the buoyantFoam (density changes due to temperature changes are important) are valid for that stuff. At least not on planet earth.

From your postings above I assume that you are talking about a mixture of air and CO2. I think that buoyantFoam is only applicable if that mixture is homogenous. Anyway: thermophysicalProperties is explained in chapter 8.1 of the user guide (it's not a long read, but very interesting). Should you have more specific questions after reading that chapter, I'm sure they will be answered here.

liugx212 June 13, 2006 20:57

Hi Bernhard Gschaider, Anyw
Hi Bernhard Gschaider,

Anyway, many thanks for your advice.

Actually, I want to modify buoyantFOAM to simulate CO2 sequestration. It means that I want to injection CO2 to coal Basin. Here, CO2 is the fluid flow but coal as the fixed media for CO2 concentration. Of course, this should be the porous flow with the CO2 concentration.

For now, I am thinking to develop a solve for this CO2 sequestration based on the buoyantFOAM or others. So I asked these questions which I am interested in.

Whatever, I also want to get your any more advice about this.

Thanks a lot.

gschaider June 14, 2006 06:46

Hi Guoxiang! I don't think
Hi Guoxiang!

I don't think that the buoyantFoam-solver is applicable for your problem at all. If I understand your problem correctly you'll want something like the twoLiquidMixingFoam (liquids are CO2 and air) extended to support porous Media (or with Navier removed and only Darcy).

liugx212 June 14, 2006 09:53

Hi Bernhard Gschaider, Man
Hi Bernhard Gschaider,

Many thanks for your advice. I will dig continuly.

Best regards.

maka July 4, 2006 10:23

Hi All, I wonder if anybody
Hi All,

I wonder if anybody has some written report or paper that discuss the compressible formulation used in buoyantFoam. Sorry for any inconvenience, but I only have experience with incompressible formulation. I want to shift to buoyantFoam so I can avoid the Boussinesq approx. I just need any doc that may help me derive the formulation my self so, I can modify the solver comfortably.

best regards,

liugx212 July 4, 2006 10:57

Hi Maka Mohu, I have a atta
Hi Maka Mohu,

I have a attachement here. But I am not sure what it is you said.
The link is:

maka July 4, 2006 12:28

Hi Guoxiang, I checked the
Hi Guoxiang,

I checked the report on the descriped link.
Thanks for the report. It includes lots of nice information about software development issues of unstructured solvers, but it seems to be targeted for multiphase formulation and I work with single phase. I posted my message here since the title of the thread was relevant and you were asking about the formulation of buoyantFoam.

What I needed is a reference that derives the formulation that is used in buoyantFoam. It is a compressible formulation that still uses a Poisson's equation for pressure. I'm also curios about the advantages, such formulation may have over the common compressible formulation where there is no Passion's equation for pressure. Thanks for you kind response.

best regards,

mss February 13, 2007 04:41

Hi, I met the error when I

I met the error when I run buoyantFoam case:

--> FOAM FATAL ERROR : Maximum number of iterations exceeded

From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.3/src/thermophysicalModels/specie/lnInclude/ specieThermoI.H at line 83.

FOAM aborting

Foam::error::printStack(Foam:: Ostream&)
Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::calculate()
Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::correct()
buoyantFoam [0x418c59]

I'm using this case for Argon gas with teperature from 300 K upto 20 000K:
Name Argon
nMoles 1.0
molecularWeight 39.948
deltaH 1.0499e+05
viscosity 3.4079e-5
Pr 0.65

Thanks so much,

mss May 9, 2007 10:09

Hi, I am using buoyancyFoa

I am using buoyancyFoam case. In the thermophysical properties I need to make some changes.
Originaly OpenFOAM is using : hThermo > pureMixture > constThermo > specieThermo > hConstThermo > perfectGas

In my case :
1. instead constThermo I should put sutherlandTransport library
2. instead hConstThermo I should put JanafThermo.

My question is how can I do the link between libraries in OpenFOAM?

Thank you for your help,

All times are GMT -4. The time now is 00:23.