Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “induction”

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.

At least 757 records · Page 42

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗

A comprehensive and fair comparison of two neural operators (with practical extensions) based on $\mathrm{FAIR}$ data

Neural operators can learn nonlinear mappings between function spaces and offer a new simulation paradigm for real-time prediction of complex dynamics for realistic diverse applications as well as for system identification in science and engineering. Herein, we investigate the performance of two neural operators, which have shown promising results so far, and we develop new practical extensions that will make them more accurate and robust and importantly more suitable for industrial-complexity applications. The first neural operator, DeepONet, was published in 2019 (Lu et al., 2019), and its original architecture was based on the universal approximation theorem of Chen & Chen (1995). The second one, named Fourier Neural Operator or FNO, was published in 2020, and it is based on parameterizing the integral kernel in the Fourier space. DeepONet is represented by a summation of products of neural networks (NNs), corresponding to the branch NN for the input function and the trunk NN for the output function; both NNs are general architectures, e.g., the branch NN can be replaced with a CNN or a ResNet. According to Kovachki et al. (2021), FNO in its continuous form can be viewed conceptually as a DeepONet with a specific architecture of the branch NN and a trunk NN represented by a trigonometric basis. In order to compare FNO with DeepONet computationally for realistic setups, we develop several extensions of FNO that can deal with complex geometric domains as well as mappings where the input and output function spaces are of different dimensions. We also develop an extended DeepONet with special features that provide inductive bias and accelerate training, and we present a faster implementation of DeepONet with cost comparable to the computational cost of FNO, which is based on the Fast Fourier Transform. Here we consider 16 different benchmarks to demonstrate the relative performance of the two neural operators, including instability wave analysis in hypersonic boundary layers, prediction of the vorticity field of a flapping airfoil, porous media simulations in complex-geometry domains, etc. We follow the guiding principles of FAIR (Findability, Accessibility, Interoperability, and Reusability) for scientific data management and stewardship. The performance of DeepONet and FNO is comparable for relatively simple settings, but for complex geometries the performance of FNO deteriorates greatly. We also compare theoretically the two neural operators and obtain similar error estimates for DeepONet and FNO under the same regularity assumptions.

42 ENGINEERING↗

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗

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↗

Investigation of impurity driven corrosion behavior in molten 2LiF-BeF 2 salt

A 1000 -h corrosion experiment of 316H stainless steel at 700 °C investigated the role of metal fluoride impurities CrF 2 and FeF 2 in LiF-BeF 2 (66-34 mole %) (FLiBe) salt. Inductively coupled plasma-optical emission spectroscopy (ICP-OES) demonstrated the ability to measure changes in the concentration of these fluorides in FLiBe salt. Cyclic voltammetry (CV) showed the potential for real-time quantitative concentration measurement of corrosion products in the molten salt. Materials characterization revealed that both impurity species influence the dissolution of Cr from the alloy. The collective results from these analysis methods advance the understanding of the influence impurities have on materials corrosion.

2LiF-BeF2↗

Fouling characterization and treatment of water reuse concentrate with membrane distillation: Do organics really matter

Membrane distillation (MD) for the treatment of concentrated brines has been limited in part by membrane fouling, resulting in subsequent flux decline and membrane wetting. This study provides new insight into the identification of fouling and scaling mechanisms and pretreatment strategies for mitigating flux decline with MD treatment of water reuse reverse osmosis concentrate (ROC). Bench-scale direct contact MD experiments were performed with untreated and pretreated ROC. Biological activated carbon (BAC), chemical water softening, or fluidized bed crystallization reactor coupled with ion exchange (FBCR-IX) were selected as pretreatment strategies to isolate the effects of organic fouling and calcium scaling. Organic and inorganic compounds were analyzed by high-performance liquid chromatography (HPLC) and inductively coupled plasma mass spectrometry (ICP-MS). Calcium ions were found to be the major contributor to flux decline despite the high organic content in the ROC. Minimal organic fouling is likely because the organic matter in the ROC is hydrophilic, limiting hydrophobic-hydrophobic interactions between the organics and the membrane. Furthermore, the water flux declined by 63 % after removing organic compounds by BAC pretreatment, with 60 % of the calcium mass precipitating from the solution. Whereas, the water flux remained constant after removing multivalent ions with fluidized bed crystallization. Cleaning the membrane by acid washing and temperature reversal recovered 73 % and 12 % of the water flux, respectively. The analyses outlined in this study can assist in selecting appropriate fouling and scaling mitigation strategies for water reuse ROC and a wide range of feed solutions used in MD applications.

