[agustinvo]

Hello,

I am using a modificated buoyantBoussinesqPimpleFoam, and when I run the simulation, I get something like this:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Reading field U

Reading field p_rgh

Reading field T

Reading/calculating face flux field phi

Reading thermophysical properties

Selecting incompressible transport model TempPolyLaw
Selecting incompressible thermal transport model Newtonian
Creating field rho

Creating field Cp

Creating field kThermal

Creating field alpha

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model LaunderSharmaKE
LaunderSharmaKECoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    sigmaEps        1.3;
}

Selecting turbulence model type RASThermalModel
Selecting RAS thermal turbulence model reynoldsAnalogy
reynoldsAnalogyCoeffs
{
    Prt             0.85;
}

Calculating field g.h

Creating derivated fields

Creating field U*U

Creating field U*T

Creating field U*p

Creating field p*T

Creating field U*U*U

Creating field U*U*T

Creating field U*U*p

Creating field Pr

Radiation model not active: radiationProperties not found
Selecting radiationModel none
No finite volume options present


SIMPLE: no convergence criteria found. Calculations will run for 12000 steps.


Starting time loop

wallShearStress wallShearStress1:
    processing wall patches: 
        heatedWall

wallShearStress wallShearStress1 output:
    min/max(heatedWall) = (-0.0003237043 0 0), (0 0 0)
wallShearStress wallShearStress1 output:
    writing field wallShearStress

calcFvcGrad wallGradT1 output:
    writing field wallGradT

Time = 2

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 3.563288e-12, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 8.45721e-12, No Iterations 6
DILUPBiCG:  Solving for T, Initial residual = 1, Final residual = 7.8688e-12, No Iterations 5
time step continuity errors : sum local = 1.96001, global = -0.01777022, cumulative = -0.01777022
DILUPBiCG:  Solving for epsilon, Initial residual = 0.9998042, Final residual = 2.899231e-12, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 1.746308e-12, No Iterations 3
Turbulent thermal diffusivity calculated
ExecutionTime = 6.55 s  ClockTime = 7 s

wallShearStress wallShearStress1 output:
    min/max(heatedWall) = (-0.3334798 -0.001289673 0), (-0.1195151 -3.818758e-07 0)
wallShearStress wallShearStress1 output:
    writing field wallShearStress

calcFvcGrad wallGradT1 output:
    writing field wallGradT

My question is why there is no Time=1? Ok, it's a steady simulation (so there is not a problem in not iterating there), but it should it in any case.

[akidess]

Did you set your deltaT to 2?

[agustinvo]

No, it is dT=1

[agustinvo]

Hi,

I found the error.

since my code came from a unsteady solver, I had something like this:

    Info<< "\nStarting time loop\n" << endl;

    while (runTime.loop())
    {
        

        // --- Pressure-velocity SIMPLE corrector loop
        while (simple.loop())
        {
        
            Info<< "Time = " << runTime.timeName() << nl << endl;
            
            #include "UEqn.H"
            #include "TEqn.H"

            // --- Pressure corrector loop
            //while (simple.correct())
            {
                #include "pEqn.H"
            }

when you actually need this:

   Info<< "\nStarting time loop\n" << endl;

    while (simple.loop())
    {
        // --- Pressure-velocity SIMPLE corrector loop
        {
        
            Info<< "Time = " << runTime.timeName() << nl << endl;
            
            #include "UEqn.H"
            #include "TEqn.H"

            // --- Pressure corrector loop
            //while (simple.correct())
            {
                #include "pEqn.H"
            }

Only one loop is necessary.