基于Uzawa算法的弹塑性扭转问题的局部微分求积法

研究了由椭圆变分不等式描述的弹塑性扭转问题,构造了基于Uzawa算法的局部微分求积法,给出了数值算例,通过与有限元方法的比较,说明了方法的有效性。     

An efficient matrix solver for finite-element analysis of non-Newtonian fluid flow problems

When the Galerkin finite-element method with a nine-node isoparametric Lagrangian element is applied to solve non-Newtonian fluid flow problems, a considerable amount of computing time is required to solve the discretized non-linear system of equations by Newton's method. A method proposed by Broyden has been modified to compute the Jacobian matrix associated with Newton's method. This modified Broyden's method can be combined with the frontal method to efficiently solve the linearized finite-element equations during the iteration. Numerical results of a sample problem concerning the determination of the pressure-drop/flow-rate relationship for power-law fluids in rectangular ducts show that the application of this new method can reduce computing time substantially.     

Performance analysis of IDEAL algorithm for three‐dimensional incompressible fluid flow and heat transfer problems

Recently, an efficient segregated algorithm for incompressible fluid flow and heat transfer problems, called inner doubly iterative efficient algorithm for linked equations (IDEAL), has been proposed by the present authors. In the algorithm there exist inner doubly iterative processes for pressure equation at each iteration level, which almost completely overcome two approximations in SIMPLE algorithm. Thus, the coupling between velocity and pressure is fully guaranteed, greatly enhancing the convergence rate and stability of solution process. However, validations have only been conducted for two‐dimensional cases. In the present paper the performance of the IDEAL algorithm for three‐dimensional incompressible fluid flow and heat transfer problems is analyzed and a systemic comparison is made between the algorithm and three other most widely used algorithms (SIMPLER, SIMPLEC and PISO). By the comparison of five application examples, it is found that the IDEAL algorithm is the most robust and the most efficient one among the four algorithms compared. For the five three‐dimensional cases studied, when each algorithm works at its own optimal under‐relaxation factor, the IDEAL algorithm can reduce the computation time by 12.9–52.7% over SIMPLER algorithm, by 45.3–73.4% over SIMPLEC algorithm and by 10.7–53.1% over PISO algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

高黏性流动有限元模拟的迭代稳定分步算法

分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法,它能在 保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实 了此特点,该方法被有效应用于充填流动过程的数值模拟.     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 John Wiley & Sons, Ltd.     

Perturbation solution to 3-D nonlinear supercavitating flow problems

A 3-D nonlinear problem of supercavitating flow past an axisymmetric body at a small angle of attack is investigated by means of the perturbation method and Fourier-cosine-expansion method. The first three order perturbation equations are derived in detail and solved numerically using the boundary integral equation method and iterative techniques. Computational results of the hydrodynamic characteristics and cavity shapes of each order are presented for nonaxisymmetric supercavitating flow past cones with various apex-angles at different cavitation numbers. The numerical results are found in good agreement with experimental data. The project supported by the National Natural Science Foundation of China     

Non-swirling axisymmetric flow of the ω-D fluid

In this paper a constitutive equation for the extra stress tensor τ is considered, which can be used for steady axisymmetric flows when u _ϕ=0 (non-swirling). It is explicit with coefficients which, in case of incompressibility, depend only upon three invariants. As such it should prove useful in numerical calculations, especially so if used in its simplest form, namely the quasi-Newtonian fluid. Here, a purely viscous code can be used. The constitutive equation is self-consistent and receives additional justification from the fact that any constitutive equation is bound to result in the form proposed if the flow is a constant stretch history flow. Received: 30 March 1998 Accepted: 6 August 1998     

Spectral methods and their implementation to solution of aerodynamic and fluid mechanic problems

Fundamental concepts underlying spectral collocation methods, especially pertaining to their use in the solution of partial differential equations, are outlined. Theoretical accuracy results are reviewed and compared with results from test problems. A number of practical aspects of the construction and use of spectral methods are detailed, along with several solution schemes which have found utility in applications of spectral methods to practical problems. Results from a few of the successful applications of spectral methods to problems of aerodynamic and fluid mechanic interest are then outlined, followed by a discussion of the problem areas in spectral methods and the current research under way to overcome these difficulties.     

