Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “evolving boundary layer”

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

Effects of Cavities and Protuberances on Transition over Hypersonic Vehicles

Surface protuberances and cavities on a hypersonic vehicle are known to cause several aerodynamic or aerothermodynamic issues. Most important of all, premature transition due to these surface irregularities can lead to a significant rise in surface heating. To help understand laminar-turbulent transition induced by protuberances or cavities on a Crew Exploration Vehicle (CEV) surface, high-fidelity numerical simulations are carried out for both types of trips on a CEV wind tunnel model. Due to the large bluntness, these surface irregularities reside in an accelerating subsonic boundary layer. For the Mach 6 wind tunnel conditions with a roughness Reynolds number Re(sub kk) of 800, it was found that a protuberance with a height to boundary layer thickness ratio of 0.73 leads to strong wake instability and spontaneous vortex shedding, while a cavity with identical geometry only causes a rather weak flow unsteadiness. The same cavity with a larger Reynolds number also leads to similar spontaneous vortex shedding and wake instability. The wake development and the formation of hairpin vortices for both protuberance and cavity were found to be qualitatively similar to that observed for an isolated hemisphere submerged in a subsonic, low speed flat-plate boundary layer. However, the shed vortices and their accompanying instability waves were found to be slightly stabilized downstream by the accelerating boundary layer along the CEV surface. Despite this stabilizing influence, it was found that the wake instability spreads substantially in both wall-normal and azimuthal directions as the flow is evolving towards a transitional state. Similarities and differences between the wake instability behind a protuberance and a cavity are investigated. Computations for the Mach 6 boundary layer over a slender cylindrical roughness element with a height to the boundary layer thickness of about 1.1 also shows spontaneous vortex shedding and strong wake instability. Comparisons of detailed flow structures associated with protuberances at subsonic and supersonic edge Mach numbers indicate distinctively different instability mechanisms.

Chang, Chau-Lyan↗

Dependence of Cloud Macrophysical Properties and Phase Distributions on Environmental Conditions Over the North Atlantic and Southern Ocean: Results From COMBLE and MARCUS

The accurate representation of Cold Air Outbreaks (CAOs) and affiliated mixed-phase boundary layer (BL) clouds in models is challenging. How BL cloud properties evolve during CAOs and their dependence on meteorological conditions is not well understood but is important for the simulation of Earth's energy budgets. Here the properties of polar BL clouds over the North Atlantic (NA) and Southern Ocean (SO) are compared using observations from the Measurements of Aerosol Radiation and CloUds over the SO (MARCUS) and CAOs in the Marine BL Experiment (COMBLE) conducted over the NA. MARCUS observations show a stronger BL inversion than COMBLE, with a higher mean EIS (estimated inversion strength)/LTS (lower tropospheric stability) of -0.03 K/13 K compared to COMBLE’s -3.2 K/9.3 K. 39% of CAOs observed during COMBLE were intense with M > 5 K, while MARCUS only had 1.3%. 78%/72% of clouds sampled in CAOs during COMBLE/MARCUS had cloud top heights <4 km. The mean BL cloud top height was over 400 m higher, and the BL was over 500 m deeper for M of 10 K compared to 0 K for both regions. MARCUS observed a 27% moister BL structure than COMBLE when M > 5 K due to stronger BL inversion trapping more moisture within the BL. Under the same LTS, EIS, and M conditions, MARCUS observed a 12% drier BL structure, and clouds were 46% more turbulent than COMBLE. During CAOs, 54% of single-layer BL clouds sampled during MARCUS had liquid-dominated bases compared to 39% during COMBLE.

54 ENVIRONMENTAL SCIENCES↗

Drag characteristics of unsteady, perturbed boundary layer flows

A series of time-dependent numerical computations have been performed for the flow of an unsteady boundary layer over an embedded cavity, which is aligned normal to the flow direction. The unsteady flow is the result of the imposition of a Stuart vortex, at the inflow boundary of the computational domain, onto an otherwise unperturbed two-dimensional laminar boundary layer flow. This produces a strongly vortical, but spatially monochromatic, motion at inflow which is allowed to evolve downstream through the time-dependent numerical solution of appropriate conservation equations. As the imposed vortex evolves over the solid boundary, a weaker induced vortical motion is produced and together these perturbing motions interact with the fluid motion in the embedded cavity. An analysis of the pertinent dynamic variables as well as the time-dependent drag characteristics for this type of non-planar geometry is performed; with particular emphasis on the relative contribution of both pressure drag and frictional drag at various times in the evolution process.

