SEARCH · Engineering Papers

Results for “criticality simulations”

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 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Ferroelectric phase-transition frustration near a tricritical composition point

Phase transition describes a mutational behavior of matter states at a critical transition temperature or external field. Despite the phase-transition orders are well sorted by classic thermodynamic theory, ambiguous situations interposed between the first- and second-order transitions were exposed one after another. Here, we report discovery of phase-transition frustration near a tricritical composition point in ferroelectric Pb(Zr 1-x Ti x )O 3 . Our multi-scale transmission electron microscopy characterization reveals a number of geometrically frustrated microstructure features such as self-assembled hierarchical domain structure, degeneracy of mesoscale domain tetragonality and decoupled polarization-strain relationship. Associated with deviation from the classic mean-field theory, dielectric critical exponent anomalies and temperature dependent birefringence data unveil that the frustrated transition order stems from intricate competition of short-range polar orders and their decoupling to long-range lattice deformation. With supports from effective Hamiltonian Monte Carlo simulations, our findings point out a potentially universal mechanism to comprehend the abnormal critical phenomena occurring in phase-transition materials.

42 ENGINEERING↗

DOE OSTI2026

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Assessment and Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY25 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., the 200 DV-1 Operable Unit), amendment particles need to be delivered a reasonable distance away from an injection well to provide cost effective in situ treatment. Field-scale particle transport models can estimate spatial deposition of amendment particles in the subsurface, which is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capability and a lack of experimental data to validate and parameterize particle transport models. During this fiscal year, the following progress has been made toward a field-scale particle transport modeling evaluation: (1) in addition to the two particle transport models implemented last FY, four additional particle transport models have been implemented within PFLOTRAN; (2) a Python-based pre-screening tool was finalized, enabling users to quickly estimate the particle radius of influence (ROI) for any given particle-amendment system; and (3) an initial compatibility assessment was completed using both particle transport simulations deployed through the pre-screening tool, in conjunction with general guidelines to (a) identify the most suitable amendment particle sizes for various Hanford sediments and (b) evaluate amendment-delivery fluid compatibility. The preliminary compatibility assessment revealed that for amendment delivery success to the various Hanford target formations, amendment particle sizes will likely need to be smaller than the amendment sizes tested in the DV-1 treatability study. It is recommended that amendment particles be decreased in size, or alternative smaller size amendments be obtained from the manufacture, prior to any further experimental testing. Also, preliminary testing suggests that xanthan gum may be the most broadly compatible delivery fluid. Planned laboratory experiments will be instrumental in validating and refining the PFLOTRAN particulate transport model formulations, ultimately enabling predictive capabilities to facilitate the design of field-scale amendment delivery systems. This work consists of acquiring new theoretical or experimental knowledge. The information associated with this report should not be used as design input or operating parameters without additional qualification.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Determining Spatial Scales of Soil Moisture—Cloud Coupling Pathways Using Semi-Idealized Simulations

The spatial scales of land-atmosphere interactions over heterogeneous, realistic soil moisture patterns are critical to improve our understanding of land-atmosphere interactions. However, these scales are poorly quantified. We use high-resolution numerical experiments to quantify the spatial scales for a selected day during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) campaign in Oklahoma. A semi-idealized approach (removing background winds and surface topography while retaining the real-world data for other boundary conditions) enables examination of the local influence from the surface. Simulations that apply uniform surface, variable land cover with uniform soil moisture, and variable land cover with variable soil moisture are then compared. Land cover variations, especially water and urban surfaces, create two variance peaks in the surface sensible heat flux in the scales >10 km (meso-β) and <6 km (meso-γ). Both peaks are enhanced by the observation-based soil moisture distribution that has most of its variance in the meso-β scale. Atmospheric response in the meso-β scale is hydrostatic with weak vertical motion. The land cover variations induce secondary circulations in the scales of 2.5–6 km. The additional soil moisture variability intensifies extreme circulations over dry soil patches and sharp sensible heat flux gradients, slightly extending the characteristic scales to 2.5–9 km. While their areal coverage is small, soil moisture variations on such “hotspots” lead to substantially larger convective cloud structure and domain-mean rainfall, calling for a consideration of their non-linear effect in large-scale models.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2026

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗

DOE OSTI2022

Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

