Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY LAYER - PLATES, FLAT - LAMINAR”

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 73 records · Page 4

Hubble Space Telescope observations of the dwarf Nova Z Chamaeleontis through two eruption cycles

We have obtained the first high-speed photometry of the eclipsing dwarf nova Z Cha at ultraviolet wavelengths with the Hubble Space Telescope (HST). We observed the eclipse roughly every 4 days over two cycles of the normal eruptions of Z Cha, giving a uniquely complete coverage of its outburst cycle. The accretion disk dominated the ultraviolet light curve of Z Cha at the peak of an eruption; the white dwarf, the bright spot on the edge of the disk, and the boundary layer were all invisible. We were able to obtain an axisymmetric map of the accretion disk at this time only by adopting a flared disk with an opening angle of approximately 8 deg. The run of brightness temperature with radius in the disk at the peak of the eruption was too flat to be consistent with a steady state, optically thick accretion disk. The local rate of mass flow through the disk was approximately 5 x 10(exp -10) solar masses/yr near the center of the disk and approximately 5 x 10(exp -9) solar masses/yr near the outer edge. The white dwarf, the accretion disk, and the boundary layer were all significant contributors to the ultraviolet flux on the descending branches of the eruptions. The temperature of the white dwarf during decline was 18,300 K less than T(sub wd) less than 21,800 K, which is significantly greater than at minimum light. Six days after the maximum of an eruption Z Cha has faded to near minimum light at ultraviolet wavelenghts, but was still approximately 70% brighter at minimum light in the B band. About one-quarter of the excess flux in the B band came from the accretion disk. Thus, the accretion disk faded and became invisible at ultraviolet wavelengths before it faded at optical wavelenghts. The disk did, however, remain optically thick and obscured the lower half of the white dwarf at ultraviolet and possibly at optical wavelenghts for 2 weeks after the eruption ended. By the third week after eruptiuons the eclipse looked like a simple occultation of an unobscured, spherical white dwarf by a dark secondary star. The center of the accretion disk was, therfore, optically thin at ultraviolet wavelenghts and the boundary layer was too faint to be visible.

Robinson, E. L.↗

Pressure-fluctuation inputs and response of panels underlying attached and separated supersonic turbulent boundary layers

Results are presented for an investigation of surface pressure fluctuations and response of panels underlying attached and separated turbulent boundary layers and shock waves. Tests conducted at transonic and supersonic Mach numbers to 3.6 to study the pressure fields. Assorted fixed edge flat panels were tested at Mach numbers from 1.6 to 3.6 in attached and completely separated flow fields, and also in mixed flow with a step-induced shock wave oscillating on the panels. The surface pressure fluctuations are described in terms of broadband rms, spectral density, and spatial correlation information. The effectiveness of parameters for scaling the pressure fluctuations is also illustrated. Measurements of the amplitude and strain response of the panels are compared with response computations by the normal mode method of analysis.

Coe, C. F.↗

Effects of Vegetation and Topography on the Boundary Layer Structure above the Amazon Forest

Observational data from two field campaigns in the Amazon forest were used to study the vertical structure of turbulence above the forest. The analysis was performed using the reduced turbulent kinetic energy (TKE) budget and its associated two-dimensional phase space. Results revealed the existence of two regions within the roughness sublayer in which the TKE budget cannot be explained by the canonical flat-terrain TKE budgets in the canopy roughness sublayer or in the lower portion of the convective ABL. Data analysis also suggested that deviations from horizontal homogeneity have a large contribution to the TKE budget. Results from LES of a model canopy over idealized topography presented similar features, leading to the conclusion that flow distortions caused by topography are responsible for the observed features in the TKE budget. These results support the conclusion that the boundary layer above the Amazon forest is strongly impacted by the gentle topography underneath.

54 ENVIRONMENTAL SCIENCES↗

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of idealized large-eddy simulations performed with the Weather Research and Forecasting model using turbulence measurements from a 250 m meteorological mast

