Engineering topics

Shariff, K.

Publications and source records attributed to Shariff, K..

10 records

NASA NTRS2004

Direct Numerical Simulation of Turbulent Condensation in Clouds

In this brief, we investigate the turbulent condensation of a population of droplets by means of a direct numerical simulation. To that end, a coupled Navier-Stokes/Lagrangian solver is used where each particle is tracked and its growth by water vapor condensation is monitored exactly. The main goals of the study are to find out whether turbulence broadens the droplet size distribution, as observed in in situ measurements. The second issue is to understand if and for how long a correlation between the droplet radius and the local supersaturation exists for the purpose of modeling sub-grid scale microphysics in cloud-resolving codes. This brief is organized as follows. In Section 2 the governing equations are presented, including the droplet condensation model. The implementation of the forcing procedure is described in Section 3. The simulation results are presented in Section 4 together with a sketch of a simple stochastic model for turbulent condensation. Conclusions and the main outcomes of the study are given in Section 5.

Shariff, K.↗

NASA NTRS2000

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.↗

NASA NTRS1997

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↗

NASA NTRS1996

Simulation and modeling of the elliptic streamline flow

Direct numerical simulations are performed for the elliptic streamline flow, which is a homogeneous turbulent flow that combines the effects of solid body rotation and strain. Simulations are run over a range of parameters in order to determine the effect of changing rotation and strain separately. For early times the nonlinear cascade is suppressed, but then is re-established at later times. The growth rate of turbulent kinetic energy agrees at early times with the trends from linear theory, but at later times the flow seems to approach an asymptotic state that is independent of the ratio of mean flow rotation rate to strain rate. A comparison with standard Reynolds stress turbulence models is made. It is found that for strong rotation rates, the models predict decay of the turbulence, while the simulations show exponential growth. Close examination of the simulation results shows that they are affected by excessively low Reynolds numbers. Suggestions for reducing low Reynolds number effects in future simulations is given.

Blaisdell, G. A.↗

NASA NTRS1995

Numerical Experiments on Homogeneous Strained Turbulence Subjected to Coriolis Force

Homogeneous turbulent flows with various combinations of strain-rate, rotation rate and coriolis force capture some important aspects of more complex flows with streamline curvature and rotation. Presently, a situation is considered in which as a box of turbulence rotates, strain axes rotate with it. This is to be contrasted with the elliptic streamline flow in which the strain axes are fixed in an inertial frame. The elliptic flow is known to exhibit (inviscid) growth of turbulent energy and one might expect even more rapid growth with the strain-axes following the box. Instead, it is found that the sign of the Reynolds shear stress is reversed leading to a negative production term for turbulent energy. Partial understanding of the phenomenon is obtained from a consideration of the rotation of inertial waves relative to the strain axes as well as the "pressure-less" RDT argument put forward by Cambon etal. [J. Fluid Mech, 278, 175]. Some comparisons with the predictions of second-order closure models will be presented.

Shariff, K.↗

NASA NTRS1994

Vortex ring instability and its sound

This work carries earlier finite-difference calculations of the Widnall instability of vortex rings into the late non-linear stage. Plots of energy in azimuthal Fourier modes indicate that low-order modes dominate at large times; their structure and dynamics remain unexplored, however. An attempt was made to calculate the acoustic signal using the theory of Mohring (1978), valid for unbounded flow. This theory shows that only low-order azimuthal modes contribute to the sound. As a check on the effects of axial periodicity and a slip wall at large radius imposed by the numerical scheme, the acoustic integrals were also computed in a truncated region. Half of the terms contributing to the sound have large differences between the two regions, and the results are therefore unreliable. The error is less severe for a contribution involving only the m = 2 mode, and its low frequency is consistent with a free elliptic bending wave on a thin ring.

Verzicco, R.↗

NASA NTRS1994

Homogeneous turbulence subjected to mean flow with elliptic streamlines

Direct numerical simulations are performed for homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. Qualitative comparisons are made with linear theory for cases with high Rossby number. The nonlinear transfer process is monitored using a generalized skewness. In general, rotation turns off the nonlinear cascade; however, for moderate ellipticities and rotation rates the nonlinear cascade is turned off and then reestablished. Turbulence statistics of interest in turbulence modeling are calculated, including full Reynolds stress budgets.

Blaisdell, G. A.↗

NASA NTRS1993

Instabilities of two-dimensional inviscid compressible vortices

An investigation of the linear stability and subsequent nonlinear evolution and acoustic radiation of a planar inviscid compressible vortex is presented. The effects of the entropy gradient are investigated, and for the particular entropy profile chosen, the positive average entropy in the vortex core is destabilizing, while the opposite is true for the negative average entropy gradient. Finite-difference methods are used to study the linear initial value problem. These methods are found to be capable of accurately computing the frequencies and weak growth rates of the normal modes. When the initial condition consists of random perturbations, the long-time behavior is found to correspond to the most unstable normal mode in all cases. The numerical solution of the Euler equations is used to study the nonlinear development of an elliptic-mode perturbation.

Chan, W. M.↗

NASA NTRS1988

Head-on collision of viscous vortex rings

The head-on collision of two identical axisymmetric viscous vortex rings is studied through direct simulations of the incompressible Navier-Stokes equations. The initial vorticity distributions considered are those of Hill's spherical vortex and of rings with circular Gaussian cores, each at Reynolds numbers of about 350 and 1000. The Reynolds number is defined by Gamma/Nu, the ratio of circulation to viscosity. As the vortices approach each other by self-induction, the radii increase by mutual induction, and vorticy cancels through viscous cross-diffusion across the collision plane. Following contact, the vorticity distribution in the core forms a head-tail structure (for the cases considered). The characteristic time of vorticity annihilation is compared with that of a 3D collision experiment and 3D numerical simulations. It is found that the annihilation time is somewhat longer in the axisymmetric case than it is in the symmetry plane of the experiment and 3D numerical simulation. By comparing the annihilatiom time with a viscous timescale and a circulation timescale, it is deduced that both the strain rate due to local effects and to 3D vorticity realignment are important.

Stanaway, S.↗

NASA NTRS1987

Evaluation of a theory for pressure-strain rate

A theoretical expression for the slow part (the nonlinear fluctuation part) of the pressure-strain rate is compared with simulations of anisotropic homogeneous flows. The objective is to determine the quantitative accuracy of the theory and to test its prediction that the generalized Rotta coefficient, a non-dimensionalized ratio of slow term to the Reynolds stress anisotropy, varies with direction and can be negative. Comparisons are made between theoretical and simulation values of the slow term itself and of the generalized Rotta coefficients. The implications of the comparison for two-point closure theories and for Reynolds stress modeling are pointed out.

Weinstock, J.↗