Search NASASearch

Engineering topics

Shariff, Karim

Publications and source records attributed to Shariff, Karim.

Detached-Eddy Simulation Based on the v2-f Model

Detached eddy simulation (DES) based on the v2-f RANS model is proposed. This RANS model incorporates the anisotropy of near-wall turbulence which is absent in other RANS models commonly used in the DES community. In LES mode, the proposed DES formulation reduces to a transport equation for the subgrid-scale kinetic energy. The constant, CDES, required by this model was calibrated by simulating isotropic turbulence. In the final paper, DES simulations of canonical separated flows will be presented.

Jee, Sol Keun

LES of High-Reynolds-Number Coanda Flow Separating from a Rounded Trailing Edge of a Circulation Control Airfoil

This slide presentation reviews the Large Eddy Simulation of a high reynolds number Coanda flow that is separated from a round trailing edge of a ciruclation control airfoil. The objectives of the study are: (1) To investigate detailed physics (flow structures and statistics) of the fully turbulent Coanda jet applied to a CC airfoil, by using LES (2) To compare LES and RANS results to figure out how to improve the performance of existing RANS models for this type of flow.

Nichino, Takafumi

The force exerted by the membrane potential during protein import into the mitochondrial matrix

The force exerted on a targeting sequence by the electrical potential across the inner mitochondrial membrane is calculated on the basis of continuum electrostatics. The force is found to vary from 3.0 pN to 2.2 pN (per unit elementary charge) as the radius of the inner membrane pore (assumed aqueous) is varied from 6.5 to 12 A, its measured range. In the present model, the decrease in force with increasing pore width arises from the shielding effect of water. Since the pore is not very much wider than the distance between water molecules, the full shielding effect of water may not be present; the extreme case of a purely membranous pore without water gives a force of 3.2 pN per unit charge, which should represent an upper limit. When applied to mitochondrial import experiments on the protein barnase, these results imply that forces between 11 +/- 2 pN and 13.5 +/- 2.5 pN catalyze the unfolding of barnase in those experiments. A comparison of these results with unfolding forces measured using atomic force microscopy is made.

Models, Molecular

Particle Size Distributions in Atmospheric Clouds

In this note, we derive a transport equation for a spatially integrated distribution function of particles size that is suitable for sparse particle systems, such as in atmospheric clouds. This is done by integrating a Boltzmann equation for a (local) distribution function over an arbitrary but finite volume. A methodology for evolving the moments of the integrated distribution is presented. These moments can be either tracked for a finite number of discrete populations ('clusters') or treated as continuum variables.

Paoli, Roberto

The Speed of Axial Propagation of a Cylindrical Bubble Through a Cylindrical Vortex

Inspired by the rapid elongation of air columns injected into vortices by dolphins, we present an exact inviscid solution for the axial speed (assumed steady) of propagation of the tip of a semi-infinite cylindrical bubble along the axis of a cylindrical vortex. The bubble is assumed to be held at constant pressure by being connected to a reservoir, the lungs of the dolphin, say. For a given bubble pressure, there is a modest critical rotation rate above which steadily propagating bubbles exist. For a bubble at ambient pressure, the propagation speed of the bubble (relative to axial velocity within the vortex) varies between 0.5 and 0.6 of the maximum rotational speed of the vortex. Surprisingly, the bubble tip can propagate (almost as rapidly) even when the pressure minimum in the vortex core is greater than the bubble pressure; in this case, solutions exhibit a dimple on the nose of the bubble. A situation important for incipient vortex cavitation, and one which dolphins also demonstrate, is elongation of a free bubble, i.e., one whose internal pressure may vary. Under the assumption that the acceleration term is small (checked a posteriori), the steady solution is applied at each instant during the elongation. Three types of behavior are then possible depending on physical parameters and initial conditions: (A) Unabated elongation with slowly increasing bubble pressure, and nearly constant volume. Volume begins to decrease in the late stages. (B1) Elongation with decreasing bubble pressure. A limit point of the steady solution is encountered at a finite bubble length. (B2) Unabated elongation with decreasing bubble pressure and indefinite creation of volume. This is made possible by the existence of propagating solutions at bubble pressures below the minimum vortex pressure. As the bubble stretches, its radius initially decreases but then becomes constant; this is also observed in experiments on incipient vortex cavitation.

