SEARCH · Engineering Papers

Results for “species distribution modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1996

An Electrothermodynamic Model with Distributed Properties for Effective Permittivities of Sea Ice

This paper presents a model to calculate the temperature dependence of effective permittivities for sea ice, a heterogeneous medium containing multiphase scatterers. With the strong permittivity fluctuation approach the model accounts for the electrodynamic scattering effect together with the quasi-static characteristics of multiple species and subspecies of inhomogeneities with distributed orientations, sizes, and shapes. Because of a preferential direction in the orientation distribution, the medium is effectively anisotropic. The size distribution is described with a probability density function in terms of normalized volumetric sizes. Scatterer shapes are nonuniform and have a general ellipsoidal form characterized by arbitrary axial ratios of correlation lengths which are related to physical geometries of the scatterers. In this formulation, sea ice consisting of solid ice, liquid brine, and gaseous inclusions is modeled to derive effective permittivities with thermodynamic phase redistribution and structural metamorphism. Theoretical results are in good agreement with experimental data at the C band frequency of 4.8 GHz for saline ice undergoing warming and cooling cycles. A competitive effect between the increase of liquid brine and the shape rounding of ellipsoidal scatterers at increasing temperatures explains the trend observed in measured data. Sensitivities of effective permittivities to structural and physical parameters characterizing sea ice are also studied.

Nghiem, S. V.↗

NASA NTRS1976

Defining constants, equations, and abbreviated tables of the 1975 US Standard Atmosphere

The U.S. Standard Atmosphere, 1975 (COESA, 1975) is an idealized, steady-state representation of the earth's atmosphere from the surface of the earth to 1000-km altitude, as it is assumed to exist in a period of moderate solar activity. From 0 to 86 km, the atmospheric model is specified in terms of the hydrostatic equilibrium of a perfect gas, with that portion of the model from 0 to 51 geopotential kilometers being identical with that of the U.S. Standard Atmosphere, 1962 (COESA, 1962). Between 51 and 86 km, the defining temperature-height profile has been modified from that of the 1962 Standard to lower temperatures between 51 and 69.33 km, and to greater values between 69.33 and 86 km. Above 86 km, the model is defined in terms of quasi-dynamic considerations involving the vertical component of the flux of molecules of individual gas species. These conditions lead to the generation of independent number-density distributions of the major species, N2, O2, O, Ar, Ne, and H, consistent with observations. The detailed definitions of the model are presented along with graphs and abbreviated tables of the atmospheric properties of the 1975 Standard.

Minzner, R. A.↗

DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1986

The removal of soluble species by warm stratiform clouds

The development of a one-dimensional, time-dependent model to study the removal of soluble gases from a warm, precipitating stratiform cloud is reported. The model calculates the distributions of water vapor and condensed water, in the form of cloud drops and raindrops, as well as the in-cloud concentration of a soluble species in the gas and aqueous phases for a specified profile of pressure and temperature and an assumed updraft velocity. Highly soluble gases are found to be rapidly dissolved into cloud droplets and then slowly incorporated into raindrops as cloudwater is converted to rainwater. The rainout rate for highly soluble species is found to be ultimately limited by the rate at which new gaseous material can be transferred from the cloud-free air into the cloud. The model calculations indicate that turbulence represents an important mechanism by which highly soluble gases are transported into stratiform clouds.

Qin, YU↗

NASA NTRS1990

Analysis of observations of the middle atmosphere from satellites

Satellite data are being used to investigate problems in middle atmosphere chemistry and dynamics. Efforts have been focused primarily on studies to determine the quality of observed distributions of trace species and derived dynamical quantities. Those data have been used as diagnostics for model-derived constituent profiles and fields and for improving our understanding of some of the fundamental processes occurring in the middle atmosphere. Temperatures and derived winds from Nimbus 7 Limb Infrared Monitoring of the Stratosphere (LIMS) data were compared with long-time series of rawinsonde data at Invercargill, New Zealand, and Berlin, West Germany, and the results are excellent for both quantities. It was also demonstrated that more highly-derived dynamical quantities can be obtained reliably from those LIMS fields. Furthermore, both the diabatic and residual-mean circulations derived using LIMS data agree qualitatively with changes in the distribution of trace species determined independently with the Nimbus 7 SAMS and LIMS experiments. Subsequently, an examination of LIMS data at mid to high latitudes of the Southern Hemisphere has revealed a synoptic-scale, upper stratospheric instability during late autumn that is associated with the development of the stratospheric polar jet. Investigation of this phenomenon continues with Stratospheric Sounding Unit (SSU) data sets.

Remsberg, Ellis E.↗

DOE OSTI2026

