CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

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

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

Sponsored Links

Reply
 
LinkBack Thread Tools Display Modes
Old   March 11, 2017, 05:05
Default Jet Engine case, initiating a detonation wave, and reactingParcelFoam errors
  #1
New Member
 
Robert Huang
Join Date: Jan 2017
Posts: 11
Rep Power: 2
HRobertHS is on a distinguished road
Sponsored Links
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

Last edited by HRobertHS; March 11, 2017 at 15:55. Reason: added fvSchemes and fvSolutions for easier download
HRobertHS is offline   Reply With Quote
Sponsored Links

Old   March 11, 2017, 15:49
Default Tried so far:
  #2
New Member
 
Robert Huang
Join Date: Jan 2017
Posts: 11
Rep Power: 2
HRobertHS is on a distinguished road
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 is offline   Reply With Quote

Old   March 11, 2017, 21:23
Default
  #3
New Member
 
Robert Huang
Join Date: Jan 2017
Posts: 11
Rep Power: 2
HRobertHS is on a distinguished road
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:

HRobertHS is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Sponsored Links


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