Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear flow”

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 1,261 records · Page 70

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Development of buried wire gages for measurement of wall shear stress in Blastane experiments

Buried Wire Gages operated from a Constant Temperature Anemometer System are among the special types of instrumentation to be used in the Boundary Layer Apparatus for Subsonic and Transonic flow Affected by Noise Environment (BLASTANE). These Gages are of a new type and need to be adapted for specific applications. Methods were developed to fabricate Gage inserts and mount those in the BLASTANE Instrumentation Plugs. A large number of Gages were prepared and operated from a Constant Temperature Anemometer System to derive some of the calibration constants for application to fluid-flow wall shear-stress measurements. The final stage of the calibration was defined, but could not be accomplished because of non-availability of a suitable flow simulating apparatus. This report provides a description of the Buried Wire Gage technique, an explanation of the method evolved for making proper Gages and the calibration constants, namely Temperature Coefficient of Resistance and Conduction Loss Factor.

Murthy, S. V.↗

Rheological properties of nondilute suspensions of deformable particles

The creeping flow equations have been solved to account for effects of shear on drop shape when the influence of neighboring drops is considered. A cell model is used to account for finite concentrations of the dispersed phase. From these results a constitutive equation is derived and implications for several flows are discussed. The shape and orientation of drops are computed and compared with existing data. Finally, dynamic viscosity and rigidity are obtained for nondilute suspensions of viscoelastic particles in an oscillatory shear flow.-

Choi, S. J.↗

Wall Shear Stress, Wall Pressure and Near Wall Velocity Field Relationships in a Whirling Annular Seal

The mean and phase averaged pressure and wall shear stress distributions were measured on the stator wall of a 50% eccentric annular seal which was whirling in a circular orbit at the same speed as the shaft rotation. The shear stresses were measured using flush mounted hot-film probes. Four different operating conditions were considered consisting of Reynolds numbers of 12,000 and 24,000 and Taylor numbers of 3,300 and 6,600. At each of the operating conditions the axial distribution (from Z/L = -0.2 to 1.2) of the mean pressure, shear stress magnitude, and shear stress direction on the stator wall were measured. Also measured were the phase averaged pressure and shear stress. These data were combined to calculate the force distributions along the seal length. Integration of the force distributions result in the net forces and moments generated by the pressure and shear stresses. The flow field inside the seal operating at a Reynolds number of 24,000 and a Taylor number of 6,600 has been measured using a 3-D laser Doppler anemometer system. Phase averaged wall pressure and wall shear stress are presented along with phase averaged mean velocity and turbulence kinetic energy distributions located 0.16c from the stator wall where c is the seal clearance. The relationships between the velocity, turbulence, wall pressure and wall shear stress are very complex and do not follow simple bulk flow predictions.

Morrison, Gerald L.↗

Alkaline phosphatase in osteoblasts is down-regulated by pulsatile fluid flow

It is our hypothesis that interstitial fluid flow plays a role in the bone remodeling response to mechanical loading. The fluid flow-induced expression of three proteins (collagen, osteopontin, and alkaline phosphatase) involved in bone remodeling was investigated. Rat calvarial osteoblasts subjected to pulsatile fluid flow at an average shear stress of 5 dyne/cm2 showed decreased alkaline phosphatase (AP) mRNA expression after only 1 hour of flow. After 3 hours of flow, AP mRNA levels had decreased to 30% of stationary control levels and remained at this level for an additional 5 hours of flow. Steady flow (4 dyne/cm2 fluid shear stress), in contrast, resulted in a delayed and less dramatic decrease in AP mRNA expression to 63% of control levels after 8 hours of flow. The reduced AP mRNA expression under pulsatile flow conditions was followed by reduced AP enzyme activity after 24 hours. No changes in collagen or osteopontin mRNA expression were detected over 8 hours of pulsatile flow. This is the first time fluid flow has been shown to affect gene expression in osteoblasts.

Non-NASA Center↗

Turbulent boundary layer around a group of obstacles in the direction of flow

Results of an investigation of a boundary layer in a turbulent flow on the surface of a wall having a group of obstacles on the path of flow are presented with regard to the mean velocity field, velocity distribution of the two dimensional flow, wall surface shear stresses and Reynolds stresses measured in a downstream cross section where an interference of boundary layers takes place in a flow around adjacent obstacles arranged on the path of flow.

Nakamura, I.↗

Evaluation of Baldwin-Barth turbulence model with thin-layer Navier-Stokes computation of an axisymmetric afterbody-exhaust jet flowfield

Axisymmetric steady flowfield over a steep circular arc afterbody with a choked exhaust jet in a subsonic freestream is numerically investigated using overflow code, developed by NASA Ames. The accuracy of Baldwin-Barth one-equation turbulence model for the computation of separated flow and plume shear layers is assessed by a comparison of predictions with available measurements from NASA Langley. The results show that the Baldwin-Barth turbulence model performs better than the Baldwin-Lomax algebraic model in predicting the surface pressure distribution in the separated flow region of the afterbody. In the jet mixing layer, the Baldwin-Barth model considerably unpredicts the growth of the shear layer, as does the Baldwin-Lomax model. It is shown that the one-equation turbulence model needs improvement for the prediction of separated flow due to smooth adverse pressure gradient, and mixing layer growth.

