Search NASA⌕ Search

Engineering topics

Ihme, Matthias

Publications and source records attributed to Ihme, Matthias.

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the importance of species immiscibility in mixing-layer dynamics at supercritical pressures

The mixing dynamics of injected propellants is a key factor in determining the ignition performance and combustion-instability response of rocket engines, internal combustion engines, and gas turbines. A key source of uncertainty in the prediction of phase-exchange dynamics is the immiscibility of the injected propellants at supercritical pressure conditions, which induces liquid/vapor phase separation and surface-tension dynamics. While experimental observations indicate the presence of liquid/vapor interfacial structures, this interfacial dynamics is typically neglected in numerical analyses. To address this issue, the objective of the present study is to systematically evaluate the importance of species immiscibility on the phase-exchange dynamics of cryogenic LOX/GH2 mixing layers at typical rocket engine injection conditions. This is accomplished by comparing simulations of (i) a recently developed interface-capturing Regularized-Interface Method (RIM) formulation, and (ii) the commonly employed Diffuse-Interface Method (DIM) formulation. Analysis shows that the interfacial dynamics significantly impact the atomization and mixing of the propellants in the near-injector region, which is not captured by the DIM formulation. The findings of this study extend to other immiscible injection systems, such as LOX/kerosene and hydrocarbon-fuel/air, thereby highlighting the importance of resolving species immiscibility in simulating high-pressure combustion engines.

diffused-interface method↗
DOE OSTI2024

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗
DOE OSTI2024

Final Technical Report: Center for Mechanistic Control of Unconventional Formations

The overarching mission of the Center for Mechanistic Control of Unconventional Formations (CMC-UF) was to garner cross-cutting, fundamental, geoscience knowledge to achieve mechanistic control over the strongly coupled nonequilibrium physical and geochemical processes in extreme geological environments including shale, mudstone, marls, and other tight rocks with nanoscale pores. Collectively, these are referred to as unconventional formations and they often play the role of seals for other subsurface storage formations. The fundamental knowledge garnered by CMC-UF enabled science-based management of US unconventional formations for subsurface storage of carbon dioxide and TWh quantities of renewable energy as hydrogen and/or compressed air over longer timeframes as well as for natural gas production, with reduced environmental impacts, in the short term.

58 GEOSCIENCES↗
DOE OSTI2024

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗
DOE OSTI2024

Probing Supercritical Phase Transition using Ultrafast X-ray Diagnostics (Final Technical Report)

Over recent years, several hypotheses have been put forward to elucidate supercritical transition states. Despite significant progress on the fundamental understanding of fluids at these supercritical conditions, important questions concerning the microstructure and dynamical processes remain. A particular research issue is hereby the fundamental understanding of the morphology of the molecular microstructure and its effect on the macroscopic behaviour and thermodynamic response functions. It is widely believed that the supercritical state is homogeneous without structural and dynamic observables to distinguish between a liquid and a vapor. However, recent investigations have identified regions of distinct liquid-like or vapor-like properties even under supercritical conditions. Specifically, it was shown that the transition between liquid-like and vapor-like states occurs across an extension to the coexistence line, marked by almost discontinuous changes in fluid properties. This transition was first identified experimentally by Nishikawa and Tanaka and stands in contrast to the classical presentation of the supercritical state space as a featureless, homogeneous domain. Structurally, the most important properties of supercritical fluids are the dynamic heterogeneities and local density fluctuations that are present within the fluid state at the microscopic level. These inhomogeneities are associated with the formation of molecular clusters of various sizes with liquid-like properties separated by voids of unbound gas-phase molecules which continuously restructure itself at picosecond timescales. By utilizing ultrafast X-ray Photon Correlation Spectroscopy (XPCS) at the Linac Coherent Light Source (LCLS) and X-ray synchrotron measurements as an experimental probe of these structural changes at the molecular level we seek to elucidate the higher-order phase transition from a liquid-like to a vapor-like state in the region extending the critical point. Open questions we seek to address particularly with XPCS are: (i) the underlying mechanisms responsible for the density fluctuations at supercritical conditions and how these dynamical processes at the molecular level affect the thermodynamic response functions and (ii) on what timescales do dynamical processes evolve that are associated with thermalization, molecular diffusion, and intermolecular cluster transfer.

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

