Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper, a new analytical method of symplectic system, Hamiltonian system, is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain. In the system, the fundamental problem is reduced to an eigenvalue and eigensolution problem. The solution and boundary conditions can be expanded by eigensolutions using adjoint relationships of the symplectic ortho-normalization between the eigensolutions. A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space. The results show that fundamental flows can be described by zero eigenvalue eigensolutions, and local effects by nonzero eigenvalue eigensolutions. Numerical examples give various flows in a rectangular domain and show effectiveness of the method for solving a variety of problems. Meanwhile, the method can be used in solving other problems.     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper,a new analytical method of symplectic system.Hamiltonian system,is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain.In the system,the fundamental problem is reduced to all eigenvalue and eigensolution problem.The solution and boundary conditions call be expanded by eigensolutions using ad.ioint relationships of the symplectic ortho-normalization between the eigensolutions.A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space.The results show that fundamental flows can be described by zero eigenvalue eigensolutions,and local effects by nonzero eigenvalue eigensolutions.Numerical examples give various flows in a rectangular domain and show effectivenees of the method for solving a variety of problems.Meanwhile.the method can be used in solving other problems.     

A theorem for a shear stress-free plane boundary in Stokes flow

A theorem for non-axisymmetrical Stokes flow about a shear stress-free plane boundary is established by using a representation for the velocity and pressure fields for the same flow in terms of biharmonic and harmonic functions. A corollary of the theorem is derived which gives the axisymmetrical Stokes flow in terms of the Stokes function about the same boundary. The formulae for drag and torque on the boundary are also given. A few illustrative examples are presented.     

Stokes paradox for power-law flow around a cylinder

A simple analysis for power-law fluids shows that the Stokes paradox for creeping flow around a cylinder is removed for shear-thinning (n < 1) but not for shear-thickening (n 1) fluids. An approximate drag value is found for n < 1 and is compared with computed results.     

Stokes流问题中的辛本征解方法

通过引入哈密顿体系,将二维Stokes流问题归结为哈密顿体系下的本 征值和本征解问题. 利用辛本征解空间的完备性,建立一套封闭的求解问题方法. 研究结果 表明零本征值本征解描述了基本的流动,而非零本征值本征解则显示着端部效应影响特点. 数值算例给出了辛本征值和本征解的一些规律和具体例子. 这些数值例子说明了端部非规则 流动的衰减规律. 为研究其它问题提供了一条路径.     

矩形空腔内Stokes流的状态空间有限元法

基于Hellinger-Reissner二类变分原理,从平面Stokes流问题的平衡方程、连续性要求和边界条件出发,得到相应的Hamilton函数,建立Hamilton正则方程后,采用分离变量法对场变量进行离散求解:在x方向采用有限元插值,在y方向采用状态空间法给出控制坐标方向的解析解。计算过程中的指数矩阵均采用精细积分法求解,使得本文算法具有高效率、高精度、对步长不敏感的优点。通过对侧边自由液面边界条件的单板驱动矩形空腔Stokes流问题的求解,得到与文献相同的结果,从而验证了本文方法的有效性。本文旨在将弹性力学状态空间有限元法的思想引入到低雷诺数流体力学中,为Hamilton体系下研究复杂边界Stokes流问题提供新的途径。     

On Stokes slip flow through a transversely wavy channel

The Stokes flow through a wavy or corrugated channel with surface slip is studied. The correct Navier's partial slip condition is applied and perturbation solutions about the small amplitude to channel width ratio are obtained. As in Stokes slip flow over a sphere, the resistance is not zero even when slip is infinite. The resistance (due to the interaction of waviness and slip) is larger when the corrugations of the two plates are out of phase than that when they are in phase.     

Optimum volume profile in axi-symmetric Stokes flow

In this paper, the problem of optimum volume profile under constant drag constraint when uniform stream is parallel to the axis of symmetry and flow is governed by Stokes equations has been tackled. Here, we take up a class of bodies to be of given maximum cross-section with fore and aft symmetry about this section. The possible shape under the stationary value drag has been obtained by making use of method of extremals [Fox C (1950) An introduction to the calculus of variations. Oxford University Press, Oxford; Elsgolc EC (1962) Differential equations and the calculus of variations. Pergamon Press; Sagan H (1969) Introduction to the calculus of Variations. MacGrawhill]. It has been found that the body profile possesses conical front and rear ends. The numerical value of the optimized volume has also been calculated for the profile and compared with some known values.     

