粘性可压混合层时间稳定性对称紧致差分求解

基于可压扰动方程组的一阶改型 ,将高精度对称紧致格式引入边值法数值线性稳定性分析。对所获非线性离散特征值问题给出了一个通用形式二阶迭代局部算法 ,实现了时间模式和空间模式的统一求解 ,并将扰动特征值及其特征函数同时得到。据此分析了可压平面自由混合层时间稳定性 ,涉及二维 /三维扰动波、粘性 /无粘扰动波、第一 /第二模态、特征函数、伪特征值谱等。研究表明 ,压缩性效应和粘性效应对最不稳定扰动波数和增长率呈相似的减抑作用 ;在 Mc=1附近 ,从高波数段开始 ,粘性效应可强化二维不稳定扰动波由第一模态向第二模态的过渡     

Efficient computation of the eigenspectrum of viscoelastic flows using submatrix-based transformation of the linearized equations

Linear stability analysis of incompressible viscoelastic flows based on normal mode expansions of the eigenfunctions requires the numerical solution of a generalized eigenvalue problem (GEVP). The complex boundary layer structure of the leading eigenfunctions and the singular character of the continuous set of eigenvalues, necessitate the use of fine mesh sizes, leading to large algebraic GEVPs. In this paper, we present a submatrix-based transformation of the linearized equations (SubTLE) that converts the GEVP into a simple eigenvalue problem (EVP) of half the original dimension for the purely elastic isothermal and non-isothermal flows of an Oldroyd-B liquid. This leads to significant (up to an order of magnitude) reduction in the CPU time and memory required for the solution of the EVP. This is illustrated in the context of isothermal and non-isothermal shear flows.     

Stability Analysis of Diffusive Displacement in Three-Layer Hele-Shaw Cell or Porous Medium

We study the linear stability of three-layer Hele-Shaw flow, which models the secondary oil recovery by polymer flooding, in the presence of a diffusion process and a variable viscosity in the middle layer (denoted by M.L.). Then the hydrodynamic stability of the flow is related with the advection–diffusion equation of the species. The diffusion coefficient and the viscosity in M.L. are used as parameters for minimizing the Saffman–Taylor instability. This model was studied also by Daripa and Pa?a (Transp Porous Med 70(1):11–23, 2007). A particular basic solution was considered. The stabilizing effect of diffusion was proved, by using a variational formulation of the stability system. However, this analytical method was not giving sufficient conditions for improving the stability; the obtained upper bound of the growth constant (in time) of the perturbations was depending on the eigenfunctions of the stability system. In this paper, we improve the above result. We use a discretization method and obtain a classical algebraic eigenvalue problem, equivalent with the Sturm-Liouville system which governs the flow stability. A generalization of the Gerschgorin’s localization theorem is given and two estimates of the growth constant are obtained, not depending on the eigenfunctions. The new estimates are used to obtain sufficient conditions for improving the stability. These conditions are given in terms of the viscosity profile, the diffusion coefficient, the injection velocity, and the M.L. length. We conclude that a strong diffusion process improves the stability in the range of large wavenumbers. In the range of small wavenumbers, a stability improvement is obtained if the viscosity jump on the M.L.–oil interface is small enough and the length of M.L. is large enough.     

Fractional optimal control of a 2-dimensional distributed system using eigenfunctions

This paper presents an eigenfunctions expansion based scheme for Fractional Optimal Control (FOC) of a 2-dimensional distributed system. The fractional derivative is defined in the Riemann–Liouville sense. The performance index of a FOC problem is considered as a function of both state and control variables, and the dynamic constraints are expressed by a Partial Fractional Differential Equation (PFDE) containing two space parameters and one time parameter. Eigenfunctions are used to eliminate the terms containing space parameters and to define the problem in terms of a set of generalized state and control variables. For numerical computation Grünwald–Letnikov approximation is used. A direct numerical technique is proposed to obtain the state and the control variables. For a linear case, the numerical technique results into a set of algebraic equations which can be solved using a direct or an iterative scheme. The problem is solved for different number of eigenfunctions and time discretization. Numerical results show that only a few eigenfunctions are sufficient to obtain good results, and the solutions converge as the size of the time step is reduced.     

