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

chtMultiRegionFoam: Solid region not heating up

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 19, 2021, 17:18
Default chtMultiRegionFoam: Solid region not heating up
  #1
New Member
 
Edgar Alejandro Martínez Ojeda
Join Date: Jul 2019
Posts: 20
Rep Power: 5
Edgar Alejandro Martínez is on a distinguished road
Greetings,

I have a 2.5D case of a fluid past a square cylinder with heat transfer (images attached below). The mesh is done with blockMeshDict.

There's an inlet, outlet, the sides are declared as walls and the top and bottom are declared as symmetry planes (there are top_fluid and top_solid, same goes for bottom_fluid and bottom_solid).

The solid is kept at 500K (supposedly), the inlet velocity is water at 1m/s.

I used topoSetDict to create the regions and indeed you can visualize them in paraView, then I used splitMeshRegions -cellZones -overwrite to create the fluid-solid interfaces and they are created as well (although I can't see the patches in paraView).

The PROBLEM is that the solid region drops its temperature dramatically soon after the simulation begins (it reaches 300K which is the fluid's temperature). The fluid barely absorbs some heat and its temperature does not change dramatically (289K-302K).

I've tried small time steps but nothing works, the square cylinder is 10 X 10 cm and I really doubt it can be cooled to ambient temperature rapidly despite the fluid being water.

I've also tried modifying the scalarSemiImplicitSource (supposedly it has no time duration and I think it's supposed to heat permanently) but results don't change.

Here's the 0/fluid file for temperature:

Code:
 
dimensions      [ 0 0 0 1 0 0 0 ];

internalField   uniform 300;

boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }

    "(left|right)"
    {
        type            zeroGradient;
    }

    fluid_to_solid
    {
        type            compressible::turbulentTemperatureCoupledBaffleMixed;
        value           $internalField;
        Tnbr            T;
        kappaMethod     fluidThermo;
    }

	"(top_fluid|bottom_fluid)"
	{
		type			symmetryPlane;
	}
}
Here's the 0/solid file for T:

Code:
dimensions      [ 0 0 0 1 0 0 0 ];

internalField   uniform 500;

boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    solid_to_fluid
    {
        type            compressible::turbulentTemperatureCoupledBaffleMixed;
        value           $internalField;
        Tnbr            T;
        kappaMethod     solidThermo;
    }

	"(top_solid|bottom_solid)"
	{
		type			symmetryPlane;
	}
}
Here's the fvOptions file:

Code:
options
{
	mySource
	{
        type    scalarSemiImplicitSource;
        active  true;

        scalarSemiImplicitSourceCoeffs
        {
        	selectionMode   all;
            cellSet      	c0;
            volumeMode      absolute;
            injectionRateSuSp
            {
            	T	(500 0);
            }
        }

	}
}
Thank you for your help and time.
Attached Images
File Type: png 1.png (31.3 KB, 3 views)
File Type: jpg 2.jpg (41.8 KB, 4 views)
File Type: jpg 3.jpg (48.9 KB, 3 views)
Attached Files
File Type: zip square_cylinder.zip (103.1 KB, 5 views)
Edgar Alejandro Martínez is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, heat transfer, square cylinder

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
Solid-Fluid region interface mesh for chtMultiRegionFoam ishan_ae OpenFOAM Running, Solving & CFD 4 February 8, 2021 13:09
[solids4Foam] How to calculate drag coeff when using solids4Foam amuzeshi OpenFOAM CC Toolkits for Fluid-Structure Interaction 15 November 7, 2019 13:50
alphat Wien3 OpenFOAM Running, Solving & CFD 22 August 15, 2019 11:00
problem with Min/max rho tH3f0rC3 OpenFOAM 8 July 31, 2019 10:48
Induction Motor in chtMultiRegionFoam xsa OpenFOAM Running, Solving & CFD 0 December 28, 2016 07:16


All times are GMT -4. The time now is 03:19.