双剪统一弹塑性有限差分方法研究

基于拉格朗日有限差分方法,建立了双剪统一弹塑性有限差分计算格式,并利用VC++语言编写动态链接库文件将双剪统一弹塑性模型导入拉格朗日有限差分程序FLAC(Fast Lagrangian Analysis of Continua)中进行计算分析。双剪统一弹塑性有限差分方法可以模拟复杂应力状态下结构的渐进破坏,无需形成刚度矩阵,对于材料非线性问题无需进行迭代计算,因此在理论和工程应用中都有积极的意义。本文利用双剪统一弹塑性有限差分方法对拉压强度不等材料的厚壁圆筒受内压、中心带孔板条受拉压、条形基础下的地基极限分析及边坡问题进行了数值分析并与滑移线场等解析方法计算结果进行对比,结果均吻合较好。     

应用于计算气动声学的优化有限紧致格式

对于含间断的计算气动声学问题,数值计算的格式不仅要求低耗散低色散的设计,对短波具有较高的分辨率,还要求能捕捉激波.中心紧致格式具有高精度,具有无耗散和低色散特征,但不能捕捉间断和激波;WENO格式处理间断较为成功,而耗散和色散误差相对较大.有限紧致格式可以将紧致格式与WENO格式相结合构造成混合格式,利用光滑因子之间的关系对激波区域进行自动判断,将传统的全域求解的紧致格式划分为有限的局部紧致求解,间断点上的激波捕捉铜梁自动作为局部紧致求解的边界通量,在在光滑区域具有紧致格式的高精度低耗散性质,在激波附近不产生非物理振荡.本文利用有限紧致格式思想,构造了新的适合于气动声学问题的优化有限紧致格式,将其应用于计算气动声学一维标准测试问题,对相关格式的模拟性能进行了评估,显示该格式在宽频声波传播和含有间断的声波传播模拟方面具有优势.     

基于非结构/混合网格的高阶精度格式研究进展

尽管以二阶精度格式为基础的计算流体力学(CFD) 方法和软件已经在航空航天飞行器设计中发挥了重要的作用, 但是由于二阶精度格式的耗散和色散较大, 对于湍流、分离等多尺度流动现象的模拟, 现有成熟的CFD 软件仍难以给出满意的结果, 为此CFD 工作者发展了众多的高阶精度计算格式. 如果以适应的计算网格来分类, 一般可以分为基于结构网格的有限差分格式、基于非结构/混合网格的有限体积法和有限元方法,以及各种类型的混合方法. 由于非结构/混合网格具有良好的几何适应性, 基于非结构/混合网格的高阶精度格式近年来备受关注. 本文综述了近年来基于非结构/混合网格的高阶精度格式研究进展, 重点介绍了空间离散方法, 主要包括k-Exact 和ENO/WENO 等有限体积方法, 间断伽辽金(DG) 有限元方法, 有限谱体积(SV) 和有限谱差分(SD) 方法, 以及近来发展的各种DG/FV 混合算法和将各种方法统一在一个框架内的CPR (correctionprocedure via reconstruction) 方法等. 随后简要介绍了高阶精度格式应用于复杂外形流动数值模拟的一些需要关注的问题, 包括曲边界的处理方法、间断侦测和限制器、各种加速收敛技术等. 在综述过程中, 介绍了各种方法的优势与不足, 其间介绍了作者发展的基于"静动态混合重构" 的DG/FV 混合算法. 最后展望了基于非结构/混合网格的高阶精度格式的未来发展趋势及应用前景.     

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法,本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法,基于有限体积方法和迎风型矢通量分裂方法,建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件,区域内部采用隐式格式和局部时间步长等,以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较,二者基本吻合。计算结果表明采用分析“V”型多重网格方法,能提高计算效率,加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。     

FPM改进算法及在应力波传播问题中的应用

有限粒子法(FPM)是传统SPH方法的重要发展,大大提高了边界区域粒子的计算精度。然而在迭代计算过程中,高耗时和潜在的数值不稳定性是制约FPM应用的关键因素。通过对FPM基本方程进行矩阵分解,建立了一种特殊格式的FPM改进算法。该方法保持FPM方法在边界区域较高计算精度的同时,成功地规避了传统FPM方法对系数矩阵可逆性的限制,大大提高了计算效率。最后,将改进算法在一维应力波传播问题中予以实现,获得了较好的数值结果。     

FPM改进算法及在应力波传播问题中的应用A specified improved algorithm for Finite Particle Method and its application to wave propagation

有限粒子法(FPM)是传统SPH方法的重要发展,大大提高了边界区域粒子的计算精度。然而在迭代计算过程中,高耗时和潜在的数值不稳定性是制约FPM应用的关键因素。通过对FPM基本方程进行矩阵分解,建立了一种特殊格式的FPM改进算法。该方法保持FPM方法在边界区域较高计算精度的同时,成功地规避了传统FPM方法对系数矩阵可逆性的限制,大大提高了计算效率。最后,将改进算法在一维应力波传播问题中予以实现,获得了较好的数值结果。     

非理想爆源产生的爆炸场数值模拟

