stardis

Perform coupled heat transfer calculations
git clone git://git.meso-star.fr/stardis.git
Log | Files | Refs | README | LICENSE

commit 2bfda89bbafb986c71ca229b6d6ce9342ec7098e
parent be606b926fdd3890934cb16777d84603cfb645c3
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Thu, 22 Sep 2022 16:35:00 +0200

Fix: set flux on both sides for H+F boundaries

Diffstat:
Msrc/stardis-intface.c | 12++++++------
1 file changed, 6 insertions(+), 6 deletions(-)

diff --git a/src/stardis-intface.c b/src/stardis-intface.c @@ -434,7 +434,8 @@ create_intface ASSERT(intface->d.hf_boundary->imposed_flux != SDIS_FLUX_NONE); if(intface->d.hf_boundary->imposed_flux) { interface_props->imposed_flux = intface->d.hf_boundary->imposed_flux; - fluid_side_shader->flux = interface_get_flux; + interface_shader.front.flux = interface_get_flux; + interface_shader.back.flux = interface_get_flux; } if(intface->d.hf_boundary->emissivity > 0) { fluid_side_shader->emissivity = interface_get_emissivity; @@ -504,13 +505,12 @@ create_intface interface_shader.convection_coef_upper_bound = hf_bound_prog_get_hmax(intface->d.hf_boundary_prog); ASSERT(interface_shader.convection_coef_upper_bound >= 0); - interface_shader.convection_coef - = intface_prog_get_hc; - fluid_side_shader->reference_temperature - = intface_prog_get_ref_temp; + interface_shader.convection_coef = intface_prog_get_hc; + fluid_side_shader->reference_temperature = intface_prog_get_ref_temp; fluid_side_shader->emissivity = intface_prog_get_emissivity; fluid_side_shader->specular_fraction = intface_prog_get_alpha; - fluid_side_shader->flux = intface_prog_get_flux; + interface_shader.front.flux = intface_prog_get_flux; + interface_shader.back.flux = intface_prog_get_flux; interface_props->get_hc = intface->d.hf_boundary_prog->hc; interface_props->get_emissivity = intface->d.hf_boundary_prog->emissivity; interface_props->get_alpha = intface->d.hf_boundary_prog->alpha;