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

My *first* multiregion case

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 25, 2021, 15:01
Default trying to debug - running into brick wall
  #41
Senior Member
 
Alan w
Join Date: Feb 2021
Posts: 277
Rep Power: 6
boffin5 is on a distinguished road
In a very methodical way, I am trying to debug my multiRegion case. It is composed of an air region, which incorporates a simple fairing around the radiator, and a porous region, which is the radiator itself consisting of a porous zone.


My first step was to modify the case to eliminate the radiator. In the process, I cleaned up the fairing mesh; it's good now. When I ran chtMultiRegionFoam to 200 seconds, it finished normally. Partial Success!


Then, I did the same for the other zone. I deleted the radiator fairing, such that the air zone is simply a blockMesh domain. So now, the only mesh being analyzed is the radiator itself, hanging in air. When I launched the solver, it immediately failed with my old friend, the message stating: energy -> temperature conversion failed to converge.


This debug version is basically identical to the heatExchanger tutorial. I closely compared the two, which only differ only in the porous zone being modelled, but I can't find any reason why my case fails, and the tutorial doesn't. As they say in France, I am gobsmacked.


Attached are my runlog file, and a zip file of the case with the radiator only.



I am 99.9% through my OF self course, and am being held up by this mystical problem. Please, if anyone has ideas, I would be delighted to hear them.
Attached Files
File Type: zip rad-hx-only.zip (67.4 KB, 4 views)
File Type: txt runlog.txt (13.5 KB, 2 views)
boffin5 is offline   Reply With Quote

Old   December 16, 2023, 14:26
Default my *first* multiregion case - The Saga is over (I hope)! plus a question
  #42
Senior Member
 
Alan w
Join Date: Feb 2021
Posts: 277
Rep Power: 6
boffin5 is on a distinguished road
A year and half ago, I posted my first plea for help with a chtMultiRegion case. Over the ensuing months, I received many posts of help from the smart people in this forum. Finally Yann provided the assistance that broke the log jam. Attached is an image of the temperature flow that has me celebrating. (I hope there is no other error lurking like the one that that fooled me for months.)

This case concerns a cuboid shaped radiator in a duct. Many of the suggestions said, just look in the tutorials. However it turns out that my method of setting up the case has no equivalent in the heat transfer tutorials. These use toposet, but Yann said I don't need it; all I have to do is to make sure the radiator lines up properly with the duct, and define it that way in the meshes. And it worked. It seems to me that my method is much more representative of real world scenarios.

Many times, I have seen posts from students in Formula SAE or Formula Student expressing the desire to model radiator installations for their cars, I imagine that many of them had the same problems as I did. In my opinion, the tutorials, while invaluable, can be incomplete.

Here is my question: In my successful case, after having problems with stl files, I modelled the radiator with blockMeshDict. But I thought that I should also be able to do it with salome, exporting the mesh as a unv file. When I specifed the face names (e.g. inlet, outlet, etc.), ideasUnvToFoam repeatedly failed. When I made no attempt to define face names, it worked, but in paraView the image is bizarre (see attached). Ultimately, using a unv file is probably the best way, but as of now, I can't make it work.

Thanks again to Yann and all the others who ventured to help me.
Attached Images
File Type: png flow-temperature.png (28.8 KB, 11 views)
File Type: png unv-file-image.png (27.6 KB, 11 views)
Yann likes this.
boffin5 is offline   Reply With Quote

Old   December 17, 2023, 09:28
Default
  #43
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,098
Rep Power: 26
Yann will become famous soon enough
Hello Alan,

Thanks for your kind words, I'm glad you finally managed to reach your goal.

You should be able to mesh your radiator with either snappy, blockMesh or Salome, it doesn't really matter.
I won't be of any help with Salome since I don't use it, but I'm pretty sure there are people around here who can help.
Yann is offline   Reply With Quote

Old   January 2, 2024, 09:11
Default chtMultiRegionFoam kappa setting goes wrong
  #44
New Member
 
