View Issue Details

IDProjectCategoryView StatusLast Update
0002608OpenFOAM[All Projects] Bugpublic2017-07-10 12:49
ReportermiriAssigned Tohenry 
PrioritynormalSeverityminorReproducibilityalways
Status resolvedResolutionfixed 
Product Versiondev 
Fixed in Version 
Summary0002608: Function 'residuals' doesn't write Uz
DescriptionHello,
I wanted to plot the residuals as described in:
https://cfd.direct/openfoam/user-guide/graphs-monitoring/

So I copied this snipped into the controlDict:

functions
{
    #includeFunc residuals;
}

and copied this file:
/opt/openfoam-dev/etc/caseDicts/postProcessing/numerical/residuals
without changes into the system-directory of the hotBoxes-Tutorial.

Now when running reactingParcelFilmFoam the residuals.dat file gets written into postProcessing.

The problem is, that the Uz-field doesn't get written and the residuals.dat file looks like the one attached.

As a result foamMonitor writes all fields with the names data1, data2, data3...

Steps To Reproducecopy
    functions
    {
        #includeFunc residuals;
    }
into the controlDict

and copy
/opt/openfoam-dev/etc/caseDicts/postProcessing/numerical/residuals
into the system-directory of the hotBoxes-Tutorial

run foamMonitor
TagsNo tags attached.

Relationships

related to 0002547 closedhenry CLI - residual plotting 

Activities

miri

2017-07-09 18:52

reporter  

residuals.dat (867 bytes)

miri

2017-07-09 18:52

reporter  

foamMonitor.png (14,118 bytes)
foamMonitor.png (14,118 bytes)

miri

2017-07-09 18:53

reporter  

