A semi-implicit three-step method based on SUPG finite element formulation for flow in lid driven cavities with different geometries

A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) stabilization scheme is used for the formulation of the Navier-Stokes equations. For the spatial discretization, the convection term is treated explicitly, while the viscous term is treated implicitly, and for the temporal discretization, a three-step method is employed. The present method is applied to simulate the lid driven cavity problems with different geometries at low and high Reynolds numbers. The results compared with other numerical experiments are found to be feasible and satisfactory.     

Numerical investigation of flow characteristics around two side-by-side cylinders by immersed boundary-lattice Boltzmann flux solver

In this paper, to study the characteristics of the flow in a laminar regime, an immersed boundary-lattice Boltzmann flux solver (IB-LBFS) is applied to numerically simulate the unsteady viscous flows around two fixed and rotating circular cylinders in side-by-side arrangement. This method applies finite volume discretization to solve the macroscopic governing equations with the flow variables defined at cell centers. At the cell interface, numerical fluxes are physically evaluated by a local lattice Boltzmann solution. In addition, the no-slip boundary condition is accurately imposed by using the implicit boundary condition-enforced immersed boundary method. Due to the simplicity and high efficiency of IB-LBFS on non-uniform grids, it is suitable for simulating fluid flows with complex geometries and moving boundaries. Firstly, numerical simulations of laminar flow past two side-by-side cylinder are performed with different gap spacings at Reynolds numbers of 100 and 200. The simulation results show that a small gap spacing induces a biased flow and forms an irregular big wake behind two cylinders at a low Reynolds number. As the gap spacing increases, an in-phase or anti-phase flow is observed. Then, the effects of the main important parameters on flow characteristics are analyzed for flow past two side-by-side rotating cylinders, including the rotational speed, Reynolds number, and gap spacing. As the rotational speed is increased, the numerical results illustrate that unsteady wakes are suppressed and the flow becomes steady. As the gap spacing is increased, two separate vortex streets behind each cylinder are formed with a definite phase relationship and single shedding frequency.     

Inverstigation on the Separated Turbulent Flow Field in Dual Rectangular Jets

In the present paper,the flow field of dual rectangular jets was numerically simulated by solving the full Reynolds averaged Navier-Stokes equations,where the RNG κ-ε model and the finite volume method were used.The flow structure in dual rectangular jets and the effects of the velocity were investigated.The numerical results agree qualitatively with the experimental data.     

Numerical simulation of interaction between laminar flow and elastic sheet

A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.     

Negative Pressure Consolidation of Silt Inside Bucket Foundation

A series of model experiments of bucket foundations concerning suction installation and negative pressure consolidation in saturated silt were carried out in a cube steel bin at Tianjin University. The experimental results show that the silt inside the bucket has been strengthened by negative pressure, and the strengthening effect decreases with the increase of the distance from the bucket. A three-dimensional numerical model of the experiments was built by means of finite element software ABAQUS with fluid-solid coupling method. The results show that the bearing capacity of the silt inside the bucket foundation increases significantly at the former stage of negative pressure consolidation, while the increasing trend slows down over time. The rotation centers of the bucket foundation and the inner soil region tend to be closer to each other based on the consolidation. The bearing capacity of the bucket foundation is improved effectively with the increase of soil strength. The effects of negative pressure consolidation on the bearing capacity of bucket foundation were also illustrated by an actual offshore wind power project case.     

求解瞬态Naiver-Stokes问题的一种稳定化方法

在压力投影及多尺度变分法的基础上,本文针对瞬态Naiver-Stokes问题提出一种新的稳定化方法,此方法采用非协调且不满足inf-sup条件的速度及压力有限元空间。文中的稳定化方法具备一些优势:能避免高阶微分的计算,数据结构与单元边界无关;只在好的尺度上增加离散的速度项,对高雷诺数的流体很有效;不增加计算成本。     

火电厂耐热钢的非线性局部化蠕变损伤研究

针对火电厂中蒸汽管道与汽包一类承压部件常用耐热钢材 ,在经历长时间高温高压作用后 ,处于蠕应变Ⅲ阶段过程时 ,发生非线性局部化 (大变形 )蠕变损伤问题 ,给出损伤本构描述 ,提出非线性局部化蠕变损伤本构模型及其数值变分原理与有限元离散化形式 ,从而形成另一种弹塑性蠕变损伤理论与数值变分新方法     

