Skip to content

Commit ea756d9

Browse files
committed
First version for Unsteady CHT tutorial.
1 parent 5a9d8a5 commit ea756d9

File tree

1 file changed

+75
-0
lines changed

1 file changed

+75
-0
lines changed
Lines changed: 75 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,75 @@
1+
---
2+
title: Unsteady Conjugate Heat Transfer
3+
permalink: /tutorials/Inc_Heated_Cylinders_Unsteady/
4+
---
5+
6+
## Goals
7+
8+
This tutorial is a follow-up on the [heated cylinders with conjugate heat transfer tutorial](/tutorials/Inc_Heated_Cylinders/) where a steady CHT solution was computed for a problem involving multiple physical zones.
9+
The following capabilities of SU2 will be showcased in this tutorial:
10+
11+
- Time domain and time-marching config file options (plus related ones) for unsteady simulations
12+
- Use of time iterations, outer and inner iterations
13+
- Paraview multiblock output
14+
15+
The intent of this tutorial is to demonstrate how a steady CHT simulation can be turned into an unsteady one.
16+
17+
## Resources
18+
19+
The resources for this tutorial can be found in the [Inc_Heated_Cylinders_Unsteady](https://github.com/su2code/su2code.github.io/blob/unsteady_cht_tutorial/Inc_Heated_Cylinders_Unsteady) directory in the [tutorial repository](https://github.com/su2code/su2code.github.io/blob/unsteady_cht_tutorial/). You will need the configuration files for all physical zones ([flow_cylinder.cfg](../../Inc_Heated_Cylinders_Unsteady/flow_cylinder.cfg), [solid_cylinder1.cfg](../../Inc_Heated_Cylinders_Unsteady/solid_cylinder1.cfg), [solid_cylinder2.cfg](../../Inc_Heated_Cylinders_Unsteady/solid_cylinder2.cfg), [solid_cylinder3.cfg](../../Inc_Heated_Cylinders_Unsteady/solid_cylinder3.cfg)), the cofiguration file to invoke a multiphysics simulation run ([cht_2d_3cylinders.cfg](../../Inc_Heated_Cylinders_Unsteady/cht_2d_3cylinders.cfg)) and the mesh file ([mesh_cht_3cyl.su2](../../Inc_Heated_Cylinders_Unsteady/mesh_cht_3cyl.su2)).
20+
21+
## Tutorial
22+
23+
The following tutorial will walk you through the steps required when solving for an unsteady coupled CHT solution. It is assumed you have already obtained and compiled the SU2_CFD code for a serial computation. If you have yet to complete these requirements, please see the [Download](/docs/Download/) and [Installation](/docs/Installation/) pages and that make sure you have completed the [heated cylinders with conjugate heat transfer tutorial](/tutorials/Inc_Heated_Cylinders/).
24+
25+
### Background
26+
27+
For unsteady flows around walls that are transferring heat from an adjacent (solid) zone, the coupling of temperature and heat flux distributions has to be resolved for each and every time step. Both will vary over time as they depend on the current flow field.
28+
29+
### Problem Setup
30+
31+
The problem setup is the same as in the [heated cylinders with conjugate heat transfer tutorial](/tutorials/Inc_Heated_Cylinders/) except for the density. It is increased in all zones by a factor of 100 so that for the flow we obtain a Reynolds number of 4000 which will make it unsteady. Thus we set
32+
```
33+
INC_DENSITY_INIT= 0.0210322
34+
```
35+
in [flow_cylinder.cfg](../../Inc_Heated_Cylinders/flow_cylinder.cfg) and
36+
37+
```
38+
SOLID_DENSITY= 0.0210322
39+
```
40+
in [solid_cylinder1.cfg](../../Inc_Heated_Cylinders/solid_cylinder1.cfg), [solid_cylinder2.cfg](../../Inc_Heated_Cylinders/solid_cylinder2.cfg) and [solid_cylinder3.cfg](../../Inc_Heated_Cylinders/solid_cylinder3.cfg)
41+
42+
For simplicity we leave all other parameters unchanged.
43+
44+
### Mesh Description
45+
46+
The mesh is the same as in the [heated cylinders with conjugate heat transfer tutorial](/tutorials/Inc_Heated_Cylinders/).
47+
48+
### Configuration File Options
49+
50+
An unsteady simulation is set up by enabling the time domain and choosing a time marching algorithm in the [master config file](../../Inc_Heated_Cylinders/cht_2d_3cylinders.cfg):
51+
52+
```
53+
TIME_DOMAIN = YES
54+
%
55+
%
56+
TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER
57+
```
58+
59+
The time marching parameters have to match the flow physics that should be resolved. For a given inlet velocity of 3.40297 m/s at Re = 4000, the Strouhal number estimation for the most upstream cylinder is Sr = 0.21. This gives a frequency of f = Sr*v = 0.71Hz for the vortex shedding so that a time step of 0.05s is chosen in the [master config file](../../Inc_Heated_Cylinders/cht_2d_3cylinders.cfg):
60+
61+
```
62+
TIME_STEP= 0.05
63+
```
64+
65+
In order to sufficiently resolve the coupling in each time step, we set the number of outer iterations to 200 in the [master config file](../../Inc_Heated_Cylinders/cht_2d_3cylinders.cfg):
66+
67+
```
68+
OUTER_ITER = 200
69+
```
70+
71+
The number of inner (zone-internal) iterations is set to 1 by default. We do not have to touch any of the zone-specific config files for unsteady options.
72+
73+
### Running SU2
74+
75+
One time iteration will run rather quick and it is up to the user for how long the simulation should run or, equivalently, which physical time span should be covered. In the video above, 1000 time steps had been computed to generate a 50s realtime video. See the [heated cylinders with conjugate heat transfer tutorial](/tutorials/Inc_Heated_Cylinders/) how to execute SU2_CFD.

0 commit comments

Comments
 (0)