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

Steady simulation does not iterate in Time=1

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 6, 2015, 07:20
Default Steady simulation does not iterate in Time=1
  #1
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 307
Rep Power: 13
agustinvo is on a distinguished road
Hello,

I am using a modificated buoyantBoussinesqPimpleFoam, and when I run the simulation, I get something like this:
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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.
agustinvo is offline   Reply With Quote

Old   November 6, 2015, 07:53
Default
  #2
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 28
akidess will become famous soon enough
Did you set your deltaT to 2?
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   November 6, 2015, 09:09
Default
  #3
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 307
Rep Power: 13
agustinvo is on a distinguished road
No, it is dT=1
agustinvo is offline   Reply With Quote

Old   November 19, 2015, 05:57
Default
  #4
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 307
Rep Power: 13
agustinvo is on a distinguished road
Hi,

I found the error.

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

Code:
    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:
Code:
   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.
agustinvo is offline   Reply With Quote

Reply

Tags
iteration, running

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
steady state simulation diverge Ema40 Fluent Multiphase 5 July 21, 2018 15:09
Adding heat source to chtMultiRegionFoam maddalena OpenFOAM Programming & Development 61 February 17, 2018 09:33
Under-relaxation for steady state simulation in CFX Chander CFX 7 May 1, 2014 13:44
Transient features in steady state simulation? PedFr0 Main CFD Forum 16 October 10, 2013 08:17
two phase flow steady VOF simulation yoyoyona ANSYS 0 March 4, 2012 07:09


All times are GMT -4. The time now is 16:59.