An efficient adaptive finite element method algorithm with mass conservation for analysis of liquid face seals简

To improve lubrication effect and seal performance, complicated geometrical hydrodynamic grooves or patterns are often processed on end faces of liquid lubricated mechanical seals. These structures can lead to difficulties in precisely estimating the seal performance. In this study, an efficient adaptive finite element method （FEM） algorithm with mass conservation was presented, in which a streamline upwind/Petrov-Galerkin （SUPG） weighted residual FEM and a fast iteration algorithm were applied to solve the lubrication equations （Reynolds equation）. A mesh adaptation technique was utilized to refine the computation domain based on a residual posterior error estimator. Validation, applicability, and efficiency were verified by comparison among different algorithms and by case studies on seals＇ faces with different groove structures. The study investigated the influence of the order of shape function and the mesh number on the leakage balance. Mesh refinement occurred mainly in cavitation zones when cavitation happened, otherwise it occurred in regions with a high pressure gradient. Numerical experiments verified that the proposed algorithm is a fast, effective, and accurate method to simulate lubrication problems in the engineering field apart from end face seals.     

Performance investigation of 2D lattice Boltzmann simulation of forces on a circular cylinder

The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...     

Numerical investigation on the drag reduction properties of a suction controlled high-rise building

To reduce the wind-induced drag and improve the wind-resistance performance of a high-rise building, steady suction control is introduced into the building structure. Based on validation of the numerical methods by experiment with suction control over the flow separation ofa 3D backward-facing step, the Reynolds stress equation model is used to investigate the drag reduction （DR） properties of a high-rise building whose side faces are controlled by all-height suction. Effects of the orifice geometrical parameters and suction flux parameters on the DR and the separation control are analyzed, and the detailed flow fields are shown to clarify the mechanism of suction control. The results indicate that suction control is very effective in reducing the wind loads on the high-rise building model, and only the dimensionless suction flux dominates. Lastly, the power consumed and the counterforce induced by suction are discussed, the suction models become the ＂zero-drag＂ model under certain suction angles.     

A meshless method based on moving least squares for the simulation of free surface flows

In this paper, a meshless method based on moving least squares (MLS) is presented to simulate free surface flows. It is a Lagrangian particle scheme wherein the fluid domain is discretized by a finite number of particles or pointset; therefore, this meshless technique is also called the finite pointset method (FPM). FPM is a numerical approach to solving the incompressible Navier-Stokes equations by applying the projection method. The spatial derivatives appearing in the governing equations of fluid flow are obtained using MLS approximants. The pressure Poisson equation with Neumann boundary condition is handled by an iterative scheme known as the stabilized bi-conjugate gradient method. Three types of benchmark numerical tests, namely, dam-breaking flows, solitary wave propagation, and liquid sloshing of tanks, are adopted to test the accuracy and performance of the proposed meshless approach. The results show that the FPM based on MLS is able to simulate complex free surface flows more efficiently and accurately.     

Large eddy simulation of unsteady aerodynamic behavior of long-span vaulted roofs

This paper discusses the unsteady aerodynamic behavior of long-span vaulted roofs. First, a forced vibration test in a turbulent boundary layer is conducted in a wind tunnel. The models are force vibrated in the first anti-symmetric mode to investigate the effects of wind speed, rise/span ratio, and the amplitude and frequency of forced vibration on the distributions of wind pressures and unsteady aerodynamic forces. Then, a large eddy simulation (LES) is carried out to clarify the physical mechanism of wind-roof interaction as well as to investigate the influences of a roof’s vibration on the flow field around the roof. From the results of the wind tunnel experiment and the LES, we discuss the characteristics of unsteady aerodynamic forces on a long-span vaulted roof over a wide range of the reduced frequency of vibration. The effect of unsteady aerodynamic forces on the dynamic response of the roof is also discussed. A comparison between the wind tunnel experiment and the LES indicates that the LES can be used effectively to evaluate the unsteady aerodynamic behavior.     

Motion and orientation of cylindrical and cubic particles in pipe flow with high concentration and high particle to pipe size ratio

