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 955 records · Page 53

Nonmodal Growth of TravelingWaves on Blunt Cones at Hypersonic Speeds

The existing database of transition measurements in hypersonic ground facilities has established that, as the nosetip bluntness is increased, the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream. However, this trend is reversed at sufficiently large values of the nose Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nose-tip bluntness. Because modal amplification is too weak to initiate transition at moderate-to-large bluntness values, nonmodal growth has been investigated as the potential basis for a physics-based model for the frustum transition. The present analysis investigates the nonmodal growth of traveling disturbances initiated within the nose-tip vicinity that peak within the entropy layer. Results show that, with increasing nose bluntness, both planar and oblique traveling disturbances experience appreciable energy amplification up to successively higher frequencies. For moderately blunt cones, the initial nonmmodal growth is followed by a partial decay that is more than overcome by an eventual, modal growth as Mack-mode waves. For larger bluntness values, the Mack-mode waves are not amplified anywhere upstream of the experimentally measured transition location, but the traveling modes still undergo a significant amount of nonmodal growth. This finding does not provide a definitive link between optimal growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition in the absence of sufficient Mack-mode amplification implies a double peak in disturbance amplification and the appearance of transitional events above the boundary-layer edge.

Hypersonic speed↗

The Impact of Dry Intrusions on Midlatitude Cold-Air Outbreak Cloud Transitions

Cold-air outbreaks (CAOs) form marine boundary layer (MBL) clouds that undergo rapid overcast-to-broken cloud regime transitions, initiated by substantial rain. CAOs are usually accompanied by dry intrusions (DIs) that subside as free-tropospheric (FT) air into the postfrontal sector of mid-latitude storms. For an exemplary cold-air outbreak in the NW Atlantic that showed faster transitions (corresponding to reduced extents of overcast clouds) closer to the low-pressure system, we posit that varying transitions are caused by an uneven meteorological pattern imposed by the prevailing DI. We compile satellite observations, reanalysis fields, and Lagrangian large-eddy simulations (LES) translating along MERRA2-based trajectories to show that postfrontal trajectories closer to the low-pressure system are uniquely favorable to rain formation (and, thus, cloud transitions) as they show (1) weaker FT subsidence rates, (2) greater FT humidity, (3) greater MBL windspeeds, and (4) a colder MBL as well as reduced lower-tropospheric stability. We present an updated conceptual view of postfrontal cloud formation that may guide future investigations.

Cold-air outbreaks↗

Refining Planetary Boundary Layer Height Retrievals From Micropulse‐Lidar at Multiple ARM Sites Around the World

Abstract Knowledge of the planetary boundary layer height (PBLH) is crucial for various applications in atmospheric and environmental sciences. Lidar measurements are frequently used to monitor the evolution of the PBLH, providing more frequent observations than traditional radiosonde‐based methods. However, lidar‐derived PBLH estimates have substantial uncertainties, contingent upon the retrieval algorithm used. In addressing this, we applied the Different Thermo‐Dynamic Stabilities (DTDS) algorithm to establish a PBLH data set at five separate Department of Energy's Atmospheric Radiation Measurement sites across the globe. Both the PBLH methodology and the products are subject to rigorous assessments in terms of their uncertainties and constraints, juxtaposing them with other products. The DTDS‐derived product consistently aligns with radiosonde PBLH estimates, with correlation coefficients exceeding 0.77 across all sites. This study delves into a detailed examination of the strengths and limitations of PBLH data sets with respect to both radiosonde‐derived and other lidar‐based estimates of the PBLH by exploring their respective errors and uncertainties. It is found that varying techniques and definitions can lead to diverse PBLH retrievals due to the inherent intricacy and variability of the boundary layer. Our DTDS‐derived PBLH data set outperforms existing products derived from ceilometer data, offering a more precise representation of the PBLH. This extensive data set paves the way for advanced studies and an improved understanding of boundary‐layer dynamics, with valuable applications in weather forecasting, climate modeling, and environmental studies.

54 ENVIRONMENTAL SCIENCES↗

Robust, High-Performing Maize–Perovskite-Based Solar Cells with Improved Stability

