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

Residuals in buoyantBoussinesqSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 28, 2013, 05:51
Default Residuals in buoyantBoussinesqSimpleFoam
  #1
New Member
 
Join Date: Apr 2011
Location: Magdeburg, Germany
Posts: 23
Rep Power: 6
Nico A. is on a distinguished road
Dear FOAMers,

I am simulating a (laminar) natural convection boundary layer in 2D with the above-named solver. So my domain has a uniformly heated wall at the left and atmospheric inlet and outlet conditions at the top, bottom and right respectivly. The domain is attached.

The case is already running with the following BCs for T, U, p_rgh:

0/T
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [ 0 0 0 1 0 0 0 ];

internalField   uniform 288.15;

boundaryField
{
    heatedWall
    {
        type            fixedValue;
        value           uniform 298.15;
    }
    right
    {
        type            inletOutlet;
        inletValue      uniform 288.15;
        value           uniform 288.15;
    }
   
    frontAndBack
    {
        type            empty;
    }
    inlet
    {
        type            inletOutlet;
        inletValue      uniform 288.15;
        value           uniform 288.15;
    }
    outlet
    {
        type            zeroGradient;
    }
}


// ************************************************************************* //
0/U
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform ( 0 0 0 );

boundaryField
{
    heatedWall
    {
        type            fixedValue;
        value           uniform ( 0 0 0 );
    }
    
    right
    {
        type            pressureDirectedInletVelocity;
    value        uniform (0 0 0);
    inletDirection  uniform (-1 0 0);
    }
    
    frontAndBack
    {
        type            empty;
    }
    
    inlet
    {
        type            pressureInletOutletVelocity;
    value        uniform (0 0 0);
    }
    
    outlet
    {
        type            inletOutlet;
        inletValue      uniform ( 0 0 0 );
        value           uniform ( 0 0 0 );        
    }
}

// ************************************************************************* //
0/p_rgh
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    heatedWall
    {
        type            buoyantPressure;
    rho        rhok;
        value           uniform 0;
    }
    
    right
    {
        type            totalPressure;
        U               U;
        phi             phi;
        rho             rhok;
        psi             none;
        gamma           1.4;
        p0              uniform 0;
    }
    
    frontAndBack
    {
        type            empty;
    }
    
    inlet
    {
        type            totalPressure;
        U               U;
        phi             phi;
        rho             rhok;
        psi             none;
        gamma           1.4;
        p0              uniform 0; 
    }
     
    outlet
    {
    type            buoyantPressure;
    rho        rhok;
        value           uniform 0;   
    }
}

// ************************************************************************* //
For the temperature difference of 10K and other flow parameters, the Rayleigh number is Ra=10^6, so a laminar boundary layer will develop, as you can see in the attached pictures for the Uz and T distribution at the vertical wall. The results seem to be physically correct.

But when I look at the residuals, they are oscillating (see attached picture). Normally they should fall due to given tolerances/ convergence criterions of 10^-08. The results are nearly the same for upwind and linear div-schemes.

I appreciate any suggestion regarding this case.

Best regards,
Nico
Attached Images
File Type: png domain.png (12.0 KB, 11 views)
File Type: jpg Uz.jpg (22.3 KB, 9 views)
File Type: jpg T.jpg (18.2 KB, 9 views)
File Type: png residuals.png (11.0 KB, 14 views)
Nico A. 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
Tutorial of how to plot residuals ! wolle1982 OpenFOAM Running, Solving & CFD 112 June 3, 2015 07:50
under-relaxation factors -> level of residuals Zigainer FLUENT 16 April 14, 2015 23:17
what to monitor besides residuals? franzdrs FLUENT 5 March 21, 2013 04:59
judging convergence through residuals MachZero Main CFD Forum 7 December 25, 2012 13:18
Convergence - scaled vs unscaled residuals HS FLUENT 1 November 7, 2005 06:45


All times are GMT -4. The time now is 12:24.