Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Feature solid solid cht #52

Merged
merged 3 commits into from
Aug 21, 2024
Merged

Feature solid solid cht #52

Show file tree
Hide file tree
Changes from 2 commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
0f817a8
Added gmsh files for generating tutorial mesh
EvertBunschoten Aug 6, 2024
acc5cf1
Added configuration files for ss-CHT tutorial
EvertBunschoten Aug 6, 2024
779d6b5
Added mesh files from gmsh scripts
EvertBunschoten Aug 21, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
55 changes: 55 additions & 0 deletions multiphysics/contact_resistance_cht/fluid.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,55 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% SU2 configuration file %
% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer %
% with contact resistance. %
% Author: E.C.Bunschoten %
% Date: August 6, 2024 %
% File Version 8.0.1 "Harrier" %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%

SOLVER= INC_NAVIER_STOKES
INC_ENERGY_EQUATION= YES

% ------------------- INCOMPRESSIBLE FREE-STREAM DEFINITION -------------------%

FREESTREAM_TEMPERATURE= 300

INC_VELOCITY_INIT=(0.1, 0.0, 0.0)
FLUID_MODEL=INC_IDEAL_GAS

VISCOSITY_MODEL=CONSTANT_VISCOSITY
MU_CONSTANT=1e-5

REF_DIMENSIONALIZATION= DIMENSIONAL
INC_NONDIM= DIMENSIONAL

% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
%
INC_INLET_TYPE=VELOCITY_INLET
INC_OUTLET_TYPE=PRESSURE_OUTLET
MARKER_SYM= ( side_3 )
MARKER_INLET= ( inlet, 300, 0.1, 1.0, 0.0, 0.0 )
MARKER_OUTLET= ( outlet, 0.0)
MARKER_INTERNAL= (fluid)
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%

CFL_NUMBER= 100

% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
%
CONV_NUM_METHOD_FLOW= FDS
MUSCL_FLOW= YES

% --------------------------- CONVERGENCE PARAMETERS --------------------------%

CONV_RESIDUAL_MINVAL= -12

% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
%
MESH_FILENAME= fluid_mesh.su2
CONV_FILENAME= history_fluid_3
TABULAR_FORMAT= CSV
60 changes: 60 additions & 0 deletions multiphysics/contact_resistance_cht/fluid_3.geo
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,60 @@
//+
Point(1) = {0, 0, 0, 1.0};
//+
Point(2) = {-0.01, 0, 0, 1.0};
//+
Point(3) = {0.01, 0, 0, 1.0};
//+
Point(4) = {0.01, -0.005, 0, 1.0};
//+
Point(5) = {0.0, -0.005, 0, 1.0};
//+
Point(6) = {-0.01, -0.005, 0, 1.0};
//+
Line(1) = {6, 2};
//+
Line(2) = {5, 1};
//+
Line(3) = {4, 3};
//+
Line(4) = {2, 1};
//+
Line(5) = {1, 3};
//+
Line(6) = {6, 5};
//+
Line(7) = {5, 4};
//+
Curve Loop(1) = {1, 4, -2, -6};
//+
Plane Surface(1) = {1};
//+
Curve Loop(2) = {5, -3, -7, 2};
//+
Plane Surface(2) = {2};
//+
Physical Curve("inlet", 8) = {1};
//+
Physical Curve("outlet", 9) = {3};
//+
Physical Curve("side_3", 10) = {6, 7};
//+
Physical Curve("cht_interface_3_1", 11) = {4};
//+
Physical Curve("cht_interface_3_2", 12) = {5};
//+
Physical Surface("fluid", 13) = {1, 2};
//+
Transfinite Curve {1, 2, 3} = 20 Using Progression 0.8;
//+
Transfinite Curve {4, 6, 7, 5} = 80 Using Progression 1;
//+
Transfinite Surface {1} = {6, 5, 1, 2};
//+
Transfinite Surface {2} = {5, 4, 3, 1};
//+
Recombine Surface {1, 2};

Mesh 2;

Save "fluid_mesh.su2";
31 changes: 31 additions & 0 deletions multiphysics/contact_resistance_cht/master.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,31 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% SU2 configuration file %
% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer %
% with contact resistance. %
% Author: E.C.Bunschoten %
% Date: August 6, 2024 %
% File Version 8.0.1 "Harrier" %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
%
SOLVER= MULTIPHYSICS
MATH_PROBLEM= DIRECT
MULTIZONE_MESH=NO
CONFIG_LIST = (solid_1.cfg, solid_2.cfg, fluid.cfg)
MARKER_ZONE_INTERFACE= (cht_interface_1_2, cht_interface_2_1, \
cht_interface_1_3, cht_interface_3_1, \
cht_interface_2_3, cht_interface_3_2)
MARKER_CHT_INTERFACE= (cht_interface_1_2, cht_interface_2_1, \
cht_interface_1_3, cht_interface_3_1, \
cht_interface_2_3, cht_interface_3_2)
CHT_INTERFACE_CONTACT_RESISTANCE = (1e-5, 0, 0)

CHT_COUPLING_METHOD= DIRECT_TEMPERATURE_ROBIN_HEATFLUX
TIME_DOMAIN = NO
OUTER_ITER = 2000
WRT_ZONE_HIST=YES