Lattice Boltzmann method was used to numerically investigate the motion and orientation distribution of cylindrical and cubic particles in pipe flow with high concentration and high particle to pipe size ratio. The transient impulse model of 3D collisions between particles and between particle and wall is proposed. The numerical results are qualitatively in agreement with and quantitatively comparable to the experiment data. The results show that the increases of both the cylindrical particle to pipe size ratio and the particle aspect ratio decrease the rotation about all axes. All rotations of cubic particles decrease with increasing the particle concentration. The cubic particles, rotating more drastically in the flow with large Reynolds number, rotate faster than the cylindrical particles with the same size. The cylindrical particles align with the flow direction more obviously with decreasing Reynolds numbers. However, the orientations of cubic particles are spread all over the range with no significant difference in magnitude, and the Reynolds numbers have no obvious effect on the orientations of cubic particles.     

火电厂耐热钢的弹塑性蠕变损伤本构模型研究

针对火电厂蒸汽管道与汽包一类承压部件常用耐热钢材 ,由于长时间承受高温高压作用 ,处于蠕应变Ⅱ阶段过程中 ,发生弹塑性蠕变损伤问题 ,给出完整的损伤本构描述 ,建立起弹塑性蠕变损伤本构模型及其数值变分原理与有限元离散化形式 ,从而形成完整的数值变分新方法     

二维振荡机翼可压缩绕流IA型反命题变域变分有限元解

In this paper, on the basis of the variational principles developed the finite element method (FEM) is employed for numerical solution of the inverse problem of 2-D unsteady compressible flow around oscillating airfoils by incorporating the non-reflecting far-field boundary conditions and a new unsteady Kutta condition. All unknown boundary (airfoil contour) and discontinuities(shocks and free trailing vortex sheets) are determined via the functional variation with variable domain and artificial density concept. For the numerical realization of the variable-domain variation, a special finite element with self-adjusting nodes is also suggested herein. The numerical results show that the present method is effective for the design of unsteady airfoil.     

Using FEM to predict tree motion in a wind field

In this paper we propose a finite element （FE） simulation method to predict tree motion in a wind field, Two FE tree models were investigated： One model was generated based on a realistic nature-looking geometric tree model, and the other was a symmetric model to investigate the influence of asymmetric material properties on tree motion. The vortex-induced vibration （VIV） theory is introduced to estimate the fluctuating wind force being exerted on tree stems and the fluid-structure interaction （FSI） analysis is also included in the simulation. The results indicate that asymmetric material properties result in the crosswind displacement of the investigated node and the main swaying direction deviation. The simulation reveals that under wind loading, a tree with leaves has much larger swaying amplitude along the wind direction and longer swaying period than a tree without leaves. However, the crosswind swaying amplitude is mainly due to branch interaction. The numerical simulation proved that the interaction of tree branches can prevent dangerous swaying motion developing.     

Variable Domain Variational Finite Element Solution of I A Inverse Problem of 2 D Unsteady Compressible Flow around Oscillating Airfoils

１ＩｎｔｒｏｄｕｃｔｉｏｎＵｐｔｏｎｏｗ,ａｌｔｈｅｐａｐｅｒｓｂｕｔＲｅｆ．［１］ｏｎｔｈｅｉｎｖｅｒｓｅａｅｒｏｄｙｎａｍｉｃｐｒｏｂｌｅｍａｒｅｓｏｌｅｌｙｃｏｎｆｉｎｅｄｔｏｓｔｅａｄｙｆｌｏｗ．ＩｎＲｅｆ．［１］ｔｈｅａｕｔｈｏｒｃａｒｉｅ...     

盲孔法中释放系数的有限元标定研究

基于有限元方法,研究了盲孔法在测量残余应力时释放系数的标定方法。提出了一种新的对释放系数进行标定的有限元模拟办法。运用新的有限元数值模拟方法对ASTM标准中的释放系数值进行了标定。结果表明,有限元标定结果与ASTM标准推荐值有较好的一致性。     

Influence of coherent structures in the gas-particle circular cylinder wake flow

INTRODUCTION Turbulent gas-particle flows are frequently found in natural phenomena and industrial processes.Cases of cylinders in cross flows with particles occur in heat exchange equipment,including the convective zone of a fluidized-bed combustor,and in the primary superheaters,reheaters,and economizers of coal-fired boilers.Coherent structures often occur in the above-mentioned gas-particle flow,and have great effect on such different systems,while some features of tur-bulent multiphas…