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

kOmegaSSTLM: edited delta to avoid floating point exception #41

Open
wants to merge 1 commit into
base: master
Choose a base branch
from

Conversation

CodeHesed
Copy link

@CodeHesed CodeHesed commented Jun 24, 2022

In some cases (e.g. airfoil with non zero AoA), floating point exception occurs due to zero values in delta of kOmegaSSTLM.
added max function to the Omega part of delta such that delta will change zero values into a small value.

In cases of Uy != 0, floating point exception occurs due to zero values in delta of kOmegaSSTLM.
added max function to the Omega part of delta such that delta will change zero values into a small value.
@CodeHesed
Copy link
Author

Example of error message:

#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::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator/<Foam::volMesh>(Foam::DimensionedField<double, Foam::volMesh> const&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:?
#5  Foam::RASModels::kOmegaSSTLM<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::Fthetat(Foam::DimensionedField<double, Foam::volMesh> const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::DimensionedField<double, Foam::volMesh> const&) const at ??:?
	#6  Foam::RASModels::kOmegaSSTLM<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correctReThetatGammaInt() at ??:?
#7  Foam::RASModels::kOmegaSSTLM<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correct() at ??:?
#8  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
Floating point exception

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

Successfully merging this pull request may close these issues.

1 participant