Bugfix: clean-up of total boundary conditions

src/finiteVolume/fields/fvPatchFields/derived/totalPressure/totalPressureFvPatchScalarField.C

@@ -264,29 +264,6 @@ Foam::totalPressureFvPatchScalarField::snGrad() const
 }
 
 
-Foam::tmp<Foam::scalarField>
-Foam::totalPressureFvPatchScalarField::valueInternalCoeffs
-(
-    const tmp<scalarField>&
-) const
-{
-    return tmp<scalarField>
-    (
-        new scalarField(this->size(), 0.0)
-    );
-}
-
-
-Foam::tmp<Foam::scalarField>
-Foam::totalPressureFvPatchScalarField::valueBoundaryCoeffs
-(
-    const tmp<scalarField>&
-) const
-{
-    return *this;
-}
-
-
 Foam::tmp<Foam::scalarField>
 Foam::totalPressureFvPatchScalarField::gradientInternalCoeffs() const
 {

src/finiteVolume/fields/fvPatchFields/derived/totalPressure/totalPressureFvPatchScalarField.H

@@ -25,7 +25,8 @@ Class
     Foam::totalPressureFvPatchScalarField
 
 Description
-    Foam::totalPressureFvPatchScalarField
+    Total pressure boundary condition.  The snGrad is set to zero to stop
+    back-diffusion
 
 SourceFiles
     totalPressureFvPatchScalarField.C
@@ -207,31 +208,20 @@ public:
             //- Update the coefficients associated with the patch field
             virtual void updateCoeffs();
 
-            //- Return patch-normal gradient
+            //- Return patch-normal gradient: set to zero
             virtual tmp<scalarField> snGrad() const;
 
-            //- Return the matrix diagonal coefficients corresponding to the
-            //  evaluation of the value of this patchField with given weights
-            virtual tmp<scalarField> valueInternalCoeffs
-            (
-                const tmp<scalarField>&
-            ) const;
-
-            //- Return the matrix source coefficients corresponding to the
-            //  evaluation of the value of this patchField with given weights
-            virtual tmp<scalarField> valueBoundaryCoeffs
-            (
-                const tmp<scalarField>&
-            ) const;
-
             //- Return the matrix diagonal coefficients corresponding to the
             //  evaluation of the gradient of this patchField
+            // Removed gradient contribution
             virtual tmp<scalarField> gradientInternalCoeffs() const;
 
             //- Return the matrix source coefficients corresponding to the
             //  evaluation of the gradient of this patchField
+            // Removed gradient contribution
             virtual tmp<scalarField> gradientBoundaryCoeffs() const;
 
+
         //- Write
         virtual void write(Ostream&) const;
 };

src/finiteVolume/fields/fvPatchFields/derived/totalTemperature/totalTemperatureFvPatchScalarField.C

@@ -174,6 +174,36 @@ void Foam::totalTemperatureFvPatchScalarField::updateCoeffs()
 }
 
 
+Foam::tmp<Foam::scalarField>
+Foam::totalTemperatureFvPatchScalarField::snGrad() const
+{
+    return tmp<scalarField>
+    (
+        new scalarField(this->size(), 0.0)
+    );
+}
+
+
+Foam::tmp<Foam::scalarField>
+Foam::totalTemperatureFvPatchScalarField::gradientInternalCoeffs() const
+{
+    return tmp<scalarField>
+    (
+        new scalarField(this->size(), 0.0)
+    );
+}
+
+
+Foam::tmp<Foam::scalarField>
+Foam::totalTemperatureFvPatchScalarField::gradientBoundaryCoeffs() const
+{
+    return tmp<scalarField>
+    (
+        new scalarField(this->size(), 0.0)
+    );
+}
+
+
 void Foam::totalTemperatureFvPatchScalarField::write(Ostream& os) const
 {
     fvPatchScalarField::write(os);

src/finiteVolume/fields/fvPatchFields/derived/totalTemperature/totalTemperatureFvPatchScalarField.H

@@ -25,7 +25,8 @@ Class
     Foam::totalTemperatureFvPatchScalarField
 
 Description
-    Foam::totalTemperatureFvPatchScalarField
+    Total temperature boundary condition.  The snGrad is set to zero to stop
+    back-diffusion
 
 SourceFiles
     totalTemperatureFvPatchScalarField.C
@@ -180,6 +181,18 @@ public:
             //- Update the coefficients associated with the patch field
             virtual void updateCoeffs();
 
+            //- Return patch-normal gradient: set to zero
+            virtual tmp<scalarField> snGrad() const;
+
+            //- Return the matrix diagonal coefficients corresponding to the
+            //  evaluation of the gradient of this patchField
+            // Removed gradient contribution
+            virtual tmp<scalarField> gradientInternalCoeffs() const;
+
+            //- Return the matrix source coefficients corresponding to the
+            //  evaluation of the gradient of this patchField
+            virtual tmp<scalarField> gradientBoundaryCoeffs() const;
+
 
         //- Write
         virtual void write(Ostream&) const;

src/thermophysicalModels/basic/derivedFvPatchFields/isentropicTotalPressure/isentropicTotalPressureFvPatchScalarField.C

@@ -194,29 +194,6 @@ Foam::isentropicTotalPressureFvPatchScalarField::snGrad() const
 }
 
 