Gatski, T. B.↗

Direct Simulation of Evolution and Control of Three-Dimensional Instabilities in Attachment-Line Boundary Layers

The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic- source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in at-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.

Joslin, Ronald D.↗

Parallel spatial direct numerical simulations on the Intel iPSC/860 hypercube

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube is documented. The direct numerical simulation approach is used to compute spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows. The feasibility of using the PSDNS on the hypercube to perform transition studies is examined. The results indicate that the direct numerical simulation approach can effectively be parallelized on a distributed-memory parallel machine. By increasing the number of processors nearly ideal linear speedups are achieved with nonoptimized routines; slower than linear speedups are achieved with optimized (machine dependent library) routines. This slower than linear speedup results because the Fast Fourier Transform (FFT) routine dominates the computational cost and because the routine indicates less than ideal speedups. However with the machine-dependent routines the total computational cost decreases by a factor of 4 to 5 compared with standard FORTRAN routines. The computational cost increases linearly with spanwise wall-normal and streamwise grid refinements. The hypercube with 32 processors was estimated to require approximately twice the amount of Cray supercomputer single processor time to complete a comparable simulation; however it is estimated that a subgrid-scale model which reduces the required number of grid points and becomes a large-eddy simulation (PSLES) would reduce the computational cost and memory requirements by a factor of 10 over the PSDNS. This PSLES implementation would enable transition simulations on the hypercube at a reasonable computational cost.

Joslin, Ronald D.↗

Simulations of bypass transition for spatially evolving disturbances

The spatial evolution of disturbances in plane Poiseuille flow and zero pressure gradient boundary layer flow is considered. For disturbances governed by the linearized equations, potential for significant transient growth of the amplitude is demonstrated. The maximum amplification occurs for disturbances with zero or near zero frequencies. Spatial numerical simulations of the transition scenario involving a pair of oblique waves has been conducted for both flows. A fully spectral solver using a simple but efficient fringe region technique allowed the flows to be computed with high resolution into the fully turbulent domain. A modal decomposition of the simulation results indicates that non-linear excitation of the transient growth is responsible for the rapid emergence of low-frequency structures. Physically, this corresponds to streaky flow structures, as seen from the results of a numerical amplitude expansion. Thus, this spatial transition scenario has been found to be similar to the corresponding temporal one. In the boundary layer simulations the streaks are seen to break down from what appears to be a secondary instability.

Lundbladh, A.↗

Scalability of Parallel Spatial Direct Numerical Simulations on Intel Hypercube and IBM SP1 and SP2

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube and IBM SP1 and SP2 parallel computers is documented. Spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows are computed with the PSDNS code. The feasibility of using the PSDNS to perform transition studies on these computers is examined. The results indicate that PSDNS approach can effectively be parallelized on a distributed-memory parallel machine by remapping the distributed data structure during the course of the calculation. Scalability information is provided to estimate computational costs to match the actual costs relative to changes in the number of grid points. By increasing the number of processors, slower than linear speedups are achieved with optimized (machine-dependent library) routines. This slower than linear speedup results because the computational cost is dominated by FFT routine, which yields less than ideal speedups. By using appropriate compile options and optimized library routines on the SP1, the serial code achieves 52-56 M ops on a single node of the SP1 (45 percent of theoretical peak performance). The actual performance of the PSDNS code on the SP1 is evaluated with a "real world" simulation that consists of 1.7 million grid points. One time step of this simulation is calculated on eight nodes of the SP1 in the same time as required by a Cray Y/MP supercomputer. For the same simulation, 32-nodes of the SP1 and SP2 are required to reach the performance of a Cray C-90. A 32 node SP1 (SP2) configuration is 2.9 (4.6) times faster than a Cray Y/MP for this simulation, while the hypercube is roughly 2 times slower than the Y/MP for this application. KEY WORDS: Spatial direct numerical simulations; incompressible viscous flows; spectral methods; finite differences; parallel computing.

Joslin, Ronald D.↗

Impact of Clouds and Blowing Snow on Surface and Atmospheric Boundary Layer Properties over Dome C, Antarctica

