 |
 |
 |
|
|
|
[Sponsors] | ||||
January 20, 2015, 21:39
|
rhoReactingFoam error
|
#1 |
|
New Member
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 11  |
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;
}
}
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;
}
}
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;
}
}
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;
}
}
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);
}
}
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;
}
}
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;
}
}
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;
}
}
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;
}
}
|
|
|
|
|
|
January 20, 2015, 21:50
|
the constant files
|
#2 |
|
New Member
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 11 |
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";
Code:
species
(
air
gas
);
reactions
{
}
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;
}
}
Code:
chemistry off; Code:
combustionModel noCombustion<rhoThermoCombustion>;
active false;
noCombustionCoeffs
{
}
Code:
simulationType RASModel; Code:
RASModel RNGkEpsilon; turbulence on; printCoeffs on; |
|
|
|
|
|
|
January 20, 2015, 21:56
|
the system files
|
#3 |
|
New Member
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 11 |
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; 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;
}
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;
}
}
Code:
temperatureConstraints
{
type temperatureLimitsConstraint;
active true;
selectionMode all;
temperatureLimitsConstraintCoeffs
{
Tmin 273;
Tmax 3500;
}
}
|
|
|
|
|
|
|
January 20, 2015, 22:19
|
|
#4 |
|
New Member
Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 11 |
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;
}
}
|
|
|
|
|
|
 |
|
|
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 |
Linear Mode