Development of High-Fidelity and Efficient Modeling Capabilities for Enabling Co-Optimization of Fuels and Multi-Mode Engines

The objective of this closely coordinated research project is to develop improved physical sub-models and innovative numerical algorithms to enable accurate and efficient predictions of multi-mode combustion to support EERE’s Co-Optima program. To this end, a research pro-gram is proposed that seeks • to develop accurate submodels for predicting multi-mode combustion regimes, wall-heat transfer, non-equilibrium plasma ignition, and combustion-mode transition, • to develop dynamic load-rebalancing algorithms and efficient semi-implicit time-integration schemes for enabling exascale computing, • to validate the computational submodels against experiments and direct numerical simulation (DNS) data in conjunction with sensitivity analysis, and • to support the direct transition of the resulting sub-models into other engine-simulation tools through the development of self-contained modules. All sub-models will be integrated into the open-source exascale engine-simulation platform Nek5000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

A versatile pressure-cell design for studying ultrafast molecular-dynamics in supercritical fluids using coherent multi-pulse x-ray scattering

Supercritical fluids (SCFs) can be found in a variety of environmental and industrial processes. They exhibit an anomalous thermodynamic behavior, which originates from their fluctuating heterogeneous micro-structure. Characterizing the dynamics of these fluids at high temperature and high pressure with nanometer spatial and picosecond temporal resolution has been very challenging. The advent of hard x-ray free electron lasers has enabled the development of novel multi-pulse ultrafast x-ray scattering techniques, such as x-ray photon correlation spectroscopy (XPCS) and x-ray pump x-ray probe (XPXP). These techniques offer new opportunities for resolving the ultrafast microscopic behavior in SCFs at unprecedented spatiotemporal resolution, unraveling the dynamics of their micro-structure. However, harnessing these capabilities requires a bespoke high-pressure and high-temperature sample system that is optimized to maximize signal intensity and address instrument-specific challenges, such as drift in beamline components, x-ray scattering background, and multi-x-ray-beam overlap. We present a pressure cell compatible with a wide range of SCFs with built-in optical access for XPCS and XPXP and discuss critical aspects of the pressure cell design, with a particular focus on the design optimization for XPCS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Analysis of direct and indirect noise in a next-generation aviation gas turbine combustor

A large-eddy simulation (LES) of a next-generation combustor is performed to examine effects of this combustor concepts on direct and indirect combustion noise characteristics. Direct noise is computed considering the unsteady heat release predictions while entropy fluctuations in the downstream part of the combustor are used to estimate indirect noise through the transfer function of the outlet nozzle. A low-order acoustic reconstruction technique, which utilizes the Green’s function of the configuration, is developed to compute the acoustics inside the combustor. Comparisons to experimental results are reported, showing the reasonable accuracy of the LES and combustion noise computations. The direct noise spectrum peaks at frequencies around 3 kHz because of the fast timescales associated with the chemical reactions inside the combustor. In addition, the indirect noise is dominant at low frequencies, i.e., less than 400 Hz, because of the characteristics of the entropy spectrum at the outlet nozzle. Further, compared to a conventional rich-quench-lean (RQL) combustor configuration, significant differences in the acoustics are observed, and are explained by two specific technological novelties. Firstly, the direct noise spectrum peaks at significantly higher frequencies due to the combustor’s compactness. Secondly, the entropy inhomogeneties downstream of the combustion region are an order of magnitude smaller in amplitude due to the lack of dilution with cold air. This leads to a significant reduction in the indirect noise amplitude at low frequencies, where it is the dominant source of noise. These results suggest opportunities for advanced combustion technologies to achieve substantial reductions in combustion-related noise.

42 ENGINEERING↗
DOE OSTI2023

Evaluation of Electron Tomography Capabilities for Shale Imaging