Abstract. We investigate the ability of the Weather Research and Forecasting model to perform large-eddy simulation of canonical flows. This is achieved through comparison of the simulation outputs with measurements from sonic anemometers on a 250 m meteorological mast located at Østerild, in northern Denmark. Østerild is on a flat and rough area, and for the predominant wind directions, the atmospheric flow can be considered to be close to homogeneous. The idealized simulated flows aim at representing atmospheric boundary layer turbulence under unstable, neutral, and stable stability conditions at the surface, which are statistically significant conditions observed at Østerild. We found that the resolved fields from the simulations appear to have the characteristics of the three stability regimes. Vertical profiles of observed mean wind speeds and direction are well reproduced by the simulations, with the largest differences under near-neutral conditions, where the effect of the subgrid-scale model is evident on the vertical wind shear close to the surface. Vertical profiles of observed eddy fluxes are also well reproduced by the simulations, with the largest differences for the three velocity component variances under stable stability conditions, although nearly always within the observed variability. With regards to turbulent kinetic energy, we find good agreement between observations and simulations at all vertical levels. Simulated and observed velocity spectra match very well and show very similar behavior with height and with atmospheric stability within the low-frequency interval; at the effective resolution, the simulated spectra show the typical drop-off of finite differences. Our findings demonstrate that these idealized simulations reproduce the characteristics of atmospheric stability regimes often observed at a high turbulent and flat site within a direction sector, where the air flows over nearly homogeneous land.

17 WIND ENERGY↗

Lessons learned in coupling atmospheric models across scales for onshore and offshore wind energy

Abstract. The Mesoscale to Microscale Coupling team, part of the U.S. Department of Energy Atmosphere to Electrons (A2e) initiative, has studied various important challenges related to coupling mesoscale models to microscale models for the use case of wind energy development and operation. Several coupling methods and techniques for generating turbulence at the microscale that is subgrid to the mesoscale have been evaluated for a variety of cases. Case studies included flat-terrain, complex-terrain, and offshore environments. Methods were developed to bridge the terra incognita, which scales from about 100 m through the depth of the boundary layer. The team used wind-relevant metrics and archived code, case information, and assessment tools and is making those widely available. Lessons learned and discerned best practices are described in the context of the cases studied for the purpose of enabling further deployment of wind energy.

17 WIND ENERGY↗

Stagnation-point ablation of carbonaceous flat disks. I Theory

The process of ablation is calculated for the stagnation region of a flat disk in a radiation-dominated, massive-blowing environment produced in a ballistic range filled with argon. Flow environments are determined by solving the boundary-layer equations while radiative transfer is calculated through a line-by-line spectral computation. The resulting wall heat-transfer rates are coupled with an existing material's response code to determine surface recession and char thickness. The calculation is performed for six 5-cm-diam models made of carbon-phenolic and carbon-carbon composite launched in the Track-G facility at the Arnold Engineering Development Center. Significant surface recessions are predicted to occur for these models due mostly to radiative heating.

Park, C.↗

Heat and mass transfer rates during flow of dissociated hydrogen gas over graphite surface

To improve upon the performance of chemical rockets, the nuclear reactor has been applied to a rocket propulsion system using hydrogen gas as working fluid and a graphite-composite forming a part of the structure. Under the boundary layer approximation, theoretical predictions of skin friction coefficient, surface heat transfer rate and surface regression rate have been made for laminar/turbulent dissociated hydrogen gas flowing over a flat graphite surface. The external stream is assumed to be frozen. The analysis is restricted to Mach numbers low enough to deal with the situation of only surface-reaction between hydrogen and graphite. Empirical correlations of displacement thickness, local skin friction coefficient, local Nusselt number and local non-dimensional heat transfer rate have been obtained. The magnitude of the surface regression rate is found low enough to ensure the use of graphite as a linear or a component of the system over an extended period without loss of performance.

Nema, V. K.↗

Lidar observations of the fine-scale variability of marine stratocumulus clouds

