Search NASA⌕ Search

SEARCH · Search NASA

Results for “Roughness sublayer”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Aeolian Sediment Transport through Large Patches of Roughness in the Atmospheric Inertial Sublayer

Roughness influences the flux of wind driven sand transport. In this paper, we report on sediment transport measurements for four different surface roughness configurations composed of the same size (solid) roughness elements in the atmospheric inertial sublayer (ISL). Results of these tests indicate that sediment transport rates through patches of roughness in the atmospheric inertial sublayer are to a large extent controlled and scale proportionally with the roughness density (lambda = n b h/S, where n is number of elements of breadth b and height h in area S) of the surface. However, element size apparently increases the magnitude of the reduction beyond that attributable to lambda. A sediment transport model that incorporates the effect of shear stress partitioning appears to predict reasonably well the effect of roughness on sand transport in the cases where the roughness elements are less than or equal to 0.10 m in height. However, when the dimensions of the roughness itself are equivalent to or are greater than the range of saltation lengths (vertical and horizontal), additional interactions of the elements with the saltation cloud appear to reduce the transport efficiency.

Gillies, John A.↗
NASA NTRS2020

Combining MODIS LAI with ICESat-Based Canopy Heights Improves Spaceborne Estimates of Vegetation Roughness Length for Momentum

Most land-surface models require parameterization of vertical wind profiles within the atmospheric boundary layer. For vegetated surfaces, it is common to assume a logarithmic profile in the surface layer, which includes estimates of vegetation roughness length for momentum (z0) and zero-plane displacement height (d0). This study finds that remotely-sensed forest canopy heights improve estimates of aerodynamic roughness length for momentum using a previously-developed representation of the roughness sublayer (Raupach 1992; Jasinski et al. 2005). Resulting roughness products consist of two datasets: 1) 14 years of 8-day snapshots of the global land surface at a nominal spatial resolution of 500-meters for users who wish to retain full temporal resolution and interannual variability; and 2) multiyear averages of the 8-day snapshots, here referred to as "climatologies" of roughness, which retain underlying seasonality. Both products are suitable for use in data assimilation and reanalyses such as the National Climate Assessment Land Data Assimilation System (NCA-LDAS), for which these products were initially developed.

Borak, Jordan S.↗
NASA NTRS2005

Bulk Surface Momentum Parameters for Satellite-Derived Vegetation Fields

The bulk aerodynamic parameters associated with the absorption of surface momentum by vegetated landscapes are theoretically estimated within the context of Raupach's roughness sublayer formulation. The parameters include the bulk plant drag coefficient, maximum u*/U(sub h), sheltering coefficient, and canopy area density at onset of sheltering. Parameters are estimated for the four principal IGBP land cover classes within the U.S. Southern Great Plains: evergreen needleleaf forests, grasslands, croplands, and open shrublands. The estimation approach applies the Method of Moments to roughness data from several international field experiments and other published sources. The results provide the necessary land surface parameters for satellite-based estimation of momentum aerodynamic roughness length and zero-plane displacement height for seasonally variable vegetation fields employed in most terrestrial and atmospheric simulation models used today. Construction of sample displacement and roughness maps over the Southern United States using MODIS land products demonstrates the potential of this approach for regional to global applications.

Jasinski, Michael F.↗
NASA NTRS1987

A comparative study of some mathematical models of the mean wind structure and aerodynamic drag of plant canopies

A semianalytical method for describing the mean wind profile and shear stress within plant canopies and for estimating the roughness length and the displacement height is presented. This method incorporates density and vertical structure of the canopy and includes simple parameterizations of the roughness sublayer and shelter factor. Some of the wind profiles examined are consistent with first-order closure techniques while others are consistent with second-order closure techniques. Some profiles show a shearless region near the base of the canopy; however, none displays a secondary maximum there. Comparing several different analytical expressions for the canopy wind profile against observations suggests that one particular type of profile (an Airy function which is associated with the triangular foliage surface area density distribution) is superior to the others. Because of the numerical simplicity of the methods outlined, it is suggested that they may be profitably used in large-scale models of plant-atmosphere exchanges.

Massman, William↗
NASA NTRS1994

An evaluation of the regional acid deposition model surface module for ozone uptake at three sites in the San Joaquin Valley of California

