# pyrk.th_system module¶

class pyrk.th_system.THSystem(kappa, components)[source]

Bases: object

This class models:

• lumped capacitance model for slab geometry
• 1-D heat diffusion in spherical geometry, for heat generation at any radius in the sphere,
• advective heat transfer by fluid
• convective heat transfer at the solid surface to fluid.
BC_center(component, t_idx)[source]

Volumetric conductive heat flux Qconduction from the center of a sphere to the first boundary (conduction without interface) in watts/meter**3

Parameters: component (str) – name of the inner most component(mesh element) t_idx (int) – timestep at which Qconduction is calculated

:return : Qcondction :rtype:float, dimensionless

advection(component, t_idx, t_in, m_flow, cp)[source]

calculate heat transfer by advection in watts

Parameters: component (str) – name of the component t_idx (int) – time step that conduction heat is computed t_in (float) – inlet temperature m_flow (float) – mass flow rate though the control volume cp (float) – specific heat capacity of the fluid dimemsionless quantity of Qadvective float
comp_from_name(name)[source]

Returns the component with the matching name

conductionFVM(component, env, t_idx, L=<Quantity(0.0, 'meter')>, k=<Quantity(0.0, 'meter')>, A=<Quantity(0.0, 'meter ** 2')>)[source]

compute volumetric conductive heat transfer by conduction(watts/m3)

Parameters: component (str) – name of the component env (str) – name of the environment t_idx (int) – time step that conduction heat is computed L (float, units meter) – conduction distance conductivity Qond, dimemsionless quantity float
conduction_slab(component, env, t_idx, L, A)[source]

compute volumetric heat transfer by conduction(watts/m3)

Parameters: component (str) – name of the component env (str) – name of the environment t_idx (int) – time step that conduction heat is computed Qond, dimemsionless quantity float
convBoundary(component, t_b, t_env, h, R)[source]

calculate heat transfer through convective boundray condition for the mesh element at the surface of the spherical Supercomponent (watts)

Parameters: component (str) – name of the outer most solid component t_b (float) – temperature of the body t_env (float) – temperature of the environment(fluid) h (float) – convective heat transfer coefficient R (float) – outer radius of the component dimensionless quantity of Qconv float
convection(t_b, t_env, h, A)[source]

heat transfer by convection(watts)

Parameters: t_b (float.) – The temperature of the body t_env (float.) – The temperature of the environment h (float.) – the heat transfer coefficient between environment and body A (float.) – the surface area of heat transfer
custom(t_b, t_env, res)[source]
dtempdt(component, power, omegas, t_idx)[source]

compute the derivative dtemperature/dt

Parameters: component (str) – name of the component that dtemp/dt is calculated power (float) – nuclear power density omegas (list) – decay heat nuclear data t_idx (int) – the timestep that dtempdt is calculated for value of dtemp/dt float, kelvin/s
heatgen(component, power, omegas)[source]

calculate heat transfer by conduction(watts/m3)

Parameters: component (str) – name of the component power – normalized nuclear power generated in the component

i.e.: power*power_tot = power (in watts) :type power: float

mass_trans(t_b, t_inlet, H, u)[source]
Parameters: t_b (float.) – The temperature of the body t_inlet – The temperature of the flow inlet
metadata(component)[source]

a recorder function that calls down to each component. used for the th/th_params table

record(component)[source]

a recorder function that calls down to each component. used for the th/th_timeseries table