Applications of parabolized stability equation for predicting transition position in boundary layers

The phenomenon of laminar-turbulent transition exists universally in nature and various engineering practice.The prediction of transition position is one of crucial theories and practical problems in fluid mechanics due to the different characteristics of laminar flow and turbulent flow.Two types of disturbances are imposed at the entrance,i.e.,identical amplitude and wavepacket disturbances,along the spanwise direction in the incompressible boundary layers.The disturbances of identical amplitude are consisted of one two-dimensional(2D) wave and two three-dimensional(3D) waves.The parabolized stability equation(PSE) is used to research the evolution of disturbances and to predict the transition position.The results are compared with those obtained by the numerical simulation.The results show that the PSE method can investigate the evolution of disturbances and predict the transition position.At the same time,the calculation speed is much faster than that of the numerical simulation.     

On the wave-cancelling nature of boundary layer flow control

This work deals with the feedforward active control of Tollmien–Schlichting instability waves over incompressible 2D and 3D boundary layers. Through an extensive numerical study, two strategies are evaluated; the optimal linear–quadratic–Gaussian (LQG) controller, designed using the Eigensystem realization algorithm, is compared to a wave-cancellation scheme, which is obtained using the direct inversion of frequency-domain transfer functions of the system. For the evaluated cases, it is shown that LQG leads to a similar control law and presents a comparable performance to the simpler, wave-cancellation scheme, indicating that the former acts via a destructive interference of the incoming wavepacket downstream of actuation. The results allow further insight into the physics behind flow control of convectively unstable flows permitting, for instance, the optimization of the transverse position for actuation. Using concepts of linear stability theory and the derived transfer function, a more efficient actuation for flow control is chosen, leading to similar attenuation of Tollmien–Schlichting waves with only about 10% of the actuation power in the baseline case.     

Spatial direct numerical simulation of high-speed boundary-layer flows Part II: Transition on a cone in Mach 8 flow

The laminar breakdown of the boundary-layer flow of an axisymmetric sharp cone in a Mach 8 flow is simulated by a synergistic approach that combines the parabolized stability equation (PSE) method and spatial direct numerical simulation (DNS). The transitional state is triggered by a symmetric pair of oblique second-mode disturbances whose nonlinear interactions generate strong streamwise vorticity, which leads in turn to severe spanwise variations in the flow and eventual laminar breakdown. The PSE method is used to compute the weakly and moderately nonlinear initial stages of the transition process and, thereby, to derive a harmonically rich inflow condition for the DNS. The strongly nonlinear and laminar-breakdown stages of transition are then computed by well-resolved DNS, with a highly accurate algorithm that exploits spectral collocation and high-order compact-difference methods. Evolution of the flow is presented in terms of modal energies, mean quantities (e.g., skin friction), Reynolds stresses, turbulent kinetic energy, and flow visualization. The numerical test case is an approximate computational analog of one of the few stability experiments performed for hypersonic boundary-layer flows. Comparisons and contrasts are drawn between the experimental and the computational results. Rope-like waves similar to those observed in schlieren images of high-speed transitional flows are also observed in the numerical experiment and are shown to be visual manifestations of second-mode instability waves.This research was supported under NASA Contracts NAS1-19831 and NAS1-20059 for the first and second authors, respectively.     

Long wavelength approximation to peristaltic motion of a power law fluid

Peristaltic motion of a power law fluid in a two-dimensional channel is studied. Assuming that the wavelength of the peristaltic wave is large in comparison to the mean half-width of the channel, a solution for the stream function is obtained as an asymptotic expansion in terms of slope parameter. Expressions for axial pressure gradient and shear stress are derived. The effect of flow behaviour indexn on the streamline pattern and shear stress is studied and the phenomenon of trapping is discussed.     

Nonparallel stability of the flow in an axially rotating pipe

The linear stability of the developing flow in an axially rotating pipe is analyzed using parabolized stability equations (PSE). The results are compared with those obtained from a near-parallel stability approximation that only takes into account the axial variation of the basic flow. Though the PSE results obviously coincide with the near-parallel ones far downstream, when the flow has reached a Hagen-Poiseuille axial velocity profile with superimposed solid-body rotation, they differ significantly in the developing region. Therefore, the onset of instability strongly depends on the axial evolution of the perturbations. The PSE results are also compared with experimental data from Imao et al. [Exp. Fluids 12 (1992) 277], showing a good agreement in the frequencies and wavelengths of the unstable disturbances, that take the form of spiral waves. Finally, a simple method for detecting one of the conditions to characterize the onset of absolute instability using PSE is given.     

