各向异性介质中弹性波的数值模拟

提出了一种非均匀各向异性介质中弹性波传播的数值模拟算法。该方法可以灵活地运用于具有任意地表形状、内部孔洞、固液边界和不规则内部交界面的介质情况,另外,该方法自然满足复杂几何边界的自由表面条件。这种基于三角形和四边形离散网格的算法使用的是围绕每个节点的积分平衡方程,而不是其它有限差分法中使用的各个节点满足的弹性动力学的微分方程。该文工作是非均匀各向同性介质中弹性波传播格子法研究的继续。除了研究各向异性介质中波的传播以外,还给出了一种能够省时的四边形网格的格子法。     

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity field by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.     

NUMERICAL SIMULATION OF ELASTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA

A new numerical simulation algorithm is presented for the elastic wave propagation in heterogeneous anisotropic media. We make discretization of the computational domain by using triangular and quadrangular grids. The scheme is based on integral equilibrium at each node to simulate the elastic wave propagation in heterogeneous anisotropic media. The scheme is very flexible in dealing with arbitrary surface topography, inner openings, liquid-solid boundaries and irregular interfaces. Moreover, the free-surface condition of complex geometrical boundaries can be satisfied naturally. This work is an extension of the grid method for the elastic wave propagation in heterogeneous isotropic media, and a quadrangular grid with low computational cost is also introduced. Project supported by the National Natural Science Foundation of China(19672016).     

Die-swell,splashing drop and a numerical technique for solving the Oldroyd B model for axisymmetric free surface flows

This work deals with the development of a numerical method for simulating viscoelastic axisymmetric free surface flow of an Oldroyd B fluid. A novel formulation is developed for the computation of the non-Newtonian extra-stress components on rigid boundaries and on the symmetry axis. The full free surface stress conditions are employed. The resulting governing equations are solved by finite differences on a Marker-and-cell (MAC) type grid. Validation is provided by simulating a pipe flow problem. The classical die-swell problem is solved and swelling ratios are provided. The height of the splash caused by a falling liquid drop for various Reynolds and Weissenberg numbers is then studied, and the height of the splash is shown to diminish with increasing viscoelasticity.     

离子选择性表面电对流的多块LB模拟

离子选择性表面(如纳米通道、离子交换膜等)复杂的动力学现象为微纳流控技术的发展提供了新思路. 向带有离子选择性表面的电解质溶液施加电压, 通过液体的电流密度会经历复杂的非线性变化过程; 当电压超过某一临界值时会引发对流现象, 这种流动被称为第二类电渗或离子选择性表面的电对流, 关于此类问题的数值求解引发了大量的研究. 本文提出一种基于多块网格加密的格子玻尔兹曼方法(lattice Boltzmann method, LBM)的数值模型, 用于模拟第二类电渗流动. 结合该算法, 给出了求解流动、电势和离子浓度的网格信息交换方程, 较好地解决了此类问题中大浓度梯度边界对计算分辨率的要求. 利用该数值模型模拟获得的电流?电压特性曲线先随着电压升高而迅速增大, 随后达到饱和状态, 与理论解吻合良好. 此外, 模拟结果还表明, 当流动发生后, 相对低电压下的流动倾向于形成大涡且流动呈指数趋势增强; 而较大电压则会先激发多个小涡, 并逐渐合并为大涡流动, 且大涡流动有更高的离子输运效率. 此外, 除了模拟离子选择性表面的电对流现象, 本文提出的数值格式还可拓展到其他电流体动力学问题的模拟.      

基于2D网格的轴对称浸没边界法

