boundary_enthalpy_from_temperature Function

private function boundary_enthalpy_from_temperature(mesh, params, temperature, Y_point) result(h_value)

Convert a boundary temperature to enthalpy using the active thermo model. Boundary enthalpy from a boundary temperature.

When a transported boundary composition is available, the fixed-temperature boundary enthalpy is evaluated as h(T_b,Y_b,p0). This keeps hot/cold species inlets thermodynamically consistent.

Arguments

Type IntentOptional Attributes Name
type(mesh_t), intent(in) :: mesh
type(case_params_t), intent(in) :: params
real(kind=rk), intent(in) :: temperature
real(kind=rk), intent(in), optional :: Y_point(:)

Return Value real(kind=rk)


Calls

proc~~boundary_enthalpy_from_temperature~~CallsGraph proc~boundary_enthalpy_from_temperature mod_energy::boundary_enthalpy_from_temperature proc~enthalpy_from_temperature_cantera_point mod_energy::enthalpy_from_temperature_cantera_point proc~boundary_enthalpy_from_temperature->proc~enthalpy_from_temperature_cantera_point proc~enthalpy_from_temperature_value mod_energy::enthalpy_from_temperature_value proc~boundary_enthalpy_from_temperature->proc~enthalpy_from_temperature_value interface~cantera_update_thermo_c mod_energy::cantera_update_thermo_c proc~enthalpy_from_temperature_cantera_point->interface~cantera_update_thermo_c proc~build_c_species_names mod_energy::build_c_species_names proc~enthalpy_from_temperature_cantera_point->proc~build_c_species_names proc~fatal_error mod_kinds::fatal_error proc~enthalpy_from_temperature_cantera_point->proc~fatal_error

Called by

proc~~boundary_enthalpy_from_temperature~~CalledByGraph proc~boundary_enthalpy_from_temperature mod_energy::boundary_enthalpy_from_temperature proc~advance_energy_transport mod_energy::advance_energy_transport proc~advance_energy_transport->proc~boundary_enthalpy_from_temperature program~lowmach_react_hex lowmach_react_hex program~lowmach_react_hex->proc~advance_energy_transport

Source Code

   function boundary_enthalpy_from_temperature(mesh, params, temperature, Y_point) result(h_value)
      type(mesh_t), intent(in) :: mesh
      type(case_params_t), intent(in) :: params
      real(rk), intent(in) :: temperature
      real(rk), intent(in), optional :: Y_point(:)
      real(rk) :: h_value

      if (.not. params%enable_cantera_thermo) then
         h_value = enthalpy_from_temperature_value(params, temperature)
         return
      end if

      call enthalpy_from_temperature_cantera_point(params, temperature, h_value, Y_point)

      associate(dummy => mesh%ncells)
      end associate
   end function boundary_enthalpy_from_temperature