stardis

commit 8e4624192dc1a9063ddfc58591af99f861fdb874
parent c19b7f6c6f29464d3893ce55af408c383e21eeb3
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Mon, 13 Apr 2020 14:09:39 +0200

Fix imposed T boundaries

Diffstat:
M
src/stardis-intface.c
 | 
18
++++++++----------

1 file changed, 8 insertions(+), 10 deletions(-)
diff --git a/src/stardis-intface.c b/src/stardis-intface.c
@@ -254,17 +254,15 @@ create_intface
           == (connect->type == DESC_BOUND_T_FOR_FLUID ? SDIS_SOLID: SDIS_FLUID));
         connection_count++;
         boundary_count++;
-        if(front_defined) {
+        if(front_defined)
           back_med = media[ext_id];
-          /* We set the known T inside
-           * TODO: should be outside to allow contact resistances */
-          interface_shader.front.temperature = interface_get_temperature;
-        } else {
-          front_med = media[ext_id];
-          /* We set the known T inside
-           * TODO: should be outside to allow contact resistances */
-          interface_shader.back.temperature = interface_get_temperature;
-        }
+        else front_med = media[ext_id];
+        /* The imposed T is for the 2 sides (until there is contact resistnces) */
+        interface_shader.front.temperature = interface_get_temperature;
+        interface_shader.back.temperature = interface_get_temperature;
+        /* Set emissivity to 1 to allow radiative paths to 'see' the imposed T */
+        fluid_side_shader->emissivity = interface_get_emissivity;
+        interface_props->emissivity = 1;
         ASSERT(connect->d.t_boundary.imposed_temperature >= 0);
         interface_props->imposed_temperature
           = connect->d.t_boundary.imposed_temperature;