Propagation of a Tollmien-Schlichting wave over the junction between rigid and compliant surfaces

The propagation of an instability wave over the junction region between rigid and compliant panels is studied theoretically. The problem is investigated using three different methods with reference to flow in a plane channel containing sections with elastic walls. Within the framework of the first approach, using the solution of the problem of flow receptivity to local wall vibration, the problem considered is reduced to the solution of an integro-differential equation for the complex wall oscillation amplitude. It is shown that at the junction of rigid and elastic channel walls the instability-wave amplitude changes stepwise. For calculating the step value, another, analytical, method of investigating the perturbation propagation process, based on representing the solution as a superposition of modes of the locally homogeneous problem, is proposed. This method is also applied to calculating the flow stability characteristics in channels containing one or more elastic sections or consisting of periodically alternating rigid and compliant sections. The third method represents the unknown solution as the sum of a local forced solution and a superposition of orthogonal modes of flow in a channel with rigid walls. The latter method can be used for calculating the eigenvalues and eigenfunctions of the stability problem for flow in a channel with uniformly compliant walls.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 31–48. Original Russian Text Copyright © 2004 by Manuilovich.     

Consideration of longitudinal diffusion in flow within a channel

The stationary problem of convective diffusion in a channel with absorbent walls is considered. It is assumed that a Poiseuille flow exists. Two methods are employed in the solution, the method of separation of variables, and the method of expansion in eigenfunctions of the corresponding problem with piston profile (expansion method). It is established by comparison with independently obtained solutions for high Peclet number that for the first eigenfunctions and eigenvalues the expansion method gives satisfactory results over the entire Peclet-number range. For approximate calculation of subsequent eigenfunctions and eigenvalues a modification of the smooth asymptotic expansion method is used. The results are used to calculate matter flow density on the wall, to evaluate the length of the entrance region, and to obtain an analytical expression for the limiting Nusselt number in terms of the Peclet number.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 63–73, November–December, 1973.     

Solution to the orr–sommerfeld equation for liquid film flowing down an inclined plane: An optimal approach

The Orr–Sommerfeld equation is solved numerically for a layer of liquid film flowing down an inclined plane under the action of gravity using the sequential gradient-restoration algorithm (SGRA) The method consists of solving the governing equation as it is a Bolza problem in the calculus of variations. The neutral stability curves, eigenvalues and eigenfunctions to the stability problem can be determined simultaneously during the process.     

Torsion of an elastic rod whose cross-section is a parallelogram,trapezoid, or triangle,or has an arbitrary shape by the method of transformation to a rectangular domain

A coordinate transformation is used to take the domain of the rod cross-section to a rectangular domain for which the spectra of eigenfunctions and eigenvalues are known. The torsion function is represented as a generalized Fourier series to reduce the problem to solving a closed linear system of algebraic equations for the expansion coefficients. It is shown that these Fourier series converge absolutely, because the expansion coefficients decrease by a cubic law depending on the term number. We prove that the approximate solution in the form of a finite sum of the Fourier series converges to the exact solution. This theorem is generalized to the case of a rod cross-section of arbitrary shape.     

A spectral approach for the stability analysis of turbulent open-channel flows over granular beds

A novel Orr–Sommerfeld-like equation for gravity-driven turbulent open-channel flows over a granular erodible bed is here derived, and the linear stability analysis is developed. The whole spectrum of eigenvalues and eigenvectors of the complete generalized eigenvalue problem is computed and analyzed. The fourth-order eigenvalue problem presents singular non-polynomial coefficients with non-homogenous Robin-type boundary conditions that involve first and second derivatives. Furthermore, the Exner condition is imposed at an internal point. We propose a numerical discretization of spectral type based on a single-domain Galerkin scheme. In order to manage the presence of singular coefficients, some properties of Jacobi polynomials have been carefully blended with numerical integration of Gauss–Legendre type. The results show a positive agreement with the classical experimental data and allow one to relate the different types of instability to such parameters as the Froude number, wavenumber, and the roughness scale. The eigenfunctions allow two types of boundary layers to be distinguished, scaling, respectively, with the roughness height and the saltation layer for the bedload sediment transport.     

