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,081 records · Page 60

On solving the compressible Navier-Stokes equations for unsteady flows at very low Mach numbers

The properties of a preconditioned, coupled, strongly implicit finite difference scheme for solving the compressible Navier-Stokes equations in primitive variables are investigated for two unsteady flows at low speeds, namely the impulsively started driven cavity and the startup of pipe flow. For the shear-driven cavity flow, the computational effort was observed to be nearly independent of Mach number, especially at the low end of the range considered. This Mach number independence was also observed for steady pipe flow calculations; however, rather different conclusions were drawn for the unsteady calculations. In the pressure-driven pipe startup problem, the compressibility of the fluid began to significantly influence the physics of the flow development at quite low Mach numbers. The present scheme was observed to produce the expected characteristics of completely incompressible flow when the Mach number was set at very low values. Good agreement with incompressible results available in the literature was observed.

Pletcher, R. H.↗

Drop Breakup in Fixed Bed Flows as Model Stochastic Flow Fields

We examine drop breakup in a class of stochastic flow fields as a model for the flow through fixed fiber beds and to elucidate the general mechanisms whereby drops breakup in disordered, Lagrangian unsteady flows. Our study consists of two parallel streams of investigation. First, large scale numerical simulations of drop breakup in a class of anisotropic Gaussian fields will be presented. These fields are generated spectrally and have been shown in a previous publication to be exact representations of the flow in a dilute disordered bed of fibers if close interactions between the fibers and the drops are dynamically unimportant. In these simulations the drop shape is represented by second and third order small deformation theories which have been shown to be excellent for the prediction of drop breakup in steady strong flows. We show via these simulations that the mechanisms of drop breakup in these flows are quite different than in steady flows. The predominant mechanism of breakup appears to be very short lived twist breakups. Moreover, the occurrence of breakup events is poorly predicted by either the strength of the local flow in which the drop finds itself at breakup, or the degree of deformation that the drop achieves prior to breakup. It is suggested that a correlation function of both is necessary to be predictive of breakup events. In the second part of our research experiments are presented where the drop deformation and breakup in PDMS/polyisobutylene emulsions is considered. We consider very dilute emulsions such that coalescence is unimportant. The flows considered are simple shear and the flow through fixed fiber beds. Turbidity, small angle light scattering, dichroism and microscopy are used to interrogate the drop deformation process in both flows. It is demonstrated that breakup at very low capillary numbers occurs in both flows but larger drop deformation occurs in the fixed bed flow. Moreover, it is witnessed that breakup in the bed occurs continuously during flow and apparently with uniform probability through the bed length. The drop deformations witnessed in our experiments are larger than those predicted by the numerical simulations, and future plans to investigate these differences are discussed.

Shaqfeh, Eric S. G.↗

Atomic Force Microscope Mediated Chromatography

The atomic force microscope (AFM) is used to inject a sample, provide shear-driven liquid flow over a functionalized substrate, and detect separated components. This is demonstrated using lipophilic dyes and normal phase chromatography. A significant reduction in both size and separation time scales is achieved with a 25-micron-length column scale, and one-second separation times. The approach has general applications to trace chemical and microfluidic analysis. The AFM is now a common tool for ultra-microscopy and nanotechnology. It has also been demonstrated to provide a number of microfluidic functions necessary for miniaturized chromatography. These include injection of sub-femtoliter samples, fluidic switching, and sheardriven pumping. The AFM probe tip can be used to selectively remove surface layers for subsequent microchemical analysis using infrared and tip-enhanced Raman spectroscopy. With its ability to image individual atoms, the AFM is a remarkably sensitive detector that can be used to detect separated components. These diverse functional components of microfluidic manipulation have been combined in this work to demonstrate AFM mediated chromatography. AFM mediated chromatography uses channel-less, shear-driven pumping. This is demonstrated with a thin, aluminum oxide substrate and a non-polar solvent system to separate a mixture of lipophilic dyes. In conventional chromatographic terms, this is analogous to thin-layer chromatography using normal phase alumina substrate with sheardriven pumping provided by the AFM tip-cantilever mechanism. The AFM detection of separated components is accomplished by exploiting the variation in the localized friction of the separated components. The AFM tip-cantilever provides the mechanism for producing shear-induced flows and rapid pumping. Shear-driven chromatography (SDC) is a relatively new concept that overcomes the speed and miniaturization limitations of conventional liquid chromatography. SDC is based on a sliding plate system, consisting of two flat surfaces, one of which has a recessed channel. A fluid flow is produced by axially sliding one plate past another, where the fluid has mechanical shear forces imposed at each point along the channel length. The shear-induced flow rates are very reproducible, and do not have pressure or voltage gradient limitations. SDC opens up a new range of enhanced separation kinetics by permitting the sample confinement with submicron dimensions. Small, highly confined liquid is advantageous for chromatographic separation because the separation rate is known to scale according to the square of the confined sample diameter. In addition, because shear-driven flows are not limited by fluid velocity, shear-driven liquid chromatography may provide up to 100,000 plate efficiency.