Herein, we focus on improving the long-term chemical and thermomechan-ical stability of perovskite solar cells (PSCs), two major challenges currently limiting their commercial deployment. Our strategy incorporates a long-chain starch polymer into the perovskite precursor. The starch polymer confers multiple beneficial effects by forming hydrogen bonds with the methylammonium iodide precursor, templating perovskite growth that results in a compact and homogeneous film deposited in a simple one-step coating (antisolvent-free). The inclusion of starch in the methylammonium lead iodide films strongly improves their thermomechanical and environmental stability while maintaining a high photovoltaic performance. The fracture energy (G c ) of the film is increased to above 5 J/m 2 by creating a nanocomposite that provides intrinsic reinforcement at grain boundaries. Additionally, improved optoelectronic properties achieved with the starch polymer enable good photostability of the active layer and enhanced resistance to thermal cycling.

14 SOLAR ENERGY↗

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗

Causes of and Solutions to Wind Speed Bias in NREL's 2020 Offshore Wind Resource Assessment for the California Pacific Outer Continental Shelf

This report provides the results of a detailed analysis into the causes of high wind speed bias in the 20-year wind resource data set for offshore California the National Renewable Energy Laboratory (NREL) released in 2020, herein called CA20. The data set was developed using the state-of-the-art Weather Research and Forecasting (WRF) model. Notably, no floating lidars were available at the time in offshore California to validate offshore hub-height wind speeds. In late 2020, the Pacific Northwest National Laboratory (PNNL) deployed two floating lidars in the California outer continental shelf (OCS), near the Bureau of Ocean Energy Management (BOEM) call areas of Humboldt and Morro Bay. Using these observations through 2021, NREL found considerable bias in modeled hub-height winds at both locations: up to +2 m/s at Humboldt over a 6-month period, and up to +1 m/s at Morro Bay over a one-year period. Upon the discovery of this bias, the Department of Energy (DOE) and BOEM funded NREL and PNNL to investigate the causes of, impacts of, and solutions to the bias in the CA20 data set. This report summarizes the findings of this research. We first investigated whether different WRF model setups could lead to reduced bias. We found that the choice of planetary boundary layer (PBL) scheme - which controls the vertical turbulent mixing of momentum, heat, and moisture in the lowermost part of the atmosphere - greatly affected hub-height wind speeds in the region. Specifically, switching from the Mellor-Yamada-Nakanishi-Niino (MYNN) scheme used in CA20 (and widely used across a range of operational and research weather models) to the less common Yonsei University (YSU) scheme nearly eliminated the bias at both the Humboldt and Morro Bay lidar locations. The large discrepancy between the MYNN- and YSU-modeled hub-height winds pointed towards the role of atmospheric stability. In general, PBL schemes agree well in conditions of high turbulence and mixing, normally referred to as "unstable" conditions. By contrast, PBL schemes start to diverge in "stable" conditions, where turbulence is low and thermal stratification (i.e., higher temperature air sitting on top of colder air) greatly suppresses vertical mixing. Under such conditions, winds aloft can decouple from surface effects and greatly accelerate, causing high wind speeds at hub-height and frequent low-level jets (LLJs). We determined that these stable conditions are in fact dominant in offshore California. The region is characterized by moderate-to-extreme stable stratification with a LLJ on average around 200 meters above sea-level. To our knowledge, no wind energy area globally has as strongly stable stratification as offshore California. Under these extreme conditions, we determined that the MYNN scheme models higher stability than YSU, resulting in less vertical turbulent mixing than YSU, allowing for the acceleration of hub-height winds, more intense LLJs, and higher-amplitude inertial oscillations. Using surface observations, we found that MYNN overestimates near-surface stability, whereas YSU tends to model stability better. We then considered several short-term case studies to assess additional meteorological drivers of the bias at Humboldt. We found that during synoptic scale northerly flows driven by the North Pacific High and inland thermal low, a coastal warm bias in the MYNN case studies contributes to the modeled wind speed bias by altering the boundary layer thermodynamics via a thermal wind mechanism. Given the strong performance of the YSU-based runs in offshore California, NREL has produced and published an updated version of the CA20 data set with YSU as the PBL scheme. This updated data set is now part of NREL's 2023 National Offshore Wind (NOW-23) data set, which covers all the U.S. offshore waters. The development and final validation of the NOW-23 data set in offshore California is documented in this report.

17 WIND ENERGY↗

Modeling Electrokinetics of Oxygen Electrodes in Solid Oxide Electrolyzer Cells