Mechanistic Modeling of TEG Dehydrator Emissions in Oil and Gas Industry

This work presents a mechanistic modeling approach for simulating methane emissions from triethylene glycol (TEG) dehydrators used in oil & gas (O&G) operations. The model was developed as a modular component of the Mechanistic Air Emissions Simulator (MAES) tool, incorporating species-specific absorption and emission dynamics through two-level, second-order polynomial regression (PR) models trained on ProMax simulation data: (1) species-level regression models that track the transfer rates of individual gas species within the dehydrator unit streams, and (2) outlet flow stream regression models that predict the fraction of inlet gas distributed among the outlet streams of the dehydrator unit. These behaviors were characterized over a range of glycol circulation ratios, wet gas pressures, and temperatures. The model was validated using root mean square error (RMSE) analysis. The species-level PR achieved low root mean square error (RMSE) values (<0.03) for light hydrocarbon species across all dehydrator components, ranging from 0.0009 for methane to 0.029 for normal pentane. Similarly, the outlet-level PR yielded RMSE values below 0.002 for the dry gas fraction, 0.001 for the flash tank fraction, and 0.002 for the still vent fraction, demonstrating strong agreement between predicted and reference ProMax values. When deployed at field facilities, the model significantly improved MAES-simulated dehydrator emissions, revealing that gas-assisted glycol pump emissions are the dominant contributors to both dehydrator-level and site-level methane emissions under uncontrolled conditions. Further analysis of the 154 dehydrator units reported by operators under the AMI 2024 project showed that 54 units (31%) used gas-driven glycol pumps, of which 6 units (11%) operated with uncontrolled flash tanks, and 22 units (40.7%) were identified as potentially oversized. Of the six dehydrator units with uncontrolled gas-assisted pumps, pump emissions accounted for 90.25% of total dehydrator emissions and 63.10% of total site-level emissions. These findings highlight substantial opportunities for emissions mitigation through equipment upgrades.

MAES↗

DOE OSTI2023

Imaging plant metabolism in situ

Abstract Mass spectrometry imaging (MSI) has emerged as an invaluable analytical technique for investigating the spatial distribution of molecules within biological systems. In the realm of plant science, MSI is increasingly employed to explore metabolic processes across a wide array of plant tissues, including those in leaves, fruits, stems, roots, and seeds, spanning various plant systems such as model species, staple and energy crops, and medicinal plants. By generating spatial maps of metabolites, MSI has elucidated the distribution patterns of diverse metabolites and phytochemicals, encompassing lipids, carbohydrates, amino acids, organic acids, phenolics, terpenes, alkaloids, vitamins, pigments, and others, thereby providing insights into their metabolic pathways and functional roles. In this review, we present recent MSI studies that demonstrate the advances made in visualizing the plant spatial metabolome. Moreover, we emphasize the technical progress that enhances the identification and interpretation of spatial metabolite maps. Within a mere decade since the inception of plant MSI studies, this robust technology is poised to continue as a vital tool for tackling complex challenges in plant metabolism.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗

NASA NTRS1994

The tropical water and energy cycles in a cumulus ensemble model. Part 1: Equilibrium climate

A cumulus ensemble model is used to study the tropical water and energy cycles and their role in the climate system. The model includes cloud dynamics, radiative processes, and microphysics that incorporate all important production and conversion processes among water vapor and five species of hydrometeors. Radiative transfer in clouds is parameterized based on cloud contents and size distributions of each bulk hydrometeor. Several model integrations have been carried out under a variety of imposed boundary and large-scale conditions. In Part 1 of this paper, the primary focus is on the water and heat budgets of the control experiment, which is designed to simulate the convective - radiative equilibrium response of the model to an imposed vertical velocity and a fixed sea surface temperature at 28 C. The simulated atmosphere is conditionally unstable below the freezing level and close to neutral above the freezing level. The equilibrium water budget shows that the total moisture source, M(sub s), which is contributed by surface evaporation (0.24 M(sub s)) and the large-scale advection (0.76 M(sub s)), all converts to mean surface precipitation bar-P(sub s). Most of M(sub s) is transported verticaly in convective regions where much of the condensate is generated and falls to surface (0.68 bar-P(sub s)). The remaining condensate detrains at a rate of 0.48 bar-P(sub s) and constitutes 65% of the source for stratiform clouds above the melting level. The upper-level stratiform cloud dissipates into clear environment at a rate of 0.14 bar-P(sub s), which is a significant moisture source comparable to the detrained water vapor (0.15 bar-P(sub s)) to the upper troposphere from convective clouds. In the lower troposphere, stratiform clouds evaporate at a rate of 0.41 bar-P(sub s), which is a more dominant moisture source than surface evaporation (0.22 bar-P(sub s)). The precipitation falling to the surface in the stratiform region is about 0.32 bar-P(sub s). The associated latent heating in the water cycle is the dominant source in the heat budget that generates a net upward motion in convective regions, upper stratiform regions (above the freezing level), and a downward motion in the lower stratiform regions. The budgets reveal a cycle of water and energy resulted from radiation-dynamic-convection interactions that maintain equilibrium of the atmosphere.

