Interaction of a three-dimensional boundary layer with an extensive obstacle

A special variant is considered of the theory of longitudinal—transverse interaction in which the pressure field in the perturbed region of flow forms under the influence of centrifugal forces which lead to a change in the pressure across the boundary layer. This regime of flow is realized in flow of an incompressible fluid, when the two-dimensional boundary layer developing along the smooth section of the contour of a solid body enters into interaction with a three-dimensional irregularity on the surface around which flow is taking place, a projection or a depression. On the assumption that the height of the irregularity is not great, a solution is constructed for the linearized problem of interaction. It is shown that the properties of the flow of fluid in the region of interaction, in particular the possibility of penetration of perturbations into the boundary layer in front of an irregularity, depend on the sign of the curvature exhibited by the contour of the body.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 39–48, January–February, 1988.     

On the Perturbation of Hamel Flow Due to Unevenness of the Channel Walls

A method of analyzing the receptivity of longitudinally inhomogeneous flows is proposed. The process of excitation of natural oscillations is studied with reference to the simplest inhomogeneous flow: the two-dimensional flow of a viscous incompressible fluid in a channel with plane nonparallel walls. As physical factors generating perturbations, the cases of a stationary irregularity and localized vibration of the channel walls are considered. By changing the independent variables and unknown functions of the perturbed flow, the problem of the generation of stationary perturbations above an irregularity is reduced to a longitudinally homogeneous boundary-value problem which is solved using a Fourier transform in the longitudinal variable. The same problem is investigated using another method based on representing the required solution in the form of a superposition of solutions of the homogeneous problem and a forced solution calculated in the locally homogeneous approximation. As a result, the problem of calculating the longitudinal distributions of the amplitudes of the normal modes is reduced to the solution of an infinite-dimensional inhomogeneous system of ordinary differential equations. The numerical solution obtained using this method is tested by comparison with an exact calculation based on the Fourier method. Using the method proposed, the problem of flow receptivity to harmonic oscillations of parts of the channel walls is analyzed. The calculations performed show that the method is promising for investigating the receptivity of longitudinally inhomogeneous flow in a laminar boundary layer.     

On local perturbations of a boundary layer with small skin friction

This paper studies the effect of two-dimensional surface irregularities on the flow in a plane steady boundary layer with small skin friction in an incompressible fluid. A detailed analysis is carried out for the flow regime with a given pressure gradient determined on the scale of the small irregularity by its shape. It is shown that there is a critical value of the height (depth) of the irregularity at which the skin friction first becomes zero, and the nonuniqueness of the corresponding solution is established. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 87–97, November–December, 1998. The work was financially supported by the Russian Foundation for Basic Research (project No. 97-01-00354).     

GREEN QUASIFUNCTION METHOD FOR FREE VIBRATION OF SIMPLY-SUPPORTED TRAPEZOIDAL SHALLOW SPHERICAL SHELL ON WINKLER FOUNDATION

The idea of Green quasifunction method is clarified in detail by considering a free vibration problem of simply-supported trapezoidal shallow spherical shell on Winkler foundation.A Green quasifunction is established by using the fundamental solution and boundary equation of the problem.This function satisfies the homogeneous boundary condition of the problem.The mode shape differential equation of the free vibration problem of simply-supported trapezoidal shallow spherical shell on Winkler foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula.There are multiple choices for the normalized boundary equation.Based on a chosen normalized boundary equation, the irregularity of the kernel of integral equations is avoided.Finally, natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations derived from the integral equations.Numerical results show high accuracy of the Green quasifunction method.     

Green quasifunction method for vibration of simply-supported thin polygonic plates on Pasternak foundation

A new numerical method—Green quasifunction is proposed.The idea of Green quasifunction method is clarified in detail by considering a vibration problem of simply-supported thin polygonic plates on Pasternak foundation.A Green quasifunction is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the problem.The mode shape differential equation of the vibration problem of simply-supported thin plates on Pasternak foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula.There are multiple choices for the normalized boundary equation.Based on a chosen normalized boundary equation,a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome.Finally,natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations derived from the integral equations.Numerical results show high accuracy of the Green quasifunction method.     

