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

chtMultiRegionSimpleFoam temperature goes to below zero

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 1, 2019, 04:26
Default chtMultiRegionSimpleFoam temperature goes to below zero
  #1
New Member
 
shach
Join Date: Apr 2019
Posts: 26
Rep Power: 6
shach934 is on a distinguished road
Hi Fomers,
I am using chtMultiRegionSimpleFoam to simulate cooling of a rotating disc. The geometry is attached. In a wind tunnel a disc is rotating, temperature of the surface of the disc is higher than the air. Here is the boundary condition and the thermophysicalProperties of air and disc.
For air 0/air/T
Code:
dimensions      [ 0 0 0 1 0 0 0 ];

internalField   uniform 300;

boundaryField
{
    inLet
    {
        type            fixedValue;
        value           $internalField;
    }
    rr_slave
    {
        type            cyclicAMI;
    }
    air
    {
        type            fixedValue;
        value           $internalField;
    }
    outLet
    {
        type            zeroGradient;
    }
    air2disc
    {
        type            compressible::turbulentTemperatureCoupledBaffleMixed;
        value           uniform 300;
        Tnbr            T;
        kappaMethod     fluidThermo;
    }
    rr
    {
        type            cyclicAMI;
    }
	
    air2disc_heat
    {	type            fixedValue;
        value           uniform 500;
	}
}
For disc 0/disc/T
Code:
dimensions      [ 0 0 0 1 0 0 0 ];

internalField   uniform 300;

boundaryField
{
    disc2air
    {
        type            compressible::turbulentTemperatureCoupledBaffleMixed;
        Tnbr            T;
        kappaMethod     solidThermo;
        value           uniform 300;
    }
	heat2air
    {
        type            fixedValue;
        value           uniform 500;
    }
}
The thermophysicalProperties for air:
Code:
thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    // coefficients for air

    specie
    {
        molWeight       28.85;
    }
    thermodynamics
    {
		Cp				1004.4;
		Hf				2.5e6;
	}
    transport
    {
       mu				1.8e-5;
	   Pr				0.705;
    }
}
For disc:
Code:
thermoType
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        molWeight   55;
    }

    transport
    {
        kappa   80;
    }

    thermodynamics
    {
        Hf      0;
        Cp      450;
    }

    equationOfState
    {
        rho     7800;
    }
}
After the chtMultiRegionSimpleFoam started a few steps, the minimum temperature in the air domain starts to decrease and very soon below zero, the solver exit with error.
Anyone has any suggestion? Thank!
Attached Images
File Type: jpg dis.jpg (20.7 KB, 19 views)
File Type: jpg tu.jpg (22.1 KB, 18 views)
shach934 is offline   Reply With Quote

Old   July 1, 2019, 05:17
Default
  #2
Member
 
Join Date: May 2013
Posts: 34
Rep Power: 12
carye is on a distinguished road
Have you checked the pressure and velocity boundary?
I use chtMultiRegionSimpleFoam also and have met the same problem before.
After I corrected the pressure boundary, it converged normally.

Another reason maybe the schemes you used.
Sometimes when I used a not such appropriate scheme, the calculation diverged soon.
carye is offline   Reply With Quote

Old   July 1, 2019, 05:59
Default
  #3
New Member
 
shach
Join Date: Apr 2019
Posts: 26
Rep Power: 6
shach934 is on a distinguished road
Hello Carye,
Thanks for your suggestion!
Would you mind to take a look at my pressure and velocity boundary conditions? Thanks!
Here is the pressure for the disc
Code:
dimensions      [ 1 -1 -2 0 0 0 0 ];

internalField   uniform 0;

boundaryField
{
    disc2air
    {
        type            calculated;
        value           $internalField;
    }
    heat2air
    {
        type            calculated;
        value           $internalField;
    }
}
Here are the pressure p boundary for air:
Code:
dimensions      [ 1 -1 -2 0 0 0 0 ];

internalField   uniform 1e5;

