stardis

commit cfdd7b6365eb7d59f44cb9860519ee175f2772bd
parent cc907b099a120896d37d364cf4b2704b9cdfc1ac
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Tue, 19 Oct 2021 16:48:26 +0200

Send Tref to solver

Diffstat:
M
src/stardis-intface.c
 | 
17
+++++++++++++++--
M
src/stardis-intface.h
 | 
1
+

2 files changed, 16 insertions(+), 2 deletions(-)
diff --git a/src/stardis-intface.c b/src/stardis-intface.c
@@ -56,6 +56,16 @@ interface_get_flux
 }
 
 static double
+interface_get_ref_temperature
+  (const struct sdis_interface_fragment* frag,
+   struct sdis_data* data)
+{
+  const struct intface* interface_props = sdis_data_cget(data);
+  (void)frag;
+  return interface_props->ref_temperature;
+}
+
+static double
 interface_get_emissivity
   (const struct sdis_interface_fragment* frag,
    struct sdis_data* data)
@@ -249,8 +259,9 @@ create_intface
         if(connect->d.h_boundary.hc > 0) {
           interface_shader.convection_coef = interface_get_convection_coef;
         }
+        ASSERT(fluid_side_shader);
+        fluid_side_shader->reference_temperature = interface_get_ref_temperature;
         if(connect->d.h_boundary.emissivity > 0) {
-          ASSERT(fluid_side_shader);
           fluid_side_shader->emissivity = interface_get_emissivity;
           fluid_side_shader->specular_fraction = interface_get_alpha;
         }
@@ -278,6 +289,7 @@ create_intface
         /* Set emissivity to 1 to allow radiative paths comming from
          * a possible external fluid to 'see' the imposed T */
         ASSERT(fluid_side_shader);
+        fluid_side_shader->reference_temperature = interface_get_ref_temperature;
         fluid_side_shader->emissivity = interface_get_emissivity;
         interface_props->emissivity = 1;
         ASSERT(connect->d.t_boundary.imposed_temperature >= 0);
@@ -317,8 +329,9 @@ create_intface
         if(connect->d.sf_connect.hc > 0) {
           interface_shader.convection_coef = interface_get_convection_coef;
         }
+        ASSERT(fluid_side_shader);
+        fluid_side_shader->reference_temperature = interface_get_ref_temperature;
         if(connect->d.sf_connect.emissivity > 0) {
-          ASSERT(fluid_side_shader);
           fluid_side_shader->emissivity = interface_get_emissivity;
           fluid_side_shader->specular_fraction = interface_get_alpha;
         }
diff --git a/src/stardis-intface.h b/src/stardis-intface.h
@@ -29,6 +29,7 @@ struct dummies;
 struct intface {
   /* fluid - solid */
   double hc;
+  double ref_temperature;
   double emissivity;
   double alpha;
   /* solid - solid */