An Nd:YAG lidar system was flown aboard NASA's ER-2 high altitude aircraft. Observations of cloud top height were made with 70 m along-track and 7.5 m vertical-height resolution. The lidar data observed from an East Pacific stratocumulus cloud height deck revealed large cloud variability on 1-5 km scales. The cloud deck sloped upward from 700 to 1000 m in a northeast-southwest direction over a distance of 120 km. Vertical cloud top distributions were negatively skewed indicating flat-topped clouds. The dominant spectral peak of the cloud top variations was found at 4.5 km, which is 5 to 7 times the depth of the local boundary layer. No other peaks were significant in the average spectrum. The cloud layer was stable with respect to cloud top entrainment instability. The southwestern region of the study area was more prone to shear instability at cloud top than the northeastern region. The results of this study show that a lidar system is ideal to provide the topography of clouds and local boundary layer depth. This information is useful in the study of cloud top radiation and parameterization of clouds in numerical models.

Boers, Reinout↗

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

A mixing-length model for the prediction of convex curvature effects on turbulent boundary layers

A mixing-length model is developed for the prediction of turbulent boundary layers with convex streamwise curvature. For large layer thickness ratio, delta/R greater than 0.05, the model scales mixing length on the wall radius of curvature, R. For small delta/R, ordinary flat wall modeling is used for the mixing-length profile with curvature corrections, following the recommendations of Eide and Johnston (1976). Effects of streamwise change of curvature are considered; a strong lag from equilibrium is required when R increases downstream. Fifteen separate data sets were compared, including both hydrodynamic and heat transfer results. Six of these computations are presented and compared to experiment.

Adams, E. W.↗

Interference of electrical double layers: Confinement effects on structure, dynamics, and screening of ionic liquids

Ionic liquids are widely used as electrolytes in electronic devices in which they are subject to nanoconfinement within nanopores or nanofilms. Because the intrinsic width of an electrical double layer is on the order of several nanometers, nanoconfinement is expected to fundamentally alter the double layer properties. Furthermore, in confined systems, a large portion of the ions are interfacial, e.g., at the electrode interface, leading to significant deviations of electrostatic screening and ion dynamics as compared to bulk properties. In this article, we systematically investigate the interference between electrical double layers for nanoconfined ionic liquids and the resulting influence on the structure, dynamics, and screening behavior. We perform molecular dynamics simulations for the ionic liquids [BMIm + ][BF$^-_4$] and [BMIm + ][PF$^-_6$] confined between two flat electrodes at systematic separation distances between 1.5 nm and 4.5 nm for both conducting and insulating boundary conditions. We find that while ion dynamics is expectedly slower than in the bulk (by ~2 orders of magnitude), there is an unexpected non-linear trend with the confinement length that leads to a local maximum in dynamic rates at ~3.5–4.5 nm confinement. We show that this nonlinear trend is due to the ion correlation that arises from the interference between opposite double layers. We further evaluate confinement effects on the ion structure and capacitance and investigate the influence of electronic polarization of the ionic liquid on the resulting properties. This systematic evaluation of the connection between electrostatic screening and structure and dynamics of ionic liquids in confined systems is important for the fundamental understanding of electrochemical supercapacitors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

Evaluation of Atmospheric Boundary Layer Height From Wind Profiling Radar and Slab Models and Its Responses to Seasonality of Land Cover, Subsidence, and Advection

In this study, we evaluated the effect of land cover, atmospheric subsidence, and advection on the annual dynamics of atmospheric boundary layer (ABL) height from two contrasting sites. The first site is the Walker Branch forest, a deciduous forest of temperate climate, complex topography, and cloudy summers. The second site is the Sacramento-San Joaquin River Delta, a site of Mediterranean climate, flat terrain on a local scale, and clear summers. After testing a new algorithm to calculate ABL heights from 915 MHz radar wind profilers, we evaluated a hierarchy of three slab models to recreate the diurnal and annual patterns of ABL growth. We found that the lower ABL heights in the Delta, particularly during late summer, are driven by the combined effects of increased atmospheric subsidence and marine air advection. In both sites, the annual pattern of ABL height was strongly correlated to total daily incoming radiation, and in the Delta, the annual pattern of ABL height closely followed the seasonal patterns of sensible heat flux from a mosaic of different land covers. A land composite of latent and sensible heat fluxes obtained through a meso-network of eddy covariance measurements and the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission resulted in higher model skill, thus showing that land cover heterogeneity is an important driver of ABL growth. To conclude, model simulations show that in the Delta, restoring agricultural land to wetlands with large open water areas could result in a reduction of ABL height during those months with low subsidence and advection.

