SEARCH · Engineering Papers

Results for “local dynamic operator”

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 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS2002

An Exploration of Discontinuous Time Synchronous Averaging for Helicopter HUMS Using Cruise and Terminal Area Vibration Data

Recent research using NASA Ames AH-1 and OH-58C helicopters, and NASA Glenn test rigs, has shown that in-flight vibration data are typically non-stationary [l-4]. The nature and extent of this non-stationarity is most likely produced by several factors operating simultaneously. The aerodynamic flight environment and pilot commands provide continuously changing inputs, with a complex dynamic response that includes automatic feedback control from the engine regulator. It would appear that the combined effects operate primarily through an induced torque profile, which causes concomitant stress modulation at the individual internal gear meshes in the transmission. This notion is supported by several analyses, which show that upwards of 93% of the vibration signal s variance can be explained by knowledge of torque alone. That this relationship is stronger in an AH-1 than an OH-58, where measured non-stationarity is greater, suggests that the overall mass of the vehicle is an important consideration. In the lighter aircraft, the unsteady aerodynamic influences transmit relatively greater unsteady dynamic forces on the mechanical components, quite possibly contributing to its greater non-stationarity . In a recent paper using OH-58C pinion data [5], the authors have shown that in computing a time synchronous average (TSA) for various single-value metric computations, an effective trade-off can be obtained between sample size and measured stationarity by using data from only a single mesh cycle. A mesh cycle, which is defined as the number of rotations required for the gear teeth to return to their original mating position, has the property of representing all of the discrete phase angles of the opposing gears exactly once in the average. Measured stationarity is probably maximized because a single mesh cycle of the pinion gear occurs over a very short span of time, during which time-dependent non-stationary effects are kept to a minimum. Clearly, the advantage of local stationarity diminishes as the temporal duration of the cycle increases. This is most evident for a planetary mesh cycle, which can take several minutes to complete.

Huff, Edward M.↗

DOE OSTI2011

3D CFD ELECTROCHEMICAL AND HEAT TRANSFER MODEL OF AN INTERNALLY MANIFOLDED SOLID OXIDE ELECTROLYSIS CELL

A three-dimensional computational fluid dynamics (CFD) electrochemical model has been created to model high-temperature electrolysis cell performance and steam electrolysis in an internally manifolded planar solid oxide electrolysis cell (SOEC) stack. This design is being evaluated at the Idaho National Laboratory for hydrogen production from nuclear power and process heat. Mass, momentum, energy, and species conservation and transport are provided via the core features of the commercial CFD code FLUENT. A solid-oxide fuel cell (SOFC) model adds the electrochemical reactions and loss mechanisms and computation of the electric field throughout the cell. The FLUENT SOFC user-defined subroutine was modified for this work to allow for operation in the SOEC mode. Model results provide detailed profiles of temperature, operating potential, steam-electrode gas composition, oxygen-electrode gas composition, current density and hydrogen production over a range of stack operating conditions. Single-cell and five-cell results will be presented. Flow distribution through both models is discussed. Flow enters from the bottom, distributes through the inlet plenum, flows across the cells, gathers in the outlet plenum and flows downward making an upside-down ''U'' shaped flow pattern. Flow and concentration variations exist downstream of the inlet holes. Predicted mean outlet hydrogen and steam concentrations vary linearly with current density, as expected. Effects of variations in operating temperature, gas flow rate, oxygen-electrode and steam-electrode current density, and contact resistance from the base case are presented. Contour plots of local electrolyte temperature, current density, and Nernst potential indicate the effects of heat transfer, reaction cooling/heating, and change in local gas composition. Results are discussed for using this design in the electrolysis mode. Discussion of thermal neutral voltage, enthalpy of reaction, hydrogen production, cell thermal efficiency, cell electrical efficiency, and Gibbs free energy are discussed and reported herein.

08 HYDROGEN↗

DOE OSTI2025

Understanding the Dissolution and Passivation of an Aluminum Electrode during Electrocoagulation of Groundwater Using Neutron and X-ray Reflectometry

