Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transition path”

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 73 records · Page 4

Thermodynamic signatures of an antiferromagnetic quantum critical point inside a superconducting dome

Recent experiments in unconventional superconductors, and in particular iron-based materials, have reported evidence of an antiferromagnetic quantum critical point (AFM-QCP) emerging inside the superconducting dome of the phase diagram. Fluctuations associated with such an AFM-QCP are expected to promote unusual temperature dependencies of thermodynamic quantities. Here, we compute the T dependence of the specific heat C(T) deep inside a fully gapped s +– superconducting state as the AFM-QCP is approached. We find that at the AFM-QCP, the specific heat C(T) vanishes quadratically with temperature, as opposed to the typical exponential suppression seen in fully gapped BCS superconductors. This robust result is due to a nonanalytic contribution to the free energy arising from the general form of the bosonic (AFM) propagator in the SC state. Away from the AFM-QCP, as temperature is lowered, C(T) shows a crossover from a T2 behavior to an exponential behavior, with the crossover temperature scale set by the value of the superconducting gap and the distance to the QCP. As a result, we argue that these features in the specific heat can be used to unambiguously determine the existence of AFM-QCPs inside the superconducting domes of iron-based and other fully gapped unconventional superconductors.

36 MATERIALS SCIENCE↗

Mechanistic study of a CO-free pathway in the methanol oxidation reaction over oxygen vacancies in NiOOH

The methanol oxidation reaction is a key reaction in direct methanol fuel cells. Prior research indicates that if oxygen vacancies in NiOOH serve as the active sites, the methanol oxidation mainly proceeds through the formate-involving pathway, which is a CO-free pathway, distinct from the conventional path over transition metal catalysts, though the fundamental reason for this suppressed CO formation is unclear. Herein, we report density functional theory calculations, through which we uncover the underlying reasons for this alternate path of methanol oxidation over the oxygen vacancies in NiOOH. We find that the existence of oxygen vacancies in NiOOH affects the adsorption configuration of adsorbates and that the interfacial charge transfer is minimal for CHO and CO intermediates. In addition, CHO, a key intermediate to form CO, adsorbs at the oxygen vacancy through the oxygen atom, leading to low stability due to the incomplete valence saturation of the carbon atom. This weak electronic interaction and instability effectively inhibit CHO formation and, consequently, CO formation. These insights provide valuable guidance for the development of efficient and CO-tolerant catalysts for methanol oxidation.

