Jet Engine case, initiating a detonation wave, and reactingParcelFoam errors

HRobertHS · March 11, 2017, 05:05

Hi OpenFOAM community,

I am trying to simulate a rotational detonation engine, except I am getting many errors and need help with starting the detonation wave.

CheckMesh output:

Code:

checkMesh
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : checkMesh
Date   : Mar 11 2017
Time   : 00:42:57
Host   : "robert-Lenovo-G580"
PID    : 7901
Case   : /home/robert/RDESim
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           44976
    faces:            488140
    internal faces:   474684
    cells:            240706
    faces per cell:   4
    boundary patches: 1
    point zones:      0
    face zones:       2
    cell zones:       3

Overall number of cells of each type:
    hexahedra:     0
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    240706
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    defaultFaces        13456    6730     ok (closed singly connected)      

Checking geometry...
    Overall domain bounding box (-65 -64.99995 0) (65 64.99994 80)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (-4.958124e-18 -2.418597e-19 2.611782e-16) OK.
    Max cell openness = 4.504071e-16 OK.
    Max aspect ratio = 20.56824 OK.
    Minimum face area = 0.1640221. Maximum face area = 131.0312.  Face area magnitudes OK.
    Min volume = 0.04628868. Max volume = 199.7307.  Total volume = 1052052.  Cell volumes OK.
    Mesh non-orthogonality Max: 80.69549 average: 19.35907
   *Number of severely non-orthogonal (> 70 degrees) faces: 141.
    Non-orthogonality check OK.
  <<Writing 141 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 1.820535 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

Paraview:

Error #1:

Code:

 DILUPBiCG:  Solving for Ux, Initial residual = 0.01593041, Final residual = 3.321646e-05, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.01461341, Final residual = 1.355258e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.03112335, Final residual = 1.287348e-06, 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  Foam::EulerImplicit<Foam::chemistryModel<Foam::rhoChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::solve(Foam::Field<double>&, double&, double&, double&, double&) const at ??:?
#4  double Foam::chemistryModel<Foam::rhoChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::solve<Foam::UniformField<double> >(Foam::UniformField<double> const&) at ??:?
#5  Foam::chemistryModel<Foam::rhoChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::solve(double) at ??:?
#6  Foam::combustionModels::laminar<Foam::combustionModels::rhoChemistryCombustion>::correct() at ??:?
#7  Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion>::correct() at ??:?
#8  ? at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? at ??:?

Warning #1:

Code:

FOAM Warning : 
 From function void Foam::CentredFitSnGradData<Polynomial>::calcFit(Foam::scalarList&, const Foam::List<Foam::Vector<double> >&, Foam::scalar, Foam::scalar, Foam::label) [with Polynomial = Foam::linearFitPolynomial; Foam::scalarList = Foam::List<double>; Foam::scalar = double; Foam::label = int]
 in file lnInclude/CentredFitSnGradData.C at line 158

Floating point exception (core dumped)

Code:

Warning #2: smoothSolver:  Solving for Ux, Initial residual = 0.03141642, Final residual = 3.558397e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.01473086, Final residual = 1.896657e-09, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.0291702, Final residual = 4.048035e-09, No Iterations 1
--> FOAM Warning : 
    From function Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::RASModels::SpalartAllmaras<BasicTurbulenceModel>::epsilon() const [with BasicTurbulenceModel = Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > >]
    in file ../turbulenceModels/lnInclude/SpalartAllmaras.C at line 325
    Turbulence kinetic energy dissipation rate not defined for Spalart-Allmaras model. Returning zero field

Warning #3:

Code:

--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 105
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 6000;  T = 189.5443
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 105
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 6000;  T = 188.6775
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 105
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 6000;  T = 196.8531
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 105
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 6000;  T = 196.468

Also, I have two injectors, but neither work in injecting parcels. How am I supposed to create a subpatch for a detonation wave. I have looked online, but could not find any adequate explanations for myself. Finally, I was wondering how I can keep track of the thrust and mass flow rate to calculate specific impulse. I have a forces file in postprocessing and thought that that was needed. Any help would be appreciated, and I need to get this done very soon! The files can be downloaded here and is a relatively large case: https://www.dropbox.com/s/4kh5wkdsde...im.tar.gz?dl=0

fvSolutions: https://www.dropbox.com/s/93i2qx11ae...vSolution?dl=0
fvSchemes: https://www.dropbox.com/s/volce020cz...fvSchemes?dl=0
Sincerely,
HRobertHS

HRobertHS · March 11, 2017, 15:49

I tried switching my thermoType model to

Code:

thermoType
{
	type	hePsiThermo;
	mixture	reactingMixture;
	transport	sutherland;
	thermo	econst;
	equationOfState	perfectGas;
	specie	specie;
	energy	sensibleInternalEnergy;
}

from

Code:

thermoType
{
	type	heRhoThermo;
	mixture	reactingMixture;
	transport	sutherland;
	thermo	janaf;
	equationOfState	perfectGas;
	specie	specie;
	energy	sensibleEnthalpy;
}

and a bunch of different types of thermoTypes, and could not obtain a result without an error.
For this one:

Code:

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#1  Foam::sigSegv::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#2  ? in "/lib64/libc.so.6"
#3  std::string::swap(std::string&) in "/opt/OpenFOAM/ThirdParty-v1612+/platforms/linux64/gcc-4.8.5/lib64/libstdc++.so.6"
#4  Foam::word::operator=(std::string const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so"
#5  Foam::basicThermo::splitThermoName(Foam::word const&, int) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so"
#6  Foam::autoPtr<Foam::rhoChemistryModel> Foam::basicChemistryModel::New<Foam::rhoChemistryModel>(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libchemistryModel.so"
#7  Foam::rhoChemistryModel::New(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libchemistryModel.so"
#8  Foam::combustionModels::rhoChemistryCombustion::rhoChemistryCombustion(Foam::word const&, Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libcombustionModels.so"
#9  Foam::combustionModels::laminar<Foam::combustionModels::rhoChemistryCombustion>::laminar(Foam::word const&, Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libcombustionModels.so"
#10  Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion>::PaSR(Foam::word const&, Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libcombustionModels.so"
#11  Foam::combustionModels::rhoCombustionModel::adddictionaryConstructorToTable<Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion> >::New(Foam::word const&, Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libcombustionModels.so"
#12  Foam::combustionModels::rhoCombustionModel::New(Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/lib/libcombustionModels.so"
#13  ? in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam"
#14  __libc_start_main in "/lib64/libc.so.6"
#15  ? in "/opt/OpenFOAM/OpenFOAM-v1612+/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam"
Segmentation fault (core dumped)

I can't obtain any other thermoType to solve the problem besides the original.... Should I also start from a higher initial temperature and make T_react above 0?

HRobertHS · March 11, 2017, 21:23

Please help.... I tried creating a new mesh and forming a patch from a partition in SALOME. The detonation patch doesn't register and OpenFOAM creates a default patch for it along with many other cells. I tried changing some parameters around, but still get the same thing... I then tried to use simFlow and obtain this:

March 11, 2017, 21:23		#3
HRobertHS New Member Robert Huang Join Date: Jan 2017 Posts: 11 Rep Power: 9	Please help.... I tried creating a new mesh and forming a patch from a partition in SALOME. The detonation patch doesn't register and OpenFOAM creates a default patch for it along with many other cells. I tried changing some parameters around, but still get the same thing... I then tried to use simFlow and obtain this: