Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemically reacting flows”

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.

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

608 records · Page 34

The sensitivity of tropospheric chemistry to cloud interactions

Clouds, although only occupying a relatively small fraction of the troposphere volume, can have a substantial impact on the chemistry of the troposphere. In newly formed clouds, or in clouds with air rapidly flowing through, the chemistry is expected to be far more active than in aged clouds with stagnant air. Thus, frequent cycling of air through shortlived clouds, i.e. cumulus clouds, is likely to be a much more efficient media for altering the composition of the atmosphere than an extensive cloud cover i.e. frontal cloud systems. The impact of clouds is tested out in a 2-D channel model encircling the globe in a latitudinal belt from 30 to 60 deg N. The model contains a detailed gas phase chemistry. In addition physiochemical interactions between the gas and aqueous phases are included. For species as H2O2, CH2O, O3, and SO2, Henry's law equilibria are assumed, whereas HNO3 and H2SO4 are regarded as completed dissolved in the aqueous phase. Absorption of HO2 and OH is assumed to be mass-transport limited. The chemistry of the aqueous phase is characterized by rapid cycling of odd hydrogen, (H2O2, HO2, and OH). O2(-) (produced through dissociation of HO2) reacting with dissolved O3 is a major source of OH in the aqueous phase. This reaction can be a significant sink for O3 in the troposphere. In the interstitial cloud air, odd hydrogen is depleted, whereas NO(x) remains in the gas phase, thus reducing ozone production due to the reaction between NO and HO2. Our calculations give markedly lower ozone levels when cloud interactions are included. This may in part explain the overpredictions of ozone levels often experienced in models neglecting cloud chemical interactions. In the present study, the existence of clouds, cloud types, and their lifetimes are modeled as pseudo random variables. Such pseudo random sequences are in reality deterministic and may, given the same starting values, be reproduced. The effects of cloud interactions on the overall chemistry of the troposphere are discussed. In particular, tests are performed to determine the sensitivity of cloud frequencies and cloud types.

Jonson, Jan E.↗

An Adaptive-Mesh-Refinement Based Computational Tool for Simulating Catalysis at Mesoscale

In this work, we present a computational tool for mesoscale applications using open-source exascale- computing compatible adaptive-mesh-refinement (AMR) library, AMReX [2]. AMReX is software library that enables development of application solvers with block-structured Cartesian AMR. Our tool has capabilities to include realistic geometry representation, chemical species transport, reactions and thermodynamics that are critical for capturing mesoscale physics. A significant achievement is the ability of our solver to automatically import electron microscopy data in the form of a stereolithography (STL) or pixelated file format (mrc, tiff) without undergoing the tedious task of unstructured mesh generation. This feature allows for rapid simulation of catalyst particles with complex morphologies using an immersed-boundary formulation. The use of AMR allows for higher resolutions at catalyst surface interfaces, which in turn provides an accurate description of surface reactions and transport. Our solver uses a hybrid distributed and shared memory parallelism (OpenMP/GPU-based) with which strong scaling up to 10,000 processors for realistic catalyst particle simulations have been demonstrated.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Corrosion Prevention in Copper Combustion Chamber Liners of Liquid Oxygen/Methane Booster Engines