Kandula, M.↗

Temperature-Corrected Model of Turbulence in Hot Jet Flows

An improved correction has been developed to increase the accuracy with which certain formulations of computational fluid dynamics predict mixing in shear layers of hot jet flows. The CFD formulations in question are those derived from the Reynolds-averaged Navier-Stokes equations closed by means of a two-equation model of turbulence, known as the k-epsilon model, wherein effects of turbulence are summarized by means of an eddy viscosity. The need for a correction arises because it is well known among specialists in CFD that two-equation turbulence models, which were developed and calibrated for room-temperature, low Mach-number, plane-mixing-layer flows, underpredict mixing in shear layers of hot jet flows. The present correction represents an attempt to account for increased mixing that takes place in jet flows characterized by high gradients of total temperature. This correction also incorporates a commonly accepted, previously developed correction for the effect of compressibility on mixing.

Abdol-Hamid, Khaled S.↗

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗

Solvent-Dependent Dynamics of Cellulose Nanocrystals in Process-Relevant Flow Fields

Flow-assisted alignment of anisotropic nanoparticles is a promising route for the bottom-up assembly of advanced materials with tunable properties. While aligning processes could be optimized by controlling factors such as solvent viscosity, flow deformation, and the structure of the particles themselves, it is necessary to understand the relationship between these factors and their effect on the final orientation. In this study, we investigated the flow of surface-charged cellulose nanocrystals (CNCs) with the shape of a rigid rod dispersed in water and propylene glycol (PG) in an isotropic tactoid state. In situ scanning small-angle X-ray scattering (SAXS) and rheo-optical flow-stop experiments were used to quantify the dynamics, orientation, and structure of the assigned system at the nanometer scale. The effects of both shear and extensional flow fields were revealed in a single experiment by using a flow-focusing channel geometry, which was used as a model flow for nanomaterial assembly. Due to the higher solvent viscosity, CNCs in PG showed much slower Brownian dynamics than CNCs in water and thus could be aligned at lower deformation rates. Moreover, CNCs in PG also formed a characteristic tactoid structure but with less ordering than CNCs in water owing to weaker electrostatic interactions. The results indicate that CNCs in water stay assembled in the mesoscale structure at moderate deformation rates but are broken up at higher flow rates, enhancing rotary diffusion and leading to lower overall alignment. Albeit being a study of cellulose nanoparticles, the fundamental interplay between imposed flow fields, Brownian motion, and electrostatic interactions likely apply to many other anisotropic colloidal systems.

36 MATERIALS SCIENCE↗

Tolerable ranges of fluid shear for early life-stage fishes: implications for safe fish passage at hydropower and irrigation infrastructure

Egg and larval fish drifting downstream are likely to encounter river infrastructure, increasing mortality. Elevated fluid shear is one likely cause of this. To confirm this and determine tolerable levels, egg and larvae of three Australian species were exposed to elevated shear in a laboratory flume. Mortality was modelled over a broad range of strain rates, allowing critical thresholds to be estimated. Eggs were very susceptible to mortality at low strain rates and 100% of golden and silver perch died once strain rate exceeded 629 and 148 cm s-1 cm-1, respectively. Larvae were less vulnerable than eggs, but mortality increased at higher strain rates and at younger ages. Murray cod larvae were the most susceptible, particularly around 9 days post hatch. Younger golden perch and silver perch were also susceptible but little mortality occurred once maturity reached ~25 days post hatch (nearing juvenile metamorphosis). The results indicate that fluid shear will contribute to the large mortalities of drifting eggs and larvae entrained at river infrastructure, with the impact reducing with age. Minimising turbulent flows and fluid shear should be the objective when designing and operating structures like spillways, turbines and irrigation weirs in areas of significant egg and larval drift.

Navarro, Anna↗

Toward large eddy simulation of turbulent flow over an airfoil

The flow field over an airfoil contains several distinct flow characteristics, e.g. laminar, transitional, turbulent boundary layer flow, flow separation, unstable free shear layers, and a wake. This diversity of flow regimes taxes the presently available Reynolds averaged turbulence models. Such models are generally tuned to predict a particular flow regime, and adjustments are necessary for the prediction of a different flow regime. Similar difficulties are likely to emerge when the large eddy simulation technique is applied with the widely used Smagorinsky model. This model has not been successful in correctly representing different turbulent flow fields with a single universal constant and has an incorrect near-wall behavior. Germano et al. (1991) and Ghosal, Lund & Moin have developed a new subgrid-scale model, the dynamic model, which is very promising in alleviating many of the persistent inadequacies of the Smagorinsky model: the model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model has been remarkably successful in prediction of several turbulent and transitional flows. We plan to simulate turbulent flow over a '2D' airfoil using the large eddy simulation technique. Our primary objective is to assess the performance of the newly developed dynamic subgrid-scale model for computation of complex flows about aircraft components and to compare the results with those obtained using the Reynolds average approach and experiments. The present computation represents the first application of large eddy simulation to a flow of aeronautical interest and a key demonstration of the capabilities of the large eddy simulation technique.

