CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Error in compressibleInterFoam with Viscoelasticity and dynamic mesh (https://www.cfd-online.com/Forums/openfoam-solving/232461-error-compressibleinterfoam-viscoelasticity-dynamic-mesh.html)

rezvani December 15, 2020 04:23
Error in compressibleInterFoam with Viscoelasticity and dynamic mesh
 
5 Attachment(s)
hello my model is rotating pin in 2 phase (1 phase is non isothermal & Viscoelastic) in compressibleInterFoam solver. and i get this error:
Code:
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: solidBody
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone innerCylinder

PIMPLE: No convergence criteria found


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


PIMPLE: Operating solver in transient mode with 2 outer correctors


Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture        pureMixture;
    transport      const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture        pureMixture;
    transport      const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading thermophysical properties


Reading g

Reading hRef
Calculating field g.h

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating field kinetic energy K

No MRF models present

Creating finite volume options from "system/fvOptions"

Selecting finite volume options model type limitTemperature
    Source: limitT
    - selecting all cells
    - selected 1044384 cell(s) with volume 0.0003577383686
Courant Number mean: 0.001733236415 max: 0.1129887195
Constructing face velocity Uf


Starting time loop

Courant Number mean: 0.001733236415 max: 0.1129887195
Interface Courant Number mean: 0 max: 0
deltaT = 0.001199040767
Time = 0.00119904

PIMPLE: Iteration 1
smoothSolver:  Solving for alpha.polymer, Initial residual = 3.79535123e-06, Final residual = 5.564230094e-12, No Iterations 1
Phase-1 volume fraction = 0.3016640131  Min(alpha.polymer) = 0  Max(alpha.polymer) = 1
MULES: Correcting alpha.polymer
MULES: Correcting alpha.polymer
Phase-1 volume fraction = 0.3016640131  Min(alpha.polymer) = -1.12080956e-10  Max(alpha.polymer) = 1
DILUPBiCGStab:  Solving for T, Initial residual = 1, Final residual = 9.759238949e-07, No Iterations 164


--> FOAM FATAL ERROR:
Negative initial temperature T0: -3294.499179

    From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const, bool) const [with Thermo = Foam::hConstThermo<Foam::rhoConst<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-8/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 56.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double) const, bool) const at ??:?
#3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5  Foam::twoPhaseMixtureThermo::correctThermo() at ??:?
#6  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
Aborted (core dumped)
alpha.polymer
Code:
dimensions      [1 -1 -2 0 0 0 0];

internalField  uniform 1e5;

boundaryField
{

    #includeEtc "caseDicts/setConstraintTypes"

polyinlet
    {
        type            calculated;
        value          $internalField;
    }
    outlet
    {
          type            calculated;
        value          $internalField;
    }

    vessel
    {
        type            calculated;
      value          $internalField;
    }

    propeller
    {
        type            calculated;
        value          $internalField;
    }
    allBoundary
      {
        type            calculated;
        value          $internalField;
    }
}
p
Code:
dimensions      [1 -1 -2 0 0 0 0];

internalField  uniform 1e5;

boundaryField
{

    #includeEtc "caseDicts/setConstraintTypes"

polyinlet
    {
        type            calculated;
        value          $internalField;
    }
    outlet
    {
          type            calculated;
        value          $internalField;
    }

    vessel
    {
        type            calculated;
      value          $internalField;
    }

    propeller
    {
        type            calculated;
        value          $internalField;
    }
    allBoundary
      {
        type            calculated;
        value          $internalField;
    }
}
p_rgh
Code:
dimensions      [1 -1 -2 0 0];

internalField  uniform 1e5;

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

polyinlet
    {
        type            fixedValue;
        value          $internalField;
    }
    outlet
    {
        type            totalPressure;
        p0              $internalField;
    }

    vessel
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value          uniform 1e5;
  }

    propeller
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value          uniform 1e5;
    }
    allBoundary
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value          uniform 1e5;
    }
}
sigma.polymer
Code:
dimensions      [0 2 -2 0 0 0 0];

internalField  uniform (0 0 0 0 0 0);

boundaryField
{

    #includeEtc "caseDicts/setConstraintTypes"

polyinlet
    {
        type            inletOutlet;
        inletValue      $internalField;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
    }

  vessel
    {
        type            zeroGradient;
  }

        propeller
    {
        type            zeroGradient;
    }
    allBoundary
    {
        type            zeroGradient;
    }

}
T
Code:
dimensions      [0 0 0 1 0 0 0];

internalField  uniform 360;

