Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixing”

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 451 records · Page 25

Determination of stress intensity factors for interface cracks under mixed-mode loading

A simple technique was developed using conventional finite element analysis to determine stress intensity factors, K1 and K2, for interface cracks under mixed-mode loading. This technique involves the calculation of crack tip stresses using non-singular finite elements. These stresses are then combined and used in a linear regression procedure to calculate K1 and K2. The technique was demonstrated by calculating three different bimaterial combinations. For the normal loading case, the K's were within 2.6 percent of an exact solution. The normalized K's under shear loading were shown to be related to the normalized K's under normal loading. Based on these relations, a simple equation was derived for calculating K1 and K2 for mixed-mode loading from knowledge of the K's under normal loading. The equation was verified by computing the K's for a mixed-mode case with equal and normal shear loading. The correlation between exact and finite element solutions is within 3.7 percent. This study provides a simple procedure to compute K2/K1 ratio which has been used to characterize the stress state at the crack tip for various combinations of materials and loadings. Tests conducted over a range of K2/K1 ratios could be used to fully characterize interface fracture toughness.

Naik, Rajiv A.↗

An experimental study of secondary vortex structure in mixing layers

This report covers the first eight months of an experimental research project on the secondary vortex structure in plane mixing layers. The aim of the project is to obtain quantitative data on the behavior of the secondary structure in a turbulent mixing layer at reasonable reynolds numbers (Re(sub delta(sub w)) approx. 50,000). In particular, we hope to resolve the questions of how the scale of the secondary vortex structure changes with the scale of the mixing layer, and whether the structures are fixed in space, or whether they 'meander' in the spanwise direction.

Bell, J. H.↗

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗

Lidar observation of marine mixed layer

Marine mixed layer is known to play an important role in the transportation of pollution exiting ship funnels. The application of a diffusion model is critically dependent upon a reliable estimate of a lid. However, the processes that form lids are not well understood, though considerable progress toward marine boundary layer has been achieved. This report describes observations of the marine mixed layer from the course Ise-wan to Nii-jima with the intention of gaining a better understanding of their structure by a shipboard lidar. These observations were made in the summer of 1991. One interesting feature of the observations was that the multiple layers of aerosols, which is rarely numerically modeled, was encountered. No attempt is yet made to present a systematic analysis of all the data collected. Instead we focus on observations that seem to be directly relevant to the structure of the mixed layer.

Yamagishi, Susumu↗

Sensitivity of the structure of untripped mixing layers to small changes in initial conditions

An experimental study was conducted concerning the influence of small changes in initial conditions on the near- and far-field evolution of the three-dimensional structure of a plan mixing layer. A two-stream mixing layer with a velocity ratio of 0.6 was generated with the initial boundary layers on the splitter plate laminar and was nominally two-dimensional. The initial conditions were changed slightly by interchanging the high- and low-speed sides of the wind tunnel, while maintaining the same velocities, and hence velocity ratio. This resulted in small changes in the initial boundary layer properties, and the perturbations present in the boundary layers were interchanged between the high- and low-speed sides for the two cases. The results indicate that, even with this relatively minor change in initial conditions, the near-field regions of the two cases differ significantly. The peak Reynolds stress levels in the near-field differ by up to 100 percent, and this is attributed to a difference in the location of the initial spanwise vortex roll-up. In addition, the positions and shapes of the individual streamwise vortical structures differ for the two cases, although the overall structures differ for the two cases, although the overall qualitative description of these structures is comparable. The subsequent reorganization and decay of the streamwise vortical structures is very similar for the two cases. As a result, in the far field, both mixing layers achieve similar structure, yielding comparable growth rates, Reynolds stress, distribution, and spectral content.

Plesniak, M. W.↗

Effects of unsteady shock impingement on high-speed gaseous mixing

Tests are performed to determine the effects of impinging oscillating shocks of different frequencies on a 15 deg downstream angled, underexpanded, sonic helium jet injected into a supersonic airflow. Information on mixing, penetration, total pressure loss, and turbulence structure from these tests is employed to estimate mixing control achieved by adding an oscillating shock to the helium injection flow field. The principal result of this study is that impingement of an oscillating shock on a high-speed shear layer can be utilized to control the rate of mixing.

Wood, C. W.↗

An assumed pdf approach for the calculation of supersonic mixing layers

In an effort to predict the effect that turbulent mixing has on the extent of combustion, a one-equation turbulence model is added to an existing Navier-Stokes solver with finite-rate chemistry. To average the chemical-source terms appearing in the species-continuity equations, an assumed pdf approach is also used. This code was used to analyze the mixing and combustion caused by the mixing layer formed by supersonic coaxial H2-air streams. The chemistry model employed allows for the formation of H2O2 and HO2. Comparisons are made with recent measurements using laser Raman diagnostics. Comparisons include temperature and its rms, and concentrations of H2, O2, N2, H2O, and OH. In general, good agreement with experiment was noted.

