黏弹流体流动的数值模拟研究进展

综述了黏弹流体流动数值模拟的研究进展,突出介绍近十年来有限元法在黏弹流体流动数值模拟研究中取得的成果,通过动量方程的适当变形和本构方程离散权函数的合理选择,可以显著增强数值计算的稳定性。得到较高Weissenberg数下的解,同时文中对黏弹流体流动数值模拟中本构方程的应用、非等温情况和三维空间下的研究进行了介绍。     

Bingham流体在旋转圆盘上流动的数值解

导出了描述Ｂｉｎｇｈａｍ流体在旋转圆盘上流动的基本方程，用差分方法数值解薄膜厚度分布方程，得到二种类型的厚度分布。数值解分别和计算机磁盘的厚度分布，Ｊｅｎｅｋｈｅ等的实验结果定性一致。     

自由液面大晃动的流固耦合数值分析方法研究进展

本文较为全面地评述了国内外对具有自由液面大晃动的流固耦合问题的研究。在第2节,以两种典型的工程问题（石油化工储液罐抗震性和提高机理问题和液态金属快中子增殖堆主容器流固耦合问题）为背景全面评述了这类流固耦合问题的有关研究进展;文章的第3节较为深入讨论了这类流固耦合问题中所采用的分析方法（包括解析方法、半解析方法和数值方法以及问题描述的ALE格式等）和分析模型（包括位移－位移模型和多种的位移－势模型等）;第4节首先从问题求解方面,说明了具有自由液面大晃动的流固耦合问题的性质和特点,然后讨论了交替求解方法及其对求解此类问题的特别适应性;最后,本文给出了这方面今后应进一步研究的问题。     

贮腔类三维自由液面动力学问题数值研究

讨论了贮腔类三维自由液面动力学问题的数值研究,将任意的拉格朗日-欧拉运动学描述关系引入到系统的控制方程中,采用任意的拉格朗日-欧拉描述跟踪自由液面,推导了自由面上结点的法向矢量计算公式。采用Galerkin余量法推导了Navier-Stokes方程的空间离散有限元方程,采用三维自由液面上微分几何理论推导了表面张力计算公式。数值研究中考虑了接触角效应,最后进行了三维数值算例分析。     

纤维悬浮液搅拌流动的数值模拟

由于缺乏适当的本构方程，对纤维悬浮液流动的研究一直局限于纤维的牛顿流体悬浮液。本文采用ＭＵＣＭ模型对作者最近提出的纤维Ｏｌｄｒｏｙｄ－Ｂ流体悬浮液的本构方程作了改进，并对锚式桨搅拌槽的二维Ｏｌｄｒｏｙｄ－Ｂ流体和牛顿流体纤维悬浮液搅拌流动作了数值模拟。模拟的结果表明，本文所用的模型和方法能有效地抑制过大局部应力的影响并合理地处理流体的记忆效应。     

黏性不可压缩流体流动前沿的数值模拟

提出了模拟注射成型中黏性、不可压缩流体流动前沿的新方法. 将Hele-Shaw流动应用于非 等温条件下的黏性、不可压缩流体，建立了流动分析模型，用充填因子的输运方程描述流动 前沿. 应用高阶Taylor展开式计算每一时间步长的充填因子，用Galerkin方法导出了计算 充填因子各阶导数的递推公式. 给出了时间增量的选取方法，证明了它的稳定性. 针对Han 设计的试验模具，用相同的材料及工艺条件模拟充填过程，比较了传统方法和该方法的模 拟结果与实验结果的差异. 算例分析表明，该方法可以有效地提高注射成型中流动前沿的 模拟精度和计算效率.     

离心泵叶轮内宾汉流体湍流流场的数值计算

考虑宾汉流体本构关系特点,建立了任意曲线坐标系（ζ,η）下宾汉流体湍流流动的基本方程,应用压力加权校正算法,实现了速度场和压力场的关联,采用交错网格技术,解决了非物理压力振荡问题,在此基础上,对离心泵叶轮内回转面上宾汉流体湍流流动进行了数值模拟,并分析探讨了离心泵叶轮内宾汉流体湍流流动机理。     

粘弹性流体的入口收敛流动

本文以聚合物流体为研究对象,对其在入口收敛流动中产生的粘弹效应及机理进行了初步的讨论和分析,并就近10年来国内外有关粘弹性流体入口流动研究及进展作了简要的评述。      

粘弹流体轴对称流动的格子Boltzmann方法

基于BGK碰撞模型，通过在迁移方程中引入作用力项，建立了粘弹流体的轴对称格子Boltzmann模型．通过Chapman-Enskog展开，获得了准确的柱坐标下轴对称宏观流动方程．采用双分布函数对运动方程和本构方程进行迭代求解，模拟分析了粘弹流体管道流动，获得了流场中的速度和构型张量的分布，通过与解析解进行比较，验证了模型的准确性．研究了作为粘弹流体流动基准问题的收敛流动，对涡旋位置进行了定量分析，将回转长度的计算结果与有限体积法进行了比较，两种数值结果十分吻合．研究结果表明，模型能够准确表征粘弹流体的轴对称流动，具有较广阔的应用前景．     

幂率流体圆管非定常流动的数值分析

本文讨论水平圆管中幂率流体的起动问题。用显式和隐式两种格式,我们得到了问题的数值解并从而得到流动建立时间的近似公式。     

Numerical approach for generic three-phase flow based on cut-cell and ghost fluid methods

