Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Crossflow Stability and Transition Experiments in Swept-Wing Flow

An experimental examination of crossflow instability and transition on a 45deg swept wing was conducted in the Arizona State University Unsteady Wind Tunnel. The stationary-vortex pattern and transition location are visualized by using both sublimating chemical and liquid-crystal coatings. Extensive hot-wire measurements were obtained at several measurement stations across a single vortex track. The mean and travelling wave disturbances were measured simultaneously. Stationary crossflow disturbance profiles were determined by subtracting either a reference or a span-averaged velocity profile from the mean velocity data. Mean, stationary crossflow, and traveling wave velocity data were presented as local boundary layer profiles and contour plots across a single stationary crossflow vortex track. Disturbance mode profiles and growth rates were determined. The experimental data are compared with predictions from linear stability theory.

Dagenhart, J. Ray↗
NASA NTRS2014

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Nielsen, Tanner B.↗
NASA NTRS1990

Experiments on swept-wing boundary layers

Three-dimensional boundary-layer experiments are currently being conducted on a 45-deg swept wing in the Arizona State University Unsteady Wind Tunnel. Crossflow-dominated transition is produced via a model with contoured end liners to simulate infinite swept-wing flow. Fixed-wavelength stationary and traveling crossflow vortices are observed. The stationary vortex wavelengths vary with Reynolds number as predicted by linear-stability theory, but with observed wavelengths which are about 25 percent smaller than theoretically predicted. The frequencies of the most amplified moving waves are in agreement with linear stability theory; however, traveling waves at higher frequencies than predicted are also observed. These higher-frequency waves may be harmonics of the primary crossflow waves generated by a parametric resonance phenomena. Boundary-layer profiles measured at several spanwise locations show streamwise disturbance profiles characteristic of the crossflow instability.

Dagenhart, J. R.↗
NASA NTRS1992

Prediction of high frequency combustion instability in liquid propellant rocket engines

The present use of a numerical model developed for the prediction of high-frequency combustion stabilities in liquid propellant rocket engines focuses on (1) the overall behavior of nonlinear combustion instabilities (2) the effects of acoustic oscillations on the fuel-droplet vaporization and combustion process in stable and unstable engine operating conditions, oscillating flowfields, and liquid-fuel trajectories during combustion instability, and (3) the effects of such design parameters as inlet boundary conditions, initial spray conditions, and baffle length. The numerical model has yielded predictions of the tangential-mode combustion instability; baffle length and droplet size variations are noted to have significant effects on engine stability.

Kim, Y. M.↗
DOE OSTI2022

A Converse to Lieb–Robinson Bounds in One Dimension Using Index Theory

Abstract Unitary dynamics with a strict causal cone (or “light cone”) have been studied extensively, under the name of quantum cellular automata (QCAs). In particular, QCAs in one dimension have been completely classified by an index theory. Physical systems often exhibit only approximate causal cones; Hamiltonian evolutions on the lattice satisfy Lieb–Robinson bounds rather than strict locality. This motivates us to study approximately locality preserving unitaries (ALPUs). We show that the index theory is robust and completely extends to one-dimensional ALPUs. As a consequence, we achieve a converse to the Lieb–Robinson bounds: any ALPU of index zero can be exactly generated by some time-dependent, quasi-local Hamiltonian in constant time. For the special case of finite chains with open boundaries, any unitary satisfying the Lieb–Robinson bound may be generated by such a Hamiltonian. We also discuss some results on the stability of operator algebras which may be of independent interest.

Physics↗
DOE OSTI2023

Heuristics for Molten-Salt Synthesis of Single-Crystalline Ultrahigh-Nickel Layered Oxide Cathodes

