Merge pull request #1226 from su2code/feature_heatcoeff_BC

Heat Transfer boundary condition

Author: TobiKattmann

.github/workflows/regression.yml

@@ -87,6 +87,7 @@ jobs:
       - name: Run Tests in Container
         uses: docker://su2code/test-su2:20200303
         with:
+          # -t <Tutorials-branch> -c <Testcases-branch>
           args: -b ${{github.ref}} -t develop -c develop -s ${{matrix.testscript}}
   unit_tests:
     runs-on: ubuntu-latest

Common/include/CConfig.hpp

@@ -1086,7 +1086,8 @@ enum BC_TYPE {
   ENGINE_EXHAUST = 22,        /*!< \brief Boundary nacelle exhaust. */
   RIEMANN_BOUNDARY= 24,       /*!< \brief Riemann Boundary definition. */
   ISOTHERMAL = 25,            /*!< \brief No slip isothermal wall boundary condition. */
-  HEAT_FLUX  = 26,            /*!< \brief No slip constant heat flux wall boundary condition. */
+  HEAT_FLUX = 26,             /*!< \brief No slip constant heat flux wall boundary condition. */
+  HEAT_TRANSFER = 27,         /*!< \brief No slip heat transfer boundary condition. */
   ACTDISK_INLET = 32,         /*!< \brief Actuator disk inlet boundary definition. */
   ACTDISK_OUTLET = 33,        /*!< \brief Actuator disk outlet boundary definition. */
   CLAMPED_BOUNDARY = 34,      /*!< \brief Clamped Boundary definition. */

Common/include/option_structure.hpp

@@ -38,7 +38,7 @@
 class CIncNSSolver final : public CIncEulerSolver {
 
   /*!
-   * \brief Generic implementation of the isothermal and heatflux walls.
+   * \brief Generic implementation of the isothermal, heatflux and heat-transfer/convection walls.
    */
   void BC_Wall_Generic(const CGeometry *geometry,
                        const CConfig *config,
@@ -157,6 +157,16 @@ class CIncNSSolver final : public CIncEulerSolver {
                           CConfig *config,
                           unsigned short val_marker) override;
 
+  /*!
+   * \brief Impose a heat flux by prescribing a heat transfer coefficient and a temperature at infinity.
+   * \param[in] geometry - Geometrical definition of the problem.
+   * \param[in] config - Definition of the particular problem.
+   * \param[in] val_marker - Surface marker where the boundary condition is applied.
+   */
+  void BC_HeatTransfer_Wall(const CGeometry *geometry,
+                            const CConfig *config,
+                            const unsigned short val_marker) override;
+
   /*!
    * \brief Impose the (received) conjugate heat variables.
    * \param[in] geometry - Geometrical definition of the problem.

Common/src/CConfig.cpp

@@ -65,7 +65,7 @@ class CNSSolver final : public CEulerSolver {
    */
   void AddDynamicGridResidualContribution(unsigned long iPoint,
                                           unsigned long Point_Normal,
-                                          CGeometry* geometry,
+                                          const CGeometry* geometry,
                                           const su2double* UnitNormal,
                                           su2double Area,
                                           const su2double* GridVel,
@@ -96,6 +96,14 @@ class CNSSolver final : public CEulerSolver {
                                   unsigned short val_marker,
                                   bool cht_mode = false);
 
+  /*!
+   * \brief Generic implementation of the heatflux and heat-transfer/convection walls.
+   */
+  void BC_HeatFlux_Wall_Generic(const CGeometry *geometry,
+                                const CConfig *config,
+                                unsigned short val_marker,
+                                unsigned short kind_boundary);
+
   /*!
    * \brief Compute the viscous contribution for a particular edge.
    * \param[in] iEdge - Edge for which the flux and Jacobians are to be computed.
@@ -185,6 +193,16 @@ class CNSSolver final : public CEulerSolver {
                         CConfig *config,
                         unsigned short val_marker) override;
 
+  /*!
+   * \brief Impose a heat flux by prescribing a heat transfer coefficient and a temperature at infinity.
+   * \param[in] geometry - Geometrical definition of the problem.
+   * \param[in] config - Definition of the particular problem.
+   * \param[in] val_marker - Surface marker where the boundary condition is applied.
+   */
+  void BC_HeatTransfer_Wall(const CGeometry *geometry,
+                            const CConfig *config,
+                            const unsigned short val_marker) override;
+
   /*!
    * \brief Impose the Navier-Stokes boundary condition (strong).
    * \param[in] geometry - Geometrical definition of the problem.

SU2_CFD/include/solvers/CIncNSSolver.hpp

@@ -1061,6 +1061,16 @@ class CSolver {
                                        CConfig *config,
                                        unsigned short val_marker) { }
 
+  /*!
+   * \brief Impose a heat flux by prescribing a heat transfer coefficient and a temperature at infinity.
+   * \param[in] geometry - Geometrical definition of the problem.
+   * \param[in] config - Definition of the particular problem.
+   * \param[in] val_marker - Surface marker where the boundary condition is applied.
+   */
+  inline virtual void BC_HeatTransfer_Wall(const CGeometry *geometry,
+                                           const CConfig *config,
+                                           const unsigned short val_marker) { }
+
   /*!
    * \brief A virtual member.
    * \param[in] geometry - Geometrical definition of the problem.

SU2_CFD/include/solvers/CNSSolver.hpp

@@ -59,9 +59,9 @@ void CFluidModel::SetThermalConductivityModel(const CConfig* config) {
     case CONDUCTIVITYMODEL::CONSTANT:
       if (config->GetKind_ConductivityModel_Turb() == CONDUCTIVITYMODEL_TURB::CONSTANT_PRANDTL) {
         ThermalConductivity = unique_ptr<CConstantConductivityRANS>(
-            new CConstantConductivityRANS(config->GetKt_ConstantND(), config->GetPrandtl_Turb()));
+            new CConstantConductivityRANS(config->GetThermal_Conductivity_ConstantND(), config->GetPrandtl_Turb()));
       } else {
-        ThermalConductivity = unique_ptr<CConstantConductivity>(new CConstantConductivity(config->GetKt_ConstantND()));
+        ThermalConductivity = unique_ptr<CConstantConductivity>(new CConstantConductivity(config->GetThermal_Conductivity_ConstantND()));
       }
       break;
     case CONDUCTIVITYMODEL::CONSTANT_PRANDTL:

SU2_CFD/include/solvers/CSolver.hpp

@@ -153,6 +153,9 @@ void CIntegration::Space_Integration(CGeometry *geometry,
       case HEAT_FLUX:
         solver_container[MainSolver]->BC_HeatFlux_Wall(geometry, solver_container, conv_bound_numerics, visc_bound_numerics, config, iMarker);
         break;
+      case HEAT_TRANSFER:
+        solver_container[MainSolver]->BC_HeatTransfer_Wall(geometry, config, iMarker);
+        break;
       case CUSTOM_BOUNDARY:
         solver_container[MainSolver]->BC_Custom(geometry, solver_container, conv_bound_numerics, visc_bound_numerics, config, iMarker);
         break;

SU2_CFD/src/fluid/CFluidModel.cpp

@@ -124,13 +124,14 @@ void CConjugateHeatInterface::GetDonor_Variable(CSolver *donor_solution, CGeomet
 
     /*--- Heat solver stand-alone case ---*/
 
-    const su2double thermal_diffusivity = donor_config->GetThermalDiffusivity_Solid();
+    const su2double thermal_diffusivity = donor_config->GetThermalDiffusivity();
     heat_flux_density = thermal_diffusivity*dTdn;
 
+
     if ((donor_config->GetKind_CHT_Coupling() == CHT_COUPLING::DIRECT_TEMPERATURE_ROBIN_HEATFLUX) ||
         (donor_config->GetKind_CHT_Coupling() == CHT_COUPLING::AVERAGED_TEMPERATURE_ROBIN_HEATFLUX)) {
 
-      const su2double rho_cp_solid = donor_config->GetSpecific_Heat_Cp()*donor_config->GetDensity_Solid();
+      const su2double rho_cp_solid = donor_config->GetSpecific_Heat_Cp()*donor_config->GetMaterialDensity(0);
       conductivity_over_dist  = thermal_diffusivity*rho_cp_solid/dist;
     }
   }