Search NASA⌕ Search

Engineering topics

Xu, Chao

Publications and source records attributed to Xu, Chao.

The Tribomechadynamics Research Challenge: Confronting blind predictions for the linear and nonlinear dynamics of a thin-walled jointed structure with measurement results

The present article summarizes the submissions to the Tribomechadynamics Research Challenge announced in 2021. The task was a blind prediction of the vibration behavior of a system comprising a thin plate clamped on two sides via bolted joints. Both geometric and frictional contact nonlinearities are expected to be relevant. Provided were the CAD models and technical drawings of all parts as well as assembly instructions. The main objective was to predict the frequency and damping ratio of the lowest-frequency mode as function of the amplitude. Many different prediction approaches were pursued, ranging from well-known methods to very recently developed ones. After the submission deadline, the system has been fabricated and tested. The aim of this article is to evaluate the current state of the art in modeling and vibration prediction, and to provide directions for future methodological advancements.

42 ENGINEERING↗

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗

Large eddy simulation of lean blow-off in swirl-stabilized flame with the subgrid dissipation concept

The lean blow-off mechanism of the premixed swirl flame is numerically investigated by large eddy simulation (LES) with the subgrid dissipation concept (SDC) combustion model. Three simulated cases cover stable, near blow-off, and transient conditions. Compared with the experiment, the LES-SDC approach captures the flow and combustion features for stable and near blow-off conditions. More importantly, the predictions of the blowoff procedure and duration agree satisfactorily with the experiment, indicating that the LES-SDC approach is a promising tool for predicting strong, unsteady turbulent combustion processes. Further, the numerical results are used to investigate the blow-off mechanism. Two stages in the blow-off procedure are specified. The first is the necking and extinction of the downstream flame surface, and the second is the shrinking of the upstream flame surface. The blow-off mechanism is well explained by the theory of stretched flame extinction. At the end of the recirculation zone, the large negative radial velocity pushes the flame to the central line. The combustion process here can be abstracted as the stretched counter-flame of the reactant-reactant configuration. The excessive flame stretch dominates the flame extinction and triggers the blow-off event. The upstream flame resists the intense stretch with the help of hot product recirculation, and the combustion here can be idealized as the counter-flame of the reactant-product configuration. The alignment of the temperature gradient and flow velocity, together with the excessive stretch, clearly indicates the tendency of flame local extinction. A Damköhler number-based blow-off criterion is raised from the mean flow strain rate and laminar flame bulk extinction strain rate.

42 ENGINEERING↗

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗

Independent Fuel Property Effects of Fuel Volatility on Low Temperature Heat Release and Fuel Autoignition (Final Report)

This Cooperative Research and Development Agreement (CRADA) project between Argonne National Laboratory (ANL), Oak Ridge National Laboratory (ORNL), and Shell Global Solutions (Shell) was initiated as part of a Directed Funding Opportunity (DFO) call for proposals from the Co-Optimization of Fuels and Engines (Co-Optima) initiative. Shell had observed that volatile fuels suppress low temperature heat release (LTHR) more than expected based on conventional gasoline autoignition metrics: research octane number (RON) and motor octane number (MON). The role of LTHR contributes to autoignition phenomena for both boosted spark ignition (BSI) and advanced compression ignition (ACI) combustion modes. ACI combustion modes are applicable to large engines in the hard-to-electrify applications such as off-road, rail, and marine. Thus, having a reliable understanding of autoignition phenomena, including being able to accurately account for the effects of fuel volatility, is particularly important as new synthetic and bio-fuel compositions are considered.This CRADA project aimed to test the hypothesis that the decreased LTHR is due to preferential evaporation of multicomponent fuels when using direct injection (DI) fueling technology, creating composition and reactivity stratification. A custom set of fuels was designed and blended to test this hypothesis by Shell, with experimental engine studies at ORNL and engine combustion modeling by ANL. However, the initial experimental findings did not show the expected effect of fuel volatility suppressing LTHR. Instead, the LTHR propensity observed was independent of the fuel volatility. Due to the unexpected experimental result, the remainder of the experimental effort was redirected to study the effect of fuel volatility on emissions under spark-ignited cold-start conditions. However, as with the LTHR experiments, the cold start effort did not show a meaningful effect of fuel volatility on cold start emissions. Meanwhile, improved engine CFD models have been developed for both LTHR and cold start operations for the Shell fuels with different volatilities. While the simulation efforts were not pursued further due to the insignificant effects of fuel volatility as shown in experiments, the models developed can be easily retooled for off-road, rail, and marine applications.

09 BIOMASS FUELS↗

Simulations of n -dodecane/oxygen/nitrogen cellular detonations