浸没边界法是处理颗粒两相流中运动边界问题的一种常用数值模拟方法. 当研究的物理问题的无量纲参数满足一定要求时, 该流场结构呈现轴对称状态. 为此本文提出了一种基于2D笛卡尔网格和柱坐标系的轴对称浸没边界法. 该算法采用有限体积法(FVM)对动量方程进行空间离散, 并通过阶梯状锐利界面替代真实的固体浸没边界来封闭控制方程. 为了提高计算效率, 本文采用自适应网格加密技术提高浸没边界附近网格分辨率. 由于柱坐标系的使用, 使得动量方程中的黏性项产生多余的源项, 我们对其作隐式处理. 此外, 在对小球匀速近壁运动进行直接数值模拟时, 由于球壁间隙很小, 间隙内的压力变化比较剧烈. 因此想要精确地解析流场需要很高的网格分辨率. 此时, 需要在一个时间步内多次实施投影步来保证计算的稳定性. 而在小球自由碰壁运动中, 我们通过引入一个润滑力模型使得低网格分辨率下也能模拟小球近壁处的运动. 最后通过小球和圆盘绕流、Stokes流小球近壁运动以及小球自由下落碰壁弹跳算例验证本算法对于轴对称流的静边界和动边界问题均是适用和准确的.      

Axisymmetric shells and plates on tensionless elastic foundations

This paper first describes a finite element method for the large deflection analysis of axisymmetric shells and plates on a nonlinear tensionless elastic foundation. Through the use of discrete data points, any form of nonlinear elastic foundation behaviour can be easily modelled. The analysis is then validated by comparison with existing results for circular plates and beams as the only existing results for shells on tensionless foundations are found to be in error. Following this verification, the analysis is applied to investigate the behaviour of shallow spherical shells subject to a central concentrated load on tensionless linear elastic foundations. A number of insightful conclusions regarding the behaviour of such structure-foundation systems are drawn. The numerical results for shells are believed to be the first correct results, which may be useful in benchmarking results from other sources in the future.     

Determination of the Axisymmetric Geometrically Nonlinear Thermoviscoelastoplastic State of Laminated Medium-Thickness Shells

A method is developed to determine the axisymmetric geometrically nonlinear thermoelastoviscoplastic stress–strain state of branched laminated medium-thickness shells of revolution. The method is based on the hypotheses of a rectilinear element for the whole set of layers. The shells are subject to loads that cause a meridional stress state and torsion. They can consist of isotropic layers, which deform beyond the elastic limit, and elastic orthotropic layers. The relations of thermoviscoplastic theory, which describe simple processes of loading, are employed as the equations of state for the isotropic layers. The solution of the problem is reduced to numerical integration of systems of differential equations. The geometrically nonlinear elastoplastic state of a two-layer corrugated shell of medium thickness is calculated as an example     

Flow of elasto-viscoplastic liquids through an axisymmetric expansion–contraction

In this work, we used a recently proposed constitutive model that addresses the elastic behavior of viscoplastic liquids. We used the equation in an expansion–contraction axisymmetric flow, and compared the results with experimental data found in the literature. Steady, inertialess numerical solutions were obtained by solving the balance equations of mass and momentum via the finite element method. We obtained the velocity and stress fields for different combinations of the governing parameters, and observed that elasticity leads to significant changes in the shape and position of the yield surface, affecting both the extra-stress and the rate of deformation fields in the cavity. The trends observed are in qualitative agreement with visualization results available in the literature.     

The effect of an insoluble surfactant on the skin friction of a bubble

In this paper, we develop a novel moving mesh method suitable for solving axisymmetric free-boundary problems, including the Marangoni effect induced by surfactant or temperature variation. This method employs a body-fitted grid system where the gas–liquid interface is one line of the grid system. We model the surfactant equation of state with a non-linear Langmuir law, and, for simplicity, we limit ourselves to the situation of an insoluble surfactant. We solve complicated dynamic boundary conditions accurately on the gas–liquid interface in the framework of finite-volume methods. Our method is used to study the effect of a surfactant on the skin friction of a bubble in a uniaxial flow. For the limiting case where the surface diffusivity is zero, the effect of a tangential stress generated by the surface tension gradient, allows us to explain a new phenomenon in high concentration regimes: larger surface tension, but also larger deformation. Furthermore, this condition leads to the formation of boundary layers and flow separation at high Reynolds numbers. The influence of these complex flow patterns is examined.     

集中载荷作用下具有光滑中心波纹膜片的非线性分析

