59 research outputs found
Multiple paths to subharmonic laminar breakdown in a boundary layer
Numerical simulations demonstrate that laminar breakdown in a boundary layer induced by the secondary instability of two-dimensional Tollmien-Schlichting waves to three-dimensional subharmonic disturbances need not take the conventional lambda vortex/high-shear layer path
On the linear stability of compressible plane Couette flow
The linear stability of compressible plane Couette flow is investigated. The correct and proper basic velocity and temperature distributions are perturbed by a small amplitude normal mode disturbance. The full small amplitude disturbance equations are solved numerically at finite Reynolds numbers, and the inviscid limit of these equations is then investigated in some detail. It is found that instability can occur, although the stability characteristics of the flow are quite different from unbounded flows. The effects of viscosity are also calculated, asymptotically, and shown to have a stabilizing role in all the cases investigated. Exceptional regimes to the problem occur when the wavespeed of the disturbances approaches the velocity of either of the walls, and these regimes are also analyzed in some detail. Finally, the effect of imposing radiation-type boundary conditions on the upper (moving) wall (in place of impermeability) is investigated, and shown to yield results common to both bounded and unbounded flows
A spectral collocation solution to the compressible stability eigenvalue problem
A newly developed spectral compressible linear stability code (SPECLS) (staggered pressure mesh) is presented for analysis of shear flow stability, and applied to high speed boundary layers and free shear flows. The formulation utilizes the first application of a staggered mesh for a compressible flow analysis by a spectral technique. An order of magnitude less number of points is needed for equivalent accuracy of growth rates compared to those calculated by a finite difference formulation. Supersonic disturbances which are found to have oscillatory structures were resolved by a spectral multi-domain discretization, which requires a factor of three fewer points than the single domain spectral stability code. It is indicated, as expected, that stability of mixing layers is enhanced by viscosity and increasing Mach number. The mean flow involves a jet being injected into a quiescent gas. Higher temperatures of the injected gas is also found to enhance stability characteristics of the free shear layer
Numerical simulation of a controlled boundary layer
The problem of interest is the boundary layer over a flat plate. The three standard laminar flow control (LFC) techniques are pressure gradient, suction, and heating. The parameters used to describe the amount of control in the context of the boundary layer equations are introduced. The numerical method required to find the mean flow, the linear eigenvalues of the Orr-Sommerfeld equation, and the full, nonlinear, 3-D solution of the Navier-Stokes equations are outlined. A secondary instability exists for the parallel boundary subject to uniform pressure gradient, suction, or heating. Selective control of the spanwise mode reduces the secondary instability in the parallel boundary layer at low Reynolds number
Spectral methods for CFD
One of the objectives of these notes is to provide a basic introduction to spectral methods with a particular emphasis on applications to computational fluid dynamics. Another objective is to summarize some of the most important developments in spectral methods in the last two years. The fundamentals of spectral methods for simple problems will be covered in depth, and the essential elements of several fluid dynamical applications will be sketched
Active control of instabilities in laminar boundary-layer flow. Part 1: An overview
This paper (the first in a series) focuses on using active-control methods to maintain laminar flow in a region of the flow in which the natural instabilities, if left unattended, lead to turbulent flow. The authors review previous studies that examine wave cancellation (currently the most prominent method) and solve the unsteady, nonlinear Navier-Stokes equations to evaluate this method of controlling instabilities. It is definitely shown that instabilities are controlled by the linear summation of waves (i.e., wave cancellation). Although a mathematically complete method for controlling arbitrary instabilities has been developed (but not yet tested), the review, duplication, and physical explanation of previous studies are important steps for providing an independent verification of those studies, for establishing a framework for subsequent work which will involve automated transition control, and for detailing the phenomena by which the automated studies can be used to expand knowledge of flow control
Towards Perfectly Absorbing Boundary Conditions for Euler Equations
In this paper, we examine the effectiveness of absorbing layers as non-reflecting computational boundaries for the Euler equations. The absorbing-layer equations are simply obtained by splitting the governing equations in the coordinate directions and introducing absorption coefficients in each split equation. This methodology is similar to that used by Berenger for the numerical solutions of Maxwell's equations. Specifically, we apply this methodology to three physical problems shock-vortex interactions, a plane free shear flow and an axisymmetric jet- with emphasis on acoustic wave propagation. Our numerical results indicate that the use of absorbing layers effectively minimizes numerical reflection in all three problems considered
Fast orthogonal derivatives on the star
AbstractIn many numerical problems there is the need for obtaining derivatives in the X and Y directions of m variables at each point on an n×n plane. We consider the case where these derivatives are obtained using spectral methods (i.e. n fast Fourier transforms of length n are taken for each component, multiplied by the wave numbers and reverse transformed).On the CDC STAR-100 all data points corresponding to a plane must be stored in contiguous locations if advantage is to be taken of the powerful pipeline hardware of the machine. This means that derivatives in one direction are obtained very efficiently while derivatives in the orthogonal direction require either the substantial overhead of transposition or the use of scalar operations with no benefits of pipelining.An algorithm is described that overcomes this problem by taking derivatives of all components simultaneously. This is made possible by perfect shuffling of data to effect a pseudo-transposition that permits the FFT routine to take transforms of all m components on a plane at one time. Practical experience with this algorithm for m=5 and n=32 shows a 10% speedup for X-derivatives and a 32% speedup for Y-derivatives over the conventional algorithms (in which X and Y derivatives are taken one component at a time and Y derivatives require transposition of data).A theoretical analysis based on available STAR-100 vector instruction timing data predicts that this algorithm is superior to the conventional algorithm for M ≥ 2, n ≤ 128 (problem sizes of practical interest). We show how further improvement in running time may be obtained if derivatives of several components on more than one plane are required.This analysis is applicable to the new generation of STAR computers (the CDC Cyber 203s) since vector instruction timings are essentially unchanged in the new machines
Optimal unified combination rule in application of Dempster-Shafer theory to lung cancer radiotherapy dose response outcome analysis
- …