Search NASASearch

Engineering topics

Merriam, Marshal

Publications and source records attributed to Merriam, Marshal.

Structure of the Compressible Reacting Mixing Layer: A Linear Stability Analysis

Previous investigations have demonstrated that a mixing layer in compressible reacting conditions can develop two instability modes in addition to the more common central mode, which exists unaccompanied in incompressible non-reacting flows. These 'outer' modes were shown in the investigations of Planche and Reynolds to have a detrimental effect on mixing and combustion efficiency and therefore their presence has important consequences for applications of supersonic combustion. Results from a parametric study of the compressible reacting mixing layer's regime space using a linear stability technique will be discussed. The purpose is to develop a more complete understanding behind the effects of compressibility, heat release and the density, equivalence and velocity ratios on the amplification rate and structure of each instability mode.

Day, M. J.

Numerical Methods Using B-Splines

The seminar will discuss (1) The current range of applications for which B-spline schemes may be appropriate (2) The property of high-resolution and the relationship between B-spline and compact schemes (3) Comparison between finite-element, Hermite finite element and B-spline schemes (4) Mesh embedding using B-splines (5) A method for the incompressible Navier-Stokes equations in curvilinear coordinates using divergence-free expansions.

Shariff, Karim

The Use of DNS in Turbulence Modeling

The use of Direct numerical simulations (DNS) data in developing and testing turbulence models is reviewed. The data is used to test turbulence models at all levels: algebraic, one-equation, two-equation and full Reynolds stress models were tested. Particular examples on the development of models for the dissipation rate equation are presented. Homogeneous flows are used to test new scaling arguments for the various terms in the dissipation rate equation. The channel flow data is used to develop modifications to the equation model that take into account near-wall effects. DNS of compressible flows under mean compression are used in testing new compressible modifications to the two-equation models.

Mansour, Nagi N.

Direct Numerical Simulation of Supersonic Turbulent Boundary Layers

The talk will present some initial results from the direct numerical simulation (DNS) of compressible turbulent boundary layers. We solve numerically the compressible Navier-Stokes equations using a method based on Spalart's transformation for the incompressible turbulent boundary layer. This allows the spatially developing boundary layer to be transformed to a calculation with periodic boundary conditions in the streamwise and spanwise directions. The equations are solved using Fourier expansions in the horizontal directions and B-splines in the wall-normal direction. The first simulation is at Mach 2.5 with a momentum thickness Reynolds number based on wall viscosity of R(sub theta(sup 1)) = 825. We are examining the physics of the compressible boundary layer using turbulence statistics and budget equations. The turbulence statistics include: rms (root mean square) and mean profiles, energy spectra, and two-point correlations. It is found that there are large density gradients which require significantly more resolution than the incompressible case.

Guarini, Stephen

A Model For the Limiting Time in Vortex Ring Formation

In another presentation, Gharib et at provide experimental results to show that when a vortex ring is created from a pipe by a piston, there is a limiting time or piston stroke length beyond which multiple rings appear. This time appeared to be insensitive to piston velocity history and Reynolds number. Nature might exploit such a. limit in different contexts to coherently deliver mass or momentum flux with the least number of strokes. Here, a simple hypothesis is considered: the limiting time occurs when the apparatus is no longer able to deliver energy at a rate compatible with the requirement, due to Kelvin, that a steady vortex ring have maximum energy given circulation and impulse. More specifically, the limit is expected to occur when the quantity alpha = E/square root of Gamma(sup 3)I delivered by the piston drops below the value, alpha(sub lim) for a limiting steady vortex ring solution. The resulting predictions agree very well with the experiments (after using alpha(sub lim) measured using the experimental flow fields). The insensitivity to piston history also emerges from the model. Finally, piston histories are designed that may extend the limiting time somewhat.

Shariff, Karim

Reynolds Stress Balance in Plane Wakes Subjected to Irrotational Strains

Direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of various irrotational plane strains have been generated. A pseudospectral numerical method with up to 25 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, R(sub e), of about 2,000. All the terms in the equations governing the evolution of the Reynolds stresses have been calculated. The relative importance of the various terms is examined for the different strain geometries and the behavior of the individual terms is used to better assess whether the strained wakes are evolving self-similarly.

Rogers, Miichael M.

The Effect of Strain Rate on the Evolution of Plane Wakes Subjected to Irrotational Strains

Direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of irrotational plane strain applied at three different strain rates have been generated. The strain geometry is such that the flow is compressed in the streamwise direction and expanded in the cross-stream direction with the spanwise direction being unstrained. This geometry is the temporally evolving analogue of a spatially evolving wake in an adverse pressure gradient. A pseudospectral numerical method with up to 16 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, Re, of about 2,000. Although the evolutions of many statistics are nearly collapsed when plotted against total strain, there are some differences owing to the different strain rate histories. The impact of strain-rate on the wake spreading rate, the peak velocity deficit, the Reynolds stress profiles, and the flow structure is examined.

Rogers, Michael M.

Uncertainty Analysis for Oil-Film Interferometry Skin-Friction Measurement Techniques

Over the past 20 years, the use of oil-film interferometry to measure the skin friction coefficient (C(sub f) = tau/q where tau is the surface shear stress and q is the dynamic pressure) has increased. Different forms of this oil-film technique with various levels of accuracy and ease of use have been successfully applied in a wide range of flows. The method's popularity is growing due to its relative ease of implementation and minimal intrusiveness as well as an increased demand for C(sub f) measurements. Nonetheless, the accuracy of these methods has not been rigorously addressed to date. Most researchers have simply shown that the skin-friction measurements made using these techniques compare favorably with other measurements and theory, most of which are only accurate to within 5-20%. The use of skin-friction data in the design of commercial aircraft, whose drag at cruise is 50% skin-friction drag, and in the validation of computational fluid dynamics programs warrants better uncertainty estimates. Additional information is contained in the original extended abstract.