Fluid flow of incompressible viscous fluid through a non-linear elastic tube

The study of viscous flow in tubes with deformable walls is of specific interest in industry and biomedical technology and in understanding various phenomena in medicine and biology (atherosclerosis, artery replacement by a graft, etc) as well. The present work describes numerically the behavior of a viscous incompressible fluid through a tube with a non-linear elastic membrane insertion. The membrane insertion in the solid tube is composed by non-linear elastic material, following Fung’s (Biomechanics: mechanical properties of living tissue, 2nd edn. Springer, New York, 1993) type strain–energy density function. The fluid is described through a Navier–Stokes code coupled with a system of non linear equations, governing the interaction with the membrane deformation. The objective of this work is the study of the deformation of a non-linear elastic membrane insertion interacting with the fluid flow. The case of the linear elastic material of the membrane is also considered. These two cases are compared and the results are evaluated. The advantages of considering membrane nonlinear elastic material are well established. Finally, the case of an axisymmetric elastic tube with variable stiffness along the tube and membrane sections is studied, trying to substitute the solid tube with a membrane of high stiffness, exhibiting more realistic response.     

An efficient algorithm for solving the incompressible fluid flow equations

The present paper reports on a modified pressure implicit predictor corrector type scheme for solving the flow governing equations, in which a consistent formulation is combined with a multi-grid solver for the pressure correction. In addition a parabolic sublayer (PSL) approach for the treatment of the flow in the vicinity of solid walls is critically evaluated in terms of accuracy and computational efficiency. The lid-driven cavity flow is chosen as the test case and results are presented for Reynolds numbers ranging from 100 to 1000. Predictions with the proposed scheme indicate substantial computational savings and fairly good agreement when compared with previous work. The PSL approach reduces the computing time, but with increasing Reynolds numbers the accuracy of the solutions tends to deteriorate.     

An Uzawa-type algorithm for the coupled Stokes equations

An Uzawa-type algorithm is designed for the coupled Stokes equations discretized by the mixed finite element method. The velocity solved by the presented algorithm is weakly divergence-free, which is different from the one solved by the common Uzawa method. Besides, an optimal relaxation parameter of the presented algorithm is provided.     

An iterative parallel algorithm of finite element method

In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.     

风沙两相流PIV测量算法研究

在风沙两相流图像特征的基础上，提出了一种基于模式识别动态聚类方法中$K$-均值 算法的数字面具(Digital Mask)自动生成算法来求解风沙两相流动. 简化了传统生成Digital Mask过程中手动设置参数的操作，减少了人为误差，为批量处理风沙两相流PIV图像提供 了一种安全快捷的方法. 并将该算法应用于风沙两相流的PIV实际测量，分别得到了不同流 动状态下的气流、沙粒以及风沙气固两相流的速度场.     

Preconditioned finite element algorithms for 3D Stokes flows

Preconditioned conjugate gradient algorithms for solving 3D Stokes problems by stable piecewise discontinuous pressure finite elements are presented. The emphasis is on the preconditioning schemes and their numerical implementation for use with Hermitian based discontinuous pressure elements. For the piecewise constant discontinuous pressure elements, a variant implementation of the preconditioner proposed by Cahouet and Chabard for the continuous pressure elements is employed. For the piecewise linear discontinuous pressure elements, a new preconditioner is presented. Numerical examples are presented for the cubic lid-driven cavity problem with two representative elements, i.e. the Q2-PO and the Q2-P1 brick elements. Numerical results show that the preconditioning schemes are very effective in reducing the number of pressure iterations at very reasonable costs. It is also shown that they are insensitive to the mesh Reynolds number except for nearly steady flows (Re_m → 0) and are almost independent of mesh sizes. It is demonstrated that the schemes perform reasonably well on non-uniform meshes.     

A finite-element approach in stream-tube method for solving fluid and solid mechanics problems

