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

rhoReactingFoam error

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 20, 2015, 21:39
Default rhoReactingFoam error
  #1
New Member
 
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 7
shenzhou1987 is on a distinguished road
Hi foamers! I'm still working on my case which simulates the jet flow of a nozzle. Two species combustion gas and air should be modeled. So I choose rhoReactingFoam, because the jet flow is also supersoinc and compressible. But there is no reaction between gas and air, so I switch off the combustion model and reaction model in rhoReactingFoam. Now the case can be runned, but just after several iterations it always collapses. The error message is "Floating point exception (core dumped)". I have tried many times, using different fvschemes and linear solvers, but the problem is still here. So I post my case files in this thread, hope somebody can give me some advise. Thank you!
Ps: The model is axisymmetry. I have executed command 'checkMesh', and no error message is reported. But I still don't know whether the BCs of wedge plans are right.
BCs:
air
Code:
boundaryField
{
    wedge-front
    {
        type wedge;
    }
    wedge-back
    {
        type wedge;
    }
    pout-side
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }
    wall-penguan
    {
        type            zeroGradient;
    }
    wall-chamber
    {
        type            zeroGradient;
    }
    pin-highp
    {
        type            fixedValue;
        value           uniform 0;
    }
}
gas
Code:
boundaryField
{
    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    wall-penguan
    {
        type            zeroGradient;
    }
    wall-chamber
    {
        type            zeroGradient;
    }
    pin-highp
    {
        type            fixedValue;
        value           uniform 1;
    }
}
p
Code:
boundaryField
{
    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            inletOutlet;
        inletValue      uniform 101325;
        value           uniform 101325;
    }
    wall-penguan
    {
        type            zeroGradient;
    }
    wall-chamber
    {
        type            zeroGradient;
    }
    pin-highp
    {
        type            totalPressure;
        gamma           1.2;
        p0              uniform 1e6;
    }
}
T
Code:
boundaryField
{

    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }
    wall-penguan
    {
        type            zeroGradient;
    }
    wall-chamber
    {
        type             zeroGradient;
    }
    pin-highp
    {
        type            fixedValue;
        value           uniform 3500;
    }
}
U
Code:
boundaryField
{
    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            zeroGradient;
    }   
    wall-penguan
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    wall-chamber
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    pin-highp
    {
        type            fixedValue;
        value           uniform (10 0 0);
    }
}
k
Code:
boundaryField
{
    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }
    wall-penguan
    {
        type            compressible::kqRWallFunction;
        value           uniform 1;
    }
    wall-chamber
    {
        type            compressible::kqRWallFunction;
        value           uniform 1;
    }
    pin-highp
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;
        value           uniform 1;
    }
}
epsilon
Code:
boundaryField
{
    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            inletOutlet;
        inletValue      uniform 200;
        value           uniform 200;
    }
    wall-penguan
    {
         type            compressible::epsilonWallFunction;
         value           uniform 200;
    }
    wall-chamber
    {
         type            compressible::epsilonWallFunction;
         value           uniform 200;
    }
    pin-highp
    {
         type            compressible::turbulentMixingLengthDissipationRateInlet;
         mixingLength    0.005;
         value           uniform 200;
    }
}
mut
Code:
boundaryField
{
    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            calculated;
        value           uniform 0;
    }
    wall-penguan
    {
        type            mutkWallFunction;
        value           uniform 0;
    }
    wall-chamber
    {
        type            mutkWallFunction;
        value           uniform 0;
    }
    pin-highp
    {
        type            calculated;
        value           uniform 0;
    }
}
alphat
Code:
boundaryField
{
    wedge-front
    {
        type            wedge;
    }
    wedge-back
    {
        type            wedge;
    }
    pout-side
    {
        type            calculated;
        value           uniform 0;
    }
    wall-penguan
    {
        type            compressible::alphatWallFunction;
        value           uniform 0;
    }
    wall-chamber
    {
        type            compressible::alphatWallFunction;
        value           uniform 0;
    }
    pin-highp
    {
        type            calculated;
        value           uniform 0;
    }
}
shenzhou1987 is offline   Reply With Quote

Old   January 20, 2015, 21:50
Default the constant files
  #2
New Member
 
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 7
shenzhou1987 is on a distinguished road
thermophysicalProperties
Code:
thermoType
{
    type            heRhoThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}
inertSpecie    air;
chemistryReader foamChemistryReader;
foamChemistryFile "./constant/reactions";
foamChemistryThermoFile "./constant/thermo.compressibleGas";
reactions
Code:
species
(
    air
    gas
);
reactions
{
}
thermo.compressibleGas
Code:
air
{
    specie
    {
        nMoles      1;
        molWeight   28.9;
    }
    thermodynamics
    {
        Tlow         200;
        Thigh        5000;
        Tcommon      300;
        highCpCoeffs (3.49344 0 0 0 0 0 0);
        lowCpCoeffs  (3.49344 0 0 0 0 0 0);
    }
    transport
    {
        As          1.4792e-06;
        Ts          116;
    }
}
gas
{
    specie
    {
        nMoles       1;
        molWeight    27.7133;
    }
    thermodynamics
    {
        Tlow         200;
        Thigh        5000;
        Tcommon      300;
        highCpCoeffs (6 0 0 0 0 0 0);
        lowCpCoeffs  (6 0 0 0 0 0 0);
    }
    transport
    {
        As           1.5e-6;
        Ts           120;
    }
}
chemistryProperties
Code:
chemistry           off;
combustionProperties
Code:
combustionModel  noCombustion<rhoThermoCombustion>;
active  false;
noCombustionCoeffs
{
}
turbulenceProperties
Code:
simulationType  RASModel;
RASProperties
Code:
RASModel        RNGkEpsilon;
turbulence      on;
printCoeffs     on;
shenzhou1987 is offline   Reply With Quote