Shock tunnel flow visualization using planar laser-induced fluorescence imaging of NO and OH

Temporal sequences of planar laser-induced fluorescence (PLIF) images of several high-speed, transient flowfields created in a reflection-type shock tunnel facility were acquired. In each case, the test gas contained either nitric oxide or the hydroxyl radical, the fluorescent species. The processes of shock reflection from an endwall with a converging nozzle and of underexpanded free jet formation were examined. A comparison was also made between PLIF imaging and shadow photography. The investigation demonstrated some of the capabilities of PLIF imaging diagnostics in complex, transient, hypersonic flowfields, including those with combustion.Nomenclature A spontaneous emission rate - A _las cross sectional area of laser sheet - B laser absorption rate - C _opt constant dependent on optical arrangement, collection efficiency, etc. - D nozzle throat diameter - E _p laser pulse energy - f _J Boltzmann fraction of absorbing state - g spectral convolution of laser and absorption lineshapes - k Boltzmann constant - M _s incident shock Mach number - N noise, root-mean-square signal fluctuation - P static pressure - P ₁ initial pressure of test gas in shock tube - P _a free jet ambient pressure - P _s stagnation pressure - Q electronic quenching rate of excited state - S PLIF signal - t time between shock reflection and image acquisition - T static temperature - T _s stagnation temperature - _a mole fraction of absorbing species     

不可压缩Stokes流动的PSPG无网格法

将应用于有限元法的Pressure-Stabilizing/Petrov-Galerkin(PSPG)稳定化机制引入到无网格法中,有效消除了由于速度和压力的插值模式违反LBB条件而导致的压力场的虚假振荡。采用与有限元法耦合的连续掺混法(Continuous Blending Method)施加本质边界条件,使得边界条件不仅在边界节点上而且在整条边界上都得到严格满足。给出了两个典型算例的数值模拟结果,表明了所建议无网格法模拟不可压缩Stokes流动的有效性。     

Slow visco-elastic flow past a cylinder in a rectangular channel

This paper is concerned with the flow of a visco-elastic liquid through a rectangular channel containing a cylindrical obstruction placed either in a symmetric or asymmetric position with respect to the centre of the channel. Numerical predictions of the flow are obtained using a well established finite element Galerkin mixed formulation. The influence of elasticity on the streamline pattern is found to be negligible, and one only observes changes due to different geometries, i.e. relative positions of cylinder and channel. However, both elasticity and a variable viscosity are found to have significant effects on the forces exerted on the cylinder.     

A domain decomposition method for modelling Stokes flow in porous materials

An algorithm is presented for solving the Stokes equation in large disordered two‐dimensional porous domains. In this work, it is applied to random packings of discs, but the geometry can be essentially arbitrary. The approach includes the subdivision of the domain and a subsequent application of boundary integral equations to the subdomains. This gives a block diagonal matrix with sparse off‐block components that arise from shared variables on internal subdomain boundaries. The global problem is solved using a biconjugate gradient routine with preconditioning. Results show that the effectiveness of the preconditioner is strongly affected by the subdomain structure, from which a methodology is proposed for the domain decomposition step. A minimum is observed in the solution time versus subdomain size, which is governed by the time required for preconditioning, the time for vector multiplications in the biconjugate gradient routine, the iterative convergence rate and issues related to memory allocation. The method is demonstrated on various domains including a random 1000‐particle domain. The solution can be used for efficient recovery of point velocities, which is discussed in the context of stochastic modelling of solute transport. Copyright © 2002 John Wiley & Sons, Ltd.     

Stokes flow due to fundamental singularities before a plane boundary

A representation for the velocity and pressure fields in three-dimensional Stokes flow was presented in terms of a biharmonic function A and a harmonic function B. This representation was used to establish a general theorem for the calculation of Stokes flow due to fundamental singularities in a region bounded by a stationary no-slip plane boundary. Collins‘s theorem for axisymmetric Stokes flow before a rigid plane follows as a special case of the theorem. A few illustrative examples are given to show its usefulness.     

