Engineering Papers⌕ Search

Engineering topics

Mansour, Nagi N.

Publications and source records attributed to Mansour, Nagi N..

At least 91 records · Page 5

Performance of Several High Order Numerical Methods for Supersonic Combustion

The performance of two recently developed numerical methods by Yee et al. and Sjoegreen and Yee using postprocessing nonlinear filters is examined for a 2-D multiscale viscous supersonic react-live flow. These nonlinear filters can improve nonlinear instabilities and at the same time can capture shock/shear waves accurately. They do not, belong to the class of TVD, ENO or WENO schemes. Nevertheless, they combine stable behavior at discontinuities and detonation without smearing the smooth parts of the flow field. For the present study, we employ a fourth-order Runge-Kutta in time and a sixth-order non-dissipative spatial base scheme for the convection and viscous terms. We denote the resulting nonlinear filter schemes ACM466-RK4 and WAV66-RK4.

Sjoegreen, Bjoern↗

A New Look at Self-Similarity in Strained Plane Wakes

In early experiments, A. J. Reynolds and J. F. Keffer sought to determine whether plane wakes of circular cylinders, when strained by a wind tunnel of varying cross-section, evolved in accordance with an analytically derived self-similar solution. As pointed out by Reynolds, for the strain geometry considered this self-similar solution indicated exponential growth of the viscous term in the mean momentum equation, a result which he interpreted as suggesting that such wakes would eventually relaminarize. The experimental results were found not to agree with the similarity theory and recent direct numerical simulations confirm this. However, a more general self-similar analysis of the kind suggested by W. K. George is found to lead to families of possible similarity solutions, some of which do indeed describe the observed flaw behavior. These equilibrium similarity solutions result from creating a balance in the governing equations by grouping certain terms. For these solutions the viscous terms can be retained in the analysis.

Rogers, Michael M.↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

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

Compressing an Elliptic Vortex: Transition to Turbulence by Tumble Breakdown

Stability of the elliptic vortex attracted interest in the past decade. Cambon (1982), and Cambon, Teissedre and Jeandel (1985) have studied the stability of such flows with spatially uniform velocity gradient, and have provided RDT solutions for a wide range of the parameter S/Omega (where the strain rate S and the vorticity 2 Omega define the velocity gradient matrix). The range studied included those of hyperbolic streamlines (strain dominated, S/Omega is greater than 1), linear streamlines (simple shear, S/Omega = 1), and elliptical streamlines (vorticity dominated, S/Omega is less than 1). The latter class has more recently attracted interest and several studies appeared (Pierrehumbert 1986, Bayly 1986, Craik and Criminale 1986). These studies will be collectively referred to as PBCC.

Godeferd, Fabien S.↗

Anisotropic Structure of Rotating Homogeneous Turbulence at High Reynolds Numbers

Large eddy simulation is used to investigate the development of anisotropies and the evolution towards a quasi two-dimensional state in rotating homogeneous turbulence at high Reynolds number. The present study demonstrates the existence of two transitions in the development of anisotropy. The first transition marks the onset of anisotropy and occurs when a macro-Rossby number (based on a longitudinal integral lengthscale) has decreased to near unity while the second transition occurs when a micro-Rossby number (defined in this work as the ratio of the rms fluctuating vorticity to background vorticity) has decreased to unity. The anisotropy marked by the first transition corresponds to a reduction in dimensionality while the second transition corresponds to a polarization of the flow, i.e., relative dominance of the velocity components in the plane normal to the rotation axis. Polarization is reflected by emergence of anisotropy measures based on the two-dimensional component of the turbulence. Investigation of the vorticity structure shows that the second transition is also characterized by an increasing tendency for alignment between the fluctuating vorticity vector and the background angular velocity vector with a preference for corrotative vorticity.

Cambon, Claude↗

Energy Transfer in Rotating Turbulence