Shariff, Karim

The Force Exerted by the Membrane Potential During Protein Import into the Mitochondrial Matrix

The electrostatic force exerted on a targeting sequence by the electrical potential across the inner mitochondrial membrane is calculated and found to vary from 1.4 pN to 2.2 pN (per unit elementary charge) as the radius of the inner membrane pore (assumed aqueous) is varied from 12 to 6.5 Angstroms, its measured range. Since the pore is not very much wider than the distance between water molecules, the full shielding effect of water may not be present; the extreme case of a nonaqueous pore gives a force of 3.1 pN per unit charge, which represents an upper limit. When applied to mitochondrial import experiments on the protein harness, these results imply that a force of 11 plus or minus 4 pN is sufficient to catalyze the unfolding of harness during import. Comparison of these results with unfolding forces measured using atomic force microscopy suggests that the two are not inconsistent.

Shariff, Karim

B-spline Method in Fluid Dynamics

B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Botella, Olivier

Non-Linear Lessons from Axisymmetric Vortex Rings

For presentation at Northwestern University, May 14-28, 2001, the talk will present two types of phenomena, both recognizable to students of nonlinearity, that are exhibited by axisymmetric vortex rings in numerical and laboratory experiments. (1) The first type of phenomenon is reminiscent of inelastic solitons.

Shariff, Karim

Analysis of the Radar Reflectivity of Aircraft Vortex Wakes

Radar has been proposed as a way to track wake vortices to reduce aircraft spacing and tests have revealed radar echoes from aircraft wakes in clear air. The results are always interpreted qualitatively using Tatarski's theory of weak scattering by isotropic atmospheric turbulence. The goal of the present work was to predict the value of the radar cross-section (RCS) using simpler models. This is accomplished in two steps. First, the refractive index is obtained. Since the structure of the aircraft wakes is different from atmospheric turbulence, three simple mechanisms specific to vortex wakes are considered: (1) Radial density gradient in a two-dimensional vortex, (2) three-dimensional fluctuations in the vortex cores, and (3) Adiabatic transport of the atmospheric fluid in a two-dimensional oval surrounding the pair of vortices. The index of refraction is obtained more precisely for the two-dimensional mechanisms than for the three-dimensional ones. In the second step, knowing the index of refraction, a scattering analysis is performed. Tatarski's weak scattering approximation is kept but the usual assumptions of a far-field and a uniform incident wave are dropped. Neither assumption is generally valid for a wake that is coherent across the radar beam. For analytical insight, a simpler approximation that invokes, in addition to weak scattering, the far-field and wide cylindrical beam assumptions, is also developed and compared with the more general analysis. The predicted RCS values for the oval surround the vortices (mechanism C) agree with the experiments of Bilson conducted over a wide range of frequencies. However, the predictions have a cut-off away from normal incidence which is not present in the measurements. Estimates suggest that this is due to turbulence in the baroclinic vorticity generated at the boundary of the oval. The reflectivity of a vortex itself (mechanism A) is comparable to that of the oval (mechanism C) but cuts-off at frequencies lower than those considered in all the experiments to date. The RCS of a vortex happens to peak at the frequency (about 49 MHz) where atmospheric radars (known as ST radars) operate and so the present prediction could be verified in the future. Finally , we suggest that hot engine exhaust could increase RCE by 40 db and reveal vortex circulation, provided its mixing with the surroundings is prevented in the laminarising flow of the vortices.

Shariff, Karim

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

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

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

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

Two-dimensional mesh embedding for Galerkin B-spline methods

