From 76dc6453054872cda0b4c5c581d4213e413c71b9 Mon Sep 17 00:00:00 2001 From: Christophe Geuzaine <cgeuzaine@uliege.be> Date: Wed, 21 Apr 2021 19:51:25 +0200 Subject: [PATCH] uncomment conv+rad bc --- Magnetometer/Jacobian_Lib.pro | 4 ++++ Magnetometer/Thermal.pro | 24 +++++++++++++++--------- 2 files changed, 19 insertions(+), 9 deletions(-) diff --git a/Magnetometer/Jacobian_Lib.pro b/Magnetometer/Jacobian_Lib.pro index 4cd3cb1..3880822 100644 --- a/Magnetometer/Jacobian_Lib.pro +++ b/Magnetometer/Jacobian_Lib.pro @@ -42,6 +42,10 @@ Jacobian { { Region All ; Jacobian Vol ; } } } + { Name Sur ; + Case { { Region All ; Jacobian Sur ; } + } + } } diff --git a/Magnetometer/Thermal.pro b/Magnetometer/Thermal.pro index d22aa14..4bfab1a 100644 --- a/Magnetometer/Thermal.pro +++ b/Magnetometer/Thermal.pro @@ -10,6 +10,14 @@ FunctionSpace { } } +Group { + DefineGroup [ SurfaceConv_The, SurfaceRad_The ]; +} + +Function { + DefineFunction [ h, T0, rad ]; +} + Formulation { { Name Thermal_T ; Type FemEquation ; Quantity { @@ -26,18 +34,16 @@ Formulation { Galerkin { [ - /* 0.5 * */ 1/sigma[] * SquNorm[sigma[]*{d v}] , {T} ] ; In Domain_The ; Integration GradGrad ; Jacobian Vol ; } - /* + Galerkin { [ h[] * Dof{T} , {T} ] ; // Convection boundary condition In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; } - Galerkin { [ -h[] * TemperatureConv[] , {T} ] ; - In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; } - */ - /* - Galerkin { [ rad[] * {T}^4 , {T} ] ; // Convection boundary condition + Galerkin { [ -h[] * T0[] , {T} ] ; In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; } - Galerkin { [ -rad[] * TemperatureConv[]^4 , {T} ] ; - In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; } - */ + + Galerkin { [ rad[] * {T}^4 , {T} ] ; // Radiation boundary condition + In SurfaceRad_The ; Integration GradGrad ; Jacobian Sur ; } + Galerkin { [ -rad[] * T0[]^4 , {T} ] ; + In SurfaceRad_The ; Integration GradGrad ; Jacobian Sur ; } } } } -- GitLab