Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Genesis of Pre-Hurricane Felix (2007): The Role of the Easterly Wave Critical Layer - Part 1

The formation of pre Hurricane Felix (2007) in a tropical easterly wave is examined in a two-part study using the Weather Research and Forecasting (WRF) model with a high-resolution nested grid configuration that permits the representation of cloud system processes. The simulation commences during the wave stage of the precursor African easterly-wave disturbance. Here the simulated and observed developments are compared, while in Part II of the study various large-scale analyses, physical parameterizations, and initialization times are explored to document model sensitivities. In this first part the authors focus on the wave/vortex morphology, its interaction with the adjacent intertropical convergence zone complex, and the vorticity balance in the neighborhood of the developing storm. Analysis of the model simulation points to a bottom-up development process within the wave critical layer and supports the three new hypotheses of tropical cyclone formation proposed recently by Dunkerton, Montgomery, and Wang. It is shown also that low-level convergence associated with the ITCZ helps to enhance the wave signal and extend the "wave pouch" from the jet level to the top of the atmospheric boundary layer. The region of a quasi-closed Lagrangian circulation within the wave pouch provides a focal point for diabatic merger of convective vortices and their vortical remnants. The wave pouch serves also to protect the moist air inside from dry air intrusion, providing a favorable environment for sustained deep convection. Consistent with the authors' earlier findings, the tropical storm forms near the center of the wave pouch via system-scale convergence in the lower troposphere and vorticity aggregation. Components of the vorticity balance are shown to be scale dependent, with the immediate effects of cloud processes confined more closely to the storm center than the overturning Eliassen circulation induced by diabatic heating, the influence of which extends to larger radii.

Wang, Zhuo↗

Remote Sensing Product Verification and Validation at the NASA Stennis Space Center

Remote sensing data product verification and validation (V&V) is critical to successful science research and applications development. People who use remote sensing products to make policy, economic, or scientific decisions require confidence in and an understanding of the products' characteristics to make informed decisions about the products' use. NASA data products of coarse to moderate spatial resolution are validated by NASA science teams. NASA's Stennis Space Center (SSC) serves as the science validation team lead for validating commercial data products of moderate to high spatial resolution. At SSC, the Applications Research Toolbox simulates sensors and targets, and the Instrument Validation Laboratory validates critical sensors. The SSC V&V Site consists of radiometric tarps, a network of ground control points, a water surface temperature sensor, an atmospheric measurement system, painted concrete radial target and edge targets, and other instrumentation. NASA's Applied Sciences Directorate participates in the Joint Agency Commercial Imagery Evaluation (JACIE) team formed by NASA, the U.S. Geological Survey, and the National Geospatial-Intelligence Agency to characterize commercial systems and imagery.

Stanley, Thomas M.↗

Large-scale simulations of Floquet physics on near-term quantum computers

Abstract Periodically driven quantum systems exhibit a diverse set of phenomena but are more challenging to simulate than their equilibrium counterparts. Here, we introduce the Quantum High-Frequency Floquet Simulation (QHiFFS) algorithm as a method to simulate fast-driven quantum systems on quantum hardware. Central to QHiFFS is the concept of a kick operator which transforms the system into a basis where the dynamics is governed by a time-independent effective Hamiltonian. This allows prior methods for time-independent simulation to be lifted to simulate Floquet systems. We use the periodically driven biaxial next-nearest neighbor Ising (BNNNI) model, a natural test bed for quantum frustrated magnetism and criticality, as a case study to illustrate our algorithm. We implemented a 20-qubit simulation of the driven two-dimensional BNNNI model on Quantinuum’s trapped ion quantum computer. Our error analysis shows that QHiFFS exhibits not only a cubic advantage in driving frequency ω but also a linear advantage in simulation time t compared to Trotterization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validation of highly reliable, real-time knowledge-based systems

Knowledge-based systems have the potential to greatly increase the capabilities of future aircraft and spacecraft and to significantly reduce support manpower needed for the space station and other space missions. However, a credible validation methodology must be developed before knowledge-based systems can be used for life- or mission-critical applications. Experience with conventional software has shown that the use of good software engineering techniques and static analysis tools can greatly reduce the time needed for testing and simulation of a system. Since exhaustive testing is infeasible, reliability must be built into the software during the design and implementation phases. Unfortunately, many of the software engineering techniques and tools used for conventional software are of little use in the development of knowledge-based systems. Therefore, research at Langley is focused on developing a set of guidelines, methods, and prototype validation tools for building highly reliable, knowledge-based systems. The use of a comprehensive methodology for building highly reliable, knowledge-based systems should significantly decrease the time needed for testing and simulation. A proven record of delivering reliable systems at the beginning of the highly visible testing and simulation phases is crucial to the acceptance of knowledge-based systems in critical applications.