A microscale model is presented in this study to simulate electrode kinetics of the oxygen electrode in a solid oxide electrolyzer cell (SOEC). Two mixed ionic/electronic conducting structures are examined for the oxygen producing electrode in this work: single layer porous lanthanum strontium cobalt ferrite (LSCF), and bilayer LSCF/SCT (strontium cobalt tantalum oxide) structures. A yttrium-stabilized zirconia (YSZ) electrolyte separates the hydrogen and oxygen electrodes, as well as a gadolinium doped-ceria (GDC) buffer layer on the oxygen electrode side. Electrochemical reactions occurring at the two-phase boundaries (2PBs) and three-phase boundaries (3PBs) of single-layer LSCF and bilayer LSCF/SCT oxygen electrodes are modeled under various SOEC voltages with lattice oxygen stoichiometry as the key output. The results reveal that there exists a competition in electrode kinetics between 2PBs and 3PBs, but 3PBs are the primary reactive sites for single-layer LSCF oxygen electrode under high voltages. These locations experience the greatest oxygen stoichiometry variations and are therefore the most likely locations for dimensional changes. By applying an active SCT layer over LSCF, the 2PBs become activated to compete with the 3PBs, thus alleviating oxygen stoichiometry variations and reducing the likelihood of dimensional change. This strategy could reduce lattice structural expansion, proving to be valuable for electrode-electrolyte delamination prevention and will be the focus of future work.

30 DIRECT ENERGY CONVERSION↗

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.↗

Numerical models of thermally and mechanically coupled two-layer convection of highly viscous fluids

Thermal convection in the earth mantle is investigated by means of numerical simulations. The mantle models comprise two horizontal layers of viscous incompressible fluid (with identical or differing properties) separated by a fixed horizontal interface; an isothermal, fixed-heat-flux, or insulating lower boundary; and heat supplied either from below or internally (in the lower layer only, equally in both layers, or primarily in the upper layer). The mathematical formulation of the models is explained, and results for linear stability and finite-amplitude convection are presented in extensive tables and graphs and discussed in detail. Particular attention is given to the presence of thermal coupling without interface distortion in many cases, the predominance of long-wavelength cells in the finite-amplitude models, and the large (50 percent) temperature difference across the interface in all cases.

Ellsworth, Kirk↗

Vapor transport and crystal growth of GeSe under normal and high acceleration

Physical vapor transport experiments on GeSe in the presence of 2 atm xenon and for a nominal temperature difference of 600-500 C were performed under 1 g, 5 g, and 10 g acceleration conditions. Under high acceleration and destabilizing conditions, the GeSe crystals are generally larger than those under 1 g, stabilizing, and up to three orders of magnitude larger in surface area than those under 1 g, destabilizing conditions. The mass transport rates of the 5 g and 10 g destabilizing experiments are considerably greater than those of the 10 g, stabilizing, and 1 g experiments. The observed increase in mass flux (under destabilizing conditions) with acceleration is significantly greater than the anticipated dependence (mass flux proportional to g exp 1/4) for laminar, boundary-layer driven free convection. In view of the considerable convection under high acceleration, destabilizing conditions, the surface morphology and bulk crystallinity of the large crystal platelets are unexpectedly good.

Wiedemeier, H.↗

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗

Linear instability of curved free shear layers

The linear inviscid hydrodynamic stability of slightly curved free mixing layers is studied in this paper. The disturbance equation is solved numerically using a shooting technique. Two mean velocity profiles that represent stably and unstably curved free mixing layers are considered. Results are shown for cases of five curvature Richardson numbers. The stability characteristics of the shear layer are found to vary significantly with the introduction of the curvature effects. The results also indicate that, in a manner similar to the Goertler vortices observed in a boundary layer along a concave wall, instability modes of spatially developing streamwise vortex pairs may appear in unstable curved mixing layers.

Liou, William W.↗

Linear instability of curved free shear layers

The linear inviscid hydrodynamic stability of slightly curved free mixing layers is studied in this paper. The disturbance equation is solved numerically using a shooting technique. Two mean velocity profiles that represent stably and unstably curved free mixing layers are considered. Results are shown for cases of five curvature Richardson numbers. The stability characteristics of the shear layer are found to vary significantly with the introduction of the curvature effects. The results also indicate that, in a manner similar to the Goertler vortices observed in a boundary layer along a concave wall, instability modes of spatially developing streamwise vortex pairs may appear in centrifugally unstable curved mixing layers.

Liou, William W.↗

The winter central Arctic surface energy budget: A model evaluation using observations from the MOSAiC campaign

