SEARCH · Engineering Papers

Results for “configuration mixing”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2022

The Earth Model Column Collaboratory (EMC 2 ) v1.1: an open-source ground-based lidar and radar instrument simulator and subcolumn generator for large-scale models

Abstract. Climate models are essential for our comprehensive understanding of Earth's atmosphere and can provide critical insights on future changes decades ahead. Because of these critical roles, today's climate models are continuously being developed and evaluated using constraining observations and measurements obtained by satellites, airborne, and ground-based instruments. Instrument simulators can provide a bridge between the measured or retrieved quantities and their sampling in models and field observations while considering instrument sensitivity limitations. Here we present the Earth Model Column Collaboratory (EMC2), an open-source ground-based lidar and radar instrument simulator and subcolumn generator, specifically designed for large-scale models, in particular climate models, but also applicable to high-resolution model output. EMC2 provides a flexible framework enabling direct comparison of model output with ground-based observations, including generation of subcolumns that may statistically represent finer model spatial resolutions. In addition, EMC2 emulates ground-based (and air- or space-borne) measurements while remaining faithful to large-scale models' physical assumptions implemented in their cloud or radiation schemes. The simulator uses either single particle or bulk particle size distribution lookup tables, depending on the selected scheme approach, to perform the forward calculations. To facilitate model evaluation, EMC2 also includes three hydrometeor classification methods, namely, radar- and sounding-based cloud and precipitation detection and classification, lidar-based phase classification, and a Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) lidar simulator emulator. The software is written in Python, is easy to use, and can be straightforwardly customized for different models, radars, and lidars. Following the description of the logic, functionality, features, and software structure of EMC2, we present a case study of highly supercooled mixed-phase cloud based on measurements from the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE). We compare observations with the application of EMC2 to outputs from four configurations of the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE3) in single-column model (SCM) mode and from a large-eddy simulation (LES) model. We show that two of the four ModelE3 configurations can form and maintain highly supercooled precipitating cloud for several hours, consistent with observations and LES. While our focus is on one of these ModelE3 configurations, which performed slightly better in this case study, both of these configurations and the LES results post-processed with EMC2 generally provide reasonable agreement with observed lidar and radar variables. As briefly demonstrated here, EMC2 can provide a lightweight and flexible framework for comparing the results of both large-scale and high-resolution models directly with observations, with relatively little overhead and multiple options for achieving consistency with model microphysical or radiation scheme physics.

58 GEOSCIENCES↗

DOE OSTI2020

Comparing Cetane Number Measurement Methods

Cetane number is one of the most important fuel performance metrics for mixing controlled compression-ignition “diesel” engines, quantifying a fuel’s propensity for autoignition when injected into end-of-compression-type temperature and pressure conditions. The historical default and referee method on a Cooperative Fuel Research (CFR) engine configured with indirect fuel injection and variable compression ratio is cetane number (CN) rating. A subject fuel is evaluated against primary reference fuel blends, with heptamethylnonane defining a low-reactivity endpoint of CN = 15 and hexadecane defining a high-reactivity endpoint of CN = 100. While the CN scale covers the range from zero (0) to 100, typical testing is in the range of 30 to 65 CN. Alternatively, several constant-volume combustion chamber (CVCC)-based cetane rating devices have been developed to rate fuels with an equivalent derived cetane number (DCN) or indicated cetane number (ICN). These devices measure ignition delay for fuel injected into a fixed volume of high-temperature and high-pressure air to simulate end-of-compression-type conditions. In this study, a range of novel fuel compounds are evaluated across three CVCC methods: the Ignition Quality Tester (IQT), Fuel Ignition Tester (FIT), and Advanced Fuel Ignition Delay Analyzer (AFIDA). Resulting DCNs and ICNs are compared for fuels within the normal diesel fuel range of reactivity, as well as very high (~100) and very low DCNs/ICNs (~5). Distinct differences between results from various devices are discussed. This is important to consider because some new, high-efficiency advanced compression-ignition (CI) engine combustion strategies operate with more kinetically controlled distributed combustion as opposed to mixing controlled diffusion flames. These advanced combustion strategies may benefit from new fuel chemistries, but current rating methods of CN, DCN, and ICN may not fully describe their performance. In addition, recent evidence suggests ignition delay in modern on-road diesel engines with high-pressure common rail fuel injection systems may no longer directly correlate to traditional CN fuel ratings. Simulated end-of-compression conditions are compared for CN, DCN, and ICN and discussed in the context of modern diesel engines to provide additional insight. Results highlight the potential need for revised and/or multiple fuel test conditions to measure fuel performance for advanced CI strategies.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Refueling infrastructure planning in intercity networks considering route choice and travel time delay for mixed fleet of electric and conventional vehicles