Quasi-Green’s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green’s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi-Green’s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob-lem. The mode shape differential equations of the free vibration problem of a simply-supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa-tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green’s function method.     

Quasi-Green＇s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green＇s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi- Green＇s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob- lem. The mode shape differential equations of the free vibration problem of a simply- supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa- tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green＇s function method.     

Quasi-Green’s function method for free vibration of clamped thin plates on Winkler foundation

The quasi-Green’s function method is used to solve the free vibration problem of clamped thin plates on the Winkler foundation. Quasi-Green’s function is established by the fundamental solution and the boundary equation of the problem. The function satisfies the homogeneous boundary condition of the problem. The mode-shape differential equation of the free vibration problem of clamped thin plates on the Winkler foundation is reduced to the Fredholm integral equation of the second kind by Green’s formula. The irregularity of the kernel of the integral equation is overcome by choosing a suitable form of the normalized boundary equation. The numerical results show the high accuracy of the proposed method.     

分形理论在煤层气非平衡吸附模型中的应用

传统的煤层气动力学模型均是建立在欧几里得几何基础上的,难以描述煤层孔隙结构的复杂性及形状的不规则性。本文以分形理论为基础,通过引入分形维数来刻画煤层孔隙结构的复杂性并考虑煤层的吸附特性、双重介质特征及介质的变形,建立基于Fick第二定律的分形介质煤层气非稳态渗流数学模型。由于流动方程的强非线性,结合各类边界条件用正则摄动法和Laplace变换得到模型在拉氏空间上的近似解析解,再利用Laplace数值反演求得实空间上的数值解。对参数进行敏感性分析并绘制了典型压力曲线,这些结果为煤层气开采提供了理论依据和试井方法。     

Unsteady flow diagnostics using weak perturbations

The examination of the motion of very weak waves generated from small sources on the boundary of a flow domain gives information on how the domain shape influences the flow, both from spatial and temporal perspectives. In the study of shock wave dynamics the shock itself generates a weak wave when passing over a small step, or surface irregularity. The basic principle is that if a particle produces a series of point disturbances in a flow field, the induced perturbations will propagate outwards at the local sonic velocity whilst at the same time being convected along with the local flow velocity. A number of issues may be identified for an unsteady flow. Firstly, the flow field at later times may be influenced by perturbations produced at earlier times. Secondly, if the positions of the perturbations can be monitored as a function of time, then the trajectory and velocity of the particle may be deduced. Thirdly, if a perturbation arises from a point on a boundary, then its influence, if any, on any particular part of the flow can be established. A number of examples are presented to illustrate the value of the technique and its potential to uncover the mechanisms responsible for the formation of certain flow patterns in high-speed compressible flows.     

Optical studies of local flow behaviour of a non-Newtonian fluid inside a porous medium

In this study the local flow of a polymeric solution inside a porous medium is studied visually. While the flow is quite uniform for low volume flow rates it shows pronounced nonuniformity for higher volume flow rates. That is to say, only certain preferred passages are taken and these passages change in time (they fluctuate). This flow irregularity is the reason for increased resistance. Received: 28 April 1997 Accepted: 30 December 1997     

Supersonic boundary layer transition induced by roughness on the attachment line of a yawed cylinder

The effect of a single two-dimensional irregularity and sandy roughness on boundary layer transition in supersonic flow over a yawed cylinder (M = 6)-has been experimentally investigated. The characteristic flow regimes beyond the roughness are identified, and their limits are determined as a function of the Reynolds number and the ratio of the height of the roughness to the characteristic thickness of the boundary layer. A qualitative comparison is made with the flow regimes induced by roughness on the attachment line in incompressible flow over a cylinder [1–3]. The thermal indicator coating method is used to measure the heat fluxes along the attachment line and a comparison is made with calculations carried out in accordance with the methods of other authors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 28–35, November–December, 1991.The authors are grateful to A. F. Kiselev for helping to calculate the heat flux in the turbulent boundary layer.     

低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

Flow irregularity and wear optimization in epitrochoidal gerotor pumps

