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

Natural ventilation in a room with windows

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

Like Tree3Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 27, 2014, 04:12
Default Natural ventilation in a room with windows
  #1
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Hi everybody,
I am looking for the best solution (solver and BC) to simulate the natural ventilation of a room.
I would like to take into consideration a room with opening surface on the outdoor, so an indoor temperature and an outdoor temperature at t=0.
In addition, I would also simulate an internal load like the radiation of a human body (in W/m2 ?).
(you can see the enclosed figure)
The wished result could be the evolution of the temperature fields and speed fields into the room.
Can anyone have already coded this kind of case ?? Maybe you can send me please JJ

Donít hesitate if you have any questions,
Thanks a lot,

Matthieu
Attached Images
File Type: jpg room_ventilation.jpg (40.7 KB, 94 views)
mathieu_135 is offline   Reply With Quote

Old   June 28, 2014, 02:06
Default Natural ventilation in a room with windows
  #2
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 106
Rep Power: 9
tian is on a distinguished road
Dear Matthieu,

you can try out the: buoyantBoussinesqSimpleFoam Solver. Should be enough for you case.

Bye
Thomas
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian is offline   Reply With Quote

Old   June 28, 2014, 04:15
Default
  #3
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Dear Thomas

thanks for your answer but do you know what kind of boundary conditions for the opening on outdoor (windows of the room) ?
Inlet / outlet ?
Pressure ?
...

thanks a lot in advance,

regards,

Mat
mathieu_135 is offline   Reply With Quote

Old   June 28, 2014, 05:43
Default Natural ventilation in a room with windows
  #4
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 106
Rep Power: 9
tian is on a distinguished road
Hi,

I am not sure about your window situation. Do you want open it a Little bit so the aiflow can come in or only set a surface temperature?

Regards,
Thomas
Attached Images
File Type: jpg BIMHVACTool-Window.jpg (97.3 KB, 106 views)
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian is offline   Reply With Quote

Old   June 28, 2014, 06:33
Default
  #5
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
in fact, I would like to consider a fully open window so, yes the airflow can come from outdoor to indoor.
In addition, I would like to take into account the radiation of an human body in the room for example (radiation in W/m2 ?).

Do you understand my wish ? Have you ever seen this case ??

Thanks for your help.

Mat
mathieu_135 is offline   Reply With Quote

Old   June 28, 2014, 09:58
Default
  #6
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 106
Rep Power: 9
tian is on a distinguished road
Dear Mat,

I upload one case. Maybe you can start from this case and adjust it by yourself...

It used buoyantBoussinesqSimpleFoam with k-Omega SST and radiations

https://copy.com/YukaGKRnnBds7kSY

Bye
Thomas
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian is offline   Reply With Quote

Old   June 30, 2014, 06:40
Default
  #7
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Hi Thomas,

thanks for your help

I don't manage to run your case with buoyantBoussinesqSimplefoam ... because I am not used to working with snappyhexmesh, I use Salome.


Can you explain me please ?
Thanks again,
Regards,

Mat

Last edited by mathieu_135; June 30, 2014 at 09:39.
mathieu_135 is offline   Reply With Quote

Old   July 3, 2014, 03:00
Default
  #8
Senior Member
 
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 104
Rep Power: 9
matejfor is on a distinguished road
snappy is just the mesher, you may prepare your mesh in any software.
What you need is to run your buoyantBoussinesqSimplefoam on your geometry. Have a look at the tutorial coming with the case.
Just be aware that the solver is using Boussinesq approximation to buoyant forces.

What will be the tricky part are the open windows. Unless you want to prescribe the pressure gradient between doors and window, you need the BC on window to take into account the pressure difference across the window. It would need some testing for a suitable BC settings. Maybe even use the buoyant solver without the Boussinesq approx. calculating the real density difference.

good luck
matej
matejfor is offline   Reply With Quote

Old   July 3, 2014, 06:03
Default
  #9
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Hi,

in fact, I don't manage to run this case, I have several error.

Do you manage to run this case on your computer ?

Regards,

Mat
mathieu_135 is offline   Reply With Quote

Old   July 3, 2014, 10:43
Default
  #10
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 106
Rep Power: 9
tian is on a distinguished road
Dear all,

I used OpenFOAM 2.2x. Let me check the case again.... What kind of error?

Bye
Thomas

[EDIT] I upload the case again:
https://copy.com/4HXgDQjMNCPkJHl9
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.

Last edited by tian; July 3, 2014 at 16:00.
tian is offline   Reply With Quote

Old   July 4, 2014, 05:53
Default
  #11
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Dear all,

I copy your files on my computer, I run buoyantBoussinesqSimpleFoam and I have this error :


--> FOAM FATAL IO ERROR:
keyword Cp is undefined in dictionary "/home/etu/OpenFOAM/etu-2.0.0/run/TestRoom_brut/0/T::boundaryField::body_People_body_10"


file: /home/etu/OpenFOAM/etu2.0.0/run/TestRoom_brut/0/T::boundaryField::body_People_body_10 from line 66 to line 70.

From function dictionary::lookupEntry(const word&, bool, bool) const in file db/dictionary/dictionary.C at line 400.

FOAM exiting


Thanks for your help

Mat
mathieu_135 is offline   Reply With Quote

Old   July 9, 2014, 01:48
Default
  #12
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Do you know this error ?

Thanks,

