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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Aerosol Source Plume Physical Characteristics from Space-based Multiangle Imaging

Models that assess aerosol effects on regional air quality and global climate parameterize aerosol sources in terms of amount, type, and injection height. The multiangle imaging spectroradiometer (MISR) aboard NASA's Terra satellite retrieves total column aerosol optical thickness (AOT), and aerosol type over cloud-free land and water. A stereo-matching algorithm automatically retrieves reflecting-layer altitude wherever clouds or aerosol plumes have discernable spatial contrast, with about 500-m accuracy, at 1.1-km horizontal resolution. Near-source biomass burning smoke, volcanic effluent, and desert dust plumes are observed routinely, providing information about aerosol amount, particle type, and injection height useful for modeling applications. Compared to background aerosols, the plumes sampled have higher AOT, contain particles having expected differences in Angstrom exponent, size, single-scattering albedo, and for volcanic plume and dust cloud cases, particle shape. As basic thermodynamics predicts, thin aerosol plumes lifted only by regional winds or less intense heat sources are confined to the boundary layer. However, when sources have sufficient buoyancy, the representative plumes studied tend to concentrate within discrete, high-elevation layers of local stability; the aerosol is not uniformly distributed up to a peak altitude, as is sometimes assumed in modeling. MISR-derived plume heights, along with meteorological profile data from other sources, make it possible to relate radiant energy flux observed by the moderate resolution imaging spectroradiometer (MODIS), also aboard the Terra spacecraft, to convective heat flux that plays a major role in buoyant plume dynamics. A MISR climatology of plume behavior based on these results is being developed.

aerosol plumes↗

DOE OSTI2022

Modeling plasma response to non-axisymmetric magnetic field perturbations in tokamak boundaries

The primary goal of this project is to provide physics insight into models that elucidate the mechanisms responsible for controlling the performance of H-mode plasmas and edge localized mode (ELM) stability when subjected to intrinsic and applied non-axisymmetric magnetic perturbations. The overall goal is to facilitate the development of physics-based criteria for the design of magnetic coils to be used for ELM suppression and H-mode pedestal control in future burning plasma devices such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1995

On Supersonic-Inlet Boundary-Layer Bleed Flow

Boundary-layer bleed in supersonic inlets is typically used to avoid separation from adverse shock-wave/boundary-layer interactions and subsequent total pressure losses in the subsonic diffuser and to improve normal shock stability. Methodologies used to determine bleed requirements are reviewed. Empirical sonic flow coefficients are currently used to determine the bleed hole pattern. These coefficients depend on local Mach number, pressure ratio, hole geometry, etc. A new analytical bleed method is presented to compute sonic flow coefficients for holes and narrow slots and predictions are compared with published data to illustrate the accuracy of the model. The model can be used by inlet designers and as a bleed boundary condition for computational fluid dynamic studies.

Harloff, Gary J.↗

NASA NTRS1996

Development of Three-Dimensional Flow Code Package to Predict Performance and Stability of Aircraft with Leading Edge Ice Contamination

In the work reported herein, a simplified, uncoupled, zonal procedure is utilized to assess the capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The computational approach combines potential flow plus boundary layer simulations by VSAERO for the un-iced aircraft forces and moments with Navier-Stokes simulations by NPARC for the incremental forces and moments due to iced components. These are compared with wind tunnel force and moment data, supplied by the Boeing Company, examining longitudinal flight characteristics. Grid refinement improved the local flow features over previously reported work with no appreciable difference in the incremental ice effect. The computed lift curve slope with and without empennage ice matches the experimental value to within 1%, and the zero lift angle agrees to within 0.2 of a degree. The computed slope of the un-iced and iced aircraft longitudinal stability curve is within about 2% of the test data. This work demonstrates the feasibility of a zonal method for the icing analysis of complete aircraft or isolated components within the linear angle of attack range. In fact, this zonal technique has allowed for the viscous analysis of a complete aircraft with ice which is currently not otherwise considered tractable.

Aircraft icing↗

NASA NTRS2002

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗

NASA NTRS1989

An Eigenvalue Analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L(sub 2) stability on a finite domain.

Warming, Robert F.↗

NASA NTRS1990

An eigenvalue analysis of finite-difference approximations for hyperbolic IBVPs

Warming, Robert F.↗

NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

NASA NTRS1992

Crossflow Stability and Transition Experiments in a Swept-Wing Flow

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

Dagenhart, John Ray↗

DOE OSTI2020

Tuning the Anode–Electrolyte Interface Chemistry for Garnet-Based Solid-State Li Metal Batteries