Anderson, Mark S.↗

A Quasi-3-D Theory for Impedance Eduction in Uniform Grazing Flow

A 2-D impedance eduction methodology is extended to quasi-3-D sound fields in uniform or shearing mean flow. We introduce a nonlocal, nonreflecting boundary condition to terminate the duct and then educe the impedance by minimizing an objective function. The introduction of a parallel, sparse, equation solver significantly reduces the wall clock time for educing the impedance when compared to that of the sequential band solver used in the 2-D methodology. The accuracy, efficiency, and robustness of the methodology is demonstrated using two examples. In the first example, we show that the method reproduces the known impedance of a ceramic tubular test liner. In the second example, we illustrate that the approach educes the impedance of a four-segment liner where the first, second, and fourth segments consist of a perforated face sheet bonded to honeycomb, and the third segment is a cut from the ceramic tubular test liner. The ability of the method to educe the impedances of multisegmented liners has the potential to significantly reduce the amount of time and cost required to determine the impedance of several uniform liners by allowing them to be placed in series in the test section and to educe the impedance of each segment using a single numerical experiment. Finally, we probe the objective function in great detail and show that it contains a single minimum. Thus, our objective function is ideal for use with local, inexpensive, gradient-based optimizers.

Watson, W. R.↗

Viscometric flow of dense granular materials under controlled pressure and shear stress

This study examines the flow of dense granular materials under external shear stress and pressure using discrete element method simulations. In this method, the material is allowed to strain along all periodic directions and adapt its solid volume fraction in response to an imbalance between the internal state of stress and the external applied stress. By systematically varying the external shear stress and pressure, the steady rheological response is simulated for: (1) rate-independent quasi-static flow; and (2) rate-dependent inertial flow. The simulated flow is viscometric with non-negligible first and second normal stress differences. While both normal stress differences are negative in inertial flows, the first normal stress difference switches from negative to slightly positive, and second normal stress difference tends to zero in quasi-static flows. The first normal stress difference emerges from a lack of coaxiality between a second-rank contact fabric tensor and strain rate tensor in the flow plane, while the second normal stress difference is linked to an excess of contacts in the shear plane compared with the vorticity direction. A general rheological model of second order (in terms of strain rate tensor) is proposed to describe the two types of flow, and the model is calibrated for various values of interparticle friction from simulations on nearly monodisperse spheres. The model incorporates normal stress differences in both regimes of flow and provides a complete viscometric description of steady dense granular flows.

42 ENGINEERING↗

Early time evolution of negative ion clouds and electron density depletions produced during electron attachment chemical release experiments

Two-dimensional electrostatic particle-in-cell simulations are used to study the early time evolution of electron depletions and negative ion clouds produced during electron attachment chemical releases in the ionosphere. The simulation model considers the evolution in the plane perpendicular to the magnetic field and a three-species plasma that contains electrons, positive ions, and also heavy negative ions that result as a by-product of the electron attachment reaction. The early time evolution (less than the negative ion cyclotron period) of the system shows that a negative charge surplus initially develops outside of the depletion boundary as the heavy negative ions move across the boundary. The electrons are initially restricted from moving into the depletion due to the magnetic field. An inhomogenous electric field develops across the boundary layer due to this charge separation. A highly sheared electron flow velocity develops in the depletion boundary due to E x B and Delta-N x B drifts that result from electron density gradients and this inhomogenous electric field. Structure eventually develops in the depletion boundary layer due to low-frequency electrostatic waves that have growth times shorter than the negative ion cyclotron period. It is proposed that these waves are most likely produced by the electron-ion hybrid instability that results from sufficiently large shears in the electron flow velocity.

