Two-level stabilized finite element method for Stokes eigenvalue problem

A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered.This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h = O(H 2),which can still maintain the asymptotically optimal accuracy.It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution,which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h.Hence,the two-level stabilized finite element method can save a large amount of computational time.Moreover,numerical tests confirm the theoretical results of the present method.     

On the stability of multi-layered fluids. Part I: analysis

The linear stability of the Poiseuille flow of multi-layered different fluids, described mathematically by a system of Orr-Somerfeld differential equations, is investigated. A spectral method is used to rewrite this system into a generalized eigenvalue problem, which can be solved with the QZ-algorithm. Special attention is paid to the tractibility of the interfacial conditions of the stability problem. Since we will limit ourselves to a linear stability analysis, the analytical treatment of the interfacial conditions is simplified. Some results related to simple flow configurations are presented. The origin of certain regions of interfacial instability is explained by simple analytical reasoning.     

A new general-purpose least-squares finite element model for steady incompressible low-viscosity laminar flow using isoparametric C1-continuous Hermite elements

This paper deals with a critical evaluation of various finite element models for low-viscosity laminar incompressible flow in geometrically complex domains. These models use Galerkin weighted residuals UVP, continuous penalty, discrete penalty and least-squares procedures. The model evaluations are based on the use of appropriate tensor product Lagrange and simplex quadratic triangular elements and a newly developed isoparametric Hermite element. All of the described models produce very accurate results for horizontal flows. In vertical flow domains, however, two different cases can be recognized. Downward flows, i.e. when the gravitational force is in the direction of the flow, usually do not present any special problem. In contrast, laminar flow of low-viscosity Newtonian fluids where the gravitational force is acting in the direction opposite to the flow presents a difficult case. We show that only by using the least-squares method in conjunction with C¹-continuous Hermite elements can this type of laminar flow be modelled accurately. The problem of smooth isoparametric mapping of C¹ Hermite elements, which is necessary in dealing with geometrically complicated domains, is tackled by means of an auxiliary optimization procedure. We conclude that the least-squares method in combination with isoparmetric Hermite elements offers a new general-purpose modelling technique which can accurately simulate all types of low-viscosity incompressible laminar flow in complex domains.     

Stability and asymmetric vibrations of pressurized compressible hyperelastic cylindrical shells

Cylindrical shells of arbitrary wall thickness subjected to uniform radial tensile or compressive dead-load traction are investigated. The material of the shell is assumed to be homogeneous, isotropic, compressible and hyperelastic. The stability of the finitely deformed state and small, free, radial vibrations about this state are investigated using the theory of small deformations superposed on large elastic deformations. The governing equations are solved numerically using both the multiple shooting method and the finite element method. For the finite element method the commercial program ABAQUS is used.¹ The loss of stability occurs when the motions cease to be periodic. The effects of several geometric and material properties on the stress and the deformation fields are investigated.     

Three‐dimensional simulation of planar contraction viscoelastic flow by penalty finite element method

The planar contraction flow is a benchmark problem for the numerical investigation of viscoelastic flow. The mathematical model of three‐dimensional viscoelastic fluids flow is established and the numerical simulation of its planar contraction flow is conducted by using the penalty finite element method with a differential Phan‐Thien–Tanner constitutive model. The discrete elastic viscous split stress formulation in cooperating with the inconsistent streamline upwind scheme is employed to improve the computation stability. The distributions of velocity and stress obtained by simulation are compared with that of Quinzani's experimental results detected by laser–doppler velocimetry and flow‐induced birefringence technologies. It shows that the numerical results agree well with the experimental results. The numerical methods proposed in the study can be well used to predict complex flow patterns of viscoelastic fluids. Copyright © 2009 John Wiley & Sons, Ltd.     

无穷域势流问题的有限元／差分线法混合求解

提出了一种将有限元和差分线法相结合求解无穷域势流问题的算法。用两同心圆将求解域划分为存在重叠的有限和无限两个区域,在有限和无限域上分别用有限元和差分线法求解Laplace方程边值问题。用差分线法推导出的关系式修正有限元方程,求解该方程组从而得到原问题的解。本算法将求解无穷域问题转化为代数特征值问题和有限域内线性方程组的...     

Comparative study of forced convection of a power-law fluid in a channel with a built-in square cylinder