Prior research conducted by United Technologies Research Center and Rockwell International Rocketdyne Division encountered severe copper corrosion while flowing hydrocarbon fuels through copper cooling tubes. These experimental results have very important implications for the development of low-cost, regeneratively cooled, oxygen/hydrocarbon bipropellant booster engines. Aerojet TechSystems undertook a program sponsored by NASA/Lewis Research Center in 1986 with two objectives: (1) to define the corrosive interaction process that occurs between hydrocarbon fuels and copper combustion chamber liner materials, using both static and dynamic test methods, and (2) to identify and demonstrate protective measures against this corrosive process. The results of the first program task demonstrated that the most damaging corrosive process was caused by trace amounts of sulfur-containing impurities in the hydrocarbon fuel which react with the copper chamber liner material to form cuprous sulfide (Cu2S), as reported in Paper No. AlAA 88-3215. Experiments were conducted which demonstrated this corrosive process with Mil-Spec RP-1, propane, and methane fuels in contact with OFHC, NASA-Z, and ZrCu coppers. The preliminary results of the second program task demonstrated that electrodeposited gold and platinum coatings held promise of preventing corrosion of copper chamber liners by methane contaminated with sulfur-containing impurities, as reported in Paper No. AIAA 89-2738. Experiments were conducted which demonstrated greatly reduced corrosion of copper by hydrogen sulfide (H2S) and methyl mercaptan (CH3SH) as measured impurities in methane. This paper presents the concluding, confirmatory results from the second program task. The efficacy of a protective gold coating on NASA-Z copper cooling channel walls was conclusively demonstrated. Tests were conducted in the Aerojet Carbothermal Test Facility which provide realistic simulations of booster engine cooling channel conditions such as temperature, pressure, flow velocity, and heat flux. Dynamic test series were performed with methane containing 5 and 10 ppm (by volume) hydrogen sulfide. Gold-plated and unplated copper alloy specimens were evaluated. The tests demonstrated that gold coatings were effective in preventing corrosion of the copper. Posttest metallographic examinations of the unplated specimens showed severe corrosion as a result of reaction with the sulfur containing contaminant in the fuel. In contrast, posttest metallographic examination of the gold-coated specimens showed no corrosion under similar operating conditions. The findings from metallurgical and chemical analyses of the copper specimens, chemical analysis of the methane, and description of all test conditions are provided. Conclusions are reached as to the efficacy of the plating methods and the effectiveness of the metal coatings in the prevention of corrosion in hydrocarbon-fueled booster engine combustion chamber liners.

S D Rosenberg↗

Thermodynamic constraints on the textural evolution of eucrite, EET 90020

