Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Layer Stability”

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 973 records · Page 54

Phase-field modeling of crack growth and mitigation in solid oxide cells

Fracture and crack growth is one of the main degradation mechanisms in solid oxide cells (SOCs). However, the modeling of crack growth in SOCs is challenging due to their complex microstructures and possible plasticity development within the Ni particles in Ni-based SOC electrodes. In this study, a phase-field fracture model is developed, which incorporates the SOC microstructures and phase-dependent material properties, including yield strength, fracture toughness in the bulk and at the interphase boundaries. The model is employed to study crack initiation and growth under thermal and redox cycling on the hydrogen electrode side of SOCs. The simulation results demonstrate that under thermal cycling, work-zone cracking dominates in electrolyte-supported SOCs with cracks initiated at the triple-phase boundaries, while only minor mechanical degradation occurs in hydrogen-electrode-supported SOCs after hundreds of thermal cycles. Under redox cycling, through-cracking of yttria-stabilized zirconia (YSZ) in the hydrogen electrode and electrolyte layers dominates. Finally, the simulation results suggest several crack-mitigation strategies, including decreasing the porosity in the hydrogen electrode support layer and synchronizing thermal strain to balance oxidation strain.

08 HYDROGEN↗

Natural laminar flow airfoil design considerations for winglets on low-speed airplanes

Winglet airfoil section characteristics which significantly influence cruise performance and handling qualities of an airplane are discussed. A good winglet design requires an airfoil section with a low cruise drag coefficient, a high maximum lift coefficient, and a gradual and steady movement of the boundary layer transition location with angle of attack. The first design requirement provides a low crossover lift coefficient of airplane drag polars with winglets off and on. The other requirements prevent nonlinear changes in airplane lateral/directional stability and control characteristics. These requirements are considered in the design of a natural laminar flow airfoil section for winglet applications and chord Reynolds number of 1 to 4 million.

Vandam, C. P.↗

Stability, transition and turbulence

A glimpse is provided of the research program in stability, transition and turbulence based on numerical simulations. This program includes both the so-called abrupt and the restrained transition processes. Attention is confined to the prototype problems of channel flow and the parallel boundary layer in the former category and the Taylor-Couette flow in the latter category. It covers both incompressible flows and supersonic flows. Some representative results are presented.

Hussaini, M. Y.↗

Stability, transition and turbulence

A glimpse is provided of the research program in stability, transition, and turbulence based on numerical simulations. This program includes both the so-called abrupt and the restrained transition processes. Attention is confined to the prototype problems of channel flow and the parallel boundary layer in the former category and the Taylor-Couette flow in the latter category. It covers both incompressible flows and supersonic flows. Some representative results are presented.

Hussaini, M. Y.↗

Predicted and hot-film measured Tollmien-Schlichting wave characteristics

The Tollmien-Schlichting (TS) instability is a time-dependence instability which can lead to transition of laminar boundary layers on airfoils. A comparison of theoretical predictions and experimental observations of the TS instability on the NLF(1)-0414F airfoil designed by Viken and Pfenninger. The theoretical predictions were obtained using the SALLY stability code. Test results, from the same hot films that were used to detect transition, revealed that TS waves could be detected by the hot films if the hot-film signal was adequately modified.

Stack, John P.↗

Effect of Target-type Thrust Reverser on Transonic Aerodynamic Characteristics of a Single-engine Fighter Model

A brief investigation of a target-type thrust reverser on a single-engine fighter model has been conducted in the Langley 16-foot transonic tunnel at Mach numbers from 0.20 to 1.05.At Mach numbers of 0.80, 0.92, and 1.05, a hydrogen peroxide turbojet-engine simulator was operated with the thrust reverser extended. The angle of attack was varied from 0 degrees to 5 degrees at these Mach numbers. The Reynolds number of the free stream, based on the mean aerodynamic chord, was about 5 x 10(6). It was estimated that reversed jet operations separated the model boundary-layer flow over the upper surface of the horizontal tail and upper part of the afterbody. This resulted in a positive pitch increment due to reversed jet operation. Jet-on operation also tended to stabilize the severe lateral oscillations which occurred with the reverser extended and the jet off. It appeared that these jet-off oscillations were the result of an alternating separation and reattachment of the flow on the rearmost portions of the fuselage afterbody.

Swihert, John M↗

Instability of a Supersonic Boundary-Layer with Localized Roughness

A localized 3-D roughness causes boundary-layer separation and (weak) shocks. Most importantly, streamwise vortices occur which induce streamwise (low U, high T) streaks. Immersed boundary method (volume force) suitable to represent roughness element in DNS. Favorable comparison between bi-global stability theory and DNS for a "y-mode" Outlook: Understand the flow physics (investigate "z-modes" in DNS through sinuous spanwise forcing, study origin of the beat in DNS).

Marxen, Olaf↗

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

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↗

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗

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↗

The influence of wind shear on aerodynamic coefficients

The purpose of this study is to investigate the influence of wind shear on the lift, drag, roll and yaw moments of a wing in a horizontal wind gradient at various elevations and roll angles. The models of wind shear considered are those proposed by Leurs (1973) for atmospheric flow over horizontally homogeneous and uniform terrain. A general series solution for the distribution of lift along the wing span following the method of Houghton and Brock (1970) is employed to compute the loads and moments on the airfoil. Results of the computations indicate that wind shear can have a significant effect on the rolling and yawing moments of the wing of an aircraft flying with one wing low in the atmospheric boundary layer. These moments are directly influenced by the magnitude of the wind gradient in the neutral atmosphere reflected by the magnitude of the friction velocity. The effect of wind shear on lift and drag is negligible. Stability of the atmosphere tends to increase the magnitude of the rolling moment while decreasing its variation with elevation.

