Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transitional flow”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Observation of Double Excitations in the Resonant Inelastic X-ray Scattering of Nitric Oxide

The nitrogen K-edge resonant inelastic X-ray scattering (RIXS) map of nitric oxide (NO) has been measured and simulated to provide a detailed analysis of the observed features. High-resolution experimental RIXS maps were collected using an in situ gas flow cell and a high-transmission soft X-ray spectrometer. Accurate descriptions of the ground, excited, and core-excited states are based upon restricted active space self-consistent-field calculations using second order multiconfigurational perturbation theory. The nitrogen K-edge RIXS map of NO shows a range of features that can be assigned to intermediate states arising from 1s → π* and 1s → Rydberg excitations; additional bands are attributed to doubly excited intermediate states comprising 1s → π* and π → π* excitations. These results provide a detailed picture of RIXS for an open-shell molecule and an extensive description of the core-excited electronic structure of NO, an important molecule in many chemical and biological processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Magnetic skyrmion diode: Unidirectional skyrmion motion via symmetry breaking of potential energy barriers

In this work, we realize a magnetic skyrmion diode operated by a unidirectional skyrmion transport that flows in only one direction, which is highly significant for information processing in spintronic and nanoelectronic devices. We easily control the skyrmion transport by engineering asymmetric shapes of geometric structures. Here, we present a simple method to describe the underlying mechanism behind the unidirectional skyrmion transport by characterizing the topography of potential energy surfaces from a purely geometric perspective. Our approach enables a deeper physical insight into skyrmion transport manipulation and efficient design of skyrmion-based devices in geometric structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Evaluation of thermal energy storage for central receiver concentrating solar power plants under charging cycles through structural assessment

Electricity production by concentrated solar power (CSP) systems has stood out among energy transition alternatives due to their large generation capacity (between 30 MW and 400 MW) and high capacity factor (80%) when incorporating thermal storage systems (TES). Although TES tanks in operation worldwide are made of austenitic steels with remarkable mechanical and anticorrosive properties, operational failures have been reported due to low-cycle fatigue, creep, and the abrupt release of residual stresses generated during tank manufacturing. A critical operation for salt tanks is the initial daily charging operation, given the low fluid level and the thermal gradients between the inventory and the salt entering the tank. This study presents a structural simulation of a 39.6 m-diameter molten salt tank for central receiver CSP plants, accounting for temperature and pressure variations during charging and different sparger ring inlet configurations. The initial pre-stressed condition of the floor, resulting from the manufacture of welded plates, is used as a boundary condition. The fatigue and creep life analyses were performed in accordance with the ASME BPVC and API standards. It was found that the tank's floor is in a less favourable mechanical condition than the wall, due to its initial manufacturing condition and the maximum temperature differences during charging. A critical zone was identified in an area affected by residual stress on the floor and under the sparger ring. The baseline sparger ring configuration, with 52 2-inch-diameter orifices and a flow direction of 90 degrees , results in an average floor stress of approximately 25 MPa and creep damage after 1.2 years of tank operation.

14 SOLAR ENERGY↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Resource Characterization to Estimate Potential for Electricity Co-Production at Blackburn Oil Field, Nevada

The U.S. Department of Energy estimates that an annual average of 25 billion barrels of hot water are produced from oil and gas wells within the United States. The thermal energy available in the co-produced water stream is usually discarded, as the produced waters are considered an inconvenience by the operators and are disposed of using injection wells. However, utilizing organic Rankine cycle (ORC) generators, a vast amount of thermal energy can be captured and converted into electricity (albeit at relatively low efficiency due to the low temperatures). The National Renewable Energy Laboratory (NREL), in collaboration with Transitional Energy and Grant Canyon Oil & Gas, evaluated the feasibility of geothermal co-production of electricity by utilizing existing oil wells in Blackburn oil field in Nevada. The once prolific Blackburn oil field is located in Pine Valley, approximately 45 miles east-southeast of Elko, Nevada. Currently, the wells targeting the highly fractured Devonian Nevada dolomite reservoir are operating at a water cut ratio of more than 99%, with individual fluid (oil and water) production rates reaching 7.4 L/s (4,021 BBL/day). Analysis of publicly available data showed that the combination of the suitable wells' maximum historical production rates reached 22.90 L/s. The production from these wells occurs naturally and the wells are choked (and even shut down) by the operator to mitigate excessive water production, indicating a strong reservoir recharge and future opportunity to increase the water production for geothermal electricity generation. The main goal of this study was to evaluate the productivity of the existing wells, the performance of the reservoir, the surface network, and the operational constraints in order to achieve 1 MWe of electricity production from the field's water production. Utilizing the GEOPHIRES tool, we have determined that a twofold to threefold increase in the total fluid production, compared to the historical production under artificial restraint (choke), is required to reach a 1-MWe net target output for a low-temperature ORC system with air-cooled condensers. Lower flow rates would be required when utilizing water-based condensers instead of air-cooled condensers. However, that would require a constant supply of cold water, which may be challenging given the arid environment of the project site.

