An h4 accurate vorticity-velocity formulation for calculating flow past a cylinder

A method of solution for the two-dimensional Navier-Stokes equations for incompressible flow past a cylinder is given in which the euquation of continuity is solved by a step-by-step integration procedure at each stage of an interative process. Thus the formulation involves the solution of one first-order and one second-order equation for the velocity components, together with the vorticity transport equation. the equations are solved numerically by h⁴-accurate methods in the case of steady flow past a circular cylinder in the Reynolds number range 10–100. Results are in satisfactory agreement with recent h⁴-accurate calculations. An improved approximation to the boundary conditions at large distance is also considered.     

不可压Navier-Stokes方程基于误差估算的网格自适应解法

首先导出了广义Stokes方程Petrov—Galerkin有限元数值解的当地事后误差估算公式;以非连续二阶鼓包（bump）函数空间为速度、压强误差的近似空间,该估算基于求解当地单元上的广义Stokes问题。然后,证明了误差估算值与精确误差之间的等价性。最后,将误差估算方法应用于Navier—Stokes环境,以进行不可压粘流计算中的网格自适应处理。数值实验中成功地捕获了多强度物理现象,验证了本文所发展的方法。     

Discontinuous element pressure gradient stabilizations for compressible Navier-Stokes equations based on local projections

A pressure gradient discontinuous finite element formulation for the compressible Navier-Stokes equations is derived based on local projections. The resulting finite element formulation is stable and uniquely solvable without requiring a B-B stability condition. An error estimate is Obtained.     

A two-order and two-scale computation method for nonselfadjoint elliptic problems with rapidly oscillatory coefficients

The purpose of this paper is to solve nonselfadjoint elliptic problems with rapidly oscillatory coefficients. A two-order and two-scale approximate solution expression for nonselfadjoint elliptic problems is considered, and the error estimation of the twoorder and two-scale approximate solution is derived. The numerical result shows that the presented approximation solution is effective.     

A class of Wilson arbitrary quadrilateral elements for an axisymmetric problem

A class of modified Wilson arbitrary quadrilateral nonconforming elements for an axisymmetric problem is proposed. Their convergence is proven by means of the strong patch test. The structure of this finite element class is investigated. Thus a general method of axisymmetric nonconforming elements with convergence properties is presented. Project supported by the Foundation of the Chinese Ministry of Machine-Building Industry     

A parallel two-level finite element method for the Navier-Stokes equations

Based on domain decomposition, a parallel two-level finite element method for the stationary Navier-Stokes equations is proposed and analyzed. The basic idea of the method is first to solve the Navier-Stokes equations on a coarse grid, then to solve the resulted residual equations in parallel on a fine grid. This method has low communication complexity. It can be implemented easily. By local a priori error estimate for finite element discretizations, error bounds of the approximate solution are derived. Numerical results are also given to illustrate the high efficiency of the method.     

A kind of distortion of mean velocity profile in pipe Poiseuille flow and its stability behaviour

This paper presents a kind of distortion of Hagen-Poiseuille velocity profile in pipe Poiseuille flow. This distortion can be regarded as a general expression of the influence on the mean flow by nonlinear interaction of various components of axisymmetric perturbations. Through the investigation of the stability behaviour of this velocity profile, this paper obtains unstable result induced by axisymmetric perturbations for the first time, and thus presents a new possible approach which leads to instability of Hagen-Poiseuille flow.     

A projection scheme for Navier-Stokes with variable viscosity and natural boundary condition

Combiningthe Navier-Stokes systems with Neumann (or natural) boundary condition to characterize a fluid flow is frequent. The popular projection (or pressure correction) methods inspired by Chorin and Temam are not well adapted to such boundary condition, which translate in loss of accuracy. If some alternative projection methods have been proposed to reduce the accuracy loss due to the Neumann condition in case of Newtonian fluids, little has been proposed for generalized Newtonian fluids. In this work, we propose two methods derived from the incremental pressure correction projection that can be used for fluids with inhomogeneous or variable viscosity with natural boundary condition. Both time and space accuracy of the methods will be illustrated using a manufactured solution.     

液浮转子微陀螺筒间流动阻力矩参数影响

建立简化计算模型,将一种具有空间指数收敛特性的谱元-Fourier法引入复杂的液浮转子微陀螺筒间流动的直接数值模拟。对于柱坐标系下不可压缩Navier-Stokes方程,在子午面的径向和轴向采用Galerkin谱元法,进行标准化的Gauss-Lobatto-Legendre节点展开,在周向采用Fourier级数展开法。设定基准算例,计算分析了高径比、旋转速度和黏性系数等参数对内筒侧面和内筒端面阻力矩的影响。结果表明,两种阻力矩都随这些参数的增大而增大,且对总阻力矩的贡献量级相当。     

A streamline diffusion nonconforming finite element method for the time-dependent linearized Navier-Stokes equations

A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.     

Numerical simulation of fluid flow of two rotating side-by-side circular cylinders by Lattice Boltzmann method

