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

How to change simulated fluid from air to C02

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 27, 2011, 06:03
Default How to change simulated fluid from air to C02
  #1
New Member
 
Juanito
Join Date: Oct 2010
Posts: 27
Rep Power: 5
Aerospace is on a distinguished road
Hello everybody!

I would like to make a heat transfer simulation with BuoyantBoussinesqSimpleFoam but I need to employ CO2 as simulated fluid.
Which is the easiest way to change from air to CO2? Notice that BuoyantBoussinesqSimpleFoam doesn't employ thermophysicalProperties file.

Thank you very much in advance!
Aerospace is offline   Reply With Quote

Old   May 27, 2011, 06:47
Default
  #2
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 9
vkrastev is on a distinguished road
I'm not an expert of the buoyantBoussinesqSimpleFoam solver, but I think that you simply have to set properly your transportProperties file, changing the air properties with the CO2 properties (such as kinematic viscosity, thermal expansion coefficient and Prandtl number) at a reference temperature and pressure suitable for your case.

Regards

V.
vkrastev is offline   Reply With Quote

Old   May 27, 2011, 07:38
Default
  #3
New Member
 
Juanito
Join Date: Oct 2010
Posts: 27
Rep Power: 5
Aerospace is on a distinguished road
Quote:
Originally Posted by vkrastev View Post
I'm not an expert of the buoyantBoussinesqSimpleFoam solver, but I think that you simply have to set properly your transportProperties file, changing the air properties with the CO2 properties (such as kinematic viscosity, thermal expansion coefficient and Prandtl number) at a reference temperature and pressure suitable for your case.

Regards

V.

First, thank you very much for your answer. I agree with you but, what about the heat capacity (Cp)? I think that you can't change it in the transportProperties because there you can only define nu, beta, Tref and Pradtl numbers.
Aerospace is offline   Reply With Quote

Old   May 27, 2011, 08:40
Default
  #4
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,066
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
look for definition of prandtl
http://en.wikipedia.org/wiki/Prandtl_number

besides in source code, you have:
fvScalarMatrix TEqn
(
fvm::div(phi, T)
- fvm::Sp(fvc::div(phi), T)
- fvm::laplacian(kappaEff, T)
);

which kappaEff is calculated like below:

volScalarField kappaEff
(
"kappaEff",
turbulence->nu()/Pr + turbulence->nut()/Prt
);

so you just need to have pr and nu for energy equation calculation
nimasam is offline   Reply With Quote

Old   May 27, 2011, 08:44
Default
  #5
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 9
vkrastev is on a distinguished road
Quote:
Originally Posted by Aerospace View Post
First, thank you very much for your answer. I agree with you but, what about the heat capacity (Cp)? I think that you can't change it in the transportProperties because there you can only define nu, beta, Tref and Pradtl numbers.
Well, actually I think that the "information" carried by the Cp is included in the Prandtl numbers (laminar and turbulent). If you take a look into the TEqn.H file (see below) you will see that the thermal diffusivity employed in the temperature equation is simply derived as nu/Pr (and the same is done for its turbulent counterpart, which is added to the laminar one to obtain the "effective" diffusivity). The thermal diffusivity is also defined as mu*Cp/K, but if you have sufficiently accurate information about nu and Pr you don't need to know anything about Cp, mu or K to let the solver work properly.

Hope this helps

V.

{
volScalarField kappaEff
(
"kappaEff",
turbulence->nu()/Pr + turbulence->nut()/Prt

);

fvScalarMatrix TEqn
(
fvm::div(phi, T)
- fvm::Sp(fvc::div(phi), T)
- fvm::laplacian(kappaEff, T)
);

TEqn.relax();

eqnResidual = TEqn.solve().initialResidual();
maxResidual = max(eqnResidual, maxResidual);

rhok = 1.0 - beta*(T - TRef);
}
vkrastev is offline   Reply With Quote

Old   May 27, 2011, 08:45
Default
  #6
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 9
vkrastev is on a distinguished road
Oops, someone's been faster than me!
vkrastev is offline   Reply With Quote

Old   May 27, 2011, 12:26
Default
  #7
New Member
 
Juanito
Join Date: Oct 2010
Posts: 27
Rep Power: 5
Aerospace is on a distinguished road
Thank you very much, both are perfect answers to my question!!!!
Aerospace is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
problem with Min/max rho tH3f0rC3 OpenFOAM 7 February 23, 2013 05:37
air change effectiveness overhead HVAC system guillaume Main CFD Forum 0 February 27, 2007 20:20
Water vapour condensation in CFX-5.7.1 hdj CFX 1 November 27, 2005 07:15
Air Volume Change with Pressure Walt Muller Main CFD Forum 3 December 12, 2000 07:09
Simulated air velocity near ceiling Apple L S Chan Main CFD Forum 1 July 1, 1999 21:13


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