recover_temperature_and_update_thermo_cantera Subroutine

private subroutine recover_temperature_and_update_thermo_cantera(mesh, params, energy, species_Y)

Recover T from h and refresh cp/lambda/rho_thermo in one Cantera sync.

This is the preferred energy-step thermo sync routine: (T, cp, lambda, rho_thermo) = sync(h, Y, p0). It preserves the transported sensible enthalpy field and updates only derived thermo state: T, cp, lambda, and diagnostic rho_thermo.

Arguments

Type IntentOptional Attributes Name
type(mesh_t), intent(in) :: mesh
type(case_params_t), intent(in) :: params
type(energy_fields_t), intent(inout) :: energy
real(kind=rk), intent(in), optional :: species_Y(:,:)

Calls

proc~~recover_temperature_and_update_thermo_cantera~~CallsGraph proc~recover_temperature_and_update_thermo_cantera mod_energy::recover_temperature_and_update_thermo_cantera interface~cantera_recover_temperature_and_update_thermo_c mod_energy::cantera_recover_temperature_and_update_thermo_c proc~recover_temperature_and_update_thermo_cantera->interface~cantera_recover_temperature_and_update_thermo_c proc~build_c_species_names mod_energy::build_c_species_names proc~recover_temperature_and_update_thermo_cantera->proc~build_c_species_names proc~build_thermo_y mod_energy::build_thermo_Y proc~recover_temperature_and_update_thermo_cantera->proc~build_thermo_y proc~fatal_error mod_kinds::fatal_error proc~recover_temperature_and_update_thermo_cantera->proc~fatal_error proc~build_thermo_y->proc~fatal_error

Called by

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

Source Code

   subroutine recover_temperature_and_update_thermo_cantera(mesh, params, energy, species_Y)
      type(mesh_t), intent(in) :: mesh
      type(case_params_t), intent(in) :: params
      type(energy_fields_t), intent(inout) :: energy
      real(rk), intent(in), optional :: species_Y(:,:)

      integer :: n_len
      real(rk), allocatable :: P_arr(:), Y_local(:,:)
      character(kind=c_char), allocatable :: c_names_flat(:)

      if (.not. params%enable_cantera_thermo) return

      if (params%nspecies <= 0) then
         call fatal_error('energy', 'enable_cantera_thermo requires at least one inert thermo species')
      end if

      allocate(P_arr(mesh%ncells))
      P_arr = params%background_press
      call build_thermo_Y(params, mesh%ncells, Y_local, species_Y)
      call build_c_species_names(params, c_names_flat, n_len)

      call cantera_recover_temperature_and_update_thermo_c(mesh%ncells, &
                                   energy%h, &
                                   P_arr, &
                                   params%nspecies, &
                                   Y_local, &
                                   energy%T, &
                                   energy%cp, &
                                   energy%lambda, &
                                   energy%rho_thermo, &
                                   params%energy_reference_T, &
                                   c_names_flat, &
                                   n_len)

      deallocate(P_arr)
      deallocate(Y_local)
      deallocate(c_names_flat)
   end subroutine recover_temperature_and_update_thermo_cantera