|
[Sponsors] |
chtMultiRegionFoam "thermoI.H" Not Converging |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Charlie Marshall
Join Date: Jan 2020
Posts: 17
Rep Power: 3 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Carlos Rubio Abujas
Join Date: Jan 2018
Location: Spain
Posts: 100
Rep Power: 6 ![]() |
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! |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Ashley Melvin
Join Date: Aug 2019
Posts: 1
Rep Power: 0 ![]() |
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.
|
|
![]() |
![]() |
![]() |
Tags |
chtmultiregionfoam, convergence, user error |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
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 |