CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

chtMultiRegionFoam "thermoI.H" Not Converging

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 13, 2020, 17:09
Post chtMultiRegionFoam "thermoI.H" Not Converging
New Member
Charlie Marshall
Join Date: Jan 2020
Posts: 17
Rep Power: 3
charlliemarshalll is on a distinguished road
OpenFOAM 7 v1912

I was trying to replicate the OpenFOAM validation case in this paper:

"Conjugate heat transfer of backward-facing step flow: A benchmark problem revisited" by Matjazˇ Ramšak

Although I have double-checked everything, I cannot get the case converge. I have checked all the tutorial cases to borrow the configurations, but none worked.

The case is uploaded. Just run the Allrun file.

The error log produced by chtMultiRegionFoam is as follows:

Solving for fluid region water
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for h, Initial residual = 1, Final residual = 1.2406221e-11, No Iterations 1

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const, bool) const [with Thermo = Foam::hConstThermo<Foam::rhoConst<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::rho Const<Foam::specie> >, Foam::sensibleEnthalpy>]
in file /home/ubuntu/OpenFOAM/OpenFOAM-7/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 70
Energy -> temperature conversion failed to converge:
iter Test e/h Cv/p Tnew
0 0 0 4181 -1.4382206e-22
1 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
2 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
3 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
4 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
5 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
6 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
7 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
8 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
9 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
10 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
11 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
12 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
13 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
14 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22
15 -1.4382206e-22 -6.0132005e-19 4181 -1.4382206e-22


Maximum number of iterations exceeded: 100


in file /home/ubuntu/OpenFOAM/OpenFOAM-7/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 106.
Attached Files
File Type: gz back1.tar.gz (4.6 KB, 5 views)
charlliemarshalll is offline   Reply With Quote

Old   May 21, 2020, 13:03
Default chtMultiRegionFoam not valid for adimensional numbers
Senior Member
Carlos Rubio Abujas
Join Date: Jan 2018
Location: Spain
Posts: 100
Rep Power: 6
crubio.abujas is on a distinguished road
Hi charlliemarshalll,

It seems that you are applying unrealistic temperatures to the problem, varying the T field between 0 and 1. I guess that the paper have adimensional number and you're trying to supply the same conditions directly. All the thermal properties are apparently introduced as S.I variables: solid_cp 450 J/kg? solid_rho=8000 kg/m3? It seems that they are rather S.I values or too extreme adimensional ones.

I've tried using more reasonable T values [300K for the water an 600K for the wall] and not found any problem. You can put numbers into the paper formulas and solve for an specific case or try to convert the thermo properties into adminesional ones.

Hope it helps!
crubio.abujas is offline   Reply With Quote

Old   February 21, 2021, 17:54
New Member
Ashley Melvin
Join Date: Aug 2019
Posts: 1
Rep Power: 0
ashansellum is on a distinguished road
I might be late, but I think the error might be because of the initial temperature field. You could try with a non-zero internalField for temperature for both the heater & water. Try to keep the internalField value between 0 & 1 for faster convergence.
ashansellum is offline   Reply With Quote


chtmultiregionfoam, convergence, user error

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
Error in thermophysical properties (chtMultiRegionFoam) mukut OpenFOAM Pre-Processing 27 October 5, 2017 06:48
Error in chtMultiRegionFoam michael157 OpenFOAM Running, Solving & CFD 17 May 22, 2017 03:32
FOAM FATAL IO ERROR for chtMultiRegionFoam xiaoyoyo OpenFOAM Running, Solving & CFD 0 May 8, 2012 16:49
Embed explicitSetValue in chtMultiRegionFoam samiam1000 OpenFOAM 2 April 18, 2012 05:14
chtMultiRegionFoam laminar heat exchanger non converging Pierpaolo OpenFOAM 7 October 26, 2009 09:06

All times are GMT -4. The time now is 10:53.