The influence or rotation on the spectral energy transfer of homogeneous turbulence is investigated in this paper. Given the fact that linear dynamics, e.g. the inertial waves regime tackled in an RDT (Rapid Distortion Theory) fashion, cannot Affect st homogeneous isotropic turbulent flow, the study of nonlinear dynamics is of prime importance in the case of rotating flows. Previous theoretical (including both weakly nonlinear and EDQNM theories), experimental and DNS (Direct Numerical Simulation) results are gathered here and compared in order to give a self-consistent picture of the nonlinear effects of rotation on tile turbulence. The inhibition of the energy cascade, which is linked to a reduction of the dissipation rate, is shown to be related to a damping due to rotation of the energy transfer. A model for this effect is quantified by a model equation for the derivative-skewness factor, which only involves a micro-Rossby number Ro(sup omega) = omega'/(2(OMEGA))-ratio of rms vorticity and background vorticity as the relevant rotation parameter, in accordance with DNS and EDQNM results fit addition, anisotropy is shown also to develop through nonlinear interactions modified by rotation, in an intermediate range of Rossby numbers (Ro(omega) = (omega)' and Ro(omega)w greater than 1), which is characterized by a marco-Rossby number Ro(sup L) less than 1 and Ro(omega) greater than 1 which is characterized by a macro-Rossby number based on an integral lengthscale L and the micro-Rossby number previously defined. This anisotropy is mainly an angular drain of spectral energy which tends to concentrate energy in tile wave-plane normal to the rotation axis, which is exactly both the slow and the two-dimensional manifold. In Addition, a polarization of the energy distribution in this slow 2D manifold enhances horizontal (normal to the rotation axis) velocity components, and underlies the anisotropic structure of the integral lengthscales. Finally is demonstrated the ability of a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) model to predict the underlying spectral transfer structure and all the subsequent developments of classic anisotropy indicators in physical space, when compared to recent LES results. Even if the applications mainly concern developed strong turbulence, a particular emphasis is given to the strong formal analogy of this EDQNM2 model with recent weakly nonlinear approaches to wave-turbulence.

Cambon, Claude↗

The asymptotic state of rotating homogeneous turbulence at high Reynolds numbers

The long-time, asymptotic state of rotating homogeneous turbulence at high Reynolds numbers has been examined using large-eddy simulation of the incompressible Navier-Stokes equations. The simulations were carried out using 128 x 128 x 512 collocation points in a computational domain that is four times longer along the rotation axis than in the other directions. Subgrid-scale motions in the simulations were parameterized using a spectral eddy viscosity modified for system rotation. Simulation results show that in the asymptotic state the turbulence kinetic energy undergoes a power-law decay with an exponent which is independent of rotation rate, depending only on the low-wavenumber form of the initial energy spectrum. Integral lengthscale growth in the simulations is also characterized by power-law growth; the correlation length of transverse velocities exhibiting much more rapid growth than observed in non-rotating turbulence.

Squires, Kyle D.↗

The effect of initial conditions on the development of temporally evolving planar three dimensional incompressible wakes

A set of direct numerical simulations of temporally evolving incompressible plane wakes started from a variety of initial conditions have been performed. The intent is to explore the effect of initial conditions on the development of three dimensionality. It was found that the existence of an oblique disturbance at a streamwise wavelength equal to the two-dimensional subharmonic wavelength plays a very important role in the development of the wake. The Reynolds number and the presence or absence of longer wavelength disturbances were also found to be important at late times in the development of the flow. Disturbance phase was found to have an effect on the details of the structure of the wake, but has relatively little effect on the wake growth.

Sondergaard, Rolf↗

A comparison of turbulence models in computing multi-element airfoil flows

Four different turbulence models are used to compute the flow over a three-element airfoil configuration. These models are the one-equation Baldwin-Barth model, the one-equation Spalart-Allmaras model, a two-equation k-omega model, and a new one-equation Durbin-Mansour model. The flow is computed using the INS2D two-dimensional incompressible Navier-Stokes solver. An overset Chimera grid approach is utilized. Grid resolution tests are presented, and manual solution-adaptation of the grid was performed. The performance of each of the models is evaluated for test cases involving different angles-of-attack, Reynolds numbers, and flap riggings. The resulting surface pressure coefficients, skin friction, velocity profiles, and lift, drag, and moment coefficients are compared with experimental data. The models produce very similar results in most cases. Excellent agreement between computational and experimental surface pressures was observed, but only moderately good agreement was seen in the velocity profile data. In general, the difference between the predictions of the different models was less than the difference between the computational and experimental data.

Rogers, Stuart E.↗

Large Eddy Simulation of Homogeneous Rotating Turbulence