This study evaluates the simulation of wintertime (15 October, 2019, to 15 March, 2020) statistics of the central Arctic near-surface atmosphere and surface energy budget observed during the MOSAiC campaign with short-term forecasts from 7 state-of-the-art operational and experimental forecast systems. Five of these systems are fully coupled ocean-sea ice-atmosphere models. Forecast systems need to simultaneously simulate the impact of radiative effects, turbulence, and precipitation processes on the surface energy budget and near-surface atmospheric conditions in order to produce useful forecasts of the Arctic system. This study focuses on processes unique to the Arctic, such as, the representation of liquid-bearing clouds at cold temperatures and the representation of a persistent stable boundary layer. It is found that contemporary models still struggle to maintain liquid water in clouds at cold temperatures. Given the simple balance between net longwave radiation, sensible heat flux, and conductive ground flux in the wintertime Arctic surface energy balance, a bias in one of these components manifests as a compensating bias in other terms. This study highlights the different manifestations of model bias and the potential implications on other terms. Three general types of challenges are found within the models evaluated: representing the radiative impact of clouds, representing the interaction of atmospheric heat fluxes with sub-surface fluxes (i.e., snow and ice properties), and representing the relationship between stability and turbulent heat fluxes.

54 ENVIRONMENTAL SCIENCES↗

On the internal structure of the major satellites of the outer planets

The structures and thermal evolutions of the large icy-satellites of the outer solar system are considered. It is shown (for bodies comparable in size and mass to the Galilean satellites, having sizeable mass fraction of H2O, and with meteoritic abundances of radioactive materials contained within their silicate fractions) that the crust of solid ice over a liquid mantle predicted by conductive heat-transfer calculations is unstable to large-scale solid-state convection. For appropriate material parameters, convective heat-transfer rates are sufficient to freeze a large liquid mantle on a time scale that is short compared to the lifetime of the body. It is also concluded that the ice layer is convecting at the present time. A reevaluation of previous work, using improved values for material parameters and boundary conditions, reverses earlier conclusions and implies a rigid outer crust with resulting long-term stability of surface features to creep deformation. The combination of a rigid crust with active internal convection presents the additional possibility of surface features that are produced and maintained by dynamic internal processes.

Reynolds, R. T.↗

Factors Influencing Lava-Substrate Heat Transfer and Implications for Thermomechanical Erosion

We develop numerical simulations of basaltic lava flowing laminarly over a basalt substrate in order to examine the details of the lava dynamics and thermal boundary layers and to understand the implications for substrate heating. As the initial stage of a larger study of thermomechanical erosion in different planetary environments, we aim to understand why erosion occurs on Earth, why erosion features are not ubiquitous given the high temperatures involved, and whether it is a plausible mechanism for the formation of planetary channels such as lunar sinuous rilles and Venusian canali. Here we confine our attention to terrestrial lavas with well-known properties and eruption parameters. With relatively simple computational fluid dynamic simulations, most closely representing tube-fed hawaiian basalts (for which erosion has been documented), we demonstrate the importance of incorporating several key factors in models of lava flow/ substrate heat transfer, which have commonly been neglected in previous treatments. By addressing the interaction of the flow dynamics and heat transfer in the lava, our work suggests that the development of a temperature gradient in the base of the lava, even for undeveloped flow, has a significant influence on substrate temperature. The sensitivity of the lava-substrate interface temperature to the thermophysical properties of the lava and substrate suggests that a delicate balance is required for partial melting to occur. Thus, it might take weeks of continuous flow to initiate partial melting of the substrate at distances of several kilometers from the vent. These durations exceed the periods of stability typical of lava flowing in tubes; pauses, blockages, surges, and break-outs frequently disrupt the flow. However, natural irregularities in the flow dynamics or substrate topography might help to initiate and maintain substrate melting on shorter timescales by disturbing the intimately coupled dynamic and thermal boundary layers. Although a purely thermal mechanism cannot be ruled out, our findings support the premise that mechanical erosion may play a key role in reports of erosion based on field evidence.

Fagents, S. A.↗

Aerodynamic characteristics of NACA RM-10 missile in 8- by 6-foot supersonic wind tunnel at Mach numbers from 1.49 to 1.98 I : presentation and analysis of pressure measurements (stabilizing fins removed)

Experimental investigation of flow about a slender body of revolution (NACA RM-10 missile) aligned and inclined to a supersonic stream was conducted at Mach numbers from 1.49 to 1.98 at a Reynolds number of approximately 30,000,000. Boundary-layer measurements at zero angle of attack are correlated with subsonic formulations for predicting boundary-layer thickness and profile. Comparison of pressure coefficients predicted by theory with experimental values showed close agreement at zero angle of attack and angle of attack except over the aft leeward side of body. At angle of attack, pitot pressure measurements in plane of model base indicated a pair of symmetrically disposed vortices on leeward side of body.

LOADS - FUSELAGE, NACELLES, AND CANOPIES↗

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↗