This manuscript presents an updated open-source version of the Hypersonics Task-based Research (HTR) solver. The solver, whose main features are presented in Di Renzo et al. (2020) [9] and Di Renzo & Pirozzoli (2021) [10], is designed for direct numerical simulation of reacting flows at high Reynolds numbers. This new version extends the applications of the HTR solver to turbulent combustion in the presence of external electric fields. In particular, a new distributed Poisson solver compatible with heterogeneous architectures has been incorporated in the algorithm to compute the electric potential distribution in bi-periodic configurations. The drift fluxes of the electrically charged species are now included in the transport equations using a targeted essentially non-oscillatory scheme. A verification of these new features of the solver is provided using one-dimensional burner stabilized flames, whereas a three dimensional turbulent flame is utilized to discuss the scalability of the proposed numerical tool. Program summaryProgram Title: Hypersonics Task-based Research solverCPC Library link to program files:https://doi.org/10.17632/9zsxjtzfr7.3Developer's repository link:https://github.com/stanfordhpccenter/HTR-solver.gitLicensing provisions: BSD 2-clauseProgramming language: Regent, C++, and CUDAJournal reference of previous versions:•M. Di Renzo, L. Fu, J. Urzay, HTR solver: an open-source exascale-oriented task-based multi-GPU high-order code for hypersonic aerothermodynamics, Comput. Phys. Commun. 255 (2020) 107262.•M. Di Renzo, S. Pirozzoli, HTR-1.2 solver: hypersonic task-based research solver version 1.2, Comput. Phys. Commun. 261 (2021) 107733.Does the new version supersede the previous version?: YesReasons for the new version: New features of the solverSummary of revisions:•normalization of the multicomponent mixture•Lu & Law [1] 30 species chemical scheme added to the solver as a mixture model for combustion of methane in air•TENO-LAD [2] implemented in the solver as an option for the inviscid fluxes discretization•new FFT-based Poisson solver for bi-periodic problems•calculation of transport properties using the (n,6,4) interaction theory for charge-neutral species interactions•new optional module for ion-wind effects on charged species transport•FFCM-1 [3] provided as a mixture model•Boivin et al. [4] added as a mixture model for hydrogen combustionNature of problem: This code solves the compressible multicomponent Navier-Stokes equations at high Mach numbers including finite-rate chemistry and complex chemical species transport. The implemented numerical methods are designed for direct numerical simulations (DNS) of high Reynolds number reacting flows, such as transitional and turbulent hypersonic boundary layers at high enthalpy and turbulent flames.Solution method: This code uses low-dissipation finite difference schemes for the spatial discretization of the conservation equations on Cartesian stretched grids. The time advancement is performed either with an explicit method, when the chemistry is slow and therefore does not introduce additional stiffness in the integration, or with an operator-splitting method that integrates the chemical production rates with an implicit discretization.Additional comments including restrictions and unusual features: The HTR solver builds on the runtime Legion [5] and is written in the programming language Regent [6] recently developed at Stanford University. Instructions for the installation of the components are provided in the README.md file enclosed with the HTR solver repository.