Skip to content

Commit 3ddf707

Browse files
davidscnuekermanBenjaminRodenberg
authored
Restructure flow over plate (#154)
* Add a first draft version * Add cleaning script and calculix draft * Adjust precice config and openfoam cases * Remove non-existing mesh tag * Tune cases to new structure * Make Nutils case consistent (one mesh, 2D) * Update README for flow-over-plate * Remove Calculix case * Clarify Vynnycky reference description Co-authored-by: Benjamin Rodenberg <[email protected]> * Rename heat.py to solid.py for fenics * Update precice config in order to run in parallel * Remove calculix from the README * Add physical parameters to text * Remove converToMeters in OpenFOAM blockMesh * Use correct inlet velocity * Clearify on Prandtl number Co-authored-by: Benjamin Uekermann <[email protected]> Co-authored-by: Benjamin Rodenberg <[email protected]>
1 parent 5d78546 commit 3ddf707

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

64 files changed

+442
-146361
lines changed

flow-over-heated-plate/README.md

+59-2
Original file line numberDiff line numberDiff line change
@@ -1,10 +1,67 @@
11
---
22
title: Flow over heated plate
33
permalink: tutorials-flow-over-heated-plate.html
4-
keywords: tutorial
5-
summary:
4+
keywords: tutorial, CHT, conjugate-heat transfer, OpenFOAM, FEniCS, Nutils
5+
summary: This tutorial describes how to run a conjugate heat transfer coupled simulation using preCICE and any fluid-solid solver combination of our [officially provided adapter codes](adapters-overview.html).
66
---
77

88
{% include important.html content="We have not yet ported the documentation of the preCICE tutorials from the preCICE wiki to here. Please go to the [preCICE wiki](https://github.com/precice/precice/wiki#2-getting-started---tutorials)" %}
99

10+
11+
## Setup
12+
13+
This scenario consists of one fluid and one solid participant and it is inspired by Vynnycky et al. [1]. A fluid enters a channel with temperature `T_\infty`, where it comes in contact with a solid plate. The plate is heated at its bottom and has a constant temperature of `T_hot`.
14+
1015
![img](images/tutorials-flow-over-heated-plate-example.png)
16+
17+
The test case is two-dimensional and a serial-implicit coupling with Aitken underrelaxation is used for the coupling.
18+
19+
The inlet velocity is `u_\infty = 0.1 m/s`, the inlet temperature is `T_\infty = 300K`. The fluid and solid have the same thermal conductivities `k_S = k_F = 100 W/m/K`. Further material properties of the fluid are its viscosity `mu = 0.0002 kg/m/s` and specific heat capacity `c_p = 5000 J/kg/K`. The Prandtl number `Pr = 0.01` follows from thermal conductivity, viscosity and specific heat capacity. The solid has the thermal diffusivity `\alpha = 1 m^2/s`. The gravitational acceleration is `g = 9.81 m/s^2`.
20+
21+
## Available solvers
22+
23+
By default, the fluid participant reads heat-flux values and the solid participant reads temperature values for the coupled simulation. The following participants are available:
24+
25+
Fluid participant:
26+
27+
* OpenFOAM (buoyantPimpleFoam). For more information, have a look at the [OpenFOAM adapter documentation](adapter-openfoam-overview.html).
28+
29+
Solid participant:
30+
31+
* OpenFOAM (laplacianFoam). For more information, have a look at the [OpenFOAM adapter documentation](adapter-openfoam-overview.html).
32+
33+
* FEniCS. For more information, have a look at the [FeniCS adapter documentation](adapter-fenics.html).
34+
35+
## Running the Simulation
36+
37+
All listed solvers can be used in order to run the simulation. Open two separate terminals and start the desired fluid and solid participant by calling the respective run script `run.sh` located in the participant directory. For example:
38+
39+
```
40+
cd fluid-openfoam
41+
./run.sh
42+
```
43+
and
44+
```
45+
cd solid-fenics
46+
./run.sh
47+
```
48+
in order to use OpenFOAM and FEniCS for this test case.
49+
50+
## Post-processing
51+
52+
How to visualize the simulation results depends on the selected solvers. Most of the solvers generate `vtk` files which can visualized using, e.g., ParaView.
53+
An example of the visualized expected results looks as follows:
54+
55+
![result](images/result-openfoam.png)
56+
57+
Observe that the temperature at the bottom of the plate is 310K and at the inlet 300K. On the interface, the temperature is between these values. An area of higher temperature is formed above the plate, which is shifted towards the front, driven by the flow.
58+
59+
You may use additional filters, such as the Calculator and the Plot Over Line, to obtain the distribution of the non-dimensional temperature (T-T_inlet)/(T_solid-T_inlet):
60+
61+
![graph](images/graph-result.png)
62+
63+
## References
64+
65+
[1] M. Vynnycky, S. Kimura, K. Kanev, and I. Pop. Forced convection heat transfer from a flat plate: the conjugate problem. International Journal of Heat and Mass Transfer, 41(1):45 – 59, 1998.
66+
67+
{% include disclaimer.html content="This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software via www.openfoam.com, and owner of the OPENFOAM® and OpenCFD® trade marks." %}

flow-over-heated-plate/buoyantPimpleFoam-fenics/.gitignore

-9
This file was deleted.

flow-over-heated-plate/buoyantPimpleFoam-fenics/Allclean

-61
This file was deleted.

flow-over-heated-plate/buoyantPimpleFoam-fenics/Allrun

-47
This file was deleted.

flow-over-heated-plate/buoyantPimpleFoam-fenics/Allrun_parallel

-58
This file was deleted.

flow-over-heated-plate/buoyantPimpleFoam-fenics/Fluid/constant/polyMesh/boundary

-71
This file was deleted.

0 commit comments

Comments
 (0)