boundaryField
{
    inLet
    {
        type            zeroGradient;
    }
    rr_slave
    {
        type            cyclicAMI;
    }
    air
    {
        type            zeroGradient;
    }
    outLet
    {
        type            fixedValue;
        value           $internalField;
    }
    air2disc
    {
        type            zeroGradient;
    }
    rr
    {
        type            cyclicAMI;
    }
    air2disc_heat
    {
        type            zeroGradient;
    }
}
and p_rgh for air:
Code:
dimensions      [ 1 -1 -2 0 0 0 0 ];

internalField   uniform 1e5;

boundaryField
{
    inLet
    {
        type            zeroGradient;
    }
    rr_slave
    {
        type            cyclicAMI;
    }
    air
    {
        type            fixedFluxPressure;
        value           $internalField;
    }
    outLet
    {
        type            fixedValue;
        value           $internalField;
    }
    air2disc
    {
        type            fixedFluxPressure;
        value           $internalField;
    }
    rr
    {
        type            cyclicAMI;
    }
    air2disc_heat
    {
        type            fixedFluxPressure;
        value           $internalField;
    }
}
The velocity boundary conditions for the air is:
Code:
dimensions      [ 0 1 -1 0 0 0 0 ];

internalField   uniform ( 0 -6.944 0 );

boundaryField
{
    inLet
    {
        type            fixedValue;
        value           $internalField;
    }
    rr_slave
    {
        type            cyclicAMI;
    }
    air
    {
        type            noSlip;
    }
    outLet
    {
        type            inletOutlet;
        inletValue      uniform ( 0 0 0 );
        value           $internalField;
    }
    air2disc
    {
        type            noSlip;
    }
    rr
    {
        type            cyclicAMI;
    }
    air2disc_heat
    {
        type            noSlip;
    }
}
Well, there is no velocity for the disc, since i am using the MRF to simulate the rotation. Any suggestion? Thanks!
shach934 is offline   Reply With Quote

Old   July 1, 2019, 06:42
Default
  #4
Member
 
Join Date: May 2013
Posts: 34
Rep Power: 12
carye is on a distinguished road
I'm sorry that I'm not an expert.

But how about set the fixedFluxPressure in the p_rgh as zeroGradient?
I remember in my calculation I also tried fixedFluxPressure but it was unstable, so I used zeroGradient and finally got the result.

Another is that the p boundary is not used in the chtMultiRegionSimpleFoam because p is calculated from p_rgh.
So the pressure boundary will affect the calculation is p_rgh.

Hope this may help you.
carye is offline   Reply With Quote

Old   July 1, 2019, 08:12
Default
  #5
New Member
 
shach
Join Date: Apr 2019
Posts: 26
Rep Power: 6
shach934 is on a distinguished road
Hi Carye,
Thanks for your suggestion. I tried with zeroGradient boundary condition for p_rgh, unfortunately it is still not working. I will try with a simpler case now and check if this is caused by cyclicAMI boundary.
shach934 is offline   Reply With Quote

Old   July 1, 2019, 17:11
Default
  #6
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Just one hint... why do you use cyclicAMI if there is no motion?
Furthermore, if you get a few iterations, you can save them and check the results to see what field makes problems and especially at which location.


Can you let us know your fvSchemes and solver settings?
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   July 1, 2019, 17:54
Default
  #7
New Member
 
shach
Join Date: Apr 2019
Posts: 26
Rep Power: 6
shach934 is on a distinguished road
Hi Tobi,
The cyclicAMI is for the nonconformal mesh between the rotating region and the stator region. It is not the source of the problem since I tried without the MRF.
In fact, I already solved the problem. There is something wrong with the fvSolution, the pRefCell and pRefValue was set wrong, I changed it a little bit and it works....
Thanks.
shach934 is offline   Reply With Quote

Reply

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
Shadow Wall and temperature norger FLUENT 10 September 28, 2019 12:43
Static temperature and Total temperature Jong-Yoon FLUENT 0 September 25, 2016 02:42
How to get free stream temperature in boundary condition saharesobh FLUENT 0 October 9, 2012 18:12
difference of temperature HKH FLUENT 0 January 28, 2010 03:45
monitoring point of total temperature rogbrito FLUENT 0 June 21, 2009 18:31


All times are GMT -4. The time now is 09:18.