This object performs the variation of temperature in a given region for a specific system. More...

#include <ThermalCoupling.hxx>

Public Member Functions
Constructor &destructor
	ThermalCoupling (int NumberOfLevels, int ZoneNumber=0, int IterationNumber=0)
	Default constructor.

	~ThermalCoupling ()
	Delete standard.

void	Run ()
	start the calculation

Miscellaneous methods
int	fi
	Iteration step number.

int	fLevels
	Number of axial levels.

int	fRadialZone
	Radial zone number.

double	fFuelRadius
	Radius of fuel rod cylinder (= radius of cladding internal crown)

double	fInnerFuelRadius
	Inner Radius of fuel rod cylinder (annular pellet)

double	fCladRadius
	Radius of exterior cylinder (cladding radius)

double	fCellHeight
	Height of each cell cylinder.

double	fCellSurface
	Surface of the cell (coolant+cladd+fuel)

string	fCoolant
	Name of the coolant (h2o, d2o, sodium, ...)

string	fFuel
	Name of the fuel (uox, mox, thpu, ...)

string	fCladding
	Name of the cladding (zircaloy, steel, ...)

double	fMassSpeed
	Mass speed that verifies the conservation of the mass (rho*V=cte=fMassSpeed in kg/m2s)

double	fRodAveragePower
	Rod average power (W)

bool	fCosinePower
	default value = false ; if true : power deposits type cosinus

bool	fSigmaPhiUpdated
	Used to check if the update of sigma phi is already done.

bool	fHexLattice
	Used for flow area calculation.

bool	fHeatTransferGazSpace
	Used for a heat transfer exchange inside the gaz space.

double	fDeltaTempOfGazSpace
	Value of temperature difference between cladding and fuel inside the gaz space (given by user)

double	fHGap
	Value of the heat exchange coefficient inside the gap (W/m2/K)

double	fFuelPorosity
	default value = 0 ; if set 0 < < 1 by user, will be given to ThermalDataReader for density correction

double	fPassageSurface
	Surface of passage (surface of coolant in an under shannel)

double	fDh
	Hydraulic diameter (4*S/P)

vector< double >	fCp
	Heat capacity value of each coolant cell [J/kg.K].

vector< double >	fVis
	Viscosities values of each coolant cell [Pa*s].

vector< double >	fk
	Thermal conductivities values of each coolant cell [W/(m.K)].

vector< double >	fFuelk
	Thermal conductivities values of each fuel cell [W/(m.K)].

vector< double >	fPr
	Prandtl Numbers of each coolant cell.

vector< double >	fRe
	Reynolds Numbers of each coolant cell.

vector< double >	fNu
	Nusselt Numbers of each coolant cell.

vector< double >	fCoolantEnterTemp
	Coolant entrance temperature in cells.

vector< double >	fCoolantExitTemp
	Coolant exit temperature in cells.

vector< double >	fPower
	Power deposited in each fuel cell [W].

vector< double >	fPowerDensities
	Power densities of each fuel cell [W/m3].

vector< double >	fFuelTemps
	Temperatures of fuel cells.

vector< double >	fFuelTp
	Parietal Temperature of fuel cells.

vector< double >	fFuelTi
	Interior Temperature of fuel cells.

vector< double >	fCladdingTemps
	Temperatures of cladding cells.

vector< double >	fCladdingTp
	Parietal Temperature of cladding cells.

vector< double >	fCladdingTi
	Interior Temperature of cladding cells.

vector< double >	fCoolantTemps
	Average temperatures of coolant cells.

vector< double >	fCoolantEntryTemps
	Temperatures of entrance in coolant cells.

vector< double >	fCoolantExitTemps
	Temperatures of exit in coolant cells.

vector< double >	fCoolantDensities
	Densities of coolant cells.

vector< double >	fParietalStream
	Parietal Stream see by each cell [W/m2].

vector< double >	fPressureLosses
	Losses of pressure.

vector< Cell * >	fFuelCells
	Fuel cells implicated in the calculation.

vector< Cell * >	fCladdingCells
	Claddings cells implicated in the calculation.

vector< Cell * >	fCoolantCells
	Coolants cells implicated in the calculation.