In pursuit of Li-ion batteries with higher energy density, ultrahigh-nickel layered oxides are a leading candidate for next-generation cathode materials. Single-crystalline morphology offers a neat solution to the poor stability of ultrahigh-Ni cathodes; a lower active surface area mitigates electrolyte decomposition at high voltages, and the elimination of grain boundaries improves mechanical resilience and increases volumetric energy density. However, single-crystal cathodes possess their own challenges, several of which originate from synthesis at elevated temperatures meant to induce grain growth. Molten-salt synthesis is an alternative method for obtaining single crystals, accelerating grain growth through the presence of a molten flux without the need for increased temperature. Herein, we offer heuristic guidelines for molten-salt synthesis, discussing key factors for designing reaction mixtures and the necessary exploratory research for novel molten salt/cathode systems. Here, the influence of different salts and synthesis conditions on the morphology and properties of single-crystal LiNiO 2 is presented. It is found that oxidative salts, such as Li 2 O 2 and LiNO 3 , are crucial to supplementing dissolution of gaseous oxygen into the molten phase. Through these discussions, this work aims to provide a set of overarching principles for obtaining higher-quality single-crystal layered oxide cathodes and engender more rigorous and impactful investigation into their fundamental nature and applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamics of Three-Dimensional Shock-Wave/Boundary-Layer Interactions

Advances in measuring and understanding separated, nominally two-dimensional (2D) shock-wave/turbulent-boundary-layer interactions (STBLI) have triggered recent campaigns focused on three-dimensional (3D) STBLI, which display far greater configuration diversity. Nonetheless, unifying properties emerge for semi-infinite interactions, taking the form of conical asymptotic behavior where shock-generator specifics become insignificant. The contrast between 2D and 3D separation is substantial; the skewed vortical structure of 3D STBLI reflects the essentially 2D influence of the boundary layer on the 3D character of the swept shock. As with 2D STBLI, conical interactions engender prominent spectral content below that of the turbulent boundary layer. However, the uniform separation length scale, which is crucial to normalizing the lowest-frequency dynamics in 2D STBLI, is absent. Comparatively, the spectra of 3D STBLI are more representative of the mid-frequency, convective, shear-layer dynamics in 2D, while phenomena associated with 2D separation-shock breathing are muted. Asymptotic behavior breaks down in many regions important to 3D-STBLI dynamics, occurring in a configuration-dependent manner. Aspects of inceptive regions near shock generators and symmetry planes are reviewed. Focused efforts toward 3D modal and nonmodalanalyses, moving-shock/boundary-layer interactions, fluid/structure interactions, and flow control are suggested as directions for future work.

42 ENGINEERING↗
DOE OSTI2021

Temperature controlled friction extrusion process

Processing temperature is decisive to the properties and performance of the product in friction extrusion. However, the constantly varying thermal boundary conditions result in an unstable processing temperature even fixed die rotational speed and feed rate were implemented, which extensively disturb the stability of the process and the quality of the product. It is key to maintain the process temperature to ensure the properties of the product is uniform. What's more, a controlled processing temperature enables the properties of the product to be manipulatable between parts and between different locations of one part

Li, Xiao↗
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method was developed that could predict the pressure distribution of a ducted source from far field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) was solved using explicit marching techniques. The Von Neumann method was used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach was successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method is developed that could predict the pressure distribution of a ducted source from far-field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) is solved using explicit marching techniques. The Von Neumann method is used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach is successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗
NASA NTRS1989

Numerical Study of Three-dimensional Spatial Instability of a Supersonic Flat Plate Boundary Layer

The behavior of spatially growing three-dimensional waves in a supersonic boundary layer was studied numerically by solving the complete Navier-Stokes equations. Satisfactory comparison with linear parallel and non-parallel stability theories, and experiment are obtained when a small amplitude inflow disturbance is used. The three-dimensional unsteady Navier-Stokes equations are solved by a finite difference method which is fourth-order and second-order accurate in the convection and viscous terms respectively, and second-order accurate in time. Spanwise periodicity is assumed. The inflow disturbance is composed of eigenfunctions from linear stability theory. By increasing the amplitude of the inflow disturbance, nonlinear effects in the form of a relaxation type oscillation of the time signal of rho(u) are observed.

Maestrello, Lucio↗
NASA NTRS2016

Regional Changes in Earths Color and Texture as Observed From Space Over a 15-Year Period

