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
Returns:

dimemsionless quantity of Qadvective
Return type:

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
Pram k:

conductivity
Returns:

Qond, dimemsionless quantity
Return type:

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
Returns:

Qond, dimemsionless quantity
Return type:

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
Returns:

dimensionless quantity of Qconv
Return type:

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
Returns:

value of dtemp/dt
Return type:

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