residuals (846 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | Web:      www.OpenFOAM.org
     \\/     M anipulation  |
-------------------------------------------------------------------------------
Description
    For specified fields, writes out the initial residuals for the first
    solution of each time step; for non-scalar fields (e.g. vectors), writes
    the largest of the residuals for each component (e.g. x, y, z).

\*---------------------------------------------------------------------------*/

#includeEtc "caseDicts/postProcessing/numerical/residuals.cfg"

fields (p U);

// ************************************************************************* //
residuals (846 bytes)

miri

2017-07-09 18:53

reporter  

controlDict (1,368 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  dev                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     reactingParcelFilmFoam;

startForm       latestTime;

startTime       0;

stopAt          endTime;

endTime         2;

deltaT          1e-04;

writeControl    adjustableRunTime;

writeInterval   0.1;

purgeWrite      0;

writeFormat     ascii; // binary;

writePrecision  10;

writeCompression uncompressed;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

adjustTimeStep  yes;

maxCo           0.5;

maxDeltaT       0.002;

functions
{
    #includeFunc residuals;
}

// ************************************************************************* //
controlDict (1,368 bytes)

miri

2017-07-09 18:56

reporter  

log (51,300 bytes)
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  dev                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : dev-1888b122cc34
Exec   : reactingParcelFilmFoam
Date   : Jul 09 2017
Time   : 19:26:45
Host   : "test-VirtualBox"
PID    : 2962
I/O    : uncollated
Case   : /home/test/OpenFOAM/test-dev/run/hotBoxes
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Creating combustion model

Selecting combustion model laminar<psiChemistryCombustion>
Selecting chemistry type 
{
    chemistrySolver noChemistrySolver;
    chemistryThermo psi;
}

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

Selecting chemistryReader foamChemistryReader
    elements not defined in "/home/test/OpenFOAM/test-dev/run/hotBoxes/constant/foam.inp"
chemistryModel: Number of species = 3 and reactions = 0
    using integrated reaction rate
Creating component thermo properties:
    multi-component carrier - 3 species
    liquids - 1 components
    solids - 0 components
Creating field rho


Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
Creating field dpdt

Creating field kinetic energy K


Reading g

Reading hRef
Calculating field g.h

No MRF models present


Constructing reacting cloud
Constructing particle forces
    Selecting particle force sphereDrag
    Selecting particle force gravity
Constructing cloud functions
    none
Constructing particle injection models
Creating injector: model1
Selecting injection model coneInjection
    Constructing 3-D injection
Selecting distribution model RosinRammler
Selecting dispersion model none
Selecting patch interaction model standardWallInteraction
Selecting stochastic collision model none
Selecting surface film model thermoSurfaceFilm
    Applying splashBai interaction model
Selecting U integration scheme Euler
Selecting heat transfer model none
Selecting T integration scheme analytical
Selecting composition model singlePhaseMixture
Selecting phase change model none
Selecting radiationModel none

Constructing surface film model
Selecting surfaceFilmModel thermoSingleLayer
    Selecting region model functions
    none
    Selecting filmThermoModel liquid
    Selecting film injection models
        curvatureSeparation
        drippingInjection
Selecting distribution model RosinRammler
    Selecting film transfer models
        none
    Selecting filmTurbulenceModel laminar
    Selecting film force models
        thermocapillary
        distributionContactAngle
Selecting distribution model normal
    Selecting filmViscosityModel liquid
    Selecting heatTransferModel mappedConvectiveHeatTransfer
    Selecting heatTransferModel constant
    Selecting phaseChangeModel standardPhaseChange
    Selecting radiationModel none
No finite volume options present

Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Film max Courant number: 0
deltaT = 0.0001436781609
Time = 0.000143678


Solving 3-D cloud reactingCloud1
Cloud: reactingCloud1
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    model1:
        number of parcels added     = 0
        mass introduced             = 0
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 0, 0
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.9912919068, Final residual = 1.793914086e-17, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.9904653317, Final residual = 1.986328704e-17, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.9973363737, Final residual = 1.832228392e-17, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.000125581383, Final residual = 6.782269813e-07, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 0, 3.286543803e-19
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.999999712, Final residual = 3.010999246e-08, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.9999997108, Final residual = 3.029136101e-08, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 1, Final residual = 2.355493116e-09, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 4.387856388e-05, Final residual = 6.518119935e-13, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.9851939041, Final residual = 2.635852857e-09, No Iterations 1
T gas min/max   = 300, 300
GAMG:  Solving for p_rgh, Initial residual = 0.003766716204, Final residual = 0.0002246622948, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.908922321e-08, global = -1.507838742e-11, cumulative = -1.507838742e-11
GAMG:  Solving for p_rgh, Initial residual = 0.0002243899932, Final residual = 4.088293376e-07, No Iterations 4
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.462470582e-11, global = -2.39276453e-14, cumulative = -1.510231506e-11
DILUPBiCGStab:  Solving for epsilon, Initial residual = 2.161907638e-05, Final residual = 3.153210653e-13, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.9999999395, Final residual = 2.261030041e-08, No Iterations 1
ExecutionTime = 31.48 s  ClockTime = 33 s

Courant Number mean: 1.524538238e-05 max: 2.745249783e-05
Film max Courant number: 0
deltaT = 0.0002063147145
Time = 0.000349993


Solving 3-D cloud reactingCloud1
Cloud: reactingCloud1
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    model1:
        number of parcels added     = 0
        mass introduced             = 0
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 0, 0
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.9999997751, Final residual = 1.730429965e-17, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.9999997479, Final residual = 1.573065806e-17, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.9999999928, Final residual = 3.80081344e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0002007270347, Final residual = 1.004836608e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 1.860799133e-20, 3.571948712e-12
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.4179561026, Final residual = 2.270798364e-08, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.4178584889, Final residual = 2.267262541e-08, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.9131622376, Final residual = 1.280773484e-08, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 1.751449989e-05, Final residual = 1.009220634e-12, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.9999531508, Final residual = 4.656763504e-08, No Iterations 1
T gas min/max   = 299.9999882, 300.0000115
GAMG:  Solving for p_rgh, Initial residual = 0.01209393665, Final residual = 0.0008557312455, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.215045997e-08, global = 1.971139215e-11, cumulative = 4.609077086e-12
GAMG:  Solving for p_rgh, Initial residual = 0.0008617705561, Final residual = 3.407065647e-07, No Iterations 6
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.909864338e-11, global = -4.857900982e-16, cumulative = 4.608591296e-12
DILUPBiCGStab:  Solving for epsilon, Initial residual = 1.881621782e-05, Final residual = 6.25164579e-13, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.340644112, Final residual = 1.024222128e-08, No Iterations 1
ExecutionTime = 41.88 s  ClockTime = 53 s

Courant Number mean: 5.083992092e-05 max: 0.0001141046605
Film max Courant number: 0
deltaT = 0.0002965774022
Time = 0.00064657


Solving 3-D cloud reactingCloud1
Cloud: reactingCloud1
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    model1:
        number of parcels added     = 0
        mass introduced             = 0
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 0, 0
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.853064311, Final residual = 1.577545174e-17, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.845556154, Final residual = 1.625203592e-17, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.7120700499, Final residual = 2.158519236e-17, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0002086864523, Final residual = 1.069739111e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 2.930842408e-18, 3.002108484e-11
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.3654801555, Final residual = 3.462562813e-08, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.3557411888, Final residual = 3.618928438e-08, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.7432654016, Final residual = 2.27566734e-08, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 1.288750158e-05, Final residual = 1.65323269e-12, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.9904586946, Final residual = 8.905303121e-08, No Iterations 1
T gas min/max   = 299.9991509, 300.0008496
GAMG:  Solving for p_rgh, Initial residual = 0.03632635406, Final residual = 0.002477847607, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.022394361e-07, global = 1.615295688e-10, cumulative = 1.6613816e-10
GAMG:  Solving for p_rgh, Initial residual = 0.002558362649, Final residual = 6.251797752e-07, No Iterations 7
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.218177075e-11, global = -1.42025442e-13, cumulative = 1.659961346e-10
DILUPBiCGStab:  Solving for epsilon, Initial residual = 2.715939298e-05, Final residual = 2.079673225e-12, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.2300674946, Final residual = 1.666429968e-08, No Iterations 1
ExecutionTime = 52.22 s  ClockTime = 76 s

Courant Number mean: 0.000124275936 max: 0.0003372829596
Film max Courant number: 0
deltaT = 0.0004245873065
Time = 0.00107116


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 2 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 2
    Current mass in system          = 0.0005355787921
    Linear momentum                 = (0.0003164004759 -0.0004857027585 -0.001450756416)
   |Linear momentum|                = 0.001562277379
    Linear kinetic energy           = 0.002402968041
    model1:
        number of parcels added     = 2
        mass introduced             = 0.0005355787921
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.6771630271, Final residual = 1.342035268e-17, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.6601497499, Final residual = 1.350425899e-17, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.5117201214, Final residual = 1.144018566e-17, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001937301813, Final residual = 1.374591894e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 3.035133984e-18, 1.35208614e-10
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.3266887609, Final residual = 1.518173882e-07, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.3078447775, Final residual = 1.555738102e-07, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.5532380832, Final residual = 9.727082612e-08, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 1.767973301e-05, Final residual = 5.078108947e-12, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.6017018428, Final residual = 9.23124487e-08, No Iterations 1
T gas min/max   = 299.9971475, 300.0022663
GAMG:  Solving for p_rgh, Initial residual = 0.09972176151, Final residual = 0.006699941019, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.96729416e-07, global = 2.159394746e-09, cumulative = 2.32539088e-09
GAMG:  Solving for p_rgh, Initial residual = 0.007370393307, Final residual = 8.602005625e-07, No Iterations 8
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.403166964e-11, global = 2.665988243e-12, cumulative = 2.328056869e-09
DILUPBiCGStab:  Solving for epsilon, Initial residual = 4.76390648e-05, Final residual = 6.100029229e-10, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.5467419413, Final residual = 7.385612475e-11, No Iterations 1
ExecutionTime = 62.89 s  ClockTime = 99 s

Courant Number mean: 0.0002571648565 max: 0.02723408886
Film max Courant number: 0
deltaT = 0.0006069254136
Time = 0.00167808


Solving 3-D cloud reactingCloud1
Cloud: reactingCloud1
    Current number of parcels       = 2
    Current mass in system          = 0.0005355787921
    Linear momentum                 = (0.0003154366317 -0.0004844374015 -0.001449733857)
   |Linear momentum|                = 0.001560739607
    Linear kinetic energy           = 0.00239776779
    model1:
        number of parcels added     = 2
        mass introduced             = 0.0005355787921
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.5167803697, Final residual = 1.202437035e-17, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.5058052084, Final residual = 1.54346951e-17, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.3799194506, Final residual = 8.574946698e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001912458769, Final residual = 1.483033371e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 1.781356341e-16, 4.259605346e-10
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.2866046139, Final residual = 4.196626304e-06, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.2781591761, Final residual = 5.656293095e-10, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.3893171839, Final residual = 4.505135113e-07, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 1.318873363e-05, Final residual = 8.320296031e-12, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.4723370681, Final residual = 1.373500842e-07, No Iterations 1
T gas min/max   = 299.990842, 300.0042705
GAMG:  Solving for p_rgh, Initial residual = 0.2786317296, Final residual = 0.01299865219, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.358531006e-07, global = 1.501288462e-08, cumulative = 1.734094149e-08
GAMG:  Solving for p_rgh, Initial residual = 0.01780210003, Final residual = 6.923372135e-07, No Iterations 12
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.359023717e-11, global = 6.151874981e-12, cumulative = 1.734709337e-08
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.0003684006135, Final residual = 5.180936195e-08, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.8680268641, Final residual = 5.401165254e-09, No Iterations 1
ExecutionTime = 74.49 s  ClockTime = 116 s

Courant Number mean: 0.0004494954439 max: 0.05882150995
Film max Courant number: 0
deltaT = 0.0008701054602
Time = 0.00254819


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 2 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 4
    Current mass in system          = 0.001274094229
    Linear momentum                 = (-0.0006356994802 -0.000442263073 -0.0031742132)
   |Linear momentum|                = 0.003267313865
    Linear kinetic energy           = 0.005695832398
    model1:
        number of parcels added     = 4
        mass introduced             = 0.001274094229
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.3743465194, Final residual = 8.033617082e-18, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.3709747016, Final residual = 1.182570289e-17, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.2663496403, Final residual = 9.526895854e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001912677875, Final residual = 1.722740304e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 7.032790538e-16, 9.63641482e-10
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.2184096655, Final residual = 3.786899622e-08, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.2457234271, Final residual = 1.000283553e-07, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.2317827989, Final residual = 3.640388859e-08, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 1.821380058e-05, Final residual = 2.704289121e-11, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.3942902936, Final residual = 3.054617149e-07, No Iterations 1
T gas min/max   = 299.9862001, 300.0071402
GAMG:  Solving for p_rgh, Initial residual = 0.9416154982, Final residual = 0.05055889382, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.575094562e-06, global = 1.217028726e-07, cumulative = 1.39049966e-07
GAMG:  Solving for p_rgh, Initial residual = 0.1137086347, Final residual = 9.494671855e-07, No Iterations 18
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.215854705e-11, global = -3.885717525e-12, cumulative = 1.390460802e-07
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.02170932406, Final residual = 7.124762374e-08, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.9156937898, Final residual = 5.814984233e-06, No Iterations 1
ExecutionTime = 87.78 s  ClockTime = 132 s

Courant Number mean: 0.0005950820354 max: 0.1652498025
Film max Courant number: 0
deltaT = 0.001233567235
Time = 0.00378176


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 2 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 6
    Current mass in system          = 0.001890877846
    Linear momentum                 = (-0.0009194844489 -0.0006004942013 -0.004725088832)
   |Linear momentum|                = 0.004851031788
    Linear kinetic energy           = 0.008467733285
    model1:
        number of parcels added     = 6
        mass introduced             = 0.001890877846
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.2016375836, Final residual = 1.056986874e-17, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.1992194727, Final residual = 1.153646267e-17, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.1357433112, Final residual = 4.658988095e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001908415269, Final residual = 1.199849563e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 1.074222246e-15, 1.455326632e-09
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.2207326822, Final residual = 3.229329252e-06, No Iterations 3
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.3304007722, Final residual = 1.17503265e-06, No Iterations 3
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.1652203261, Final residual = 8.078769707e-06, No Iterations 2
DILUPBiCGStab:  Solving for O2, Initial residual = 1.338029602e-05, Final residual = 1.13110927e-09, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.2975589734, Final residual = 5.244055867e-06, No Iterations 1
T gas min/max   = 299.9852851, 300.0109699
GAMG:  Solving for p_rgh, Initial residual = 0.7017982406, Final residual = 0.05358238757, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.468345257e-06, global = 1.989094897e-07, cumulative = 3.3795557e-07
GAMG:  Solving for p_rgh, Initial residual = 0.03626031201, Final residual = 5.680620663e-07, No Iterations 21
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.126720868e-11, global = -3.147579283e-11, cumulative = 3.379240942e-07
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.177289353, Final residual = 2.630929777e-06, No Iterations 3
DILUPBiCGStab:  Solving for k, Initial residual = 0.6676722838, Final residual = 4.067290019e-06, No Iterations 3
ExecutionTime = 102.33 s  ClockTime = 149 s

Courant Number mean: 0.0003630463357 max: 0.1020527478
Film max Courant number: 0
deltaT = 0.001718182934
Time = 0.00549994


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 4 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 10
    Current mass in system          = 0.002749969313
    Linear momentum                 = (0.0001238453605 -0.0005202793152 -0.006827845738)
   |Linear momentum|                = 0.006848759425
    Linear kinetic energy           = 0.01226637279
    model1:
        number of parcels added     = 10
        mass introduced             = 0.002749969313
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.02960168823, Final residual = 8.929440322e-18, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.03001242313, Final residual = 8.596503971e-18, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.01876854453, Final residual = 3.789202812e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001931329698, Final residual = 1.398247745e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 1.025362006e-15, 1.538927607e-09
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.7063234912, Final residual = 1.898172439e-06, No Iterations 4
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.5523807914, Final residual = 7.650142041e-06, No Iterations 3
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.8075216818, Final residual = 2.755654e-06, No Iterations 3
DILUPBiCGStab:  Solving for O2, Initial residual = 1.797637276e-05, Final residual = 3.547962454e-09, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.1568392338, Final residual = 3.162840351e-06, No Iterations 2
T gas min/max   = 299.9852535, 300.0142316
GAMG:  Solving for p_rgh, Initial residual = 0.1092976544, Final residual = 0.005376141373, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.034328402e-07, global = 8.325786448e-09, cumulative = 3.462498806e-07
GAMG:  Solving for p_rgh, Initial residual = 0.01133531387, Final residual = 9.749352355e-07, No Iterations 21
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.727358483e-11, global = -6.722949517e-11, cumulative = 3.461826511e-07
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.2053350499, Final residual = 7.638619578e-06, No Iterations 3
DILUPBiCGStab:  Solving for k, Initial residual = 0.4374804583, Final residual = 6.558653358e-06, No Iterations 4
ExecutionTime = 117.38 s  ClockTime = 165 s

Courant Number mean: 0.0004003326454 max: 0.1209326295
Film max Courant number: 0
deltaT = 0.001968751279
Time = 0.00746869


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 4 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 14
    Current mass in system          = 0.003734344953
    Linear momentum                 = (0.0004782709281 -0.001496928448 -0.009409415578)
   |Linear momentum|                = 0.009539740006
    Linear kinetic energy           = 0.01659357267
    model1:
        number of parcels added     = 14
        mass introduced             = 0.003734344953
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.01628035342, Final residual = 8.508540245e-18, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.01534590802, Final residual = 8.322771625e-18, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.01031921346, Final residual = 3.822908842e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001944748044, Final residual = 9.180557252e-07, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 9.475311507e-16, 1.492811515e-09
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.4323204804, Final residual = 8.345363091e-06, No Iterations 3
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.459363298, Final residual = 7.435280613e-06, No Iterations 4
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.8116034889, Final residual = 5.128385326e-06, No Iterations 4
DILUPBiCGStab:  Solving for O2, Initial residual = 1.56341456e-05, Final residual = 4.379710537e-09, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.0160919992, Final residual = 1.062637011e-06, No Iterations 2
T gas min/max   = 299.9855401, 300.013905
GAMG:  Solving for p_rgh, Initial residual = 0.7139985876, Final residual = 0.03419010994, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.149008028e-06, global = 7.799035168e-08, cumulative = 4.241730028e-07
GAMG:  Solving for p_rgh, Initial residual = 0.1094455548, Final residual = 8.511017492e-07, No Iterations 34
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.517846981e-11, global = 1.32790135e-11, cumulative = 4.241862818e-07
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.1876173135, Final residual = 2.799151792e-06, No Iterations 4
DILUPBiCGStab:  Solving for k, Initial residual = 0.2415347191, Final residual = 8.007128523e-06, No Iterations 4
ExecutionTime = 136.01 s  ClockTime = 184 s

Courant Number mean: 0.0006771370084 max: 0.1771335188
Film max Courant number: 0
deltaT = 0.001968751279
Time = 0.00943744


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 4 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 18
    Current mass in system          = 0.004718720592
    Linear momentum                 = (0.0007108337238 -0.001856449729 -0.01197241308)
   |Linear momentum|                = 0.0121363242
    Linear kinetic energy           = 0.02092610733
    model1:
        number of parcels added     = 18
        mass introduced             = 0.004718720592
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.03574893192, Final residual = 9.730905859e-18, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.03420715854, Final residual = 8.377491508e-18, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.02262464733, Final residual = 4.808408907e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001955047204, Final residual = 1.267382786e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 8.892371221e-16, 1.398943581e-09
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.214837119, Final residual = 3.117352949e-06, No Iterations 4
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.2689653003, Final residual = 5.650176939e-06, No Iterations 4
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.3118950253, Final residual = 2.734786004e-06, No Iterations 4
DILUPBiCGStab:  Solving for O2, Initial residual = 1.570758351e-05, Final residual = 2.784293583e-08, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.1014844272, Final residual = 1.023929514e-06, No Iterations 2
T gas min/max   = 299.9879522, 300.0116033
GAMG:  Solving for p_rgh, Initial residual = 0.8792101132, Final residual = 0.03838567975, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.277461879e-06, global = 4.846273154e-08, cumulative = 4.726490134e-07
GAMG:  Solving for p_rgh, Initial residual = 0.05256280098, Final residual = 8.679862221e-07, No Iterations 31
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.03596089e-10, global = 3.791173844e-11, cumulative = 4.726869251e-07
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.08273113257, Final residual = 1.245751872e-06, No Iterations 4
DILUPBiCGStab:  Solving for k, Initial residual = 0.1867730679, Final residual = 9.531053723e-06, No Iterations 4
ExecutionTime = 154.17 s  ClockTime = 203 s

Courant Number mean: 0.0003127922614 max: 0.1535993523
Film max Courant number: 0
deltaT = 0.001968751279
Time = 0.0114062


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 4 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 22
    Current mass in system          = 0.005703096231
    Linear momentum                 = (0.0007795546931 -0.001307647942 -0.0147449412)
   |Linear momentum|                = 0.01482332418
    Linear kinetic energy           = 0.02524631132
    model1:
        number of parcels added     = 22
        mass introduced             = 0.005703096231
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.004764941147, Final residual = 1.00052226e-17, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.00464327307, Final residual = 8.987383096e-18, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.002945199598, Final residual = 4.91186358e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001981299824, Final residual = 8.313758097e-07, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 8.846262043e-16, 1.387673319e-09
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.4393719856, Final residual = 2.934842393e-06, No Iterations 7
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.4033347183, Final residual = 1.410795316e-06, No Iterations 7
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.738394929, Final residual = 9.291964057e-06, No Iterations 6
DILUPBiCGStab:  Solving for O2, Initial residual = 1.577808676e-05, Final residual = 1.785714971e-07, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.07500484785, Final residual = 4.970899758e-06, No Iterations 2
T gas min/max   = 299.9859934, 300.0101506
GAMG:  Solving for p_rgh, Initial residual = 0.2066373205, Final residual = 0.01352917403, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.784229325e-07, global = -6.887242966e-09, cumulative = 4.657996821e-07
GAMG:  Solving for p_rgh, Initial residual = 0.04555389435, Final residual = 8.295865495e-07, No Iterations 27
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.812567447e-11, global = 3.809574032e-11, cumulative = 4.658377779e-07
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.06368126834, Final residual = 6.772523588e-06, No Iterations 3
DILUPBiCGStab:  Solving for k, Initial residual = 0.2736861211, Final residual = 4.598812498e-06, No Iterations 7
ExecutionTime = 172.71 s  ClockTime = 226 s

Courant Number mean: 0.0002683984753 max: 0.186715495
Film max Courant number: 0
deltaT = 0.001968751279
Time = 0.0133749


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 4 new parcels

Cloud: reactingCloud1
    Current number of parcels       = 26
    Current mass in system          = 0.006687471871
    Linear momentum                 = (0.001441432611 -0.0007975278285 -0.01746003534)
   |Linear momentum|                = 0.01753757716
    Linear kinetic energy           = 0.02958050067
    model1:
        number of parcels added     = 26
        mass introduced             = 0.006687471871
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    New film splash parcels         = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.001858599765, Final residual = 9.895607892e-18, No Iterations 1
smoothSolver:  Solving for Ufy, Initial residual = 0.001847191566, Final residual = 8.499616814e-18, No Iterations 1
smoothSolver:  Solving for Ufz, Initial residual = 0.001260126205, Final residual = 5.115656128e-18, No Iterations 1
smoothSolver:  Solving for hf, Initial residual = 0.0001960784919, Final residual = 1.07254366e-06, No Iterations 100
smoothSolver:  Solving for deltaf, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
    added mass         = 0
    current mass       = 0
    min/max(mag(U))    = 8.817350071e-16, 1.392951418e-09
    min/max(delta)     = 0, 0
    coverage           = 0
    injected mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    transferred mass      = 0
      - patch: region0_to_wallFilmRegion_cubeFaces: 0
      - patch: region0_to_wallFilmRegion_floorFaces: 0
    min/mean/max(T)    = 300, 300, 300
    mass phase change  = 0
    vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.4319352657, Final residual = 3.503566308e-06, No Iterations 8
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.4169593406, Final residual = 3.72562284e-06, No Iterations 8
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.6601804208, Final residual = 7.663666336e-07, No Iterations 8
DILUPBiCGStab:  Solving for O2, Initial residual = 1.582359533e-05, Final residual = 2.484976279e-07, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.005734836191, Final residual = 7.856249973e-06, No Iterations 2
T gas min/max   = 299.9858623, 300.0102451
GAMG:  Solving for p_rgh, Initial residual = 0.6884617283, Final residual = 0.06257313158, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.208157059e-06, global = 1.233555493e-08, cumulative = 4.781733328e-07
log (51,300 bytes)

wyldckat

2017-07-09 22:43

updater   ~0008369

Actually it's the "p" entry that doesn't have any values for it. If you look at the log file, there doesn't seem to be any "p" equation being solved.

This is similar to what happened in report #2547.

From one perspective, this is a user error; from another perspective, if the entry "p" does not exist, then an explicit entry should be placed on the outputted file that indicates that it wasn't found, such as "N/A" or "-".

wyldckat

2017-07-09 23:11

updater  

proposition_v1.patch (953 bytes)
diff --git a/src/functionObjects/utilities/residuals/residualsTemplates.C b/src/functionObjects/utilities/residuals/residualsTemplates.C
index a221085..e1022bc 100644
--- a/src/functionObjects/utilities/residuals/residualsTemplates.C
+++ b/src/functionObjects/utilities/residuals/residualsTemplates.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2015-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2015-2017 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -87,6 +87,10 @@ void Foam::functionObjects::residuals::writeResidual(const word& fieldName)
                 }
             }
         }