vector< ofstream * >	fFiles
	Generates output files for thermal data (Axial averages and radial distribution in cells.

vector< ofstream * >	fFuelOuter
	Generates output files for outer pellet temperature.

vector< ofstream * >	fFuelInner
	Generates output files for inner pellet temperature.

vector< ofstream * >	fCoolantDensity
	Generates output files for coolant density.

ThermalDataReader *	fThermalData

void	SetCosinePowerDeposits ()

void	SetRodRadius (double RodRadius)

void	SetInnerRodRadius (double InnerRodRadius)

void	SetCladRadius (double CladRadius)

void	SetCellsHeight (double ElementaryCellHeight)

void	SetCellsSurface (double ElementaryCellSurface)

void	SetEnterCoolantTemp (double EnterCoolantTemp)

void	SetCoolantType (string Coolant)

void	SetFuelType (string Fuel)

void	SetCladdingType (string Cladding)

void	SetRodAveragePower (double Power)

void	SetMassSpeed (double MassSpeed)

void	SetHexLattice ()

void	SetFuelPorosity (double FuelPorosity)

void	CalculateHeatTransferInGazSpace ()

void	SetDeltaTempOfGazSpace (double DeltaTemp)

void	SetHeatExchangeCoeffOfGap (double Hgap)

void	CheckConsistent ()
	Checks the consistent with all data.

void	GetAllCells ()
	Gets fuel, cladding and coolants cells.

void	LoadingData ()
	Update of each cell data of the previous step (temperature, densities, ...)

void	CalculateRodPowerDeposits ()
	Calculate power densities of each fuel cell in a SINGLE rod.

void	CalculateCosinusPowerDeposits ()
	Calculates each power density with a cosine distribution.

void	PrintPowerResults ()
	Prints results of power deposits.

void	FluidHeating ()
	Calculates the warm-up of the fluid.

void	DimensionlessNumbersCalculation ()
	Calculate the dimensionless numbers (Re, Nu)

void	Convection ()
	Calculates the convective heat transfert between rod and coolant.

void	CladConduction ()
	Resolves Fourier equation of conduction in the claddind.

void	ThermalRadiation ()

void	FuelConduction ()
	Resolves Fourier equation of conduction in the fuel.

void	PressureLosses ()
	Calculates losses of pressure due to linear phenomena (no consideration of the term due to the peculiarities)

void	PrintThermalResults ()
	Prints radial thermal results for each cell.

void	CloseFiles ()
	Closes output files created.

void	UploadingOfNewProperties ()
	Updates news properties of cells and materials.

void	ForcedDestroy ()
	Explicit destruction of this thermalcoupling object.

Detailed Description

This object performs the variation of temperature in a given region for a specific system.

At present, the geometry must be cylindrical for the fuel and the cladding It solve the heat equation in the case of a stationnary state Estimates done : heat conduction in homogeneous isotropic media -> T(r) only no thermal radiation -> T(ext fuel) = T(int cladd) The flag SetEvolution() is needed to the calculation of power deposits The calculation is done on an average cell : no crossflow. The first cell MUST be at the bottom of the fuel rod

Author: Nico*

Version: 2.0

Constructor & Destructor Documentation

◆ ThermalCoupling()

ThermalCoupling::ThermalCoupling	(	int	NumberOfLevels,
		int	ZoneNumber = `0`,
		int	IterationNumber = `0`
	)

Default constructor.

Parameters

NumberOfLevels	Number of axial levels
ZoneNumber	Radial zone number if differentiation by region is required (different universes)
IterationNumber	Iteration step

◆ ~ThermalCoupling()

ThermalCoupling::~ThermalCoupling ( )

default

Delete standard.

Member Function Documentation

◆ CalculateCosinusPowerDeposits()

void ThermalCoupling::CalculateCosinusPowerDeposits ( )

protected

Calculates each power density with a cosine distribution.

◆ CalculateHeatTransferInGazSpace()

void ThermalCoupling::CalculateHeatTransferInGazSpace ( )

inline

< Force a calculation of the heat transfers in the gaz space

◆ CalculateRodPowerDeposits()

void ThermalCoupling::CalculateRodPowerDeposits ( )

protected

Calculate power densities of each fuel cell in a SINGLE rod.

◆ CheckConsistent()

void ThermalCoupling::CheckConsistent ( )

protected

Checks the consistent with all data.

◆ CladConduction()

void ThermalCoupling::CladConduction ( )

protected

Resolves Fourier equation of conduction in the claddind.

◆ CloseFiles()

void ThermalCoupling::CloseFiles ( )

protected

Closes output files created.

◆ Convection()

void ThermalCoupling::Convection ( )

protected

Calculates the convective heat transfert between rod and coolant.

◆ DimensionlessNumbersCalculation()

void ThermalCoupling::DimensionlessNumbersCalculation ( )

protected

Calculate the dimensionless numbers (Re, Nu)

◆ FluidHeating()

void ThermalCoupling::FluidHeating ( )

protected

Calculates the warm-up of the fluid.

◆ ForcedDestroy()

void ThermalCoupling::ForcedDestroy ( )

protected

Explicit destruction of this thermalcoupling object.

◆ FuelConduction()

void ThermalCoupling::FuelConduction ( )

protected

Resolves Fourier equation of conduction in the fuel.

◆ GetAllCells()

void ThermalCoupling::GetAllCells ( )

protected

Gets fuel, cladding and coolants cells.

◆ LoadingData()

void ThermalCoupling::LoadingData ( )

protected

Update of each cell data of the previous step (temperature, densities, ...)

◆ PressureLosses()

void ThermalCoupling::PressureLosses ( )

protected

Calculates losses of pressure due to linear phenomena (no consideration of the term due to the peculiarities)

◆ PrintPowerResults()

void ThermalCoupling::PrintPowerResults ( )

protected

Prints results of power deposits.

◆ PrintThermalResults()

void ThermalCoupling::PrintThermalResults ( )

protected

Prints radial thermal results for each cell.

◆ Run()

void ThermalCoupling::Run ( )

start the calculation

◆ SetCellsHeight()

void ThermalCoupling::SetCellsHeight ( double ElementaryCellHeight )

inline

< Sets height of cells (the same for each one)

◆ SetCellsSurface()

void ThermalCoupling::SetCellsSurface ( double ElementaryCellSurface )

inline

< Sets the total surface of the elementary cell (coolant+cladding+fuel)

◆ SetCladdingType()

void ThermalCoupling::SetCladdingType ( string Cladding )

inline

< Sets cladding type by giving its name

◆ SetCladRadius()

void ThermalCoupling::SetCladRadius ( double CladRadius )

inline

< Sets cladding radius

◆ SetCoolantType()

void ThermalCoupling::SetCoolantType ( string Coolant )

inline

< Sets coolant type by giving its name

◆ SetCosinePowerDeposits()

void ThermalCoupling::SetCosinePowerDeposits ( )

inline

< Allows a calculation without real coupling -> cosinus distribution deposits

◆ SetDeltaTempOfGazSpace()

void ThermalCoupling::SetDeltaTempOfGazSpace ( double DeltaTemp )

inline

< Fix an arbitrary value of the temperature difference inside the gaz space

◆ SetEnterCoolantTemp()

void ThermalCoupling::SetEnterCoolantTemp ( double EnterCoolantTemp )

inline

< Sets the entrance coolant temperature (K)

◆ SetFuelPorosity()

void ThermalCoupling::SetFuelPorosity ( double FuelPorosity )

inline

< sets fuel porosity to be given to ThermalDataReader for density correction

◆ SetFuelType()

void ThermalCoupling::SetFuelType ( string Fuel )

inline

< Sets fuel type by giving its name

◆ SetHeatExchangeCoeffOfGap()

void ThermalCoupling::SetHeatExchangeCoeffOfGap ( double Hgap )

inline

< Fix the heat exchange coefficient inside the gap space (W/m2/K)

◆ SetHexLattice()

void ThermalCoupling::SetHexLattice ( )

inline

< Information of hexagonal lattice (used for flow area calculation)

◆ SetInnerRodRadius()

void ThermalCoupling::SetInnerRodRadius ( double InnerRodRadius )

inline

< Sets fuel pellet inner radius (annular pellet)

◆ SetMassSpeed()

void ThermalCoupling::SetMassSpeed ( double MassSpeed )

inline

< Mass speed [kg/m2s] (e.g. in PWR mass speed = 3900 kg/m2s)

◆ SetRodAveragePower()

void ThermalCoupling::SetRodAveragePower ( double Power )

inline

< Sets the average power of ONE rod (Watt)

◆ SetRodRadius()

void ThermalCoupling::SetRodRadius ( double RodRadius )

inline

< Sets fuel pellet radius

◆ ThermalRadiation()

void ThermalCoupling::ThermalRadiation ( )

protected

◆ UploadingOfNewProperties()

void ThermalCoupling::UploadingOfNewProperties ( )

protected

Updates news properties of cells and materials.

Member Data Documentation

◆ fCellHeight

double ThermalCoupling::fCellHeight

private

Height of each cell cylinder.

◆ fCellSurface

double ThermalCoupling::fCellSurface

private

Surface of the cell (coolant+cladd+fuel)

◆ fCladding

string ThermalCoupling::fCladding

private

Name of the cladding (zircaloy, steel, ...)

◆ fCladdingCells

vector< Cell *> ThermalCoupling::fCladdingCells

private

Claddings cells implicated in the calculation.

◆ fCladdingTemps

vector< double > ThermalCoupling::fCladdingTemps

private

Temperatures of cladding cells.

◆ fCladdingTi

vector< double > ThermalCoupling::fCladdingTi

private

Interior Temperature of cladding cells.

◆ fCladdingTp

vector< double > ThermalCoupling::fCladdingTp

private

Parietal Temperature of cladding cells.

◆ fCladRadius

double ThermalCoupling::fCladRadius

private

Radius of exterior cylinder (cladding radius)

◆ fCoolant

string ThermalCoupling::fCoolant

private

Name of the coolant (h2o, d2o, sodium, ...)

◆ fCoolantCells

vector< Cell *> ThermalCoupling::fCoolantCells

private

Coolants cells implicated in the calculation.

◆ fCoolantDensities

vector< double > ThermalCoupling::fCoolantDensities

private

Densities of coolant cells.

◆ fCoolantDensity

vector< ofstream *> ThermalCoupling::fCoolantDensity

private

Generates output files for coolant density.

◆ fCoolantEnterTemp

vector< double > ThermalCoupling::fCoolantEnterTemp

private

Coolant entrance temperature in cells.

◆ fCoolantEntryTemps

vector< double > ThermalCoupling::fCoolantEntryTemps

private

Temperatures of entrance in coolant cells.

◆ fCoolantExitTemp

vector< double > ThermalCoupling::fCoolantExitTemp

private

Coolant exit temperature in cells.

◆ fCoolantExitTemps

vector< double > ThermalCoupling::fCoolantExitTemps

private

Temperatures of exit in coolant cells.

◆ fCoolantTemps

vector< double > ThermalCoupling::fCoolantTemps

private

Average temperatures of coolant cells.

◆ fCosinePower

bool ThermalCoupling::fCosinePower

private

default value = false ; if true : power deposits type cosinus

◆ fCp

vector< double > ThermalCoupling::fCp

private

Heat capacity value of each coolant cell [J/kg.K].

◆ fDeltaTempOfGazSpace

double ThermalCoupling::fDeltaTempOfGazSpace

private

Value of temperature difference between cladding and fuel inside the gaz space (given by user)

◆ fDh

double ThermalCoupling::fDh

private

Hydraulic diameter (4*S/P)

◆ fFiles

vector< ofstream *> ThermalCoupling::fFiles

private

Generates output files for thermal data (Axial averages and radial distribution in cells.

◆ fFuel

string ThermalCoupling::fFuel

private

Name of the fuel (uox, mox, thpu, ...)

◆ fFuelCells

vector< Cell *> ThermalCoupling::fFuelCells

private

Fuel cells implicated in the calculation.

◆ fFuelInner

vector< ofstream *> ThermalCoupling::fFuelInner

private

Generates output files for inner pellet temperature.

◆ fFuelk

vector< double > ThermalCoupling::fFuelk

private

Thermal conductivities values of each fuel cell [W/(m.K)].

◆ fFuelOuter

vector< ofstream *> ThermalCoupling::fFuelOuter

private

Generates output files for outer pellet temperature.

◆ fFuelPorosity

double ThermalCoupling::fFuelPorosity

private

default value = 0 ; if set 0 < < 1 by user, will be given to ThermalDataReader for density correction

◆ fFuelRadius

double ThermalCoupling::fFuelRadius

private

Radius of fuel rod cylinder (= radius of cladding internal crown)

◆ fFuelTemps

vector< double > ThermalCoupling::fFuelTemps

private

Temperatures of fuel cells.

◆ fFuelTi

vector< double > ThermalCoupling::fFuelTi

private

Interior Temperature of fuel cells.

◆ fFuelTp

vector< double > ThermalCoupling::fFuelTp

private

Parietal Temperature of fuel cells.

◆ fHeatTransferGazSpace

bool ThermalCoupling::fHeatTransferGazSpace

private

Used for a heat transfer exchange inside the gaz space.

◆ fHexLattice

bool ThermalCoupling::fHexLattice

private

Used for flow area calculation.

◆ fHGap

double ThermalCoupling::fHGap

private

Value of the heat exchange coefficient inside the gap (W/m2/K)

◆ fi

int ThermalCoupling::fi

private

Iteration step number.

◆ fInnerFuelRadius

double ThermalCoupling::fInnerFuelRadius

private

Inner Radius of fuel rod cylinder (annular pellet)

◆ fk

vector< double > ThermalCoupling::fk

private

Thermal conductivities values of each coolant cell [W/(m.K)].

◆ fLevels

int ThermalCoupling::fLevels

private

Number of axial levels.

◆ fMassSpeed

double ThermalCoupling::fMassSpeed

private

Mass speed that verifies the conservation of the mass (rho*V=cte=fMassSpeed in kg/m2s)

◆ fNu

vector< double > ThermalCoupling::fNu

private

Nusselt Numbers of each coolant cell.

◆ fParietalStream

vector< double > ThermalCoupling::fParietalStream

private

Parietal Stream see by each cell [W/m2].

◆ fPassageSurface

double ThermalCoupling::fPassageSurface

private

Surface of passage (surface of coolant in an under shannel)

◆ fPower

vector< double > ThermalCoupling::fPower

private

Power deposited in each fuel cell [W].

◆ fPowerDensities

vector< double > ThermalCoupling::fPowerDensities

private

Power densities of each fuel cell [W/m3].

◆ fPr

vector< double > ThermalCoupling::fPr

private

Prandtl Numbers of each coolant cell.

◆ fPressureLosses

vector< double > ThermalCoupling::fPressureLosses

private

Losses of pressure.

◆ fRadialZone

int ThermalCoupling::fRadialZone

private

Radial zone number.

◆ fRe

vector< double > ThermalCoupling::fRe

private

Reynolds Numbers of each coolant cell.

◆ fRodAveragePower

double ThermalCoupling::fRodAveragePower

private

Rod average power (W)

◆ fSigmaPhiUpdated

bool ThermalCoupling::fSigmaPhiUpdated

private

Used to check if the update of sigma phi is already done.

◆ fThermalData

ThermalDataReader* ThermalCoupling::fThermalData

private

◆ fVis

vector< double > ThermalCoupling::fVis

private

Viscosities values of each coolant cell [Pa*s].

The documentation for this class was generated from the following files:

/gpr/meplan/SMURE/source/include/ThermalCoupling.hxx
/gpr/meplan/SMURE/source/src/ThermalCoupling.cxx

Public Member Functions

Miscellaneous methods

Detailed Description

Constructor & Destructor Documentation

◆ ThermalCoupling()

◆ ~ThermalCoupling()

Member Function Documentation

◆ CalculateCosinusPowerDeposits()

◆ CalculateHeatTransferInGazSpace()

◆ CalculateRodPowerDeposits()

◆ CheckConsistent()

◆ CladConduction()

◆ CloseFiles()

◆ Convection()

◆ DimensionlessNumbersCalculation()

◆ FluidHeating()

◆ ForcedDestroy()

◆ FuelConduction()

◆ GetAllCells()

◆ LoadingData()

◆ PressureLosses()

◆ PrintPowerResults()

◆ PrintThermalResults()

◆ Run()

◆ SetCellsHeight()

◆ SetCellsSurface()

◆ SetCladdingType()

◆ SetCladRadius()

◆ SetCoolantType()

◆ SetCosinePowerDeposits()

◆ SetDeltaTempOfGazSpace()

◆ SetEnterCoolantTemp()

◆ SetFuelPorosity()

◆ SetFuelType()

◆ SetHeatExchangeCoeffOfGap()

◆ SetHexLattice()

◆ SetInnerRodRadius()

◆ SetMassSpeed()

◆ SetRodAveragePower()

◆ SetRodRadius()

◆ ThermalRadiation()

◆ UploadingOfNewProperties()

Member Data Documentation

◆ fCellHeight

◆ fCellSurface

◆ fCladding

◆ fCladdingCells

◆ fCladdingTemps

◆ fCladdingTi

◆ fCladdingTp

◆ fCladRadius

◆ fCoolant

◆ fCoolantCells

◆ fCoolantDensities

◆ fCoolantDensity

◆ fCoolantEnterTemp

◆ fCoolantEntryTemps

◆ fCoolantExitTemp

◆ fCoolantExitTemps

◆ fCoolantTemps

◆ fCosinePower

◆ fCp

◆ fDeltaTempOfGazSpace

◆ fDh

◆ fFiles

◆ fFuel

◆ fFuelCells

◆ fFuelInner

◆ fFuelk

◆ fFuelOuter

◆ fFuelPorosity

◆ fFuelRadius

◆ fFuelTemps

◆ fFuelTi

◆ fFuelTp

◆ fHeatTransferGazSpace

◆ fHexLattice

◆ fHGap

◆ fi

◆ fInnerFuelRadius