Baurle, R. A.↗

Model free simulations of a high speed reacting mixing layer

The effects of compressibility, chemical reaction exothermicity and non-equilibrium chemical modeling in a combusting plane mixing layer were investigated by means of two-dimensional model free numerical simulations. It was shown that increased compressibility generally had a stabilizing effect, resulting in reduced mixing and chemical reaction conversion rate. The appearance of 'eddy shocklets' in the flow was observed at high convective Mach numbers. Reaction exothermicity was found to enhance mixing at the initial stages of the layer's growth, but had a stabilizing effect at later times. Calculations were performed for a constant rate chemical rate kinetics model and an Arrhenius type kinetics prototype. The Arrhenius model was found to cause a greater temperature increase due to reaction than the constant kinetics model. This had the same stabilizing effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗

Structure of a swirling turbulent mixing layer

A single-stream swirling turbulent mixing layer was produced in which the angular momentum instability acts over most of the layer for most of its development length. The instability is associated with a large increase in the level of the conventional Reynolds stresses over those in an unswirled mixing layer. This paper presents measurements of the three-dimensionality of the mean flow at one axial location well downstream of the origin. The output of a fixed hot-wire probe was ensemble-averaged to produce an approximation to the mean flow field that would be seen by a rotating observer. The main goals of this study were to establish the presence or absence of Taylor-Goertler vortices found in unstable laminar boundary layers. The variations in the mean velocities and Reynolds stresses are generally consistent with a bodily 'wrinkling' of the mixing layer, without a significant variation in the layer thickness. The wrinkling is most probably a result of circumferential variations in the initial boundary layer, which in turn originate at the swirl generator.

Wood, D. H.↗

Measurements of the streamwise vortical structures in a plane mixing layer

The 3D structure of a plane two-stream mixing layer of velocity ratio 0.6 and originating from laminar initial boundary layers was investigated through direct measurements made in a specially constructed mixing-layer wind tunnel. The main objective of the study was to establish quantitatively the presence and the role of the secondary streamwise vortex structure (of the kind that has been shown in past flow visualization investigations to ride among the primary spanwise vortices) in the development of a plane turbulent mixing layer at relatively high Reynolds numbers. Results indicate that the instability leading to the formation of streamwise vortices is initially amplified just downstream of the first spanwise roll-up. The streamwise vortices, which first appear in clusters containing vorticity of both signs, realign further downstream to form counterrotating pairs. Due to the amalgamation of like-sign vortices, the streamwise vortex spacing increases in a stepwise fashion.

Bell, James H.↗

Vorticity generation mechanisms in parallel injection schemes for supersonic mixing

A numerical study is reported of the three-dimensional nonreacting supersonic flow field produced by three parallel fuel injection schemes. Such injections are being considered as a means of enhancing fuel mixing in future Scramjet engine designs. The strength of vortex structures produced by ramp injectors and by rocket/jet interactions are quantified and their relative effectiveness in enhancing the fuel mixing process is addressed. An experimental setup for validation of the above numerical results has been constructed and preliminary results are presented. For the flow field parameters chosen, the ramp generated vorticity is found to be considerably larger than that generated by the shock. The unsteady recirculating flow field in the wake of the injector appears also to give rise to important additional vorticity and mixing.

Donohue, James M.↗

A numerical study of Li-SF6 wick combustion - Forced and mixed convective burning

A numerical study is conducted to study Li-SF6 wick diffusion flames under mixed convective burning conditions at a pressure of 0.01 MPa. Both planar and cylindrical wicks are considered. The model is based on a conserved scalar approach. The objective of this study is to assess the effects of particular parameters on the burning rate and heat transfer. The flat-plate solution yields a fuel mass burning rate per unit surface area following the x exp -1/2 dependence of the classical similarity solution, where x is the streamwise distance. Cylindrical wick geometries yield enhanced burning rates over planar wicks. For the case of mixed convective burning, the burning rate results approach either the forced or natural convective burning limits as ambient streamwise velocity is changed. Critical Richardson numbers specifying these burning limits are determined for a given condition. Reducing gravity results in a lower burning rate because the influence of natural convection is diminished. Under reduced gravity of 1/1000 of the sea-level value, mixed convective burning nearly resembles forced convection.

Damaso, R. C.↗

An equivalent domain integral method for three-dimensional mixed-mode fracture problems