Plants and soils act as major sinks for the destruction of tropospheric ozone, especially during daylight hours when plant stomata open and are thought to provide the dominant pathway for the uptake of ozone. The present study, part of the California Ozone Deposition Experiment, compares predictions of the regional acid deposition model ozone surface conductance module with surface conductance data derived from eddy covariance measurements of ozone flux taken at a grape, a cotton, and a grassland site in the San Joaquin Valley of California during the summer of 1991. Results indicate that the model (which was developed to provide long-term large-area estimates for the eastern United States) significantly overpredicts the surface conductance at all times of the day for at least two important types of plant cover of the San Joaquin Valley and that it incorrectly partitions the ozone flux between transpiring and nontranspiring components of the surface at the third site. Consequently, the model either overpredicts or inaccurately represents the observed deposition velocities. Other results indicate that the presence of dew does not reduce the rate of ozone deposition, contradicting to model assumptions, and that model assumptions involving the dependency of stomata upon environmental temperature are unnecessary. The effects of measurement errors and biases, arising from the presence of the roughness sublayer and possible photochemical reactions, are also discussed. A simpler model for ozone surface deposition (at least for the San Joaquin Valley) is proposed and evaluated.

Massman, W. J.↗
NASA NTRS2005

Effect of Stability on Mixing in Open Canopies

In open canopies, the within-canopy flux from the ground surface and understory can account for a significant fraction of the total flux above the canopy. This study incorporates the important influence of within-canopy stability on turbulent mixing and subcanopy fluxes into a first-order closure scheme. Toward this goal, we analyze within-canopy eddy-correlation data from the old aspen site in the Boreal Ecosystem - Atmosphere Study (BOREAS) and a mature ponderosa pine site in Central Oregon, USA. A formulation of within-canopy transport is framed in terms of a stability- dependent mixing length, which approaches Monin-Obukhov similarity theory above the canopy roughness sublayer. The new simple formulation is an improvement upon the usual neglect of the influence of within-canopy stability in simple models. However, frequent well-defined cold air drainage within the pine subcanopy inversion reduces the utility of simple models for nocturnal transport. Other shortcomings of the formulation are discussed.

Lee, Young-Hee↗
NASA NTRS2007

Shear Stress Partitioning in Large Patches of Roughness in the Atmospheric Inertial Sublayer

Drag partition measurements were made in the atmospheric inertial sublayer for six roughness configurations made up of solid elements in staggered arrays of different roughness densities. The roughness was in the form of a patch within a large open area and in the shape of an equilateral triangle with 60 m long sides. Measurements were obtained of the total shear stress (tau) acting on the surfaces, the surface shear stress on the ground between the elements (tau(sub S)) and the drag force on the elements for each roughness array. The measurements indicated that tau(sub S) quickly reduced near the leading edge of the roughness compared with tau, and a tau(sub S) minimum occurs at a normalized distance (x/h, where h is element height) of approx. -42 (downwind of the roughness leading edge is negative), then recovers to a relatively stable value. The location of the minimum appears to scale with element height and not roughness density. The force on the elements decreases exponentially with normalized downwind distance and this rate of change scales with the roughness density, with the rate of change increasing as roughness density increases. Average tau(sub S): tau values for the six roughness surfaces scale predictably as a function of roughness density and in accordance with a shear stress partitioning model. The shear stress partitioning model performed very well in predicting the amount of surface shear stress, given knowledge of the stated input parameters for these patches of roughness. As the shear stress partitioning relationship within the roughness appears to come into equilibrium faster for smaller roughness element sizes it would also appear the shear stress partitioning model can be applied with confidence for smaller patches of smaller roughness elements than those used in this experiment.

Gillies, John A.↗
NASA NTRS1991

Characteristics of the Martian atmosphere surface layer