CO2 reduction↗

DOE OSTI2022

Co-simulation Framework for Community-scale Building-grid Integration [SWR-21-75]

Distributed energy resources (DERs), including rooftop solar, energy storage, and flexible loads, are gaining popularity as costs decline and as building owners and utilities realize their benefits. DERs can improve distribution system efficiency, help prevent the need for expensive grid upgrades, and increase the resilience of local communities. However, they can also cause difficulties in grid operations and can require controls to achieve their benefits. To address this challenge, NREL researchers have developed a community-scale solution that assesses the impacts of DERs and their control strategies on a distribution system. The framework has been shown to reduce solar photovoltaic (PV) curtailment to 0%, mitigate the adverse impact of solar variability on the distribution voltage, and provide up to 5-day critical load support during emergency events. Utilizing 5 different modules representing the feeder, buildings, home energy management systems, an aggregator, and a utility controller, NREL expects this simulation technology to play a critical role in the continued integration of DERs. According to the Energy Information Administration (EIA), solar curtailments accounted for 94% of the total energy curtailed in the California Independent System Operator (CAISO) in 2020. By enabling Independent System Operators (ISOs) and utility operators to bring solar curtailments to 0%, the electrical grid can become less dependent on fossil-fueled power generation sources. NREL's co-simulation framework contains five major components: Distribution Feeder Model: describes the distribution feeder topology using OpenDSS, including the locations of all DERs. Residential Building Model: simulates a large number of buildings at a high resolution using OCHRETM. The model is equipped to control equipment based on signals from an external module. The model includes major household appliances such as HVAC and a water heater, non-dispatchable load models, a distributed PV system, and a home battery system. Home Energy Management System: optimizes the controls for the devices in a home using foreseeTM. The control can adjust based on the user preferences including cost, comfort, and convenience. In hierarchical control scenarios, where the houses follow signals from an aggregator, the home energy management system provides a flexibility band with a range of power and follows the dispatch signals received from aggregator. Community-Level Aggregator: solves for optimal energy dispatch based on the flexibility bands received from each home and the grid service signal received from the utility controller. Utility-Level Controller: provides grid signals for voltage control using Distributed Energy Resources (DERs), such as solar systems, in the community.

Balamurugan, Sivasathya Pradha↗

DOE OSTI2025

Criticality and Magnetic Phases of Ising Shastry–Sutherland Candidate Holmium Tetraboride

Frustrated magnetic systems arising in geometrically constrained lattices represent rich platforms for exploring unconventional phases of matter, including fractional magnetization plateaus, incommensurate orders and complex domain dynamics. However, determining the microscopic spin configurations that stabilize such phases is a key challenge, especially when in-plane and out-of-plane spin components coexist and compete. Here, we combine neutron scattering and magnetic susceptibility experiments with simulations to investigate the emergence of field-induced fractional plateaus and the related criticality in a frustrated magnet holmium tetraboride (HoB 4 ) that represents the family of rare earth tetraborides that crystalize in a Shastry–Sutherland lattice in the ab plane. We focus on the interplay between classical and quantum criticality near phase boundaries, as well as the role of material defects in the stabilization of the ordered phases. We find that simulations using classical annealing can explain certain observed features in the experimental Laue diffraction and the origin of multiple magnetization plateaus. Our results show that defects and out-of-plane interactions play an important role and can guide the route towards resolving microscopic spin textures in highly frustrated magnets.

Frustration↗

NASA NTRS1995

Managing Complexity in Multidisciplinary Visualization

As high performance computing technology progresses, computational simulations are becoming more advanced in their capabilities. In the computational aerosciences domain, single discipline steady-state simulations computed on a single grid are far from the state-of-the-art. In their place are complex, time-dependent multidisciplinary simulations that attempt to model a given geometry more realistically. The product of these multidisciplinary simulations is a massive amount of data stored in different formats, grid topologies, units of measure, etc., as a result of the differences in the simulated physical domains. In addition to the challenges posed by setting up and performing the simulation, additional challenges exist in analyzing computational results. Visualization plays an important role in the advancement of multidisciplinary simulations. To date, visualization has been used to aid in the interpretation of large amounts of simulation data. Because the human visual system is effective in digesting a large amount of information presented graphically, visualization has helped simulation scientists to understand complex simulation results. As these simulations become even more complex, integrating several different physical domains, visualization will be critical to digest the massive amount of information. Another important role for visualization is to provide a common communication medium from which the domain scientists can use to develop, debug, and analyze their work. Multidisciplinary analyses are the next step in simulation technology, not only in computational aerosciences, but in many other areas such as global climate modeling. Visualization researchers must understand and work towards the challenges posed by multidisciplinary simulation scenarios. This paper addresses some of these challenges, describing technologies that must be investigated to create a useful visualization analysis tool for domain scientists.