Clouds and blowing snow (BLSN) occur frequently over Antarctica, where it is critical to understand their feedbacks to surface and atmospheric boundary layer processes. Dome C, an elevated East Antarctic station, dominated by lengthy periods of surface longwave (LW) radiative cooling, is selected to reveal cloud and BLSN impacts within a largely stable environment. The sky condition is classified as clear, cloudy, or BLSN, using 3 years of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations satellite data. Co-located and contemporaneous in situ observations are used to investigate the relationship of sky condition with surface and atmospheric boundary layer thermal structure, focusing on seasonal variability. Results show that increased downwelling LW radiation from clouds abate surface radiative cooling losses, contributing to warming during all seasons. An increase of 3°C in the mean surface air temperature is observed during spring, whereas, a more dramatic rise (around 10°C), due to accompanying large-scale subsidence, is observed during fall and winter in association with clouds. For all seasons, the wind speed and wind speed shear are strongest during BLSN events, and the surface-based inversion is weakened by cooling which peaks in a shallow above-surface turbulent layer. The stronger background stability during fall and winter seasons, restricts turbulence and BLSN depths generally to the lowest tens of meters. The Earth's cryosphere is among the most rapidly evolving yet least well-observed regions, and knowledge of clouds and BLSN interactions with the typical stable atmospheric boundary layer can help further understand energy and moisture exchanges.

CALIPSO↗

Factors which influence the development of a low-level jet and coastal cyclogenesis

Mesoscale model simulations were run to examine the mechanisms which generate a low-level jet (LLJ) and the sea-level pressure decrease (SLPD) associated with secondary cyclogenesis along the East Coast of the U.S. Data collected during the Presidents' Day cyclone of February 18-19, 1979 are reviewed, including the behavior of the LLJ preceding cyclogenesis. The simulations covered adiabatic conditions, the absence and presence of latent heating, and the inclusion of all physical parameters with and without computations of boundary layer phenomena, 60-km grid-scale precipitation, and convective precipitation. The results indicate that synergistic reactions among the LLJ, latent heat release, jet-induced circulation, and boundary layer processes are necessary to account for secondary cyclogenesis and the accompanying rapidly evolving mass, momentum and moisture fields.

Uccellini, Louis W.↗

Boundary conditions for direct computation of aerodynamic sound generation

A numerical scheme suitable for the computation of both the near field acoustic sources and the far field sound produced by turbulent free shear flows utilizing the Navier-Stokes equations is presented. To produce stable numerical schemes in the presence of shear, damping terms must be added to the boundary conditions. The numerical technique and boundary conditions are found to give stable results for computations of spatially evolving mixing layers.

Colonius, Tim↗

Stereo image visualization of numerically simulated turbulence

Stereo imaging techniques are used in the current study to analyze the results of a direct Navier-Stokes simulation of a turbulent boundary layer. Coherent structures embedded within the unsteady, highly three-dimensional turbulent flow fields are significantly enhanced through the use of stereo computer graphic renderings. In the paper, the geometrical foundations of two different methods for generating stereo images are reviewed. Sufficient detail is provided to enable readers to code stereo display algorithms from scratch. The effects of varying the geometric parameters of the stereo algorithm are demonstrated. Sample image pairs from animations of the turbulence numerical simulation are presented in color. Stereo imagery, especially when used with time-evolving data, provides new insight into the physical processes responsible for the generation an maintenance of turbulence within boundary layers.

Robinson, Stephen K.↗

Water-tunnel experiments on an oscillating airfoil at RE equals 21,000

Flow visualization experiments were performed in a water tunnel on a modified NACA 0012 airfoil undergoing large amplitude harmonic oscillations in pitch. Hydrogen bubbles were used to: (1) create a conveniently striated and well preserved set of inviscid flow markers; and (2) to expose the succession of events occurring within the viscous domain during the onset of dynamic stall. Unsteady effects were shown to have an important influence on the progression of flow reversal along the airfoil surface prior to stall. A region of reversed flow underlying a free shear layer was found to momentarily exist over the entire upper surface without any appreciable disturbance of the viscous-inviscid boundary. A flow protuberance was observed to develop near the leading edge, while minor vortices evolve from an expanding instability of the free shear layer over the rear portion of the airfoil. The complete breakdown of this shear layer culminates in the successive formation of two dominant vortices.

Mcalister, K. W.↗

Linear and nonlinear stability of the Blasius boundary layer

Two new techniques for the study of the linear and nonlinear instability in growing boundary layers are presented. The first technique employs partial differential equations of parabolic type exploiting the slow change of the mean flow, disturbance velocity profiles, wavelengths, and growth rates in the streamwise direction. The second technique solves the Navier-Stokes equation for spatially evolving disturbances using buffer zones adjacent to the inflow and outflow boundaries. Results of both techniques are in excellent agreement. The linear and nonlinear development of Tollmien-Schlichting (TS) waves in the Blasius boundary layer is investigated with both techniques and with a local procedure based on a system of ordinary differential equations. The results are compared with previous work and the effects of non-parallelism and nonlinearity are clarified. The effect of nonparallelism is confirmed to be weak and, consequently, not responsible for the discrepancies between measurements and theoretical results for parallel flow.