Despite the advantageous resolution of electron tomography (ET), reconstruction of three-dimensional (3D) images from multiple two-dimensional (2D) projections presents several challenges, including small signal-to-noise ratios, and a limited projection range. This study evaluates the capabilities of ET for thin sections of shale, a complex nanoporous medium. A numerical phantom with 1.24 nm pixel size is constructed based on the tomographic reconstruction of a Barnett shale. A dataset of 2D projection images is numerically generated from the 3D phantom and studied over a range of conditions. First, common reconstruction techniques are used to reconstruct the shale structure. The reconstruction uncertainty is quantified by comparing overall values of storage and transport metrics, as well as the misclassification of pore voxels compared to the phantom. We then select the most robust reconstruction technique and we vary the acquisition conditions to quantify the effect of artifacts. We find a strong agreement for large pores over the different acquisition workflows, while a wider variability exists for nanometer-scale features. Furthermore, the limited projection range and reconstruction are identified as the main experimental bottlenecks, thereby suggesting that sample thinning, advanced holders, and advanced reconstruction algorithms offer opportunities for improvement.

3D scanning transmission electron microscopy↗
DOE OSTI2023

A high-resolution large-eddy simulation framework for wildland fire predictions using TensorFlow

Background: Wildfires are becoming more severe, so we need improved tools to predict them over a wide range of conditions and scales. One approach towards this goal entails the use of coupled fire/atmosphere modelling tools. Although significant progress has been made in advancing their physical fidelity, existing tools have not taken full advantage of emerging programming paradigms and computing architectures to enable high-resolution wildfire simulations. Aims: The aim of this study was to present a new framework that enables landscape-scale wildfire simulations with physical representation of combustion at an affordable cost. Methods: We developed a coupled fire/atmosphere simulation framework using TensorFlow, which enables efficient and scalable computations on Tensor Processing Units. Key Results: Simulation results for a prescribed fire were compared with experimental data. Predicted fire behavior and statistical analysis for fire spread rate, scar area, and intermittency showed overall reasonable agreement. Scalability analysis was performed, showing close to linear scaling. Conclusions: While mesh refinement was shown to have less impact on global quantities, such as fire scar area and spread rate, it benefits predictions of intermittent fire behavior, buoyancy-driven dynamics, and small-scale turbulent motion. Implications: This new simulation framework is efficient in capturing both global quantities and unsteady dynamics of wildfires at high spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Local Rearrangement in Adsorption Layers of Nanoconfined Ethane

Fluids under nanoconfinement control many environmental and engineering processes. In particular, nanoconfined ethane commonly occurs in hydrocarbon production and gas separation applications. While previous research has measured the pore critical properties and self-diffusivity of nanoconfined ethane, the relationship between molecular organization and phase behavior in nanoporous media persists as a knowledge gap. Through experimental adsorption measurements and molecular dynamics simulations, we report the impact of pore wall proximity, pore size, and temperature on the structure of nanoconfined ethane. Finally, we offer evidence for an energetically favorable rearrangement of ethane molecules in the adsorption layer that may contribute to more accurate gas-in-place estimates and more efficient design of gas separation membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2020

Modeling Heatshield Erosion due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant↗
NASA NTRS2020

Modeling Heatshield Erosion Due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant↗
NASA NTRS2019

Numerical Predictions of Dust-Induced Heat Flux Augmentation in Hypersonic Blunt-Body Flows Using a Discontinuous Galerkin Multiphase Flow Solver

Recent interest in human-scale missions to Mars has motivated the need for high-fidelity simulations of reentry flows. During a dust storm, there can be high levels of suspended dust in the Martian atmosphere, which cannot only enhance erosion of thermal protection systems but also transfer energy and momentum to the shock layer, thereby significantly augmenting the surface heat flux. Second-order finite-volume schemes are typically employed for hypersonic flow simulations, but such schemes suffer from a number of disadvantages. An attractive alternative is discontinuous Galerkin methods, which benefit from arbitrarily high spatial order of accuracy, geometric flexibility, and other properties. To enable accurate computations of high-speed particle-laden flows, an Euler-Lagrange methodology was developed in which the Eulerian field of the carrier gas is calculated using a discontinuous Galerkin scheme while the disperse phase is treated with Lagrangian particle tracking. We discuss challenges associated with coupling these two formulations and how to handle them. Momentum and energy transfer between the carrier gas and the particle phase is considered, and the importance of accounting for interparticle collisions is assessed. In addition, we describe the physical model of the particle phase and examine effects of its uncertainties on the numerical solution. We demonstrate the performance of the Euler-Lagrange method in representative testcases, with focus on the accurate prediction of particle trajectories and heating augmentation. Quantitative comparisons with experiments are provided.

Ching, Eric J.↗