Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parabolic”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗
DOE OSTI2025

Direct Discontinuous Galerkin methods for the reacting multi-component flow equations

The Direct Discontinuous Galerkin (DDG (Liu and Yan, 2008)) method and a counterpart with Interface Correction (DDGIC (Danis and Yan, 2022)) are extended to compute diffusion terms that arise when solving the compressible multi-component flow equations in thermochemical nonequilibrium. Thermodynamic properties, transport properties, chemical reaction rates, and energy exchange terms are computed using Mutation++ (Scoggins et al., 2020). The DG method is applied on unstructured grids, where the accuracy and convergence rates can be sensitive to the numerical method chosen for parabolic terms. A method for determining the homogeneity tensor of the flow equations required for DDGIC is shown. The convergence properties of the DDG methods are studied and compared to the Interior Penalty (IP) method. A number of numerical experiments are conducted to assess the accuracy and performance of the method. The numerical results and convergence studies indicate that DDG and DDGIC provide accurate solutions and perform well for general flows in thermochemical nonequilibrium.

Diffusion↗
DOE OSTI2021

Water vapor oxidation of SiC layer in surrogate TRISO fuel particles

Under accidental conditions for high temperature gas-cooled reactors (HTGR), the SiC layer in tri-structural-isotropic (TRISO) fuel particles can be exposed to water vapor. In this study, oxidation behaviors of surrogate TRISO fuel particles were investigated in a He-20 vol% water vapor mixed atmosphere at temperatures up to 1600 °C. The growth of the crystalline oxide passivation layer with cracks and pores followed a parabolic law with time, where the maximum was 2.3 μm at 1600 °C. The oxide layer and the SiC surface under the oxide became flattened with increasing temperature, as a function of the silica viscosity and the diffusion path of water vapor. Volatilization of the oxide layer was analyzed using a mechanistic model that an inert gas in oxidizing atmospheres could influence the magnitude of volatilization. The fracture load and strength of the oxidized and thinned SiC layer were numerically estimated to decrease from 2.27 to 1.69 N and from 317 to 299 MPa, respectively, with the SiC thickness decrease from 35 to 32 μm. This prediction indicates that the oxidized SiC layer should retain fission products. Additionally, the mechanical integrity of each layer in the TRISO fuel particle after oxidation was evaluated. The results in this work provide important data for the safety analysis of accidental scenarios in HTGRs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High temperature oxidation behaviors of bulk SiC with low partial pressures of air and water vapor in argon

Oxidation behaviors of bulk SiC at 1000–1400 °C in Ar-20 vol% H 2 O-10 vol% air were investigated. The thickness of the passivation oxide layer showed a parabolic increase with time. Dense and uniform layers formed below 1200 °C, whereas pores and cracks were observed at higher temperatures. An analytical model was proposed to estimate the oxide layer growth and SiC mass loss rate. The SiC layer should be effective in retaining the fission products after oxidation by excessive water and air. This work is the first detailed study of SiC layer oxidation during loss of coolant accidents (LOCA) for HTGRs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of exposure temperature on the corrosion behavior of a FeNiCrCuAl high entropy alloy in supercritical water

The corrosion behavior of a FeNiCrCuAl high entropy alloy exposed to supercritical water at 380–650 °C and 25 MPa was investigated. The weight gain exhibited a nearly parabolic kinetics. At 380 °C, Al 2 O 3 precipitated on Fe 3 O 4 particles in FeCr-rich regions and intensified as the temperature increased to 650 °C, eventually resulting in an Al 2 O 3 particle. A mixed Fe 3 O 4 and NiAl 2 O 4 outer layer developed on NiAl-rich regions at 380 °C, whereas an Al 2 O 3 layer formed at 650 °C. Cu 2 O formed along the boundaries of NiAl-rich regions, regardless of exposure temperatures. In conclusion, the corrosion mechanism at different temperatures was elucidated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The microscopic mechanisms of high temperature oxidation of Haynes 282

Nickel-based superalloy finds widespread applications in aerospace and extreme environments. They are known for their high temperature oxidation resistance under extended periods. However, the oxide formation and evolution which sets the stages for the later parabolic oxidation kinetics is not fully understood. This paper aims to provide new insights into the transient stage oxidation mechanism and kinetics of a typical Ni-based superalloy (Haynes 282). While tracking the chemical and microstructural evolution under micrometer scale at 800°C, we show that the oxide scale and its grain boundary species change significantly during the initial stages of oxidation and can have a profound impact on the oxidation kinetics. Cr 2 O 3 forms initially at the grain boundary along with minor amount of Al 2 O 3 . Then, the grain boundary region is enriched with copious amounts TiO 2 while Cr migrates away from the grain boundary. A noticeable change in the isothermal oxidation kinetics observed likely results from a mechanistic change from uniform surface oxidation to preferential outward diffusion of Ti 4+ ions through the grain boundary. Through first-principles calculations combined with energy dispersive spectroscopy, we confirmed the preferential outward diffusion of Ti 4+ ions as the diffusion barrier is lower for Ti than Cr along the grain boundaries. In conclusion, these findings highlight the critical role of early-stage grain boundary oxidation dynamics in dictating the long-term oxidation resistance of Ni-based superalloys and provide a foundation for future strategies to enhance their performance in extreme environments.

