|
| 1 | +--- |
| 2 | +title: Unsteady Conjugate Heat Transfer |
| 3 | +permalink: /tutorials/Inc_Heated_Cylinders_Unsteady/ |
| 4 | +written_by: oleburghardt |
| 5 | +for_version: 7.1.0 |
| 6 | +revised_by: TobiKattmann |
| 7 | +revision_date: 2021-03-04 |
| 8 | +revised_version: 7.1.1 |
| 9 | +solver: INC_NAVIER_STOKES, HEAT_EQUATION |
| 10 | +requires: SU2_CFD |
| 11 | +complexity: advanced |
| 12 | +follows: Static_CHT |
| 13 | +userguide: Multizone |
| 14 | +--- |
| 15 | + |
| 16 | +<br> |
| 17 | +<iframe width="560" height="315" src="https://www.youtube.com/embed/MqN8GalSyzk" frameborder="0" allow="accelerometer; autoplay; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe> |
| 18 | + |
| 19 | +## Goals |
| 20 | + |
| 21 | +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. |
| 22 | +The following capabilities of SU2 will be showcased in this tutorial: |
| 23 | + |
| 24 | +- Time domain and time-marching config file options (plus related ones) for unsteady simulations |
| 25 | +- Use of time iterations, outer and inner iterations |
| 26 | +- Paraview multiblock output |
| 27 | + |
| 28 | +The intent of this tutorial is to demonstrate how a steady CHT simulation can be turned into an unsteady one. |
| 29 | + |
| 30 | +## Resources |
| 31 | + |
| 32 | +The resources for this tutorial can be found in the [Inc_Heated_Cylinders_Unsteady](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht) directory in the [tutorial repository](https://github.com/su2code/Tutorials). You will need the configuration files for all physical zones ([flow_cylinder.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/flow_cylinder.cfg), [solid_cylinder1.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/solid_cylinder1.cfg), [solid_cylinder2.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/solid_cylinder2.cfg), [solid_cylinder3.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/solid_cylinder3.cfg)), the cofiguration file to invoke a multiphysics simulation run ([cht_2d_3cylinders.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/cht_2d_3cylinders.cfg)) and the mesh file ([mesh_cht_3cyl.su2](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/mesh_cht_3cyl.su2)). |
| 33 | + |
| 34 | +## Tutorial |
| 35 | + |
| 36 | +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/). |
| 37 | + |
| 38 | +### Background |
| 39 | + |
| 40 | +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. |
| 41 | + |
| 42 | +### Problem Setup |
| 43 | + |
| 44 | +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 |
| 45 | +``` |
| 46 | +INC_DENSITY_INIT= 0.0210322 |
| 47 | +``` |
| 48 | +in [flow_cylinder.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/flow_cylinder.cfg) and |
| 49 | + |
| 50 | +``` |
| 51 | +MATERIAL_DENSITY= 0.0210322 |
| 52 | +``` |
| 53 | +in [solid_cylinder1.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/solid_cylinder1.cfg), [solid_cylinder2.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/solid_cylinder2.cfg) and [solid_cylinder3.cfg](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/solid_cylinder3.cfg) |
| 54 | + |
| 55 | +For simplicity we leave all other parameters unchanged. |
| 56 | + |
| 57 | +### Mesh Description |
| 58 | + |
| 59 | +The [mesh](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/mesh_cht_3cyl.su2) is the same as in the [heated cylinders with conjugate heat transfer tutorial](/tutorials/Inc_Heated_Cylinders/). |
| 60 | + |
| 61 | +### Configuration File Options |
| 62 | + |
| 63 | +An unsteady simulation is set up by enabling the time domain and choosing a time marching algorithm in the [master config file](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/cht_2d_3cylinders.cfg): |
| 64 | + |
| 65 | +``` |
| 66 | +TIME_DOMAIN = YES |
| 67 | +% |
| 68 | +% |
| 69 | +TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER |
| 70 | +``` |
| 71 | + |
| 72 | +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](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/cht_2d_3cylinders.cfg): |
| 73 | + |
| 74 | +``` |
| 75 | +TIME_STEP= 0.05 |
| 76 | +``` |
| 77 | + |
| 78 | +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](https://github.com/su2code/Tutorials/tree/master/multiphysics/unsteady_cht/cht_2d_3cylinders.cfg): |
| 79 | + |
| 80 | +``` |
| 81 | +OUTER_ITER = 200 |
| 82 | +``` |
| 83 | + |
| 84 | +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. |
| 85 | + |
| 86 | +### Running SU2 |
| 87 | + |
| 88 | +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