42 ENGINEERING↗

Electron and proton conductivity of Fe-N-C cathodes for PEM fuel cells: A model-based electrochemical impedance spectroscopy measurement

Impedance spectra of a PEM fuel cell with three Fe-N-C cathodes have been measured under the H 2 /N 2 testing regime. The spectra have been fitted using a recently developed physics-based impedance model, which takes into account finite proton ($σ_p$) and electron ($σ_e$) conductivity of the oxygen-free cathode catalyst layer. Fitting allowed to extract numerical data for $σ_p$, $σ_e$, the double layer capacitance, and the inductance of cables used for measuring impedance spectra. The values of σp and σe are close to what previously found for standard Pt/C electrodes, which is found for the first time using PGM-free catalysts. The method enables simultaneous measurement of reference proton and electron conductivity of PEMFC cathode.

25 ENERGY STORAGE↗

Active electrochemical interfaces stabilized through self-organized potential oscillations

Some electrochemical systems are known to have higher average efficiency when operated under an oscillatory regime. Given the compromise between activity and stability, the stability of electrochemical interfaces in a self-organized, oscillatory state must be taken into account. Here we evaluate the electro-oxidation of methanol and formic acid on platinum under regular and oscillatory conditions, and study the stability by following the Pt dissolution rates in situ with a stationary probe rotating disk electrode (SPRDE) coupled to an inductively coupled plasma mass spectrometer (ICP-MS). Generally speaking, as the electro-oxidation reaction proceeds, the platinum dissolution rate increases considerably. To guarantee Pt stability, the potential must be kept below 1.0 V vs. RHE. Interestingly, no dissolution is detectable when the electrode potential undergoes temporary self-organization, ensuring a stable and active interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of solvent formulations in electrolytes on fast charging of Li-ion cells

Improving the fast charging performance of lithium ion batteries (LIBs) has the promise to increase the widespread adoption of electric vehicles (EVs). Electrolyte development plays an important role in enabling fast charging. In this study, fast charging performance of LIBs is studied with different electrolytes of 1.2 M LiPF 6 in Ethylene Carbonate (EC)/Ethyl Methyl Carbonate (EMC)/co-solvents at 30/50/20 wt%. The co-solvents are methyl acetate (MA), ethyl acetate (EA), ethyl formate (EF), dimethyl carbonate (DMC) and EMC. Long term cycling performance under fast charging shows different capacity retention behaviors for different co-solvents. The structural changes in the electrode material are studied by X-ray absorption spectroscopy (XAS) and X-ray diffraction (XRD). The molarity changes in electrolyte is investigated by inductively coupled plasma-optical emission spectroscopy (ICP-OES). The electrode/electrolyte interfaces before and after fast charging are analyzed by X-ray photoemission spectroscopy (XPS). The characterization results are in good agreement with the long-term cycling performance. DMC shows the highest fast-charging capability among the five studied co-solvents due to its increased conductivity, improved electrode/electrolyte interface and stable electrode structural integrity.

25 ENERGY STORAGE↗

Aqueous solution-based synthesis approach for carbon-disordered rocksalt composite cathode development and its limitations

Disordered rocksalt cathodes exhibit high specific capacities and high energy density; however, their low electronic conductivity poses a great challenge. Herein, we explored an aqueous-solution-based synthesis route that involves controlling the surface charges of Li 1.2 Mn 0.6 Ti 0.2 O 1.8 F 0.2 (LMTOF) to be anchored by a few-layer reduced graphene oxide (rGO) for the first time. The uniform rGO wrapping on the surface of the LMTOF particles is achieved by electrostatic attraction between the negatively charged rGO and positively charged LMTOF particles. Although the initial specific capacity of rGO-LMTOF composite increased by 58 % compared to the pristine LMTOF, the composite experienced a severe capacity fade over cycling. The synthesis process in an aqueous medium resulted in Li + /H + exchange and TM dissolution as evidenced from inductively coupled plasmon analysis and X-ray diffraction analysis. Therefore, this work suggests the search for alternative media or conditions for the synthesis of carbon-disordered rock salt cathode composite.