DFT↗
DOE OSTI2023

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SEM petrography of dispersed organic matter in black shales: A review

Organic matter (OM)-hosted pores are important constituents of the pore system of black shales and play a crucial role in determining their methane adsorption capacity and porosity. OM-hosted pores are generally observed and described with scanning electron microscope (SEM) on Ar ion-milled surfaces. However, SEM imaging is not able to reliably distinguish OM types and relate the observed pores to specific macerals. Partly because of this inability to relate organic pores to macerals, the evolution of organic porosity during thermal maturation remains poorly understood. In this paper, we review the petrographic characteristics of dispersed organic matter (DOM) in black shales under the SEM. Organic petrographic classification of DOM developed for reflected-light microscopy is so far the most practical method when describing DOM in black shales under the SEM because this classification has information on the origin of DOM. Therefore, correlative microscopy (combination of reflected-light and electron microscopy) is the most effective method to identify both OM types and OM-hosted pores. This method, however, is not readily available to most researchers. Although identifying OM on the basis of SEM observations is a challenging task, it is achievable provided there is a good understanding of the studied shales, especially their thermal maturity and original OM composition. Therefore, the overall objective of this paper is to review petrographic characteristics of DOM in black shales under the SEM to provide some guidelines for identifying DOM from SEM observations. We also review factors that control the formation and preservation of OM-hosted pores. OM-hosted pores consist of primary and secondary organic pores. Primary organic pores are pores inherited from the biological structure of the original OM. Secondary organic pores develop during hydrocarbon generation and expulsion from oil-prone OM and are hosted by solid bitumen or pyrobitumen. The development of secondary organic pores is controlled by thermal maturity and OM type, and their preservation is subject to thermal maturity, OM content, and mineralogical composition. The presented view of the evolution of micropore and mesopore characteristics of OM with thermal maturity is based on data from the literature. The specific surface area and pore volume of OM in black shales follow a parabolic pattern with increasing thermal maturity (quantified via vitrinite reflectance, R o ). The initial increase reflects development of OM-hosted pores, and the subsequent decrease is due to denser stacking of aromatic units in the macromolecular structure of OM, with maximum values (specific surface area ~ 300 m 2 /g and pore volume ~ 0.3 cm 3 /g) reached at R o values in the 2.5-3.5% range. The contribution of OM-hosted pores to the pore characteristics of black shales depends on OM content, OM type, and thermal maturity.

03 NATURAL GAS↗
DOE OSTI2021

Design and characterization of a single-channel microwave interferometer for the Helicon Physics Prototype eXperiment

A heterodyne interferometer for measuring the plasma electron density in Helicon Physics Prototype eXperiment (HPPX) has been developed and characterized in the laboratory. The heterodyne-type interferometer uses a probing beam of 85 GHz and a reference beam of 85.7 GHz by using two Gunn oscillators for facilitating frequency modulation. The cutoff density for this frequency is 8.96 x 10 19 m –3 that is much higher than the nominal electron density of the HPPX plasmas. The maximum deflection angle due to the refraction effect is estimated to be 0.63° assuming a parabolic electron density profile, and the impact on the measurement is negligible for the plasmas with a size of 100 mm. Two identical HDPE lenses act as the vacuum windows, and the lenses working together with the transmitter and receiver horns, form a beam with the beam waist of about 20 mm located at the plasma center. The characterization results of the lens and the critical microwave components show good performance. The plasma electron density measurement results are in line with the expectation. Furthermore, the Signal to Noise Ratio (SNR) of the measurement is defined as SNR = 20lg($\bar{n}_{e}$/S), where $\bar{n}_{e}$ and S are respectively the average electron density and the standard deviation of the signal, and the SNR is 80 dB. Finally, this demonstrates the system's capability of providing a clear and detailed density resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Alluvial fan morphology: A self-similar free boundary problem description

In this work, we examine approximate geometrically self-similar solutions to a parabolic free boundary value problem applied to alluvial fan surface morphology and growth. Alluvial fans are fan- or cone-shaped sedimentary deposits caused by the rapid deposition of sediment from a canyon discharging onto a flatter plain. Longitudinal, topographic profiles of fans can be readily described by a seemingly time independent dimensionless profile (DeChant et al., 1999). However, because an alluvial fan can be expected to grow over time, it is clear that this “steady” profile is certainly time dependent and can be described using a space-time self-similar solution. In an experimental and theory-based study, Guerit et al. (2014) developed a self-similar (or as they describe it a self-affine) linear solution based upon an approximate first order small parameter expansion solution for a 1-d homogeneous nonlinear diffusion equation. Direct substitution of this result into a linear diffusion equation suggests that this first order expression may not fully satisfy the associated governing equation. In contrast, we develop a more complete solution based upon a modeled approximation for the axi-symmetric formulation such that the associated temporal behavior is consistent with a 1/3 time power-law as described by Reitz and Jerolmack (2014). The resulting expression is an exact solution to a linear heat equation. Furthermore, we emphasize that a small parameter is not inherent to the resulting profile result and is not included in our model development. Though developed using rather different approaches, the formal solution developed here is in good agreement with the simple polynomial described by DeChant et al. (1999) suggesting that this self-similar solution is a suitable time dependent representation of alluvial fan longitudinal profile form and improves on earlier work.

