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

air circulation with buoyantBoussinesqSimpleFoam solver

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 5, 2018, 14:52
Post air circulation with buoyantBoussinesqSimpleFoam solver
  #1
New Member
 
Aurélien Mingot
Join Date: Jun 2018
Location: London
Posts: 6
Rep Power: 8
Datsusong is on a distinguished road
Hi,
I am new on the forum and I will introduce myself.
My name is Aurélien and I am french student working in London in engineers building domain.
I work on dissertation. I want to have visualisation of air ventilation in a room with openfoam tool. The goal is to improve placement of supply and extract.

-First step was to create geometry of a room. This step is done. I've had some difficulties but it's OK. So BlockMeshDict file is correct and I can open on Paraview.

-I defined a air flowrate for output and velocity for input but I can't see any evolution on paraview.

-I use buoyantBoussinesqSimpleFoam solver. Is it the right one for air ventilation ?

-What kind of inout I missed ? I mean, basicly what kind of input we nee to run a simple air flow ?

I will send my files this week end
Thanks for reading me and sorry for my bad english ^^
Aurélien
Datsusong is offline   Reply With Quote

Old   June 5, 2018, 15:23
Default
  #2
Senior Member
 
tian's Avatar
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 122
Rep Power: 17
tian is on a distinguished road
Hi,

buoyantBoussinesqSimpleFoam is ok for this case. All other questions need to wait for your case.

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

Old   June 6, 2018, 09:07
Default
  #3
New Member
 
Aurélien Mingot
Join Date: Jun 2018
Location: London
Posts: 6
Rep Power: 8
Datsusong is on a distinguished road
Quote:
Originally Posted by tian View Post
Hi,

buoyantBoussinesqSimpleFoam is ok for this case. All other questions need to wait for your case.

Bye
Thomas
Hi, thanks for your reply
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices
(
(0 0 0)
(2.2 0 0)
(2.2 3 0)
(0 3 0)
(0 0 1.1)
(2.2 0 1.1)
(2.2 3 1.1)
(0 3 1.1)
(5.3 0 0)
(5.3 3 0)
(2.2 0 4.7)
(5.3 0 4.7)
(5.3 3 4.7)
(2.2 3 4.7)
(8.1 0 0)
(8.1 3 0)
(5.3 0 2.8)
(8.1 0 2.8)
(8.1 3 2.8)
(5.3 3 2.8)
(2.2 2.98 -0.1)
(5.3 2.98 -0.1)
(5.3 3 -0.1)
(2.2 3 -0.1)
(2.2 2.98 0)
(5.3 2.98 0)
(2.2 3 2.8)
(5.3 3.1 2.8)
(2.2 3.1 2.8)
(5.3 3.1 4.7)
(2.2 3.1 4.7)
);

blocks
(
hex (0 1 2 3 4 5 6 7) (50 50 50) simpleGrading (1 1 1)
hex (1 8 9 2 10 11 12 13) (50 50 50) simpleGrading (1 1 1)
hex (8 14 15 9 16 17 18 19) (50 50 50) simpleGrading (1 1 1)
hex (20 21 22 23 24 25 9 2) (10 10 1) simpleGrading (1 1 1)
hex (26 19 27 28 13 12 29 30) (50 50 50) simpleGrading (1 1 1)
);

edges
(
);

boundary

(
fixedWalls
{
type wall;
faces
(
(0 3 2 1)
(3 7 6 2)
(4 5 6 7)
(0 1 5 4)
(1 2 9 8)
(2 13 12 9)
(10 11 12 13)
(1 8 11 10)
(1 10 13 2)
(8 9 12 11)
(8 9 15 14)
(16 17 18 19)
(8 14 17 16)
(15 9 19 18)
(2 23 20 24)
(2 9 22 23)
(25 21 22 9)
(24 20 21 25)
(13 26 19 12)
(13 12 29 30)
(29 12 19 27)
(28 27 19 26)
(28 26 13 30)
);
}
supply
{
type inletOutlet;
faces
(
(20 21 22 23)
);
}
extract
{
type inletOutlet;
faces
(
(28 30 29 27)
);
}
);
mergePatchPairs
(
);