A number of advantages result from using B-splines as basis functions in a Galerkin method for solving partial differential equations. Among them are arbitrary order of accuracy and high resolution similar to that of compact schemes but without the aliasing error. This work develops another property, namely, the ability to treat semi-structured embedded or zonal meshes for two-dimensional geometries. This can drastically reduce the number of grid points in many applications. Both integer and non-integer refinement ratios are allowed. The report begins by developing an algorithm for choosing basis functions that yield the desired mesh resolution. These functions are suitable products of one-dimensional B-splines. Finally, test cases for linear scalar equations such as the Poisson and advection equation are presented. The scheme is conservative and has uniformly high order of accuracy throughout the domain.

Shariff, Karim

A Study of Two-Equation Turbulence Models on the Elliptic Streamline Flow

Several two-equation turbulence models are compared to data from direct numerical simulations (DNS) of the homogeneous elliptic streamline flow, which combines rotation and strain. The models considered include standard two-equation models and models with corrections for rotational effects. Most of the rotational corrections modify the dissipation rate equation to account for the reduced dissipation rate in rotating turbulent flows, however, the DNS data shows that the production term in the turbulent kinetic energy equation is not modeled correctly by these models. Nonlinear relations for the Reynolds stresses are considered as a means of modifying the production term. Implications for the modeling of turbulent vortices will be discussed.

Blaisdell, Gregory A.

Vortex rings

The vortex-ring problem in fluid mechanics is examined generally in terms of formation, the steady state, the duration of the rings, and vortex interactions. The formation is studied by examining the generation of laminar and turbulent vortex rings and their resulting structures with attention given to the three stages of laminar ring development. Inviscid dynamics is addressed to show how core dynamics affects overall ring motion, and laminar vortex structures are described in two dimensions. Viscous and inviscid structures are related in terms of 'leapfrogging', head-on collisions, and collisions with a no-slip wall. Linear instability theory is shown to successfully describe observational data, although late stages in the breakdown are not completely understood. This study of vortex rings has important implications for key aerodynamic issues including sound generation, transport and mixing, and vortex interactions.

Shariff, Karim

Dynamics of a class of vortex rings

The contour dynamics method is extended to vortex rings with vorticity varying linearly from the symmetry axis. An elliptic core model is also developed to explain some of the basic physics. Passage and collisions of two identical rings are studied focusing on core deformation, sound generation and stirring of fluid elements. With respect to core deformation, not only the strain rate but how rapidly it varies is important and accounts for greater susceptibility to vortex tearing than in two dimensions. For slow strain, as a passage interaction is completed and the strain relaxes, the cores return to their original shape while permanent deformations remain for rapidly varying strain. For collisions, if the strain changes slowly the core shapes migrate through a known family of two-dimensional steady vortex pairs up to the limiting member of the family. Thereafter energy conservation does not allow the cores to maintain a constant shape. For rapidly varying strain, core deformation is severe and a head-tail structure in good agreement with experiments is formed. With respect to sound generation, good agreement with the measured acoustic signal for colliding rings is obtained and a feature previously thought to be due to viscous effects is shown to be an effect of inviscid core deformation alone. For passage interactions, a component of high frequency is present. Evidence for the importance of this noise source in jet noise spectra is provided. Finally, processes of fluid engulfment and rejection for an unsteady vortex ring are studied using the stable and unstable manifolds. The unstable manifold shows excellent agreement with flow visualization experiments for leapfrogging rings suggesting that it may be a good tool for numerical flow visualization in other time periodic flows.

Shariff, Karim

Acoustics and dynamics of coaxial interacting vortex rings

Using a contour dynamics method for inviscid axisymmetric flow we examine the effects of core deformation on the dynamics and acoustic signatures of coaxial interacting vortex rings. Both 'passage' and 'collision' (head-on) interactions are studied for initially identical vortices. Good correspondence with experiments is obtained. A simple model which retains only the elliptic degree of freedom in the core shape is used to explain some of the calculated features.

Shariff, Karim