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

Error in compressibleInterFoam with Viscoelasticity and dynamic mesh

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

Like Tree1Likes
  • 1 Post By CFD-HSNR

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 15, 2020, 05:23
Default Error in compressibleInterFoam with Viscoelasticity and dynamic mesh
  #1
New Member
 
yahya rezvani
Join Date: Aug 2019
Posts: 17
Rep Power: 4
rezvani is on a distinguished road
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.
Attached Images
File Type: png porpeller1.png (100.1 KB, 5 views)
Attached Files
File Type: txt snappyHexMesh.txt (11.1 KB, 1 views)
File Type: txt sigmapolymer.txt (1.4 KB, 1 views)
File Type: txt fvSchemes.txt (1.7 KB, 1 views)
File Type: txt fvSolution.txt (1.8 KB, 2 views)

Last edited by rezvani; December 18, 2020 at 12:03.
rezvani is offline   Reply With Quote

Old   December 16, 2020, 10:21
Default
  #2
New Member
 
Max
Join Date: Apr 2020
Location: Germany
Posts: 8
Rep Power: 3
CFD-HSNR is on a distinguished road
Quote:
Originally Posted by rezvani View Post
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 likes this.
CFD-HSNR is offline   Reply With Quote

Old   December 18, 2020, 08:59
Default
  #3
New Member
 
yahya rezvani
Join Date: Aug 2019
Posts: 17
Rep Power: 4
rezvani is on a distinguished road
Quote:
Originally Posted by CFD-HSNR View Post
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
Attached Files
File Type: zip viscoelasticmixer.zip (24.9 KB, 2 views)
rezvani is offline   Reply With Quote

Reply

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



All times are GMT -4. The time now is 07:18.