ThermalCoupling.hxx

Go to the documentation of this file.

/*
  This file is part of MURE,
  Copyright (C) 2007-2021 MURE developers.
 
  MURE is free software: you can redistribute it and/or modify
  it under the terms of the GNU Lesser General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  MURE is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU Lesser General Public License for more details.
 
  You should have received a copy of the GNU Lesser General Public License
  along with MURE.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef _THERMAL_
#define _THERMAL_
 
#include <vector>
#include <string>
#include <fstream>
#include <algorithm>
 
#include "ThermalDataReader.hxx"
#include "Cell.hxx"
 
class Cell;
class ThermalDataReader;
 
 
class ThermalCoupling
{
  public :
    ThermalCoupling(
        int NumberOfLevels,     
        int ZoneNumber = 0,         
        int IterationNumber = 0     
    );
    ~ThermalCoupling(); 
    void Run(); 
 
 
    void SetCosinePowerDeposits()
    {
        fCosinePower = true;   
    }
    void SetRodRadius(double RodRadius)
    {
        fFuelRadius = RodRadius;   
    }
    void SetInnerRodRadius(double InnerRodRadius)
    {
        fInnerFuelRadius = InnerRodRadius;   
    }
    void SetCladRadius(double CladRadius)
    {
        fCladRadius = CladRadius;   
    }
    void SetCellsHeight(double ElementaryCellHeight)
    {
        fCellHeight = ElementaryCellHeight;   
    }
    void SetCellsSurface(double ElementaryCellSurface)
    {
        fCellSurface = ElementaryCellSurface;   
    }
    void SetEnterCoolantTemp(double EnterCoolantTemp)
    {
        fCoolantEnterTemp.push_back(EnterCoolantTemp);   
    }
    void SetCoolantType(string Coolant)
    {
        fCoolant = Coolant;   
    }
    void SetFuelType(string Fuel)
    {
        fFuel = Fuel;   
    }
    void SetCladdingType(string Cladding)
    {
        fCladding = Cladding;   
    }
    void SetRodAveragePower(double Power)
    {
        fRodAveragePower = Power;   
    }
    void SetMassSpeed(double MassSpeed)
    {
        fMassSpeed = MassSpeed;   
    }
    void SetHexLattice()
    {
        fHexLattice = true;   
    }
    void SetFuelPorosity(double FuelPorosity)
    {
        fFuelPorosity = FuelPorosity;   
    }
    void CalculateHeatTransferInGazSpace()
    {
        fHeatTransferGazSpace = true;   
    }
    void SetDeltaTempOfGazSpace(double DeltaTemp)
    {
        fDeltaTempOfGazSpace = DeltaTemp;    
        CalculateHeatTransferInGazSpace();
    }
    void SetHeatExchangeCoeffOfGap(double Hgap)
    {
        fHGap = Hgap;    
        CalculateHeatTransferInGazSpace();
    }
 
 
  protected :
    void CheckConsistent(); 
    void GetAllCells(); 
    void LoadingData(); 
    void CalculateRodPowerDeposits(); 
    void CalculateCosinusPowerDeposits();
    void PrintPowerResults(); 
    void FluidHeating(); 
    void DimensionlessNumbersCalculation(); 
    void Convection(); 
    void CladConduction(); 
    void ThermalRadiation(); // Resolves Thermal radiation phenomena between cladding and fuel
    void FuelConduction(); 
    void PressureLosses(); 
    void PrintThermalResults(); 
    void CloseFiles(); 
    void UploadingOfNewProperties(); 
    void ForcedDestroy(); 
    //double Dittus(double, double, double); //!< The Dittus-Boelter equation to calculate the heat transfer coefficient
 
  private :
    int             fi; 
    int             fLevels; 
    int             fRadialZone; 
    double          fFuelRadius; 
    double          fInnerFuelRadius;   
    double          fCladRadius; 
    // double           fGazSpace; //!< Void space between fuel rod and cladding : Fission Products space
    double          fCellHeight; 
    double          fCellSurface; 
    string          fCoolant; 
    string          fFuel; 
    string          fCladding; 
    double          fMassSpeed; 
    double          fRodAveragePower; 
 
    bool            fCosinePower; 
    bool            fSigmaPhiUpdated; 
    bool            fHexLattice; 
    bool            fHeatTransferGazSpace; 
    double          fDeltaTempOfGazSpace; 
    double          fHGap; 
    double          fFuelPorosity; 
    double          fPassageSurface; 
    double          fDh; 
    vector < double >   fCp; 
    vector < double >   fVis; 
    vector < double >   fk; 
    vector < double >   fFuelk; 
    vector < double >   fPr; 
    vector < double >   fRe; 
    vector < double >   fNu; 
    vector < double >   fCoolantEnterTemp;  
    vector < double >   fCoolantExitTemp;   
    vector < double >   fPower; 
    vector < double >   fPowerDensities; 
    vector < double >   fFuelTemps; 
    vector < double >   fFuelTp; 
    vector < double >   fFuelTi; 
    vector < double >   fCladdingTemps; 
    vector < double >   fCladdingTp; 
    vector < double >   fCladdingTi; 
    vector < double >   fCoolantTemps; 
    vector < double >   fCoolantEntryTemps; 
    vector < double >   fCoolantExitTemps; 
    vector < double >   fCoolantDensities; 
    vector < double >   fParietalStream; 
    vector < double >   fPressureLosses; 
 
    vector < Cell *>    fFuelCells; 
    vector < Cell *>    fCladdingCells; 
    vector < Cell *>    fCoolantCells; 
 
    vector < ofstream *> fFiles;        
    vector < ofstream *> fFuelOuter;    
    vector < ofstream *> fFuelInner;    
    vector < ofstream *> fCoolantDensity; 
 
    ThermalDataReader *fThermalData;
};
 
#endif