Basaltic eucrites, which are possibly parts of the Vestan crust, formed as either lava flows or intrusions, and can offer insights into early crust formation. However most eucrites have experienced thermal metamorphism, which exacerbates the challenges of understanding these samples. A variety of petrogenetic models such as, partial melting of a primitive source, fractional crystallization of magmas emplaced in a primitive crust, and/or partial melting of a eucritic source coupled with melt mixing and assimilation have been proposed in order to explain observed petrologic and chemical characteristics of eucrites. EET 90020 is an unbrecciated eucrite that has experienced significant thermal metamorphism, with temperatures of metamorphic equilibration calculated to be ~840-1040°C. Its petrologic history remains contentious, in part due to differences in trace element analyses [10,17, this work], which has led to multiple metamorphic interpretations [8,10-11]. Here, we combined petrologic observations, chemical analyses, and thermodynamic modeling, to interpret micro-domain textures identified in EET 90020 and better constrain its petrologic history. Ultimately, the development of such textures is a direct result of the geologic processes operating on a young Vesta or similar asteroid. Sample description: We identified three textural domains in EET 90020 (Figure 1); (I) A coarse grain domain with granoblastic plagioclase and pyroxene, (II) a medium grain domain with curved grain boundaries between plagioclase and pyroxene, and (III) a fine grain domain dominated by tridymite that fills in interstitial space between spherical plagioclase. Major element chemistry in phases is consistent across domains. Pyroxene have distinctive high Ca lamellae (Wo40.92En28.1Fs31.0) forming from the low-Ca pigeonite host (Wo3.5En33.7Fs62.8). There are small amounts of fayalitic olivine along pyroxene/oxide grain boundaries, where pyroxenes have reacted with ilmenite/chromite clasts. Plagioclase is anorthitic (Avg ~An88), with little variation (An86-92.5). Trace element bulk rock analyses were collected via ICP-MS and individual mineral analyses with LA-ICP-MS (Figure 2). The bulk sample has an observable enrichment in light rare earth elements, with a notable depletion in Eu. Plagioclase is enriched in LREEs, depleted in HREEs, and has a positive Eu anomaly while pyroxene is enriched in HREEs with a negative Eu anomaly. Variations in trace element analyses is likely influenced by the presence of phosphates because are highly enriched in REEs and have a dramatic effect on the bulk domain compositions. Thermodynamic modeling: The software package, Perple_X, was used to calculate mineral phase equilibria over a range of conditions using a Gibbs free energy minimization approach [12]. Models used all major/minor elements except P2O5. Thermodynamic properties from [13] were used to constrain endmember phase stabilities. Activity models were used to describe mixing in phases with solid solution. Isochemical P-T phase diagrams were constructed for the bulk thin section, and the coarse, medium and fine grain domain compositions. Domain compositions were determined using microprobe analyses of phases identified in each domain and observed modal mineral abundancies. Thermodynamically stable pyroxene compositions were extracted from model results and compared to pyroxene compositions collected via EMPA, allowing us to calculate temperatures for metamorphic equilibrium. Results: Results that use the coarse domain composition (Fig. 3a) demonstrate overlap at T ~ 1020°C indicating high T equilibration. There was no overlap for the medium and fine grain domains, or the bulk composition (e.g., Fig. 3b) indicating that these are not equilibrium assemblages. These results were replicated using data from [8,20]. Discussion: The disequilibrium results from the bulk rock, fine grain and medium grain models imply that EET 90020 was modified during and/or after peak metamorphism (open system behavior). However, this is in contrast to the coarse domain where it appears as if equilibrium was maintained. We suggest that this conflict is explained by highly localized equilibrium occurring at millimeter length-scales, which has been observed in terrestrial metamorphism where aqueous fluid infiltrate [14] or melt loss occurs [15]. It is unlikely that aqueous alteration was significant in modifying the bulk composition of EET 90020, because it lacks hydrated minerals. However, the presence of melt is supported by the presence of an abundance of curved grain boundaries and spherical grains in the medium and fine grain domains This contrasts the coarse domain, where grain boundaries are angular. Given the evidence for disequilibrium in the medium and fine grain domains, it would be inappropriate to draw additional conclusions about textural evolution using thermodynamic models of those domains. However, modeling results of the coarse grain domain, where equilibrium was maintained, can provide useful insights into the geologic evolution of EET 90020. We approximate maximum metamorphic temperatures, using the Ca component in pyroxene, at ~1020°C (Fig. 4), which is consistent with previous two-pyroxene thermometry [11]. At this temperature, up to 10 vol % melt can be produced in the coarse grain domain, which is at the boundary of minimum melt needed in order to segregate and form a melt network [16]. Tridymite and plagioclase would be the first minerals to melt out during heating starting ~990°C. This is consistent with the observed melt textures formed by tridymite and plagioclase in the fine grain domain. We suggest that a partial melting model best explains the textural development of EET 90020. Upon heating, melt generated was either trapped, resulting in no net change in the localized rock composition, or migrated, resulting in disequilibrium between pyroxene and the surrounding matrix. A positive Eu anomaly in plagioclase indicates that all of the plagioclase was melted previously. This is consistent with previous studies that have suggested that EET 90020 represents a residual eucrite that experienced partial melting and subsequent melt loss [8]. We speculate that partial melting could have occurred when EET 90020 was heated in the asteroid’s lower crust or adjacent to an intrusion [7]. Future work will study other, texturally heterogeneous samples to see if EET 90020 presents a unique circumstance or whether crustal partial melting was ubiquitous on the eucrite asteroid.

Eucrites↗

Finite-rate entrainment effects on nitrogen oxide (NO x ) emissions in staged combustors

