From 8d2e5ecae86bd935148d4a1ef5aefb62bb47c0c2 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Wed, 31 May 2017 10:37:11 +0100 Subject: [PATCH] CourantNo: Removed check for nInternalFaces which is not correct in general in parallel --- src/dynamicFvMesh/include/meshCourantNo.H | 3 +-- .../cfdTools/compressible/compressibleCourantNo.H | 6 ++---- src/finiteVolume/cfdTools/incompressible/CourantNo.H | 3 +-- 3 files changed, 4 insertions(+), 8 deletions(-) diff --git a/src/dynamicFvMesh/include/meshCourantNo.H b/src/dynamicFvMesh/include/meshCourantNo.H index 7e014529a..13a8c48d0 100644 --- a/src/dynamicFvMesh/include/meshCourantNo.H +++ b/src/dynamicFvMesh/include/meshCourantNo.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -32,7 +32,6 @@ Description scalar meshCoNum = 0.0; scalar meanMeshCoNum = 0.0; -if (mesh.nInternalFaces()) { scalarField sumPhi ( diff --git a/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H b/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H index 765115e9b..bdfe7063c 100644 --- a/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H +++ b/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -32,12 +32,10 @@ Description scalar CoNum = 0.0; scalar meanCoNum = 0.0; -if (mesh.nInternalFaces()) { scalarField sumPhi ( - fvc::surfaceSum(mag(phi))().primitiveField() - / rho.primitiveField() + fvc::surfaceSum(mag(phi))().primitiveField()/rho.primitiveField() ); CoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue(); diff --git a/src/finiteVolume/cfdTools/incompressible/CourantNo.H b/src/finiteVolume/cfdTools/incompressible/CourantNo.H index 23c48e2ae..9bd9b348f 100644 --- a/src/finiteVolume/cfdTools/incompressible/CourantNo.H +++ b/src/finiteVolume/cfdTools/incompressible/CourantNo.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -32,7 +32,6 @@ Description scalar CoNum = 0.0; scalar meanCoNum = 0.0; -if (mesh.nInternalFaces()) { scalarField sumPhi (