The range anxiety has been a major factor that affects the market acceptance of electric vehicles. Even with the recent development of battery technologies, a lack of charging stations and range anxiety are still significant concerns, specifically for intercity trips. This calls for more investments in building charging stations and advancing battery technologies to increase the market share of electric vehicles and improve sustainability. This study suggests a configuration for plug-in electric vehicle charging infrastructure to support long-distance intercity trips of electric vehicles at the network level. A model is proposed to minimize the total system cost including infrastructure investment (building charging stations/spots) and travel time delays (charging time, waiting time in the queue, and detour time to access charging stations). This study fills existing gaps in the literature by capturing realistic patterns of travel demand and considering flow-dependent charging delays at charging stations. Furthermore, the proposed model, which is formulated as a mixed-integer program with nonlinear constraints, solves the optimization problem at the network level. At the network level, impacts of charging station locations on the traffic assignment problem with a mixed fleet of electric and conventional vehicles need to be considered. To this end, a traffic assignment module is integrated with a simulated annealing algorithm. The numerical experiments show a satisfactory application of the model for a full-scale case study (intercity network in Michigan). The solution quality and efficiency of the proposed solution algorithm are evaluated against those of an enumeration approach for a small case study. The results suggest that even for the current market share and charging stations’ setting, a significant investment is needed to support intercity trips without range anxiety issues and with acceptable delays. Additionally, through sensitivity analyses, the required infrastructure and battery investments to support intercity trips with acceptable delays are established for hypothetical increased market shares and battery size in the future.

42 ENGINEERING↗

DOE OSTI2025

Mechanisms and stability of Li dynamics in amorphous Li-Ti-P-S-based mixed ionic–electronic conductors: A machine learning molecular dynamics study

Mixed ionic–electronic conductors (MIECs) exhibit both high ionic and electronic conductivity to improve the battery performance. In this work, we investigate the mechanism and stability of transport channels in our recently developed MIEC material, amorphous Ti-doped lithium phosphorus sulfide (LPS), using molecular dynamics (MD) simulations with a 99% accurate machine-learning force field (MLFF) trained on ab initio MD data. The achieved MLFF helps efficient large-scale MD simulations on LPS with three Ti concentrations (10%, 20%, and 30%) and six temperatures (25°C to 225°C) to calculate ionic conductivity, activation energy, Li-ion transport mechanism, and configurational entropy. Results show that ionic conductivities and activation energies are consistent with our recent experimental values. Moreover, Li-ion transport occurs via free-volume diffusion facilitated by the formation of disordered Li-S polyhedra. The enhanced stability of transport channels at 10% and 20% Ti doping, compared to 0% and 30%, is observed by analyzing the vibrational and configurational entropy of these disordered Li-S polyhedra. Overall, this study highlights the utility of MLFF-based large-scale MD simulations in explaining the transport mechanism and the stability of Li-ion in Ti-doped LPS electrolyte with significant computational efficiency.

And configuration entropy↗

DOE OSTI2024

Optimizing dynamic wireless charging for electric buses: A data-driven approach to infrastructure planning

The network configuration significantly impacts the performance of dynamic wireless charging (DWC) technology for electric buses. Here, this study presents a novel approach to planning charging infrastructure for public transit using data-driven nonconvex mixed-integer optimization. Integrating DWC and charging station technologies reveals a trade-off between enroute and stationary charging times. Our framework optimizes bus frequency settings and transmitter coil arrangements to minimize operational and infrastructure costs. A case study in Chattanooga, Tennessee, demonstrates the method's effectiveness in mitigating range anxiety and reducing charging expenses. This research implies that integrating DWC technology into public transit systems can enhance the feasibility and cost-effectiveness of electric bus operations, promoting sustainable urban mobility.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Pb:Sn Ratio-Driven Polytypism and Band Modulation in Photoresponsive Hexagonal Perovskitoids

Three-dimensional (3D) perovskites of the formula AMX 3 are known for their excellent optoelectronic properties, but their design is limited by the narrow range of A-site cations that can template the 3D corner-sharing structure, and many of the viable options have already been explored. These materials also face structural instability under environmental conditions. In contrast, 3D hexagonal perovskitoids, with the same chemical formula, may offer enhanced stability and richer structural diversity through a range of corner- and face-sharing octahedral configuration options, providing greater opportunities for structural design; however, challenges such as synthesis complexity and wide band gaps have hindered optoelectronic performance achieved to date. Herein, we synthesized a structural homologous series of mixed-metal hexagonal perovskitoids with the formula APb 1-x Sn x I 3 (x = 0, 0.25, 0.50, 0.75, 1; A = ethylammonium, guanidinium) and identified three polytypes (9R, 12R, 6H) using single-crystal X-ray diffraction (SCXRD). These structures exhibited increasing corner-sharing connectivity with higher Sn content, revealing a previously unobserved relationship between metal composition and structural evolution in perovskitoid materials. The incorporation of Sn reduced the band gap (tunable from 2.51 eV to 1.87 eV) and drove structural transformations, a trend seen also in density functional theory (DFT) calculations, which suggest a thermodynamic preference for Pb at face-sharing sites and Sn at corner-sharing sites. DFT band structure calculations and optical spectroscopy also reveal an anomalous behavior in these materials, which we term polytypic band modulation. This phenomenon combines conventional band bowing with the structural transformations that occur as Sn content increases, as observed in the band gap and in photoluminescence spectra. Photodiodes fabricated from thin films of these materials exhibited stable and pronounced photoresponses across various light intensities over time. Furthermore, the combination of templating 3D perovskitoids with multiple cations and alloying Pb and Sn, suggests a vastly underexplored phase space that offers new parameters to tune perovskitoids.