In this work, two-dimensional n-dodecane/air/nitrogen cellular detonations are simulated with various equivalence ratios (ERs). A skeletal mechanism consisting of 54 species and 269 reactions is used. The lower and upper equivalence ratio boundaries for self-sustained detonation are 0.3 and 2.2, respectively. Detonation with different regimes characterized by the detonation cell patterns is observed, which aligns well with the category based on the stability parameter, i.e., weakly and highly unstable detonations, and extinction. Further, in terms of the frontal structure, non-negligible effect of diffusion on the cellular detonation is revealed, especially in the vicinity of the leading shock front. In highly unstable and quenched detonations, the alternation in reaction pathway within the induction zone accounts for the changes of detonation dynamics, such as the absence or extended sequence of important radical formation, e.g., OH. In addition, the composition of the unburned pockets depends on both pocket location from the leading shock front and the ER in the fresh mixture, because the former determines the residence time, whilst the latter affects the pocket reaction rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A numerical study of ignition and flame development characteristics in GCI combustion using large eddy simulations and chemical explosive mode analysis

Here, this work investigates the ignition and flame development processes of low reactivity fuel combustion under compression ignition conditions based on the large eddy simulation approach. The chemical explosive mode analysis (CEMA) is employed to characterize the local combustion features, including gas-liquid fuel zone, auto-ignition, diffusion-assisted, extinction, cool flame and post-ignition zone, among which auto-ignition and post-ignition are found to play a key role in the overall heat release process. The local flame propagation modes in gasoline compression ignition (GCI) are determined by quantifying the relative magnitude of diffusion/chemistry at a representative progress variable in the pre-ignition zone. The results show that autoignition fronts and deflagration waves exist simultaneously in the ignition and intense high temperature heat release (HTHR) stages, but autoignition fronts dominate. In addition, the chemical kinetic processes of four heat release periods are analyzed. The heat release during the ignition period is found to be dominated by the reactions CH 3 + H (+M) <=> CH 4 (+M) and CH 3 CHO + H <=> CH 2 CHO + H 2 . The reaction CH 2 OH + OH <=> CH 2 O + H 2 O always plays an important role in the heat releases during the other three combustion stages including intense HTHR, moderate HTHR and post-combustion.

33 ADVANCED PROPULSION SYSTEMS↗

Grid resolution requirement of chemical explosive mode analysis for large eddy simulations of premixed turbulent combustion

Full Article Figures & data References Citations Metrics Reprints & Permissions Read this article Abstract The grid resolution requirement for trustworthy Chemical Explosive Mode Analysis (CEMA) in Large Eddy Simulation (LES) of premixed turbulent combustion is proposed. Explicit filtering, to emulate the effect of the LES filter, is applied to one-dimensional laminar flame and three-dimensional planar turbulent flames across a wide range of Karlovitz numbers (5 - 239). The identification of the flame front by CEMA is found relatively insensitive to the cell size (Δ), while the combustion mode identification shows more significant sensitivity. Specifically, increasing Δ falsely enhances the auto-ignition and local extinction modes and suppresses the diffusion-assisted mode. Limited dependence of the CEMA performance on the turbulent combustion regime (Karlovitz number) is observed. A simple grid size criterion for reliable CEMA mode identification in LES is proposed as Δ ≲ δ L /2; The criterion can be relaxed to Δ ≲ δ L in the laminar flame limit. Furthermore, theoretical analysis is conducted on an idealised chemistry-diffusion system. The effects of the filtering process and turbulence on the local combustion mode are demonstrated, which is consistent with the numerical observations. Further, by incorporating turbulent combustion models in CEMA, potential improvement in identifying local combustion modes can be expected.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafiltration separation of Am(VI)-polyoxometalate from lanthanides

Partitioning of americium from lanthanides (Ln) present in used nuclear fuel plays a key role in the sustainable development of nuclear energy. This task is extremely challenging because thermodynamically stable Am(III) and Ln(III) ions have nearly identical ionic radii and coordination chemistry. Oxidization of Am(III) to Am(VI) produces AmO 2 2+ ions distinct with Ln(III) ions, which has the potential to facilitate separations in principle. However, the rapid reduction of Am(VI) back to Am(III) by radiolysis products and organic reagents required for the traditional separation protocols including solvent and solid extractions hampers practical redox-based separations. Herein, we report a nanoscale polyoxometalate (POM) cluster with a vacancy site compatible with the selective coordination of hexavalent actinides ( 238 U, 237 Np, 242 Pu and 243 Am) over trivalent lanthanides in nitric acid media. To our knowledge, this cluster is the most stable Am(VI) species in aqueous media observed so far. Ultrafiltration-based separation of nanoscale Am(VI)-POM clusters from hydrated lanthanide ions by commercially available, fine-pored membranes enables the development of a once-through americium/lanthanide separation strategy that is highly efficient and rapid, does not involve any organic components and requires minimal energy input.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