Sui, C. H.↗

DOE OSTI2021

Development of a reduced model for energetic particle transport by sawteeth in tokamaks

We report that the sawtooth instability is known for inducing transport and loss of energetic particles (EPs), and for generating seed magnetic islands that can trigger tearing modes. Both effects degrade the overall plasma performance. Several theories and numerical models have been previously developed to quantify the expected EP transport caused by sawteeth, with various degrees of sophistication to differentiate the response of EPs at different energies and on different orbits (e.g. passing vs. trapped), although the analysis is frequently limited to a single time slice during a tokamak discharge. This work describes the development and initial benchmark of a framework that enables a reduced model for EP transport by sawteeth retaining the full EP phase-space information. The model, implemented in the ORBIT hamiltonian particle-following code, can be used either as a standalone post-processor taking input data from codes such as TRANSP, or as a pre-processor to compute transport coefficients that can be fed back to TRANSP for time-dependent simulations including the effects of sawteeth on EPs. The advantage of the latter approach is that the evolution of the EP distribution can be simulated quantitatively for sawtoothing discharges, thus enabling a more accurate modeling of sources, sinks and overall transport properties of EP and thermal plasma species for comprehensive physics studies that require detailed information of the fast-ion distribution function and its evolution over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1998

Spatial Aspects of Interspecific Competition

Using several variants of a stochastic spatial model introduced by Silvertown et al., we investigate the effect of spatial distribution of individuals on the outcome of competition. First, we prove rigorously that if one species has a competitive advantage over each of the others, then eventually it takes over all the sites in the system. Second, we examine tradeoffs between competition and dispersal distance in a two-species system. Third, we consider a cyclic competitive relationship between three types. In this case, a nonspatial treatment leads to densities that follow neutrally stable cycles or even unstable spiral solutions, while a spatial model yields a stationary distribution with an interesting spatial structure.

Durrett, Rick↗

NASA NTRS1999

Trace Chemistry

The goals of the trace chemistry group were to identify the processes relevant to aerosol and aerosol precursor formation occurring within aircraft gas turbine engines; that is, within the combustor, turbine, and nozzle. The topics of discussion focused on whether the chemistry of aerosol formation is homogeneous or heterogeneous; what species are important for aerosol and aerosol precursor formation; what modeling/theoretical activities to pursue; what experiments to carry out that both support modeling activities and elucidate fundamental processes; and the role of particulates in aerosol and aerosol precursor formation. The consensus of the group was that attention should be focused on SO2, SO3, and aerosols. Of immediate concern is the measurement of the concentration of the species SO3, SO2, H2SO4 OH, HO2, H2O2, O, NO, NO2, HONO, HNO3, CO, and CO2 and particulates in various engines, both those currently in use and those in development. The recommendation was that concentration measurements should be made at both the combustor exit and the engine exit. At each location the above species were classified into one of four categories of decreasing importance, Priority I through IV, as follows: Combustor exit: Priority I species - SO3:SO2 ratio, SO3, SO2, and particulates; Priority II species: OH and O; Priority III species - NO and NO2; and Priority IV species - CO and CO2. For the Engine exit: Priority I species - SO3:SO2 ratio, SO3, SO2,H2SO4, and particulates; Priority II species: OH,HO2, H2O2, and O; Priority III species - NO, NO2, HONO, and HNO3; and Priority IV species - CO and CO2. Table I summarizes the anticipated concentration range of each of these species. For particulate matter, the quantities of interest are the number density, size distribution, and composition. In order to provide data for validating multidimensional reacting flow models, it would be desirable to make 2-D, time-resolved measurements of the concentrations of the above species and, in addition, of the pressure, temperature, and velocity. A near term goal of the experimental program should be to confirm the nonlinear effects of sulfur speciation, and if present, to provide an explanation for them. It is also desirable to examine if the particulate matter retains any sulfur. The recommendation is to examine the effects on SOx production of variations in fuel-bound sulfur and aromatic content (which may affect the amount of particulates formed). These experiments should help us to understand if there is a coupling between particulate formation and SO, concentration. Similarly, any coupling with NOx can be examined either by introducing NOx into the combustion air or by using fuel-bound nitrogen. Also of immediate urgency is the need to establish and validate a detailed mechanism for sulfur oxidation/aerosol formation, whose chemistry is concluded to be homogeneous, because there is not enough surface area for heterogeneous effects. It is envisaged that this work will involve both experimental and theoretical programs. The experimental work will require, in addition to the measurements described above, fundamental studies in devices such as flow reactors and shock tubes. Complementing this effort should be modeling and theoretical activities. One impediment to the successful modeling of sulfur oxidation is the lack of reliable data for thermodynamic and transport properties for several species, such as aqueous nitric acid, sulfur oxides, and sulfuric acid. Quantum mechanical calculations are recommended as a convenient means of deriving values for these properties. Such calculations would also help establish rate constants for several important reactions for which experimental measurements are inherently fraught with uncertainty. Efforts to implement sufficiently detailed chemistry into computational fluid dynamic codes should be continued. Zero- and one-dimensional flow models are also useful vehicles for elucidating the minimal set of species and reactions that must be included in two- and three-dimensional modeling studies.