本文采用有限差分方法TVD(Total Variation Diminishing,A. Harten,1983)格式对非理想爆源在地面上形成的二维轴对称爆炸场进行了数值模拟,数值计算值与试验结果吻合得相当好。     

非理想爆源产生的爆炸场数值模拟

本文采用有限差分方法——TVD(Total Variation Diminishing,A. Harten,1983)格式对非理想爆源在地面上形成的二维轴对称爆炸场进行了数值模拟,数值计算值与试验结果吻合得相当好。     

基于制造解的非结构二阶有限体积离散格式的精度测试与验证

随着计算机技术的飞速进步,计算流体力学得到迅猛发展,数值计算虽能够快速得到离散结果,但是数值结果的正确性与精度则需要通过严谨的方法来进行验证和确认.制造解方法和网格收敛性研究作为验证与确认的重要手段已经广泛应用于计算流体力学代码验证、精度分析、边界条件验证等方面.本文在实现标量制造解和分量制造解方法的基础上,通过将制造解方法精度测试结果与经典精确解(二维无黏等熵涡)精度测试结果进行对比,进一步证实了制造解精度测试方法的有效性,并将两种制造解方法应用于非结构网格二阶精度有限体积离散格式的精度测试与验证,对各种常用的梯度重构方法、对流通量格式、扩散通量格式进行了网格收敛性精度测试.结果显示,基于Green-Gauss公式的梯度重构方法在不规则网格上会出现精度降阶的情况,导致流动模拟精度严重下降,而基于最小二乘(least squares)的梯度重构方法对网格是否规则并不敏感.对流通量格式的精度测试显示,所测试的各种对流通量格式均能达到二阶精度,且各方法精度几乎相同;而扩散通量离散中界面梯度求解方法的选择对流动模拟精度有显著影响.     

多维迎风方法在非结构/混合网格热流计算中的应用研究

非结构/混合网格具有极强的几何灵活性,在复杂外形飞行器的气动力特性数值模拟中已得到广泛应用,但目前还难以准确地预测气动热环境。本文从非结构/混合网格热流计算的三个需求出发,选取了多维迎风方法,并与其他方法进行了对比研究。以二维圆柱高超声速绕流这一Benchmark典型问题为例,对比研究了多维迎风方法和几种广泛使用的无粘通量格式(Roe格式、Van Leer格式和AUSMDV格式)对混合网格热流计算精度的影响。结果表明,多维迎风方法在热流计算精度、鲁棒性以及收敛性方面表现良好。最后,将多维迎风方法应用于常规混合网格上的圆柱和钝双锥绕流问题,均得到了较好的热流计算结果,为非结构/混合网格热流计算在复杂高超飞行器中的应用奠定了基础。     

An application of the capacitance matrix method to accommodate masked land areas and island circulations in a primitive equation ocean model

The capacitance matrix method has been implemented in a primitive equation ocean model to accommodate islands and portions of irregular coastal boundaries that cannot be treated adequately by boundary-fitted orthogonal curvilinear co-ordinates. The algorithm preserves the ability to solve the streamfunction equation using fast and accurate elliptic solvers that require a rectangular computational domain. By superposition of a set of island Green functions, the solution is adjusted to ensure continuity of pressure around each island. The implementation is tested by comparison with an analytic solution for wind-driven flow in a closed basin similar to the southwest Pacific Ocean.     

Motion of the interface of two fluids in a porous medium

In this study we obtain the Integro-differential equation for the motion of the interface of two incompressible fluids in various well areal arrangement systems. The solution of the equation is presented for a five-point system in the form of a power series with respect to time. Formulas are assumed which describe the motion of the particles belonging to the interface along invariant streamlines for five-point, seven-point, and nine-point well arrangement systems. The stratum sweeping coefficients for the fluid which is displacing the stratum oil are calculated (under conditions of the five-point system) at the instant when the fluid breaks through into the operation wells. The results of the calculations are compared with experimental data [1].The author wishes to thank V. L. Danilov for valuable counsel and comments.     

A VARIATIONAL INEQUALITY APPROACH FOR NON-STEADY SEEPAGE FLOW THROUGH THREE-DIMENSIONAL FRACTURE NETWORK

为了求解裂隙岩体有自由面非稳定渗流问题,将Darcy定律延拓至整个研究区域,使得潜在溢出边界条件满足Signorini型边界条件,建立了三维裂隙网络非稳定渗流问题的抛物型变分不等式（parabolic variational inequality,PVI）提法,并证明其与偏微分方程（partial differential equation,PDE）提法的等价性,从而将自由面上的流量条件以及潜在溢出边界上的互补条件转化成自然边界条件,降低该问题求解难度。同时给出了基于PVI提法的有限元数值求解方法,通过与交叉裂隙模型理论解的对比分析,证明了该方法的正确性。最后将该方法对含复杂三维裂隙网络的边坡进行非稳定渗流分析,计算结果表明该方法对于复杂裂隙网络求解具有较强的可靠性和适应性。     

River and reservoir flow modelling using the transformed shallow water equations