Researchers extend elements of various terrestrial boundary layer models to Mars in order to estimate sensible heat, latent heat, and momentum fluxes within the Martian atmospheric surface layer. To estimate the molecular viscosity and thermal conductivity of a CO2-H2O gas mixture under Martian conditions, parameterizations were developed. Parameterizations for specific heat and and binary diffusivity were also determined. The Prandtl and Schmidt numbers derived from these thermophysical properties were found to range from 0.78 - 1.0 and 0.47 - 0.70, respectively, for Mars. Brutsaert's model for sensible and latent heat transport within the interfacial sublayer for both aerodynamically smooth and rough airflow was experimentally tested under similar conditions, validating its application to Martian conditions. For the surface sublayer, the researchers modified the definition of the Monin-Obukhov length to properly account for the buoyancy forces arising from water vapor gradients in the Martian atmospheric boundary layer. This length scale was then utilized with similarity theory turbulent flux profiles with the same form as those used by Businger et al. and others. It was found that under most Martian conditions, the interfacial and surface sublayers offer roughly comparable resistance to sensible heat and water vapor transport and are thus both important in determining the associated fluxes.

Clow, G. D.↗
NASA NTRS1990

Characteristics of the Martian atmosphere surface layer

Elements of various terrestrial boundary layer models are extended to Mars in order to estimate sensible heat, latent heat, and momentum fluxes within the Martian atmospheric surface ('constant flux') layer. The atmospheric surface layer consists of an interfacial sublayer immediately adjacent to the ground and an overlying fully turbulent surface sublayer where wind-shear production of turbulence dominates buoyancy production. Within the interfacial sublayer, sensible and latent heat are transported by non-steady molecular diffusion into small-scale eddies which intermittently burst through this zone. Both the thickness of the interfacial sublayer and the characteristics of the turbulent eddies penetrating through it depend on whether airflow is aerodynamically smooth or aerodynamically rough, as determined by the Roughness Reynold's number. Within the overlying surface sublayer, similarity theory can be used to express the mean vertical windspeed, temperature, and water vapor profiles in terms of a single parameter, the Monin-Obukhov stability parameter. To estimate the molecular viscosity and thermal conductivity of a CO2-H2O gas mixture under Martian conditions, parameterizations were developed using data from the TPRC Data Series and the first-order Chapman-Cowling expressions; the required collision integrals were approximated using the Lenard-Jones potential. Parameterizations for specific heat and binary diffusivity were also determined. The Brutsart model for sensible and latent heat transport within the interfacial sublayer for both aerodynamically smooth and rough airflow was experimentally tested under similar conditions, validating its application to Martian conditions. For the surface sublayer, the definition of the Monin-Obukhov length was modified to properly account for the buoyancy forces arising from water vapor gradients in the Martian atmospheric boundary layer. It was found that under most Martian conditions, the interfacial and surface sublayers offer roughly comparable resistance to sensible heat and water vapor transport and are thus both important in determining the associated fluxes.

Clow, G. D.↗
NASA NTRS1990

Free convection in the Matian atmosphere

The 'free convective' regime for the Martian atmospheric boundary layer (ABL) was investigated. This state occurs when the mean windspeed at the top of the ABL drops below some critical value U(sub c) and positive buoyant forces are present. Such forces can arise either from vertical temperature or water vapor gradients across the atmospheric surface layer. During free convection, buoyant forces drive narrow plumes that ascend to the inversion height with a return circulation consisting of broad slower-moving downdraughts. Horizontal pressure, temperature, windspeed, and water vapor fluctuations resulting form this circulation pattern can be quite large adjacent to the ground (within the surface layer). The local turbulent fluctuations cause non-zero mean surface stresses, sensible heat fluxes, and latent heat fluxes, even when the mean regional windspeed is zero. Although motions above the surface layer are insensitive to the nature of the surface, the sensible and latent heat fluxes are primarily controlled by processes within the interfacial sublayer immediately adjacent to the ground during free convection. Thus the distinction between aerodynamically smooth and rough airflow within the interfacial sublayer is more important than for the more typical situation where the mean regional windspeed is greater than U(sub c). Buoyant forces associated with water vapor gradients are particularly large on Mars at low pressures and high temperatures when the surface relative humidity is 100 percent, enhancing the likelihood of free convection under these conditions. On this basis, Ingersol postulated the evaporative heat losses from an icy surface on Mars at 237 K and current pressures would exceed the available net radiative flux at the surface, thus prohibiting ice from melting at low atmospheric pressures. Schumann has developed equations describing the horizontal fluctuations and mean vertical gradients occurring during free convection. Schumann's model was generalized to include convection driven by water vapor gradients and to include the effects of circulation above both aerodynamically smooth and rough surfaces.

Clow, G. D.↗
NASA NTRS1976