The control of a permanent rotational motion and equilibrium position of a spacecraft

This article has adopted an analytical method to obtain a non-linear control law to reach the exponential asymptotic stablity of the permanent rotational motion of a spacecraft. The control moments achieving this rotational motion are obtained. The control moments to establish exponential asymptotic stablity of the mentioned motion are obtained as non-linear functions of the phase coordinates of the spacecraft. The general solution of the equations of perturbed motion is derived. Furthermore, analysis and numerical simulation study of this solution are presented. For numerical examples the time needed for control is calculated. An equilibrium position of the spacecraft is proved to be exponentially asymptotically stable as a special case of the above-studied problem.     

Vortex structure simulation for supersonic mixing layers using nonlinear PSE method

The method of nonlinear parabolized stability equations (PSE) is applied in the simulation of vortex structures in compressible mixing layer. The spatially-evolving unstable waves, which dominate the vortex structure, are investigated through spatial marching method. The instantaneous flow field is obtained by adding the harmonic waves to basic flow. The results show that T-S waves do not keep growing exponentially as the linear evolution, the energy transfer to high order harmonic modes, and that finally all harmonic modes get saturated due to nonlinear interaction. The mean flow distortion induced by the nonlinear interaction between the harmonic modes and their conjugate harmonic ones, makes great change of the average flow and increases the thickness of mixing layer. PSE methods can well capture the two- and three-dimensional large scale nonlinear vortex structures in mixing layers such as vortex roll-up, vortex pairing, and Λ vortex.     

Sensitivity Analysis Using Adjoint Parabolized Stability Equations for Compressible Flows

An input/output framework is used to analyze the sensitivity of two- and three-dimensional disturbances in a compressible boundary layer for changes in wall and momentum forcing. The sensitivity is defined as the gradient of the kinetic disturbance energy at a given downstream position with respect to the forcing. The gradients are derived using the parabolized stability equations (PSE) and their adjoint (APSE). The adjoint equations are derived in a consistent way for a quasi-two-dimensional compressible flow in an orthogonal curvilinear coordinate system. The input/output framework provides a basis for optimal control studies. Analysis of two-dimensional boundary layers for Mach numbers between 0 and 1.2 show that wall and momentum forcing close to branch I of the neutral stability curve give the maximum magnitude of the gradient. Forcing at the wall gives the largest magnitude using the wall normal velocity component. In case of incompressible flow, the two-dimensional disturbances are the most sensitive ones to wall inhomogeneity. For compressible flow, the three-dimensional disturbances are the most sensitive ones. Further, it is shown that momentum forcing is most effectively done in the vicinity of the critical layer. This revised version was published online in July 2006 with corrections to the Cover Date.     

An evolve‐then‐filter regularized reduced order model for convection‐dominated flows

In this paper, we propose a new evolve‐then‐filter reduced order model (EF‐ROM). This is a regularized ROM (Reg‐ROM), which aims to add numerical stabilization to proper orthogonal decomposition (POD) ROMs for convection‐dominated flows. We also consider the Leray ROM (L‐ROM). These two Reg‐ROMs use explicit ROM spatial filtering to smooth (regularize) various terms in the ROMs. Two spatial filters are used: a POD projection onto a POD subspace (Proj) and a POD differential filter (DF). The four Reg‐ROM/filter combinations are tested in the numerical simulation of the three‐dimensional flow past a circular cylinder at a Reynolds number Re=1000. Overall, the most accurate Reg‐ROM/filter combination is EF‐ROM‐DF. Furthermore, the spatial filter has a higher impact on the Reg‐ROM than the regularization used. Indeed, the DF generally yields better results than Proj for both the EF‐ROM and L‐ROM. Finally, the CPU times of the four Reg‐ROM/filter combinations are orders of magnitude lower than the CPU time of the DNS. Copyright © 2017 John Wiley & Sons, Ltd.     

Control structure interaction in the nonlinear analysis of flexible mechanical systems