25 ENERGY STORAGE↗

Scaling of energy delivered through an electrostatic discharge to a small series load

We study the energy delivered through a small-resistance series “victim” load during electrostatic discharge events in air. For gap lengths over 1 mm, the fraction of the stored energy delivered is mostly gap-length independent, with a slight decrease at larger gaps due to electrode geometry. The energy to the victim scales linearly with circuit capacitance and victim load resistance but is not strongly dependent on circuit inductance. This scaling leads to a simple approach to predicting the maximum energy that will be delivered to a series resistance for the case where the victim load resistance is lower than the spark resistance.

42 ENGINEERING↗

Cooperative fault management for resilient integration of renewable energy

Cooperative fault management (CFM) is designed herein to control different types of renewable energy resources cooperatively during electrical faults. This paper studies systems with a high penetration of photovoltaic (PV) energy and wind energy. First, CFM leverages power converters of PV farms to boost the ride-through capability of nearby doubly-fed induction generators (DFIGs). By controlling PV farms’ output voltages to change smoothly during both fault initiation and fault clearance, the widely used crowbar in DFIGs is less likely to be activated. Crowbar activation adversely makes DFIGs lose controllability and absorb reactive power. The second contribution is the development of a software-defined CFM controller and a controller in-the-loop demonstration of the real-time performance of this optimization-based CFM. CFM capitalizes on distributed optimization formulation to enable flexibility, plug-and-play, and privacy-preserving. Computation time, however, is a major concern for optimization-based dynamics control. Here, real-time controller-in-the-loop simulation results show optimization-based CFM can output reference values around 60 ms and is quick enough for dynamic control.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Correlating armature and needle dynamics with voltage waveforms of solenoid-actuated GDI injector

The injector voltage hump that appears near the needle closing has been used for the real-time monitoring and feedback control of fuel injection duration in modern engines. This voltage hump has been thought to result from the abrupt change in electromagnetic induction by the stoppage of needle motion but detailed electromagnetic processes and associated armature and needle dynamics during the needle closing have not been thoroughly investigated in a wide range of injection conditions, which knowledge is crucial for the delicate control of fuel injection based on the voltage hump. Here, the current study analyzes the transient armature and needle dynamics of a solenoid-actuated gasoline direct injection injector using an X-ray phase-contrast imaging technique. Then, the results are correlated with voltage waveforms during the needle closing transient under various injection pressures, injection pulse durations, and dwell times of split injections. The time derivatives of voltage waveforms showed lower and upper peaks in order in the regime of the voltage hump. Inconsistent with conventional understandings, the lower peak timing of the voltage derivative did not match with the timing of needle closing (end of injection) but rather matched with the abrupt descent timing of the armature and needle. The inflection timing and upper peak timing of the voltage derivative matched with the timings of actual needle closing and armature closing respectively. The amplitude of the voltage hump was near linearly dependent on the needle closing speed. The needle closing speed decreased upon the decrease of injection pulse duration and injection pressure which made it difficult to detect the voltage humps in ballistic injection regimes and low injection pressures. In split injection conditions, the voltage hump of the first injection was not detectable if the dwell time was shorter than the needle closing delay, the time from the current cut-off to the actual needle closing.

42 ENGINEERING↗

Linking chemical precursors to the synthesis of erythritol tetranitrate

