Merge branch 'wip-support-eulersolver-heat-transfer-FO' into 'develop'

support eulersolver heat transfer fo

See merge request !333

applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionTwoPhaseEulerFoam/Make/options

@@ -36,5 +36,6 @@ EXE_LIBS = \
     -lradiationModels \
     -lregionModels \
     -lsampling \
-    -lreactingPhaseSystem \
-    -lreactingTwoPhaseSystem
+    -lreactingTwoPhaseSystem \
+    -lreactingPhaseSystem
+

src/functionObjects/field/Make/options

@@ -18,7 +18,8 @@ EXE_INC = \
     -I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
-    -I$(LIB_SRC)/fvAgglomerationMethods/pairPatchAgglomeration/lnInclude
+    -I$(LIB_SRC)/fvAgglomerationMethods/pairPatchAgglomeration/lnInclude \
+    -I$(LIB_SRC)/phaseSystemModels/reactingEulerFoam/phaseSystems/lnInclude
 
 LIB_LIBS = \
     -lfiniteVolume \
@@ -37,4 +38,5 @@ LIB_LIBS = \
     -lcompressibleTurbulenceModels \
     -lchemistryModel \
     -lreactionThermophysicalModels \
-    -lpairPatchAgglomeration
+    -lpairPatchAgglomeration \
+    -lreactingPhaseSystem

src/functionObjects/field/heatTransferCoeff/heatTransferCoeffModels/heatTransferCoeffModel/heatTransferCoeffModel.C

@@ -30,6 +30,7 @@ License
 #include "fluidThermo.H"
 #include "turbulentTransportModel.H"
 #include "turbulentFluidThermoModel.H"
+#include "phaseSystem.H"
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
@@ -90,6 +91,28 @@ Foam::heatTransferCoeffModel::q() const
             q[patchi] = alphabf[patchi]*hebf[patchi].snGrad();
         }
     }
+    else if (mesh_.foundObject<phaseSystem>("phaseProperties"))
+    {
+        const phaseSystem& fluid =
+        (
+            mesh_.lookupObject<phaseSystem>("phaseProperties")
+        );
+
+        for (label patchi : patchSet_)
+        {
+            forAll(fluid.phases(), phasei)
+            {
+                const phaseModel& phase = fluid.phases()[phasei];
+                const fvPatchScalarField& alpha =
+                    phase.boundaryField()[patchi];
+                const volScalarField& he = phase.thermo().he();
+                const volScalarField::Boundary& hebf = he.boundaryField();
+
+                q[patchi] +=
+                    alpha*phase.alphaEff(patchi)()*hebf[patchi].snGrad();
+            }
+        }
+    }
     else
     {
         FatalErrorInFunction

src/phaseSystemModels/Allwmake

@@ -4,9 +4,7 @@ cd "${0%/*}" || exit                                # Run from this directory
 
 #------------------------------------------------------------------------------
 
-wmakeLnInclude reactingEulerFoam/interfacialCompositionModels
-wmakeLnInclude reactingEulerFoam/interfacialModels
-wmakeLnInclude reactingEulerFoam/derivedFvPatchFields
+wmakeLnInclude reactingEulerFoam/phaseSystems
 
 case "$WM_COMPILER" in
 Mingw*)
@@ -18,7 +16,6 @@ esac
 
 #------------------------------------------------------------------------------
 
-wmake $targetType  reactingEulerFoam/phaseSystems
 wmake $targetType  reactingEulerFoam/interfacialModels
 wmake $targetType  reactingEulerFoam/interfacialCompositionModels
 wmake $targetType  reactingEulerFoam/derivedFvPatchFields
@@ -29,4 +26,6 @@ wmake $targetType  reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseCompr
 wmake $targetType  reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseSystem
 wmake $targetType  reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseCompressibleTurbulenceModels
 
+wmake $targetType  reactingEulerFoam/phaseSystems
+
 #------------------------------------------------------------------------------