An aluminum (Al)-based electrocoagulation (EC) system can effectively remove dissolved silica and hardness in groundwater. The effectiveness of Al-EC in terms of pollutant removal, Faradaic efficiency, and energy consumption depends on the interfacial electrolysis or passivation of the electrode in water. Thus, understanding the electrolysis reaction at the liquid/electrode interface during operation is important for sustainable EC deployment. Here, a continuous flow-through Al-EC system was tested with various groundwater simulants, i.e., chloride (Cl – )-based, sulfate (SO 4 2– )-based, and mixed solutions. High pollutant removal with low energy consumption was observed in Cl – -based groundwater treatment, while low pollutant removal with high energy consumption was observed in SO 4 2– -based groundwater. For example, the required energy per unit mass of Al dosing in SO 4 2– -based groundwater is three times higher than that in Cl – -based groundwater at 10 mA/cm 2 . However, increasing the Cl – concentration significantly reduces this energy demand. In SO 4 2– -based groundwater, the silicate removal efficiency drops from 85.1% to 24.0% compared to that for Cl – -based groundwater, while Mg 2+ and Ca 2+ removal efficiencies decrease to 0.6% from 15.8% and 5.7% from 44.8%, respectively. To better understand this EC performance, we used in situ neutron reflectometry (NR) to examine the interfacial dynamics of Al dissolution and passivation at a 100 nm scale occurring on the surface of the sacrificial Al electrodes during EC. Ex situ X-ray reflectometry (XRR) was also used to support the in situ NR results. Both NR and XRR results revealed that Al dissolution is influenced by the presence of Cl – in the simulants, while a passivating layer forms on the electrode in a SO 4 2– -based solution. In the Cl – -based solution, anodic Al dissolution occurred locally and inhomogeneously across the surface of the Al anode film, resulting in a localized thickness reduction over time. In the SO 4 2– -based solution, no apparent dissolution of the Al anode was identified. Instead, Al underwent oxidation, forming an amorphous Al 2 O 3 surface layer within the Al electrode film that increased in thickness over time. In the mixed solution, both anodic Al dissolution and surface Al 2 O 3 layer formation occurred, indicating that Al dissolution and surface Al 2 O 3 layer formation are attributable to the Cl – and SO 4 2– ions, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Investigation on the Core Bypass Flow in a Very High Temperature Reactor

Uncertainties associated with the core bypass flow are some of the key issues that directly influence the coolant mass flow distribution and magnitude, and thus the operational core temperature profiles, in the very high-temperature reactor (VHTR). Designers will attempt to configure the core geometry so the core cooling flow rate magnitude and distribution conform to the design values. The objective of this project is to study the bypass flow both experimentally and computationally. Researchers will develop experimental data using state-of-the-art particle image velocimetry in a small test facility. The team will attempt to obtain full field temperature distribution using racks of thermocouples. The experimental data are intended to benchmark computational fluid dynamics (CFD) codes by providing detailed information. These experimental data are urgently needed for validation of the CFD codes. The following are the project tasks: (1) Construct a small-scale bench-top experiment to resemble the bypass flow between the graphite blocks, varying parameters to address their impact on bypass flow. Wall roughness of the graphite block walls, spacing between the blocks, and temperature of the blocks are some of the parameters to be tested. (2) Perform CFD to evaluate pre- and post-test calculations and turbulence models, including sensitivity studies to achieve high accuracy. (3) Develop the state-of-the art large eddy simulation (LES) using appropriate subgrid modeling. (4) Develop models to be used in systems thermal hydraulics codes to account and estimate the bypass flows. These computer programs include, among others, RELAP3D, MELCOR, GAMMA, and GAS-NET. Actual core bypass flow rate may vary considerably from the design value. Although the uncertainty of the bypass flow rate is not known, some sources have stated that the bypass flow rates in the Fort St. Vrain reactor were between 8 and 25 percent of the total reactor mass flow rate. If bypass flow rates are on the high side, the quantity of cooling flow through the core may be considerably less than the nominal design value, causing some regions of the core to operate at temperatures in excess of the design values. These effects are postulated to lead to localized hot regions in the core that must be considered when evaluating the VHTR operational and accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRSDate not supplied

