[rezvani]

Hello. i want use openfoam to solve laminar viscos fluid with rotating wall. my problem solved in fluent by turn off viscous heating and use simple scheme for Coupling scheme & use heat source in rotating wall instead of viscus heating. but in openfoam i get error of inital temprarure.
this is my setting:

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(U)         cellLimited Gauss linear 1;
}

divSchemes
{
div(phi,alpha)     Gauss interfaceCompression vanLeer 1;
div((((alpha.air*thermo:rho.air)*nuEff.air)*dev2(T(grad(U)))))    Gauss linear;
div((((alpha.polymer*thermo:rho.polymer)*(thermo:mu.polymer|thermo:rho.polymer))*dev2(T(grad(U)))))   Gauss linear;
div(((alpha.polymer*thermo:rho.polymer)*sigma.polymer))    Gauss linear;
div((((alpha.polymer*thermo:rho.polymer)*nuM)*grad(U)))  Gauss linear;
div(rhoPhi,U) Gauss linearUpwindV grad(U);
div(rhoPhi,K)     Gauss upwind;
div((phi+meshPhi),p)   Gauss upwind;
div(rhoPhi,T)  Gauss linearUpwind grad(T);
div(alphaRhoPhi.polymer,sigma.polymer)   Gauss linearUpwind grad(sigma);
div(((rho*nuEff)*dev2(T(grad(U)))))   Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited corrected 0.33;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited corrected 0.33;
}

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       WLF;
    thermo           hConst;
    equationOfState  rPolynomial;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        molWeight   291.3;
    }
    equationOfState
    {
        C (0.005218 -3.1055e-06 4.8689e-09 5.7772e-13 -3.2225e-16);
    }
    thermodynamics
    {
      Cp            1800;
        Hf          0;
    }
    transport
    {
        mu0         53763179097.64641453;
        Tr          115;
        C1          31.24;
        C2          107.53;
        Pr        6.15;
    }
}

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 propellerStem

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       WLF;
    thermo          hConst;
    equationOfState rPolynomial;
    specie          specie;
    energy          sensibleEnthalpy;
}

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 1131591 cell(s) with volume 0.0005971852225
Courant Number mean: 0.001553287087 max: 0.06698099517
Constructing face velocity Uf


Starting time loop

Courant Number mean: 0.001553287087 max: 0.06698099517
Interface Courant Number mean: 0 max: 0
deltaT = 0.001199040767
Time = 0.00119904

PIMPLE: Iteration 1
smoothSolver:  Solving for alpha.polymer, Initial residual = 4.148307397e-07, Final residual = 4.148307397e-07, No Iterations 0
Phase-1 volume fraction = 0.2935706537  Min(alpha.polymer) = 0  Max(alpha.polymer) = 1
MULES: Correcting alpha.polymer
MULES: Correcting alpha.polymer
Phase-1 volume fraction = 0.2935706537  Min(alpha.polymer) = 0  Max(alpha.polymer) = 1
DILUPBiCGStab:  Solving for T, Initial residual = 0.7538694576, Final residual = 7.468251956e-07, No Iterations 59
DICPCG:  Solving for p_rgh, Initial residual = 0.9999999577, Final residual = 14.02207328, No Iterations 1000
PIMPLE: Iteration 2
smoothSolver:  Solving for alpha.polymer, Initial residual = 0.0001202851214, Final residual = 9.627369821e-07, No Iterations 18
Phase-1 volume fraction = 0.2935690002  Min(alpha.polymer) = 0  Max(alpha.polymer) = 1
MULES: Correcting alpha.polymer
MULES: Correcting alpha.polymer
Phase-1 volume fraction = 0.2935690002  Min(alpha.polymer) = -0.1114075347  Max(alpha.polymer) = 1.245386626
DILUPBiCGStab:  Solving for T, Initial residual = 0.2150890754, Final residual = 3.474027855e-07, No Iterations 169
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4  exp in "/lib/x86_64-linux-gnu/libm.so.6"
#5  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::WLFTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rPolynomial<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#6  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::WLFTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rPolynomial<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#7  Foam::twoPhaseMixtureThermo::correctThermo() at ??:?
#8  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
Floating point exception (core dumped)

thank you.
best regards.

[galoisati]

Hi,



It could be a bit late but I had a similar problem and solved it with below change in the thermophysicalProperties dictionary:



thermoType
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}


changed to:



thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}