SEARCH · Engineering Papers

Results for “Shear Layers”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗

DOE OSTI2020

Effects of Tethered Polymers on Dynamics of Nanoparticles in Unentangled Polymer Melts

Polymer-tethered nanoparticles (NPs) are commonly added to a polymer matrix to improve the material properties. Critical to the fabrication and processing of such composites is the mobility of the tethered NPs. Here, we study the motion of tethered NPs in unentangled polymer melts using molecular dynamics simulations, which offer a precise control of the grafted chain length N g and the number z of grafted chains per particle. As Ng increases, there is a crossover from particle-dominated to tethered-chain-dominated terminal diffusion of NPs with the same z. The mean squared displacement of loosely tethered NPs in the case of tethered-chain-dominated terminal diffusion exhibits two subdiffusive regimes at intermediate time scales for small z. The first one at shorter time scales arises from the dynamical coupling of the particle and matrix chains, while the one at longer time scales is due to the participation of the particle in the dynamics of the tethered chains. The friction of loosely grafted chains in unentangled melts scales linearly with the total number of monomers in the chains, as the friction of individual monomers is additive in the absence of hydrodynamic coupling. As more chains are grafted to a particle, hydrodynamic interactions between grafted chains emerge. As a result, there is a nondraining layer of hydrodynamically coupled chain segments surrounding the bare particle. Outside the nondraining layer is a free-draining layer of grafted chain segments with no hydrodynamic coupling. The boundary of the two layers is the stick surface where the shear stress due to the relative melt flow is balanced by the friction between the grafted and melt chains in the interpenetration layer. The stick surface is located further away from the bare surface of the particle with higher grafting density.

42 ENGINEERING↗

DOE OSTI2025

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗

DOE OSTI2022

The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

Distortion of a flat-plate boundary layer by free-stream vorticity normal to the plate

A nominally uniform flow over a semiinfinite flat plate is considered. The analysis shows how a small streamwise disturbance in the otherwise uniform flow ahead of the plate is amplified by leading-edge bluntness effects and eventually leads to a small-amplitude but nonlinear spanwise motion far downstream from the leading edge of the plate. This spanwise motion is then imposed on the viscous boundary-layer flow at the surface of the plate - causing an order-one change in its profile shape. This ultimately reduces the wall shear stress to zero, causing the boundary layer to undergo a localized separation, which may be characterized as a kind of bursting phenomenon that could be related to the turbulent bursts observed in some flat-plate boundary-layer experiments.

Goldstein, M. E.↗

NASA NTRS1983

A miniature, directional surface-fence gage for three-dimensional turbulent boundary layer measurements

A new miniature (3.175 mm overall diameter) surface-fence gage was developed to measure simultaneously the magnitude and direction of wall shear stress under three-dimensional turbulent boundary layers. Two fence elements were placed orthogonal to each other on the end surface of the gage. The gage was calibrated for its directional and magnitude sensitivities. It was then tested under an axisymmetric swirling boundary layer, subject to a sudden transverse strain. The advantage of the gage was demonstrated, and excellent agreement was obtained with the results from surface oil flow interferometry and velocity field measurements. The results also compared well with numerical computations.

Higuchi, H.↗

NASA NTRS1975

Experimental results for the transpired turbulent boundary layer in an adverse pressure gradient

The fluid mechanics of transpired incompressible turbulent boundary layers under zero and adverse pressure gradient conditions is investigated using an open-ended wind tunnel with a porous floor in the test section and a secondary air system for supply and metering of the transpiration air. All velocity profiles and turbulence profiles are obtained by linearized constant-temperature hot-wire anemometry. The wall shear stress is determined by measuring the shear stress away from the wall and extrapolating to the wall by integrating the boundary layer equations for the shear-stress profile. Equilibrium boundary layers are obtained when the transpiration velocity is varied such that the blowing parameter and the Clauser pressure gradient parameter are held constant. The experimental results obtained are presented in tabular and graphical forms.

Andersen, P. S.↗

NASA NTRS1959

The Interaction of an Oblique Shock Wave with a Laminar Boundary Layer

The results of some experimental and theoretical studies of the interaction of oblique shock waves with laminar boundary layers are presented. Detailed measurements of pressure distribution, shear distribution, and velocity profiles were made during the interaction of oblique shock waves with laminar boundary layers on a flat plate. From these measurements a model was derived to predict the pressure levels characteristic of separation and the length of the separated region.

Hakkinen, R. J.↗

DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗

NASA NTRS1973

Contact problem for an elastic reinforcement bonded to an elastic plate

The stiffening layer is treated as an elastic membrane and the base plate is assumed to be an elastic continuum. The bonding between the two materials is assumed to be either one of direct adhesion ro through a thin adhesive layer which is treated as a shear spring. The solution for the simple case in which both the stiffener and the base plate are treated as membranes is also given. The contact stress is obtained for a series of numerical examples. In the direct adhesion case the contact stress becomes infinite at the stiffener ends with a typical square root singularity for the continuum model, and behaving as a delta function for the membrane model. In the case of bonding through an adhesive layer the contact stress becomes finite and continuous along the entire contact area.

Erdogan, F.↗

NASA NTRS1974

Contact problem for an elastic reinforcement bonded to an elastic plate