yingting tang
Join Date: Aug 2023
Posts: 8
Rep Power: 2
tyting is on a distinguished road
I'm going to simulate electronicBoard(material:Al) cooling process with water,when I setting the metal's thermalphysicalProperties,it remind me that:
[0] --> FOAM FATAL IO ERROR:
[0] Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 43 the label 205
[0]
[0] file: /home/zh_zhanghpc/OpenFOAM/ytTang/eletricalHT/constant/metal/thermophysicalProperties/mixture/transport/kappa at line 43.
[0]
[0] From function [4]
[4]
[4] --> FOAM FATAL IO ERROR:
[4] Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 0 the label 205
[4]
[4] file: IOstream/mixture/transport/kappa at line 0.[5]
[5]
[5] --> FOAM FATAL IO ERROR:
[5] Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 0 the label 205
[5]
[5] file: IOstream/mixture/transport/kappa at line 0.
[5]
[5] From function [6]
[6]
[6] --> FOAM FATAL IO ERROR:
[6] Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 0 the label 205
[6]
[6] file: IOstream/mixture/transport/kappa at line 0.
[6]
[6] From function Foam::Istream& Foam::Istream::readBegin(const char*)

[4]
[4] From function Foam::Istream& Foam::Istream::readBegin(const char*)
[4] in file db/IOstreams/IOstreams/Istream.C at line 92.
[4]
FOAM parallel run exiting
[4]
Foam::Istream& Foam::Istream::readBegin(const char*)
[5] in file db/IOstreams/IOstreams/Istream.C at line 92.
[5]
FOAM parallel run exiting
[5]
[6] in file db/IOstreams/IOstreams/Istream.C at line 92.
[6]
FOAM parallel run exiting
[6]
[7]
[7] file: IOstream/mixture/transport/kappa at line 0.
[7]
[7] From function Foam::Istream& Foam::Istream::readBegin(const char*)
[7] in file db/IOstreams/IOstreams/Istream.C at line 92.
[7]
FOAM parallel run exiting
[7]
Foam::Istream& Foam::Istream::readBegin(const char*)
[0] in file db/IOstreams/IOstreams/Istream.C at line 92.
[0]
FOAM parallel run exiting
[0]
[1]
[1]
[1] --> FOAM FATAL IO ERROR:
[1] Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 0 the label 205
[1]
[1] [2]
[2]
[2] --> FOAM FATAL IO ERROR:
[2] Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 0 the label 205
[2]
[2] file: IOstream/mixture/transport/kappa at line file: IOstream/mixture/transport/kappa at line 0.
[1]
[1] 0.
[2]
[2] From function Foam::Istream& Foam::Istream::readBegin(const char*)
[2] in file db/IOstreams/IOstreams/Istream.C at line 92 From function Foam::Istream& Foam::Istream::readBegin(const char*)
[1] in file db/IOstreams/IOstreams/Istream.C at line 92.
[1]
FOAM parallel run exiting
[1]
.
[2]
FOAM parallel run exiting
[2]
[3]
[3]
[3] --> FOAM FATAL IO ERROR:
[3] Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 0 the label 205
[3] --------------------------------------------------------------------------
MPI_ABORT was invoked on rank 6 in communicator MPI COMMUNICATOR 3 SPLIT FROM 0
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

[3] file: IOstream/mixture/transport/kappa at line 0.
[3]
[3] From function Foam::Istream& Foam::Istream::readBegin(const char*)
[3] in file db/IOstreams/IOstreams/Istream.C at line 92.
[3]
FOAM parallel run exiting
[3]
[DESKTOP-91GJLA5:28832] 7 more processes have sent help message help-mpi-api.txt / mpi-abort
[DESKTOP-91GJLA5:28832] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

my setup for metal:
thermoType
{
type heSolidThermo;
mixture pureMixture;
transport constAnIsoSolid;
thermo eConst;
equationOfState rhoConst;
specie specie;
energy sensibleInternalEnergy;
}

mixture
{
specie
{
nMoles 1;
molWeight 26.98;

}
thermodynamics
{
Cv 700;
Hf 0;
}
transport
{
kappa 205;
}
equationOfState
{
rho 2.7;
}
}

Is anyone know how can I fixed it? Any help will be appreciated!!Thanks!!
tyting is offline   Reply With Quote

Reply


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
OF 4.0 multiregion case calculate yPlus pbnuclex OpenFOAM Post-Processing 6 July 16, 2020 04:27
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 14:53
[OpenFOAM] ParaView 4.10 and OpenFOAM 2.3.0 Multiregion and decomposed case romant ParaView 3 April 7, 2014 15:42
Transient case running with a super computer microfin FLUENT 0 March 31, 2009 11:20
Turbulent Flat Plate Validation Case Jonas Larsson Main CFD Forum 0 April 2, 2004 10:25


All times are GMT -4. The time now is 01:30.