Johnson, Sally C.↗

Modeling of turbulent separated flows for aerodynamic applications

Steady, high speed, compressible separated flows modeled through numerical simulations resulting from solutions of the mass-averaged Navier-Stokes equations are reviewed. Emphasis is placed on benchmark flows that represent simplified (but realistic) aerodynamic phenomena. These include impinging shock waves, compression corners, glancing shock waves, trailing edge regions, and supersonic high angle of attack flows. A critical assessment of modeling capabilities is provided by comparing the numerical simulations with experiment. The importance of combining experiment, numerical algorithm, grid, and turbulence model to effectively develop this potentially powerful simulation technique is stressed.

Marvin, J. G.↗

Modeling of turbulent separated flows for aerodynamic applications

Steady, high speed, compressible separated flows modeled through numerical simulations resulting from solutions of the mass-averaged Navier-Stokes equations are reviewed. Emphasis is placed on benchmark flows that represent simplified (but realistic) aerodynamic phenomena. These include impinging shock waves, compression corners, glancing shock waves, trailing edge regions, and supersonic high angle of attack flows. A critical assessment of modeling capabilities is provided by comparing the numerical simulations with experiment. The importance of combining experiment, numerical algorithm, grid, and turbulence model to effectively develop this potentially powerful simulation technique is stressed.

Marvin, J. G.↗

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

Friction stir welding (FSW) has found increased applications in automotive and aerospace industries due to its advantages of solid-state bonding, no fusion and melting, and versatility in various working conditions and material combinations. However, the relationship among processing parameters, material properties, and bonding extent and fidelity remains largely empirical, primarily because of the lack of the mechanistic understanding of the tool-workpiece frictional behavior that affects our subsequent understanding of microstructural evolution and interface bonding formation. While the tool-workpiece stick-slip condition is believed to dictate the resulting torque and heat generation rate during the welding process, it remains rare and elusive to conduct a quantitative experimental measurement of such interfacial field. On the other hand, numerical simulations based on Computational Fluid Dynamics (CFD) rely on ad hoc assumptions of interfacial pressure and shear-stress conditions, but predictions can only be validated via the medium- and far-range temperature field which is known to be insensitive to the interfacial frictional behavior. This work first presents a comparison among two CFD-based simulation methodologies and the Coupled Eulerian Lagrangian (CEL) model in finite element method, the last of which uses the Coulomb friction so that the stick-slip is naturally developed. Based on the Hill-Bower similarity relationship in the contact analysis, an analytical model is developed here to prove why a constant stick-slip fraction will be developed in the steady state, to correlate the stick-slip fraction to processing parameters such as the tool spin rate, and further to derive dimensionless functions for torque and heat-generation-rate predictions. Pros and cons of various numerical approaches in predicting stick-slip are discussed, and our analytical model has been found to agree well with our numerical simulation and literature experimental results. These analyses provide the critical strain-rate and temperature fields that are needed for the bonding analysis in our future work.

36 MATERIALS SCIENCE↗

On the origin of internal obstacles to dislocation glide in single-phase NiFe random alloys

Over the past decade, a new class of unconventional alloys, where several constitutive elements are presented in close to equal concentrations, is considered as promising structural materials in many applications. Properties and behavior of many concentrated solid solution alloys (CSSAs) under different conditions, including mechanical loading, differs from that of conventional alloys and has been the subject of intensive studies by different techniques. A key feature of these materials is that there is no clear distinction between solute and solvent species of atoms, which is the conventional approach to considering dislocation motion in alloys. To understand this regime, we have recently reported on the glide of a screw dislocation glide in a simple equiatomic NiFe alloy and reported the two-mode character. At low applied stress, the dislocation motion is similar to a glide through a field of obstacles. In this work, we report the results of a detailed study on the nature of these internal obstacles. We demonstrate that these are not conventional localized “obstacles” but is the result of significantly different energies of instantaneous dislocation configurations over the alloy bulk. The strong pinning positions are associated with composition changes in the partial dislocation cores, while observing no significant composition changes along the stacking faults. Local energies can vary significantly, as demonstrated using simulations of short-periodicity dislocation lines. Therefore, the critical stress to move a dislocation under such conditions is not a constant Peierls stress but represents a wide spectrum of barriers between different non-straight valleys that the dislocation line accommodates in the bulk at particular stress/strain conditions.

36 MATERIALS SCIENCE↗

Influence of Ag metal dispersion on the thermal conversion of ethanol to butadiene over Ag-ZrO 2 /SiO 2 catalysts