58 GEOSCIENCES↗
DOE OSTI2026

Multi-arrival infrasound from meteoroids: Fragmentation signatures versus propagation effects in a fine-scale layered atmosphere

Infrasonic signatures of meteoroid fragmentation are frequently ambiguous: do multiple arrivals signify a complex breakup or merely the distorting effects of a layered atmosphere? Resolving this ambiguity is critical for accurate energy estimates and source reconstruction. In this study, we address this challenge by analyzing a unique regional dataset of well-constrained meteoroid events observed by the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. We employ pseudo-differential parabolic equation (PPE) simulations to quantify how fine-scale gravity-wave structures in the stratosphere and lower thermosphere modify acoustic waveforms at ranges <300 km. Our modeling reveals that while fine-scale layering can stretch signals and generate diffuse oscillatory tails, it does not produce discrete, high-amplitude pulse splitting at ranges below ∼140 km. By applying these results to the rare multi-arrival event 20060305, we demonstrate that its distinct double arrival at 100 km range is inconsistent with atmospheric multipathing and provides definitive evidence of separate fragmentation episodes. These findings establish new diagnostic criteria for separating source physics from propagation artifacts, improving the reliability of infrasound as a monitoring tool for natural bolides, space debris re-entries, and catastrophic launch failures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Gas nitriding behavior of refractory metals and implications for multi-principal element alloy design

Multi-principal element alloys (MPEAs) comprise a large, flexible compositional space that enables tuning of their chemistry, structure, and properties. Here, to facilitate the development of nitriding-based surface-enhancement strategies that harness a broad compositional space, this study examined the gas nitriding behavior of Hf, Mo, Nb, Ta, Ti, and Zr as a function of time, temperature (750 and 1000 °C), and nitriding potential (i.e. ammonia-to-hydrogen ratio). These metals were selected because they have a strong driving force to form nitrides, and appear in many promising refractory MPEA compositions. The nitriding temperatures were selected based on the phase transformation temperature of Ti and Zr, and the nitriding potentials were chosen such that all elements are expected to form nitrides. Mass gain measurements indicate that all six elements follow parabolic kinetics. The microstructure observations and quantitative microchemical analysis show formation of dense and well-adhered compound layers for Mo, Nb, and Ta. Thick diffusion zones appear in Hf, Ta, Ti, and Zr, and diffusion coefficients were fit to the composition profiles. Partial delamination of the compound layer occurred for Ti and Zr. Peak hardness values above 30 GPa are obtained in the dense compound layers, and the solute hardening of the underlying alloy is correlated with the nitrogen content. The results provide insight into the dynamics of nitride compound formation relative to interstitial dissolution of nitrogen, and are discussed in the context of MPEA composition and processing design.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A unified formulation of splitting-based implicit time integration schemes

Splitting-based time integration approaches such as fractional step, alternating direction implicit, operator splitting, and locally one dimensional methods partition the system of interest into components, and solve individual components implicitly in a cost-effective way. Here this work proposes a unified formulation of splitting time integration schemes in the framework of general-structure additive Runge–Kutta (GARK) methods. Specifically, we develop implicit-implicit (IMIM) GARK schemes, provide the order conditions for this class, and explain their application to partitioned systems of ordinary differential equations. We show that classical splitting methods belong to the IMIM GARK family, and therefore can be studied in this unified framework. New IMIM-GARK splitting methods are developed and tested using parabolic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Neural operator prediction of linear instability waves in high-speed boundary layers

Here, we investigate if neural operators can predict the linear evolution of instability waves in high-speed boundary layers. To this end, we extend the design of the DeepOnet to ensure accurate and robust predictions, and also to perform data assimilation. In particular, we train DeepONet to take as inputs an upstream disturbance and a downstream location of interest, and to provide as output the perturbation field downstream in the boundary layer. DeepONet thus approximates the linearized and parabolized Navier-Stokes operator for this flow. For successful application to the high-speed boundary layer problem, we add sample weighting and Fourier input features to the regular DeepONet formulation. Once trained, the DeepOnet can perform fast and accurate predictions of the downstream disturbances within the range of training frequencies (inside the distribution). In addition, we show that DeepONet can solve the inverse problem, where downstream wall measurements are adopted as input, and a trained network can predict the upstream disturbances that led to these observations. This capability, along with the forward predictions, allows us to perform a full data assimilation cycle efficiently: starting from wall-pressure data, we predict the upstream disturbance using the inverse DeepONet and its evolution using the forward DeepONet. Finally, we introduce three new metrics to benchmark the training, evaluation and break-even cost of neural operators.

42 ENGINEERING↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