Frost, W.↗

Relaminarization of turbulent flow on a flat plate by localized surface heating

Relaminarization of a turbulent boundary layer in air on a flat plate is demonstrated experimentally using localized surface heating. Localized heating is achieved by electrically heating a wire embedded in a thermally insulated substrate (Space Shuttle Tile) on the surface. The stability of the flow downstream of the applied control point increases with decreasing stream temperature in the flow direction. The mean and perturbation velocity profiles without control show that the flow is turbulent. With control, these profiles sequentially change from intermittently turbulent to a fully laminar state. In the relaminarization stage, the turbulent energy is dissipated by molecular transport due to viscous and conductivity mechanisms. The new profile adjusts to a lower Reynolds number based on the momentum thickness than that of the previous turbulent state.

Maestrello, Lucio↗

Blunt-nosed swept supercritical LFC wings without nose flaps

To simplify the design of swept supercritical (SC) Laminar Flow Control (LFC) wings and maintain satisfactory low speed characteristics, blunt-nosed swept SC LFC wings without nose flaps and lower wing loadings were studied. Their boundary layer crossflow in the leading edge area is optimally controlled (1) by compensating the boundary layer crossflow of the front acceleration zone by an opposite crossflow in a downstream pressure rise area, (2) by maintaining a neutrally stable boundary layer crossflow by suction within a narrow spanwise suction strip located close to the wing attachment line in the front acceleration zone. The required suction massflow and power are then very small, especially considering the strongly stabilizing effect of surface and streamline curvature on crossflow stability.

Pfenninger, W.↗

Deposition of ozone to tundra

Eddy correlation measurements of O3 deposition fluxes to tundra during the Arctic Boundary Layer Expedition (ABLE 3A) are reported. The mean O3 deposition velocity was 0.24 cm/s in the daytime and 0.12 cm/s at night. The day-to-day difference in deposition velocity was driven by both atmospheric stability and surface reactivity. The mean surface resistance to O3 deposition was 2.6 s/cm in the daytime and 3.4 s/cm at night. The relatively low surface resistance at night is attributed to light-insensitive uptake of O3 at dry upland tundra surfaces. The small day-tonight difference in surface resistance is attributed to additional stomatal uptake by wet meadow tundra plants in the daytime. The mean O3 deposition flux to the world north of 60 deg N in July-August is estimated at 8.2 x 10 exp 10 molecules/sq cm/s. Suppression of photochemical loss by small anthropogenic inputs of nitrogen oxides could have a major effect on O3 concentrations in the summertime Arctic troposphere.

Jacob, D. J.↗

An Experimental and Numerical Study of Roughness-Induced Instabilities in a Mach 3.5 Boundary Layer

Progress on a joint experimental and numerical study of laminar-to-turbulent transition induced by an isolated roughness element in a high-speed laminar boundary layer is reported in this paper. The numerical analysis suggests that transition is driven by the instability of high- and low-speed streaks embedded in the wake of the isolated roughness element. In addition, spatial stability analysis revealed that the wake flow supports multiple modes (even and odd) of convective instabilities that experience strong enough growth to cause transition. The experimental measurements, which included hot-wire and pitot-probe surveys, confirmed the existence of embedded high- and low-speed streaks in the roughness wake. Furthermore, the measurements indicate the presence of both even and odd modes of instability, although their relative magnitude depends on the specifics of the roughness geometry and flow conditions (e.g., the value of Re(sub kk) or k/delta. For the two test cases considered in the measurements (Re(sub kk) values of 462 and 319), the even mode and the odd mode were respectively dominant and appear to play a primary role in the transition process.

Kegerise, Michael A.↗

Statistically Resolved Planetary Boundary Layer Height Diurnal Variability Using Spaceborne Lidar Data

The Planetary Boundary Layer Height (PBLH) significantly impacts weather, climate, and air quality. Understanding the global diurnal variation of the PBLH is particularly challenging due to the necessity of extensive observations and suitable retrieval algorithms that can adapt to diverse thermodynamic and dynamic conditions. This study utilized data from the Cloud-Aerosol Transport System (CATS) to analyze the diurnal variation of PBLH in both continental and marine regions. By leveraging CATS data and a modified version of the Different Thermo-Dynamics Stability (DTDS) algorithm, along with machine learning denoising, the study determined the diurnal variation of the PBLH in continental mid-latitude and marine regions. The CATS DTDS-PBLH closely matches ground-based lidar and radiosonde measurements at the continental sites, with correlation coefficients above 0.6 and well-aligned diurnal variability, although slightly overestimated at nighttime. In contrast, PBLH at the marine site was consistently overestimated due to the viewing geometry of CATS and complex cloud structures. The study emphasizes the importance of integrating meteorological data with lidar signals for accurate and robust PBLH estimations, which are essential for effective boundary layer assessment from satellite observations.

54 ENVIRONMENTAL SCIENCES↗

A numerical method for solving the Navier-Stokes equations with application to shock-boundary layer interactions

A numerical method for solving the compressible form of the unsteady Navier-Stokes equations is described. This method was originally presented in 1970 and has since been modified during the development of computer programs at Ames for implementing models that account for the effects of turbulence in shock-induced separated flows. Although this paper does not describe the turbulence models themselves, a complete description of the basic numerical method is given with emphasis on the choice of a computational mesh for high Reynolds number flows, finite-difference approximations for mixed partial derivatives, extension of the Courant-Friedrichs-Lewy stability condition for viscous flows, mesh boundary conditions, and numerical smoothing for strong shock-wave calculations.

Maccormack, R. W.↗