This paper presents a theoretical formulation in which the stream-tube method (STM) is examined through a variational approach for solving solid strain and fluid flow problems with finite elements. The analysis considers a reference domain, used as computational domain, related to the physical domain by an unknown transformation function to be determined numerically. Mass conservation is automatically verified by STM. The variational approach leads to eliminate the pressure in fluid problems and avoids to set up a mixed displacement–pressure procedure in the case of incompressible solids. Examples are given for fluid flows, applications and comparisons are also provided in the bending problem in elasticity.     

大体积混凝土不稳定温度场求解的迭代算法

在大体积混凝土不稳定温度场计算中，无论是向前差分法、中点差分法还是向后差分法都会产生计算误差．根据拉格朗日中值定理及各节点温差和温度变化率的互不耦合假设，提出一种迭代算法，可以在适当增加差分步长的情况下，减少计算量并保证求解精度，或在保持常规差分步长的情况下，增加少许计算量而提高求解精度，使数值计算结果逼近理论解，是一个便于工程应用的有效计算方法．     

An iterative method for solution of a boundary value problem in non-newtonian fluid flow

The boundary value problem

arises in boundary layer equations for the steady flow of a power-law fluid over an impermeable, semi-infinite flat plane. The parameter μ is equal to $\text{1}\text{n}$ where n is the exponent of the strain rate in the expression for the shear stress. We develop and prove the convergence of an iterative method for the solution of the given boundary value broblem for dilatant fluids (0 < μ <1). The iterative method can be easily implemented computationally. An added feature of our technique is that it accurately yields y(0), an important parameter which is related to the drag at the plate. The iterative method works well computationally not only for 0 < μ < 1 but for the range 1 < μ < 4 (pseudoplastic fluids with 1 > n > $\text{1}\text{4}$), as well.     

Boundary conformed co-ordinate systems for selected two-dimensional fluid flow problems. Part I: Generation of BFGs

In this paper the generation of general curvilinear co-ordinate systems for use in selected two-dimensional fluid flow problems is presented. The curvilinear co-ordinate systems are obtained from the numerical solution of a system of Poisson equations. The computational grids obtained by this technique allow for curved grid lines such that the boundary of the solution domain coincides with a grid line. Hence, these meshes are called boundary fitted grids (BFG). The physical solution area is mapped onto a set of connected rectangles in the transformed (computational) plane which form a composite mesh. All numerical calculations are performed in the transformed plane. Since the computational domain is a rectangle and a uniform grid with mesh spacings Δξ = Δη = 1 (in two-dimensions) is used, the computer programming is substantially facilitated. By means of control functions, which form the r.h.s. of the Poisson equations, the clustering of grid lines or grid points is governed. This allows a very fine resolution at certain specified locations and includes adaptive grid generation. The first two sections outline the general features of BFGs, and in section 3 the general transformation rules along with the necessary concepts of differential geometry are given. In section 4 the transformed grid generation equations are derived and control functions are specified. Expressions for grid adaptation arc also presented. Section 5 briefly discusses the numerical solution of the transformed grid generation equations using sucessive overrelaxation and shows a sample calculation where the FAS (full approximation scheme) multigrid technique was employed. In the companion paper (Part II), the application of the BFG method to selected fluid flow problems is addressed.     

A three-dimensional fictitious domain method for incompressible fluid flow problems

A new Galerkin finite element method for the solution of the Navier–Stokes equations in enclosures containing internal parts which may be moving is presented. Dubbed the virtual finite element method, it is based upon optimization techniques and belongs to the class of fictitious domain methods. Only one volumetric mesh representing the enclosure without its internal parts needs to be generated. These are rather discretized using control points on which kinematic constraints are enforced and introduced into the mathematical formulation by means of Lagrange multipliers. Consequently, the meshing of the computational domain is much easier than with classical finite element approaches. First, the methodology will be presented in detail. It will then be validated in the case of the two-dimensional Couette cylinder problem for which an analytical solution is available. Finally, the three-dimensional fluid flow inside a mechanically agitated vessel will be investigated. The accuracy of the numerical results will be assessed through a comparison with experimental data and results obtained with a standard finite element method. © 1997 John Wiley & Sons, Ltd.