+        else
+        {
+            file() << token::TAB << "'N/A'";
+        }
     }
 }
 
proposition_v1.patch (953 bytes)

wyldckat

2017-07-09 23:11

updater  

proposition_v1.tar.gz (1,188 bytes)

wyldckat

2017-07-09 23:15

updater   ~0008370

@Henry: Attached are the following files:

 - "proposition_v1.patch" - Shows the proposed change to write 'N/A' for when the entry is not found on the residuals list.

 - "proposition_v1.tar.gz" - Has the modified file "src/functionObjects/utilities/residuals/residualsTemplates.C" as indicated by the patch.

I tested it as according to the report and also used 'foamMonitor' as described in the User Guide in the reported section, which worked just fine, although it does give a warning message every time it refreshes, since it had to ignore the unusable column.

henry

2017-07-10 12:49

manager   ~0008375

Thanks for the patch Bruno, I applied it with only a small modification: I removed the single quotes from the 'N/A' to reduce clutter. Please reopen if you think the quotes are needed/useful for parsing of the file generated.

Resolved by 8b8e5b9492387a622eb28f0002529f312da3fc52

Issue History

Date Modified Username Field Change
2017-07-09 18:52 miri New Issue
2017-07-09 18:52 miri File Added: residuals.dat
2017-07-09 18:52 miri File Added: foamMonitor.png
2017-07-09 18:53 miri File Added: residuals
2017-07-09 18:53 miri File Added: controlDict
2017-07-09 18:56 miri File Added: log
2017-07-09 22:38 wyldckat Relationship added related to 0002547
2017-07-09 22:43 wyldckat Note Added: 0008369
2017-07-09 23:11 wyldckat File Added: proposition_v1.patch
2017-07-09 23:11 wyldckat File Added: proposition_v1.tar.gz
2017-07-09 23:15 wyldckat Note Added: 0008370
2017-07-09 23:15 wyldckat Assigned To => henry
2017-07-09 23:15 wyldckat Status new => assigned
2017-07-10 12:49 henry Status assigned => resolved
2017-07-10 12:49 henry Resolution open => fixed
2017-07-10 12:49 henry Note Added: 0008375