The effect of the control structure interaction on the feedforward control law as well as the dynamics of flexible mechanical systems is examined in this investigation. An inverse dynamics procedure is developed for the analysis of the dynamic motion of interconnected rigid and flexible bodies. This method is used to examine the effect of the elastic deformation on the driving forces in flexible mechanical systems. The driving forces are expressed in terms of the specified motion trajectories and the deformations of the elastic members. The system equations of motion are formulated using Lagrange's equation. A finite element discretization of the flexible bodies is used to define the deformation degrees of freedom. The algebraic constraint equations that describe the motion trajectories and joint constraints between adjacent bodies are adjoined to the system differential equations of motion using the vector of Lagrange multipliers. A unique displacement field is then identified by imposing an appropriate set of reference conditions. The effect of the nonlinear centrifugal and Coriolis forces that depend on the body displacements and velocities are taken into consideration. A direct numerical integration method coupled with a Newton-Raphson algorithm is used to solve the resulting nonlinear differential and algebraic equations of motion. The formulation obtained for the flexible mechanical system is compared with the rigid body dynamic formulation. The effect of the sampling time, number of vibration modes, the viscous damping, and the selection of the constrained modes are examined. The results presented in this numerical study demonstrate that the use of the driving forees obtained using the rigid body analysis can lead to a significant error when these forces are used as the feedforward control law for the flexible mechanical system. The analysis presented in this investigation differs significantly from previously published work in many ways. It includes the effect of the structural flexibility on the centrifugal and Coriolis forces, it accounts for all inertia nonlinearities resulting from the coupling between the rigid body and elastic displacements, it uses a precise definition of the equipollent systems of forces in flexible body dynamics, it demonstrates the use of general purpose multibody computer codes in the feedforward control of flexible mechanical systems, and it demonstrates numerically the effect of the selected set of constrained modes on the feedforward control law.     

Spatial direct numerical simulation of boundary-layer transition mechanisms: Validation of PSE theory

A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to evaluate critically the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave breakdown. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the two-dimensional test problem. This discrepancy is attributed to far-field boundary-condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.     

Large Eddy Simulation (2D) of a Reacting Plane Mixing Layer Using Filtered Density Function Closure

The filtered density function (FDF) is implemented for a two-dimensional, large eddy simulation (LES) of a gas phase, spatially developing, reacting and non-reacting, constant-density, plane mixing layer in a flow regime prior to the mixing transition where the flow is mainly two-dimensional. The unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities following the FDF approach. In the derived FDF transport equation, the effect of chemical reactions appears in a closed form. The Lagrangian Monte Carlo scheme is used to solve the FDF transport equation. The applicability and performance of the FDF for LES of a reacting plane mixing layer are assessed by comparisons with experimental measurements. In non-reacting flow, the calculated mean streamwise velocity profiles and mean mixture fraction profiles relax to self-similarity, which is in satisfactory agreement with the measurements. In reacting flow, the FDF calculation provided a satisfactory accuracy in comparison with measurements of mean reactant and product concentration. The increase in the total amount of product formation in the flip case demonstrates the asymmetric characteristics of the entrainment and mixing characteristics in the mixing layer. This revised version was published online in July 2006 with corrections to the Cover Date.     

Equation of a shock wave front

Second-order differential equations of the hyperbolic type are derived for describing the local law of shock wave propagation. The shock waves are assumed to be two-dimensional unsteady in a stationary gas flow and three-dimensional steady in a supersonic flow. The behavior of the characteristics of these equations is investigated as a function of the governing flow parameters and their relative position with respect to the typical bicharacteristics of the characteristic cone behind the shock is analyzed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 159–165, May–June, 2000.     

Stability analysis method considering non-parallelism: EPSE method and its application

The e-N method is widely used in transition prediction. The amplitude growth rate used in the e-N method is usually provided by the linear stability theory (LST) based on the local parallel hypothesis. Considering the non-parallelism effect, the parabolized stability equation (PSE) method lacks local characteristic of stability analysis. In this paper, a local stability analysis method considering non-parallelism is proposed, termed as EPSE since it may be considered as an expansion of the PSE method. The EPSE considers variation of the shape function in the streamwise direction. Its local characteristic is convenient for stability analysis. This paper uses the EPSE in a strong non-parallel flow and mode exchange problem. The results agree well with the PSE and the direct numerical simulation (DNS). In addition, it is found that the growth rate is related to the normalized method in the non-parallel flow. Different results can be obtained using different normalized methods. Therefore, the normalized method must be consistent.     