Blackburn↗
DOE OSTI2024

EPCAPE-PT-LANL Measurements: Single Particle Soot Photometer

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Single Particle Soot Photometer (Droplet Measurements Technology) Data Notes: Contact us if you want additional data products from this instrument. Reported data is the black carbon (rBC) number and mass concentration with diagnostic flags. The SP2 measures incandescence from particles that is induced with a 1064 Nd-YAG laser. Particles that absorb the laser energy and then emit radiation is assumed to contain black carbon. The amount of intensity of radiation is related to the mass of absorbing material in the particle. Here, we calibrate the incandescence intensity to size selected Regal Black (Cabot) nebulized from solution. Data provided for the EPCAPE campaign used the combined broadband high-gain + low-gain channels (BHBL). Each channel has a lower threshold of detection equivalent to 2-s of the respective channel noise. Thresholding has been applied to select between the high-gain and low-gain channels. The detection limit is 80-540 nm (Deq) or 0.36- 55 fg. It is assumed that rBC is the only aerosol type that is in significant concentration in the sampled atmosphere that absorbs the laser energy. Particle data are integrated over 10 second windows to calculate a rBC number and mass concentration. QC/QA: • Diagnostic flags that impacted measured concentration: - Sample, Sheath, and Purge Flow Rates: Despite observing fluctuations in all flows, the rBC detection and mass quantification are generally observed to be stable, although large changes in sample flow rate did impact detection efficiency. A flag was implemented for sample flow deviations >12 cm3/min from the set point averaged over 30 seconds - Laser Power: Detection efficiency decreases with laser power and deviations in laser power also affect mass quantification. Flag for laser power is set for deviations in laser current from the set value >5 mA. - Primary Detector Threshold: Primary thresholding is the signal value which determines whether a “particle” is recorded. Thresholding errors occur when the threshold value is too HIGH and real particles are ignored. • Several periods without data: - 16-18 Nov: Ultra Zero Air generator failure, flows deviated significantly from set points. rBC # conc recovered but questionable. - 20-21 Nov: Offline for calibrations for several hours each day. - 25-27 Nov: Data is missing. - 30 Nov: Power outage 3 Dec: SP2 hard-drive full. - 3 Dec: Power Outage Header: - BHBL_NumbConc[#/cc]: Refers to the number concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in particles per cubic centimeter. - BHBL_BCMass_Conc[ug/m3]: Refers to the mass concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in micrograms per cubic meter. - NoData_Flag[bool]: A boolean flag that indicates whether no data was recorded during a measurement. - NoBC_Flag[bool]: A boolean flag indicating whether no black carbon particles were recorded during the measurement. - Laser_Flag[bool]: A boolean flag indicating deviations in laser current during the measurement. - SampleFlow_Flag[bool]: A boolean flag indicating deviations from the set sample flow rate during the measurement. - PrimThresh_Flag[bool]: A boolean flag indicating deviations from the set primary threshold during the measurement, potentially ignoring real particles. - Manual_Flag[bool]: A boolean flag indicating manual intervention or adjustments during the measurement. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active or inactive during the measurement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formation of edge pressure pedestal and reversed magnetic shear due to toroidal rotation in a tokamak equilibrium

Toroidal rotation is well known to play significant roles in the edge transport and L-H transition dynamics of tokamaks. Our recent calculation finds that a sufficiently strong localized toroidal rotation can directly bring out the formation of edge pressure pedestal with reversed magnetic shear that is reminiscent of an H-mode plasma, purely through the effects of toroidal rotation on the tokamak MHD equilibrium itself. In particular, the enhanced edge toroidal rotation enables a substantial peaking of the parallel current profile near edge in higher β regimes, which leads to the flattening or reversal of the local q (safety factor) profile. Here the formation of pressure pedestal along with the reversed magnetic shear region is shown to be the natural outcome of the MHD tokamak equilibrium in a self-consistent response to the presence of a localized toroidal rotation typically observed in H-mode or QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electrochemical Ammonia Compression