A general formulation of the equivalent domain integral (EDI) method for mixed mode fracture problems in cracked solids is presented. The method is discussed in the context of a 3-D finite element analysis. The J integral consists of two parts: the volume integral of the crack front potential over a torus enclosing the crack front and the crack surface integral due to the crack front potential plus the crack face loading. In mixed mode crack problems the total J integral is split into J sub I, J sub II, and J sub III representing the severity of the crack front in three modes of deformations. The direct and decomposition methods are used to separate the modes. These two methods were applied to several mixed mode fracture problems, were analyzed, and results were found to agree well with those available in the literature. The method lends itself to be used as a post-processing subroutine in a general purpose finite element program.

Shivakumar, K. N.↗

Numerical study on mixing of sprayed liquid in an LNG storage tank

This paper presents a numerical method to simulate the mixing of heavier LNG sprayed on lighter layer. Numerical results for evolutions of flow field and density field are obtained in a rectangular computational domain which includes the vicinity of the liquid surface. At the surface boundary, uniform distributions of the fluid velocity and the density are assumed. Detail structure of flow caused by impingements of liquid drops are neglected. But, to trigger a realistic motion, a series of random numbers is employed. It is used as an initial distribution of the density near the surface. This method successfully gives a realistic simulation of the mixing process. Numerical results for mixing velocity shows good agreement with experimental data.

Uchida, Hiroyuki↗

Mixed-valence iron minerals on Venus: Fe(2+)-Fe(3+) oxides and oxy-silicates formed by surface-atmosphere interactions

Inferences from these investigations are that Fe(3+)-bearing minerals such as hematite magnesioferrite, acmite, and epidote are thermodynamically unstable, and that magnetite is the predominant mixed-valence iron oxide mineral on venus. Recently, the Fe(2+)-Fe(3+) silicate mineral laihunite was proposed to be a reaction product of olivine with the venusian atmosphere. This possibility is discussed further here. We suggest that other mixed-valence Fe(2+)-Fe(3+)-Oz-OH(-) silicates could also result from surface-atmosphere interactions on Venus. Topics discussed include the following: (1) conversion of hematite to magnetite; (2) stability of laihunite; (3) the possible existence of oxy-amphiboles and oxy-micas on Venus; and (4) other mixed-valence Fe(2+)-Fe(3+) silicates likely to exist on Venus.

Burns, Roger G.↗

Stability of mixing layers

The research program for the first year of this project (see the original research proposal) consists of developing an explicit marching scheme for solving the parabolized stability equations (PSE). Performing mathematical analysis of the computational algorithm including numerical stability analysis and the determination of the proper boundary conditions needed at the boundary of the computation domain are implicit in the task. Before one can solve the parabolized stability equations for high-speed mixing layers, the mean flow must first be found. In the past, instability analysis of high-speed mixing layer has mostly been performed on mean flow profiles calculated by the boundary layer equations. In carrying out this project, it is believed that the boundary layer equations might not give an accurate enough nonparallel, nonlinear mean flow needed for parabolized stability analysis. A more accurate mean flow can, however, be found by solving the parabolized Navier-Stokes equations. The advantage of the parabolized Navier-Stokes equations is that its accuracy is consistent with the PSE method. Furthermore, the method of solution is similar. Hence, the major part of the effort of the work of this year has been devoted to the development of an explicit numerical marching scheme for the solution of the Parabolized Navier-Stokes equation as applied to the high-seed mixing layer problem.

Tam, Christopher↗

The three-dimensional evolution of a plane mixing layer. Part 2: Pairing and transition to turbulence

The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings was simulated numerically. Initial conditions for all simulations consisted of a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found both to inhibit the growth of infinitesimal three-dimensional disturbances and to trigger the transition to turbulence in highly three dimensional flows. The mechanisms responsible for the growth of three-dimensionality as well as the initial phases of the transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of span wise vorticity, which undergo a secondary roll up. Transition also produces an increase in the degree of scalar mixing, in agreement with experimental observations of mixing transition. Simulations were also conducted to investigate changes in span wise length scale that may occur in response to the change in stream wise length scale during a pairing. The linear mechanism for this process was found to be very slow, requiring roughly three pairings to complete a doubling of the span wise scale. Stronger three-dimensionality can produce more rapid scale changes but is also likely to trigger transition to turbulence. No evidence was found for a change from an organized array of rib vortices at one span wise scale to a similar array at a larger span wise scale.

Moser, Robert D.↗

Effects of material mixing on planetary gamma ray spectroscopy

The effect of mixing of rocks and soil on the gamma-ray signal emitted by Martian surface was investigated with special consideration given to the effect of the geometry of mixing. It is shown that, under certain circumstances, this scale can have a significant effect on the gamma-ray spectrum observed. The techniques that can be used to address this problem are discussed, showing that, using these techniques, rock compositions can be inferred from the signal produced by a rock-soil mixtures of unknown rock concentration and mixing geometry.

Squyres, Steven W.↗