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

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

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

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

幂律流体圆柱绕流的格子波尔兹曼模拟

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.      

广义有限差分法求解Kirchhoff和Winkler薄板弯曲问题

本文将广义有限差分法用于数值计算Kirchhoff板和Winkler板的弯曲问题.广义有限差分法是基于最小二乘原理的一种区域型无网格方法. 相比于传统的网格类数值解法,广义有限差分法无需网格生成且无需数值积分.通过数值实验结果表明,广义有限差分法可以有效地求解两类薄板在不同横向荷载作用下的弯曲问题.     

幂律流体圆柱绕流的格子波尔兹曼模拟简

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明：插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.     

Burgers方程的小波精细积分算法

求解偏微分方程的常用方法包括有限差分法、有限元法等。近年来,小波分析在偏微分方程数值求解中的应用已引起很多学者的关注,例如采用Daubechies小波或shannon小波构造的小波配置方法已经取得较好的结果。钟万勰院士提出的偏微分方程的子域精细积分方法是一种半解析方法,方法简单,精度高。将小波方法和精细积分方法相结合应用于偏微分方程的数值求解中将有利于提高算法的精度和稳定性,为此本文以Burgers方程为例,提出了一种求解一维非线性抛物型偏微分方程的小波精积分方法。该方法用拟小波配点法对空间域进行离散,建立起对时间的常微分方程组,然后采用精细时程积分方法对该方程组求解。数值计算结果表明,该方法同其它方法相比,具有计算格式简单,数值稳定性和精度较高的优点。     

等截面梁有限变形的传递函数增量算法

本文介绍了一种计算等截面梁有限变形的新方法－传递函数增量算法，它是一种半解析数值计算方法。此算法充分利用增量失空法Ｇａｕｓｓ求积公式计算非线性有限变形的特点，并将这些特点与传递函数方法，有效地结合起来，既避免了数值方法计算量大的困难，又使得求解高阶非线性微分方程的解析解成为可能，算例分析表明，这是一种易编程，计算量小，收敛快，求解精度高的行之有效的计算方法。     

二阶双曲型方程的精细时程积分法

对于二阶双曲型偏微分方程初边值问题，可以用有限差分法进行求解。通常的有限差分法在使用过程中受到精确度和稳定性的限制，本文提出求解二阶双曲型方程的精细时程积分法。由于这种方法是半解析方法，在时间域上可以精确计算，所以这种方法不仅精确度高，而且还绝对稳定。文末的数值算例进一步验证了上述结构，而且对大的时间步长（例如△t=0.5）仍然获得精度很高的数值结果。可见，精细时程积分法是一种很实用的方法。     

含肿瘤皮肤组织传热分析的广义有限差分法

生物传热分析在低温外科手术、肿瘤热疗、病热诊断等临床医学治疗和诊断中有着广泛的应用. 由于健康皮肤组织内肿瘤的存在使得肿瘤附近区域的温度会明显升高, 这一特性常被用于检测皮肤组织内的肿瘤生长, 因此有必要开展生物传热数值分析的研究. 本文以含肿瘤的皮肤组织为研究对象, 将一种新型区域型无网格配点法——广义有限差分法应用于能描述含肿瘤皮肤组织传热过程的Pennes方程求解. 广义有限差分法利用泰勒展开式与移动最小二乘法将计算区域内的每个离散点上的物理量导数表示成其与邻近点物理量及权重系数的线性组合, 进而构建得到仅含各离散点未知物理量的线性方程组. 该方法不仅具有无需划分网格、避免数值积分等无网格配点法的优点, 同时还克服了大多数无网格配点法中插值矩阵高度病态的问题, 为此类方法在大规模工程数值计算中的应用提供了可能性. 本文首先介绍了模拟含肿瘤皮肤组织传热过程的广义有限差分法离散模型, 随后通过不含肿瘤与含规则形状肿瘤的基准算例, 检验广义有限差分法的计算精度与收敛性; 在此基础上, 通过数值模拟研究不同肿瘤形状及肿瘤位置分布对皮肤组织内温度分布的影响.      

一种低数值耗散的求解油水驱替过程的数值方法