An electrochemical (EC) compressor is a solid-state compression device. For decades, researchers have been studying EC compressors for applications for a variety of energy systems. As the world transitions away from fossil fuel energy sources, EC compression emerges as a promising technique for energy storage, specifically energy stored in the form of pressurized ammonia. Ammonia is also a commonly used refrigerant. The present studies examine the viability of EC compression for ammonia storage and refrigeration. EC ammonia compression increases the concentration of an ammonia-hydrogen mixture via the input of electrical energy and a series of chemical reactions. A polymer electrolyte membrane separates the low-concentration side from the high-concentration side of the device. On the low-pressure (anode) side hydrogen atoms oxidize and react with ammonia molecules, forming positively charge ammonium ions. The ions traverse the membrane electrolytically. The ammonium ions are reduced and revert back to hydrogen and ammonia upon reaching the high-pressure (cathode). An external circuit provides the current needed to sustain the reactions. In this project, we studied the performance of the ammonia EC compressor under a variety of different conditions. We replicated preliminary data and using a small cell with 5 cm2 of active area. We demonstrated the operation of larger cells with 100 cm2 of active area. Further, we used a commercially available hydrogen EC compressor stack to analyze the scaled-up compressor performance. While previous experiments examined only the transient EC compressor performance, we developed test facilities that allowed the compressor to reach steady state. We analyzed the effects of pressure and current on the EC compressor performance. We measured the flow rates of gas leaving the compressors and analyzed the composition using gas chromatography. We tested methods of separating ammonia from the effluent vapor, which contained hydrogen and water vapor. We found that back diffusion adversely affected the performance, especially when we maintained high pressure lifts.

25 ENERGY STORAGE↗
DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2023

Roles of Mn and Co in the Air Synthesizability of Layered Oxide Cathodes for Lithium-Based Batteries

High-nickel layered oxides (LiNi 1-x-y Mn x Co y O 2 ) are the prevailing cathode materials for high-energy-density lithium-based batteries, but they are plagued with deleterious surface air instabilities stemming from residual lithium formation. These issues severely hinder mass production as cathode calcination is limited to a flowing oxygen atmosphere, which entails high manufacturing costs as opposed to simpler and more economical air calcination. Here, while higher Ni contents are known to worsen air instabilities, the influence of Mn and Co contents on impacting these phenomena are less elucidated. We herein present the synthesis in ambient air and flowing oxygen atmospheres of three cathode variants with the same Ni contents, but varying Mn and Co contents: LiNi 0.7 Mn 0.3 O 2 , LiNi 0.7 Mn 0.15 Co 0.15 O 2 , and LiNi 0.7 Co 0.3 O 2 . It is found that the critical parameter influencing the air stability of the cathodes is the average Ni oxidation state, which is greatly dependent on the Mn and Co contents. Substitution of Mn for Ni drives down the Ni oxidation state as Mn exists as Mn4 + and reduces surface residual lithium formation, which vastly improves the overall air stability and, therefore, the synthesizability in air, but with a penalty of lowered capacity. In contrast, substitution of Co for Ni maintains Ni 3+ as Co exists as Co 3+ , offering increased initial capacity, but worsens the air stability and cyclability as the driving force for residual lithium formation and surface reactivity is increased.

25 ENERGY STORAGE↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ecosystem-Level Biomimicry for the Built Environment: Adopting Systems Ecology Principles for the Control of Heterogeneous Energy Systems

This paper presents, to our knowledge, the first system-level engineering study to bio-mimic the cybernetics and flow dynamics of energy resources in natural ecosystems for the control of heterogeneous energy infrastructures in the built environment. To this end, we introduce a novel Biomimetic Pulsing State (BPS) control that functionally mimics mature ecosystems. A preliminary Modelica-based case study features a single-family residential building with electrical and HVAC subsystems. The BPS control objective is to minimize the energy exchange between the building and the grid for the purposes of future self-supporting buildings and grid stability. The building contains PV, a wind turbine, a battery storage system, and a fan coil/heat pump HVAC system served by an ambient district energy network. Evaluating the control performance (BPS vs. constant setpoint) over several renewable energy scenarios (net importer, net zero, net exporter), simulation results show how the building's HVAC system can dynamically adjust its electrical load and temperatures to the electrical system's net energy status in real-time with BPS control. As a net importer, the heat pump consumed 29% less energy and its peak power reduced by 15% with BPS control compared to the constant setpoint case, with the zone air temperature 1 degree C lower on average. As a net exporter, the heat pump effectively consumed the same energy, but the peak power increased by 34% with BPS control, while the zone air temperature was 1 degree C higher when renewable energy was abundant, preheating the home. BPS and constant setpoint control produced comparable results under a net zero scenario. While further evaluation is essential, BPS control may help communities meet their sustainability and resiliency targets as they transition towards fully distributed and renewable energy grids.