Cations↗

DOE OSTI2022

Coupling coarse-mesh CFD with fine-mesh CFD for modeling molten-salt reactors in the Virtual Test Bed (VTB)

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS's coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS's high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Method and apparatus for in situ synthesis of SiC, SiC ceramic matrix composites, and SiC metal matrix composites during additive manufacturing

Methods and apparatuses for in situ synthesis of SiC, CMCs, and MMCs are disclosed, comprising: providing an apparatus having: an electromagnetic energy source; an autofocusing scanner; a powder system for SiC and one or more powders; a powder delivery system; a shielding gas comprising argon and/or nitrogen; and a computer coupled to and configured to control the energy source, scanner, powder system, and powder delivery system to deposit layers of the sample; programming the computer with specifications of the sample; using the computer to control electromagnetic radiation, mixing ratio, and powder deposition parameters based on the specifications of the sample; and using the autofocusing scanner to focus and scan the electromagnetic radiation onto the sample while the powders are concurrently deposited by the powder delivery system onto the sample to create a melting pool to deposit one or more layers onto the sample. Other embodiments are described and claimed.

Liu, Jian↗

DOE OSTI2022

Coupling coarse mesh CFD with fine mesh CFD for modeling for Molten Salt Reactors in the Virtual Test Bed VTB

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS’s coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS’s high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Understanding Lithium Local Environments in LiMn 0.5 Ni 0.5 O 2 Cathodes: A DFT-Supported 6 Li Solid-State NMR Study

LiMn 0.5 Ni 0.5 O 2 has been plagued by the structural failure driven by the irreversible phase change upon subsequent cycling and the structural degradation caused by the seemingly unavoidable Li/Ni exchange. This begs for a better understanding of the correlation between lithium local structures, synthesis conditions, and overall electrochemical properties of materials. In this study, we use density functional theory (DFT) to assist experimental 6 Li solid-state magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy to understand the nature and change in lithium local structures with extended annealing times (168 vs 15 h) synthesized at 900 °C. To correlate the NMR peak changes after different synthesis conditions, we calculated the NMR spectra of representative low-energy model configurations of LiMn 0.5 Ni 0.5 O 2 composition–zigzag, row, and flower structures–as well as the influence of Ni/Mn ordering and Li/Ni exchange on each spectrum using a supercell that is large enough to allow us to consider the effect of different relative configurations of adjacent transition metal (TM) layers. The analysis based on a combination of the calculated NMR and the experimental spectra suggests that most configurations contributing to the experimental spectra are in fact composed of a blend of configurations higher in energy than the classical well-ordered ground-state structures. Further, the extended annealing of LiMn 0.5 Ni 0.5 O 2 slightly enhances Ni/Mn and Li/Mn orderings while reducing the Li/Ni mixing ratio. This study shows how DFT calculations are crucial in providing a better understanding of the experimental characterization data and therefore local environments and domain structures in oxide materials, such as Li-ion battery cathodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Search for heavy neutral leptons decaying into muon-pion pairs in the MicroBooNE detector

We present upper limits on the production of heavy neutral leptons (HNLs) decaying to $\mu \pi$ pairs using data collected with the MicroBooNE liquid-argon time projection chamber (TPC) operating at Fermilab. This search is the first of its kind performed in a liquid-argon TPC. We use data collected in 2017 and 2018 corresponding to an exposure of $2.0 \times 10^{20}$ protons on target from the Fermilab Booster Neutrino Beam, which produces mainly muon neutrinos with an average energy of $\approx 800$ MeV. HNLs with higher mass are expected to have a longer time-of-flight to the liquid-argon TPC than Standard Model neutrinos. The data are therefore recorded with a dedicated trigger configured to detect HNL decays that occur after the neutrino spill reaches the detector. We set upper limits at the $90\%$ confidence level on the element $\lvert U_{\mu4}\rvert^2$ of the extended PMNS mixing matrix in the range $\lvert U_{\mu4}\rvert^2<(6.6$-$0.9)\times 10^{-7}$ for Dirac HNLs and $\lvert U_{\mu4}\rvert^2<(4.7$-$0.7)\times 10^{-7}$ for Majorana HNLs, assuming HNL masses between $260$ and $385$ MeV and $\lvert U_{e 4}\rvert^2 = \lvert U_{\tau 4}\rvert^2 = 0$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria: A DFT cluster expansion study