A complete model of turbulence

A set of constitutive equations suitable for a priori computation of turbulent shear flows has been developed. Since no properties of a given turbulent flow need be known in advance in order to obtain a solution, the equations comprise a complete model of turbulence. Perturbation analysis shows that the model predicts a composite five-layer structure for an incompressible turbulent boundary layer, viz, a defect layer, a law-of-the-wall layer, a viscous sublayer, a near-surface roughness layer, and a viscous superlayer at the boundary-layer edge. Analysis of the defect layer demonstrates the key improvement of the model over its predecessor, the Saffman-Wilcox two-equation model of turbulence. Examination of model-predicted sublayer structure yields model-parameter boundary conditions appropriate for surfaces with roughness and mass injection. Results of numerical computations of compressible and incompressible equilibrium boundary layers show that, for such flows, the model is as accurate as mixing-length theory. Applications to transitional boundary layers and to nonequilibrium relaxation of a boundary layer passing from a rough to a smooth surface indicate that the model's applicability extends far beyond that of mixing-length theory's.

Wilcox, D. C.↗
NASA NTRS1959

Experimental Investigation of the Heat-Transfer Rate to a Series of 20 deg Cones of Various Surface Finishes at a Mach Number of 4.95

The heat-transfer rates were measured on a series of cones of various surface finishes at a Mach number of 4.95 and Reynolds numbers per foot varying from 20 x 10(exp 6) to 100 x 10(exp 6). The range of surface finish was from a very smooth polish to smooth machining with no polish (65 micro inches rms). Some laminar boundary-layer data were obtained, since transition was not artificially tripped. Emphasis, however, is centered on the turbulent boundary layer. The results indicated that the turbulent heat-transfer rate for the highest roughness tested was only slightly greater than that for the smoothest surface. The laminar-sublayer thickness was calculated to be about half the roughness height for the roughest model at the highest value of unit Reynolds number tested.

Jones, Jim J.↗
NASA NTRS2007

Modeling and Grid Generation of Iced Airfoils

SmaggIce Version 2.0 is a software toolkit for geometric modeling and grid generation for two-dimensional, singleand multi-element, clean and iced airfoils. A previous version of SmaggIce was described in Preparing and Analyzing Iced Airfoils, NASA Tech Briefs, Vol. 28, No. 8 (August 2004), page 32. To recapitulate: Ice shapes make it difficult to generate quality grids around airfoils, yet these grids are essential for predicting ice-induced complex flow. This software efficiently creates high-quality structured grids with tools that are uniquely tailored for various ice shapes. SmaggIce Version 2.0 significantly enhances the previous version primarily by adding the capability to generate grids for multi-element airfoils. This version of the software is an important step in streamlining the aeronautical analysis of ice airfoils using computational fluid dynamics (CFD) tools. The user may prepare the ice shape, define the flow domain, decompose it into blocks, generate grids, modify/divide/merge blocks, and control grid density and smoothness. All these steps may be performed efficiently even for the difficult glaze and rime ice shapes. Providing the means to generate highly controlled grids near rough ice, the software includes the creation of a wrap-around block (called the "viscous sublayer block"), which is a thin, C-type block around the wake line and iced airfoil. For multi-element airfoils, the software makes use of grids that wrap around and fill in the areas between the viscous sub-layer blocks for all elements that make up the airfoil. A scripting feature records the history of interactive steps, which can be edited and replayed later to produce other grids. Using this version of SmaggIce, ice shape handling and grid generation can become a practical engineering process, rather than a laborious research effort.

Vickerman, Mary B.↗
NASA NTRS2006

SmaggIce Version 1.8

SmaggIce version 1.8 is a set of software tools for geometrical modeling of, and generation of grids that conform to, both clean and iced airfoils. Ice shapes, especially those that include rough surfaces, pose difficulty in generating high-quality grids that are essential for predicting airflows by use of computational fluid dynamics. SmaggIce version 1.8 contains software tools needed to overcome this difficulty. For a given airfoil, it allows the user to define the flow domain, decompose the domain into blocks, generate grids, merge gridded blocks, and control the density and smoothness of each grid. Among the unique features of version 1.8 is a thin C-shaped block, called a "viscous sublayer block," which is wrapped around an iced airfoil and its wake line and serves as a means to generate highly controlled grids near the rough ice surface. Users can modify block boundary shapes using control points of non-uniform rational B-spline (NURBS) curves. Concave ice regions can be smoothed during geometrical modeling or creation of the viscous sublayer block.

