Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Running hy2Foam with 5 neutral species and 19 reactions #30

Closed
iainsymthe opened this issue Jan 13, 2020 · 2 comments
Closed

Running hy2Foam with 5 neutral species and 19 reactions #30

iainsymthe opened this issue Jan 13, 2020 · 2 comments

Comments

@iainsymthe
Copy link

I am trying to run hy2Foam with all file as as in the example case cylinderReactingMach20 but with the shocktube example. This case works fine with two species N2 and N and the associated reactions in hTCReactionsQK/thermoDEM (default setting for cylinderReactingMach20). When I try to run the same case with all the neutral species on all 19 reactions, the code crashes immediately (see part of output file). I have the 0/ directory set-up to include O2, O, NO species and include/initialConditions has Y_N2 = 0.77, Y_O2=0.23. Are there any other changes to be made?

===== part of output file
Selecting thermodynamics package
{
type heRho2Thermo;
mixture reacting2Mixture;
transport BlottnerEucken;
thermo decoupledEnergyModes;
energy sensible2InternalEnergy;
equationOfState perfect2Gas;
specie advancedSpecie;
}

Selecting chemistry2Reader foam2ChemistryReader
chemistry2Model: Number of species = 5 and reactions = 19
using integrated reaction rate
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib64/libc.so.6"
#3 ? at rho2ReactionThermos.C:?
#4 Foam::species::multiThermo<Foam::decoupledEnergyModesThermo<Foam::perfect2GasFoam::advancedSpecie >, Foam::sensible2InternalEnergy>::evels(double, double) const at ??:?
#5 Foam::Specie2Mixture<Foam::reacting2Mixture<Foam::BlottnerEuckenTransport<Foam::species::multiThermo<Foam::decoupledEnergyModesThermo<Foam::perfect2GasFoam::advancedSpecie >, Foam::sensible2InternalEnergy> > > >::HEvel(int, double, double) const at ??:?
#6 Foam::rho2ReactionThermo::initialiseLight() at ??:?
#7 Foam::rho2ReactionThermo::initialise2T_Light() at ??:?
#8 ? at ??:?
#9 ? at ??:?
#10 __libc_start_main in "/lib64/libc.so.6"
#11 ? at ??:?
Floating point exception (core dumped)

@iainsymthe
Copy link
Author

I was able to run the case. Must have been some typo in some file.

@iainsymthe
Copy link
Author

iainsymthe commented Dec 20, 2020 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant