SEARCH · Search NASA
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.
The Chebyshev matrix collocation method is applied to obtain the spatial modes of the Orr-Sommerfeld equation for Poiseuille flow and the Blausius boundary layer. The problem is linearized by the companion matrix technique for semi-infinite domain using a mapping transformation. The method can be easily adapted to problems with different boundary conditions requiring different transformations.
Spectral methods have been used to great advantage in hydrodynamic stability calculations; the concepts are described in Orszag's seminal application of the Chebyshev tau method to the Orr-Sommerfeld equation for plane Poiseuille flow in 1971. Orszag discusses both the Chebyshev Galerkin and the Chebyshev tau methods, but presents results for the tau method, which is easier to implement than the Galerkin method. The tau method has the disadvantage that two unstable eigenvalues are produced that are artifacts of the discretization. An extremely simple modification to the Chebyshev tau method is presented which eliminates the spurious eigenvalues. First a simplified model of the Orr-Sommerfeld equation discussed by Gottlieb and Orszag was studied. Then the Chebyshev tau method is considered, which has two spurious eigenvalues, and then a modification which eliminates them is described. Finally, results for the Orr-Sommerfeld equation are considered where the modified tau method also eliminates the spurious eigenvalues. The simplicity of the modification makes it a convenient alternative to other approaches to the problem.
The stability characteristics of a laminar, axisymmetric jet, issuing from fully developed Poiseuille flow, are investigated. The jet preferred frequency, as inferred from surveys of u'-spectra, is found to yield a Strouhal number (St) that depends on the Reynolds number (R); St and R are based on the jet diameter (D) and the average velocity (U sub av) at the jet origin. The value of St increases with increasing R in the range 400 less than approximately R less than approximately 4000, attaining an asymptotic value of about 0.45. Flow visualization confirms that the instability is primarily in a helical mode, as predicted by stability analyses. Analyses do predict a similar St versus R variation in approximately the correct St-range. However, the R-range where this is predicted is lower than that found experimentally.
The Chebyshev matrix collocation method is applied to obtain the spatial modes of the Orr-Sommerfeld equation for Poiseuille flow and the Blasius boundary layer. The problems is linearized by the companion matrix technique for semiinfinite domain using a mapping transformation. The method can be easily adapted to problems with different boundary conditions requiring different transformations.
A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.
The applicability of active control by periodic suction blowing in spatially evolving plane Poiseuille flow is investigated by the direct simulation of the three-dimensional, incompressible Navier-Stokes equations. The results reveal that significant reductions in perturbation amplitudes can be obtained by a proper choice of the control wave amplitude and phase. The upstream influence of the control wave is shown to be confined to a region in the vicinity of the control slot with no apparent effect on the flow development.
The fourth-order finite-difference scheme with fully implicit time-marching presently used to computationally study the spatial instability of planar Poiseuille flow incorporates a novel treatment for outflow boundary conditions that renders the buffer area as short as one wavelength. A semicoarsening multigrid method accelerates convergence for the implicit scheme at each time step; a line-distributive relaxation is developed as a robust fast solver that is efficient for anisotropic grids. Computational cost is no greater than that of explicit schemes, and excellent agreement with linear theory is obtained.
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.
We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitras~y inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.
Velocity profiles of non-Newtonian fluid /polyisobutylene solution/ in helical flow by flow visualization technique, discussing Poiseville and viscometric flows
The recent hydrodynamic phonon transport theory for graphitic materials has been supported by the measurements of the second sound at temperatures up to about 100 K. When boundary scattering becomes comparable to momentum-conserving normal phonon scattering processes that are responsible for phonon hydrodynamics, Poiseuille phonon flow phenomena can emerge to give rise to unique size-dependent thermal conductivity in thin graphite. Here, we examine the thickness range for the Poiseuille phonon flow to become observable in thin graphite with the use of both deviational Monte Carlo simulation of the Peierls-Boltzmann transport equation and four-probe thermal transport measurements. As the basal-plane thermal conductivity calculated by prior first-principles theories saturates to the graphite value when the thickness is increased to five graphene layers, the phonon dispersion of graphite is used in the current calculations of thin graphite of micrometer thickness and a 23-layer thick ultrathin graphite (UTG) sample. The calculations show that diffuse surface scattering by surface defects can lead to Poiseuille phonon flow at 50 K in thin graphite with the thickness close to several micrometers but not in the 65 μm thin graphite and 23-layer UTG, where phonon scattering by the top and bottom surfaces become, respectively, much less and more frequent than the normal processes. In addition, the calculation results with the bulk graphite dispersion and diffuse surface scattering show decreased basal-plane thermal conductivity with decreasing thickness, opposite to recent thermocouple measurements of thin graphite samples. In comparison, the calculation results reveal that partially diffuse surface defect scattering can yield the four-probe measurement results of UTG samples, which are prepared here with an improved process to minimize surface contamination.
Direct numerical simulations of turbulent spots in plane Poiseuille and boundary-layer flows are performed. Mature, self-similar spots are obtained. The propagation velocities and spreading angles are found to compare well with corresponding experiments. The difference in shape of the two spots is also clearly discernible: the turbulent parts are contained within arrowhead regions that point in opposite directions for the two cases. The wing-tip region of the Poiseuille spot is also found to consist of a large-amplitude semiturbulent wave packet.
Two mechanisms of drag reduction for flow over flat plates were investigated. The first mechanism employs Bushnell's hypothesis that compliant walls produce drag reduction by interfering with the formation of the turbulent spots in a turbulent boundary layer. It is shown that the amplitudes and frequencies of compliant wall motions for drag reduction might be achieved by using slightly curved walls and the resulting large amplitude motions of snap buckling. A simple structural model of an arch is used in the analysis, and an asymptotic method is developed. The required wall motions can be obtained by using materials like mylar. In addition, the delay of transition from laminar to turbulent flow by driven walls was studied for Poiseuille channel flow. The walls are driven by a periodic traveling wave. A significant increase in the transitional Reynolds number is obtained by appropriately prescibing the wavelength and phase velocity of the wall motion. Previously developed asymptotic methods are used in the analysis.
Two mechanisms of drag reduction for flow over flat plates were investigated. The first mechanism employs Bushnell's hypothesis that compliant walls produce drag reduction by interfering with the formation of the turbulent spots in a turbulent boundary layer. It is shown that the amplitudes and frequencies of compliant wall motions for drag reduction might be achieved by using slightly curved walls and the resulting large amplitude motions of snap buckling. A simple structural model of an arch is used in the analysis, and an asymptotic method is developed. The required wall motions can be obtained by using materials like mylar. In addition, the delay of transition from laminar to turbulent flow by driven walls was studied for Poiseuille channel flow. The walls are driven by a periodic traveling wave. A significant increase in the transitional Reynolds number is obtained by appropriately prescribing the wavelength and phase velocity of the wall motion. Previously developed asymptotic methods are used in the analysis. Previously announced in STAR as N83-11061
A viewgraph presentation is made showing the capabilities of the computer code SPIRALI. Overall capabilities of SPIRALI include: computes rotor dynamic coefficients, flow, and power loss for cylindrical and face seals; treats turbulent, laminar, Couette, and Poiseuille dominated flows; fluid inertia effects are included; rotor dynamic coefficients in three (face) or four (cylindrical) degrees of freedom; includes effects of spiral grooves; user definable transverse film geometry including circular steps and grooves; independent user definable friction factor models for rotor and stator; and user definable loss coefficients for sudden expansions and contractions.
This invention relates generally to electrophoresis devices and more particularly to a moving wall, continuous flow device in which an electrophoresis chamber is angularly positionable with respect to the direction of moving belt walls. A frame with an electrophoresis chamber is rotatably supported between two synchronously driven belt walls. This allows the chamber to be angularly positionable with respect to the direction of belt travel, which compensates for electroosmotic flow within the electrophoresis chamber. Injection of a buffer solution via an opening and a homogenous sample stream via another opening is performed at the end of a chamber, and collection of buffer and the fractionated species particles is done by a conventional collection array at an opposite end of the chamber. Belts are driven at a rate which exactly matches the flow of buffer and sample through the chamber, which entrains the buffer to behave as a rigid electrophoretic medium, eliminating flow distortions (Poiseuille effect). Additionally, belt material for each belt is stored at one end of the device and is taken up by drive wheels at an opposite end. The novelty of this invention particularly lies in the electrophoresis chamber being angularly positionable between two moving belt walls in order to compensate for electroosmotic flow. Additionally, new belt material is continuously exposed within the chamber, minimizing flow distortion due to contamination of the belt material by the sample.
The refraction of acoustic waves by a moving medium layer is theoretically treated and the expressions for reflection and transmission coefficients are determined. The moving medium layer velocity is assumed to have a space dependence in one direction. A partitioning of the moving medium layer into constant-velocity sublayers is introduced and the number of sublayers is allowed to increase until the reflection and transmission coefficients converage to their respective values. Numerical results for several sublayer approximations of Poiseuille's flow are presented as functions of the moving layer velocity for several angles of incidence of the acoustic wave. The degenerate case of single constant-velocity layer is also treated, both theoretically and by a numerical analysis.
Research accomplishments are summarized and publications generated under the contract are listed. The general purpose of the research was to investigate various symmetry breaking problems in fluid mechanics by the use of structure parameters and selection principles. Although all of the nonlinear problems studied involved systems of partial differential equations, many of these problems led to the study of a single nonlinear operator equation of the form F(w, lambda, gamma) = 0, (w is an element of H), (lambda is an element of R1), (gamma is an element of R1). Instead of varying only the load parameter lambda, as is often done in the study of such equations, one of the main ideas used was to vary the structure parameter gamma in such a way that stable solutions were obtained. In this way one determines detailed stability results by making use of the structure of the model equations and the known physical parameters of the problem. The approach was carried out successfully for Benard-type convection problems, Taylor-like problems for short cylinders, rotating Couette-Poiseuille channel flows, and plane Couette flows. The main focus of the research was on wave theory of vortex breakdown in a tube. A number of preliminary results for inviscid axisymmetric flows were obtained.