A detailed comparison between the lattice Boltzmann method and the finite element method is presented for an incompressible steady laminar flow and heat transfer of a power-law fluid past a square cylinder between two parallel plates. Computations are performed for three different blockage ratios (ratios of the square side length to the channel width) and different values of the power-law index n covering both pseudo-plastic fluids (n < 1) and dilatant fluids (n > 1). The methodology is validated against the exact solution. The local and averaged Nusselt numbers are also presented. The results show that the relatively simple lattice Boltzmann method is a good alternative to the finite element method for analyzing non-Newtonian fluids.     

Stability of viscous flow in a rotating porous medium in the form of an annulus: The small-gap problem

The paper deals with the linear stability analysis of laminar flow of a viscous fluid in a rotating porous medium in the form of an annulus bounded by two concentric circular impermeable cylinders. The usual no-slip condition is imposed at both the boundaries. The resulting sixth order boundary value, eigenvalue problem has been solved numerically for the small-gap case by the Runge-Kutta-Gill method, assuming that the marginal state is stationary. The results of computation reveal that the critical Taylor number increases with decreasing permeability of the medium. The problem is found to reduce to the case of ordinary viscous flow in the annulus obtained by Chandrasekhar,¹ when the permeability parameter tends to zero.     

Comparison of spectral and finite element methods applied to the study of the core‐annular flow in an undulating tube

A Galerkin/finite element and a pseudo‐spectral method, in conjunction with the primitive (velocity‐pressure) and streamfunction‐vorticity formulations, are tested for solving the two‐phase flow in a tube, which has a periodically varying, circular cross section. Two immiscible, incompressible, Newtonian fluids are arranged so that one of them is around the axis of the tube (core fluid) and the other one surrounds it (annular fluid). The physical and flow parameters are such that the interface between the two fluids remains continuous and single‐valued. This arrangement is usually referred to as Core‐Annular flow. A non‐orthogonal mapping is used to transform the uneven tube shape and the unknown, time dependent interface to fixed, cylindrical surfaces. With both methods and formulations, steady states are calculated first using the Newton–Raphson method. The most dangerous eigenvalues of the related linear stability problem are calculated using the Arnoldi method, and dynamic simulations are carried out using the implicit Euler method. It is shown that with a smooth tube shape the pseudo‐spectral method exhibits exponential convergence, whereas the finite element method exhibits algebraic convergence, albeit of higher order than expected from the relevant theory. Thus the former method, especially when coupled with the streamfunction‐vorticity formulation, is much more efficient. The finite element method becomes more advantageous when the tube shape contains a cusp, in which case the convergence rate of the pseudo‐spectral method deteriorates exhibiting algebraic convergence with the number of the axial spectral modes, whereas the convergence rate of the finite element method remains unaffected. Copyright © 2002 John Wiley & Sons, Ltd.     

旋转薄壳稳定性的数值分析

用数值方法求解旋转薄壳稳定性问题的主要困难在于:对任意应力状态下的旋转壳进行稳定性分析,各谐波间彼此耦合,致使方程推导的复杂性和计算量都增大,本文采用[3]所提供的单元及经典的能量判据,在最一般的提法下构造了有限元列式,详细讨论了耦合情况,利用本文方案编制程序并进行了许多算例分析,结果令人满意。同时表明在某些情况下边界条件对临界载荷的影响很大,这可以部分地解释古典稳定性分析结果与实验偏差较大以及实验数据分散的原因。     

Newtonian stratified flow through an abrupt expansion

The present paper is devoted to the numerical simulations of stratified generalized Newtonian flow. The main feature of the algorithm is to include the coordinates of the interface as primitive unknowns of the flow problem together with the nodal velocity components and pressures calculated on a deformable finite element mesh. Newton-Raphson's iterative method is used for solving the non-linear problem. Special attention is given to the location of the intersection between the interface and a solid boundary, when it is unknown at the outset. The method is applied to the creeping stratified flow of two Newtonian fluids through a one-to-four abrupt expansion; the results remain valid for the contraction problem. It is found that, under appropriate conditions, the small Newtonian corner vortex can increase by a large factor, both in size and intensity; the degrees of freedom of the simulation are the ratios of flow-rate and of viscosity of the fluids.     

Modeling and simulation of three‐dimensional extrusion swelling of viscoelastic fluids with PTT,Giesekus and FENE‐P constitutive models