In-Situ Resource Utilization Modeling of a Lunar Water Processing System

A key element of achieving a sustained surface presence, such as defined in NASA’s Artemis plan, is In-Situ Resource Utilization (ISRU). ISRU is the practice of using local resources to provide mission consumables that reduce system launch mass requirements, and regenerate resources (chiefly, water and oxygen) for propulsion and life support supporting both Lunar and Martian missions. ISRU systems require multiple complex processes, such as excavation, chemical reactors, and electrolysis subsystems that must operate in harmony to optimize the overall system process from beginning to end. The Mission Analysis and Integration Tool (MAIT) was previously developed with MATLAB in FY22 to connect individual subsystem models into a customized, flexible framework for the purpose of technology downselect, optimization, and end-to-end process planning. Beginning in FY24, MAIT was updated and became the capital program in the Systems Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) project. Prior work was leveraged and evolved using MATLAB/Simulink due to its ability to communicate with a vast number of other programming languages and makes up the backbone of data flow between inputs and outputs to the subsystem models. MAIT initially evaluated a suite of ISRU-related technologies, including the water processing Lunar Auger Dryer for ISRU (LADI) system with integrated upstream excavation and downstream electrolysis subsystems. With individual models consolidated, the MAIT tool generated over 5,000 cases during its first round of parametric sweeps on the water processing architecture at multiple production targets; the system analysis produced valuable insight into the optimal LADI geometry that minimized energy demands, estimated effects to cold trap size and radiator requirements, and calculated the power dynamics of the electrolysis unit and liquid oxygen storage volume. Additional efforts are being made to demonstrate the ability to scale ISRU technologies supporting the Space Technology Mission Directorate’s (STMD) commercialization strategy and increase the MAIT software capability. Work is ongoing to handle a wide array of ISRU system models beyond the Lunar environment, e.g. production of propellant for a Martian lander.

Avery Carlson↗

NASA NTRSDate not supplied

In-Situ Resource Utilization Modeling of a Lunar Water Processing System

Avery Carlson↗

NASA NTRSDate not supplied

Astrobee's Multi-year Activities at the International Space Station's Japanese Experimental Module

The Astrobee free-flying robots recently completed their third successful year of operations, housed in the Japanese Experimental Module (JEM) on the International Space Station. We summarize the three years of operation, giving special attention to JAXA's 1st and 2nd Kibo Robot Programming Challenge (RPC) and the mapping processes and tools that make Astrobees' autonomous operation possible. The JEM is an ever changing, dynamic environment where light settings, cargo, payloads, and crew members constantly move and interact with one another. The 1st JAXA Kibo RPC event, a collaboration between JAXA and NASA, was held in 2020. Students from several countries in the Asia-Pacific region competed in programming challenges with a simulated Astrobee. The finalists were then invited to run their code on an actual Astrobee in the JEM. For the final round, students programmed Astrobee to visit three different locations to obtain data that would instruct the robot to complete a final task with the participation of ISS crew. The first competition was a tremendous success, leading to an equally successful 2nd JAXA Kibo RPC in 2021 with even larger participation. The 3rd JAXA Kibo RPC will occur in 2022 expanding further to incorporate US participants. These activities led to several firsts in Astrobee’s history: operation of an Astrobee free-flying robot without crew supervision in preparation for on-orbit operations, autonomous image acquisition towards updates of the navigation map, non-NASA code running on the robot (both from JAXA and participating students), two heterogeneous free-flying robots from two different space agencies working together (Int-Ball and Astrobee) during the final event in 2020, the first payload using Astrobee, and having Astrobee controlled from a non-NASA location (Tsukuba Space Center). The preparation towards these activities involved constant evaluation of the different components of Astrobee's systems, specially mapping and localization. The paper describes the evolution of these systems such as the improvements made in localization to reduce localization drift by using graph-based optimization instead of the extended Kalman Filter localizer. Additionally, it reports on the mapping process and analysis tools created to validate map consistency across different activities in the constantly changing JEM environment. These enhancements have enabled the Astrobee facility to successfully execute over 100 ISS activities supporting over a dozen researchers and partners around the world.

