Design and testing of a new planetary boundary layer parameterization
A new planetary boundary layer parameterization was developed. Tests included forecasting experiments.
SEARCH · Engineering Papers
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.
A new planetary boundary layer parameterization was developed. Tests included forecasting experiments.
A marine planetary boundary layer (PBL) model is presented and compared with data from sea-based experiments. The PBL model comprises two layers, the outer an Ekman-Taylor layer with stratification-dependent secondary flow, and the logarithmic surface layer corrected for stratification and humidity effects and variable surface roughness. Corrections are noted for air much warmer than water in stable conditions and for low wind speeds. The layers are analytically defined along with similarity relations and a resistance law for inclusion in a program. An additional interfacial layer correction is developed and shown to be significant for heat flux calculations. Experimental data from GOASEX were used to predict the windfield in the Gulf of Alaska, and JASIN data was used for windfields SE of Iceland. The JASIN-derived wind field predictions were accurate to within 1 m/sec and 10 deg in a 200 km triangle.
The transition of the laminar boundary layer into the turbulent state, which results in an increased drag, can be avoided by sucking of the boundary layer particles near the wall. The technically-interesting case of sucking the particles using individual slits is investigated for bodies of revolution in incompressible flow. The results of the variational calculations show that there is an optimum suction height, where the slot separations are maximum. Combined with favorable shaping of the body, it is possible to keep the boundary layer over bodies of revolution laminar at high Reynolds numbers using relatively few suction slits and small amounts of suction flow.
Exact solution of the laminar-boundary-layer equations for wedge-type flow with constant property values are presented for transpiration-cooled surfaces with variable wall temperatures. The difference between wall and stream temperature is assumed proportional to a power of the distance from the leading edge. Solutions are given for a Prandtl number of 0.7 and ranges of pressure-gradient, cooling-air-flow, and wall-temperature-gradient parameters. Boundary-layer profiles, dimensionless boundary-layer thicknesses, and convective heat-transfer coefficients are given in both tabular and graphical form. Corresponding results for constant wall temperature and for impermeable surfaces are included for comparison purposes.
Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon similar correlation techniques that have been is use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. Recent entry flight testing performed with the Orbiter Discovery now provides a means to validate this engineering correlation approach to higher confidence. These results only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, and to discuss the characteristics which are desirable in a next generation advancement. The details of the paper will review the experimental facilities and techniques that were utilized to perform the implementation of the Orbiter RTF BLT Vsn 2 prediction capability. Statistically significant results for multiple engineering correlations from a ground testing campaign will be reviewed in order to describe why only certain correlations were selected for complete implementation to support the Shuttle Program. Historical Orbiter flight data on early boundary layer transition due to protruding gap fillers will be described in relation to the selected empirical correlations. In addition, Orbiter entry flight testing results from the BLT Flight Experiment will be discussed in relation to these correlations. Applicability of such correlations to the entry design problem will be reviewed, and finally a perspective on the desirable characteristics for a next generation capability based on high fidelity physical models will be provided.
The discrete spectrum of the Orr-Sommerfeld problem of hydrodynamic stability for boundary layer flows in semi-infinite regions is examined. Related questions concerning the continuous spectrum are also addressed. Emphasis is placed on the stability problem for the Blasius boundary layer profile. A general theoretical result is given which proves that the discrete spectrum of the Orr-Sommerfeld problem for boundary layer profiles (U(y), 0,0) has only a finite number of discrete modes when U(y) has derivatives of all orders. Details are given of a highly accurate numerical technique based on collocation with splines for the calculation of stability characteristics. The technique includes replacement of 'outer' boundary conditions by asymptotic forms based on the proper large parameter in the stability problem. Implementation of the asymptotic boundary conditions is such that there is no need to make apriori distinctions between subcases of the discrete spectrum or between the discrete and continuous spectrums. Typical calculations for the usual Blasius problem are presented.
The law of the wake and the law of the wall in incompressible turbulent boundary layers formulated by Coles (1956) and their use by Mathews et al. (1970) in the development of a wall-wake representation of the velocity profile in a form applicable for isoenergetic compressible boundary layers are extended to a modified wall-wake velocity profile for turbulent compressible boundary layers. The modified wall-wake profile is shown to provide good representations of experimental velocity distributions.
Due to their fast evolution and large natural variability in macro- and microphysical properties, the accurate representation of boundary layer clouds in current climate models remains a challenge. One of the regions with large intermodel spread in the Coupled Model Intercomparison Project Phase 6 ensemble is the western North Atlantic Ocean. Here, statistically representative in situ measurements can help to develop and constrain the parameterization of clouds in global models. To this end, we performed comprehensive measurements of boundary layer clouds, aerosol, trace gases, and radiation in the western North Atlantic Ocean during the NASA Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) mission. In total, 174 research flights with 574 flight hours for cloud and precipitation measurements were performed with the HU-25 Falcon during three winter (February–March 2020, January–April 2021, and November 2021–March 2022) and three summer seasons (August–September 2020, May–June 2021, and May–June 2022). Here we present a statistical evaluation of 16 140 individual cloud events probed by the fast cloud droplet probe and the two-dimensional stereo cloud probe during 155 research flights in a representative and repetitive flight strategy allowing for robust statistical data analyses. We show that the vertical profiles of distributions of the liquid water content and the cloud droplet effective diameter (ED) increase with altitude in the marine boundary layer. Due to higher updraft speeds, higher cloud droplet number concentrations (N liquid ) were measured in winter compared to summer despite lower cloud condensation nucleus abundance. Flight cloud cover derived from statistical analysis of in situ data is reduced in summer and shows large variability. This seasonal contrast in cloud coverage is consistent with a dominance of a synoptic pattern in winter that favors conditions for the formation of stratiform clouds at the western edge of cyclones (post-cyclonic). In contrast, a dominant summer anticyclone is concomitant with the occurrence of shallow cumulus clouds and lower cloud coverage. The evaluation of boundary layer clouds and precipitation in the N liquid ED phase space sheds light on liquid, mixed-phase, and ice cloud properties and helps to categorize the cloud data. Ice and liquid precipitation, often masked in cloud statistics by a high abundance of liquid clouds, is often observed throughout the cloud. The ACTIVATE in situ cloud measurements provide a wealth of cloud information useful for assessing airborne and satellite remote-sensing products, for global climate and weather model evaluations, and for dedicated process studies that address precipitation and aerosol–cloud interactions.
This presentation summarizes the computations of a set of shock wave / turbulent boundary layer interaction (SWTBLI) test cases using the Wind-US code, as part of the 2010 American Institute of Aeronautics and Astronautics (AIAA) shock / boundary layer interaction workshop. The experiments involve supersonic flows in wind tunnels with a shock generator that directs an oblique shock wave toward the boundary layer along one of the walls of the wind tunnel. The Wind-US calculations utilized structured grid computations performed in Reynolds-averaged Navier-Stokes mode. Three turbulence models were investigated: the Spalart-Allmaras one-equation model, the Menter Shear Stress Transport wavenumber-angular frequency two-equation model, and an explicit algebraic stress wavenumber-angular frequency formulation. Effects of grid resolution and upwinding scheme were also considered. The results from the CFD calculations are compared to particle image velocimetry (PIV) data from the experiments. As expected, turbulence model effects dominated the accuracy of the solutions with upwinding scheme selection indicating minimal effects.!
The present scheme for the prediction of boundary-layer development in small, low-speed wind tunnel contraction sections proceeds by calculating the wall pressure distributions, and hence the wall velocity distributions, by means of a three-dimensional potential-flow method. For the family of contractions presently treated, the assumption of a laminar boundary layer appears to be justified; the measured boundary layer momentum thicknesses at the exit of the four contractions were found to lie within 10 percent of predicted values.
The underlying premise in this work is that large structures in turbulent shear flows are not universal and similar large structures in different flows owe their origins to basically similar instabilities. The turbulent boundary-layer flow over a backward-facing step and also the reattaching transitional boundary layer behind a large spanwise rod have been studied experimentally in parallel primarily to understand the origin of the curious large structures observed in the far downstream region (after 39 step heights) of the former flowfield. Simultaneous flow visualization and hot-wire anemometry and also plain-flow visualization have been carried out. In the backward-facing step, the three-dimensional mixing layer developing after detachment persists even after reattachment. The flow does not recover fully back to a regular turbulent boundary layer even far downstream of reattachment. The turbulence-producing instability of a turbulent boundary layer continues to compete with the wall-bounded mixing layer instability, giving rise to the curious schismatic arrays of large structures.
The method described utilizes a nonorthogonal coordinate system for boundary-layer calculations. It includes a geometry program that represents the wing analytically, and a velocity program that computes the external velocity components from a given experimental pressure distribution when the external velocity distribution is not computed theoretically. The boundary layer method is general, however, and can also be used for an external velocity distribution computed theoretically. Several test cases were computed by this method and the results were checked with other numerical calculations and with experiments when available. A typical computation time (CPU) on an IBM 370/165 computer for one surface of a wing which roughly consist of 30 spanwise stations and 25 streamwise stations, with 30 points across the boundary layer is less than 30 seconds for an incompressible flow and a little more for a compressible flow.
Two and three-dimensional shock-boundary layer interaction data were obtained from supersonic wind tunnel tests. These interactions are studied both with and without boundary layer bleed. The data verify computational fluid dynamic codes. Surface static pressure, pitot pressure, flow angularity, and bleed rates, are studied by flow visualization techniques. Surface oil flow using fluorescent dye and laser sheet using water droplets as the scattering material are used for flow visualization.
Preliminary results are presented of a finite element/finite difference method (semidiscrete Galerkin method) used to calculate compressible boundary layer flow about airfoils, in which the group finite element scheme is applied to the Dorodnitsyn formulation of the boundary layer equations. The semidiscrete Galerkin (SDG) method promises to be fast, accurate and computationally efficient. The SDG method can also be applied to any smoothly connected airfoil shape without modification and possesses the potential capability of calculating boundary layer solutions beyond flow separation. Results are presented for low speed laminar flow past a circular cylinder and past a NACA 0012 airfoil at zero angle of attack at a Mach number of 0.5. Also shown are results for compressible flow past a flat plate for a Mach number range of 0 to 10 and results for incompressible turbulent flow past a flat plate. All numerical solutions assume an attached boundary layer.
An accurate method was developed for predicting effects of streamline curvature and coordinate system rotation on turbulent boundary layers. A new two-equation model of turbulence was developed which serves as the basis of the study. In developing the new model, physical reasoning is combined with singular perturbation methods to develop a rational, physically-based set of equations which are, on the one hand, as accurate as mixing-length theory for equilibrium boundary layers and, on the other hand, suitable for computing effects of curvature and rotation. The equations are solved numerically for several boundary layer flows over plane and curved surfaces. For incompressible boundary layers, results of the computations are generally within 10% of corresponding experimental data. Somewhat larger discrepancies are noted for compressible applications.
The laminar boundary-layer flow about a circular cone at large angles of attack to a supersonic stream has been analyzed in the plane of symmetry by a method applicable in general to the flow about conical bodies. At the bottom of the cone, velocity profiles were obtained showing the expected tendency of the boundary layer to become thinner on the under side of the cone as the angle of attack is increased. At the top of the cone, the analysis failed to yield unique solutions except for small angles of attack. Beyond a certain critical angle of attack, boundary-layer flow does not exist in the plane of symmetry, thus indicating separation. This critical angle is presented as a function of Mach number and cone vertex angle.
Abstract Uncrewed aircraft systems (UAS) demonstrate significant potential for filling data gaps in the atmospheric boundary layer. However, the extent to which UAS observations—typically vertical profiles taken over 15 min—are representative of the boundary layer as a whole remains poorly characterized. Using large eddy simulations (LES) of the daytime convective boundary layer (CBL), we quantify random errors in UAS measurements that occur due to insufficient statistical convergence of the time average to the true ensemble mean. Random errors in first‐order moments increase as the CBL becomes increasingly unstable, and are largest near the surface for most quantities. Errors are on the order of 2–6 m for wind speed, 15–60 for wind direction, 0.2–3 K for potential temperature, and 0.1–1 g for specific humidity, with errors in turbulent fluxes on the order of 50%–100%. Sampling strategies that mitigate random errors are discussed in light of our results.
The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.