Magnetic contributions have the potential to significantly influence predicted phase stability within alloy and mineral mixing phase diagrams, yet have been historically challenging to incorporate due to a significant increase to phase space sampling. Here, in this work, we employ a computational protocol that includes spin orientation as an additional configurational component within multi-component cluster expansions between magnetic and non-magnetic metal oxide alloys [calculated using density functional theory (DFT) and the generalized gradient approximation]. This approach was used to determine the effect of magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria from first principles. Two-component cluster expansions of the magnetic components of eskolaite and hematite were first performed showing the ability of this method to properly calculate their respective magnetic properties. Two-component cluster expansions were then performed for non-magnetic Al(III) and ferromagnetic Cr(III) and Fe(III), and phase diagrams were calculated for later comparison. Finally, a non-magnetic Al(III) and “up” and “down” magnetic configurations for anti-ferromagnetic Cr(III) and Fe(III) were performed. Magnetic contributions to the calculated phase diagram for the corundumeskolaite system were shown to be inconsequential but are vital for accurate determination of the corundum-hematite solvus.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems

Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

CI-MBPT and Intensity-Based Lifetime Calculations for Th II

Lifetime calculations of Th II J = 1.5 and 2.5 odd states are performed with configuration–interaction many-body perturbation theory (CI-MBPT). For many J = 2.5 states, lifetimes are quite accurate, but two pairs of J = 2.5 odd states and many groups of J = 1.5 states are strongly mixed, making theoretical predictions unreliable. To solve this problem, a method based on intensities is used. To relate experimental intensities to lifetimes, two parameters, one an overall coefficient of proportionality for transition rates and one temperature of the Boltzmann distribution of populations, are introduced and fitted to minimize the deviation between theoretical and intensity-derived lifetimes. For strongly mixed groups of states, the averaged lifetimes obtained from averaged transition rates were used instead of individual lifetimes in the fit. Close agreement is obtained. Then intensity branching ratios are used to extract individual lifetimes for the strongly mixed states. The resulting lifetimes are compared to available directly measured lifetimes and reasonable agreement is found, considering limited accuracy of intensity measurements. The method of intensity-based lifetime calculations with fit to theoretical lifetimes is quite general and can be applied to many complex atoms where strong mixing between multiple states exists.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗

DOE OSTI2022

Experimental and computational investigation of the influence of iso -butanol on autoignition of n -decane and n -heptane in non-premixed flows

Experimental and computational investigation is carried out to elucidate the influence of iso-butanol on critical conditions of autoignition of n-decane and n-heptane. The counterflow configuration is employed. In this configuration an axisymmetric stream of air is directed over the surface of an evaporating pool of a liquid fuel. A stagnation plane is established. The temperature of the air is increased until autoignition is achieved in the mixing layer in the vicinity of the stagnation plane. The temperature of the air stream at autoignition, T ig , is measured at various values of the strain rate, which is defined as the axial gradient of the axial component of the flow velocity at the stagnation plane. Fuels tested are n-decane, n-heptane, iso-butanol and various mixtures of n-decane/iso-butanol and n-heptane/iso-butanol. Kinetic modeling is carried out using the recently updated version of the comprehensive CRECK chemical-kinetic mechanism. Critical conditions of autoignition are predicted for all fuels and fuel mixtures and the results are compared with the measurements. Computations show that low-temperature chemistry plays a significant role in promoting autoignition of n-decane and n-heptane, and the influence of low-temperature chemistry decreases with increasing strain rates because there is insufficient residence for the low-temperature reactions to take place. Experimental data and numerical simulations show that addition of even small amounts of iso-butanol to n-decane or n-heptane increases the value of T ig at low strain rates indicating that iso-butanol strongly inhibits the low-temperature chemistry of n-decane and n-heptane. Furthermore, the simulations show that when iso-butanol is present in the liquid fuel, only a small amount of n-decane is available in the gas phase and its low concentration cannot effectively sustain the low temperature degenerate branching oxidation process. Here, predicted flame structures show that the peak values of mole fraction of ketohydroperoxide are significantly reduced when iso-butanol is added to n-decane indicating that the kinetic pathway to low temperature ignition is blocked. This observation is confirmed from sensitivity analysis.

42 ENGINEERING↗