Astrobee↗

DOE OSTI2016

On the origin of high ionic conductivity in Na-doped SrSiO 3

Understanding the local structure and ion dynamics is at the heart of ion conductor research. This paper reports on high-resolution solid-state 29 Si, 23 Na, and 17 O NMR investigation of the structure, chemical composition, and ion dynamics of a newly discovered fast ion conductor, Na-doped SrSiO 3 , which exhibited a much higher ionic conductivity than most of current oxide ion conductors. Quantitative analyses reveal that with a small dose (<10 mol%) of Na, the doped Na integrates into the SrSiO 3 structure to form Na x Sr 1-x SiO 3-0.5x , and with >10 mol% Na doping, phase separation occurs, leading to the formation of an amorphous phase β-Na 2 Si 2 O 5 and a crystalline Sr-rich phase. Variable-temperature 23Na and 17O magic-angle-spinning NMR up to 618 °C have shown significant changes in Na ion dynamics at high temperatures but little oxide ion motion, suggesting that Na ions are responsible for the observed high ionic conductivity. In addition, β-Na 2 Si 2 O 5 starts to crystallize at temperatures higher than 480 °C with prolonged heating, resulting in reduction in Na+ motion, and thus degradation of ionic conductivity. This study has contributed critical evidence to the understanding of ionic conduction in Na-doped SrSiO 3 and demonstrated that multinuclear high-resolution and high-temperature solid-state NMR is a uniquely useful tool for investigating ion conductors at their operating conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Evolutionary Optimization of a Charge Transfer Ionic Potential Model for Ta/Ta-Oxide Heterointerfaces

Heterostructures of tantalum and its oxide are of tremendous technological interest for a myriad of technological applications, including electronics, thermal management, catalysis and biochemistry. In particular, local oxygen stoichiometry variation in TaO x memristors comprising of thermodynamically stable metallic (Ta) and insulating oxide (Ta 2 O 5 ) have been shown to result in fast switching on the subnanosecond timescale over a billion cycles. This rapid switching opens up the potential for advanced functional platforms such as stateful logic operations and neuromorphic computation. Despite its broad importance, an atomistic scale understanding of oxygen stoichiometry variation across Ta/TaO x heterointerfaces, such as during early stages of oxidation and oxide growth, is not well understood. This is mainly due to the lack of a unified interatomic potential model for tantalum oxides that can accurately describe metallic (Ta), ionic (TaO x ) as well as mixed (Ta/TaO x interfaces) bonding environments simultaneously. To address this challenge, we introduce a Charge Transfer Ionic Potential (CTIP) model for Ta/Ta-oxide system by training against lattice parameters, cohesive energies, equations of state (EOS), elastic properties, and surface energies of the various experimentally observed Ta 2 O 5 polymorphs (hexagonal, orthorhombic and monoclinic) obtained from density functional theory (DFT) calculations. The best CTIP parameters are determined by employing a global optimization scheme driven by genetic algorithms followed by local Simplex optimization. Our newly developed CTIP potential accurately predicts structure, thermodynamics, energetic ordering of polymorphs, as well as elastic and surface properties of both Ta and Ta 2 O 5 , in excellent agreement with DFT calculations and experiments. We employ our newly parameterized CTIP potential to investigate the early stages of oxidation and atomic scale mechanisms associated with oxide growth on Ta surface at various temperatures. Furthermore, the CTIP potential developed in this work is an invaluable tool to investigate atomic-scale mechanisms and transport phenomena underlying the response of Ta/TaO x interfaces to external stimuli (e.g, temperature, pressure, strain, electric field etc.), as well as other interesting dynamical phenomena including the physics of switching dynamics in TaO x based memristors and neuromorphic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Local Analysis-Test Correlation Of Tow-Steered Composite Shells With Small Cutouts