The Gerotor pumps have a vast number of applications in industries and automobiles. The flow rate irregularity and wear rate proportional factor influence on the lifetime of the Gerotor pumps. In this paper, the optimization problem cost functions and constraints have been developed according to volumetric, dynamic and geometric properties. In order to have general optimum solution and reduce number of design variables, all variables have changed to non-dimensional variables by using the outer rotor lobe center radius which causes the non-dimensional cost functions too. The multi-objective optimization problem has been changed to single objective optimization problem by using a multi-objective optimization classical method. The single objective optimization problem has been solved by using of a mixed integer nonlinear optimization algorithm. The optimization programming has been implemented for several values of number of the outer rotor teeth, non-dimensional displacement and rotors width. On the other hand the influence of varied parameter such as number of outer rotors teeth, the value of non-dimensional displacement and non-dimensional rotor width have been surveyed on each of cost functions. The results show in constant input torque, as using greater number of outer rotor teeth up to critical value improves the wear of teeth. Also optimum value of flow rate irregularity in odd values of outer rotor teeth is smaller than even ones. At last for having better comparison, two sample commercial pumps have been optimized. The results show both wear rate proportional factor and flow irregularity have been significantly improved.     

A continuous shape sensitivity equation method for unsteady laminar flows

This paper presents the application of a general shape sensitivity equation method (SEM) to unsteady laminar flows. The formulation accounts for complex parameter dependence and is suitable for a wide range of problems. The flow and sensitivity equations are solved on 3D meshes using a Streamline-Upwind Petrov Galerkin (SUPG) finite element method. In the case of shape parameters, boundary conditions for sensitivities depend on the flow gradient at the boundary. Therefore, an accurate recovery of solution gradients is crucial to the success of shape sensitivity computations. In this work, solution gradients at boundary points are extracted using the Finite Node Displacement (FiND) method on which the finite element discretization is enriched locally via the insertion of nodes close to the boundary points. The normal derivative of the solution is then determined using finite differences. This approach to evaluate shape sensitivity boundary conditions is embedded in the continuous SEM. The methodology is applied to the flow past a cylinder in ground proximity. First, the proposed method is verified on a steady state problem. The computed sensitivity is compared to the actual change in the solution when a small perturbation is imposed to the shape parameter. Then, the study investigates the ability of the SEM to anticipate the unsteady flow response to changes in the ground to cylinder gap. A reduction of the gap causes damping of the vortex shedding while an increase amplifies the unsteadiness.     

Gravity-driven laminar film flow of power-law fluids along vertical walls

A theoretical analysis is presented which brings steady laminar film flow of power-law fluids within the framework of classical boundary layer theory. The upper part of the film, which consists of a developing viscous boundary layer and an external inviscid freestream, is treated separately from the viscous dominated part of the flow, thereby taking advantage of the distinguishing features of each flow region. It is demonstrated that the film boundary layer developing along a vertical wall can be described by a generalized Falkner-Skan type equation originally developed for wedge flow. An exact similarity solution for the velocity field in the film boundary layer is thus made available.Downstream of the boundary layer flow regime the fluid flow is completely dominated by the action of viscous shear, and fairly accurate solutions are obtained by the Von Karman integral method approach. A new form of the velocity profile is assumed, which reduces to the exact analytic solution for the fully-developed film. By matching the downstream integral method solution to the upstream generalized Falkner-Skan similarity solution, accurate estimates for the hydrodynamic entrance length are obtained. It is also shown that the flow development in the upstream region predicted by the approximate integral method closely corresponds to the exact similarity solution for that flow regime. An analytical solution of the resulting integral equation for the Newtonian case is compared with previously published results.     

Solution of flow around complex‐shaped surfaces by an immersed boundary‐body conformal enrichment method

