BUG: Only use ".group" if internalField is not null

Author: jheylmun

src/thermodynamicModels/derivedFvPatchFields/dynamicThermoPressure/dynamicThermoPressureFvPatchScalarField.C

@@ -97,7 +97,12 @@ Foam::dynamicThermoPressureFvPatchScalarField::dynamicThermoPressureFvPatchScala
 )
 :
     fixedValueFvPatchScalarField(tppsf, iF),
-    thermoBasePatchField(tppsf, this->patch(), iF.group()),
+    thermoBasePatchField
+    (
+        tppsf,
+        this->patch(),
+        isNull(iF) ? tppsf.phaseName_ : iF.group()
+    ),
     p0_(tppsf.p0_)
 {}
 

src/thermodynamicModels/derivedFvPatchFields/pressureWaveTransmissive/pressureWaveTransmissiveFvPatchField.C

@@ -85,7 +85,12 @@ pressureWaveTransmissiveFvPatchField
 )
 :
     advectiveFvPatchField<Type>(ptpsf, iF),
-    thermoBasePatchField(ptpsf, ptpsf.patch(), iF.group())
+    thermoBasePatchField
+    (
+        ptpsf,
+        ptpsf.patch(),
+        isNull(iF) ? ptpsf.phaseName_ : iF.group()
+    )
 {}
 
 

src/thermodynamicModels/derivedFvPatchFields/thermoBasePatchField/thermoBasePatchField.C

@@ -42,11 +42,23 @@ Foam::thermoBasePatchField::thermoBasePatchField
 )
 :
     patch_(p),
-    determinePhase_(true),
+    determinePhase_(false),
     phaseName_(defaultGroup)
 {}
 
 
+Foam::thermoBasePatchField::thermoBasePatchField
+(
+    const thermoBasePatchField& tbpf,
+    const fvPatch& p
+)
+:
+    patch_(p),
+    determinePhase_(false),
+    phaseName_(tbpf.phaseName_)
+{}
+
+
 Foam::thermoBasePatchField::thermoBasePatchField
 (
     const thermoBasePatchField& tbpf,

src/thermodynamicModels/derivedFvPatchFields/thermoBasePatchField/thermoBasePatchField.H

@@ -87,6 +87,13 @@ public:
             const word&
         );
 
+        //- Construct from copy and patch
+        thermoBasePatchField
+        (
+            const thermoBasePatchField&,
+            const fvPatch&
+        );
+
         //- Construct from patch and internal field
         thermoBasePatchField
         (

src/thermodynamicModels/derivedFvPatchFields/totalThermoTemperature/totalThermoTemperatureFvPatchScalarField.C

@@ -113,7 +113,12 @@ totalThermoTemperatureFvPatchScalarField
 )
 :
     fixedValueFvPatchScalarField(tppsf, iF),
-    thermoBasePatchField(tppsf, this->patch(), iF.group()),
+    thermoBasePatchField
+    (
+        tppsf,
+        this->patch(),
+        isNull(iF) ? tppsf.phaseName_ : iF.group()
+    ),
     UName_(tppsf.UName_),
     phiName_(tppsf.phiName_),
     T0_(tppsf.T0_)