Erythritol tetranitrate (ETN) is a homemade explosive with high explosive performance which can easily be synthesized using widely available chemical precursors. Ascertaining the manner in which ETN is synthesized and being able to trace the ETN back to its precursor materials may aid in the forensic identification. As a proof of concept, seven different types of ETN were synthesized with laboratory grade materials as well as commercially available chemicals from the hardware store and grocery store. The seven different types of solid ETN and the reaction quenches were analyzed using high performance liquid chromatography (HPLC). Partial least squares discriminant analysis (PLS-DA) was used to probe the different classes of ETN and reaction quenches. Several PLS-DA models were produced with the seven different types of ETN and quenches having their own classification, ETN and quenches based on their acid source, and ETN and the quenches based on their nitrate source. Regardless of which classification manner was used, there were no differences found between the different types of ETN or their reaction quenches. Furthermore, it is demonstrated that there is inherent variability in the synthesis of ETN making HPLC an unreliable method to determine the source of the synthesis materials. Elemental impurities were measured in seven solid ETN samples and their reaction quenches by inductively coupled plasma mass spectrometry (ICP-MS). No differences were discovered between the solid ETN samples, but there were trace metals found in the reaction quenches which were found to be identifiers of the source material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multimodal study of the impact of stimulation pH on shale pore structure, with an emphasis on organics behavior in alkaline environments

The tight nature of shale formations calls for hydraulic fracturing techniques able of altering the pore architecture to facilitate the production of hydrocarbon resources. The composition, pH, salinity, density and chemistry of Hydraulic Fracturing Fluids (HFF) vary substantially depending on reservoir characteristics. pH is among the most important stimulation fluid properties, ranging from very acidic, for carbonate reservoirs, to basic, in the case of clay-rich formations. pH regulates what geochemical reactions take place between the stimulation fluid and shale, as well as pore architecture alteration. The dissolution of carbonates and pyrite, and the precipitation of common minerals such as barite, gypsum and iron oxides in acidic environments have been extensively documented. In contrast, alkali stimulation environments and the role of organic components have received less attention. This research provides insight into the role of both minerals and organic components during alkaline stimulation, and the resulting pore architecture alterations. Here, a set of reactive experiments are performed using three different shales, with varying organic matter (OM) content, at different pH (3, 6, 8, 10, 12). The analysis of pore architecture alteration induced by the reaction was performed via nitrogen (N2) gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Inductively Coupled Plasma - Optical Emission Spectrometry (ICP-OES) was used to measure the composition of the stimulation fluid post-reaction. The system with a larger impact on the pore architecture (pH 12) underwent a more extensive analysis through the application of Focused Ion Beam Scanning Electron Microscopy (FIB-SEM) and x-ray Diffraction (XRD). Results show that mineral dissolution contributed to the creation of secondary porosity, mainly through dissolution of clays and silicates. We show that the abundance and distribution of organic matter both play a significant role in changes in pore architecture. Most samples showed an increase in inter-organic porosity, likely due to the organic acids within kerogen reacting with the stimulation fluid. Samples with a larger amount and more widespread distribution of organic matter show the most significant alterations. Finally, we also highlight the dissolution process of iron-rich minerals, including framboidal pyrite, as another important source of secondary porosity.

04 OIL SHALES AND TAR SANDS↗

Neural network model of neutral beam injection in the EAST tokamak to enable fast transport simulations

The neutral beam injection (NBI) system in EAST produces energetic neutral particles, which collide with electrons and ions in tokamak plasmas and heat the plasmas through Coulomb collisions. Moreover, it drives a non-inductive source of current, due to the charge-exchange collision between neutral particles and ions, and injects toroidal torque, which generates a toroidal rotation of the plasma. The effect caused by the NBI system, such as plasma heating, current drive, total neutron rate, momentum transfer, and shine-through, are modeled by a comprehensive module called NUBEAM. However, NUBEAM is computationally intensive since it relies on Monte Carlo methods. In this work, a neural network model has been developed as a surrogate model for NUBEAM in EAST. The database for neural-network model training, validation and testing is generated by running TRANSP for experimental discharges from recent EAST campaigns (after the latest NBI upgrade) while using the NUBEAM module. Simulation results illustrate that the trained neural network has the capability of replicating the predictions made by NUBEAM while demanding a significantly shorter execution time. Finally, these results indicate that surrogate models like the one proposed in this work could enable fast transport simulations for EAST after integrating them into a control-oriented predictive code such as COTSIM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