波纹膜片是一种薄壳弹性体，由于它的参数很多，又相互制约，所以使得它的设计很复杂。在大多数位移式仪器仪表中，要求波纹膜片产生至少和膜片厚度是同样数量级的弹性位移。这就要求使用薄壳的几何非线性理论进行分析。大多数学者研究波纹膜片的弯曲问题，是采用扁壳理论讨论具有浅波纹的膜片。而工程实际中，经常遇到深波纹膜片，这就要求从一般壳体大挠度方程进行求解。本文采用轴对称旋转壳体的简化Re-issner方程。研究了在中心集中载荷作用下具有光滑中心波纹膜片的非线性弯曲问题。应用积分方程方法，可以获得膜片的特征关系（载荷-中心挠度关系）和应力分布。文末给出实例计算的数值结果。     

三维非均匀介质中弹性波传播的数值模拟

提出了一种三维非均匀介质中弹性波传播数值模拟的方法,文中称为三维格子法。该算法是二维格子法（一种二维非均匀介质中P－SV波传播的数值模拟算法）向三维非均匀介质情况的推广。在空间离散上该文方法与有限元方法类似,容许根据连续体的形状和介质分界面任意剖面网格,且自然满足自由表面边界条件。不同于常规有限差分法在各个节点上满足动力学微分方程,该算法通过满足各节点周围格子的整体平衡（积分平衡方程）来对问题进行求解,三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当。算例表明,该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射和绕射。     

Implementation of high‐order compact finite‐difference method to parabolized Navier–Stokes schemes

The numerical solution to the parabolized Navier–Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth‐order compact finite‐difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite‐difference algorithm of Beam and Warming type with a high‐order compact accuracy. A shock‐fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously computed with the fourth‐order accuracy. The computations are carried out for a benchmark case: hypersonic axisymmetric flow over a blunt cone at Mach 8. A sensitivity study is performed for the basic flowfield, including profiles and their derivatives obtained from the fourth‐order compact PNS and IPNS solutions, and the effects of grid size and numerical dissipation term used are discussed. The present results for the flowfield variables and also their derivatives are compared with those of other basic flow models to demonstrate the accuracy and efficiency of the proposed method. The present work represents the first known application of a high‐order compact finite‐difference method to the PNS schemes, which are computationally more efficient than Navier–Stokes solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

Computation of complex fluid flows using an adaptive grid method

Recently the concept of adaptive grid computation has received much attention in the computational fluid dynamics research community. This paper continues the previous efforts of multiple one-dimensional procedures in developing and asessing the ideas of adaptive grid computation. The focus points here are the issue of numerical stability induced by the grid distribution and the accuracy comparison with previously reported work. Two two-dimensional problems with complicated characteristics—namely, flow in a channel with a sudden expansion and natural convection in an enclosed square cavity—are used to demonstrate some salient features of the adaptive grid method. For the channel flow, by appropriate distribution of the grid points the numerical algorithm can more effectively dampen out the instabilities, especially those related to artificial boundary treatments, and hence can converge to a steady-state solution more rapidly. For a more accurate finite difference operator, which contains less undesirable numerical diffusion, the present adaptive grid method can yield a steady-state and convergent solution, while uniform grids produce non-convergent and numerically oscillating solutions. Furthermore, the grid distribution resulting from the adaptive procedure is very responsive to the different characteristics of laminar and turbulent flows. For the problem of natural convection, a combination of a multiple one-dimensional adaptive procedure and a variational formulation is found very useful. Comparisons of the solutions on uniform and adaptive grids with the reported benchmark calculations demonstrate the important role that the adaptive grid computation can play in resolving complicated flow characteristics.     

Numerical calculations for two-phase jet flow of a mixture of pulverized-coal and gas

In this paper, Euler-Lagrange type equations are used to describe the jet flow of a mixture of pulverized-coal and gas, which is an unsteady axisymmetric two-phase flow. By means of the finite-difference method, the coal particle's distribution, velocity and trajectory in the flow field are obtained. The coal particles are represented by a finite number of computational particles. Each particle's diameter is randomly assigned according to a given distribution. The states of the computational particles are different from each other. Turbulence is accounted for in a stochastic model. Explicit time-splitting scheme is used to calculate the strongly coupling interphase term. The numerical results are reasonable. The comparison between the numerical results and the experiment data for the case of the oil droplet injection shows good agreement. This numerical technique can be extended to the calculation of other two-phase flows of dilute particles or a droplet system. Mr. Mei Renwei also participated in the work of this paper.     

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