A numerical investigation of the two-dimensional laminar flow around side-by-side rotating circular cylinders using Lattice Boltzmann method is conducted. The effects of variation of rotational speed ratio β and different gap spacings g* at Reynolds number of 100 are studied. A various range of rotational speed ratio 0 ≤ β ≤ 2 for four different gap spacings of 3, 1.5, 0.7 and 0.2 are investigated. Flow conditions and its characteristics, such as lift and drag coefficients and Strouhal number, is studied. The results indicated that as β increases, the flow changes its condition from periodic to steady after a critical rotational speed. Results also indicated that variation of the gap spacing and rotational speed has significant effect on wake pattern. Wake pattern in turn has significant effect on the Strouhal number. Finally, the result is compared with experimental and other numerical data.     

Vibration modelling of structural networks using a hybrid finite element/wave and finite element approach

The vibration modelling of waveguide structures is considered. These structures comprise waveguides connected via joints. Traditionally, analytical models of the wave behaviour of such structures can be developed if they are simple (beams or rods connected at point joints, etc.). However, if the waveguides are of complicated constructions (truss-like, layered media, etc.) or the joints are complicated (e.g. of significant physical dimensions), obtaining the wave characteristics might be a formidable task. In this paper, such structures are modelled using a hybrid finite element/wave and finite element (FE/WFE) approach. The waveguides are modelled using the WFE method and thus their wave characteristics are obtained regardless of the complexity of their cross-section. The joints are modelled using standard FE, and the WFE and FE models are coupled to yield the scattering properties of the joints. The propagation and scattering models are assembled to describe the behaviour of the structure using relatively small models, while also providing information for other applications such as structure-borne sound, statistical energy analysis, etc. Numerical examples are presented to illustrate the approach.     

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

Reference values for drag and lift of a two‐dimensional time‐dependent flow around a cylinder

This paper presents a numerical study of a two‐dimensional time‐dependent flow around a cylinder. Its main objective is to provide accurate reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time. In addition, the accuracy of these values obtained with different time stepping schemes and different finite element methods is studied. Copyright © 2004 John Wiley & Sons, Ltd.     

亚临界雷诺数下圆柱受迫振动的数值研究

在Navier-Stokes方程和k-ω湍流模型的基础上,利用流线迎风有限元方法结合ALE动网格技术对亚临界雷诺数下的圆柱受迫振动问题开展了数值模拟研究。本文的数值模型成功模拟了Re=5000条件下,圆柱发生受迫振动时尾迹区内的2S,2P和P+S尾流模式;对Re=10000情况下,无量纲振幅分别为0.3,0.4,0.5的圆柱受迫振动问题开展了数值模拟,分析了给定振幅条件下圆柱受力随振动频率的变化关系以及受迫振动的锁定区间。以上数值计算结果与Gopalkrishnan (1993)的实验结果基本符合。研究结果表明,二维数值模型能够基本正确地反映出圆柱发生受迫振动时的涡激振动特性以及有关的受力变化趋势,为今后进一步开展三维数值分析工作奠定了基础。     

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the local- truncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases. The project supported by the China NKBRSF (2001CB409604) The English text was polished by Yunming Chen     

An upwind finite element scheme for high-Reynolds-number flows

A new upwind finite element scheme for the incompressible Navier-Stokes equations at high Reynolds number is presented. The idea of the upwind technique is based on the choice of upwind and downwind points. This scheme can approximate the convection term to third-order accuracy when these points are located at suitable positions. From the practical viewpoint of computation, the algorithm of the pressure Poisson equation procedure is adopted in the framework of the finite element method. Numerical results of flow problems in a cavity and past a circular cylinder show excellent dependence of the solutions on the Reynolds number. The influence of rounding errors causing Karman vortex shedding is also discussed in the latter problem.     

The dynamic plastic response characteristics of circular beam

In this paper the problem of a circular beam subjected to radial impact by a rigid mass at its tip in its own plane is investigated on the basis of rigid-perfectly plastic assumption. The analytical solution of the particle velocities is obtained as the function of travelling plastic hinge location. By analysing the solution, some special properties of circular beam problem are found.     

A problem of heat and mass transfer: Proof of the existence condition by a finite difference method

We solve by a finite difference method a system of simultaneous non-linear partial differential equations which modelizes the transfer of heat and mass when a fluid evaporates from the hot wall and condenses on the cold wall of an upright rectangular cavity. The need to verify a certain condition associating the physical parameters of the fluid for the existence of steady state solutions is proved.     

A design of residual error estimates for a high order BDF‐DGFE method applied to compressible flows

We deal with the numerical solution of the non‐stationary compressible Navier–Stokes equations with the aid of the backward difference formula – discontinuous Galerkin finite element method. This scheme is sufficiently stable, efficient and accurate with respect to the space as well as time coordinates. The nonlinear algebraic systems arising from the backward difference formula – discontinuous Galerkin finite element discretization are solved by an iterative Newton‐like method. The main benefit of this paper are residual error estimates that are able to identify the computational errors following from the space and time discretizations and from the inexact solution of the nonlinear algebraic systems. Thus, we propose an efficient algorithm where the algebraic, spatial and temporal errors are balanced. The computational performance of the proposed method is demonstrated by a list of numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.