Naughton, Jonathan W.

Some Aspects of Nonlinear Dynamics and CFD

The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with examples of spurious behavior observed in CFD computations.

Yee, Helen C.

Ring Bubbles of Dolphins

The article discusses how dolphins create and play with three types of air-filled vortices. The underlying physics is discussed. Photographs and sketches illustrating the dolphin's actions and physics are presented. The dolphins engage in this behavior on their own initiative without food reward. These behaviors are done repeatedly and with singleminded effort. The first type is the ejection of bubbles which, after some practice on the part of the dolphin, turn into toroidal vortex ring bubbles by the mechanism of baroclinic torque. These bubbles grow in radius and become thinner as they rise vertically to the surface. One dolphin would blow two in succession and guide them to fuse into one. Physicists call this a vortex reconnection. In the second type, the dolphins first create an invisible vortex ring in the water by swimming on their side and waving their tail fin (also called flukes) vigorously. This vortex ring travels horizontally in the water. The dolphin then turns around, finds the vortex and injects a stream of air into it from its blowhole. The air "fills-out" the core of the vortex ring. Often, the dolphin would knock-off a smaller ring bubble from the larger ring (this also involves vortex reconnection) and steer the smaller ring around the tank. One other dolphin employed a few other techniques for planting air into the fluke vortex. One technique included standing vertically in the water with tail-up, head-down and tail piercing the free surface. As the fluke is waved to create the vortex ring, air is entrained from above the surface. Another technique was gulping air in the mouth, diving down, releasing air bubbles from the mouth and curling them into a ring when they rose to the level of the fluke. In the third type, demonstrated by only one dolphin, the longitudinal vortex created by the dorsal fin on the back is used to produce 10-15 foot long helical bubbles. In one technique she swims in a curved path. This creates a dorsal fin vortex since centrifugal force has to be balanced by a lift-like force. She then re-traces her path and injects air into the vortex from her blowhole. She can even make a ring reconnect from the helix. In the second technique, demonstrated a few times, she again swims in a curved path, releases a cloud or group of bubbles from her blowhole and turns sharply away (Which presumably strengthens the vortex). As the bubbles encounter the vortex, they travel to the center of the vortex, merge and, in a flash, elongate along the core of the vortex. In all the three types, the air-water interface is shiny smooth and stable because the pressure gradient in the vortex flow around the bubble stabilizes it. A lot of the interesting physics still remains to be explored.

Shariff, Karim

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. we consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength Of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies.

Mansour, N. N.

A Numerical Experiment on the Role of Surface Shear Stress in the Generation of Sound

The sound generated due to a localized flow over an infinite flat surface is considered. It is known that the unsteady surface pressure, while appearing in a formal solution to the Lighthill equation, does not constitute a source of sound but rather represents the effect of image quadrupoles. The question of whether a similar surface shear stress term constitutes a true source of dipole sound is less settled. Some have boldly assumed it is a true source while others have argued that, like the surface pressure, it depends on the sound field (via an acoustic boundary layer) and is therefore not a true source. A numerical experiment based on the viscous, compressible Navier-Stokes equations was undertaken to investigate the issue. A small region of a wall was oscillated tangentially. The directly computed sound field was found to to agree with an acoustic analogy based calculation which regards the surface shear as an acoustically compact dipole source of sound.

Shariff, Karim

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise voracity is generated and ejected due to the oscillatory motion of the lid. The present simulations c Implement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University. Jacobson and Reynolds used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. More recently, Lachowiez and Wlezien are investigating the flow generated by an electro-mechanically driven lid to be used for assertion control in aerodynamic applications. We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin. The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands ol'particle's allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies. Estimates of the forces needed to drive the devices will also be presented.

Mansour, N. N.

3-D CFD in a day - The laser digitizer project

The computation of airflow over complex configurations requires a complete description of the geometry. This can be obtained from CAD data, from blueprints, or from actual models. In any case, the time required is currently estimated at 4 to 6 months. It is proposed to shorten this time by a factor of 10 to 100 through the use of automated software, a fast, highly parallel computer and a three-dimensional laser digitizer. This device can provide (x,y,z) coordinates of surface points at rates exceeding 14,500/sec. Thus, it is possible to digitize an entire model in a few minutes. The accuracy of measurement on a flat white surface is better than 0.005 inches. Higher accuracy is available at higher cost. This work discusses the challenges which remain to be addressed. In particular, the surface point data need to be converted into a surface description, the surface description needs to be made into a surface grid, and the surface grid used to make a volume grid for the flow solver. Algorithms are kept in place or in mind for all of these problems. Integration of the more mature flow solution and visualization algorithms then allows generation of solution graphics directly from a wind tunnel model.

Merriam, Marshal

Parallel unstructured grid generation

A parallel unstructured grid generation algorithm is presented and implemented on the Hypercube. Different processor hierarchies are discussed, and the appropraite hierarchies for mesh generation and mesh smoothing are selected. A domain-splitting algorithm for unstructured grids which tries to minimize the surface-to-volume ratio of each subdomain is described. This splitting algorithm is employed both for grid generation and grid smoothing. Results obtained on the Hypercube demonstrate the effectiveness of the algorithms developed.

Loehner, Rainald