PSE as applied to problems of transition in compressible boundary layers

A new idea of using the parabolized stability equation (PSE) method to predict laminar-turbulent transition is proposed. It is tested in the prediction of the location of transition for compressible boundary layers on flat plates, and the results are compared with those obtained by direct numerical simulations (DNS). The agreement is satisfactory, and the reason for this is that the PSE method faithfully reproduces the mechanism leading to the breakdown process in laminar-turbulent transition, i. e., the modification of mean flow profile leads to a remarkable change in its stability characteristics.     

弹性支撑圆柱绕流稳定性分析

基于CFD 技术,采用系统辨识方法,建立了亚临界雷诺数(Re < 47) 下绕圆柱流动的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程和降阶气动力模型,建立了弹性支撑圆柱绕流的稳定性分析模型. 算例分析了亚临界雷诺数下,结构固有频率、质量比等参数以及支撑方式对弹性系统稳定性的影响. 对于单自由度横向支撑圆柱,当结构固有频率趋近流动最不稳定模态频率时,弹性系统会在一定频率范围内失稳,这种现象最低可在Re~20 时出现. 旋转自由度的释放能够进一步降低系统的稳定性,可将临界雷诺数进一步降低至18 左右. ROM 方法不仅具有很高的效率,而且清晰地指出了弹性系统失稳的根本原因:流动模态和结构模态耦合作用导致结构模态失稳所致. 因此,失稳状态下系统振荡频率锁定于结构固有频率. 基于ROM 技术预测的失稳边界与直接CFD/CSD 仿真结果吻合,证明了该方法的正确性和精度.      

On the use of sensitivity analysis in model reduction to predict flows for varying inflow conditions

The proper orthogonal decomposition (POD)‐based model reduction method is more and more successfully used in fluid flows. However, the main drawback of this methodology rests in the robustness of these reduced order models (ROMs) beyond the reference at which POD modes have been derived. Any variation in the flow or shape parameters within the ROM fails to predict the correct dynamics of the flow field. To broaden the spectrum of these models, the POD modes should have the global characteristics of the flow field over which the predictions are required. Mixing of snapshots with varying parameters is one way to improve the global nature of the modes but is computationally demanding because it requires full‐order solutions for a number of parameter values in order to assemble atextitrich enough database on which to perform POD. Instead, we have used sensitivity analysis (SA) to include the flow and shape parameters influence during the basis selection process to develop more robust ROMs for varying viscosity (Reynolds number), changing orientation and shape definition of bodies. This study aims at extending these ideas to inflow conditions to demonstrate the effectiveness of the proposed approach in capturing the effect of varying inflow on the dynamics of the flow over an elliptic cylinder. Numerical experiments show that the newly derived models allow for a more accurate representation of the flows when exploring the parameter space. Copyright © 2011 John Wiley & Sons, Ltd.     

Investigation of the turbulent model for pressure fluctuations with spectral theory

The pressure fluctuations in turbulent shear flows are investigated with the theory of spectral analysis.An expression for pressure spectra is analytically derived in terms of velocity spectra.This derivation is based on a formal solution of the Navier-Stokes equation and quasi-normal assumption to express the third and fourth order velocity correlations in terms of double velocity correlation.Then,a turbulent model for the computation of pressure fluctuation intensity with Renolds stress and mean flow velocity gradients is established.The turbulent constants in the model are calculated from the assumptions about the general behaviour of velocity spectra in high Renolds number flows.Comparison with direct simulation of turbulent boundary layer is made.It is found that the turbulent-turbulent,cross correlation,and turbulent-shear source terms for mean square value of pressure fluctuation are about the same magnitude.     

Performance modeling and analysis of blood flow in elastic arteries

Nomenclature　　τ　　wallshearstressγshearrateτy yieldstressηc Cassonviscosityktheconsistencyindexnnon_Newtonianindexτp shearstressofthepthelementωangularvelocityRvessel’sradiusCwavespeedM　　magneticparameter (Hartmannnumber)u,w　velocitycomponentinther_andz_directions,respectivelyP　　pressureα　　unsteadinessparameter k , R meanparametersTp relaxationtimeofthepthelementρ densityIntroductionTheimportancetoatherogenesisofarterialflowphenomenasuchasflowseparation ,recirculationands…