Old   January 20, 2015, 21:56
Default the system files
  #3
New Member
 
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 7
shenzhou1987 is on a distinguished road
controlDict
Code:
application     rhoReactingFoam;
startFrom       startTime;
startTime       0;
stopAt          endTime;
endTime         1e-3;
deltaT          1e-5;
writeControl    timeStep;
writeInterval   5;
purgeWrite      0;
writeFormat     ascii;
writePrecision  8;
writeCompression off;
timeFormat      general;
timePrecision   6;
graphFormat     raw;
runTimeModifiable true;
adjustTimeStep yes;
maxCo 0.75;
fvSchemes
Code:
ddtSchemes
{
    default         Euler;
}
gradSchemes
{
    default         Gauss linear;
}
divSchemes
{
    default         none;

    div(phi,U)                       Gauss upwind;
    div(phi,Yi_h)                    Gauss upwind;
    div(phi,K)                       Gauss linear;
    div(phid,p)                      Gauss upwind;
    div(phi,epsilon)                 Gauss upwind;
    div(phi,k)                       Gauss upwind;
    div((muEff*dev2(T(grad(U)))))    Gauss linear;
}
laplacianSchemes
{
    default         Gauss linear corrected;
}
interpolationSchemes
{
    default         linear;
}
snGradSchemes
{
    default         corrected;
}
fluxRequired
{
    default         no;
    p;
}
fvSolution
Code:
solvers
{
    p
    {
        solver           PCG;
        preconditioner   DIC;
        tolerance        1e-6;
        relTol           0.0;
        maxIter          100;
    }

    pFinal
    {
        $p;
        tolerance        1e-6;
        relTol           0.0;
    }

    "rho.*"
    {
        $p;
    }

    "(U|h|k|epsilon)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0.0;
        maxIter         100;
    }

    "(U|h|k|epsilon)Final"
    {
        $U;
        relTol          0;
    }

    Yi
    {
        $hFinal;
    }
}
PIMPLE
{
    nCorrectors                 2;
    nOuterCorrectors            50;
    nNonOrthogonalCorrectors    0;

    residualControl
    {
        "(p|U|k|epsilon)"
        {
            tolerance    1e-6;
            relTol       0.0;
        }
        "(h|Yi)"
        {
            tolerance    1e-6;
            relTol       0.0;
        }
    }
}
relaxationFactors
{
    fields
    {
        p         0.3;
        pFinal    1;
    }
    equations
    {
        U                             0.7;
        k                             0.25;
        epsilon                       0.25;
        rho                           0.8;
        h                             0.8;
        Yi                            0.8;
        "(U|k|epsilon|rho|h|Yi)Final"     1;
    }
}
fvOptions
Code:
temperatureConstraints
{
    type             temperatureLimitsConstraint;
    active           true;
    selectionMode    all;

    temperatureLimitsConstraintCoeffs
    {
        Tmin    273;
        Tmax    3500;
    }
}
shenzhou1987 is offline   Reply With Quote

Old   January 20, 2015, 22:19
Default
  #4
New Member
 
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 7
shenzhou1987 is on a distinguished road
Otherwise, in the file "thermophysicalProperties" I attempted to use "multiComponentMixture" as the mixture model. But it always complained about the definition of species. If I use "reactingMixture", no error is reported. Here is the original file. What's wrong with this file?
thermophysicalProperties
Code:
thermoType
{
    type            hePsiThermo;
    mixture         multiComponentMixture;
    transport       sutherland;
    thermo          janaf;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

species
{
    gas
    air
};

inertSpecies    air;

gas
{
    specie
    {
        nMoles       1;
        molWeight    27.7133;
    }
    thermodynamics
    {
        Tlow         200;
        Thigh        5000;
        Tcommon      300;
        highCpCoeffs (6 0 0 0 0 0 0);
        lowCpCoeffs  (6 0 0 0 0 0 0);
    }
    transport
    {
        As           1.5e-6;
        Ts           120;
    }
}

air
{
    specie
    {
        nMoles      1;
        molWeight   28.9;
    }
    thermodynamics
    {
        Tlow         200;
        Thigh        5000;
        Tcommon      300;
        highCpCoeffs (3.49344 0 0 0 0 0 0);
        lowCpCoeffs  (3.49344 0 0 0 0 0 0);
    }
    transport
    {
        As          1.4792e-06;
        Ts          116;
    }
}
shenzhou1987 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


Similar Threads
Thread Thread Starter Forum Replies Last Post
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries NickG OpenFOAM Installation 3 December 30, 2019 00:21
[blockMesh] blockMesh with double grading. spwater OpenFOAM Meshing & Mesh Conversion 92 January 12, 2019 09:00
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM Community Contributions 300 October 29, 2014 18:00
OpenFOAM without MPI kokizzu OpenFOAM Installation 4 May 26, 2014 09:17
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51


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