The prebuckling and postbuckling behavior of two composite tow-steered shells with small cutouts is assessed using nonlinear dynamic finite element analyses and compared in detail with experimental measurements. The cylindrical shells were manufactured without cutouts using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the axially stiff crown and keel, to ±45 degrees on the shear-stiff sides. The first shell with overlaps has laminate thickness variations on the crown and keel that result from application of all 24 tows during each pass of the fiber placement system. The second shell without overlaps uses the fiber placement system’s tow drop/add capability to achieve a more uniform shell wall thickness. An unreinforced cutout representing a passenger door on a commercial aircraft fuselage is machined into the side of each of the two shells. These shells with cutouts were tested in axial compression and buckled elastically in previous work. Detailed nonlinear finite element analysis results are compared with their corresponding measured local load-displacement and load-strain responses in prebuckling, at global buckling, and into a stable postbuckled state. Test data from displacement transducers, strain gages, and digital image correlation are extracted at the centers of the crown and keel, and at the middles of the top and left edges of the cutout. The overall agreement between these measured and analytical responses is excellent in prebuckling through global buckling, and very good from global buckling through postbuckling. As such, the excellent correlation observed here increases confidence in applying tow-steered composites in operational vehicles.

Composites↗

NASA NTRSDate not supplied

Local Analysis-Test Correlation Of Tow-Steered Composite Shells With Small Cutouts

Composites↗

NASA NTRSDate not supplied

Multiple Probe Measurements at Uranus Motivated By Spatial Variability

Motivation: Spatial variations in the temperature field and composition of Uranus' atmosphere demonstrate a need for multiple entry probes to characterize vertical profiles in multiple locations. We will review variation of composition and temperature, which are produced by dynamical processes on a range of scales from global (polar anomalies, zonal bands) to regional (vortices, storms) [1–2]. In particular, the spatial variation of convective activity is not well understood based on existing remote sensing observations [3]. Understanding how these processes operate, and how they modulate variable composition, is key to constraining bulk atmospheric abundances. Abundances in turn provide cosmochemical constraints on planetary origins. Secondary probes at Uranus: Key measurements for secondary probes are temperature-pressure profiles, along with compositional profiles. Radio occultations are limited to shallow levels less than ~2 bar [4]. In situ temperature-pressure measurements (or atmospheric structure measurements) can unambiguously extend these results to deeper levels. Measurements with vertical resolution of at least 2 km are needed to characterize anomalies like the 1.2-bar feature from the Voyager 2 occultation at 2–6° S [4], which supports a range of temperature gradients depending on assumptions of composition. Simultaneous measurements of temperature, pressure, and composition are key to understanding this class of features. Species such as methane, hydrogen sulfide, and ammonia must be measured along a probe descent profile because they vary over several orders of magnitude due to the strong temperature dependence of their saturation vapor pressures [e.g., 5]. Results of these measurements can be interpreted to understand the potential for moist convective activity in the atmosphere. Lessons from the other giant planets: On Jupiter, the Galileo Probe's entry into a meteorologically distinct five-micron hot spot led many to interpret the local composition as column-stretched, so that well-mixed abundances were reached at deeper levels than in surrounding, unperturbed regions [6–8]. Ground-based microwave measurements and Juno data now indicate that the deep depletion of ammonia is a very widespread atmospheric characteristic not limited to 5-µm hot spots [9–11]. But in the absence of multiple probes, we do not know if the other volatiles H2S and H2O behave in the same way. The open questions for the Jupiter case strongly motivate sending multiple probes to Uranus. On Saturn, retrievals of NH3 and PH3 abundances at shallow levels vary with latitude [12], but it is unknown how deep these differences extend, which is why a Saturn probe (or probes) was a mission theme for NASA's New Frontiers 4 and 5 opportunities [13]. Challenges for secondary probes: Cost is an issue due to the perception that it involves sacrifices to other mission elements. Spacecraft trajectories may be constrained by needs for orbit insertion that limit probe deliveries to different latitudes, and additional limitations may be placed on communication windows for probe descent phases (particularly if multiple probes are released from the orbiter simultaneously) [14]. Finally, composition sensors for miniature secondary probes are not at the required technological maturity [15]. Mass spectrometers are typically too large, massive, and powerhungry, while smaller nanosensors are only beginning to be developed for planetary missions. Finally, the need for probe survival heating is most easily met by radioisotope heat sources, but these require regulatory approvals that are even more difficult to satisfy compared to standard environmental reviews [16], unless the secondary probe is designed as a core element of a mission. References: [1] Molter E.M., et al. (2021) PSJ, 2, 3. [2] Rowe-Gurney N., et al. (2021) Icar, 365, 114506. [3] Hueso R., Sánchez-Lavega A. (2019) SSRv, 215, 52. [4] Lindal G.F., et al. (1987) JGR, 92, 14987– 15001. [5] Simon A.A., et al. (2022) RemS, 14, 1518. [6] Atreya, S.K., et al. (1997) in The Three Galileos: The Man, the Spacecraft, the Telescope, pp. 249–260 (C. Barbieri et al., eds.). [7] Showman A.P., Ingersoll A.P. (1998) Icar, 132, 205–220. [8] Friedson A.J. (2005) Icar, 177, 1–17. [9] de Pater I., et al. (2001) Icar, 149, 66–78. [10] Li C., et al. (2017) GeoRL, 44, 5317–5325. [11] de Pater I., et al. (2019) Icar, 322, 168–191. [12] Fletcher L.N., et al. (2009) Icar, 202, 543–564. [13] National Research Council (2011) Vision and Voyages. [14] Sayanagi K.M., et al. (2020) SSRv, 216, 72. [15] Wong M.H. et al. (2021) BAAS, 53, 486. [16] Zide A., Mendoza-Hill A., Cheney D. (2022) COSPAR Abstracts H0.6-0012-22.

