在用自己的构型适用于reactingFoam，反应物是C7H16和空气，目前修改了边界条件和反应机理等，已解决问题1
问题1：
--> FOAM FATAL ERROR:
Sum of mass fractions is zero for species
544
(
各种组分
）
From function void Foam::basicSpecieMixture::correctMassFractions()
in file mixtures/basicSpecieMixture/basicSpecieMixture.C at line 205.
FOAM exiting
解决方法：看0文件里有无以上涉及的组分，组分之和必须为1，我是修改后0文件里原本是C7H16，但调用的反应机理里面是NC7H16，更改了文件名后错误消失。
问题2：
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 9
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 9-c96968adf6d8
Exec : reactingFoam
Date : Apr 26 2023
Time : 08:29:07
Host : "yang-Kuangshi-X-Series-GM7AG7N"
PID : 12197
I/O : uncollated
Case : /home/yang/OpenFOAM/yang-9/run/smartinjection4.24
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
PIMPLE: No convergence criteria found
PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode
Reading thermophysical properties
Selecting thermodynamics package
{
type heRhoThermo;
mixture multiComponentMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}
Reading field U
Reading/calculating face flux field phi
Creating turbulence model.
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Prt 0.85;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}
Creating thermophysical transport model
Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
Creating reaction model
Combustion model not active: combustionProperties not found
Selecting combustion model none
Creating field dpdt
Creating field kinetic energy K
No MRF models present
No fvModels present
No fvConstraints present
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::nutkWallFunctionFvPatchScalarField::nut() const at ??:?
#4 Foam::nutWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#5 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#6 Foam::RASModels::kEpsilon<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correctNut() at ??:?
#7 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#8 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception (core dumped)
问题3：c7h16单步化学反应机理的reactingfoam
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = 23.14951573
Time = 4.3e-06
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab: Solving for Ux, Initial residual = 0.8932644369, Final residual = 0.01804921196, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.7925881184, Final residual = 0.04676651581, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.8082444789, Final residual = 0.07178631378, No Iterations 1
#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 Foam::combustionModels::EDC<Foam::psiReactionThermo>::correct() at ??:?
#5 ? in /opt/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/reactingFoam
#6 ? in /lib/x86_64-linux-gnu/libc.so.6
#7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8 ? in /opt/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/reactingFoam
Floating point exception (core dumped)
问题4：用了gri3.0的reactingfoam（燃烧模型EDM）
--> FOAM FATAL ERROR: (openfoam-2206)
Inconsistent thermo package for CombustionModel<psiReactionThermo> model:
hePsiThermo<reactingMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>
Please select a thermo package based on singleStepReactingMixture
From Foam::combustionModels::singleStepCombustion<ReactionThermo, ThermoType>::singleStepCombustion(const Foam::word&, ReactionThermo&, const Foam::compressibleTurbulenceModel&, const Foam::word&) [with ReactionThermo = Foam::psiReactionThermo; ThermoType = Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >]
in file lnInclude/singleStepCombustion.C at line 75.
FOAM exiting
问题5：用了gri3.0的reactingfoam（燃烧模型EDC）
Create time
Create mesh for time = 0
PIMPLE: Operating solver in PISO mode
Using LTS
Reading thermophysical properties
Selecting thermodynamics package
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}
Selecting chemistryReader chemkinReader
Reading field U
Reading/calculating face flux field phi
Creating turbulence model.
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Prt 0.85;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}
Creating reaction model
Selecting combustion model EDC
Selecting chemistry solver
{
solver ode;
method standard;
}
StandardChemistryModel: Number of species = 544 and reactions = 2446
Selecting ODE solver seulex
using integrated reaction rate
Creating field dpdt
Creating field kinetic energy K
No MRF models present
No finite volume options present
Starting time loop
Time scales min/max:
Flow = 9.931677386e-08, 0.0001
Temperature = 1e+300, 1e+300
Composition = 1e+300, 1e+300
Overall = 9.931677386e-08, 0.0001
Time = 1e-07
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 2.322751201e-16, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 8.928355209e-17, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 1.099301832e-16, No Iterations 1
Chemistry time solved max/min : 0 / 0
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 double Foam::StandardChemistryModel<Foam::psiReactionThermo, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) at ??:?
#4 Foam::combustionModels::EDC<Foam::psiReactionThermo>::correct() at ??:?
#5 ? in /opt/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/reactingFoam
#6 ? in /lib/x86_64-linux-gnu/libc.so.6
#7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8 ? in /opt/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/reactingFoam
Floating point exception (core dumped)
问题6：改七改八不知道改出了个啥玩意
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -348.7413834
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -357.1173073
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -642.2389819
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -654.7336172
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -586.8210643
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -598.5484188
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = 191.7211225
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = 190.9258018
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = 158.245464
--> FOAM Warning :
From Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /opt/OpenFOAM-v2206/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = 156.9752934
min/max(T) = 200, 5000
DICPCG: Solving for p, Initial residual = 0.2776863156, Final residual = 6.681837376e-07, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.998430897e-07, global = 8.558466776e-08, cumulative = 8.682383884e-06
DICPCG: Solving for p, Initial residual = 0.008298322922, Final residual = 7.419768043e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.185531302e-08, global = 2.538559952e-08, cumulative = 8.707769484e-06
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.01101147019, Final residual = 4.417456436e-07, No Iterations 1
bounding epsilon, min: -2.841173778e+29 max: 6.936160996e+30 average: 2.531930651e+26
DILUPBiCGStab: Solving for k, Initial residual = 1.944555359e-12, Final residual = 1.944555359e-12, No Iterations 0
ExecutionTime = 467.37 s ClockTime = 469 s
Time scales min/max:
Flow = 2.528261534e-12, 0.0001
Temperature = 7.245868382e+83, 1e+300
Composition = 1e+300, 1e+300
Overall = 2.528261534e-12, 5.165711841e-05
Time = 5.7e-06
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab: Solving for Ux, Initial residual = 0.912797583, Final residual = 0.03907631423, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.8698878836, Final residual = 0.0217991991, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.7918226439, Final residual = 0.05245962132, No Iterations 1
#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 Foam::combustionModels::EDC<Foam::psiReactionThermo>::correct() at ??:?
#5 ? in /opt/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/reactingFoam
#6 ? in /lib/x86_64-linux-gnu/libc.so.6
#7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8 ? in /opt/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/reactingFoam
Floating point exception (core dumped)
问题7：用of9已成功调用的reactingfoam改
Create time
Create mesh for time = 0
PIMPLE: No convergence criteria found
PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode
Using LTS
Reading thermophysical properties
Selecting thermodynamics package
{
type heRhoThermo;
mixture multiComponentMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}
Reading field U
Reading/calculating face flux field phi
Creating turbulence model.
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Prt 0.85;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}
Creating thermophysical transport model
Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
Creating reaction model
Selecting combustion model EDC
Selecting chemistry solver
{
solver ode;
method standard;
}
standardChemistryModel: Number of species = 544 and reactions = 2446
Selecting ODE solver seulex
Creating field dpdt
Creating field kinetic energy K
No MRF models present
No fvModels present
No fvConstraints present
Starting time loop
Time scales min/max:
Flow = 9.728569846e-08, 0.0001
Temperature = 4.494232837e+307, 4.494232837e+307
Composition = 4.494232837e+307, 4.494232837e+307
Overall = 9.728569846e-08, 0.0001
Time = 1e-07
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.0001092756986, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.006862564884, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.0001167898495, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::standardChemistryModel<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) at ??:?
#4 Foam::combustionModels::EDC::correct() at ??:?
#5 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#6 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception (core dumped)
扫一扫
1万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