A staggered control volume scheme for unstructured triangular grids

The purpose of this work is to introduce and validate a new staggered control volume method for the simulation of 2D/axisymmetric incompressible flows. The present study introduces a numerical procedure for solving the Navier–Stokes equations using the primitive variable formulation. The proposed method is an extension of the staggered grid methodology to unstructured triangular meshes for a control volume approach which features ease of handling of irregularly shaped domains. Two alternative elements are studied: transported scalars are stored either at the sides of an element or at its vertices, while the pressure is always stored at the centre of an element. Two interpolation functions were investigated for the integration of the momentum equations: a skewed mass-weighted upwind function and a flow-oriented exponential shape function. The momentum equations are solved over the covolume of a side or of a vertex and the pressure–velocity coupling makes use of a localized linear reconstruction of the discontinuous pressure field surrounding an element in order to obtain the pressure gradient terms. The pressure equation is obtained through a discretization of the continuity equation which uses the triangular element itself as the control volume. The method is applied to the simulation of the following test cases: backward-facing step flow, flow over a two-dimensional obstacle and flow in a pipe with sudden contraction of cross-sectional area. All numerical investigations are compared with experimental data from the literature. A grid convergence and error analysis study is also carried out for flow in a driven cavity. Results compared favourably with experimental data and so the new control volume scheme is deemed well suited for the prediction of incompressible flows in complex geometries. © 1997 John Wiley & Sons, Ltd.     

A new coupled level set and volume-of-fluid method to capture free surface on an overset grid system

We present a new coupled level set and volume-of-fluid (CLSVOF) method for free surface flow simulations on an overset grid system. The coupled method takes advantages of the strengths of the level set (LS) method and the volume-of-fluid (VOF) method, and is superior to either single method. The novelty of the present method lies in that we develop the methodology for an overset grid system of embedding, overlapping and moving structured grids. The new methodology accurately captures interface and greatly preserves mass on an overset grid system by demonstrating the 3D sphere advection test. The method is coupled to a well validated Reynolds-Averaged Navier–Stokes incompressible flow solver. The method is validated with the dam-breaking flow interacting with a 3D obstacle (square structure/circular cylinder) by comparing the numerical results with available experimental and numerical studies. The water impact of a sphere case is further performed to demonstrate the capabilities of the new method on a complicated moving overset grid system.     

An extended finite element method for the simulation of particulate viscoelastic flows

We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.     

Nonlinear analysis of the deformation and breakup of viscous microjets using the method of lines

We present a numerical model for predicting the instability and breakup of viscous microjets of Newtonian fluid. We adopt a one‐dimensional slender‐jet approximation and obtain the equations of motion in the form of a pair of coupled nonlinear partial differential equations (PDEs). We solve these equations using the method of lines, wherein the PDEs are transformed to a system of ordinary differential equations for the nodal values of the jet variables on a uniform staggered grid. We use the model to predict the instability and satellite formation in infinite microthreads of fluid and continuous microjets that emanate from an orifice. For the microthread analysis, we take into account arbitrary initial perturbations of the free‐surface and jet velocity, as well as Marangoni instability that is due to an arbitrary variation in the surface tension. For the continuous nozzle‐driven jet analysis, we take into account arbitrary time‐dependent perturbations of the free‐surface, velocity and/or surface tension as boundary conditions at the nozzle orifice. We validate the model using established computational data, as well as axisymmetric, volume of fluid (VOF) computational fluid dynamic (CFD) simulations. The key advantages of the model are its ease of implementation and speed of computation, which is several orders of magnitude faster than the VOF CFD simulations. The model enables rapid parametric analysis of jet breakup and satellite formation as a function of jet dimensions, modulation parameters, and fluid rheology. Copyright © 2010 John Wiley & Sons, Ltd.