Search NASASearch

Engineering topics

Wray, A.

Publications and source records attributed to Wray, A..

High Order Numerical Methods for LES of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten's subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.

High Order

MHD Turbulence Sheared in Fixed and Rotating Frames

We consider homogeneous turbulence in a conducting fluid that is exposed to a uniform external magnetic field while being sheared in fixed and rotating frames. We take both the frame-rotation axis and the applied magnetic field to be aligned in the direction normal to the plane of the mean shear. Here a systematic parametric study is carried out in a series of Direct Numerical Simulations (DNS) in order to clarify the main effects determining the structural anisotropy and stability of the flow.

Kassinos, S. C.

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.

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

Ensemble averaged dynamic modeling

The possibility of using the information from simultaneous equivalent Large Eddy Simulations (LAS) for improving the subgrid scale modeling is investigated. An ensemble average dynamic model is proposed as an alternative to the usual spatial average versions. It is shown to be suitable independently of the existence of any homogeneity directions, and its formulation is thus universal. The ensemble average dynamic model is shown to give very encouraging results for as few as 16 simultaneous LES's.

Carati, D.

A new numerical method for the simulation of three dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction.

Leonard, A.

Optical studies of shock generated transient supersonic base flows

A shock tube employing interferometric and schlieren techniques is used to study transient base flow phenomena following shock wave passage over two plane bluff bodies: a hemicircular cylinder and a cylinder with the Galileo Jovian probe profile. An attempt is made to understand the physics of transition from transient to steady state flow, and to provide code verification for a study employing the Illiac IV computer. Transient base flow interactions include a series of shock diffraction, regular, and Mach reflections, coupled with boundary layer development, separation, and recompression. Vorticity generation and transport underlie these features. The quantitative verification of the computer code includes comparisons of transient pressure and density fields, near wake geometries, and bow shock standoff distances.

Liang, P.-Y.

A new numerical method for the simulation of three-dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction. Previously announced in STAR as N82-31644

Leonard, A.

Numerical experiments in boundary-layer stability

Numerical solution of the three-dimensional incompressible Navier-Stokes equations is used to study the instability of a flat-plate boundary layer in a manner analogous to the vibrating-ribbon experiments. Flow-field structures are observed which are very similar to those found in the vibrating-ribbon experiment to which computational initial conditions have been matched. Streamwise periodicity is assumed in the simulation so that the evolution occurs in time, but the events which constitute the instability are so similar to the spatially occurring ones of the laboratory that it seems clear the physical processes involved are the same. A spectral and finite difference numerical algorithm is employed in the simulation.

Wray, A.