Miceli, Kristina D.↗

NASA NTRS1993

Control structural interaction testbed: A model for multiple flexible body verification

Conventional end-to-end ground tests for verification of control system performance become increasingly complicated with the development of large, multiple flexible body spacecraft structures. The expense of accurately reproducing the on-orbit dynamic environment and the attendant difficulties in reducing and accounting for ground test effects limits the value of these tests. TRW has developed a building block approach whereby a combination of analysis, simulation, and test has replaced end-to-end performance verification by ground test. Tests are performed at the component, subsystem, and system level on engineering testbeds. These tests are aimed at authenticating models to be used in end-to-end performance verification simulations: component and subassembly engineering tests and analyses establish models and critical parameters, unit level engineering and acceptance tests refine models, and subsystem level tests confirm the models' overall behavior. The Precision Control of Agile Spacecraft (PCAS) project has developed a control structural interaction testbed with a multibody flexible structure to investigate new methods of precision control. This testbed is a model for TRW's approach to verifying control system performance. This approach has several advantages: (1) no allocation for test measurement errors is required, increasing flight hardware design allocations; (2) the approach permits greater latitude in investigating off-nominal conditions and parametric sensitivities; and (3) the simulation approach is cost effective, because the investment is in understanding the root behavior of the flight hardware and not in the ground test equipment and environment.

Chory, M. A.↗

DOE OSTI2022

Modelling single atom catalysts for water splitting and fuel cells: A tutorial review

Here, in this tutorial review, we report the state-of-the-art of the modeling approaches of Single Atom Catalysts (SAC) for water splitting and fuel cells reactions. The discussion applies for Hydrogen Evolution Reaction (HER), Oxygen Reduction Reaction (OER), Hydrogen Oxidation Reaction (HOR), and Oxygen Reduction Reaction (ORR). The main scope of this work is to underline the relevant aspects of SACs modelling. On the one hand, the review could help computational chemists aiming to start the study of SACs. On the other hand, experimentalists could find the critical analysis of DFT results of interest to understand better the strengths and weaknesses of simulations, and how to interpret computational results. After an introductory section of SACs, we start by briefly presenting the state-of-the-art methodologies. Then, we analyze the critical aspects for a reliable prediction of the electronic properties, and we discuss the robustness of the structural models and ways to validate them. Furthermore, we discuss the main approaches to predict catalytic activity and selectivity, which is the final goal of the computational catalysis. We conclude this review with a critical analysis of the current challenges in the field, and the main limitations of the modeling approaches that are described.

25 ENERGY STORAGE↗

NASA NTRS1981

Onboard utilization of ground control points for image correction. Volume 2: Analysis and simulation results

An approach to remote sensing that meets future mission requirements was investigated. The deterministic acquisition of data and the rapid correction of data for radiometric effects and image distortions are the most critical limitations of remote sensing. The following topics are discussed: onboard image correction systems, GCP navigation system simulation, GCP analysis, and image correction analysis measurement.

Source record↗

DOE OSTI2020

Constrained Bayesian Optimization of Criticality Experiments at LLNL [Slides]

The design of criticality experiments is typically an iterative process that employs a Monte Carlo transport code. The goal is to find a design that optimizes some variable, like the sensitivity of a response to a cross section, while simultaneously ensuring criticality. The high fidelity of the Monte Carlo code is a great asset, but it makes exploring the design space computationally expensive. Herein, we present how a constrained Bayesian optimization algorithm can be used to efficiently design a criticality experiment. It uses Gaussian processes as a surrogate model to probe the design space and to reduce the number of code executions that are needed to find the optimum. We demonstrate constrained Bayesian optimization with a Pu-239/polyethylene solution system and a TEX experiment that is designed for critical ity safety validation of a nuclear waste model at the Hanford Site. For both systems, a global optimum was found within 75 Monte Carlo simulations.