The contact problem for a thin elastic reinforcement bonded to an elastic plate is considered. The stiffening layer is treated as an elastic membrane and the base plate is assumed to be an elastic continuum. The bonding between the two materials is assumed to be either one of direct adhesion or through a thin adhesive layer which is treated as a shear spring. The solution for the simple case in which both the stiffener and the base plate are treated as membranes is also given. The contact stress is obtained for a series of numerical examples. In the direct adhesion case the contact stress becomes infinite at the stiffener ends with a typical square root singularity for the continuum model and behaving as a delta function for the membrane model. In the case of bonding through an adhesive layer the contact stress becomes finite and continuous along the entire contact area.

Erdogan, F.↗

NASA NTRS1951

Analysis of turbulent free-convection boundary layer on flat plate

With the use of Karman's integrated momentum equation for the boundary layer and data on the wall-shearing stress and heat transfer in forced-convection flow, a calculation was carried out for the flow and heat transfer in the turbulent free-convection boundary layer on a vertical flat plate. The calculation is for a fluid with a Prandtl number that is close to 1. A formula was derived for the heat-transfer coefficient that was in good agreement with experimental data in the range of Grashof numbers from 10sup10 to 10sup12. Because of the good agreement between the theoretical formula and the experimental data, the formula may be used to obtain data for high Grashof numbers. The calculation also yielded formulas for the maximum velocity in the boundary layer and for boundary-layer thickness.

Eckert, E R G↗

NASA NTRS1961

Sound Radiation from a Turbulent Boundary Layer

If the restriction of incompressibility in the turbulence problem is relaxed, the phenomenon of energy radiation in the form of sound from the turbulent zone arises. In order to calculate this radiated energy, it is shown that new statistical quantities, such as time-space correlation tensors, have to be known within the turbulent zone in addition to the conventional quantities. For the particular case of the turbulent boundary layer, indications are that the intensity of radiation becomes significant only in supersonic flows. Under these conditions, the recent work of Phillips is examined together with some experimental findings of the author. It is shown that the qualitative features of the radiation field (intensity, directionality) as predicted by the theory are consistent with the measurements; however, even for the highest Mach number flow, some of the assumptions of the asymptotic theory are not yet satisfied in the experiments. Finally, the question of turbulence damping due to radiation is discussed, with the result that in the Mach number range covered by the experiments, the energy lost from the boundary layer due to radiation is a small percentage of the work done by the wall shearing stresses.

TURBULENT BOUNDARY LAYER↗

NASA NTRS1989

Development of high resolution simulations of the atmospheric environment using the MASS model

Numerical simulations were performed with a very high resolution (7.25 km) version of the MASS model (Version 4.0) in an effort to diagnose the vertical wind shear and static stability structure during the Shuttle Challenger disaster which occurred on 28 January 1986. These meso-beta scale simulations reveal that the strongest vertical wind shears were concentrated in the 200 to 150 mb layer at 1630 GMT, i.e., at about the time of the disaster. These simulated vertical shears were the result of two primary dynamical processes. The juxtaposition of both of these processes produced a shallow (30 mb deep) region of strong vertical wind shear, and hence, low Richardson number values during the launch time period. Comparisons with the Cape Canaveral (XMR) rawinsonde indicates that the high resolution MASS 4.0 simulation more closely emulated nature than did previous simulations of the same event with the GMASS model.

Kaplan, Michael L.↗

NASA NTRS1992

The three-dimensional evolution of a plane mixing layer - The Kelvin-Helmholtz rollup

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber and Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Rogers, Michael M.↗

NASA NTRS1992

Magellan observations of extended impact crater related features on the surface of Venus

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber & Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Campbell, D. B.↗

NASA NTRS1984

Kelvin wave packets and flow acceleration - A comparison of modeling and observations

Atmospheric Kelvin waves, as revealed by temperatures obtained from the recent Limb Infrared Monitor of the Stratosphere (LIMS) experiment, commonly occur in packets. A simple two-dimensional gravity-wave model is used to study the upward propagation of these packets through different zonal mean wind profiles derived from the LIMS data. The observed prevalence of high frequency waves in the lower mesosphere and low frequency waves in the lower stratosphere can be exlained by dispersion of energy associated with the range of frequencies comprising a packet. Dominant wave frequencies at upper and lower levels are more distinctly separated if the packet propagates through a layer of westerly winds. Due to dispersion and shear effects, a packet of short temporal length at low levels will have a considerably extended impact on a layer of westerly winds at higher levels. Observed and modeled westerly accelerations resulting from packet absorption occur in the same layer, and are similar in magnitude and duration. These results support the theory that Kelvin waves are responsible for the westerly phase of the semiannual oscillation.

Coy, L.↗

NASA NTRS1988

Skin friction measurements by laser interferometry in swept shock wave/turbulent boundary-layer interactions

The laser interferometric skin friction meter was used to measure wall shear stress distributions in two interactions of fin-generated swept shock waves with turbulent boundary layers. The basic research configuration was an unswept sharp-leading-edge fin of variable angle mounted on a flatplate. The results indicate that such measurements are practical in high-speed interacting flows, and that a repeatability of + or - 6 percent or better is possible. Marked increases in wall shear were observed in both swept interactions tested.

Kim, Kwang-Soo↗