Chemical reaction in a turbulent flow field with uniform velocity gradient
Statistical transport, energy and chemical reaction rate in turbulent flow field with uniform velocity gradient, using Brownian motion theory
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.
Statistical transport, energy and chemical reaction rate in turbulent flow field with uniform velocity gradient, using Brownian motion theory
The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This modele is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence.
The turbulent, incompressible reattaching flow over a rearward-facing step has been studied by many researchers over the years. One of the principal quantities determined in these experiments has been the distance from the step to the point (or region) where the separated shear layer reattaches to the surface (x(r)). The values for x(r)/h, where h is the step height, have covered a wider range than can reasonably be attributed to experimental technique or inaccuracy. Often the reason for a largely different value of x(r)/h can be attributed to an incompletely developed turbulent layer, or a transitional or laminar boundary layer. However, for the majority of experiments where the boundary layer is believed to be fully developed and turbulent, x(r)/h still varies several step heights; generally, 5 1/2 approximately < x(r)/h approximately < 7 1/2. This observed variation has usually been attributed to such variables as l/h (step length to height, h/delta (step height to initial boundary-layer thickness), R(e)(theta)), or the experimental technique for determining reattachment location. However, there are so many different combinations of variables in the previous experiments that it was not possible to sort out the effects of particular conditions on the location of reattachment. In the present experiment velocity profiles have been measured in and around the region of separated flow. Results show a large influence of adverse pressure gradient on the reattaching flow over a rearward-facing step that has not been reported previously. Further, the many previous experiments for fully developed, turbulent flow in parallel-walled channels have shown a range of reattachment location that has not been explained by differences in initial flow conditions. Although these initial flow conditions might contribute to the observed variation of reattachment location, it appears that the pressure gradient effect can explain most of that variation.
Report describes experimental study of turbulence in two boundary-layer flows with adverse gradients of pressure. Flows produced about cylinder oriented with axis along that of low-speed wind tunnel of rectangular cross section. Fluctuations of velocities analyzed with respect to various mathematical models of turbulence.
Measurements have been made of the mean flow properties and turbulent fluctuations in adiabatic turbulent boundary layer flows subjected to distributed adverse pressure gradients. In the freestream region upstream of the adverse pressure gradient the Mach number was 3.86, the unit Reynolds number 5.3 million per foot. The boundary layer developed on the wall of an axisymmetric nozzle and straight test section. In order to avoid the effects of streamwise surface curvature the adverse pressure gradients at the test section wall were induced by contoured centerbodies mounted on the wind tunnel centerline. The flow under study simulated that which might be found in an axially symmetric engine inlet of a supersonic aircraft.
Theoretical and experimental research on fluid conductivity of lunar surface materials is summarized. Theoretical methods were developed for the analysis of transitional and free-molecular flows, and for analysis of lunar permeability probe data in general. Experimental studies of rarefied flows under conditions of a large pressure gradient show flows in the continuum regime to be responsible for the largest portion of the pressure drop between source and sink for one dimensional flow, provided the entrance Knudsen number is sufficiently small. The concept of local similarity leading to a universal nondimensional function of Knudsen number was shown to have approximate validity; flows in all regimes may be described in terms of an area fraction and a single length parameter. Synthetic porous media prepared from glass beads exhibited flow behavior similar in many regards to that of a natural sandstone; studies using artificial stones with known pore configurations may lead to new insight concerning the structure of natural materials. The experimental method involving the use of segmented specimens of large permeability is shown to be fruitful.
The space shuttle main engine (SSME) has extremely complex internal flow structure. The geometry of the flow domain is three-dimensional with complicated topology. The flow is compressible, viscous, and turbulent with large gradients in flow quantities and regions of recirculations. The analysis of the flow field in SSME involves several tedious steps. One is the geometrical modeling of the particular zone of the SSME being studied. Accessing the geometry definition, digitalizing it, and developing surface interpolations suitable for an interior grid generator require considerable amount of manual labor. There are several types of grid generators available with some general-purpose finite element programs. An efficient and robust computational scheme for solving 3D Navier-Stokes equations has to be implemented. Post processing software has to be adapted to visualize and analyze the computed 3D flow field. The progress made in a project to develop software for the analysis of the flow is discussed. The technical approach to the development of the finite element scheme and the relaxation procedure are discussed. The three dimensional finite element code for the compressible Navier-Stokes equations is listed.
A numerical study of the time variation of the angular momentum and mass capture rates for a central object accreting from a uniform medium with a velocity gradient transverse to the direction of the mean flow is presented, covering a range of velocity asymmetries and Mach numbers in the incident flow. It is found that the mass accretion rate in a given evolutionary sequence varies in an irregular manner, with the matter accreting onto the central object from either a continuously moving accretion wake or from an accretion disk. The implications of the results from the study of short-term fluctuations observed in the pulse period and luminosity of X-ray pulsars are discussed.
Transport aircraft performance is strongly influenced by the effectiveness of high-lift systems. Developing wakes generated by the airfoil elements are subjected to strong pressure gradients and can thicken very rapidly, limiting maximum lift. This paper focuses on the effects of various pressure gradients on developing symmetric wakes and on the ability of a linear eddy viscosity model and a non-linear explicit algebraic stress model to accurately predict their downstream evolution. In order to reduce the uncertainties arising from numerical issues when assessing the performance of turbulence models, three different numerical codes with the same turbulence models are used. Results are compared to available experimental data to assess the accuracy of the computational results.
Explore the source record for details and available documents.
Experimental results are presented for separating flows over a NASA GA(W)-1 airfoil with 2-percent trailing-edge thickness, with emphasis on the mean velocity data that demonstrate flowfield similarity downstream of separation. The length scale presently employed was selected because it furnishes the requisite matching point between the outer, mixing region flow and the back-flow.
Standard approaches to neoclassical theory do not extend into regions of strong gradients in tokamaks such as the pedestal and internal transport barriers. Here, we calculate the modifications to neoclassical electron physics inside strong gradient regions of large aspect ratio tokamaks in the banana regime. We show that these modifications are due to the different ion flow and the strong poloidal variation of the potential. We also provide a physical interpretation of the mechanisms that drive poloidal asymmetries and hence a poloidal electric field. We apply our model to two specific example cases of pedestal profiles, calculating the neoclassical electron flux and the bootstrap current. We find that, depending on the ion flow, weak gradient neoclassical theory overestimates or underestimates the neoclassical electron transport and the bootstrap current in regions with strong gradients. We show that the determination of the mean parallel flow is more complex than in weak gradient neoclassical theory. For vanishing turbulence, we can determine the radial electric field for a given flow profile in the pedestal.
There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests that for many proteins crystal nucleation and growth is by addition of associated species that are preformed by reversible concentration-driven self association processes in the bulk solution. We have developed a self-association model for the crystal nucleation and growth of the protein chicken egg lysozyme. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small an octamer in the 43 helix configuration (the proposed average sized growth unit) would have a M.W. approx. 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm (a 24-mer) commonly attach to the crystal. AFM results from Weichmann et al. (Ultramicroscopy 86, 159-166, 2001) suggest that associated species of up to 40-mers in size add to the (101) faces. These measurements reflect the sizes of units that both added and desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. On Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species is more likely to be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.
An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations are integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and the convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the lower-upper (L-U)-successive symmetric over-relaxation iterative scheme is more efficient than a preconditioner based on incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional line Gauss-Seidel relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.
Turbulence measurements were obtained for two cases of boundary layer flow with an adverse pressure gradient, one attached and the other separated. A three-component laser Doppler velocimeter system was used to measure three mean velocity components, all six Reynolds stress components, and all ten velocity triple product correlations. Independent measurements of skin-friction obtained with a laser oil-flow interferometer were used to examine the law of the wall in adverse pressure gradient flows where p(+) is less than 0.05. Strong similiarities were seen between the two adverse pressure gradient flows and free shear layer type flows. Eddy viscosities, dissipation rates, and pressure-strain rates were deduced from the data and compared to various turbulence modeling assumptions.
The effect of basic flow density gradients on the supersonic wall modes were investigated in Part 1 of this analysis. In that investigation only the 2-D modes were studied. Tam and Hu investigated the 3-D modes in a confined vortex sheet and reported that the first 2-D Class A mode (A01) had the highest growth rate compared to all other 2-D and 3-D modes present in the vortex sheet for that particular set of flow patterns. They also showed that this result also held true for finite thickness shear layers with delta(sub w) less than 0.125. For free shear layers, Sandham and Reynolds showed that the 3-D K-H mode became the dominant mode for M(sub c) greater than 0.6. Jackson and Grosch investigated the effect of crossflow and obliqueness on the slow and fast odes present in a M(sub c) greater than 1 environment and showed that for certain combination of crossflow and wave angles the growth rates could be increased by up to a factor of 2 with respect to the 2-D case. The case studied here is a confined shear layer shown in Part 1. All solution procedures and basic low profiles are the same as in Part 1. The effect of density gradients on the 3-D modes present in the density ratios considered in Part 1 are investigated.
A series of preliminary experiments were done to investigate solutal flow about growing lysozyme crystals and its effects. Density-gradient-driven flow was observed using a schlieren optical system. Crystals used ranged from 0.3 to 1.72 mm across the (110) face, and protein concentrations were from 3.7 to 23.7 mg/ml. The convective plume velocities were found to be from 10 to 50 microns/s, which correlated with those predicted to occur based upon a diffusive-convective model. When microcrystals of lysozyme, less than 20 microns across the (110) face were subjected to directed solution flows, the growth rate was found to rapidly decrease over the 8-20 h course of the experiment. Solution flow rates used ranged from 18 to 40 microns/s, and protein concentrations were from 7.3 to 11.7 mg/ml, conditions typical of larger (greater than 0.5 mm) crystals in the terminal phases of a typical crystal growth procedure.
We consider a viscous incompressible fluid flow driven between two parallel plates by a constant pressure gradient. The flow is at a finite Reynolds number, with an 0(l) disturbance in the form of a traveling wave. A phase equation approach is used to discuss the evolution of slowly varying fully nonlinear two dimensional wavetrains. We consider uniform wavetrains in detail, showing that the development of a wavenumber perturbation is governed by Burgers equation in most cases. The wavenumber perturbation theory, constructed using the phase equation approach for a uniform wavetrain, is shown to be distinct from an amplitude perturbation expansion about the periodic flow. In fact we show that the amplitude equation contains only linear terms and is simply the heat equation. We review, briefly, the well known dynamics of Burgers equation, which imply that both shock structures and finite time singularities of the wavenumber perturbation can occur with respect to the slow scales. Numerical computations have been performed to identify areas of the (wavenumber, Reynolds number, energy) neutral surface for which each of these possibilities can occur. We note that the evolution equations will breakdown under certain circumstances, in particular for a weakly nonlinear secondary flow. Finally we extend the theory to three dimensions and discuss the limit of a weak spanwise dependence for uniform wavetrains, showing that two functions are required to describe the evolution. These unknowns are a phase and a pressure function which satisfy a pair of linearly coupled partial differential equations. The results obtained from applying the same analysis to the fully three dimensional problem are included as an appendix.