Earth-observing satellites provide global observations of many geophysical variables. As these variables are derived from measured radiances, the underlying radiance data are the most reliable sources of information for change detection. Here, we identify statistically significant trends in the color and spatial texture of the Earth as viewed from multiple directions from the Multi-angle Imaging SpectroRadiometer (MISR), which has been sampling the angular distribution of scattered sunlight since 2000. Globally, our results show that the Earth has been appearing relatively bluer (up to 1.6 % per decade from both nadir and oblique views) and smoother (up to 1.5 % per decade only from oblique views) over the past 15 years. The magnitude of the global blueing trends is comparable to that of uncertainties in radiometric calibration stability. Regional shifts in color and texture, which are significantly larger than global means, are observed, particularly over polar regions, along the boundaries of the subtropical highs, the tropical western Pacific, Southwestern Asia, and Australia. We demonstrate that the large regional trends cannot be explained either by uncertainties in radiometric calibration or variability in total or spectral solar irradiance; hence, they reflect changes internal to the Earths climate system. The 15-year-mean true color composites and texture images of the Earth at both nadir and oblique views are also presented for the first time.

regional trends↗
NASA NTRS2025

Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics

A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.

Process Modeling↗
DOE OSTI2021

Reactive boride infusion stabilizes Ni-rich cathodes for lithium-ion batteries

Engineered polycrystalline electrodes are critical to the cycling stability and safety of lithium-ion batteries, yet it is challenging to construct high-quality coatings at both the primary- and secondary-particle levels. Here, we present a room-temperature synthesis route to achieve full surface coverage of secondary particles and facile infusion into grain boundaries, thus offering a complete “coating-plus-infusion” strategy. Cobalt boride metallic glass is successfully applied to Ni-rich layered cathode LiNi 0.8 Co 0.1 Mn 0.1 O 2 . Here, it dramatically improves the rate capability and cycling stability, including under high-discharge-rate and elevated-temperature conditions and in pouch full cells. The superior performance originates from simultaneous suppression of microstructural degradation of intergranular cracking and side reactions with electrolyte. Atomistic simulations identified the critical role of strong selective interfacial bonding, which offers not only a large chemical driving force to ensure uniform reactive wetting and facile infusion but also lowered the surface/interface oxygen activity, contributing to the exceptional mechanical and electrochemical stabilities of the infused electrode.

25 ENERGY STORAGE↗
DOE OSTI2024

Characterization of Precipitation-Strengthening Heat-Resistant Austenitic Stainless Steels for Life-Prediction Modeling

In this study, the role of minor alloying additions in 347H stainless steels (UNS34709, ASTM A240/240M) on creep-rupture properties at 650-750°C and microstructure evolution during isothermal exposure at 750°C has been investigated, aiming to provide the experimental dataset as boundary conditions of physics-based modeling for material/component life prediction. Four different 347H heats containing various amounts of boron and nitrogen additions were prepared and evaluated. The combined additions of B and N are found to stabilize the strengthening secondary M23C6 carbides and retarding the transition from M23C6 to sigma phase precipitates during thermal exposure. The observed kinetics of microstructure evolution reasonably explains the improvement of creep-rupture properties of 347H stainless steels with the B and N additions.

Yamamoto, Yukinori↗
NASA NTRS1981

Estimating ocean-air heat fluxes during cold air outbreaks by satellite

Nomograms of mean column heating due to surface sensible and latent heat fluxes were developed. Mean sensible heating of the cloud free region is related to the cloud free path (CFP, the distance from the shore to the first cloud formation) and the difference between land air and sea surface temperatures, theta sub 1 and theta sub 0, respectively. Mean latent heating is related to the CFP and the difference between land air and sea surface humidities q sub 1 and q sub 0 respectively. Results are also applicable to any path within the cloud free region. Corresponding heat fluxes may be obtained by multiplying the mean heating by the mean wind speed in the boundary layer. The sensible heating estimated by the present method is found to be in good agreement with that computed from the bulk transfer formula. The sensitivity of the solutions to the variations in the initial coastal soundings and large scale subsidence is also investigated. The results are not sensitive to divergence but are affected by the initial lapse rate of potential temperature; the greater the stability, the smaller the heating, other things being equal. Unless one knows the lapse rate at the shore, this requires another independent measurement. For this purpose the downwind slope of the square of the boundary layer height is used, the mean value of which is also directly proportional to the mean sensible heating. The height of the boundary layer should be measurable by future spaceborn lidar systems.

Chou, S. H.↗