Scales, W. A.↗

Scattering of coherent sound waves by atmospheric turbulence

An analytical study of the propagation of coherent sound waves through an atmosphere containing both mean and fluctuating flow variables is presented. The general flow problem is formulated as a time-dependent wave propagation in a half-space containing the turbulent medium. The coherent acoustic waves are analyzed by a smoothing technique, assuming that mean flow variables vary with the height only. The general equations for the coherent waves are derived, and then applied to two special cases, corresponding to uniform and shear mean flow, respectively. The results show that mean shear and turbulence introduce pronounced effects on the propagation of coherent acoustic disturbances.

Chow, P. L.↗

High-frequency sound propagation in a spatially varying mean flow

An equation for acoustic ray paths in a spatially varying mean flow was examined to determine some of the characteristics of the flow gradient effects on sound propagation. In a potential flow, the acoustic rays are deflected in the direction of increasing mean flow, and the gradient of the mean flow speed is the dominant factor causing the ray deflection. In contrast, in a sheared mean flow, the vorticity is the dominant factor in deflection of the acoustic rays.

Cho, Y. C.↗

High-frequency sound propagation in a spatially varying mean flow

An equation for acoustic ray paths in a spatially varying mean flow has been examined to determine some of the characteristics of the flow gradient effects on sound propagation. In a potential flow the acoustic rays are deflected in the direction of increasing mean flow, and the gradient of the mean flow speed is the dominant factor causing the ray deflection. In contrast, in a sheared mean flow, the vorticity is the dominant factor in deflection of the acoustic rays.

Cho, Y. C.↗

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗

Fiber orientation and porosity in large-format extrusion process: The role of processing parameters

Controlling fiber orientation and porosity in short-fiber thermoplastic composites is important for enhancing mechanical, electrical and thermal properties in large-format additive manufacturing. This study employs a factorial design of experiments (DoE) to assess the effects of nozzle diameter (5.08 mm–10.16 mm), temperature (230–250 °C), and extruder screw speed (150–280 rpm) on flow rate, shear rate, porosity, fiber orientation, fiber length and tensile strength in 20 % carbon fiber-filled acrylonitrile butadiene styrene. ANOVA results show that screw speed significantly impacts flow rate, while nozzle diameter and temperature have lesser effects. Shear rate increases with smaller nozzles and higher speeds. Porosity decreases from 5.58 % with a 10.16 mm nozzle to 3.11 % with a 5.08 mm nozzle at 150 rpm due to increased shear rates, which induce shear thinning, reducing viscosity and facilitating gas escape. Larger nozzles (10.16 mm) produce larger, more heterogeneous pores, while smaller nozzles (5.08 mm) yield smaller, uniform pores. Beads produced with the 5.08 mm nozzle exhibit longer fiber lengths due to reduced residence time, lower shear stress, and better alignment. Fiber orientation improves with smaller nozzles due to higher shear rates but decreases with higher screw speeds (280 rpm) due to shorter residence times. The highest fiber alignment (A xx ∼ 0.65) and low porosity (∼3%) were achieved with a 5.08 mm nozzle at 150 rpm, while equivalent additive manufacturing-compression molding samples exhibited better tensile strength (∼93 MPa) under these conditions. In conclusion, these findings emphasize the importance of optimizing processing parameters to enhance fiber alignment and reduce porosity for improved mechanical performance.

Design of experiments↗

Early time evolution of a chemically produced electron depletion

The early time evolution of an ionospheric electron depletion produced by a radially expanding electron attachment chemical release is studied with a two-dimensional simulation model. The model includes electron attachment chemistry, incorporates fluid electrons, particle ions and neutrals, and considers the evolution in a plane perpendicular to the geomagnetic field for a low beta plasma. Timescales considered are of the order of or less than the cyclotron period of the negative ions that result as a by-product of the electron attacment reaction. This corresponds to time periods of tenths of seconds during recent experiemts. Simulation results show that a highly sheared azimuthal electron flow velocity develops in the radially expanding depletion boundary. This sheared electron flow velocity and the steep density gradients in the boundary give rise to small-scale irregulatities in the form of electron density cavities and spikes. The nonlinear evolution of these irregularities results in trapping and ultimately turbulent heating of the negative ions.

Scales, W. A.↗