Solving the flow around objects with complex shapes may involve extensive meshing work that has to be repeated each time a change in the geometry is needed. Time consuming meshing can be avoided when the solution algorithm can tackle grids that do not fit the shape of immersed objects. This work presents applications of a recently proposed immersed boundary—body conformal enrichment method to the solution of the flow around complex shaped surfaces such as those of a metallic foam matrix. The method produces solutions of the flow satisfying accurately Dirichlet boundary conditions imposed on the immersed fluid/solid interface. The boundary of immersed objects is defined using a level‐set function, and the finite element discretization of interface elements is enriched with additional degrees of freedom, which are eliminated at element level. The method is first validated in the case of flow problems for which reference solutions on body‐conformal grids can be obtained: flow around an array of spheres and flow around periodic arrays of cylinders. Then, solutions are shown for the more complex flow inside a metallic foam matrix. A multiscale approach combining the solution at the pore level by the immersed boundary method and the macro‐scale solution with simulated permeability is used to solve actual experimental configurations. The computed pressure drop as a function of the flow rate on the macro scale configuration replicating two experimental setups is compared with the experimental data for various foam thicknesses. Copyright © 2011 National Research Council Canada     

Numerical Modeling of the Development of a Streaky Structure in a Boundary Layer at a High Freestream Turbulence Level

The disturbances generated by external turbulence in the boundary layer on a flat plate set suddenly in motion are determined. A turbulent flow calculated by direct numerical simulation is taken as the initial conditions. The solution obtained simulates the initial stage of laminar-turbulent transition in the flat-plate boundary layer at a high turbulence level in the oncoming flow. The solution makes it possible to estimate the effects of different factors, such as nonstationarity, nonlinearity, and the parameters of the freestream velocity fluctuation spectrum, on disturbance enhancement in the boundary layer.     

Uniqueness of the self-similar solutions of three-dimensional laminar boundary-layer equations

The equations of the three-dimensional laminar boundary layer on lines of flow outflow and inflow are studied for conical outer flow under the assumption that the Prandtl number and the productρμ are constant. It is shown that in the case of a positive velocity gradient of the secondary flow (α₁>0) the additional conditions which result from the physical flow pattern determine a unique solution of the system of boundary-layer equations. For a negative velocity gradient of the secondary flow (α₁≤0) these conditions are satisfied by two solutions. An approximate solution is obtained for the boundary layer equations which is in rather good agreement with the numerical integration results. Compressible gas flow in a three-dimensional laminar boundary layer is described by a system of nonlinear differential equations whose solution is not unique for given boundary conditions. Therefore additional conditions resulting from the physical pattern of the gas flow are imposed on the resulting solution. In the solution of problems with a negative pressure gradient these additional conditions are sufficient for a unique selection of the solution of the boundary-layer equations. However, in the case of a positive pressure gradient the solution of the boundary-layer equations satisfying the boundary and additional conditions may not be unique. In particular, in [1] in a study of a three-dimensional laminar boundary layer in the vicinity of the stagnation point it was shown that for $$c = {{\frac{{\partial v_e }}{{\partial y}}} \mathord{\left/ {\vphantom {{\frac{{\partial v_e }}{{\partial y}}} {\frac{{\partial u_e }}{{\partial x}}}}} \right. \kern-\nulldelimiterspace} {\frac{{\partial u_e }}{{\partial x}}}} > 0$$ the solution is unique, while for c<0 there are two solutions. In the present paper we study the question of the uniqueness of the self-similar solution of the three-dimensional laminar boundary-layer equations on lines of flow outflow and inflow for a conical outer flow.     

The two-dimensional transmission and reflection of a solitary wave

In the present paper, limited to the discussion of weak non-linear shallow water waves, the transmission and reflection of a planar soliton on a two-dimensional structure are considered. The whole flow field is divided mainly into two subfields. One is in the vicinity of the structure, called the inner field; the other is far from the structure, called the outer field. In the outer field, according to its definition, the influence of the structure on the flow is negligible; to the order O(α, β) the governing equation for the flow is replaced by the Boussinesq equation. In the inner field the effect of the structure on the flow is significant, so the full Laplace equation is adopted as the governing equation for the flow field. Then the matched asymptotic expansion method is employed to connect smoothly the inner and outer solutions. Owing to the irregularity of the bottom of the structure, the boundary element method is incorporated. As an example, the case in which the incoming wave is a solitary wave is calculated and the time histories of transmitted and reflected waves are plotted.