54 ENVIRONMENTAL SCIENCES↗

Quantifying wind plant blockage under stable atmospheric conditions

Wind plant blockage reduces the wind velocity upstream undermining turbine performance for the first row of the plant. We assess how atmospheric stability modifies the induction zone of a wind plant in flat terrain. We also explore different approaches to quantifying the magnitude and extent of the induction zone from field-like observations. To investigate the influence from atmospheric stability, we compare simulations of two stable boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. We find a faster cooling rate at the surface, which produces a stronger stably stratified boundary layer, amplifies the induction zone of both an isolated turbine and of a large wind plant. A statistical analysis on the hub-height wind speed field shows wind slowdowns only extend far upstream (up to 15D) of a wind plant in strong stable boundary layers. To evaluate different ways of measuring wind plant blockage from field-like observations, we consider various ways of estimating the freestream velocity upstream of the plant. Sampling a large area upstream is the most accurate approach to estimating the freestream conditions, and thus of measuring the blockage effect. Also, the choice of sampling method may induce errors of the same order as the velocity deficit in the induction zone.

54 ENVIRONMENTAL SCIENCES↗

LES simulations of a vacuum membrane distillation channel with geometric alterations

3D LES simulations were carried out to study the performance of a vacuum membrane distillation module. A wiggly wall profile with/without embedded stiffeners was considered to alleviate axial and radial temperature polarization, the cause of performance decrease in membrane distillation. Results of the flow field show the wiggles create unsteady vortex shedding inducing intense mixing in the channel. Vortex shedding intensity increases as the Reynolds number increases or stiffeners are added. 3D results show a 56% improvement in flux when moving from a flat sheet membrane to a wiggly membrane with stiffeners at a constant Reynolds number, corresponding to an improvement from 11.7 to 41.7 in the Nusselt number, showing that the alterations improved the flux performance by enhancing the heat transfer along the membrane surface and therefore alleviating temperature polarization, a critical bottleneck in membrane distillation systems. Further, a merit criterion was defined based on the Nusselt number and friction factor, and a 40% increase in merit was shown switching from a flat to a wiggly channel, while a 97% merit increase was seen going from flat to wiggly with stiffeners. A 42% enhancement in the flux was also seen moving from a straight channel to a wiggly channel at a higher Reynolds number which highlights the importance of strategically choosing the mass flow rate, as well as inducing flow separation and vortex shedding in the channel to promote mixing and dissipate the thermal boundary layer. The variation in results between the wiggly module with stiffeners and the wiggly module at high Reynolds numbers suggests enhancing mixing structures can be more impactful on flux for this geometry than increasing the flow rate to a turbulent/transitional regime. However, both are preferable for peak system flux performance. Furthermore, a 2D approximation was used to perform simulations on more extended channels to examine the length degradation. The modules with wiggly channels performed at the same flux level with a doubling of the length. On the other hand, the flat sheet modules experienced length degradation by temperature polarization and dropped in flux yield by around 13% of the short-channel value. This work illustrates that modeling the system and understanding how the performance decreases as the membrane surface area increases are critical for a larger module (scaling up from a lab to a prototype module) to maintain high flux performance.

42 ENGINEERING↗

Comparison of multiequation turbulence models for several shock-separated boundary-layer interaction flows

Several multiequation eddy viscosity models of turbulence are used with the Navier-Stokes equations to compute three classes of experimentally documented shock-separated turbulent boundary-layer flows. The types of flow studied are: (1) a normal shock at transonic speeds in both a circular duct and a two-dimensional channel; (2) an incident oblique shock at supersonic speeds on a flat surface; and (3) a two-dimensional compression corner at supersonic speeds. Established zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence models are each utilized to describe the Reynolds shear stress for the three classes of flows. These models are assessed by comparing the calculated values of skin friction, wall pressure distribution, velocity, Mach number, and turbulent kinetic energy profiles with experimental measurements. Of the models tested the two-equation model results gave the best overall agreement with the data.

Viegas, J. R.↗