Axially staged combustors offer the potential for achieving low NO x emissions from gas turbines at the elevated temperatures required to achieve at least 65% combined cycle efficiencies. While there is evidence that premixing between the main burner products and secondary stream reduces NO x , it is important to characterize the specific requirements needed to achieve low emissions at gas turbine conditions. This study examines the sensitivity of NO x to finite-rate, large-scale entrainment of the main and secondary streams under the simplifying assumption of infinitely-fast small-scale mixing. Essentially, this allows us to isolate the effects of limited large-scale entrainment rates by using a homogeneous/uniform condition for the entrained and reacting gases. We use a reduced-order reactor network model to examine a generic staged-combustor whose inputs are physical time scales that embody the finite-rate entrainment characteristics. With simulations conducted over a large parameter space at typical operating conditions (25 atm, 650 K air), the results show that, even when small-scale mixing is infinitely fast and the reaction zone is uniform, entrainment significantly affects NO x emissions due to its influence on the equivalence ratio — and thus the temperature and time — at which the entrained mixture ignites. Fuel-air staging is shown to be vital to NO x reduction in most practical cases where entrainment times exceed roughly 1 ms and the secondary stream finishes entraining before the main burner products. A constrained NO x minimization indicates that using the leanest possible main burner and re-routing as much air as possible to the secondary stage is key to low NO x with finite-rate entrainment; for example, less than 10 ppm NO x can be achieved if the secondary stage contains 20% of the total mass flow.

03 NATURAL GAS↗

Evaluation of thermodynamic closure models for partially reacted two-phase mixture of condensed phase explosives

One of the key fundamental issues that is crucial in the continuum modeling of reactive flow phenomena is the thermodynamically consistent description of reaction mixture properties. To define the mixture properties, thermodynamic closure rules that relate the properties of the individual reaction components to the mixture properties are required. In the context of reactive two-phase modeling approaches, various strategies to define the thermodynamic closures have been adopted such as pressure temperature (PT) equilibrium between the individual reaction components, pressure (specific) volume (PV) equilibrium, etc. The choice of closure rules determines the relative distribution of specific volume and energy across the reaction components that comprise the mixture. Therefore, depending on the choice of the closure, the mixture thermodynamic behavior can vary. The present work examines the effect of different closure approaches on the thermodynamic properties of the reaction mixture. The analysis is performed for a condensed phase HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) based plastic bonded explosive (PBX) 9501 explosive using four different thermodynamic closures, viz., PT equilibrium, PV equilibrium, volume temperature (VT) equilibrium, and pressure (P) equilibrium with reactants on an isentrope. Furthermore, the relative variations in the thermodynamic properties of the mixture are analyzed and compared under both compression and expansion loading regimes. It is shown that out of the four closure models, only PT equilibrium and P equilibrium closures lead to a thermodynamically accurate description of the mixture under both compression and expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Combustion of Metals in Reduced-Gravity and Extraterrestrial Environment

As a result of the ongoing exploration of Mars and the several unmanned and manned missions planned for the future, increased attention has been given to the use of the natural resources of the planet for rocket propellant production and energy generation. Since the atmosphere of Mars consists of approximately 95% carbon dioxide (CO2), this gas is the resource of choice to be employed for these purposes. Unfortunately, CO2 is also a final product in most combustion reactions, requiring further processing to extract useful reactants such as carbon monoxide (CO), oxygen (O2), and hydrocarbons. An exception is the use Of CO2 as an oxidizer reacting directly with metal fuel. Since many metals burn vigorously with CO2, these may be used as an energy source and as propellants for an ascent/descent vehicle in sample-collection missions on Mars. In response to NASA's Human Exploration and Development of Space (HEDS) Enterprise to search for appropriate in-situ resource utilization techniques, this investigation will study the burning characteristics of promising metal/CO2 combinations. The use of reduced gravity is essential to eliminate the intrusive buoyant flows that plague the high-temperature metal reactions, to remove the destructive effect of gravity on the shape of molten metal samples, and to study the influence of radiative heat transfer from solid oxides undisturbed by natural convection. In studies with large metal specimens, the burning process is invariably influenced by strong convective currents that accelerate the reaction and shorten the burning times. Although these currents are nearly absent from small burning particles, the high emissivity of the flames, rapid reaction, small length scales, and intermittent explosions make the gathering of any useful information on burning rates and flame structure very difficult. This investigation has the ultimate goal of providing a careful probing of flame structure and dynamics by taking advantage of large, free-floating spherical metal samples and their corresponding long burning times available in reduced gravity. The first set of experiments has been conducted with magnesium (Mg) samples burning in the low-gravity environment generated by an aircraft flying parabolic trajectories. Owing to its high adiabatic flame temperature, oxidizer/fuel ratio, and heat per unit mass of fuel, as well as its low toxicity and low ignition temperature, Mg has been identified as a promising metal fuel with CO2 as oxidizer. The experimental effort is complemented by the development of a numerical model combining gas-phase chemical kinetics and transport mechanisms.