Bertolotti, F. P.↗

Atomic Structure of the Initial Nucleation Layer in Hexagonal Perovskite BaRuO 3 Thin Films

Hexagonal perovskites are an attractive family of materials due to their various polymorphs and rich structure–property relationships. BaRuO 3 (BRO) is a prototypical hexagonal perovskite, in which the electromagnetic properties are significantly modified depending on its lattice structure. Whereas thin-film synthesis would allow the study of the diverse properties of hexagonal perovskites by epitaxially stabilizing various metastable polymorphs, the epitaxial growth of hexagonal perovskites, especially at the initial growth stage, is yet to be addressed. In this study, it is shown that an intriguing nucleation behavior takes place during the initial stabilization of a hexagonal 9R BRO thin film deposited on a (111) SrTiO 3 (STO) substrate. It is found that a nanoscale strained layer composed of different RuO 6 octahedral stacking is initially formed at the interface followed by a relaxed single crystal 9R BRO thin film. Within the interface layer, hexagonal BRO is nucleated on the (111) STO substrate by both face- and corner-sharing octahedra. More interestingly, it is found that boundaries between differently stacked nucleation layers (heterostructural boundaries) facilitate strain relaxation, in addition to the formation of conventional misfit dislocations evolved from homostructural boundaries. The observations reveal an important underlying mechanism to understand the thin-film epitaxy and strain accommodation in hexagonal perovskites.

36 MATERIALS SCIENCE↗

An experimental investigation of a low Reynolds number turbulent boundary layer subject to an adverse pressure gradient

The evolution of a low Re(sub theta) turbulent boundary layer in an adverse pressure gradient (APG) is being studied for comparison with direct numerical simulations by Spalart. A short region of favorable pressure gradient (FPG) is applied first to establish a self preserving layer with Re(sub theta) equals approx. 600, which is a suitable initial condition for simulations. The APG is then applied rapidly such that beta equals approx. 2 at Re(sub theta) equals approx. 1500. The streamwise extent of the measurements exceeds the current capabilities of direct simulation so that the results should also serve as a useful data base for Reynolds averaged boundary layer prediction methods and in the future for direct simulation schemes as computer technology evolves.

Watmuff, J. H.↗

Coherent structures in a simulated turbulent mixing layer

A direct numerical simulation of a plane turbulent mixing layer has been performed. The simulation was initialized using two turbulent velocity fields obtained from direct numerical simulation of a turbulent boundary layer at momentum thickness Reynolds number 300 (Spalart, 1988). The mixing layer is allowed to evolve long enough for self-similar linear growth to occur, with the visual thickness Reynolds number reaching 14,000. The simulated flow is examined for evidence of the coherent structures expected in a mixing layer (rollers and rib vortices). Before the onset of self-similar growth, such structures are present with properties similar to the corresponding laminar or transitional structures. In the self-similar growth regime, however, only the rollers are present with no indication of rib vortices and no indication of conventional pairing. This results in a reduction of mixing and layer growth.

Moser, Robert D.↗

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗

Stoichiometric effects on grain growth in zirconium carbide coatings for high-temperature nuclear fuel

Interest in coated particle fuel for space nuclear propulsion (SNP) has expanded in recent years due to successful demonstrations of the resiliency of the coatings to extreme environments. For SNP applications, the coating layer for the particle design needs to be able to withstand exposure to high temperature hydrogen during operating conditions. ZrC has been proposed as a protective layer, however, it is important to understand the high temperature behavior to ensure adequate protection to this fuel. In this study, surrogate ZrC coated particles were heat treated at 1900 °C up to 300 min, to examine how the microstructure evolves when exposed to high temperature. Scanning electron microscopy and electron backscatter diffraction (EBSD) were conducted to determine grain size and grain boundary character and orientation to determine the degree of change in the ZrC layer post heat treatment. Raman spectroscopy provided insight to understand how the as-fabricated stoichiometry of each sample contributed to the differences in grain growth behavior. Despite the as-fabricated samples showing a similar initial grain size and grain boundary character, the samples with a higher amount of excess carbon exhibited smaller grain areas and slower growth rates when exposed to 1900 °C. In conclusion, this investigation details the as-fabricated microstructure of the ZrC layer, specifically grain size, evolved under high temperature as this can impact the performance of the fuel under operating conditions.

EBSD↗