biomimicry↗
DOE OSTI2022

The role of deviatoric stress and dislocations on the α to ω phase transformation in Ti

Under extreme conditions, α-Ti becomes unstable and transforms either into β-Ti at high temperature or into ω-Ti at high pressure. In what concerns the α to ω phase transformation (PT), there has been a wide range of experimentally reported transition pressures from approximately 2 to 15 GPa at room temperature. Deviatoric stresses and internal defects are often assumed to be the root cause of this variation. Here, in this study, these postulates are revisited using both continuum mechanics and molecular dynamics (MD) simulations. First, a simple continuum model, assuming linear elasticity and isotropic plasticity, is developed to describe the effects of applied stress and dislocations on the stability of an ω nucleus in an infinite α domain. Second, a new MD simulation method is developed to generate an ω nucleus in the α domain utilizing the displacement field identified from the topological analysis. Results from MD simulations show that despite the fact that phase diagrams typically delineate the limits between two phases in terms of only P and T, deviatoric stress promotes the α to ω phase transformation by reducing the critical radius above which an ω nucleus is stable. Furthermore, the required deviatoric stress to nucleate and stabilize a nanoscale ω nucleus is likely emanating from the internal stress of defects such as dislocations. The MD-informed micromechanics models are used to identify favorable configurations where dislocations help favor the α to ω transformation. These configurations show that the interaction with a basal or prismatic dislocation reduces the critical radius of a ω nucleus by about 10 or 16 %, respectively. In addition, prismatic edge dislocations are found to promote the growth of ω nucleus when interacting with the ($\bar{1}$$\bar{1}$20) α //(0001) ω interface. Importantly, a simple model of the arrival of dislocation at an ω nucleus suggests that PT does not necessarily require a pile-up to be present but could alternatively be mediated by a constant rapid flow of dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

System Integration of Rationally Designed Dilute Alloy Catalysts for Energy-Efficient Electrochemical CO 2 -To-Fuel Conversion

Energy effiecient electrochemical conversion of CO 2 to high-demand chemicals and transportation fuels using renewable solar and wind energy is a key technology needed for a high-productivity, low-carbon future. However, the development of scalable, low-cost, active, selective, and stable electrocatalysts remains a key challenge that needs to be overcome to enable high-volume conversion of CO 2 to feedstock chemicals for the chemical industry. In previous work, Lawrence Livermore National Laboratory (LLNL) has developed a rational design platform for dilute alloy transition metal electrocatalysts that promise to make electrochemical CO 2 conversion more energy efficient and selective. In this project, we worked with our industrial partners, Twelve and TotalEnergies, to improve scale up, integration, and stability of LLNL’s dilute alloy catalyst technology into an industry-relevant zero gap electrolyzer platform. Through virtual experiments and data analysis, we designed efficient and cost-effective copper-based catalysts. The catalyst was specifically designed to streamline the slowest and most energy-intensive step of the electrochemical chemical transformation of CO 2 to multi-carbon products – that is making the carbon-carbon bond by dimerization of the reaction intermediate carbon monoxide - resulting in up to 10% improvement in energy efficiency for C 2 products while simultaneously increasing the selectivity towards C 2 products. We tested two different scalable catalyst coating technologies and down-selected magnetron sputtering as the technique that provided the best control over catalyst loading, composition, and morphology. Using this technology, we successfully demonstrated integration of our dilute alloy catalysts into a 100 cm 2 electrolyzer platform with Faradaic efficiencies for ethylene production reaching 40% at a current density of 200 mA/cm 2 . We also developed the technology to integrate a well-defined nanoscale porosity by depositing alloy compositions that were compatible with dealloying, that is, selective leaching of an alloy component to generate nanoscale porosity. We observed that integration of the dealloying-derived nanoporosity improved catalyst performance and stability by leading to a more hydrophobic catalyst/anion exchange membrane interface. Unsolved problems that still need to be addressed are corrosion of the Cu catalyst -specifically if the used catalyst is exposed to air - as well as long term stability do to salt formation/deposition and flooding of the catalyst/electrolyzer flow channels, especially at higher current densities. As we only worked on optimization of catalyst composition, coating thickness, and morphology, further performance optimization will require a system level approach that includes optimization of electrolyzer design and membrane technology.

36 MATERIALS SCIENCE↗