Mat
mathieu_135 is offline   Reply With Quote

Old   July 9, 2014, 02:49
Default
  #13
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 106
Rep Power: 9
tian is on a distinguished road
Hi,

yes: "keyword Cp is undefined in dictionary ". So add the "cp" value in "transportProperties" under "constant"-Folder like:

Cp Cp [0 0 0 0 0 0 0] 1000;

Bye
Thomas
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian is offline   Reply With Quote

Old   July 9, 2014, 03:09
Default
  #14
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Hi Thomas,

thanks for your answer. And this error, you know ? (sorry, i'm a newby ...)


********************************************
Create mesh for time = 0
 
Reading g
Reading thermophysical properties
Reading field T
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Creating turbulence model
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
alphaK1 0.85034;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.85616;
gamma1 0.5532;
gamma2 0.4403;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
c1 10;
}
Reading field kappat
Calculating field g.h
 
SIMPLE: convergence criteria
field p_rgh tolerance 0.01
field U tolerance 0.001
field T tolerance 0.001
field k tolerance 0.001
field omega tolerance 0.001
 
Starting time loop
Time = 1
 
 
--> FOAM FATAL IO ERROR:
Unknown convection scheme bounded
Valid convection schemes are :
1
(
Gauss
)
 
file: /home/etu/OpenFOAM/etu-2.0.0/run/TestRoom-A/system/fvSchemes::divSchemes::div(phi,U) at line 32.
From function convectionScheme<Type>::New(const fvMesh&, const surfaceScalarField&, Istream&)
in file /home/opencfd/OpenFOAM/OpenFOAM-2.0.0/src/finiteVolume/lnInclude/convectionScheme.C at line 97.
FOAM exiting
 
**********************************
mathieu_135 is offline   Reply With Quote

Old   July 9, 2014, 03:13
Default
  #15
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 106
Rep Power: 9
tian is on a distinguished road
Hi,

yes, because you used OpenFOA 2.0. Delete all "bounded" words in "fvScheme".

Bye
Thomas
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian is offline   Reply With Quote

Old   July 9, 2014, 03:45
Default
  #16
Senior Member
 
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 104
Rep Power: 9
matejfor is on a distinguished road
Hi the best thing to fight all these errors is:

a) look into tutorial case in $FOAM_TUTORIALS/heatTransfer/<solvername>/<tutorialCase>
for the settings and check your settings against it.

b) is there any reason why you are using 2.0 instead of the latest version 2.3? Where might be some improvements in the solver?

matej
matejfor is offline   Reply With Quote

Old   July 9, 2014, 03:53
Default
  #17
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
Hi,

I use VirtualBox (ubuntu), is there any solution to upgrade easily to OF 2.3 on a virtualbox ?

Thanks,
Mat
mathieu_135 is offline   Reply With Quote

Old   July 9, 2014, 04:02
Default
  #18
Senior Member
 
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 104
Rep Power: 9
matejfor is on a distinguished road
Yes, the easiest way is to follow this: http://www.openfoam.org/download/ubuntu.php from your VirtualBox ubuntu window. VirtualBox is just virtual computer, so you have ubuntu runing there just as it would be installed on your computer.

looking here: http://www.openfoam.org/download/history.php there has not been much development in solver itself, but in the BC you may find usefull for your case, like thermal baffles (2 different types), mixed thermal condition.

mat
matejfor is offline   Reply With Quote

Old   July 9, 2014, 06:13
Default
  #19
New Member
 
Join Date: Jan 2013
Posts: 13
Rep Power: 5
mathieu_135 is on a distinguished road
I have a lot of difficulties with my OF ....

Now :

************************************************** **********
Starting time loop
Time = 1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.0791183, No Iterations 161
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0792435, No Iterations 42
DILUPBiCG: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 0.0979706, No Iterations 9
 
--> FOAM FATAL ERROR:
lookup of (1|A(U)) from objectRegistry region0 successful
but it is not a volScalarField, it is a surfaceScalarField
From function objectRegistry::lookupObject<Type>(const word&) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.0.0/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 114.
FOAM aborting
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Foam::Ostream& Foam:perator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const& Foam:bjectRegistry::lookupObject<Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#4 Foam::fixedFluxPressureFvPatchScalarField::updateC oeffs() in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#5 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#6 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricFi eld<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const& Foam:bjectRegistry::lookupObject<Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#4 Foam::fixedFluxPressureFvPatchScalarField::updateC oeffs() in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#5 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#6 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricFi eld<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#8
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#9
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#10
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#11 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#12
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam"
Aborted
************************************************** **********
mathieu_135 is offline   Reply With Quote

Old   July 10, 2014, 02:15
Default
  #20
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 106
Rep Power: 9
tian is on a distinguished road
Dear Matej,

do you still try out my case or your own? Please compare your fvSchemes with a tutorial case like hotRooam. It seem your schemes are wrong for pressure?

Bye
Thomas
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian 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
Natural ventilation tichaona Main CFD Forum 4 December 9, 2011 15:34
CFX11 + Fortran compiler ? Mohan CFX 20 March 30, 2011 18:56
Natural convection + Species in a room Fran FLUENT 0 January 29, 2009 06:21
Natural Ventilation Boundary conditions lmitsc FLUENT 0 October 2, 2006 12:24
examples of simulation on room ventilation? lei FLUENT 1 June 20, 2003 22:38


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