Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-stage preconditioners for thermal–compositional–reactive flow in porous media

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.

42 ENGINEERING↗
NASA NTRS1991

Growth and properties of amorphous silicon films grown using pulsed-flow reactive plasma beam epitaxy

The growth and properties of a-Si:H films grown using a novel deposition technique, reactive plasma beam epitaxy, are discussed. In this technique, a remote H plasma produced in a microwave-ECR reactor is used to grow a-Si:H films at low pressures. The H ions react with SiH4 introduced near the substrate to produce the film. The flow of SiH4 is pulsed on or off, thereby achieving in-situ annealing of the film during growth by H ions and radicals. The films produced by this technique appear to have good electronic quality, and are more stable than the standard glow discharge films.

Dalal, Vikram L.↗
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A hybrid adaptive multiresolution approach for the efficient simulation of reactive flows

Computational studies that use block-structured adaptive mesh refinement (AMR) approaches suffer from unnecessarily high mesh resolution in regions adjacent to important solution features. This deficiency limits the performance of AMR codes. In this work a novel hybrid adaptive multiresolution (HAMR) approach to AMR-based calculations is introduced to address this issue. The multiresolution (MR) smoothness indicators are used to identify regions of smoothness on the mesh where the computational cost of individual physics solvers may be decreased by replacing direct calculations with interpolation. We suggest an approach to balance the errors due to the adaptive discretization and the interpolation of physics quantities such that the overall accuracy of the HAMR solution is consistent with that of the MR-driven AMR solution. The performance of the HAMR scheme is evaluated for a range of test problems, from pure hydrodynamics to turbulent combustion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy

A novel and efficient algorithm is presented in this paper to deal with DNS of turbulent reacting flows under the low-Mach-number assumption, with detailed chemistry and a quasi-spectral accuracy. The temporal integration of the equations relies on an operating-split strategy, where chemical reactions are solved implicitly with a stiff solver and the convection-diffusion operators are solved with a Runge-Kutta-Chebyshev method. The spatial discretisation is performed with high-order compact schemes, and a FFT based constant-coefficient spectral solver is employed to solve a variable-coefficient Poisson equation. The numerical implementation takes advantage of the 2DECOMP&FFT libraries developed by, which are based on a pencil decomposition method of the domain and are proven to be computationally very efficient. An enhanced pressure-correction method is proposed to speed up the achievement of machine precision accuracy. It is demonstrated that a second-order accuracy is reached in time, while the spatial accuracy ranges from fourth-order to sixth-order depending on the set of imposed boundary conditions. In conclusion, the software developed to implement the present algorithm is called HOLOMAC, and its numerical efficiency opens the way to deal with DNS of reacting flows to understand complex turbulent and chemical phenomena in flames.

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

Pore-scale modelling of multiphase reactive flow: application to mineral dissolution with production of CO 2

A micro-continuum approach is proposed to simulate the dissolution of solid minerals at the pore scale in the presence of multiple fluid phases. The approach employs an extended Darcy–Brinkman–Stokes formulation that accounts for the interfacial tension between the two immiscible fluid phases and the moving contact line at the mineral surface. The simulation framework is validated using an experimental microfluidic device that provides time-lapse images of the dissolution dynamics. The set-up involves a single-calcite crystal and the subsequent generation of CO 2 bubbles in the domain. The dissolution of the calcite crystal and the production of gas during the acidizing process are analysed. We then show that the production of CO 2 bubbles during the injection of acid in a carbonate formation may limit the overall dissolution rate and prevent the emergence of wormholes.

bio-inspired, mechanical behavior, carbon sequestr↗
DOE OSTI2014

A tomographic technique for the simultaneous imaging of temperature, chemical species, and pressure in reactive flows using absorption spectroscopy with frequency-agile lasers

This paper proposes a technique that can simultaneously retrieve distributions of temperature, concentration of chemical species, and pressure based on broad bandwidth, frequency-agile tomographic absorption spectroscopy. The technique holds particular promise for the study of dynamic combusting flows. A proof-of-concept numerical demonstration is presented, using representative phantoms to model conditions typically prevailing in near-atmospheric or high pressure flames. The simulations reveal both the feasibility of the proposed technique and its robustness. Our calculations indicate precisions of ∼70 K at flame temperatures and ∼0.05 bars at high pressure from reconstructions featuring as much as 5% Gaussian noise in the projections.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The reactive flow evolution of the polymer-bonded explosive PBX 9502: Experiments and model validation in extreme pressure regimes

The shock-to-detonation transition properties of the triaminotrinitrobenzene based PBX 9502 high explosive (HE) are experimentally and computationally explored in extremely high input pressure conditions. These include both slightly sub-Chapman–Jouguet and overdriven input pressure conditions, namely, ∼25 and ∼31 GPa, respectively. Our experiments capture the transient buildup of a shock-induced reaction via measurement of HE and polymethyl methacrylate window interface particle velocity profiles for a variety of sample thicknesses for this insensitive HE. These observations necessitate extremely thin explosive samples, and the high rates of reaction provide a considerable challenge to optical diagnostics. Samples at these thicknesses also provide an opportunity for evaluation of potential micro-structure effects on the resulting shock-to-detonation-transition measurements. To address this, the thin samples are also characterized via x-ray micro-computed tomography. Finally, a pair of previously established continuum-level detonation performance modeling approaches for PBX 9502 were used to analyze the experiments. The employed model variants crucially differ in their definition of each model’s empirical reaction rate functional form, utilization of shock state quantities, and local flow variable dependencies. As a result, the present experiments provide a novel platform to evaluate the quantitative and qualitative consequences stemming from these modeling choices in a challenging initiation scenario, largely beyond the chosen calibration range of either model. This new experimental information will provide a platform for both improved physics and model parameterizations for this well-studied explosive.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