Radhakrishnan, Krishnan↗

DOE OSTI2025

Atom Identification in Bilayer Moiré Materials with Gomb-Net

Moiré patterns in van der Waals bilayer materials complicate the analysis of atomic-resolution images, hindering the atomic-scale insight typically attainable with scanning transmission electron microscopy. Here, we report a method to detect the positions and identities of atoms in each of the individual layers that compose twisted bilayer heterostructures. We developed a deep learning model, Gomb-Net, which identifies the coordinates and atomic species in each layer, deconvoluting the moiré pattern. This enables layer-specific mapping of atomic positions and dopant distributions, unlike other commonly used segmentation models which struggle with moiré-induced complexity. Using this approach, we explored the Se atom substitutional site distribution in a twisted fractional Janus WS 2 -WS 2(1–x) Se 2x heterostructure and found that layer-specific implantation sites are unaffected by the moiré pattern’s local energetic or electronic modulation. In conclusion, this advancement enables atom identification within material regimes where it was not possible before, opening new insights into previously inaccessible material physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

National forest timber bids and export price interlinkages in the USA: The bounds testing approach

We examine the interrelationships of national forest timber bid and log export prices for Douglas-fir (Region 6) in 2003–2021 and loblolly pine (Region 8) using time-series stationarity and bounds-testing approaches. The analysis shows the difference in time series properties across different geographies and time scales. Bid prices for these dominant species of timber in each of two Forest Service regions are found to be stationary, while log export prices for those same species appear to be non-stationary. Bounds testing conducted with autoregressive distributed lag modeling methods provides strong evidence for a cointegrating relationship between timber bid prices and log export prices in Forest Service Region 6. In Region 8, evidence for a cointegrating relationship was weak. Here, the share of salvage wood in national forest timber sales is a statistically significant control variable in all models estimated.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗

NASA NTRS1987

Middle atmosphere trace constituents; Proceedings of Workshop 9 of the Twenty-sixth COSPAR Plenary Meeting, Toulouse, France, June 30-July 11, 1986

Papers are presented on a reference model for stratospheric NO2, a proposed reference model for nitric acid, proposed reference models for ozone, and a background stratospheric aerosol reference model. Other topics include a comparison between the observed and calculated distributions of trace species in the middle atmosphere, recent measurements of trace constituents from rocket and balloon probes, NO2 field observations in the stratosphere during the MAP/GLOBUS Campaign in 1983 and 1985, and the vertical distribution of ozone in the equatorial zone. Also considered are precision measurements of stratospheric ozone profiles by rocket-borne optical ozonesondes, vertical profiles of source gases at tropical latitudes, and measurements of stratospheric trace constituent distributions from balloon-borne far infrared observations.

Keating, G. M.↗

NASA NTRS1974

Reacting viscous-shock-layer solutions with multicomponent diffusion and mass injection

Numerical solutions of the viscous-shock-layer equation where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium are presented. Also the effects of the diffusion model, surface catalysis, and mass injection on surface transport and flow parameters are considered. The flow is treated as a mixture of five inert and thermally perfect species. The viscous-shock-layer equations are solved by using an implicit-difference scheme. All calculations are for hyperboloids with included angles of 20 and 45. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented to show the effects of the chemical models; diffusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; shock standoff distance; wall pressure distribution; and tangential velocity, temperature, and species profiles. The results show that an equilibrium analysis can substantially overpredict the heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. Moreover, at such conditions surface catalysis significantly influences heat-transfer and flow-field properties. If a binary rather than a multicomponent diffusion model is assumed, negligible errors in most flow properties result. Quantitative results are presented that show the effect of mass injection on flow properties within and downstream of the injection region.

Moss, J. N.↗

DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