This paper describes a versatile finite difference scheme for the solution of the two-dimensional shallow water equations on boundary-fitted non-orthogonal curvilinear meshes. It is believed that this is the first non-orthogonal shallow water equation model incorporating the advective acceleration terms to have been developed in the United Kingdom. The numerical scheme has been validated against the severe condition of jet-forced flow in a circular reservoir with vertical side walls, where reflections of the initial free surface waves pose major problems in achieving a stable solution. Furthermore, the validation exercises are designed to test the computer model for artificial diffusion, which may be a consequence of the numerical scheme adopted to stabilize the shallow water equations. The model is shown to be capable of simulating the flow conditions in an irregularly shaped domain typical of the geometries frequently encountered in civil engineering river basin management.     

Finite difference solution of a Newtonian jet swell problem

A finite difference technique has been developed to study the Newtonian jet swell problem. The streamfunction and vorticity were used as dependent variables to describe the jet flow. The boundary-fitted co-ordinate transformation method was adopted to map the flow geometry into a rectangular domain. The standard finite difference method was then applied for solving the flow equations. The location of the jet free surface was updated by the kinematic boundary condition, and an adjustable parameter was included in the free-surface iteration. We could obtain numerical solutions for the Reynolds number as high as 100, and the differences between the present study and previous finite element simulations on the jet swell ratio are less than 5%.     

Finite element computations of two-dimensional arterial flow in the presence of a transverse magnetic field

A finite element solution of the Navier-Stokes equations for steady flow under the magnetic effect through a double-branched two-dimensional section of a three-dimensional model of the canine aorta is discussed. The numerical scheme involves transforming the physical co-ordinates to a curvilinear boundary-fitted co-ordinate system. The shear stress at the wall is calculated for a Reynolds number of 1000 with the branch-to-main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and found to be in reasonable qualitative agreement. The steady flow, shear stress and branch flow under the effect of a magnetic field have been discussed in detail.     

一种改进格林元方法及在渗流问题中的应用

采用格林公式和基本解推导出直接边界积分方程来求解渗流问题.边界积分方程数值离散基于格林元方法(Green element methond),改进了原方法中压力和压力导数的求解方法,命名为混合边界元方法(Mixed boundary element method).相较于格林元类方法,该方法显式考虑了求解节点的外法向流量值和压力值,并使求得的数值解在求解区域上能够连续,符合实际的物理过程,在不增加额外未知数的情况下提高了计算精度.分析了不同网格类型对模拟计算结果的影响,并对稳定渗流问题、非稳定(瞬态)渗流问题和非稳态问题进行了实例计算,结果显示改进方法提高了计算精度,并对各类渗流问题有较好的适应性.     

Comparison of two numerical methods for the solution of non-Newtonian flow in ducts

Two numerical methods, the Galerkin finite element method (FEM) and the boundary-fitted co-ordinate transformation method (BFCTM), have been applied to solve inelastic non-Newtonian fluid flow in ducts of irregular cross-section. Three representative fluid models, namely the power-law, the Ellis and the Bingham models, have been analysed. The application of the FEM is straightforward, while for the BFCTM the accurate estimation of viscosity on the duct boundary and the proper mesh adjustment appear to be critical for generating convergent solutions. A detailed comparison of the two numerical methods in terms of volumetric flow rate, axial velocity, shear rate, viscosity and CPU time is given. Both methods can generate accurate solutions of velocity over a wide range of variables, but the FEM requires much less computing time to reach the same level of accuracy. Only the BFCTM can be used to approximate shear rate and viscosity with reasonable accuracy.     

A Quasilinear Based Procedure for Saturated/Unsaturated Water Movement in Soils

The quasilinear form of Richards equation for one-dimensional unsaturated flow in soils can be readily solved for a wide variety of conditions. However, it cannot explain saturated/unsaturated flow and the constant diffusivity assumption, used to linearise the transient quasilinear equation, can introduce significant error. This paper presents a quasi-analytical solution to transient saturated/unsaturated flow based on the quasilinear equation, with saturated flow explained by a transformed Darcy's equation. The procedure presented is based on the modified finite analytic method. With this approach, the problem domain is divided into elements, with the element equations being solutions to a constant coefficient form of the governing partial differential equation. While the element equations are based on a constant diffusivity assumption, transient diffusivity behaviour is incorporated by time stepping. Profile heterogeneity can be incorporated into the procedure by allowing flow properties to vary from element to element. Two procedures are presented for the temporal solution; a Laplace transform procedure and a finite difference scheme. An advantage of the Laplace transform procedure is the ability to incorporate transient boundary condition behaviour directly into the analytical solutions. The scheme is shown to work well for two different flow problems, for three soil types. The technique presented can yield results of high accuracy if the spatial discretisation is sufficient, or alternatively can produce approximate solutions with low computational overheads by using large sized elements. Error was shown to be stable, linearly related to element size.     

A Computational Method for Free-Surface Viscous Flows Through Porous Media

In this paper we have investigated free-surface unsteady viscous flow through porous media using a finite difference technique via primitive formulation. A transformation has been used to convert the irregular flow domain into a rectangular domain. Two illustrations have been presented in case of a standing wave problem and a cavity flow problem with a free surface. The method is found satisfactory for numerical solution for Reynolds number upto 200.