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

Help Required for chtMultiRegionSimpleFoam

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

Like Tree2Likes
  • 1 Post By Alczem
  • 1 Post By yips

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 30, 2023, 09:52
Default Help Required for chtMultiRegionSimpleFoam
  #1
New Member
 
Yiping Song
Join Date: May 2023
Posts: 7
Rep Power: 2
yips is on a distinguished road
Briefly, I am testing chtMultiRegionSimpleFoam for a simple case of conjugated heat flow but I could not get acceptable results. Below is a description of my test case and the problem I encountered:

The case is a simple one of an aluminium body with a coolant channel.
Top surface of the Al body has a heat input, I started with a heat flux and then I tried fixed gradient. Both resulted in the same issue.
Coolant has a fixed velocity at input and fixed pressure at output.
Simulation works fine but the results are not correct.

Below are the results and issues (see attached file for images):
o Pressure in the coolant developed OK and it looks correct.
o Flow velocity developed OK and looks correct.
o Heat flow within the Al body looks correct.
o The problem I have is the energy conservation does not work:
 Outlet coolant temperature is far too high for the energy input from the Al body.
 It feels as if in temperature calculation, coolant is treated as stationery.


So, to summaries:
Everything else looks correct in the simulation, except for temperature of the coolant. It behaves as if it is not flowing at all when temperature is calculated.

Any help and suggestions are welcome.

Below are my temperature boundary files for your reference:

(1) For the solid:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0/solid";
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{
solid_top
{
// type fixedValue;
// type externalWallHeatFluxTemperature;
// mode flux;
// q uniform 1e4; // heat flux [W/m^2]
// kappaMethod solidThermo;
// value uniform 320;

type fixedGradient;
gradient uniform 50;

}


solid_walls
{
type zeroGradient;
}

solid_to_liquid
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod fluidThermo;
thicknessLayers (1e-5);
kappaLayers (0.024);

}
}

(2) for the liquid:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{

liquid_inlet
{
type fixedValue;
value uniform 300;
}
liquid_outlet
{
type inletOutlet;
inletValue uniform 300;
value uniform 300;
}
liquid_to_solid
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod fluidThermo;
thicknessLayers (1e-5);
kappaLayers (0.024);
}

}
Attached Files
File Type: pdf chtMultiRegionSimpleFoam issue reduced size.pdf (147.0 KB, 3 views)
yips is offline   Reply With Quote

Old   May 31, 2023, 05:26
Default
  #2
Senior Member
 
Join Date: Dec 2021
Posts: 140
Rep Power: 3
Alczem is on a distinguished road
Hey!


I think there is a slight mistake in the definition of your boundary condition solid_to_liquid, where you set kappaMethod to fluidThermo. It should be solidThermo if I am not mistaken.


If you can share your whole case, it would be helpful to pinpoint other potential issue
yips likes this.
Alczem is offline   Reply With Quote

Old   May 31, 2023, 05:51
Default Sharing the whole case
  #3
New Member
 
Yiping Song
Join Date: May 2023
Posts: 7
Rep Power: 2
yips is on a distinguished road
Hi Alczem,

Thanks for your reply. I will give it a try.

When you mentioned about sharing the whole case, what are the files needed to be included? (I am quite new to this forum, so I need to learn.)

Yiping
yips is offline   Reply With Quote

Old   May 31, 2023, 06:05
Default set kappaMethod to fluidThermo
  #4
New Member
 
Yiping Song
Join Date: May 2023
Posts: 7
Rep Power: 2
yips is on a distinguished road
Hi Alczem,

To clarify your suggestion:

for solid_to_liquid, I should set kappaMethod to solidThermo
for liquid_to_solid, I should keep kappaMethod to fluidThermo

Is above understanding correct?

Yiping
yips is offline   Reply With Quote

Old   June 1, 2023, 11:12
Default
  #5
Senior Member
 
Join Date: Dec 2021
Posts: 140
Rep Power: 3
Alczem is on a distinguished road
Quote:
To clarify your suggestion:

for solid_to_liquid, I should set kappaMethod to solidThermo
for liquid_to_solid, I should keep kappaMethod to fluidThermo