油水驱替问题的压力场和饱和度场之间是存在粘性耦合的.在求解油水驱替问题时,常用的有限差分数值方法不仅会引入大量数值耗散,从而导致求解精度降低,而且可能会产生网格取向效应.为了减小数值耗散的影响,本文将能够高效率求解压力场的有限分析法和零数值耗散地求解饱和度场的值域离散网格法相结合,提出了一种低数值耗散的顺序求解方法.与解析解的对比表明,该方法能够保持间断的锐利形态,数值耗散得到有效抑制.利用该方法与有限差分法计算结果的对比,说明了流度比不利时,网格取向效应是在数值耗散作用下的流动不稳定性在有限差分法计算结果中的表现形式.     

不同壁面绕流特性的格子Boltzmann模拟研究

为了探讨不同壁面的绕流特性,针对粘性流场中,不同壁面诱导的涡脱落现象以及升阻力系数等流场特性进行了格子Boltzmann数值研究。利用基于分子动理论的格子Boltzmann方法(LBM)求解Navier-Stokes方程,实现对流体运动的描述,针对不同的壁面条件,分别采用不同的格子Boltzmann流-固壁面处理方法。采用Half-way反弹边界条件来处理平直壁面,而曲壁面则采用LBM与有限差分法相结合的形式进行处理,计入了壁面与标准网格不重合对结果造成的影响。开发相应的计算程序,计算结果与已发表文献结果吻合良好,验证了数值模型的正确性。同时,探讨了进出口边界与钝体中心的距离对结果的影响。对比分析了不同壁面的绕流模型中升阻力系数、斯托罗哈数和涡量云图等,并进一步研究了雷诺数条件的影响。结果表明,不同壁面的绕流特性具有明显差异,且同时受雷诺数的显著影响;一般地,平直壁面对于来流作出的响应更迅速。     

A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

Lattice Boltzmann method for simulating particle–fluid interactions

The lattice Boltzmann method （LBM） has gained increasing popularity in the last two decades as an alternative numerical approach for solving fluid flow problems. One of the most active research areas in the LBM is its application in particle-fluid systems, where the advantage of the LBM in efficiency and parallel scalability has made it superior to many other direct numerical simulation （DNS） techniques. This article intends to provide a brief review of the application of the LBM in particle-fluid systems. The numerical techniques in the LBM pertaining to simulations of particles are discussed, with emphasis on the advanced treatment for boundary conditions on the particle-fluid interface. Other numerical issues, such as the effect of the internal fluid, are also briefly described. Additionally, recent efforts in using the LBM to obtain closures for particle-fluid drag force are also reviewed.     

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid Past a Vertical Plate Immersed in a Porous Medium

The effects of viscous dissipation on unsteady free convection from an isothermal vertical flat plate in a fluid saturated porous medium are examined numerically. The Darcy–Brinkman–Forchheimer model is employed to describe the flow field. A new model of viscous dissipation is used for the Darcy–Brinkman–Forchheimer model of porous media. The simultaneous development of the momentum and thermal boundary layers are obtained by using a finite difference method. Boundary layer and Boussinesq approximation have been incorporated. Numerical calculations are carried out for various parameters entering into the problem. Velocity and temperature profiles as well as local friction factor and local Nusselt number are shown graphically. It is found that as time approaches infinity, the values of friction factor and heat transfer coefficient approach steady state.     

基于显式几何更新算法的地下管道自动形状优化

基于计算力学中的结构优化思想，应用一种新型的显式几何更新算法，自行编制C++程序，实现地下管道形状设计的自动优化。管道内的流体假设为牛顿不可压缩流，并考虑惯性项。优化区域主要为管道竖直方向和水平方向的过渡段。形状优化的设计变量是几何边界的有限元节点坐标，优化目标是实现流体黏性能耗散的最小化。优化过程基于形状梯度，即通过形状敏感度分析来求解目标函数相对于设计变量的偏导数。所使用的显式几何更新算法既可以通过网格清晰描述形状，也可以大范围地自动更新网格。详细介绍了地下管道自动形状优化过程的关键步骤。通过数值算例探讨了不同注入速度、密度和黏度对其最优形状的影响。     

Velocity–pressure coupling in finite difference formulations for the Navier–Stokes equations