// ************************************************** *********************** //
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
supply
{
type flowRateOutletVelocity;
volumetricFlowRate 0.5;
value uniform (0 0 0);
}

extract
{
type flowRateOutletVelocity;
volumetricFlowRate 0.2;
value uniform (0 0 0);
}

fixedWalls
{
type noSlip;
}
defaultFaces
{
type calculated;
}
}

// ************************************************** *********************** //
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
supply
{
type calculated;
value $internalField;
}

extract
{
type calculated;
value $internalField;
}

fixedWalls
{
type calculated;
value $internalField;
}
defaultFaces
{
type calculated;
value $internalField;
}
}

// ************************************************** *********************** //
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
supply
{
type calculated;
rho rhok;
value $internalField;
}

extract
{
type fixedFluxPressure;
rho rhok;
value uniform 30;
}

fixedWalls
{
type fixedFluxPressure;
rho rhok;
value uniform 0;
}
defaultFaces
{
type fixedFluxPressure;
rho rhok;
value uniform 0;
}
}

// ************************************************** *********************** //
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{
supply
{
type fixedValue;
value uniform 280;
}
extract
{ type calculated;
value $internalField;
}
fixedWalls
{
type fixedValue;
value uniform 300;
}
defaultFaces
{ type calculated;
value $internalField;
}
}

// ************************************************** *********************** //
here are my files for run a simulation alpha epsilon nut and k stay constant. To have a good vision I add a picture of the room on paraview.
https://imgur.com/vOyHfNi

What input I missed ? perhaps my exploitation data on paraview is not good to have visualisation ?

Thanks
Have a good day
Aurélien
Datsusong is offline   Reply With Quote

Old   June 6, 2018, 13:25
Default
  #4
Senior Member
 
tian's Avatar
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 122
Rep Power: 17
tian is on a distinguished road
Hi,

your blockMesh looks not ok. It seems you created internal walls which will block your fluid. Also the cells size and number is also terrible.

https://postimg.cc/image/o3lhcj2dp/

https://postimg.cc/image/5b9m94t6l/

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

Old   June 7, 2018, 05:13
Default
  #5
New Member
 
Aurélien Mingot
Join Date: Jun 2018
Location: London
Posts: 6
Rep Power: 8
Datsusong is on a distinguished road
Hi,

Thanks for your help !!
ohhhh yes internal walls, I suppose I have to add empty boundary for internal walls ?

Because in my BlockMeshDict I have written mergePatchPairs and I thought it was ok for face between blocks.

Yeah I know that cells size is terrible but I will adjust after a first run.

Bye
Aurélien
Datsusong is offline   Reply With Quote

Old   June 7, 2018, 09:36
Default
  #6
New Member
 
Aurélien Mingot
Join Date: Jun 2018
Location: London
Posts: 6
Rep Power: 8
Datsusong is on a distinguished road
Hi,
I think I understand my mistake for mesh !

In my case a part of a face in a block have two functions and it's impossible I have to create new blocks.
I add pictures to explain. In the first photo there are 2 blocks.
https://imgur.com/kOz6dtB

but if I want to have a part of wall as empty and other part as wall I have to create the third block or not ?
https://imgur.com/T8QpIrP

Aurélien
Datsusong 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
what solver for air cyclone simulation comingdaytime OpenFOAM Running, Solving & CFD 2 May 20, 2015 12:12
Working directory via command line Luiz CFX 4 March 6, 2011 20:02
Getting too many iterations by velocity solving (aborting). Changing U - Solver? suitup OpenFOAM Running, Solving & CFD 0 January 20, 2010 07:45
air bubble is disappear increasing time using vof xujjun CFX 9 June 9, 2009 07:59
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


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