The investigation of the extrusion swelling mechanism of viscoelastic fluids has both scientific and industrial interest. However, it has been traditionally difficult to afford theoretical and experimental researches to this problem. The numerical methodology based on the penalty finite element method with a decoupled algorithm is presented in the study to simulate three‐dimensional extrusion swelling of viscoelastic fluids flowing through out of a circular die. The rheological responses of viscoelastic fluids are described by using three kinds of differential constitutive models including the Phan‐Thien Tanner model, the Giesekus model, and the finite extensible nonlinear elastic dumbbell with a Peterlin closure approximation model. A streamface‐streamline method is introduced to adjust the swelling free surface. The calculation stability is improved by using the discrete elastic‐viscous split stress algorithm with the inconsistent streamline‐upwind scheme. The essential flow characteristics of viscoelastic fluids are predicted by using the proposed numerical method, and the mechanism of swelling phenomenon is further discussed.Copyright © 2013 John Wiley & Sons, Ltd.     

Propagation of long waves of finite amplitude at the interface of two viscous fluids

The propagation of long waves of finite amplitude at the interface of two viscous fluids has been studied theoretically. For plane Couette-Poiseuille flow of two superposed layers of fluids of different viscosity, an equation is derived to determine the development in time of the shape of these finite amplitude waves. The influence of the viscosity ratio, the density difference of the fluids and an imposed pressure gradient have been investigated.     

Variable porosity and thermal dispersion effects on vortex instability of a horizontal natural convection flow in a saturated porous medium

A numerical analysis is made to analyze the variable porosity and thermal dispersion effects on the vortex mode of instability of a horizontal natural convection boundary layer flow in a saturated porous medium. The porosity of the medium is assumed to vary exponentially with distance from the wall. In the base flow, the governing equations are solved by using a suitable variable transformation and employing an implicit finite difference Keller Box method. The stability analysis is based on the linear stability theory and the resulting eigenvalue problem is solved by the local similarity approximations. The results indicate that both effects increase the heat transfer rate. In addition, the thermal dispersion effect stabilizes the flow to the vortex mode of disturbance, while the variable porosity effect destabilizes it.

     

A new stable space–time formulation for two‐dimensional and three‐dimensional incompressible viscous flow

A space–time finite element method for the incompressible Navier–Stokes equations in a bounded domain in ?^d (with d=2 or 3) is presented. The method is based on the time‐discontinuous Galerkin method with the use of simplex‐type meshes together with the requirement that the space–time finite element discretization for the velocity and the pressure satisfy the inf–sup stability condition of Brezzi and Babu?ka. The finite element discretization for the pressure consists of piecewise linear functions, while piecewise linear functions enriched with a bubble function are used for the velocity. The stability proof and numerical results for some two‐dimensional problems are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

激光辐照下圆柱壳结构的稳定性分析

针对航天工程中常用的承受轴压作用的薄壁圆柱壳,分别采用解析方法与特征值屈曲有限元方法分析了圆柱壳结构在均匀轴压作用下的稳定性能,得到了屈曲承载力.并进行了对比,验证了有限元模型的合理性.采用线性屈曲特征值有限元方法分析了强激光辐照作用下圆柱壳的稳定性能,分析表明激光辐照导致壳体局部温度上升并由此带来材料参数的改变,是因为激光辐照在壳体中引起了热应力与热应变,使轴压圆柱壳的屈曲承载力明显降低.论文还着重对加筋圆柱壳结构的稳定性进行了研究,数值分析结果表明,加筋能有效提高圆柱壳结构的抗压承载力,激光辐照作用下,加筋对轴压作用下圆柱壳的屈曲承载力的提高作用更为明显.     

One dimensional finite element solutions for a Maxwell fluid

The one-dimensional flow of a Maxwell fluid is used to study several aspects of the numerical simulation of viscoelastic fluids. The numerical results are produced via a Galerkin finite element method of the mixed type and are compared with an analytic solution. The accuracy of the numerical solutions are evaluated with respect to variations in mesh refinement and changes in the interpolation order for the stress variables. Various convergence and accuracy problems associated with the chosen numerical schemes are noted and a tentative explanation for these difficulties is provided through use of an eigenvalue analysis.     

Investigation of the stability for inviscid compressible swirling flow between concentric cylinders

The temporal stability on inviscid compressible swirling flow between two concentric cylinders is investigated. First, a linearized differential equation is derived. Two stability criteria are derived for compressible swirling flow by an analytic method analogous to Ludwieg ’s method. A finite-difference numerical method is then used to solve the eigenvalue problem of this differential equation, to get temporal growth rate and to check these stabilitv criteria derived. Finally.The effect of compressibility for stability is disscused.