Linear and Nonlinear PSE for Stability Analysis of the Blasius Boundary Layer Using Compact Scheme

2001 ◽  
Vol 123 (3) ◽  
pp. 545-550 ◽  
Author(s):  
V. Esfahanian ◽  
K. Hejranfar ◽  
F. Sabetghadam

A highly accurate finite-difference PSE code has been developed to investigate the stability analysis of incompressible boundary layers over a flat plate. The PSE equations are derived in terms of primitive variables and are solved numerically by using compact method. In these formulations, both nonparallel as well as nonlinear effects are accounted for. The validity of present numerical scheme is demonstrated using spatial simulations of two cases; two-dimensional (linear and nonlinear) Tollmien-Schlichting wave propagation and three-dimensional subharmonic instability breakdown. The PSE solutions have been compared with previous numerical investigations and experimental results and show good agreement.

2001 ◽  
Vol 432 ◽  
pp. 69-90 ◽  
Author(s):  
RUDOLPH A. KING ◽  
KENNETH S. BREUER

An experimental investigation was conducted to examine acoustic receptivity and subsequent boundary-layer instability evolution for a Blasius boundary layer formed on a flat plate in the presence of two-dimensional and oblique (three-dimensional) surface waviness. The effect of the non-localized surface roughness geometry and acoustic wave amplitude on the receptivity process was explored. The surface roughness had a well-defined wavenumber spectrum with fundamental wavenumber kw. A planar downstream-travelling acoustic wave was created to temporally excite the flow near the resonance frequency of an unstable eigenmode corresponding to kts = kw. The range of acoustic forcing levels, ε, and roughness heights, Δh, examined resulted in a linear dependence of receptivity coefficients; however, the larger values of the forcing combination εΔh resulted in subsequent nonlinear development of the Tollmien–Schlichting (T–S) wave. This study provides the first experimental evidence of a marked increase in the receptivity coefficient with increasing obliqueness of the surface waviness in excellent agreement with theory. Detuning of the two-dimensional and oblique disturbances was investigated by varying the streamwise wall-roughness wavenumber αw and measuring the T–S response. For the configuration where laminar-to-turbulent breakdown occurred, the breakdown process was found to be dominated by energy at the fundamental and harmonic frequencies, indicative of K-type breakdown.


2006 ◽  
Vol 17 (01) ◽  
pp. 65-73 ◽  
Author(s):  
SHIRO SAWADA

The optimal velocity model which depends not only on the headway but also on the relative velocity is analyzed in detail. We investigate the effect of considering the relative velocity based on the linear and nonlinear analysis of the model. The linear stability analysis shows that the improvement in the stability of the traffic flow is obtained by taking into account the relative velocity. From the nonlinear analysis, the relative velocity dependence of the propagating kink solution for traffic jam is obtained. The relation between the headway and the velocity and the fundamental diagram are examined by numerical simulation. We find that the results by the linear and nonlinear analysis of the model are in good agreement with the numerical results.


2011 ◽  
Vol 137 ◽  
pp. 72-76
Author(s):  
Wei Zhang ◽  
Xian Wen ◽  
Yan Qun Jiang

A proper orthogonal decomposition (POD) method is applied to study the global stability analysis for flow past a stationary circular cylinder. The flow database at Re=100 is obtained by CFD software, i.e. FLUENT, with which POD bases are constructed by a snapshot method. Based on the POD bases, a low-dimensional model is established for solving the two-dimensional incompressible NS equations. The stability of the flow solution is evaluated by a POD-Chiba method in the way of the eigensystem analysis for the velocity disturbance. The linear stability analysis shows that the first Hopf bifurcation takes place at Re=46.9, which is in good agreement with available results by other high-order accurate stability analysis methods. However, the calculated amount of POD is little, which shows the availability and advantage of the POD method.


2008 ◽  
Vol 3 (3) ◽  
pp. 34-38
Author(s):  
Sergey A. Gaponov ◽  
Yuri G. Yermolaev ◽  
Aleksandr D. Kosinov ◽  
Nikolay V. Semionov ◽  
Boris V. Smorodsky

Theoretical and an experimental research results of the disturbances development in a swept wing boundary layer are presented at Mach number М = 2. In experiments development of natural and small amplitude controllable disturbances downstream was studied. Experiments were carried out on a swept wing model with a lenticular profile at a zero attack angle. The swept angle of a leading edge was 40°. Wave parameters of moving disturbances were determined. In frames of the linear theory and an approach of the local self-similar mean flow the stability of a compressible three-dimensional boundary layer is studied. Good agreement of the theory with experimental results for transversal scales of unstable vertices of the secondary flow was obtained. However the calculated amplification rates differ from measured values considerably. This disagreement is explained by the nonlinear processes observed in experiment


2018 ◽  
Vol 841 ◽  
pp. 636-653
Author(s):  
Ting-Yueh Chang ◽  
Falin Chen ◽  
Min-Hsing Chang

A three-dimensional linear stability analysis is carried out for a convecting layer in which both the temperature and solute distributions are linear in the horizontal direction. The three-dimensional results show that, for $Le=3$ and 100, the most unstable mode occurs invariably as the longitudinal mode, a vortex roll with its axis perpendicular to the longitudinal plane, suggesting that the two-dimensional results are sufficient to illustrate the stability characteristics of the convecting layer. Two-dimensional results show that the stability boundaries of the transverse mode (a vortex roll with its axis perpendicular to the transverse plane) and the longitudinal modes are virtually overlapped in the regime dominated by thermal diffusion and the regime dominated by solute diffusion, while these two modes hold a significant difference in the regime the salt-finger instability prevails. More precisely, the instability area in terms of thermal Grashof number $Gr$ and solute Grashof number $Gs$ is larger for the longitudinal mode than the transverse mode, implying that, under any circumstance, the longitudinal mode is always more unstable than the transverse mode.


Spikes and spots are discussed mostly for incompressible boundary layers, with the emphasis towards strong nonlinearity. The distinction between forced and free disturbances then becomes blurred, as spikes and spots reproduce each other. First, the forced case is concentrated on the start of spikes. The theory used is that of the two- or three-dimensional interacting boundary layer, capturing nonlinear Tollmien-Schlichting waves, for example, or following a vortex-wave interaction. Finite-time breakup produces shortened time and length scales, yielding agreement with computations and experiments on the first spike in transition, with subsequent spot formation. After the breakup, normal pressure gradients and vortex wind-up become significant locally. Second, the free case concerns initial-value problems for spots containing a wide band of three-dimensional nonlinear disturbances. The theory points to successive nonlinear stages starting at the wing tips near the spot trailing edge but gradually entering the middle as the amplitudes increase downstream. This effect combined with shortening scales produces a spread angle near 11°, very close to the experimental observations. The overall spot structure is described briefly, including also the leading edge. Viscosity arises later in two ways; for the case mentioned above with spikes originating near the surface and also through a novel interaction influencing the global spot.


2002 ◽  
Author(s):  
Leslie Ng ◽  
Richard Rand

We investigate the effect of nonlinearites on a parametrically excited ordinary differential equation whose linearization exhibits the phenomena of coexistence. The differential equation studied governs the stability mode of vibration in an unforced conservative two degree of freedom system used to model the free vibrations of a thin elastica. Using perturbation methods, we show that at parameter values corresponding to coexistence, nonlinear terms can cause the origin to become nonlinearly unstable, even though linear stability analysis predicts the origin to be stable. We also investigate the bifurcations associated with this instability.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


Sign in / Sign up

Export Citation Format

Share Document