Abbud-Madrid, A.↗

Study of Buoyancy Effects in Diffusion Flames Using Rainbow Schlieren Deflectometry

Diffusion flames are extensively encountered in many domestic and industrial processes. Even after many decades of research, a complete understanding of the diffusion flame structure is not available. The structure and properties of the flames are governed by the mixing (laminar or turbulent), chemical kinetics, radiation and soot processes. Another important phenomenon that affects flame structure in normal gravity is buoyancy. The presence of buoyancy has long hindered the rational understanding of many combustion processes. In gas jet diffusion flames, buoyancy affects the structure of the shear layer, the development of fluid instabilities, and formation of the coherent structures in the near nozzle region of the gas jets. The buoyancy driven instabilities generate vorticial structures outside the flame resulting in flame flicker. The vortices also strongly interact with the small-scale structures in the jet shear layer. This affects the transitional and turbulence characteristics of the flame. For a fundamental understanding of diffusion flames it is essential to isolate the effects of buoyancy. This is the primary goal of the experiments conducted in microgravity. Previous investigations, have shown dramatic differences between the jet flames in microgravity and normal gravity. It has been observed that flames in microgravity are taller and more sooty than in normal gravity. The fuels used in these experiments were primarily hydrocarbons. In the absence of buoyancy the soot resides near the flame region, which adversely affects the entrainment of reactants. It is very important to eliminate the interference of soot on flame characteristics in microgravity. The present work, therefore, focuses on the changes in the flame structure due to buoyancy without the added complexities of heterogeneous reactions. Clean burning hydrogen is used as the fuel to avoid soot formation and minimize radiative losses. Because of the low luminosity of hydrogen flames, we use rainbow schlieren deflectometry for visualization. The visualized images are digitized for quantification.The work reported here is divided into three sections; rainbow schlieren deflectometry (RSD), microgravity experiments and sub-atmospheric pressure experiments. The first section demonstrates the application of RSD for quantitative measurements in non-reacting and reacting flow systems. A computational effort to complement the experimental work is also included. In the second section, the experiments conducted at the 2.2s NASA Lewis Drop tower facility are described. The experiments were conducted to study the behavior of laminar, transitional and turbulent hydrogen flames in microgravity. The ability of RSD technique to provide quantitative data is highlighted. The final section deals with the sub-atmospheric pressure tests, which demonstrate that buoyancy in hydrogen diffusion flames can be scaled with pressure at normal gravity.

Agrawal, Ajay K.↗

Trade-Offs of Simplified Versus Comprehensive Representation of Mineralogy When Studying Dust Impacts on Earth’s Climate Systems

The intensity and direction of dust impacts on Earth’s climate systems depend on mineral composition. For example, the presence or absence of a few percent of iron oxides in dust will determine if dust is warming or cooling the atmosphere. Similarly, feldspar will enhance ice cloud formation, while acid gases in the atmosphere will react on the surface of dust calcite limiting acid rain. Still, most climate models use a simplified representation of dust mineralogy. They assume a fixed composition at emission which stays invariant during transport and removal. Such simplification assumes spatially and temporally constant physical and chemical properties of dust, and appears to provide satisfactory results when comparing some properties with observations. The trade-off is their lack of spatial gradients, which will fail to induce circulation, cloud and precipitation changes. The two reasons to omit mineral variations are the uncertainty of current atlases of soil mineral composition in arid regions, and, more practically, an improved runtime efficiency. The former reason is losing ground with the recent launch (July 2022) of a dedicated mission (NASA/JPL EMIT) to retrieve global soil mineralogy of dust sources at high spatial resolution. While the EMIT science team is finalizing a satisfactory global map of mineral composition of dust sources, we analyzed the interaction of dust mineralogy on radiation and its impact on the fast temperature response using different representations of mineral composition from detailed and spatially varying to simplified and globally uniform, assuming different hematite contents and methods to calculate optical properties. Our results show that resolving dust mineralogy reduces dust absorption, and results in improved agreement with observation-based single scattering albedo (SSA), radiative fluxes from CERES (the Clouds and the Earth’s Radiant Energy System), and land surface temperature from CRU (Climatic Research Unit), compared to the baseline bulk dust model version. It also results in distinct radiative impacts on Earth’s climate over North Africa. From our 19-year simulation, we will show that it leads to a reduction of over 50% in net downward radiation at top of atmosphere (TOA) across the Sahara and an approximately 20% reduction over the Sahel. We will explain how the surface temperature response affects the monsoon flow from the Gulf of Guinea. Interestingly, we find similar results by simply fixing the hematite content of dust to a globally uniform value of 0.9% by volume. We will discuss the underlying reasons for such results and show that they may be unrelated to the distribution of soil mineralogy. Still, an accurate representation of soil mineralogy is necessary to better understand dust impacts on the Earth’s climate systems.