src/phaseSystemModels/reactingEulerFoam/derivedFvPatchFields/Make/options

@@ -15,7 +15,6 @@ LIB_LIBS = \
     -lfiniteVolume \
     -lfvOptions \
     -lmeshTools \
-    -lreactingPhaseSystem \
     -lreactingEulerianInterfacialModels \
     -lreactingEulerianInterfacialCompositionModels \
     -lturbulenceModels \

src/phaseSystemModels/reactingEulerFoam/interfacialModels/Make/options

@@ -15,6 +15,5 @@ LIB_LIBS = \
     -lfluidThermophysicalModels \
     -lcompressibleTransportModels \
     -lspecie \
-    -lreactingPhaseSystem \
     -lturbulenceModels \
     -lcompressibleTurbulenceModels

src/phaseSystemModels/reactingEulerFoam/phaseSystems/Make/options

@@ -24,4 +24,9 @@ LIB_LIBS = \
     -lcompressibleTransportModels \
     -lturbulenceModels \
     -lcompressibleTurbulenceModels \
-    -lcombustionModels
+    -lcombustionModels \
+    -lreactingEulerianInterfacialCompositionModels \
+    -lreactingEulerianInterfacialModels \
+    -lreactingEulerianFvPatchFields \
+    -lmultiphaseReactingTurbulenceModels \
+    -ltwoPhaseReactingTurbulenceModels

src/phaseSystemModels/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseCompressibleTurbulenceModels/Make/options

@@ -19,7 +19,6 @@ LIB_LIBS = \
     -lfiniteVolume \
     -lfvOptions \
     -lmeshTools \
-    -lreactingPhaseSystem \
     -lreactingMultiphaseSystem \
     -lreactingEulerianInterfacialModels \
     -lreactingEulerianInterfacialCompositionModels \

src/phaseSystemModels/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/Make/options

@@ -22,5 +22,4 @@ LIB_LIBS = \
     -lfvOptions \
     -lmeshTools \
     -lsampling \
-    -lcombustionModels \
-    -lreactingPhaseSystem
+    -lcombustionModels

src/phaseSystemModels/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseCompressibleTurbulenceModels/Make/options

@@ -17,7 +17,6 @@ LIB_LIBS = \
     -lfiniteVolume \
     -lfvOptions \
     -lmeshTools \
-    -lreactingPhaseSystem \
     -lreactingTwoPhaseSystem \
     -lreactingEulerianInterfacialModels \
     -lreactingEulerianInterfacialCompositionModels \

src/phaseSystemModels/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseSystem/Make/options

@@ -21,6 +21,5 @@ LIB_LIBS = \
     -lmeshTools \
     -lsampling \
     -lcombustionModels \
-    -lreactingPhaseSystem \
     -lturbulenceModels \
     -lcompressibleTurbulenceModels

tutorials/heatTransfer/chtMultiRegionTwoPhaseEulerFoam/solidQuenching2D/system/controlDict

@@ -61,6 +61,41 @@ adjustTimeStep  yes;
 
 functions
 {
+    htc
+    {
+        type            heatTransferCoeff;
+        libs            (fieldFunctionObjects);
+        region          water;
+        field           T.liquid;
+        writeControl    outputTime;
+        writeInterval   1;
+        htcModel        fixedReferenceTemperature;
+        patches         (water_to_solid);
+        TRef            373;
+    }
+
+    htcSurfaceFieldValue
+    {
+        type            surfaceFieldValue;
+        libs            (fieldFunctionObjects);
+        log             on;
+        enabled         true;
+
+        writeControl    outputTime;
+        writeInterval   1;
+
+        writeFields     false;
+        surfaceFormat   none;
+        writeArea       false;
+
+        region          water;
+        regionType      patch;
+        name            water_to_solid;
+
+        operation       areaAverage;
+
+        fields          ( heatTransferCoeff(T.liquid) );
+    }
 }
 
 // ************************************************************************* //