==76523== Memcheck, a memory error detector
==76523== Copyright (C) 2002-2017, and GNU GPL'd, by Julian Seward et al.
==76523== Using Valgrind-3.15.0 and LibVEX; rerun with -h for copyright info
==76523== Command: reactingTwoPhaseEulerFoam
==76523==
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 7
Exec   : reactingTwoPhaseEulerFoam
Date   : Dec 19 2019
Time   : 18:34:05
Host   : "..."
PID    : 76523
I/O    : uncollated
Case   : ...
nProcs : 1
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: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calculations will do 3 corrections


PIMPLE: Operating solver in transient mode with 3 outer correctors



Reading g

Reading hRef
Creating phaseSystem

Selecting twoPhaseSystem basicTwoPhaseSystem
Selecting phaseModel for air: purePhaseModel
Selecting diameterModel for phase air: isothermal
Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.air
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Selecting phaseModel for water: purePhaseModel
Selecting diameterModel for phase water: constant
Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          eConst;
    equationOfState perfectFluid;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.water
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
No MRF models present

Selecting default blending method: linear
Selecting drag blending method: linear
Selecting surfaceTensionModel for (air and water): constant
Selecting aspectRatioModel for (air in water): constant
Selecting aspectRatioModel for (water in air): constant
Selecting dragModel for (air and water): segregated
Selecting swarmCorrection for (air and water): none
Selecting dragModel for (air in water): SchillerNaumann
Selecting swarmCorrection for (air in water): none
Selecting dragModel for (water in air): SchillerNaumann
Selecting swarmCorrection for (water in air): none
Selecting virtualMassModel for (air in water): constantCoefficient
Selecting virtualMassModel for (water in air): constantCoefficient
Selecting heatTransferModel for (air in water): RanzMarshall
Selecting heatTransferModel for (water in air): RanzMarshall
Calculating field g.h

Reading field p_rgh

Courant Number mean: 0.011125 max: 0.037

Starting time loop

fieldAverage fieldAverage1:
    Starting averaging at time 0

Courant Number mean: 0.011125 max: 0.037
Max Ur Courant Number = 0.037
Time = 0.005

PIMPLE: Iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.297297  Min(alpha1) = 0  Max(alpha1) = 1
Constructing momentum equations
==76523== Invalid read of size 16
==76523==    at 0x9AF75F4: Foam::sqr(Foam::Field<Foam::SymmTensor<double> >&, Foam::UList<Foam::Vector<double> > const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so)
==76523==    by 0x9AF796C: Foam::sqr(Foam::tmp<Foam::Field<Foam::Vector<double> > > const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so)
==76523==    by 0x6F2FDB9: Foam::pressureInletOutletVelocityFvPatchVectorField::updateCoeffs() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x4645C9: Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==    by 0x43A965: main (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==  Address 0x146fab48 is 568 bytes inside a block of size 576 alloc'd
==76523==    at 0x4C31C63: operator new[](unsigned long) (vg_replace_malloc.c:433)
==76523==    by 0x6C8EDF4: Foam::tmp<Foam::Field<Foam::Vector<double> > > Foam::operator/<Foam::Vector<double> >(Foam::UList<Foam::Vector<double> > const&, Foam::UList<double> const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6C8E44C: Foam::fvPatch::nf() const (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6F2FDAE: Foam::pressureInletOutletVelocityFvPatchVectorField::updateCoeffs() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x4645C9: Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==    by 0x43A965: main (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==
smoothSolver:  Solving for e.air, Initial residual = 1, Final residual = 6.61277e-09, No Iterations 2
smoothSolver:  Solving for e.water, Initial residual = 1, Final residual = 9.84031e-11, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 0.00018724, Final residual = 2.39268e-09, No Iterations 10
==76523== Invalid read of size 16
==76523==    at 0x9AF75F4: Foam::sqr(Foam::Field<Foam::SymmTensor<double> >&, Foam::UList<Foam::Vector<double> > const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so)
==76523==    by 0x9AF796C: Foam::sqr(Foam::tmp<Foam::Field<Foam::Vector<double> > > const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so)
==76523==    by 0x6F2FDB9: Foam::pressureInletOutletVelocityFvPatchVectorField::updateCoeffs() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6D34D01: Foam::directionMixedFvPatchField<Foam::Vector<double> >::evaluate(Foam::UPstream::commsTypes) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x45653A: Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==    by 0x43E231: main (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==  Address 0x14016598 is 568 bytes inside a block of size 576 alloc'd
==76523==    at 0x4C31C63: operator new[](unsigned long) (vg_replace_malloc.c:433)
==76523==    by 0x6C8EDF4: Foam::tmp<Foam::Field<Foam::Vector<double> > > Foam::operator/<Foam::Vector<double> >(Foam::UList<Foam::Vector<double> > const&, Foam::UList<double> const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6C8E44C: Foam::fvPatch::nf() const (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6F2FDAE: Foam::pressureInletOutletVelocityFvPatchVectorField::updateCoeffs() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6D34D01: Foam::directionMixedFvPatchField<Foam::Vector<double> >::evaluate(Foam::UPstream::commsTypes) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x45653A: Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==    by 0x43E231: main (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==
==76523== Invalid read of size 16
==76523==    at 0x9AF75F4: Foam::sqr(Foam::Field<Foam::SymmTensor<double> >&, Foam::UList<Foam::Vector<double> > const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so)
==76523==    by 0x9AF796C: Foam::sqr(Foam::tmp<Foam::Field<Foam::Vector<double> > > const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so)
==76523==    by 0x6F2FDB9: Foam::pressureInletOutletVelocityFvPatchVectorField::updateCoeffs() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6D34CFC: Foam::directionMixedFvPatchField<Foam::Vector<double> >::evaluate(Foam::UPstream::commsTypes) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x45653A: Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==    by 0x43E231: main (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==  Address 0x14017a28 is 568 bytes inside a block of size 576 alloc'd
==76523==    at 0x4C31C63: operator new[](unsigned long) (vg_replace_malloc.c:433)
==76523==    by 0x6C8EDF4: Foam::tmp<Foam::Field<Foam::Vector<double> > > Foam::operator/<Foam::Vector<double> >(Foam::UList<Foam::Vector<double> > const&, Foam::UList<double> const&) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6C8E44C: Foam::fvPatch::nf() const (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6F2FDAE: Foam::pressureInletOutletVelocityFvPatchVectorField::updateCoeffs() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x6D34CFC: Foam::directionMixedFvPatchField<Foam::Vector<double> >::evaluate(Foam::UPstream::commsTypes) (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so)
==76523==    by 0x45653A: Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==    by 0x43E231: main (in /software/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam)
==76523==
PIMPLE: Iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.297302  Min(alpha1) = 0  Max(alpha1) = 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 0.0294606, Final residual = 3.72835e-09, No Iterations 3
smoothSolver:  Solving for e.water, Initial residual = 5.07225e-05, Final residual = 1.8718e-10, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.97072e-05, Final residual = 5.83036e-09, No Iterations 7
PIMPLE: Iteration 3
MULES: Solving for alpha.air
MULES: Solving for alpha.air

...