42 ENGINEERING↗

NASA NTRS1984

Turbulent heating of heavy ions on auroral field lines

Electrostatic ion cyclotron turbulence and the associated acceleration of ions on auroral field lines are investigated analytically and by plasma simulations. The auroral plasma is assumed to consist of drifting electrons and stationary hydrogen and oxygen ions. It is found that for a given critical drift, the maximum ion perpendicular heating is generally larger for oxygen ions than for hydrogen ions. Simulation results show that unless the oxygen ions are a minority species, oxygen transverse heating generally exceeds that of hydrogen ions. Theory and numerical simulations are in good agreement.

Ashour-Abdalla, M.↗

NASA NTRS2011

Impact of Optimized land Surface Parameters on the Land-Atmosphere Coupling in WRF Simulations of Dry and Wet Extremes

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry (2006) and wet (2007) conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through the use of a new optimization and uncertainty module in NASA's Land Information System (LIS-OPT), whereby parameter sets are calibrated in the Noah land surface model and classified according to the land cover and soil type mapping of the observations and the full domain. The impact of the calibrated parameters on the a) spin up of land surface states used as initial conditions, and b) heat and moisture fluxes of the coupled (LIS-WRF) simulations are then assessed in terms of ambient weather, PBL budgets, and precipitation along with L-A coupling diagnostics. In addition, the sensitivity of this approach to the period of calibration (dry, wet, normal) is investigated. Finally, tradeoffs of computational tractability and scientific validity (e.g.,. relating to the representation of the spatial dependence of parameters) and the feasibility of calibrating to multiple observational datasets are also discussed.

Kumar, Sujay↗

DOE OSTI2024

Simulated nitrogen deposition and precipitation events alter microbial carbon cycling during early stages of litter decomposition

Plant litter decomposition is a major nutrient input to terrestrial ecosystems that is primarily driven by microorganisms. Litter quality is considered a key drive of decomposition; however, human-induced global disturbance like nitrogen deposition and increasing extreme precipitation events will shift nutrient availability during litter decomposition. Little is known about how shifting nutrient availability will impact dissolved organic matter concentrations and microbially driven carbon cycling that are critical to soil organic matter formation. This study investigated the effect of simulated nitrogen deposition and repeated precipitation events on microbially driven carbon flow during short-term litter decomposition using a ‘common garden’ experiment with microcosms containing sand and blue grama grass litter inoculated with different microbial communities. Overall, nitrogen deposition decoupled respiration and dissolved organic carbon (DOC) by increasing respiration and not affecting DOC concentrations. Moreover, nitrogen deposition had no effect on microbial carbon use efficiency (CUE). Repeated simulated precipitation events significantly increased DOC concentrations, decreased microbial CUE, increased the microbial metabolic quotient (qCO2), and altered microbial composition and diversity. These findings highlight the complex interactions and responses of surface litter decomposers to shifting nutrient availability and contradicts previous findings that nitrogen deposition will increase soil carbon sequestration from a larger supply of DOC and reduced respiration.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRSDate not supplied

Training Astronauts using Hardware-in-the-Loop Simulations and Virtual Reality

The commercial market has recently started giving significant attention to virtual and augmented reality, even though the technology has been around for many years. The Virtual Reality Training Lab (VRL) at the NASA Johnson Space Center has been using virtual reality to train astronauts for decades. This paper will focus on describing three major Hardware-in-the-Loop VR simulation systems, the Simplified Aid for EVA Rescue (SAFER) system known as the "jetpack", the Mass Handling System nicknamed Charlotte, and a simulated robotics environment for collaborative mission evaluation. Two of these systems are critical for astronaut training. Crew must certify on SAFER and go through the Charlotte Mass Handling training prior to flying to the International Space Station (ISS). Typically, they also complete at least one collaborative visualization session to review any planned Extra Vehicular Activities (EVAs), or spacewalks, before an assigned flight. Given the volatility of new technologies, the graphics and simulation environments used are maintained to be hardware agnostic to preserve a high level of fidelity. Utilizing VR for astronaut training has proved to be effective and essential for these specific systems.

Angelica D. Garcia↗

DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