Vickerman, Mary B.↗
NASA NTRS1991

Free convection in the Martian atmosphere

Researchers investigated the free convective regime for the Martian atmospheric boundary layer (ABL). Researchers generalized Schumann's model describing horizontal fluctuations and mean vertical gradients occurring during free convection to include convection driven by water vapor gradients and to include the effects of circulation above both aerodynamically smooth and rough surfaces. Applying the model to Mars, researchers found that nearly all the resistance to sensible and latent heat transfer in the ABL occurs within the thin interfacial sublayer at the surface. Free convection is found to readily occur at low pressures and high temperatures when surface ice is present. At 7 mb, the ABL should freely convect whenever the mean windspeed at the top of the surface layer drops below about 2.5 m s(-1) and surface temperatures exceed 250 K. Mean horizontal fluctuations within the surface layer are found to be as high as 3 m (-1) for windspeed, 0.5 K for temperature, and 10 (-4) kg m (-3) for water vapor density. Airflow over surfaces similar to the Antarctic Polar Plateau was found to be aerodynamically smooth on Mars during free convection for all pressures between 6 and 1000 mb, while surfaces with z sub o approx. equals 1 cm are aerodynamically rough over this pressure range.

Clow, G. D.↗
NASA NTRS1999

A Generalized Wall Function

The asymptotic solutions, described by Tennekes and Lumley (1972), for surface flows in a channel, pipe or boundary layer at large Reynolds numbers are revisited. These solutions can be extended to more complex flows such as the flows with various pressure gradients, zero wall stress and rough surfaces, etc. In computational fluid dynamics (CFD), these solutions can be used as the boundary conditions to bridge the near-wall region of turbulent flows so that there is no need to have the fine grids near the wall unless the near-wall flow structures are required to resolve. These solutions are referred to as the wall functions. Furthermore, a generalized and unified law of the wall which is valid for whole surface layer (including viscous sublayer, buffer layer and inertial sublayer) is analytically constructed. The generalized law of the wall shows that the effect of both adverse and favorable pressure gradients on the surface flow is very significant. Such as unified wall function will be useful not only in deriving analytic expressions for surface flow properties but also bringing a great convenience for CFD methods to place accurate boundary conditions at any location away from the wall. The extended wall functions introduced in this paper can be used for complex flows with acceleration, deceleration, separation, recirculation and rough surfaces.

Shih, Tsan-Hsing↗
NASA NTRS1988

The effects of longitudinal roughness elements upon the turbulent boundary layer

Velocity profiles and bursting statistics were measured in the presence of cylindrical longitudinal roughness elements (LREs) and compared to similar work by Johansen and Smith (1986). An algorithm was devised to detect low speed streaks (LSS) from the hot-wire rake data, allowing an estimation of the LREs' ability to reduce spatial randomness of the sublayer structure. The velocity profiles agree with Johansen and Smith's result that the effects of the LREs are felt only for y(+) less than about 15. The analysis of the bursting statistics and LSS algorithm output shows that the flow relaxes back to a flat plate boundary layer between the LREs. A crude first approximation is that the LREs simply impose the no-slip boundary condition at an elevation equivalent to their diameter. While the LREs might serve as a nucleation site for the LSS, spatial randomness of the structure is affected only very close to the LREs.

Blackwelder, R. F.↗
NASA NTRS1989

Roughening Surfaces Of Solar Cells

Proposed treatment involving bombardment by ions gives silicon solar cells rough surfaces and increases amount of light absorbed by cells. First step of treatment, wafer of single-crystal silicon held at temperature of 70 degree C or higher while bombarded with argon or hydrogen ions at kinetic energies between 100 and 2,000 eV. This high dose produces damaged layer on surface, consisting of outer polycrystalline sublayer and underlying sublayer containing extended defects. Bombarded wafer then etched chemically. Photovoltaic conversion efficiency increases with number of reflections at surface of silicon cell because large part of light absorbed at each reflection. Deeply textured surface provides opportunities for multiple reflections.

Singh, Ranbir↗