Heat transfer between multiple solids and fluids (snappyHexMesh)

Update: June 1, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater

Summary

We consider an analysis domain where the space is divided by a solid domain. Assuming that a heating element is set at the bottom of the solid domain, we calculate the temperature and flow of the entire system from 0 sec to 75 sec. The domain consists of five subdomains as shown below.

Model geometry (region heater) Model geometry (region "heater")
Model geometry (region leftSolid) Model geometry (region "leftSolid")
Model geometry (region rightSolid) Model geometry (region "rightSolid")
Model geometry (region bottomAir) Model geometry (region "bottomAir")
Model geometry (region topAir) Model geometry (region "topAir")

In the region "topAir", the air flows in from the X minimum plane at 0.1 m/s and out from the X maximum plane. The region "bottomAir" is assumed to be filled with air and there is no inflow or outflow of air.

The bottom of the region "heater" is fixed at a temperature of 500 K. The boundary temperature between the heater and the other regions is set to 300 K.

After the calculation, the utility "paraFoam" is used to generate .OpenFOAM files for each region. To visualize them, open each .OpenFOAM file from [File]-[Open] in the menu after starting ParaView.

Display all regions Display all regions

The meshes are as follows, and the number of mesh is 27404.

Meshes (region heater) Meshes (region "heater")
Meshes (region leftSolid) Meshes (region "leftSolid")
Meshes (region rightSolid) Meshes (region "rightSolid")
Meshes (region bottomAir) Meshes (region "bottomAir")
Meshes (region topAir) Meshes (region "topAir")

The calculation result is as follows.

Temperature in solid region (T) Temperature in solid region (T)
Temperature on XY-plane (T) Temperature on XY-plane (T)
Flow velocity on XY-plane (U) Flow velocity on XY-plane (U)

Commands

cp -r $FOAM_TUTORIALS/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater snappyMultiRegionHeater
cd snappyMultiRegionHeater

blockMesh
surfaceFeatureExtract
snappyHexMesh -overwrite
splitMeshRegions -cellZones -overwrite

rm -f 0*/heater/{nut,alphat,epsilon,k,U,p_rgh}
rm -f 0*/leftSolid/{nut,alphat,epsilon,k,U,p_rgh}
rm -f 0*/rightSolid/{nut,alphat,epsilon,k,U,p_rgh}

changeDictionary -region bottomAir
changeDictionary -region topAir
changeDictionary -region heater
changeDictionary -region leftSolid
changeDictionary -region rightSolid

chtMultiRegionFoam

paraFoam -touchAll
paraFoam

After ParaView has been started by the paraFoam command, each .OpenFOAM file must be opened from the menu [File]-[Open].

Calculation time

11 minutes 7.3 seconds *Single, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz