bloc.design_optim#

Functions#

find_2nd_inj_mass_flow_rate(gas_torch, qm_torch, ...)

Find the mass flow rate (in kg/s) of the second injection to reach the desired temperature in the reactor.

find_2nd_inj_mfr_Tin(gas_torch, qm_torch, gas_2nd_inj, ...)

Find the mass flow rate (in kg/s) of the 2nd injection to reach the desired temperature at the reactor entrance.

find_2nd_inj_mfr_Tout(gas_torch, qm_torch, ...[, ...])

Find the mass flow rate (in kg/s) of the second injection to reach the desired specific enthalpy at the reactor output.

find_quench_flow(qm_reac, gas_reac, X_quench, ...[, ...])

Compute the required quench mass flow for a desired mixed outlet temperature.

Module Contents#

bloc.design_optim.find_2nd_inj_mass_flow_rate(gas_torch, qm_torch, gas_2nd_inj, qm_2nd_inj_guess, T_reac_C, verbose=False)#

Find the mass flow rate (in kg/s) of the second injection to reach the desired temperature in the reactor.

This function does not take into account the chemical energy contained in the radicals produced by the torch. Thus, for high enthalpy plasma, the mass flow rate of the second injection will be underestimated.

Parameters:

  • gas_torch (-) – Gas object of the torch

  • qm_torch (-) – Mass flow rate of the torch in kg/s

  • gas_2nd_inj (-) – Gas object of the reactor input

  • qm_2nd_inj_guess (-) – Initial guess for the mass flow rate of the second injection in kg/s

  • T_reac_C (-) – Desired temperature in the reactor in °C

  • verbose (-) – If True, print verbose output.

Returns:

- qm_2nd_inj – Mass flow rate of the second injection in kg/s

Return type:

float

bloc.design_optim.find_2nd_inj_mfr_Tin(gas_torch, qm_torch, gas_2nd_inj, qm_2nd_inj_guess, T_in_target, t_res_PSR, verbose=False)#

Find the mass flow rate (in kg/s) of the 2nd injection to reach the desired temperature at the reactor entrance.

Simulate a PSR for t_res_PSR, and adapt qm_2nd_inj to get T_in_target after t_res_PSR in the PSR.

Parameters:

  • gas_torch (-) – Gas object of the torch

  • qm_torch (-) – Mass flow rate of the torch in kg/s

  • gas_2nd_inj (-) – Gas object of the reactor input

  • qm_2nd_inj_guess (-) – Initial guess for the mass flow rate of the second injection in kg/s

  • T_in_target (-) – Desired temperature at the reactor entrance in °C

  • t_res_PSR (-) – Residence time in the PSR in s.

  • verbose (-) – If True, print verbose output.

Returns:

- qm_2nd_inj – Mass flow rate of the second injection in kg/s

Return type:

float

bloc.design_optim.find_2nd_inj_mfr_Tout(gas_torch, qm_torch, gas_2nd_inj, qm_2nd_inj_guess, T_out_target, C_solid='C(s)', rescale_factor=1000.0, verbose=False)#

Find the mass flow rate (in kg/s) of the second injection to reach the desired specific enthalpy at the reactor output.

This enthalpy target is defined as the enthalpy of fully converted CH4 at T_out_target. Since the specific enthalpy depends on the initial CH4/H2 ratio, we define hp = ( h - Y_H2*h_H2(T_target) ) / Y_CH4 to remove the dependance of the enthalpy target with the initial CH4/H2 ratio.

Note: this function assumes that only H2 and CH4 are the only species injected in the torch and reactor.

Parameters:

  • gas_torch (-) – Gas object of the torch

  • qm_torch (-) – Mass flow rate of the torch in kg/s

  • gas_2nd_inj (-) – Gas object of the reactor input

  • qm_2nd_inj_guess (-) – Initial guess for the mass flow rate of the second injection in kg/s

  • T_out_target (-) – Desired temperature at the reactor exit in °C

  • verbose (-) – If True, print verbose output.

Returns:

- qm_2nd_inj – Mass flow rate of the second injection in kg/s

Return type:

float

bloc.design_optim.find_quench_flow(qm_reac, gas_reac, X_quench, T_quench, T_target, verbose=False)#

Compute the required quench mass flow for a desired mixed outlet temperature.

Assumptions - Ideal-gas enthalpy model (true for Cantera ideal-gas phases): h depends

only on temperature and species; mixture h(T) = sum(Y_k h_k(T)).

  • Adiabatic, no shaft work, no pressure change during mixing.

Closed-form solution:

qm_quench = qm_reac * (h1(T1) - h1(Tt)) / (h2(Tt) - h2(T2))

Parameters:

  • qm_reac (-) – Mass flow rate of the reactor in kg/s

  • gas_reac (-) – Gas object of the reactor output

  • X_quench (-) – Composition of the quench flow in the reactor input

  • T_quench (-) – Temperature of the quench flow in °C

  • T_target (-) – Desired temperature after the quench in °C

  • verbose (-) – If True, print verbose output.

Returns:

- qm_quench – Mass flow rate of the quench flow in kg/s

Return type:

float