Study of turbulent flows in rotating reference frames has proven to be one of the more challenging areas of turbulence research. The large number of theoretical, experimental, and computational studies performed over the years have demonstrated that the effect of solid-body rotation on turbulent flows is subtle and remains exceedingly difficult to predict. Because of the complexities associated with non-homogeneous turbulence, it is worthwhile to examine the effect of steady system rotation on the evolution of an initially isotropic turbulent flow. The assumption of statistical homogeneity considerably simplifies analysis and computation; calculation of homogeneous turbulence is further motivated since it possesses the essential physics found in more complex rotating flows. The principal objectives of the present study have therefore been to increase our fundamental understanding of turbulent flows in rotating reference frames through an examination of the asymptotic state of homogeneous rotating turbulence; particularly as to the existence of an asymptotic state which is self similar. Knowledge of an asymptotic similarity state permits prediction of the ultimate statistical evolution of the flow without requiring detailed knowledge of the complex, and not well understood, non-linear transfer processes. Aside from examination of possible similarity states in rotating turbulence, of further interest in this study has been an examination of the degree to which solid-body rotation induces a two-dimensional state in an initially isotropic flow.

Squires, Kyle D.↗

Investigation of the asymptotic state of rotating turbulence using large-eddy simulation

Study of turbulent flows in rotating reference frames has long been an area of considerable scientific and engineering interest. Because of its importance, the subject of turbulence in rotating reference frames has motivated over the years a large number of theoretical, experimental, and computational studies. The bulk of these previous works has served to demonstrate that the effect of system rotation on turbulence is subtle and remains exceedingly difficult to predict. A rotating flow of particular interest in many studies, including the present work, is examination of the effect of solid-body rotation on an initially isotropic turbulent flow. One of the principal reasons for the interest in this flow is that it represents the most basic turbulent flow whose structure is altered by system rotation but without the complicating effects introduced by mean strains or flow inhomogeneities. The assumption of statistical homogeneity considerably simplifies analysis and computation. The principal objective of the present study has been to examine the asymptotic state of solid-body rotation applied to an initially isotropic, high Reynolds number turbulent flow. Of particular interest has been to determine the degree of two-dimensionalization and the existence of asymptotic self-similar states in homogeneous rotating turbulence.

Squires, Kyle D.↗

Computations Of Breakup Of A Capillary Jet

Report describes computations of breakup of round jet of liquid in gas or vacuum. Breakup occurs both spontaneously and when stimulated by vibrations giving rise to surface waves, some of which grow until drops pinch off. Study focusses on formation of satellite drops, taking advantage of advanced computers and computational methods to carry theoretical analysis beyond small-perturbation, linear regime.

Mansour, Nagi N.↗

Direct numerical simulations of a spatially developing plane wake

In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Maekawa, Hiroshi↗

Instability mode interactions in a spatially developing plane wake

Direct numerical simulations are used to study the development of various instability modes in a spatially developing 2D wake. Five types of forcing of the inlet are investigated: fundamental mode, fundamental and one or two subharmonics, fundamental mode and random noise, and random noise only. Statistical analyses are carried out, and some numerical results are compared with experimental measurements.

Maekawa, Hiroshi↗

The effects of rotation on initially anisotropic homogeneous flows

Rotation of initially anisotropic homogeneous flows is studied using a model spectral tensor. It is shown that the anisotropy changes because of the influence of rotation through phase scrambling. Phase scrambling causes the Reynolds stresses to develop with damped oscillations. The final Reynolds stress anisotropy is found to be proportional to the initial structural tensor anisotropy. Closure models for the rapid pressure strain terms should reflect this change in anisotropy, and should drive the anisotropy to reach its final predicted state. Finally, it is shown that long-time integration using direct numerical simulations should be treated with care because phase scrambling effects on a discrete wave space can cause loss of resolution when time becomes large.

Mansour, Nagi N.↗

Modeling the rapid spherical compression of isotropic turbulence

Rapid distortion analysis is used to modify the form of the closure model for the dissipation rate of the turbulent kinetic energy. The modification is such that the evolution of the dissipation rate during a rapid compression is predicted exactly; good agreement between the model prediction and direct simulation data is obtained. Previous closure proposals fail to properly predict the rapid compression case. The reason for the difference between the present and previous models is traced to the fact that previous workers neglected variations of kinematic viscosity.

Coleman, Gary N.↗