dust impacts↗

Detonation Waves in High Explosives

A material at high temperature can react or decompose. For an energetic material, the reaction is exothermic and releases chemical energy that would further increase the temperature. Under some circumstances, when a reaction is triggered, such a reaction can propagate and the material rapidly releases a large amount of energy giving rise to an explosion. Examples of such materials are aerosols, suspensions of solid particles or liquid droplets in a gas; such as coal dust, grain dust and fuel-air explosions. Frequently, explosions are due to accidents. A spectacularly destructive example is the recent explosion of a large quantity of ammonium nitrate (thousands of tons) in Beirut, Lebanon (August 2020); see for example Beirut explosion. Ammonium nitrate is used as a fertilizer. It and the aerosols are not considered to be explosives due to the limited conditions for which an explosion can occur. An aerosol gets the oxidizer from the surrounding air. Burning requires diffusion of the oxidizer to the particle surface where the reaction occurs. A large density of small particles is required for a fast enough reaction to support an explosion. In contrast, an explosive is an energetic material with both fuel and oxidizer mixed on a molecular scale (either premixed gases or within molecules of a solid). This allows fast enough reactions over a wide range of conditions to support a self-propagating reactive wave known as a detonation wave. A detonation wave can be controlled and an explosive used for useful purposes such as in mining, construction, demolition, explosive welding, argon flash lamp, pulsed power using a magnetic flux generator [see also Goforth et al., 2015], jet cutter with shaped charge, explosive art, and generating conditions to study the response of materials at high strain rates and high pressures [see for example, Marsh, 1980]. Explosives are also used in conventional munitions and nuclear weapons. The focus of this book is on the theory and phenomenology of solid high explosives (HEs); in particular, plastic-bonded explosives (PBXs). Some aspects of detonation wave theory are needed to interpret explosive data. Hence, the theory is presented before the detonation wave phenomenology. A familiarity with fluid flow, specifically the notion of shock waves and the shock loci are assumed. In the remainder of this chapter we give a brief overview on the basic properties of detonation waves and PBXs.

36 MATERIALS SCIENCE↗

Lipopolysaccharide modulation of a CD14-like molecule on porcine alveolar macrophages

Cluster of differentiation antigen 14 (CD14) functions as a receptor for lipopolysaccharide (LPS) LPS-binding protein (LBP) complexes. Because LPS has varying effects on CD14 expression in vitro, we evaluated CD14 expression in response to LPS with a fully differentiated macrophage phenotype, the alveolar macrophage. By using flow microfluorometric analysis and a radioimmunoassay with an anti-human CD14 monoclonal antibody (My4) that cross-reacts with porcine CD14, we found that macrophages stimulated with LPS for 24 h exhibited a two- to fivefold increase in CD14-like antigen compared with unstimulated cells. At low concentrations of LPS, up-regulation of the CD14-like antigen was dependent on serum; at higher concentrations of LPS, serum was not required. In the absence of serum a 10-fold higher dose of LPS (10 ng/ml) was required to increase CD14-like expression. In addition, LPS-induced CD14-like up-regulation correlated with secretion of tumor necrosis factor-alpha, regardless of serum concentration. Blockade with My4 antibody significantly inhibited LPS-induced tumor necrosis factor-alpha secretion at 1 ng/ml of LPS. However, inhibition decreased as we increased the LPS concentration, suggesting the existence of CD14-independent pathways of macrophage activation in response to LPS. Alternatively, My4 may have a lower affinity for the porcine CD14 antigen than LPS, which may have only partially blocked the LPS-LBP binding site at high concentrations of LPS. Therefore, these data suggest that LPS activation of porcine alveolar macrophages for 24 h increased CD14-like receptor expression. The degree of CD14-like up-regulation was related to LPS concentration, however, activation did not require the presence of serum at high concentrations of LPS.