A singular finite element for Stokes flow: The stick–slip problem

Abrupt changes in boundary conditions in viscous flow problems give rise to stress singularities. Ordinary finite element methods account effectively for the global solution but perform poorly near the singularity. In this paper we develop singular finite elements, similar in principle to the crack tip elements used in fracture mechanics, to improve the solution accuracy in the vicinity of the singular point and to speed up the rate of convergence. These special elements surround the singular point, and the corresponding field shape functions embody the form of the singularity. Because the pressure is singular, there is no pressure node at the singular point. The method performs well when applied to the stick–slip problem and gives more accurate results than those from refined ordinary finite element meshes.     

Stokes flow before plane boundary with mixed stick-slip boundary conditions

A general theorem for the Stokes flow over a plane boundary with mixed stick-slip boundary conditions is established. This is done by using a representation for the velocity and pressure fields in the three-dimensional Stokes flow in terms of a biharmonic function and a harmonic function. The earlier theorem for the Stokes flow due to fundamental singularities before a no-slip plane boundary is shown to be a special case of the present theorem. Furthermore, in terms of the Stokes stream function, a corollary of the theorem is also derived, providing a solution to the problem of the axisymmetric Stokes flow along a rigid plane with stick-slip boundary conditions. The formulae for the drag and torque exerted by the fluid on the boundary are established. An illustrative example is given.     

The numerical solution of Stokes flow in a domain with re-entrant boundaries by the boundary element method

Numerical solutions are presented for two-dimensional low Reynolds number flow in a rotating tank with stationary barriers. The boundary element method is employed, assuming straight panels and quadratic source distribution. The feasibility of repositioning the nodes as a way to minimize the error is explored. A stretching parameter places smaller elements near the re-entrant regions. Elementary error analysis shows uniform improvement in the solution with stretching. The changing eddy pattern for different numbers and sizes of the barriers is compared with experimental results.     

Three-dimensional spectral approximations to Stokes flow between eccentrically rotating cylinders

A three-dimensional spectral algorithm for the solution of Stokes flow between eccentrically rotating cylinders is described. Included in the model are pressure boundary conditions at the two ends of the finite length cylinders and the effect of a fluid line source on the inner cylinder. A comparison of results for the load and couple on the inner cylinder is made with those available from lubrication theory in the absence of a line source. Good agreement is shown for long, short and finite journal bearings when the various geometrical assumptions inherent in the lubrication analysis are satisfied.     

Reliable validation based on optical flow visualization for CFD simulations

A reliable validation based on the optical flow visualization for numerical simulations of complex flowfields is addressed in this paper. Several test cases, including two-dimensional, axisymmetric and three-dimensional flowfields, were presented to demonstrate the effectiveness of the validation and gain credibility of numerical solutions of complex flowfields. In the validation, images of these flowfields were constructed from numerical results based on the principle of the optical flow visualization, and compared directly with experimental interferograms. Because both experimental and numerical results are of identical physical representation, the agreement between them can be evaluated effectively by examining flow structures as well as checking discrepancies in density. The study shows that the reliable validation can be achieved by using the direct comparison between numerical and experiment results without any loss of accuracy in either of them.     

Stokes flow of a rotating sphere around the axis of a circular orifice or a circular disk

This paper presents an infinite series solution to the creeping flow equations for the axisymmetric motion of a sphere of arbitrary size rotating in a quiescent fluid around the axis of a circular orifice or a circular disk whose diameters are either larger or smaller than that of the sphere. Numerical tests of the convergence are passed and the comparison with the exact solution and other computational results shows an agreement to five significant figures for the torque coefficients in both cases. The torque coefficients are obtained for the sphere located up to a position tangent to the wall plane containing either the orifice or the disk. It is concluded that the torque coefficients of the sphere and the disk are monotonically increasing with the decrease of the distance from the disk or the orifice plane in both cases.     

Hamilton体系下环扇形域的Stokes流动问题

基于极坐标下Stokes流的基本方程,将径向坐标模拟为时间坐标,推导了Hamilton体系下Stokes流动问题的对偶方程,采用本征向量展开法对环扇形域Stokes流动问题进行了分析,并给出了相应的实际算例,其结果说明了本文方法的有效性。