From d3c6f40cece274770d0ff4ef972a9cccaea0f249 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?R=C3=A9mi=20Chassagne?= Date: Fri, 20 Oct 2023 10:29:08 +0200 Subject: [PATCH] Replace nuFraExp to nuFra in UEqns.H --- solver/pU/UEqns.H | 9 ++------- solver/pUf/UEqns.H | 10 ++-------- 2 files changed, 4 insertions(+), 15 deletions(-) diff --git a/solver/pU/UEqns.H b/solver/pU/UEqns.H index 8779559c..b79a95dc 100644 --- a/solver/pU/UEqns.H +++ b/solver/pU/UEqns.H @@ -49,19 +49,14 @@ if (SUSaniso) susb = SUS*(K/rhob)*turbulenceb->nut() *((SUS_I*iso-SUS_A*aniso) & fvc::grad(alpha)); } -// limiter for the explicit stress terms to avoid numerical instability and -// too strong limitation on the time step -volScalarField nuFraExp(nuFra); -nuFraExp.min(nuMaxExp); -nuFraExp.correctBoundaryConditions(); // Compute explicit stress term for phase a volTensorField Rca ( "Rca", - ((2.0/3.0)*I)*((nuEffa+nuFraExp/(alpha+alphaSmall))*tr(gradUaT) + ((2.0/3.0)*I)*((nuEffa+nuFra/(alpha+alphaSmall))*tr(gradUaT) + spherSigmaSGSa) - - (nuEffa+nuFraExp/(alpha+alphaSmall))*gradUaT + - (nuEffa+nuFra/(alpha+alphaSmall))*gradUaT ); // Compute concentration gradient stress contribution on cell faces for phase a diff --git a/solver/pUf/UEqns.H b/solver/pUf/UEqns.H index 299ccf57..0910a29a 100644 --- a/solver/pUf/UEqns.H +++ b/solver/pUf/UEqns.H @@ -50,19 +50,13 @@ if (SUSaniso) *((SUS_I*iso-SUS_A*aniso) & fvc::grad(alpha)); } -// limiter for the explicit stress terms to avoid numerical instability and -// too strong limitation on the time step -volScalarField nuFraExp(nuFra); -nuFraExp.min(nuMaxExp); -nuFraExp.correctBoundaryConditions(); - // Compute explicit stress term for phase a volTensorField Rca ( "Rca", - ((2.0/3.0)*I)*((nuEffa+nuFraExp/(alpha+alphaSmall))*tr(gradUaT) + ((2.0/3.0)*I)*((nuEffa+nuFra/(alpha+alphaSmall))*tr(gradUaT) + spherSigmaSGSa) - - (nuEffa+nuFraExp/(alpha+alphaSmall))*gradUaT + - (nuEffa+nuFra/(alpha+alphaSmall))*gradUaT ); // Compute concentration gradient stress contribution on cell faces for phase a