Atomistic scale models were developed and coupled with experimental investigation to deliver a functional understanding of catalytic activity and selectivity in the conversion of ethanol to 1,3-butadiene over Ag/ZrO 2 /SiO 2 . A detailed evaluation of the structural and electronic properties of the resultant catalyst models led to the identification of critical active sites of the catalyst. More importantly, the extent of Ag dispersion on the SiO 2 support and relative proximity to ZrO 2 were found to vary with the oxidation state of Ag and local coordination environment (Ag-O SiO2 ) allowing for critical control of ethanol conversion towards butadiene or ethylene. Simulations revealed that less dispersed or clustered Ag contain predominantly Ag 0 charge state and promote conversion of ethanol to ethylene. The well-dispersed Ag/ZrO 2 /SiO 2 catalyst instead contain a larger fraction of cationic Ag δ+ , and predominantly promote ethanol dehydrogenation and subsequent production of butadiene. The theoretical insights drawn were validated and confirmed experimentally using TEM, XRD and reactivity measurements demonstrating the effect of Ag dispersion on the selectivity of ethanol conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predictions of Boron Phase Stability Using an Efficient Bayesian Machine Learning Interatomic Potential

Thermodynamic phase stability of three elemental boron allotropes, i.e., α-B, β-B, and γ-B, was investigated using a Bayesian interatomic potential trained via a sparse Gaussian process (SGP). SGP potentials trained with datasets from on-the-fly active learning achieve quantum mechanical level accuracy when employed in molecular dynamics simulations to predict wide-ranging thermodynamic, structural, and vibrational properties. The simulated phase diagram (500~1400 K and 0~16 GPa) agrees with experimental measurements. The SGP-based MD simulations also successfully predicted that the B13 defect is critical in stabilizing β-B below 700 K. At higher temperatures, the entropy becomes the dominant factor, making β-B the more stable phase over α-B. Furthermore, this Letter demonstrates that SGP potentials based on a training set consisting of defect-free-only systems could make correct predictions of defect-related phenomena in solid-state crystals, paving the path to investigate crystal phase stability and transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗

Crossover from Rouse to Reptation Dynamics in Salt-Free Polyelectrolyte Complex Coacervates

Considerable interest in the dynamics and rheology of polyelectrolyte complex coacervates has been motivated by their industrial application as viscosity modifiers. A central question is the extent to which classical Rouse and reptation models can be applied to systems where electrostatic interactions play a critical role on the thermodynamics. By relying on molecular simulations, we present a direct analysis of the crossover from Rouse to reptation dynamics in salt-free complex coacervates as a function of chain length. This crossover shifts to shorter chain lengths as electrostatic interactions become stronger, which corresponds to the formation of denser coacervates. To distinguish the roles of Coulomb interactions and density, here we compare the dynamics of coacervates to those of neutral, semidilute solutions at the same density. Both systems exhibit a universal dynamical behavior in the connectivity-dominated (subdiffusion and normal diffusion) regimes, but the monomer relaxation time in coacervates is much longer and increases with increasing Bjerrum length. This is similar to the cage effect observed in glass-forming polymers, but the local dynamical slowdown is caused here by strong Coulomb attractions (ion pairing) between oppositely charged monomers. Our findings provide a microscopic framework for the quantitative understanding of coacervate dynamics and rheology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wedge-Shaped Hopper Design for Milled Woody Biomass Flow

The bioenergy industry has been challenged by unstable flow and transport of milled biomass in material handling operations. Handling issues such as hopper clogging and auger jamming are attributed to knowledge gaps between existing handling units designed for bulk solids and their suitability for milled biomass with high compressibility. This work investigates various flow behaviors of granular woody biomass in wedge-shaped hoppers. Herein, hopper flow physical experiments and numerical simulations are conducted to study the influence of the critical material attributes and critical processing parameters on the flow pattern, arching, and throughput. The results show that (1) the preferred flow pattern, mass flow, can be achieved by controlling the material’s internal friction angle, hopper inclination, and hopper wall friction; (2) hopper arching, governed by the competing gravity-driven force against flow resistance from material internal friction and material–wall friction, can be controlled by the hopper wall friction angle and the inclination angle; and (3) flow throughput can be accurately estimated from our empirical equation with inputs of hopper outlet geometry and particle-scale to bulk-scale material attributes. This study elucidates woody biomass flow physics and provides guidance for industrial equipment design.

09 BIOMASS FUELS↗

Crystalline Peptoid Nanofibers with a Single-Unit Cell Cross Section