A new numerical procedure for solving the two‐dimensional, steady, incompressible, viscous flow equations on a staggered Cartesian grid is presented in this paper. The proposed methodology is finite difference based, but essentially takes advantage of the best features of two well‐established numerical formulations, the finite difference and finite volume methods. Some weaknesses of the finite difference approach are removed by exploiting the strengths of the finite volume method. In particular, the issue of velocity–pressure coupling is dealt with in the proposed finite difference formulation by developing a pressure correction equation using the SIMPLE approach commonly used in finite volume formulations. However, since this is purely a finite difference formulation, numerical approximation of fluxes is not required. Results presented in this paper are based on first‐ and second‐order upwind schemes for the convective terms. This new formulation is validated against experimental and other numerical data for well‐known benchmark problems, namely developing laminar flow in a straight duct, flow over a backward‐facing step, and lid‐driven cavity flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical study of the properties of the central moment lattice Boltzmann method

Central moment lattice Boltzmann method (LBM) is one of the more recent developments among the lattice kinetic schemes for computational fluid dynamics. A key element in this approach is the use of central moments to specify the collision process and forcing, and thereby naturally maintaining Galilean invariance, an important characteristic of fluid flows. When the different central moments are relaxed at different rates like in a standard multiple relaxation time (MRT) formulation based on raw moments, it is endowed with a number of desirable physical and numerical features. Because the collision operator exhibits a cascaded structure, this approach is also known as the cascaded LBM. While the cascaded LBM has been developed sometime ago, a systematic study of its numerical properties, such as the accuracy, grid convergence, and stability for well‐defined canonical problems is lacking, and the present work is intended to fulfill this need. We perform a quantitative study of the performance of the cascaded LBM for a set of benchmark problems of differing complexity, viz., Poiseuille flow, decaying Taylor–Green vortex flow, and lid‐driven cavity flow. We first establish its grid convergence and demonstrate second‐order accuracy under diffusive scaling for both the velocity field and its derivatives, that is, the components of the strain rate tensor, as well. The method is shown to quantitatively reproduce steady/unsteady analytical solutions or other numerical results with excellent accuracy. The cascaded MRT LBM based on the central moments is found to be of similar accuracy when compared with the standard MRT LBM based on the raw moments, when a detailed comparison of the flow fields are made, with both reproducing even the small scale vortical features well. Numerical experiments further demonstrate that the central moment MRT LBM results in significant stability improvements when compared with certain existing collision models at moderate additional computational cost. Copyright © 2015 John Wiley & Sons, Ltd.     

Finite Element Error Analysis Approach for Three-Dimensional Incompressible Viscous Fluid Flow Analysis

In our previous research, the modified Galerkin method was proposed as one of the most efficient methods for the analyses of convection-diffusion problems and two-dimensional viscous fluid flow problems. In this modified Galerkin method, the inertia term is considered explicitly, so only the symmetrical matrixes appear. Then an artificial viscosity is introduced through an error analysis approach to improve its accuracy and stability. In this paper, we proposed a new finite element formulation for three-dimensional incompressible viscous fluid flow analysis. This formulation (‘MS’ algorithm and ‘MSR’ algorithm) is based on the modified Galerkin method coupled with the Semi-Implicit Method for Pressure-Linked Equations. The cubic cavity flow problems were investigated for the Reynolds number of 400, 1,000, 2,000 and 3,200 using non-uniform meshes. Finally, we confirmed the effectiveness of our proposed method through the comparison with other research works.     

Comparative study of lattice‐Boltzmann and finite volume methods for the simulation of laminar flow through a 4:1 planar contraction

In the present paper, a comparative study of numerical solutions for Newtonian fluids based on the lattice‐Boltzmann method (LBM) and the classical finite volume method (FVM) is presented for the laminar flow through a 4:1 planar contraction at a Reynolds number of value one, Re=1. In this study, the stress field for LBM is directly obtained from the distribution function. The calculations of the stress based on the FVM‐data use the evaluations of velocity gradients with finite differences. The stress field for both LBM and FVM is expressed in the present study in terms of the shear stress and the first normal stress difference. The lateral and axial profiles of the velocity, the shear stress and the first normal stress difference for both methods are investigated. It is shown that the LBM results for the velocity and the stresses are in excellent agreement with the FVM results. Copyright © 2004 John Wiley & Sons, Ltd.     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.