Transmission of sound through nonuniform circular ducts with compressible mean flows

An acoustic theory is developed to determine the sound transmission and attenuation through an infinite, hard-walled or lined, circular duct carrying compressible, sheared, mean flows and having a variable cross section. The theory is applicable to large as well as small axial variations, as long as the mean flow does not separate. Although the theory is described for circular ducts, it is applicable to other duct configurations - annular, two dimensional, and rectangular. The theory is described for the linear problem, but the technique is general and has the advantage of being applicable to the nonlinear case as well as the linear case. The technique is based on solving for the envelopes of the quasi-parallel acoustic modes that exist in the duct instead of solving for the actual wave. A computer program was developed. The mean flow model consists of a one dimensional flow in the core and a quarter-sine profile in the boundary layer. Results are presented for the reflection and transmission coefficients in ducts with varying slopes and carrying different mean flows.

Nayfeh, A. H.↗

A wave-envelope of sound propagation in nonuniform circular ducts with compressible mean flows

An acoustic theory is developed to determine the sound transmission and attenuation through an infinite, hard-walled or lined circular duct carrying compressible, sheared, mean flows and having a variable cross section. The theory is applicable to large as well as small axial variations, as long as the mean flow does not separate. The technique is based on solving for the envelopes of the quasi-parallel acoustic modes that exist in the duct instead of solving for the actual wave, thereby reducing the computation time and the round-off error encountered in purely numerical techniques. The solution recovers the solution based on the method of multiple scales for slowly varying duct geometry. A computer program was developed based on the wave-envelope analysis for general mean flows. Results are presented for the reflection and transmission coefficients as well as the acoustic pressure distributions for a number of conditions: both straight and variable area ducts with and without liners and mean flows from very low to high subsonic speeds are considered.

Nayfeh, A. H.↗

Transmission of sound through nonuniform circular ducts with compressible mean flows

An acoustic theory is developed to determine the sound transmission and attenuation through an infinite hard-walled or lined circular duct carrying compressible, sheared mean flows and having a variable cross section. The theory is applicable to large as well as small axial variations, as long as the mean flow does not separate. The technique is based on solving for the envelopes of the quasi-parallel acoustic modes that exist in the duct instead of solving for the actual wave, thereby reducing the computation time and the round-off error encountered in purely numerical techniques. A number of test cases that demonstrate the flexibility of the program are included. Convergence of the transmission coefficients and the acoustic pressure profiles with an increasing number of modes is illustrated.

Nayfeh, A. H.↗

Dynamical evolution of a solar coronal magnetic field arcade

Calculations of the long-term dynamical evolution of a solar coronal magnetic field arcade which is subjected to shearing photospheric flows are presented. The evolution is obtained by numerical solution of a subset of the resistive magnetohydrodynamic equations. For a simplified model of the bipolar magnetic field observed in the solar corona, it is found that photospheric flow produces a slow evolution of the magnetic field, with a buildup of magnetic energy. For certain photospheric shear profiles, the field configuration produced is linearly unstable to an ideal magnetohydrodynamic mode when the shear exceeds a critical value. The nonlinear evolution of this instability shows the spontaneous formation of current sheets. Reconnection of the magnetic field produces a rapid release of magnetic energy. The major fraction of the energy is dissipated resistively, while a small fraction is converted into kinetic energy of an ejected plasmoid. The relevance of these results to two-ribbon flares is discussed.

Mikic, Z.↗

Second-order closure analysis of turbulent reacting flows

A reacting shear layer program for a complete second-order closure of turbulent reacting flows is developed. The program uses a typical eddy box model and many of the fluid mechanical turbulence models developed in earlier studies. The computer program actually solves a total of 30 equations, where the numerical integration of the equations is performed by a forward-time-centered space quasi-implicit upwind finite-difference scheme. Results for nonreacting and simple isothermal reacting flows are presented along with preliminary results for propane-air flames. It is shown that the program predictions for simple nonreacting flows are in good agreement with available data. The calculations for isothermal reacting flows demonstrate the importance of including a mixedness correlation in finite rate chemistry calculations, where significant errors in the flowfield properties can result from neglect of the interaction between turbulence and the chemistry. Results for propane-air diffusion flames further support this conclusion.

Varma, A. K.↗

Screening of resonant magnetic perturbation fields assuming various plasma flow models

A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