OUTPUT_FILES=(RESTART, PARAVIEW_MULTIBLOCK)
47 changes: 47 additions & 0 deletions multiphysics/contact_resistance_cht/solid_1.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,47 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% SU2 configuration file %
% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer %
% with contact resistance. %
% Author: E.C.Bunschoten %
% Date: August 6, 2024 %
% File Version 8.0.1 "Harrier" %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
%
SOLVER= HEAT_EQUATION


% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
%
MARKER_ISOTHERMAL= ( isothermal_wall_1, 500.0 )
MARKER_SYM=(side_1)
MARKER_MONITORING= ( NONE )
MARKER_INTERNAL=(solid_1)

% ---------------- (SOLIDS) CONDUCTION CONDITION DEFINITION -------------------%

% We should keep the dimensionalization of the coupled flow solver
INC_NONDIM= DIMENSIONAL
FREESTREAM_TEMPERATURE= 300.0

% Properties of stainless steel
MATERIAL_DENSITY= 8935
SPECIFIC_HEAT_CP= 3850
THERMAL_CONDUCTIVITY_CONSTANT=26

% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%

CFL_NUMBER= 100.0


% -------------------- HEAT NUMERICAL METHOD DEFINITION -----------------------%

TIME_DISCRE_HEAT= EULER_IMPLICIT

% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
TABULAR_FORMAT= CSV
CONV_FILENAME= history_solid_1
MESH_FILENAME=solid_mesh_1.su2
42 changes: 42 additions & 0 deletions multiphysics/contact_resistance_cht/solid_1.geo
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,42 @@
//+
Point(1) = {0, 0, 0, 1.0};
//+
Point(2) = {-0.01, 0, 0, 1.0};
//+
Point(3) = {-0.01, 0.01, 0, 1.0};
//+
Point(4) = {0, 0.01, 0, 1.0};
//+
Line(1) = {2, 3};
//+
Line(2) = {3, 4};
//+
Line(3) = {4, 1};
//+
Line(4) = {1, 2};
//+
Curve Loop(1) = {1, 2, 3, 4};
//+
Plane Surface(1) = {1};
//+
Physical Curve("isothermal_wall_1", 5) = {1};
//+
Physical Curve("cht_interface_1_2", 6) = {3};
//+
Physical Curve("side_1", 7) = {2};
//+
Physical Curve("cht_interface_1_3", 8) = {4};
//+
Physical Surface("solid_1", 9) = {1};
//+
Transfinite Curve {1, 2, 3, 4} = 40 Using Progression 1;
//+
Transfinite Surface {1} = {2, 1, 4, 3};
//+
Recombine Surface {1};
//+
Recombine Surface {1};

Mesh 2;

Save "solid_mesh_1.su2";
41 changes: 41 additions & 0 deletions multiphysics/contact_resistance_cht/solid_2.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,41 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% SU2 configuration file %
% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer %
% with contact resistance. %
% Author: E.C.Bunschoten %
% Date: August 6, 2024 %
% File Version 8.0.1 "Harrier" %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
%
SOLVER= HEAT_EQUATION

% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
%
MARKER_ISOTHERMAL= ( isothermal_wall_2, 300.0 )
MARKER_SYM=(side_2)
MARKER_INTERNAL=(solid_2)

% ---------------- (SOLIDS) CONDUCTION CONDITION DEFINITION -------------------%

% We should keep the dimensionalization of the coupled flow solver
INC_NONDIM= DIMENSIONAL
FREESTREAM_TEMPERATURE= 300.0

MATERIAL_DENSITY= 8000
SPECIFIC_HEAT_CP= 4420
THERMAL_CONDUCTIVITY_CONSTANT= 61

% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%

CFL_NUMBER= 100.0

% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%

MESH_FILENAME= solid_mesh_2.su2
TABULAR_FORMAT= CSV
CONV_FILENAME= history_solid_2

40 changes: 40 additions & 0 deletions multiphysics/contact_resistance_cht/solid_2.geo
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,40 @@
//+
Point(1) = {0, 0, 0, 1.0};
//+
Point(2) = {0.01, 0, 0, 1.0};
//+
Point(3) = {0.01, 0.01, 0, 1.0};
//+
Point(4) = {0.0, 0.01, 0, 1.0};
//+
Line(1) = {1, 4};
//+
Line(2) = {4, 3};
//+
Line(3) = {3, 2};
//+
Line(4) = {2, 1};
//+
Curve Loop(1) = {1, 2, 3, 4};
//+
Surface(1) = {1};
//+
Physical Curve("cht_interface_2_1", 5) = {1};
//+
Physical Curve("cht_interface_2_3", 6) = {4};
//+
Physical Curve("isothermal_wall_2", 7) = {3};
//+
Physical Curve("side_2", 8) = {2};
//+
Physical Surface("solid_2", 9) = {1};
//+
Transfinite Curve {1, 4, 3, 2} = 40 Using Progression 1;
//+
Transfinite Surface {1} = {1, 2, 3, 4};
//+
Recombine Surface {1};

Mesh 2;

Save "solid_mesh_2.su2";