Lithium (Li) metal anode is a promising candidate for high energy density solid-state batteries. However, the interface issues, including large interfacial resistance and Li dendrite generation, have always frustrated the attempt to commercialize solid-state Li metal batteries (SSLBs). Here, we report that infusing the garnet-type solid electrolytes (GSEs) with air-stable electrolyte Li 3 PO 4 (LPO) dramatically reduces the interfacial resistance to ~1 Ω cm 2 and achieves a high critical current density of 2.2 mA cm -2 under ambient condition due to enhanced interfacial stability to Li metal anode. The coated and infused LPO electrolytes not only improve the mechanical strength and Li-ion conductivity of grain boundaries, but also form a stable Li-ion conductive but electron-insulating LPO-derived solid-electrolyte interphase between Li metal and GSEs, and consequently eliminates the Li dendrites growth and prevents the direct reduction of GSEs by Li metal over a long cycle life. This interface engineering approach together with grain-boundary modification on GSEs represents a promising strategy to revolutionize the anode-electrolyte interface chemistry for SSLBs and provides new design strategy for other-types of solid-state batteries.

25 ENERGY STORAGE↗

DOE OSTI2025

Post-build stress-relief optimization for laser powder bed fusion 316H stainless steel

Nuclear energy remains a critical component of a diversified and efficient energy portfolio, offering reliable, high-capacity, and low-carbon power. However, in the U.S., aging infrastructure and the slow qualification and deployment of advanced materials and manufacturing techniques hinder progress in next-generation reactor technologies. This study explores the application of laser powder bed fusion (LPBF) additive manufacturing for stainless steel 316H, with a focus on optimizing post-build heat treatments to enhance material properties for high-temperature nuclear applications. The research targets the optimization of stress-relief temperatures to alleviate postbuild residual stresses, ensuring improvements in the microstructural corelated properties. A series of microstructural and mechanical evaluations were performed on LPBF-printed SS-316H samples which were subjected to annealing at temperatures varying between 650 °C and 850 °C. X-ray diffraction, scanning electron microscopy, and transmission electron microscopy analyses revealed that increasing the heattreatment temperature accelerated dislocation recovery. Vickers microhardness measurements showed an initial reduction in values, followed by stabilization over extended durations at all the temperatures. While higher temperatures facilitated faster recovery, they also promoted carbide precipitation along grain and solidification cell boundaries, narrowing the safe processing window. In contrast, heat treatment at 650°C preserved the cellular substructure and enabled controlled carbide precipitation over time. In conclusion, these findings highlight the importance of time–temperature optimization and suggest that 650°C for up to 2 h provides the most favorable balance between recovery and carbide control for a stress-relief treatment.

316 stainless steel↗

DOE OSTI2022

Investigation of the Self-Assembly Behavior of Statistical Bottlebrush Copolymers via Self-Consistent Field Theory Simulations

Bottlebrush copolymers, a class of high-density side-chain-grafted copolymers with high molecular weights, have recently attracted extensive research attention due to the dramatic reduction of chain entanglement in these copolymers, unique self-assembly behaviors, and distinctive physical properties. In this study, we employ self-consistent field theory simulations to systemically investigate the phase behavior of statistical bottlebrush copolymers A-stat-B, an intriguing subclass of bottlebrush copolymers with random or alternating sequences of A and B side chains. We find that a broader variety of ordered mesophases can be stabilized in the melts of these copolymers relative to the more extensively studied bottlebrush block copolymers. In particular, we observe that these statistical copolymers can stabilize sphere phases over a larger range of species volume fractions than their diblock and tetrablock counterparts but require much higher segregation strength to achieve this. The sphere and hexagonally packed cylinder phases formed by conformationally symmetric statistical bottlebrush copolymers possess interesting discrete core–shell structures with the backbone and lower volume fraction side-chain species forming shells and cores, respectively. Moreover, we uncover a new strategy of stabilizing the Frank–Kasper phase A15 by introducing conformational asymmetry through different side-chain lengths and also observe deflection of order–order phase boundaries and significant mixing between the backbones and short side chains as conformational asymmetry is introduced.

36 MATERIALS SCIENCE↗

DOE OSTI2022

High Entropy Oxide Relaxor Ferroelectrics