On Nonlinear Instability and Stability for Stratified Shear Flow

An example of stratified shear flow is presented in which an explicit construction is given for unstable eigenvalues with smooth eigenfunctions for the Taylor--Goldstein equation. It is proved for any stratified, plane parallel shear flow that the unstable spectrum of the linear operator is purely discrete. A general theorem is then invoked to prove that the specific example is nonlinearly unstable. A sufficient condition for nonlinear stability for stratified shear flow is discussed.     

On the stability of multi-layered fluids. Part I: analysis

The linear stability of the Poiseuille flow of multi-layered different fluids, described mathematically by a system of Orr-Somerfeld differential equations, is investigated. A spectral method is used to rewrite this system into a generalized eigenvalue problem, which can be solved with the QZ-algorithm. Special attention is paid to the tractibility of the interfacial conditions of the stability problem. Since we will limit ourselves to a linear stability analysis, the analytical treatment of the interfacial conditions is simplified. Some results related to simple flow configurations are presented. The origin of certain regions of interfacial instability is explained by simple analytical reasoning.     

Laminar heat transfer in parallel plate flow

An exact solution for the fluid temperature due to laminar heat transfer in parallel plate flow is found. The formulas obtained are valid for an arbitrary velocity profile. The basic problem encountered involves finding certain expansion coefficients in a series of nonorthogonal eigenfunctions. This problem is solved by passing to a vector system of equations having orthogonal eigenvectors. The method is applicable to more general problems.     

Eigenoscillations of a fluid in a canonical domain and functional difference equations

In this work we construct and discuss special solutions of a homogeneous problem for the Laplace equation in a domain with cone-shaped boundaries. The problem at hand is interpreted as that describing oscillatory linear wave movement of a fluid under gravity in such a domain. These solutions are found in terms of the Mellin transform and by means of the reduction to some new functional-difference equations solved in an explicit form (by quadrature). The behavior of the solutions at large distances is studied by use of the saddle point technique. The corresponding eigenoscillations of a fluid are then interpreted as generalized eigenfunctions of the continuous spectrum.     

Analysis of superposed fluids by the finite element method: Linear stability and flow development

A Galerkin finite element method is described for studying the stability of two superposed immiscible Newtonian fluids in plane Poiseuille flow. The formulation results in an algebraic eigenvalue problem of the form Aλ² + Bλ + C = 0 which, after transforming to a standard generalized eigenvalue problem, is solved by the QR algorithm. The numerical results are in good agreement with previous asymptotic results. Additional results show that the finite element method is ideally suited for studying linear stability of superposed fluids when parameters characterizing the flow fall outside the range amenable to perturbation methods. The applicability of the finite element method to similar eigenvalue problems is demonstrated by analysing the steady-state spatial development of two superposed fluids in a channel.     

Some Flows in Shape Optimization

Geometric flows related to shape optimization problems of the Bernoulli type are investigated. The evolution law is the sum of a curvature term and a nonlocal term of Hele–Shaw type. We introduce generalized set solutions, the definition of which is widely inspired by viscosity solutions. The main result is an inclusion preservation principle for generalized solutions. As a consequence, we obtain existence, uniqueness and stability of solutions. Asymptotic behavior for the flow is discussed:we prove that the solutions converge to a generalized Bernoulli exterior free-boundary problem.     

Calculation of neutral monotonic instability curves in the problem of concentration convection in a mixture with an infinite number of components

The Oberbeck-Boussinesq approximation for concentration convection in a mixture with an infinite number of components is constructed. The features of the formulation of the problem are described in detail. The large-parameter asymptotics are constructed for the linear problem of hydrodynamic stability. The problem is substantially simplified and equations not previously encountered in hydrodynamic stability theory are obtained. In the case of the non self-adjoint problem the asymptotics of the eigenvalues and eigenfunctions are obtained. Numerical results which, in particular, show that the spectrum of the boundary value problem is not connected are presented. The critical values obtained make it possible to solve the important practical problem of improving the process of mixture separation by the isoelectric focusing method.Rostov-on-Don. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 11–20, September–October, 1995.     