boundaryField
{

    #includeEtc "caseDicts/setConstraintTypes"
        polyinlet
    {
        type            inletOutlet;
        inletValue      $internalField;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
    }

    vessel
    {
        type            zeroGradient;
    }

    "propeller.*"
    {
        type            zeroGradient;
    }
    allBoundary
    {
        type            zeroGradient;
    }

}
T.air
Code:
dimensions      [0 0 0 1 0 0 0];

internalField  uniform 360;

boundaryField
{
    propeller
    {
        type            calculated;
        value          uniform 360;
    }
    vessel
    {
        type            calculated;
        value          uniform 360;
    }
    polyinlet
    {
        type            calculated;
        value          uniform 360;
    }
    outlet
    {
        type            calculated;
        value          uniform 360;
    }
}
T.polymer
Code:
dimensions      [0 0 0 1 0 0 0];

internalField  uniform 360;

boundaryField
{
    propeller
    {
        type            calculated;
        value          uniform 360;
    }
    vessel
    {
        type            calculated;
        value          uniform 360;
    }
    polyinlet
    {
        type            calculated;
        value          uniform 360;
    }
    outlet
    {
        type            calculated;
        value          uniform 360;
    }
}
U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField  uniform (0.0025 0 0);

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

    polyinlet
    {
        type            uniformFixedValue;
        uniformValue    constant;
        value          (0.0025 0 0);
    }
 
    outlet
    {
        type            pressureInletOutletVelocity;
        value            $internalField;
    }
 
    vessel
    {
        type            fixedValue;
        value          $internalField;
  }

    "propeller.*"
    {
        type            movingWallVelocity;
        value          uniform (0 0 0);
    }
    allBoundary
        {
        type            fixedValue;
        value          $internalField;
    }

}

please help me.
thank you.

CFD-HSNR December 16, 2020 09:21
Quote:
Originally Posted by rezvani (Post 790693)
hello my model is rotating pin in 2 phase (1 phase is non isothermal & Viscoelastic) in compressibleInterFoam solver. and i get this error:
Code:
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: solidBody
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone innerCylinder

PIMPLE: No convergence criteria found


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


PIMPLE: Operating solver in transient mode with 2 outer correctors


Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture        pureMixture;
    transport      const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture        pureMixture;
    transport      const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading thermophysical properties


Reading g

Reading hRef
Calculating field g.h

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating field kinetic energy K

No MRF models present

Creating finite volume options from "system/fvOptions"

Selecting finite volume options model type limitTemperature
    Source: limitT
    - selecting all cells
    - selected 1044384 cell(s) with volume 0.0003577383686
Courant Number mean: 0.001733236415 max: 0.1129887195
Constructing face velocity Uf


Starting time loop

Courant Number mean: 0.001733236415 max: 0.1129887195
Interface Courant Number mean: 0 max: 0
deltaT = 0.001199040767
Time = 0.00119904

PIMPLE: Iteration 1
smoothSolver:  Solving for alpha.polymer, Initial residual = 3.79535123e-06, Final residual = 5.564230094e-12, No Iterations 1
Phase-1 volume fraction = 0.3016640131  Min(alpha.polymer) = 0  Max(alpha.polymer) = 1
MULES: Correcting alpha.polymer
MULES: Correcting alpha.polymer
Phase-1 volume fraction = 0.3016640131  Min(alpha.polymer) = -1.12080956e-10  Max(alpha.polymer) = 1
DILUPBiCGStab:  Solving for T, Initial residual = 1, Final residual = 9.759238949e-07, No Iterations 164


--> FOAM FATAL ERROR:
Negative initial temperature T0: -3294.499179

    From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const, bool) const [with Thermo = Foam::hConstThermo<Foam::rhoConst<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-8/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 56.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double) const, bool) const at ??:?
#3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5  Foam::twoPhaseMixtureThermo::correctThermo() at ??:?
#6  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
Aborted (core dumped)
please help me.
thank you.

Hello,

As you can read from the error message, your initial temperature is responsible for the crash of your simulation:

Quote:
--> FOAM FATAL ERROR:
Negative initial temperature T0: -3294.499179
I would recommend you to check your boundary conditions. The files you attached do not allow us to identify the problem.

rezvani December 18, 2020 07:59
1 Attachment(s)
Quote:
Originally Posted by CFD-HSNR (Post 790832)
Hello,

As you can read from the error message, your initial temperature is responsible for the crash of your simulation:



I would recommend you to check your boundary conditions. The files you attached do not allow us to identify the problem.
thank you. i try to chenge boundary condition but same error happened. please help.my case file is


All times are GMT -4. The time now is 11:46.