Is above understanding correct?

Yes!




To share your case, if it is not confidential, you can create a zip archive containing the 0, constant and system folders and the scripts you may be using to run the case. Don't include the mesh if it is too big.
Alczem is offline   Reply With Quote

Old   June 5, 2023, 08:49
Default Zip file attached
  #6
New Member
 
Yiping Song
Join Date: May 2023
Posts: 7
Rep Power: 2
yips is on a distinguished road
I have uploaded the case file. Sorry I have to remove all mesh data due to size limitation.
Attached Files
File Type: zip CFD-Online_case_file.zip (19.3 KB, 1 views)
yips is offline   Reply With Quote

Old   June 5, 2023, 10:46
Default liquid temperature is still too high
  #7
New Member
 
Yiping Song
Join Date: May 2023
Posts: 7
Rep Power: 2
yips is on a distinguished road
Hi Alczem, I changed the T boundary file in the 0 directory for the solid/liquid interface to solidThermo:

solid_to_liquid
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod solidThermo;
thicknessLayers (1e-5);
kappaLayers (0.024);

}

This does not look to have made a lot of difference. Only after 6000 iteration, liquid temperature has become too high and feels like it is not flowing. See attached. On the other hand, velocity results show that the liquid is flowing nicely.
Attached Images
File Type: jpg Image 05062023.JPG (101.5 KB, 8 views)
yips is offline   Reply With Quote

Old   June 6, 2023, 04:57
Default
  #8
Senior Member
 
Join Date: Dec 2021
Posts: 140
Rep Power: 3
Alczem is on a distinguished road
Hey


I took a quick look at the case. A couple of things:

  • Relaxation factor for h for both the fluid and the solid is very low, try to get as close to 1 as possible (try 0.9, then 0.99 etc). It made a clear difference in several of my cases to achieve convergence.
  • for p file, you should have calculated as the type for every boundary.
  • you could try turbulentTemperatureRadCoupledMixed instead of turbulentTemperatureCoupledBaffleMixed for the coupled patches. I have used almost exclusively this one, the setups are similar.
  • Are you sure your flow is laminar? Turbulence can play a significant role for heat exchanges.
Other than that, I cannot see a blatant mistake, especially if the case runs ok despite the unconclusive results. Try to set a flowrate rather than a fixedValue for the velocity at the inlet, play with different numerical schemes, and double check your mesh.


Last thing you could do is to run a transient case rather than a steadystate case. It is usually more stable, and to speed up the computation, you can divide the heat capacity of your materials by 10 or 100 to reach a steady state quicker. Can't help you much more, sorry! Keep the thread updated if you manage to solve your issues
Alczem is offline   Reply With Quote

Old   June 6, 2023, 05:09
Default
  #9
New Member
 
Yiping Song
Join Date: May 2023
Posts: 7
Rep Power: 2
yips is on a distinguished road
Thanks for your suggestions. I will try them and see if I can get something sensible. May take some time before I get somewhere but I will post it if anything works.

For relaxation factor, I took a value from some examples I found. I do not really know what does it do, so I made no attempt to change it. I will try your suggestion. What is it by the way?
yips is offline   Reply With Quote

Old   June 6, 2023, 09:41
Default
  #10
New Member
 
Yiping Song
Join Date: May 2023
Posts: 7
Rep Power: 2
yips is on a distinguished road
Sorry for asking question about relaxation factor prematurely. I have found information from the OpenFoam web site now.

Thanks.
Alczem likes this.
yips is offline   Reply With Quote

Reply

Tags
chtmultiregionsimplefoam

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
Install openFOAM from with proxy server Dhruval OpenFOAM Installation 3 October 18, 2014 15:38
paraview installation woes vex OpenFOAM Installation 15 January 30, 2011 07:11
[OpenFOAM] Problem with paraFoam on a linux-64 bit bunni ParaView 4 April 14, 2010 20:55
Version 15 on Mac OS X gschaider OpenFOAM Installation 113 December 2, 2009 10:23
paraFoam reader for OpenFOAM 1.6 smart OpenFOAM Installation 13 November 16, 2009 21:41


All times are GMT -4. The time now is 23:28.