Non-NASA Center↗

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↗

Ammonia Permeation Through Inorganic Membranes

Moisture probes are used in the tritium facility to measure the moisture content of the gas stream. When tritium reacts with nitrogen to form ammonia, the ammonia gas damages the sensitive moisture probe. The gas stream tested during this experiment for moisture probe protection is composed of nitrogen and argon gas in conjunction with ammonia and water vapor. The gas stream is filtered through an inorganic membrane due to the high reactivity of tritium and a residual gas analyzer (RGA) is used to measure the ammonia and water vapor in the gas stream, as well as the permeability of nitrogen and argon through the inorganic membranes. The membranes tested were boron nitride and silver. The permeability was compared for each membrane at different pressures and flow rates. Fourier-transformed infrared spectroscopy (FT-IR) was used to compare the membranes before and after ammonia permeation. Ammonia gas is a byproduct formed from tritium gas reacting with nitrogen in the gas stream, similar to the Haber Process which occurs when hydrogen reacts with nitrogen [2]. The corrosive nature of ammonia damages metals and sensitive instruments such as moisture probes. Moisture probes are essential in various gas and liquid processes for measuring the moisture content of a gas. Membrane research has found that graphene oxide membranes successfully uptake ammonia from the gas stream at a humidity of approximately 30% at which point the pores pin closed, preventing gas permeation. Further research is being conducted to find an inorganic membrane that will prevent ammonia permeation thus circumventing damage to the sensitive moisture probe. The permeation measurements for argon and nitrogen gas and water and ammonia vapor were measured at various pressures of 800, 1000, 1200, 1500, 1800, and 2000 torr. The permeability average (PA) was calculated for each filter using a Barrer calculation. The average permeations for each filter were then compared. The gas composition was measured using the RGA Nitrogen was permeated through the membranes after ammonia permeation and a FT-IR spectrum was obtained. Boron Nitride: As the nitrogen pressure increased, the permeability of the membrane decreased from 800-1200 torr and remained relatively unchanged from 1200-2000 torr. The membrane exhibited a slight decrease in argon permeability from 800-1200 torr and a slight increase from 1500-2000 torr. The boron nitride membrane when permeated with ammonia and nitrogen prevented less ammonia permeation compared to the blank filter. The boron nitride filter demonstrated a decrease in permeability after ammonia permeation. Silver Membrane: When permeated with nitrogen alone, the membrane permeability decreased from 800 torr to 1200 torr. The permeability of the membrane showed little difference as pressure increased from 1200 torr. When permeated with argon, the same phenomenon occurred but on a much smaller scale. When compared to the blank and the boron nitride membranes the silver membrane exhibited significantly lower permeability. This may be in part due to the coating on the silver membrane and differences in pore size of the membranes. During permeation with ammonium hydroxide vapor and nitrogen, the partial pressure of the ammonia increased as compared to the blank membrane. The permeability before and after permeation with ammonia and nitrogen were compared. Permeability values suggested that after being permeated with the ammonia that the permeability for lower pressures increased slightly. There was no significant change in the FT-IR spectrums before and after permeation. Future work: Continue testing various inorganic membranes for ammonia permeation. Further research chemical separation of ammonia from the gas stream, which may prove difficult due to the high reactivity of tritium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