CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

chtMultiRegionFoam "thermoI.H" Not Converging

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

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

Background:
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

Problem:
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.

Info:
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

....

--> FOAM FATAL ERROR:
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, 2 views)
charlliemarshalll is offline   Reply With Quote

Old   May 21, 2020, 13:03
Default chtMultiRegionFoam not valid for adimensional numbers
  #2
Member
 
Carlos Rubio Abujas
Join Date: Jan 2018
Location: Spain
Posts: 34
Rep Power: 4
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

Reply

Tags
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 08:29.