From f388fa7faee59477d7a5fc046f885816b30b6ad0 Mon Sep 17 00:00:00 2001 From: "code.hesed" Date: Fri, 24 Jun 2022 22:59:54 +0900 Subject: [PATCH] kOmegaSSTLM: edited delta to avoid floating point exception 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. --- .../momentumTransportModels/RAS/kOmegaSSTLM/kOmegaSSTLM.C | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/MomentumTransportModels/momentumTransportModels/RAS/kOmegaSSTLM/kOmegaSSTLM.C b/src/MomentumTransportModels/momentumTransportModels/RAS/kOmegaSSTLM/kOmegaSSTLM.C index ff9bb0d70f..64e2bd0fd1 100644 --- a/src/MomentumTransportModels/momentumTransportModels/RAS/kOmegaSSTLM/kOmegaSSTLM.C +++ b/src/MomentumTransportModels/momentumTransportModels/RAS/kOmegaSSTLM/kOmegaSSTLM.C @@ -82,7 +82,7 @@ tmp kOmegaSSTLM::Fthetat const volScalarField::Internal& omega = this->omega_(); const volScalarField::Internal& y = this->y_(); - const volScalarField::Internal delta(375*Omega*nu*ReThetat_()*y/sqr(Us)); + const volScalarField::Internal delta(375*max(Omega, dimensionedScalar("omegaSmall", dimless/dimTime, small))*nu*ReThetat_()*y/sqr(Us)); const volScalarField::Internal ReOmega(sqr(y)*omega/nu); const volScalarField::Internal Fwake(exp(-sqr(ReOmega/1e5)));