Relaxor ferroelectrics are important in technological applications due to strong electromechanical response, energy storage capacity, electrocaloric effect, and pyroelectric energy conversion properties. Current efforts to discover and design materials in this class generally rely on substitutional doping as slight changes to local compositional order can significantly affect the Curie temperature, morphotropic phase boundary, and electromechanical responses. In this work, we demonstrate that moving to the strong limit of compositional complexity in an ABO 3 perovskite allows stabilization of relaxor responses that do not rely on a single narrow phase transition region. Entropy-assisted synthesis approaches are utilized to synthesize single-crystal Ba(Ti 0.2 Sn 0.2 Zr 0.2 Hf 0.2 Nb 0.2 )O 3 [Ba(5B)O] films. The high levels of configurational disorder present in this system are found to influence dielectric relaxation, phase transitions, nanopolar domain formation, and Curie temperature. Temperature-dependent dielectric, Raman spectroscopy, and second-harmonic generation measurements reveal multiple phase transitions, a high Curie temperature of 570 K, and the relaxor ferroelectric nature of Ba(5B)O films. The first-principles theory calculations are used to predict possible combinations of cations to design relaxor ferroelectrics and quantify the relative feasibility of synthesizing these highly disordered single-phase perovskite systems. Further, the ability to stabilize single-phase perovskites with various cations on the B-sites offers possibilities for designing high-performance relaxor ferroelectric materials for piezoelectric, pyroelectric, and electrocaloric applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Visualizing Millisecond Atomic Dynamics of Nanocrystals in Liquid

Atomic structures of nanomaterials are inherently dynamic and continuously reshaped through interactions with chemical species and external stimuli. Such dynamics are further amplified as the size and dimensionality of nanomaterials decrease. Despite advances in analytical methods, it remains challenging to capture the structural dynamics of nanomaterials in reactive environments with both atomic spatial resolution and commensurate temporal resolution. Here, in this study, we directly visualize atomic-scale dynamics of gold (Au) nanocrystals in reactive liquid environments with millisecond-speed liquid-cell electron microscopy (EM) and deep-learning denoising. We uncover reversible fluctuations in the local crystallinity of Au nanocrystals dependent on the surrounding chemical environment. These transient fluctuations, driven by interactions at nanocrystal–liquid interfaces, critically influence the dissolution kinetics and grain boundary relaxation. By overcoming the spatiotemporal limitations in conventional liquid-cell EM, our findings provide insights into how transient nanoscale structures dictate the stability and reactivity of nanomaterials.

Kang, Sungsu [University of Chicago, IL (United St↗

DOE OSTI2023

The Flux-Differencing Discontinuous Galerkin Method Applied to an Idealized Fully Compressible Nonhydrostatic Dry Atmosphere

Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite-volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux-Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high-order accuracy and stability for large-eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge-layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

The Role of Non-Local Effects on Surface Sensible Heat Flux Under Different Types of Thermal Structures Over the Arctic Sea-Ice Surface

The effects of atmospheric thermal structure on the surface energy flux are poorly understood over the Arctic sea-ice surface. Here, we explore the mechanism of sensible heat exchange under different types of thermal structures over the Arctic sea-ice surface by using data collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. The quadrant analysis indicates that strong surface temperature inversions below 100 m enhance non-local effects on the positive (upward) sensible heat flux ($\overline{w'θ'}$) through entrainment of large eddies from the convective boundary layer aloft. However, strong surface inversions restrict the contributions of large eddies to the negative (downward) $\overline{w'θ'}$ due to intensified surface stability. By inspecting the existing parameterization schemes, we found that the European Center for Medium-Range Weather Forecasts Integrated Forecasting System scheme fails to predict the impacts of non-local processes on the positive $\overline{w'θ'}$, and an adjustment term to correct the bias of parameterized $\overline{w'θ'}$ is proposed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Non-volatile electric-field control of inversion symmetry

Competition between ground states at phase boundaries can lead to significant changes in properties under stimuli, particularly when these ground states have different crystal symmetries. A key challenge is to stabilize and control the coexistence of symmetry-distinct phases. Using BiFeO 3 layers confined between layers of dielectric TbScO 3 as a model system, we stabilize the mixed-phase coexistence of centrosymmetric and non-centrosymmetric BiFeO 3 phases at room temperature with antipolar, insulating and polar semiconducting behaviour, respectively. Application of orthogonal in-plane electric (polar) fields results in reversible non-volatile interconversion between the two phases, hence removing and introducing centrosymmetry. Counterintuitively, we find that an electric field 'erases' polarization, resulting from the anisotropy in octahedral tilts introduced by the interweaving TbScO 3 layers. Consequently, this interconversion between centrosymmetric and non-centrosymmetric phases generates changes in the non-linear optical response of over three orders of magnitude, resistivity of over five orders of magnitude and control of microscopic polar order. Here our work establishes a platform for cross-functional devices that take advantage of changes in optical, electrical and ferroic responses, and demonstrates octahedral tilts as an important order parameter in materials interface design.

36 MATERIALS SCIENCE↗