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.

34 records · Page 2
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2022

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗
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.↗
DOE OSTI2023

Recent Developments in Revealing the Impact of Complex Metal Oxide Reconstruction on Catalysis

Surface reconstruction of complex metal oxides refers to the alterations of the topmost catalyst surface, relative to the bulk structure, following exposure to pretreatment or reaction conditions during heterogeneous catalysis. (1) This structural and electronic transformation may involve (1) changes of the exposed crystallographic facets, (2) changes in chemical composition, such as the enrichment of specific elements, and, (3) morphological variations, including surface roughness and defects (Scheme 1). Beyond relaxation, surface reconstruction here also encompasses disruptions to the periodicity of the sublayers and changes to its stoichiometry. This change significantly influences the chemical properties of the catalyst surface, such as redox sites, (2) acid/base pairs, hydroxyl groups, (3) and surface defects. Surface reconstruction occurs due to the thermodynamic drive to lower the surface energy of the catalyst under the surrounding environment (e.g., temperature, chemical potential of species present). In addition to the treatment or reaction conditions, the extent of reconstruction also depends on the elements in the metal oxide structure. For instance, annealing ABO 3 perovskites (SrTiO 3 , BaTiO 3 , and BaZrO 3 ) at high temperatures (500 °C) in O 2 has a very different impact on the A/B ratio at the surface (ranging from 0.9 to 2.5), depending on the identity of A and B. (4) Metal oxide surface reconstruction under an electric field and exposure to acidic and alkaline media has recently gained increased attention in electrocatalysis, such as in the oxygen evolution reaction (OER); (5,6) however, in this Viewpoint we focus on the impact of surface reconstruction on thermal catalysis, a critical enduring research theme with ongoing challenges. Here, we start showcasing approaches to characterize surface reconstruction at the topmost surface layer, such as low-energy ion scattering (LEIS), discussing advantages and limitations. Next, we connect characterization of surface reconstruction to site-specific kinetic analysis, enabled via steady-state isotopic transient kinetic analysis (SSITKA), and we express our viewpoint on building structure–performance relationships by deeply understanding reconstructed surfaces. Later, we comment on the potential of tuning surface reconstruction to enhance catalytic performance. Finally, we explore the interplay between surface reconstruction and the “intelligent behavior”.

alcohol conversion↗
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↗
NASA NTRS1993

A comparison of airborne eddy correlation and bulk aerodynamic methods for ocean-air turbulent fluxes during cold-air outbreaks

The viscous interfacial-sublayer model of Liu et al. (1979) is used to derive four bulk schemes (LKB, FG, D, and DB), with the flux-profile relationships of Lie et al., Francey and Garratt (1981), Dyer (1974), and Dyer and Bradley (1982). These schemes, with stability-dependent transfer coefficients, are tested against the eddy-correlation fluxes measured at the 50-m flight level above the western Atlantic Ocean during cold-air outbreaks. The bulk fluxes of momentum (tau), sensible heat (H), and latent heat (E) are found to increase with various von Karman constants. The dependence of transfer coefficients on wind speeds and roughness lengths is discussed. The transfer coefficients for tau and E agree excellently between LKB and FG. The ratio of the coefficent for H of LKB to that of FG, increasing with decreasing stability, is very sensitive to stability at low winds, but approaches the neutral value of 1.25 at high winds.

Chou, Shu-Hsien↗