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

Buoyancy of the buoyantFOAM case

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 13, 2006, 00:58
Default Hi friends, I want to use t
  #1
Senior Member
 
Guoxiang
Join Date: Mar 2009
Posts: 109
Rep Power: 17
liugx212 is on a distinguished road
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.
liugx212 is offline   Reply With Quote

Old   June 13, 2006, 09:29
Default The energy equation, "hEqn.H",
  #2
Member
 
E. David Huckaby
Join Date: Mar 2009
Posts: 57
Rep Power: 17
david_h is on a distinguished road
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:
"<case>/constant/environmentalProperties"

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

good luck
Dave
david_h is offline   Reply With Quote

Old   June 13, 2006, 12:48
Default Hi Dave: Thank you very muc
  #3
Senior Member
 
Guoxiang
Join Date: Mar 2009
Posts: 109
Rep Power: 17
liugx212 is on a distinguished road
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.
liugx212 is offline   Reply With Quote

Old   June 13, 2006, 16:56
Default Hi Guoxiang! When I think o
  #4
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51
gschaider will become famous soon enoughgschaider will become famous soon enough
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.
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   June 13, 2006, 21:57
Default Hi Bernhard Gschaider, Anyw
  #5
Senior Member
 
Guoxiang
Join Date: Mar 2009
Posts: 109
Rep Power: 17
liugx212 is on a distinguished road
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.
liugx212 is offline   Reply With Quote

Old   June 14, 2006, 07:46
Default Hi Guoxiang! I don't think
  #6
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51
gschaider will become famous soon enoughgschaider will become famous soon enough
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).
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   June 14, 2006, 10:53
Default Hi Bernhard Gschaider, Man
  #7
Senior Member
 
Guoxiang
Join Date: Mar 2009
Posts: 109
Rep Power: 17
liugx212 is on a distinguished road
Hi Bernhard Gschaider,

Many thanks for your advice. I will dig continuly.

Best regards.
liugx212 is offline   Reply With Quote

Old   July 4, 2006, 11:23
Default Hi All, I wonder if anybody
  #8
Senior Member
 
Maka Mohu
Join Date: Mar 2009
Posts: 305
Rep Power: 18
maka is on a distinguished road
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,
Maka.
maka is offline   Reply With Quote

Old   July 4, 2006, 11:57
Default Hi Maka Mohu, I have a atta
  #9
Senior Member
 
Guoxiang
Join Date: Mar 2009
Posts: 109
Rep Power: 17
liugx212 is on a distinguished road
Hi Maka Mohu,

I have a attachement here. But I am not sure what it is you said.
The link is: hal.iwr.uni-heidelberg.de/~peter/Papers/habil.pdf
liugx212 is offline   Reply With Quote

Old   July 4, 2006, 13:28
Default Hi Guoxiang, I checked the
  #10
Senior Member
 
Maka Mohu
Join Date: Mar 2009
Posts: 305
Rep Power: 18
maka is on a distinguished road
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,
Maka
maka is offline   Reply With Quote

Old   February 13, 2007, 04:41
Default Hi, I met the error when I
  #11
mss
Guest
 
Posts: n/a
Hi,

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::error::abort()
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]
__libc_start_main
__gxx_personality_v0
------------------------------------------------

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,
Rita
  Reply With Quote

Old   May 9, 2007, 11:09
Default Hi, I am using buoyancyFoa
  #12
mss
Guest
 
Posts: n/a
Hi,

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,
Rita
  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
TotalPressure and buoyantFoam ariorus OpenFOAM Running, Solving & CFD 1 January 22, 2008 09:41
ThermoPhysicalProperties in buoyantFoam prashant24983 OpenFOAM Running, Solving & CFD 0 October 6, 2007 10:40
BuoyantFoam braennstroem OpenFOAM Running, Solving & CFD 22 September 19, 2007 17:55
BuoyantFoam in parallel sylvester OpenFOAM Running, Solving & CFD 0 July 19, 2007 06:38
Mixing elbow case water heat transfer calculation buoyantFoam benyamin1 OpenFOAM Running, Solving & CFD 0 January 14, 2006 10:25


All times are GMT -4. The time now is 04:15.