Choi, Haecheon↗

Interactive method for computation of viscous flow with recirculation

An interactive method is proposed for the solution of two-dimensional, laminar flow fields with identifiable regions of recirculation, such as the shear-layer-driven cavity flow. The method treats the flow field as composed of two regions, with an appropriate mathematical model adopted for each region. The shear layer is computed by the compressible boundary layer equations, and the slowly recirculating flow by the incompressible Navier-Stokes equations. The flow field is solved iteratively by matching the local solutions in the two regions. For this purpose a new matching method utilizing an overlap between the two computational regions is developed, and shown to be most satisfactory. Matching of the two velocity components, as well as the change in velocity with respect to depth is amply accomplished using the present approach, and the stagnation points corresponding to separation and reattachment of the dividing streamline are computed as part of the interactive solution. The interactive method is applied to the test problem of a shear layer driven cavity. The computational results are used to show the validity and applicability of the present approach.

Brandeis, J.↗

Turbulence measurements in hypersonic shock-wave boundary-layer interaction flows

Turbulent intensity and Reynolds shear stress measurements are presented for two nonadiabatic hypersonic shock-wave boundary-layer interaction flows, one with and one without separation. These measurements were obtained using a new hot-wire probe specially designed for heated flows. Comparison of the separated and attached flows shows a significant increase above equilibrium values in the turbulent intensity and shear stress downstream of the interaction region for the attached case, while for the separated case, the turbulent fluxes remain close to equilibrium values. This effect results in substantial differences in turbulence lifetimes for the two flows. It is proposed that these differences are due to a coupling between the turbulent energy and separation bubble unsteadiness, a hypothesis supported by the statistical properties of the turbulent fluctuations.

Mikulla, V.↗

Turbulence measurements in hypersonic shock-wave boundary-layer interaction flows

Turbulent intensity and Reynolds shear stress measurements are presented for two nonadiabatic hypersonic shock-wave boundary-layer interaction flows, one with and one without separation. These measurements were obtained using a new hot-wire probe specially designed for heated flows. Comparison of the separated and attached flows shows a significant increase above equilibrium values in the turbulent intensity and shear stress downstream of the interaction region for the attached case, while for the separated case, the turbulent fluxes remain close to equilibrium values. This effect results in substantial differences in turbulence lifetime for the two flows. We propose that these differences are due to a coupling between the turbulent energy and separation bubble unsteadiness, a hypothesis supported by the statistical properties of the turbulent fluctuations.

Mikulla, V.↗

Observations of intense velocity shear and associated electrostatic waves near an auroral arc

Simultaneous measurements of energetic particles and ac electric fields made by the javelin sounding rocket NASA 8:56 during the late expansion phase of a magnetic storm have revealed an intense shear in plasma flow of magnitude 20 (m/s)/m at the edge of an auroral arc. Structure with two characteristic scales sizes is displayed in the region of shear. Larger structures are of the order of several kilometers in size. Intense irregularities with characteristic wavelengths smaller than the scale size of the shear have also been detected. The large-scale changes in the orientation of the charge sheet at the edge of the arc may be due to the Kelvin-Helmholtz branch; shorter-wavelength modes may be related to the shear driven resistive drift wave. Observations are consistent with the suggestion that velocity shear instabilities may play a role in the formation of high-latitude irregularities.

Kelley, M. C.↗

Probing Nanorod Assembly and Dynamics in Polymer Nanocomposites in Equilibrium and Shear

Coarse-grained molecular dynamics simulations are used to examine the structure and dynamics of nanorod assemblies in polymer melts under equilibrium and simple shear. We show that as the concentration of nanorods increases, there is a transition from an isotropic phase to a two-phase region in which the nanorods phase separate into a dilute phase and dense bundles of hexagonally packed nanorods. The onset of the two-phase region is below that predicted by Onsager theory, which we attribute to an effective increase in the diameter of nanorods due to a layer of polymer bound to the rod surfaces. Equilibrium simulations show that increasing polymer chain length N enhances nanorod bundling at fixed nanorod concentration. Simulations of systems undergoing simple shear show that flow enhances nanorod alignment and bundling relative to those of equilibrium systems with similar properties. Finally, simulations reveal that increasing nanorod length enhances nanorod alignment under shear at equivalent shear rates but not at equivalent Péclet numbers. Overall, our simulations highlight that polymer matrix-nanorod attraction (i.e., chemistry), shear rate, and matrix chain length are desirable design variables to control the structure of nanorod-containing soft materials under simple shear.

diffusion↗