Ultranarrow crystalline one-dimensional nanostructures formed from soft materials facilitate precise structural control in nanomaterial design, which is essential for biomedicine and nanotechnology applications. Systematic control of their hierarchical structure is challenging due to the complexities of simultaneously manipulating multiple noncovalent interactions at such small scales. We employed a polypeptoid crystal motif as a supramolecular synthon to engineer ultranarrow crystalline nanofibers constrained to a single lattice axis by incorporating a single ionizable side chain into the hydrophobic core of a nanosheet-forming peptoid. Cryogenic transmission electron microscopy of the nanofibers revealed detailed molecular arrangements of a unit-cell cross-section and the presence of distinct pH-dependent lattice isoforms that resulted in morphological transformations. Molecular dynamics simulations demonstrated that the ionizable side chain plays a critical role in changing the local conformation of the unit cell, which further impacts the dimensionality of hierarchical structures. Moreover, these fibers were readily functionalized with biological ligands to afford one-dimensional (1D) protein arrays. This approach for the high-precision bottom-up assembly of ultranarrow 1D nanostructures offers significant potential for developing novel biomimetic nanostructures.

Lee, Yen Jea↗

Giant Cloud Condensation Nuclei Facilitate Drizzle Formation in Stratocumulus—Insights From a Combined Observation‐Modeling Framework

The mechanism for initiating drizzle drop remains a gap in the current understanding of warm rain formation. One prevalent hypothesis suggests that the presence of Giant Cloud Condensation Nuclei (GCCN) generates drizzle‐sized drops necessary to trigger the Collision‐Coalescence (C‐C) process. Here, in this study, this hypothesis is investigated using a novel framework that integrates in situ observations, remote sensing measurements, and idealized models. Results show that GCCN can efficiently generate drizzle drops through condensation, producing a broad Droplet Size Distribution (DSD) comparable to in situ observations. The large drizzle drop and broad DSD strongly facilitate C‐C, further accelerating drizzle initiation. To compare with observation, the model‐generated DSDs are used to generate radar Doppler spectra where radar reflectivity and Doppler skewness is estimated. The simulated radar quantities correspond well with radar observations, providing critical evidence for the GCCN‐induced drizzle initiation mechanism.

54 ENVIRONMENTAL SCIENCES↗

Probing three-dimensional mesoscopic interfacial structures in a single view using multibeam X-ray coherent surface scattering and holography imaging

Abstract Visualizing surface-supported and buried planar mesoscale structures, such as nanoelectronics, ultrathin-film quantum dots, photovoltaics, and heterogeneous catalysts, often requires high-resolution X-ray imaging and scattering. Here, we discovered that multibeam scattering in grazing-incident reflection geometry is sensitive to three-dimensional (3D) structures in a single view, which is difficult in conventional scattering or imaging approaches. We developed a 3D finite-element-based multibeam-scattering analysis to decode the heterogeneous electric-field distribution and to faithfully reproduce the complex scattering and surface features. This approach further leads to the demonstration of hard-X-ray Lloyd’s mirror interference of scattering waves, resembling dark-field, high-contrast surface holography under the grazing-angle scattering conditions. A first-principles calculation of the single-view holographic images resolves the surface patterns’ 3D morphology with nanometer resolutions, which is critical for ultrafine nanocircuit metrology. The holographic method and simulations pave the way for single-shot structural characterization for visualizing irreversible and morphology-transforming physical and chemical processes in situ or operando .

47 OTHER INSTRUMENTATION↗

Bottom-up synthesis of mesoscale nanomeshes of graphene nanoribbons on germanium

The synthesis of functional graphene nanostructures on Ge(001) provides an attractive route toward integrating graphene-based electronic devices onto complementary metal oxide semiconductor-compatible platforms. In this study, we leverage the phenomenon of the anisotropic growth of graphene nanoribbons from rationally placed graphene nanoseeds and their rotational self-alignment during chemical vapor deposition to synthesize mesoscale graphene nanomeshes over areas spanning several hundred square micrometers. Lithographically patterned nanoseeds are defined on a Ge(001) surface at pitches ranging from 50 to 100 nm, which serve as starting sites for subsequent nanoribbon growth. Rotational self-alignment of the nanoseeds followed by anisotropic growth kinetics causes the resulting nanoribbons to be oriented along each of the equivalent, orthogonal Ge < 110 > directions with equal probability. As the nanoribbons grow, they fuse, creating a continuous nanomesh. In contrast to nanomesh synthesis via top-down approaches, this technique yields nanomeshes with atomically faceted edges and covalently bonded junctions, which are important for maximizing charge transport properties. Additionally, we simulate the electrical characteristics of nanomeshes synthesized from different initial nanoseed-sizes, size-polydispersities, pitches, and device channel lengths to identify a parameter-space for acceptable on/off ratios and on-conductance in semiconductor electronics. The simulations show that decreasing seed diameter and pitch are critical to increasing nanomesh on/off ratio and on-conductance, respectively. With further refinements in lithography, nanomeshes obtained via seeded synthesis and anisotropic growth are likely to have superior electronic properties with tremendous potential in a multitude of applications, such as radio frequency communications, sensing, thin-film electronics, and plasmonics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