Ho, Thy [University of Oklahoma, Norman, OK (Unite↗

Dynamic strength of iron under pressure-temperature conditions of Earth’s inner core

Iron (Fe) is a primary constituent of terrestrial planetary cores, yet its rheological properties under extreme conditions remain uncertain. Here we present direct measurements of Fe strength at 310-430 GPa pressures and 3700-5800 K temperatures, obtained using Rayleigh-Taylor (RT) instability experiments at the National Ignition Facility. Single-crystal α-Fe samples with [001] and [111] orientations are shock-ramp compressed past the α-ε transition along paths approaching Earth’s inner core conditions. We find that ε-Fe derived from [001] α-Fe is consistently stronger (11-20 GPa) than that from [111] α-Fe (8-18 GPa), contrary to the trend at ambient conditions. Large-scale molecular dynamics simulations reproduce this atypical strength anisotropy and attribute it to microstructural evolution during the phase transition and subsequent ε-phase plasticity. Ripple growth analysis further constrains viscosities of 100-170 Pa·s under the driven conditions. Our results provide experimental benchmarks for Fe rheology at inner-core conditions, with implications for seismic anisotropy and the geodynamo.

Condensed-matter physics↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Redesigning large-scale multimodal transit networks with shared autonomous mobility services

Here, this study addresses a large-scale multimodal transit network design problem, with Shared Autonomous Mobility Services (SAMS) as both transit feeders and an origin-to-destination mode. The framework captures spatial demand and modal characteristics, considers intermodal transfers and express services, determines transit infrastructure investment and path flows, and generates transit routes. A system-optimal multimodal transit network is designed with minimum total door-to-door generalized costs of users and operators, satisfying transit origin-destination demand within a pre-set infrastructure budget. Firstly, the geography, demand, and modes in each zone are characterized with continuous approximation. The decisions of network link investment and multimodal path flows in zonal connection optimization are formulated as a minimum-cost multi-commodity network flow (MCNF) problem and solved efficiently with a mixed-integer linear programming (MILP) solver. Subsequently, the route generation problem is solved by expanding the MCNF formulation to minimize intramodal transfers. The model is illustrated through a set of experiments with the Chicago network comprised of 50 zones and seven modes, under three scenarios. The computational results present savings in traveler journey time and operator cost demonstrating the potential benefits of collaboration between multimodal transit systems and SAMS.

Autonomous vehicles↗

Toward Drilling the Perfect Geothermal Well: An International Research Coordination Network for Geothermal Drilling Optimization Supported by Deep Machine Learning and Cloud Based Data Aggregation

The EDGE project, supported by the U.S. Department of Energy Geothermal Technologies Office under award DE-EE0008793, established a data-driven framework for improving the efficiency, cost-effectiveness, and reliability of geothermal well drilling. The project focused on developing scalable data infrastructure, advanced machine learning and probabilistic models, and integrated analytics tools to support continuous drilling optimization. A central objective was to reduce geothermal drilling costs by up to seventy percent while minimizing the risk of well failure through predictive diagnostics and adaptive planning. Over the project period, a comprehensive data repository was designed and deployed, incorporating records from over one hundred geothermal wells across varied geological settings. This repository supported both structured and unstructured data and adhered to FAIR data principles, enabling provenance tracking, quality control, and standardized metadata. The project introduced automated ingestion pipelines and a cloud-hosted platform that facilitated access to raw, processed, and derived datasets. This infrastructure served as the foundation for model development and analysis. Machine learning workflows were developed to predict key drilling metrics including rate of penetration, non-productive time, and total drilling costs. Self-organizing maps and dimensionality reduction methods were used to uncover operational patterns and outliers, while supervised learning algorithms such as random forests and deep neural networks were applied to forecast performance outcomes. The models were validated on heterogeneous datasets from both U.S. and Icelandic fields, demonstrating variable but significant predictive accuracy. The results indicated that finer temporal resolution, inclusion of lithological data, and consistency in operational annotations could substantially improve model performance. The project also implemented process mining techniques to reconstruct state-transition models from drilling event logs. These models enabled the identification of deviations from optimal workflows and provided insights into recurring failure modes. Analysis of non-productive time highlighted the impact of equipment failures, geological challenges, and human factors, offering opportunities for targeted mitigation strategies. The EDGE Dashboard was developed as a web-based expert system integrating data visualization, model outputs, and user-driven queries. It provided an accessible interface for operators to explore historical data, evaluate predicted outcomes, and compare drilling scenarios. Initial feedback from project partners suggested that the dashboard could serve as a foundation for more advanced advisory and optimization tools. Overall, the EDGE project demonstrated the feasibility and value of applying modern data science techniques to geothermal drilling. It delivered a set of interoperable tools and models that can support more efficient, lower-risk well development. The findings point toward a viable path for transitioning from advisory analytics to semi-autonomous drilling systems, contingent on continued collaboration, expanded datasets, and field validation. The project results have immediate relevance for drilling operations, data management practices, and future geothermal R&D efforts aimed at achieving reliable, cost-competitive geothermal energy at scale.

15 GEOTHERMAL ENERGY↗

Sustained Oxygen Spectral Gaps and Their Dynamic Evolution in the Inner Magnetosphere

Van Allen Probes observations of ion spectra often show a sustained gap within a very narrow energy range throughout the full orbit. To understand their formation mechanism, we statistically investigate the characteristics of the narrow gaps for oxygen ions and find that they are most frequently observed near the noon sector with a peak occurrence rate of over 30%. The magnetic moment (μ) of the oxygen ions in the gap shows a strong dependence on magnetic local time (MLT), with higher and lower μ values in the morning and afternoon sectors, respectively. Moreover, we find through superposed epoch analysis that the gap formation also depends on geomagnetic conditions. Those gaps formed at lower magnetic moments (μ < 3,000 keV/G) are associated with stable convection electric fields, which enable magnetospheric ions to follow a steady drift pattern that facilitates the gap formation by corotational drift resonance. On the other hand, gaps with higher μ values are statistically preceded by a gradual increase of geomagnetic activity. Here, we suggest that ions within the gap were originally located inside the Alfven layer following closed drift paths, before they were transitioned into open drift paths as the convection electric field was enhanced. The sunward drift of these ions, with very low fluxes, forms a drainage void in the dayside magnetosphere manifested as the sustained gap in the oxygen spectrum. This scenario is supported by particle-tracing simulations, which reproduce most of the observed characteristics and therefore provide new insights into inner magnetospheric dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

Designing polar textures with ultrafast neuromorphic features from atomistic simulations

This review summarizes recent works, all using a specific atomistic approach, that predict and explain the occurrence of key features for neuromorphic computing in three archetypical dipolar materials, when they are subject to THz excitations. The main ideas behind such atomistic approach are provided, and illustration of model relaxor ferroelectrics, antiferroelectrics, and normal ferroelectrics are given, highlighting the important potential of polar materials as candidates for neuromorphic computing. Some peculiar emphases are made in this Review, such as the connection between neuromorphic features and percolation theory, local minima in energy path, topological transitions and/or anharmonic oscillator model, depending on the material under investigation. By considering three different and main polar material families, this work provides a complete and innovative toolbox for designing polar-based neuromorphic systems.

THz pulses↗

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the nonlinear stability of a high-speed, axisymmetric boundary layer

The stability of a high-speed, axisymmetric boundary layer is investigated using secondary instability theory and direct numerical simulation. Parametric studies based on the temporal secondary instability theory identify subharmonic secondary instability as a likely path to transition on a cylinder at Mach 4.5. The theoretical predictions are validated by direct numerical simulation at temporally-evolving primary and secondary disturbances in an axisymmetric boundary-layer flow. At small amplitudes of the secondary disturbance, predicted growth rates agree to several significant digits with values obtained from the spectrally-accurate solution of the compressible Navier-Stokes equations. Qualitative agreement persists to large amplitudes of the secondary disturbance. Moderate transverse curvature is shown to significantly affect the growth rate of axisymmetric second mode disturbances, the likely candidates of primary instability. The influence of curvature on secondary instability is largely indirect but most probably significant, through modulation of the primary disturbance amplitude. Subharmonic secondary instability is shown to be predominantly inviscid in nature, and to account for spikes in the Reynolds stress components at or near the critical layer.

Pruett, C. David↗

On the nonlinear stability of a high-speed, axisymmetric boundary layer

The stability of a high-speed, axisymmetric boundary layer is investigated using secondary instability theory and direct numerical simulation. Parametric studies based on the temporal secondary instability theory identify subharmonic secondary instability as a likely path to transition on a cylinder at Mach 4.5. The theoretical predictions are validated by direct numerical simulation at temporally-evolving primary and secondary disturbances in an axisymmetric boundary-layer flow. At small amplitudes of the secondary disturbance, predicted growth rates agree to several significant digits with values obtained from the spectrally-accurate solution of the compressible Navier-Stokes equations. Qualitative agreement persists to large amplitudes of the secondary disturbance. Moderate transverse curvature is shown to significantly affect the growth rate of axisymmetric second mode disturbances, the likely candidates of primary instability. The influence of curvature on secondary instability is largely indirect but most probably significant, through modulation of the primary disturbance amplitude. Subharamonic secondary instability is shown to be predominantly inviscid in nature, and to account for spikes in the Reynolds stress components at or near the critical layer.

Pruett, C. D.↗

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

Temporal direct numerical simulation of laminar breakdown via subharmonic secondary instability in high-speed axisymmetric boundary layers has been accomplished using a highly accurate, fully explicit algorithm which combines spectral collocation and high-order compact-difference techniques. Numerical test cases confirm that subharmonic secondary instability is confirmed to be a viable path to transition in high-speed boundary-layer flow. Secondary instability is shown to account for peaks in the Reynolds stresses at or near the critical layer which are not possible from the second-mode primary instability alone. Reynolds stresses spatially reconstructed from the temporal model via the Gaster transformation show a 'spreading angle' of about 12 deg, in qualitative agreement with experimental findings. The rate of broadening of the Reynolds stress peak is a strongly nonlinear phenomenon which cannot be reproduced by secondary instability theory.

Pruett, C. D.↗

Nonlinear Evolution and Breakdown of Azimuthally Compact Crossflow Vortex Pattern over a Yawed Cone

Hypersonic boundary-layer flows over a circular cone at moderate angle of incidence can support strong crossflow instability in between the windward and leeward rays on the plane of symmetry. Due to the more efficient excitation of stationary crossflow vortices by surface roughness, a possible path to transition in such flows corresponds to rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. In the present paper, the previous analyses of crossflow instability over a 7- degree half-angle, yawed circular cone in a Mach 6 free stream have been extended to the nonlinear evolution of azimuthally localized crossflow vortex packets and the amplification characteristics and nonlinear breakdown of high-frequency secondary instabilities associated with those packets. A comparison between plane marching PSE and direct Navier-Stokes simulations (DNS) reveals favorable agreement in regard to mode shapes, most amplified disturbance frequencies, and N-factor evolution. In contrast, the quasi-parallel predictions are found to result in severe underprediction of the N-factors. The direct numerical simulations also indicate that the breakdown of secondary instabilities in a 3D hypersonic boundary layer shares certain common features with the previous computations of crossflow transition over subsonic swept wings.

Vortex structure↗

Voltage decay and redox asymmetry mitigation by reversible cation migration in lithium-rich layered oxide electrodes

Despite the high energy density of lithium-rich layered-oxide electrodes, their real-world implementation in batteries is hindered by the substantial voltage decay on cycling. This voltage decay is widely accepted to mainly originate from progressive structural rearrangements involving irreversible transition-metal migration. As prevention of this spontaneous cation migration has proven difficult, a paradigm shift toward management of its reversibility is needed. Herein, we demonstrate that the reversibility of the cation migration of lithium-rich nickel manganese oxides can be remarkably improved by altering the oxygen stacking sequences in the layered structure and thereby dramatically reducing the voltage decay. The preeminent intra-cycle reversibility of the cation migration is experimentally visualized, and first-principles calculations reveal that an O2-type structure restricts the movements of transition metals within the Li layer, which effectively streamlines the returning migration path of the transition metals. Furthermore, we propose that the enhanced reversibility mitigates the asymmetry of the anionic redox in conventional lithium-rich electrodes, promoting the high-potential anionic reduction, thereby reducing the subsequent voltage hysteresis. Our findings demonstrate that regulating the reversibility of the cation migration is a practical strategy to reduce voltage decay and hysteresis in lithium-rich layered materials.

25 ENERGY STORAGE↗

On the structure of cellular solutions in Rayleigh-Benard-Marangoni flows in small-aspect-ratio containers

Multiple steady flow patterns occur in surface-tension/buoyancy-driven convection in a liquid layer heated from below (Rayleigh-Benard-Marangoni flows). Techniques of numerical bifurcation theory are used to study the multiplicity and stability of two-dimensional steady flow patterns (rolls) in rectangular small-aspect-ratio containers as the aspect ratio is varied. For pure Marangoni flows at moderate Biot and Prandtl number, the transitions occurring when paths of codimension 1 singularities intersect determine to a large extent the multiplicity of stable patterns. These transitions also lead, for example, to Hopf bifurcations and stable periodic flows for a small range in aspect ratio. The influence of the type of lateral walls on the multiplicity of steady states is considered. 'No-slip' lateral walls lead to hysteresis effects and typically restrict the number of stable flow patterns (with respect to 'slippery' sidewalls) through the occurrence of saddle node bifurcations. In this way 'no-slip' sidewalls induce a selection of certain patterns, which typically have the largest Nusselt number, through secondary bifurcation.

Dijkstra, Henk A.↗

An Alternative Humans to Mars Approach: Reducing Mission Mass with Multiple Mars Flyby Trajectories and Minimal Capability Investments

Mars flyby trajectories and Earth return trajectories have the potential to enable lower- cost and sustainable human exploration of Mars. Flyby and return trajectories are true minimum energy paths with low to zero post-Earth departure maneuvers. By emplacing the large crew vehicles required for human transit on these paths, the total fuel cost can be reduced. The traditional full-up repeating Earth-Mars-Earth cycler concept requires significant infrastructure, but a Mars only flyby approach minimizes mission mass and maximizes opportunities to build-up missions in a stepwise manner. In this paper multiple strategies for sending a crew of 4 to Mars orbit and back are examined. With pre-emplaced assets in Mars orbit, a transit habitat and a minimally functional Mars taxi, a complete Mars mission can be accomplished in 3 SLS launches and 2 Mars Flyby's, including Orion. While some years are better than others, ample opportunities exist within a given 15-year Earth-Mars alignment cycle. Building up a mission cadence over time, this approach can translate to Mars surface access. Risk reduction, which is always a concern for human missions, is mitigated by the use of flybys with Earth return (some of which are true free returns) capability.

Whitley, Ryan J.↗

Tracing structural optima as a function of available resources by a homotopy method

Optimization problems are typically solved by starting with an initial estimate and proceeding iteratively to improve it until the optimum is found. The design points along the path from the initial estimate to the optimum are usually of no value. The present work proposes a strategy for tracing a path of optimum solutions parameterized by the amount of available resources. The paper specifically treats the optimum design of a structure to maximize its buckling load. Equations for the optimum path are obtained using Lagrange multipliers, and solved by a homotopy method. The solution path has several transitions from unimodal to bimodal solutions. The Lagrange multipliers and second-order optimality conditions are used to detect branching points and to switch to the optimum solution path. The procedure is applied to the design of a foundation which supports a column for maximum buckling load. Using the total available foundation stiffness as a homotopy parameter, a set of optimum foundation designs is obtained.

Haftka, Raphael T.↗