In this paper, we introduce numerical methods that can simulate complex multiphase flows. The finite volume method, applying Cartesian cut-cell is used in the computational domain, containing fluid and solid, to conserve mass and momentum. With this method, flows in and around any geometry can be simulated without complex and time consuming meshing. For the fluid region, which involves liquid and gas, the ghost fluid method is employed to handle the stiffness of the interface discontinuity problem. The interaction between each phase is treated simply by wall function models or jump conditions of pressure, velocity and shear stress at the interface. The sharp interface method “coupled level set (LS) and volume of fluid (VOF)” is used to represent the interface between the two fluid phases. This approach will combine some advantages of both interface tracking/capturing methods, such as the excellent mass conservation from the VOF method and good accuracy of interface normal computation from the LS function. The first coupled LS and VOF will be generated to reconstruct the interface between solid and the other materials. The second will represent the interface between liquid and gas.     

Numerical methods for a fluid mixture model

In this paper, we firstly apply generalized difference methods to solve a fluid mixture model. The model is usually used to describe the tissue deformations and contains a nonlinear hyperbolic equation and an elliptic equation. Most people have used finite difference methods for solving the elliptic equation and other schemes for solving the hyperbolic equation. It is well known that the accuracy of traditional finite difference method is not high. This may be a serious disadvantage in the fluid mixture model, which describes cell movements and tissue deformations. The numerical methods we propose to improve accuracy are based on generalized Galerkin methods and dual decomposition. By choosing suitable trial function space and test function space, our generalized upwind difference schemes exhibit second‐order convergence in space for smooth problems and can eliminate numerical oscillations for discontinuous problems. Some numerical results are presented to demonstrate the advantages of our methods. Copyright © 2012 John Wiley & Sons, Ltd.     

Asymptotic analysis of the finite cone-and-plate flow of a non-newtonian fluid

Consider the shearing flow of a viscoelastic fluid trapped by surface tension between a cone and a plate. An asymptotic analysis of this problem in the limit of small gap angle has been done. This limit is realized in many practical situations. It is assumed that the Deborah number De, the Reynolds number Re, and the retardation parameter β are all order unity and that the shape of the free surface is very nearly spherical. Closed form analytic expressions are obtained for the leading terms of the primary and weak secondary motion of the fluid as well as the meniscus shape. It is found that the velocity field is bounded and continuous if and only if . There is a family of curves in the De-β plane on which the velocity field has a removable singularity at the origin. The secondary flow is made up of either one or two toroidal vortices. The meniscus has a bulge near the rotating cone and a trough near the stationary plate.     

Numerical analysis of funnel formation during unsteady fluid outflow from a rotating cylindrical vessel

The main features of the unsteady outflow of a fluid from a cylindrical vessel rotating together with it at a constant angular velocity through a central drain hole in the bottom are studied. The software package STAR-CD tested on the results of experiments with water in the absence of rotation is used. Certain important features of the phenomenon under consideration related with the formation of vortex funnels in the fluid are established. The effect of the main control parameters of the problem, namely, angular velocity, viscosity, initial depth, etc., is analyzed.     

Possible classification of the self-oscillatory spouting regimes of plane vertical submerged jets of a heavy fluid

New results of an experimental investigation of self-oscillatory regimes of plane vertical jet spouting from beneath the free surface of a heavy incompressible fluid are discussed. The experiments were performed on a setup with discharge over a weir. The range of dimensionless jet submergence values on which bifurcation change of spouting regime is observable is studied. It is established that on the Froude number and dimensionless jet submergence ranges considered in the study six characteristic spouting regimes differing in free surface shape and self-oscillation frequency can exist. It is shown that these regimes can be subdivided into three typical groups with respect to the dependence of the self-oscillation period on the jet flow rate. A dimensionless parameter that makes it possible to identify the boundaries of the bifurcation change in spouting regimes is obtained for each of these groups. For certain spouting regimes without the formation of free jets numerical calculations are carried out using the STAR-CD software package; the calculated results are in good agreement with experimental data.     

Assessment of volume of fluid and immersed boundary methods for droplet computations

The volume of fluid (VOF) and immersed boundary (IB) methods are two popular computational techniques for multi‐fluid dynamics. To help shed light on the performance of both techniques, we present accuracy assessment, which includes interfacial geometry, detailed and global fluid flow characteristics, and computational robustness. The investigation includes the simulations of a droplet under static equilibrium as a limiting test case and a droplet rising due to gravity for Re?1000. Surface tension force models are key issues in both VOF and IB and alternative treatments are examined resulting in improved solution accuracy. A refined curvature model for VOF is also presented. With the newly developed interfacial treatments incorporated, both IB and VOF perform comparably well for the droplet dynamics under different flow parameters and fluid properties. Copyright © 2004 John Wiley & Sons, Ltd.     

An orthogonal mapping technique for the computation of a viscous free-surface flow

In this paper we describe finite element computations of the free-surface flow of a viscous fluid down an undulating inclined plane. The technique developed here employs an orthogonal mapping that is computed along with the velocity and pressure. This is allied to a technique to compute symbolically the Jacobian and other derivatives required for numerical continuation methods. The solutions obtained are compared with laboratory experiments and finite element computations reported by Pritchard and co-workers. The finite element computational method used by these authors employs spines to represent the free surface. An excellent agreement is shown to exist between the new computations and the laboratory experiments, and with the numerical solutions of Pritchard and co-workers.     

Formation of the free surface of a viscous fluid volume inside a rotating horizontal cylinder

Two-dimensional viscous flow with a free surface in a horizontal cylinder rotating at a constant speed is investigated numerically using the boundary element method. It is shown that in the initial stage of rotation of the cylinder four different variants of the behavior of the free surface can be realized in the stage of transition from horizontal to steady-state form.