Michael H Wong↗

DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗

DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

DOE OSTI2016

Co-benefits of global and regional greenhouse gas mitigation for US air quality in 2050

Policies to mitigate greenhouse gas (GHG) emissions will not only slow climate change but can also have ancillary benefits of improved air quality. Here we examine the co-benefits of both global and regional GHG mitigation for US air quality in 2050 at fine resolution, using dynamical downscaling methods, building on a previous global co-benefits study (West et al., 2013). The co-benefits for US air quality are quantified via two mechanisms: through reductions in co-emitted air pollutants from the same sources and by slowing climate change and its influence on air quality, following West et al. (2013). Additionally, we separate the total co-benefits into contributions from domestic GHG mitigation vs. mitigation in foreign countries. We use the Weather Research and Forecasting (WRF) model to dynamically downscale future global climate to the regional scale and the Sparse Matrix Operator Kernel Emissions (SMOKE) program to directly process global anthropogenic emissions to the regional domain, and we provide dynamical boundary conditions from global simulations to the regional Community Multi-scale Air Quality (CMAQ) model. The total co-benefits of global GHG mitigation from the RCP4.5 scenario compared with its reference are estimated to be higher in the eastern US (ranging from 0.6 to 1.0 µg m -3 ) than the west (0–0.4 µg m -3 ) for fine particulate matter (PM 2.5 ), with an average of 0.47 µg m -3 over the US; for O 3 , the total co-benefits are more uniform at 2–5 ppb, with a US average of 3.55 ppb. Comparing the two mechanisms of co-benefits, we find that reductions in co-emitted air pollutants have a much greater influence on both PM 2.5 (96 % of the total co-benefits) and O 3 (89 % of the total) than the second co-benefits mechanism via slowing climate change, consistent with West et al. (2013). GHG mitigation from foreign countries contributes more to the US O 3 reduction (76 % of the total) than that from domestic GHG mitigation only (24 %), highlighting the importance of global methane reductions and the intercontinental transport of air pollutants. For PM 2.5 , the benefits of domestic GHG control are greater (74 % of total). Since foreign contributions to co-benefits can be substantial, with foreign O 3 benefits much larger than those from domestic reductions, previous studies that focus on local or regional co-benefits may greatly underestimate the total co-benefits of global GHG reductions. We conclude that the US can gain significantly greater domestic air quality co-benefits by engaging with other nations to control GHGs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗

DOE OSTI2023

HFTS-1 (Final Report)

Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.

02 PETROLEUM↗