Creeping flow analyses of free surface cavity flows

Two industrially important free surface flows arising in polymer processing and thin film coating applications are modelled as lid-driven cavity problems to which a creeping flow analysis is applied. Each is formulated as a biharmonic boundary-value problem and solved both analytically and numerically. The analytical solutions take the form of a truncated biharmonic series of eigenfunctions for the streamfunction, while numerical results are obtained using a linear, finite-element formulation of the governing equations written in terms of both the streamfunction and vorticity. A key feature of the latter is that problems associated with singularities are alleviated by expanding the solution there in a series of separated eigenfunctions. Both sets of results are found to be in extremely good agreement and reveal distinctive flow transformations that occur as the operating parameters are varied. They also compare well with other published work and experimental observation.     

Dynamical eigenfunction decomposition of turbulent channel flow

The results of an analysis of low-Reynolds-number turbulent channel flow based on the Karhunen-Loéve(K-L) expansion are presented. The turbulent flow field is generated by a direct numerical simulation of the Navier-Stokes equations at a Reynolds number Re,= 80 (based on the wall shear velocity and channel half-width). The K-L procedure is then applied to determine the eigenvalues and eigenfunctions for this flow. The random coefficients of the K-L expansion are subsequently found by projecting the numerical flow field onto these eigenfunctions. The resulting expansion captures 90% of the turbulent energy with significantly fewer modes than the original trigonometric expansion. The eigenfunctions, which appear either as rolls or shearing motions, posses viscous boundary layers at the walls and are much richer in harmonics than the original basis functions. Chaotic temporal behaviour is observed in all modes and increases for higher-order eigenfunctions. The structure and dynamical behaviour of the eigenmodes are discussed as well as their use in the representation of the turbulent flow.     

Buckling of a flush-mounted plate in simple shear flow

The buckling of an elastic plate with arbitrary shape flush-mounted on a rigid wall and deforming under the action of a uniform tangential load due to an overpassing simple shear flow is considered. Working under the auspices of the theory of elastic instability of plates governed by the linear von Kármán equation, an eigenvalue problem is formulated for the buckled state resulting in a fourth-order partial differential equation with position-dependent coefficients parameterized by the Poisson ratio. The governing equation also describes the deformation of a plate clamped around the edges on a vertical wall and buckling under the action of its own weight. Solutions are computed analytically for a circular plate by applying a Fourier series expansion to derive an infinite system of coupled ordinary differential equations and then implementing orthogonal collocation, and numerically for elliptical and rectangular plates by using a finite-element method. The eigenvalues of the resulting generalized algebraic eigenvalue problem are bifurcation points in the solution space, physically representing critical thresholds of the uniform tangential load above which the plate buckles and wrinkles due to the partially compressive developing stresses. The associated eigenfunctions representing possible modes of deformation are illustrated, and the effect of the Poisson ratio and plate shape is discussed.     

Stability of Hartmann flows

The stability of Hartmann flows for arbitrary magnetic Reynolds numbers is investigated in the framework of linear theory. The initial three-dimensional problem reduces to the equivalent two-dimensional problem. Perturbation theory is used to find asymptotic expressions for the eigenvalues. Distinguishing two types of disturbances — magnetic and hydrodynamic — is shown to be advantageous in a number of cases. Simple features of the stability are considered for particular cases. The well-know Lundquist result is generalized. An energy approach is applied to the problem of stability. The results of simulations involving the solution of the linear stability problem are described. A distinctive picture of stability is developed. There are several types of instability and they can develop simultaneously. The hydrodynamic and magnetic phenomena interact with each other in a very complex fashion. The magnetic field can either enhance flow stability or reduce it.Novosibirsk. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 6, pp. 17–31, November–December, 1972.