-Foam::tmp<Foam::scalarField>
-Foam::isentropicTotalPressureFvPatchScalarField::valueInternalCoeffs
-(
-    const tmp<scalarField>&
-) const
-{
-    return tmp<scalarField>
-    (
-        new scalarField(this->size(), 0.0)
-    );
-}
-
-
-Foam::tmp<Foam::scalarField>
-Foam::isentropicTotalPressureFvPatchScalarField::valueBoundaryCoeffs
-(
-    const tmp<scalarField>&
-) const
-{
-    return *this;
-}
-
-
 Foam::tmp<Foam::scalarField>
 Foam::isentropicTotalPressureFvPatchScalarField::gradientInternalCoeffs() const
 {

src/thermophysicalModels/basic/derivedFvPatchFields/isentropicTotalPressure/isentropicTotalPressureFvPatchScalarField.H

@@ -25,7 +25,8 @@ Class
     Foam::isentropicTotalPressureFvPatchScalarField
 
 Description
-    Foam::isentropicTotalPressureFvPatchScalarField
+    Isentropic total pressure boundary condition.  The snGrad is set to
+    zero to stop back-diffusion
 
 SourceFiles
     isentropicTotalPressureFvPatchScalarField.C
@@ -147,18 +148,6 @@ public:
                 return UName_;
             }
 
-            //- Return the heat capacity ratio
-            scalar gamma() const
-            {
-                return gamma_;
-            }
-
-            //- Return reference to the heat capacity ratio to allow adjustment
-            scalar& gamma()
-            {
-                return gamma_;
-            }
-
             //- Return the total pressure
             const scalarField& p0() const
             {
@@ -197,6 +186,18 @@ public:
             //- Update the coefficients associated with the patch field
             virtual void updateCoeffs();
 
+            //- Return patch-normal gradient: set to zero
+            virtual tmp<scalarField> snGrad() const;
+
+            //- Return the matrix diagonal coefficients corresponding to the
+            //  evaluation of the gradient of this patchField
+            //  Removed gradient contribution
+            virtual tmp<scalarField> gradientInternalCoeffs() const;
+
+            //- Return the matrix source coefficients corresponding to the
+            //  evaluation of the gradient of this patchField
+            virtual tmp<scalarField> gradientBoundaryCoeffs() const;
+
 
         //- Write
         virtual void write(Ostream&) const;

src/thermophysicalModels/basic/derivedFvPatchFields/isentropicTotalTemperature/isentropicTotalTemperatureFvPatchScalarField.C

@@ -174,6 +174,38 @@ void Foam::isentropicTotalTemperatureFvPatchScalarField::updateCoeffs
 }
 
 
+Foam::tmp<Foam::scalarField>
+Foam::isentropicTotalTemperatureFvPatchScalarField::snGrad() const
+{
+    return tmp<scalarField>
+    (
+        new scalarField(this->size(), 0.0)
+    );
+}
+
+
+Foam::tmp<Foam::scalarField>
+Foam::isentropicTotalTemperatureFvPatchScalarField::
+gradientInternalCoeffs() const
+{
+    return tmp<scalarField>
+    (
+        new scalarField(this->size(), 0.0)
+    );
+}
+
+
+Foam::tmp<Foam::scalarField>
+Foam::isentropicTotalTemperatureFvPatchScalarField::
+gradientBoundaryCoeffs() const
+{
+    return tmp<scalarField>
+    (
+        new scalarField(this->size(), 0.0)
+    );
+}
+
+
 void Foam::isentropicTotalTemperatureFvPatchScalarField::write
 (
     Ostream& os

src/thermophysicalModels/basic/derivedFvPatchFields/isentropicTotalTemperature/isentropicTotalTemperatureFvPatchScalarField.H

@@ -25,7 +25,8 @@ Class
     Foam::isentropicTotalTemperatureFvPatchScalarField
 
 Description
-    Foam::isentropicTotalTemperatureFvPatchScalarField
+    Isentropic total temperature boundary condition.  The snGrad is set to
+    zero to stop back-diffusion
 
 SourceFiles
     isentropicTotalTemperatureFvPatchScalarField.C
@@ -169,6 +170,18 @@ public:
             //- Update the coefficients associated with the patch field
             virtual void updateCoeffs();
 
+            //- Return patch-normal gradient: set to zero
+            virtual tmp<scalarField> snGrad() const;
+
+            //- Return the matrix diagonal coefficients corresponding to the
+            //  evaluation of the gradient of this patchField
+            //  Removed gradient contribution
+            virtual tmp<scalarField> gradientInternalCoeffs() const;
+
+            //- Return the matrix source coefficients corresponding to the
+            //  evaluation of the gradient of this patchField
+            virtual tmp<scalarField> gradientBoundaryCoeffs() const;
+
 
         //- Write
         virtual void write(Ostream&) const;