DEV: Fiel for enthalpyJump

applications/solvers/heatTransfer/buoyantPimpleFoam/Make/options

@@ -10,4 +10,6 @@ EXE_LIBS = \
     -lcompressibleTurbulenceModel \
     -lcompressibleRASModels \
     -lcompressibleLESModels \
-    -lfiniteVolume
+    -lfiniteVolume \
+    -L$(FOAM_USER_LIBBIN) \
+        -lFvPatchScalarFieldEnthalpyJump

src/finiteVolume/fields/fvPatchFields/fvPatchField/fvPatchFieldNew.C

@@ -65,8 +65,10 @@ Foam::tmp<Foam::fvPatchField<Type> > Foam::fvPatchField<Type>::New
      || actualPatchType != p.type()
     )
     {
+        Info << actualPatchType << endl;
+        Info << p.type() << endl;
         typename patchConstructorTable::iterator patchTypeCstrIter =
-            patchConstructorTablePtr_->find(p.type());
+            patchConstructorTablePtr_->find(actualPatchType);//p.type());
 
         if (patchTypeCstrIter != patchConstructorTablePtr_->end())
         {

src/thermophysicalModels/basic/Make/options

@@ -1,6 +1,8 @@
 EXE_INC = \
     -I$(LIB_SRC)/finiteVolume/lnInclude \
-    -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude
+    -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
+    -I/home/tom3/sergio/development/chtMultiRegionCoupledFoam/lnInclude
 
 LIB_LIBS = \
-    -lfiniteVolume
+    -lfiniteVolume \
+    -L$(FOAM_USER_LIBBIN)/FvPatchScalarFieldEnthalpyJump

src/thermophysicalModels/basic/basicThermo/basicThermo.C

@@ -33,6 +33,7 @@ License
 #include "fixedInternalEnergyFvPatchScalarField.H"
 #include "gradientInternalEnergyFvPatchScalarField.H"
 #include "mixedInternalEnergyFvPatchScalarField.H"
+#include "enthalpyJumpFvPatchScalarField.H"
 
 /* * * * * * * * * * * * * * * private static data * * * * * * * * * * * * * */
 
@@ -64,10 +65,14 @@ Foam::wordList Foam::basicThermo::hBoundaryTypes()
         {
             hbt[patchi] = gradientEnthalpyFvPatchScalarField::typeName;
         }
-        else if (isA<mixedFvPatchScalarField>(tbf[patchi]))
+        else if(isA<mixedFvPatchScalarField>(tbf[patchi]))
         {
             hbt[patchi] = mixedEnthalpyFvPatchScalarField::typeName;
         }
+        else if (isA<enthalpyJumpFvPatchScalarField>(tbf[patchi]))
+        {
+            hbt[patchi] = enthalpyJumpFvPatchScalarField::typeName;
+        }
     }
 
     return hbt;
@@ -85,7 +90,7 @@ void Foam::basicThermo::hBoundaryCorrection(volScalarField& h)
             refCast<gradientEnthalpyFvPatchScalarField>(hbf[patchi]).gradient()
                 = hbf[patchi].fvPatchField::snGrad();
         }
-        else if (isA<mixedEnthalpyFvPatchScalarField>(hbf[patchi]))
+        else if(isA<mixedEnthalpyFvPatchScalarField>(hbf[patchi]))
         {
             refCast<mixedEnthalpyFvPatchScalarField>(hbf[patchi]).refGrad()
                 = hbf[patchi].fvPatchField::snGrad();

src/thermophysicalModels/basic/rhoThermo/hRhoThermo/hRhoThermo.C

@@ -118,9 +118,11 @@ Foam::hRhoThermo<MixtureType>::hRhoThermo(const fvMesh& mesh)
         ),
         mesh,
         dimEnergy/dimMass,
+        this->hBoundaryTypes(),
         this->hBoundaryTypes()
     )
 {
+    Info << "hBounda" << this->hBoundaryTypes() << endl;
     scalarField& hCells = h_.internalField();
     const scalarField& TCells = this->T_.internalField();
 
@@ -136,8 +138,9 @@ Foam::hRhoThermo<MixtureType>::hRhoThermo(const fvMesh& mesh)
     }
 
     hBoundaryCorrection(h_);
-
+    Info << h_ << endl;
     calculate();
+
 }