CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

convergence problem of chtMultiRegionFoam in simulation of free convection

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   May 1, 2017, 12:50
Default convergence problem of chtMultiRegionFoam in simulation of free convection
Join Date: Sep 2016
Posts: 63
Rep Power: 8
sitajeje is on a distinguished road
Dear Foamers,

I simulated free convection of air in a vertical channel. You can find simulation results of the thermal and viscous boundary layers in the attached pictures.

The simulated domain is a 125mm x 40 mm 2D rectangle with:
The front and back boundary as "empty",
The left and right boundary are "symmetryPlane"
The upper boundary is "pressureInletOutletVelocity" for velocity, "zeroGradient" for p_rgh and temperature,
The lower boundary has "inletOutlet" for velocity and temperature, and "totalPressure" for pressure.

A heater of 85mm tall in the vertical direction and 0.25mm thick in the horizontal direction is attached to the left boundary. fvOption was applied to add a heat source to the heater.

The simulation ran well before the maximum temperature was achieved, as shown in the pictures. I would expect the temperature to reach a steady state, but the temperature began to decrease, while the Courant number continued to increase to far above 1. As you can see in the attached animation, the thickness of the boundary layer shrank rapidly to almost zero.

The original mesh was 0.25 to 1.6mm/cell the horizontal direction, and 1.25mm/cell in the vertical direction. I double the mesh density, and the maximum temperature was achieved at the same runTime and had the same value. Afterwards the temperature decreased, and the boundary layer shrank in a similar manner although with a much slower rate.

I am not sure whether this computation deliver the correct results. Does anyone know why the temperature reduces after the maximal was achieved? I would appreciate for any comments and suggestions!

Thank you very much in advance!

Best regards,
Attached Images
File Type: jpg BLT.jpg (28.6 KB, 18 views)
File Type: jpg BLU.jpg (31.5 KB, 14 views)
Attached Files
File Type: zip (134.2 KB, 4 views)

Last edited by sitajeje; May 1, 2017 at 16:39.
sitajeje is offline   Reply With Quote

Old   May 1, 2017, 13:29
Arvind Jay
Join Date: Sep 2012
Posts: 96
Rep Power: 13
arvindpj is on a distinguished road
Is part of the left boundary wall? If so, the contours look correct.
arvindpj is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Natural convection energy convergence problem sanooj FLUENT 5 October 29, 2015 07:23
SimpleFoam convergence problem with really simple simulation mayank.dce2k7 OpenFOAM Running, Solving & CFD 2 November 19, 2013 05:28
Some problems with simulation in a Forced Convection problem Jordi P. Main CFD Forum 0 August 15, 2013 17:18
Urgent; convergence problem in MRF simulation Mansureh ANSYS 4 February 2, 2011 06:00
free convection heat transfer from a heated horizontal surface through a liquid to a thin cooled fin Kaushik FLUENT 1 May 8, 2000 06:47

All times are GMT -4. The time now is 05:29.