[Sausi]

Hi,

I am using the chtMultiRegionFoam with OpenFoam 2.1.x and I realized that if I have a normal pipe flow the velocity at the lower edge is wrong [attachment]. My case is 2D and I defined a velocity of 3m/s at my inlet.
For my case the buoyant pressure is important because my real case has a heat source and I want to investigate this buoyant driven effects.



Pressure 0 folder

/*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
 |  \\    /   O peration     | Version:  2.1.x                                 |
 |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
 {
     version     2.0;
     format      binary;
     class       volScalarField;
     location    "0";
     object      p_rgh;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 

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

 value           101325;
 internalField   uniform $value;
 

 boundaryField
 {
     air_in
     {
         type            buoyantPressure;
         value           uniform $value;
     }
     air_out
     {
         type            fixedValue;
         value           uniform $value;
     }
     wall
     {
         type            buoyantPressure;
         value           uniform $value;
     }
     air_sym_left
     {
         type            symmetryPlane;
     }
     air_sym_right
     {
         type            symmetryPlane;
     }
 }
 

 

 // ************************************************************************* //

Velocity 0 folder

/*--------------------------------*- C++ -*----------------------------------*\ 
 | =========                 |                                                 | 
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           | 
 |  \\    /   O peration     | Version:  2.1.x                                 | 
 |   \\  /    A nd           | Web:      www.OpenFOAM.org                      | 
 |    \\/     M anipulation  |                                                 | 
 \*---------------------------------------------------------------------------*/ 
 FoamFile 
 { 
     version     2.0; 
     format      binary; 
     class       volVectorField; 
     location    "0"; 
     object      U; 
 } 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 
  
 dimensions      [ 0 1 -1 0 0 0 0 ]; 
  
 value           (3 0 0); 
 internalField   uniform $value; 
  
 boundaryField 
 { 
     air_in 
     { 
         type        fixedValue; 
         value       uniform $value; 
     } 
     air_out 
     { 
         type            inletOutlet; 
         inletValue      uniform (0 0 0); 
         value           uniform $value; 
     } 
     "wall" 
     { 
         type            fixedValue; 
         value           uniform (0 0 0); 
     } 
     air_sym_left 
     { 
         type            symmetryPlane; 
     } 
     air_sym_right 
     { 
         type            symmetryPlane; 
     } 
 } 
  
  
 // ************************************************************************* //

Did anyone have similar problems? Thanks for your help.

[jherb]

Hello Sausi,

could you solve your problem? Is gravity considered in your simulation?