间断解问题的有限体积法

有限体积法（FVM）是集有限差分法和有限元方法之优点而发展起来的一种新的数值方法,首先利用守恒问题给出FVM的基本概念和基本过程,接下来就间断解问题的求解进一步介绍了FVM的一些比较活跃的构造方法,并给出了用其中几个方法模拟二维溃坝问题的数值结果,然后按问题分类介绍了FVM数值方法和实际应用的进展,最后给出了FVM理论分析的一些进展。     

数值求解瞬态温度场的有限体积法

本文介绍了求解三维瞬态热传导的有限体积法,给出了控制方程和数值方法。通过一些算例对该方法进行了验证,数值试验结果表明该方法计算稳定、可靠。     

基于样条逼近有限体积法的通量分裂算法

研究了任意曲线坐标系中求解双曲型守恒律的高精度、无波动样条逼近有限体积方法,比较了三种不同的通量分裂技术在这种方法中的应用。数值实验表明,在不同的通量分裂技术下,该方法对流场中的激波和接触间断都有很高的分辨率。     

摄动有限体积法重构近似高精度的意义

研讨有限体积(FV)方法重构近似高精度的作用问题.FV方法中积分近似采用中点规则为二阶精度时,重构近似高精度(精度高于二阶)的意义和作用是一个有争议的问题.利用数值摄动技术[1,2]构造了标量输运方程的积分近似为二阶精度、重构近似为任意阶精度的迎风型和中心型摄动有限体积(PFV)格式.迎风PFV格式无条件满足对流有界准则(CBC),中心型PFV格式为正型格式,两者均不会产生数值振荡解.利用PFV格式求解模型方程的数值结果表明:与一阶迎风和二阶中心格式相比,PFV格式精度高、对解的间断分辨率高、稳定性好、雷诺数的适用范围大,数值地"证实"重构近似高精度和PFV格式的实际意义和好处.     

激光熔蚀金属平板的三维有限体积法数值分析

为了研究气流条件下强激光对金属靶的熔蚀效应,采用有限体积方法建立了数值模型,并开发了三维Fortran计算程序。综合考虑强激光与材料耦合规律、光束能量空间分布、材料高温热物理性能以及熔蚀界面移动等关键影响因素,模拟了激光辐照下金属靶板升温、熔化和剥蚀的复杂物理过程。最后,将计算结果与试验数据进行了比较,验证了计算模型和程序的有效性。结果表明,计算模型能够反映强激光熔蚀金属平板的基本规律,熔蚀深度和后表面温度计算值与试验吻合较好,并且自编计算程序简单高效。     

激光熔蚀金属平板的三维有限体积法数值分析

为了研究气流条件下强激光对金属靶的熔蚀效应,采用有限体积方法建立了数值模型,并开发了三维Fortran计算程序。综合考虑强激光与材料耦合规律、光束能量空间分布、材料高温热物理性能以及熔蚀界面移动等关键影响因素,模拟了激光辐照下金属靶板升温、熔化和剥蚀的复杂物理过程。最后,将计算结果与试验数据进行了比较,验证了计算模型和程序的有效性。结果表明,计算模型能够反映强激光熔蚀金属平板的基本规律,熔蚀深度和后表面温度计算值与试验吻合较好,并且自编计算程序简单高效。     

空间差分格式对有限体积法数值散射的影响

数值散射是辐射传递方程近似算法中最常见的离散误差。本文主要讨论空间差分格式对有限体积法数值散射的影响。构造激光平行及倾斜入射的物理模型,验证和比较阶梯格式、中心差分格式及指数格式下温度场的计算精度及数值散射特性。计算结果表明,在激光平行入射与倾斜入射两种情况下,阶梯格式引起的的数值散射比菱形格式及指数格式要多,但其计算精度高于菱形及指数格式。不同激光入射条件下,各种差分格式表现出的数值散射分布有明显的差异。     

基于有限体积法数值模拟双撕裂模非线性演化

以磁流体理论为基础,采用基于有限体积法的通量差分分裂格式数值求解具有双曲保守律形式的电阻磁流体方程组。编写C++程序对平板几何位形下的等离子体双撕裂模进行了长时间数值模拟,得到双撕裂模不稳定性的演化图景,捕捉到了双撕裂模非线性发展过程中磁场重联的几个典型阶段,讨论了等离子体电阻和两个有理面之间的距离对双撕裂模不稳定性非线性发展的影响。为研究磁流体动力学提供了一种可行的高精度数值算法。     

有限体积法求解圆柱形散射介质内辐射与导热耦合换热

将谱带模型与有限体积解法相结合;求吸收、发射、散射性非灰介质圆柱体内辐射传递方程。考虑辐射强度场与热扩散温度场的耦合,将控制容积法与有限体积法结合,求解辐射与导热耦合换热。经与光线踪迹法、离散传递法的计算结果比较表明,谱带模型与有限体积解法相结合能处理多场耦合下非灰介质内的辐射换热。     

一种基于多矩的有限体积浸入边界法

基于多矩VSIAM3格式及浸入边界法,提出一套在复杂计算区域内求解不可压缩流动的数值格式.不可压N-S方程使用VSIAM3格式进行离散,引入浸入边界法处理复杂、移动边界,使用虚拟网格方法计算动量方程修正项,同时还考虑了对连续方程的修正.使用标准算例对数值模式进行验证.     

A hierarchical fracture model for the iterative multiscale finite volume method

An iterative multiscale finite volume (i-MSFV) method is devised for the simulation of multiphase flow in fractured porous media in the context of a hierarchical fracture modeling framework. Motivated by the small pressure change inside highly conductive fractures, the fully coupled system is split into smaller systems, which are then sequentially solved. This splitting technique results in only one additional degree of freedom for each connected fracture network appearing in the matrix system. It can be interpreted as an agglomeration of highly connected cells; similar as in algebraic multigrid methods. For the solution of the resulting algebraic system, an i-MSFV method is introduced. In addition to the local basis and correction functions, which were previously developed in this framework, local fracture functions are introduced to accurately capture the fractures at the coarse scale. In this multiscale approach there exists one fracture function per network and local domain, and in the coarse scale problem there appears only one additional degree of freedom per connected fracture network. Numerical results are presented for validation and verification of this new iterative multiscale approach for fractured porous media, and to investigate its computational efficiency. Finally, it is demonstrated that the new method is an effective multiscale approach for simulations of realistic multiphase flows in fractured heterogeneous porous media.     

谐振管内非线性驻波的有限体积数值算法

为了扩展谐振管内非线性驻波在工程中的应用, 以及克服现有数值计算方法仅局限于求解直圆柱形和指数形谐振管内非线性驻波的问题. 根据变截面的非稳态可压缩热黏性流体Navier-Stokes方程和空间守恒方程, 并基于求解压力速度耦合方程的半隐式算法和交错网格技术, 构建一种能够计算任意形状轴对称谐振管受活塞驱动时内部非线性驻波的有限体积算法. 分别对圆柱形、指数形和圆锥形谐振管内的非线性驻波进行仿真计算. 通过与现有试验结果以及数值仿真结果的对比, 验证了该方法的正确性.并获得除驻波声压之外的另外一些新的物理结果, 包括速度、密度、温度的瞬时变化.在直圆柱形谐振管内产生冲击声压波, 速度波形中出现钉状结构.而在指数形和圆锥形谐振管内产生高声压幅值的驻波, 没有出现冲击波, 速度波形中均未发现钉状结构. 计算结果表明谐振管内非线性驻波的物理属性与谐振管形状之间有密切关系.     

耦合势有限体积法高效模拟各向异性地层中海洋可控源的三维电磁响应

本文基于电场矢势与标势分解的耦合势有限体积法研究建立一套各向异性地层中海洋可控源电磁法的三维响应的高效数值模拟技术.首先引入电场的矢势和标势,将电场分解为无散场和无旋场之和,Maxwell方程转换为关于矢势与标势的混合Helmholtz方程,克服低感应数问题.在此基础上,借助Yee氏交错网格和有限体积法以及非均质单元中等效电导率公式,建立混合Helmholtz方程的离散方程.并采用直接法求解器PARDISO求解离散方程,有效保证在大的求解空间中仍然能够获得电磁场稳定可靠的数值解.此外,在数值模拟中利用差异场技术,克服源的奇异性问题,尽可能提高近场的计算精度.与解析解的对比证明了该算法的有效性.数值模拟结果表明,海洋可控源电磁法沿测线方向的电场,对油气藏的纵向电阻率敏感,对横向电阻率不敏感;对油气藏上方的覆盖层的纵向电阻率和横向电阻率都敏感.     

Integrating 1D and 2D hydrodynamic,sediment transport model for dam-break flow using finite volume method

The purpose of this study is to set up a dynamically linked 1D and 2D hydrodynamic and sediment transport models for dam break flow.The 1D-2D coupling model solves the generalized shallow water equations,the non-equilibrium sediment transport and bed change equations in a coupled fashion using an explicit finite volume method.It considers interactions among transient flow,strong sediment transport and rapid bed change by including bed change and variable flow density in the flow continuity and momentum equations.An unstructured Quadtree rectangular grid with local refinement is used in the 2D model.The intercell flux is computed by the HLL approximate Riemann solver with shock captured capability for computing the dry-to-wet interface for all models.The effects of pressure and gravity are included in source term in this coupling model which can simplify the computation and eliminate numerical imbalance between source and flux terms.The developed model has been tested against experimental and real-life case of dam-break flow over fix bed and movable bed.The results are compared with analytical solution and measured data with good agreement.The simulation results demonstrate that the coupling model is capable of calculating the flow,erosion and deposition for dam break flows in complicated natural domains.     

���������������ֵģ��˫˺��ģ�������ݻ�

以磁流体理论为基础,采用基于有限体积法的通量差分分裂格式数值求解具有双曲保守律形式的电阻磁流体方程组.编写C++程序对平板几何位形下的等离子体双撕裂模进行了长时间数值模拟,得到双撕裂模不稳定性的演化图景,捕捉到了双撕裂模非线性发展过程中磁场重联的几个典型阶段,讨论了等离子体电阻和两个有理面之间的距离对双撕裂模不稳定性非线性发展的影响.为研究磁流体动力学提供了一种可行的高精度数值算法.     

A coupled edge-/face-based smoothed finite element method for structural-acoustic problems

The edge-based smoothed finite element method (ES-FEM) and the face-based smoothed finite element method (FS-FEM) developed recently have shown great efficiency in solving solid mechanics problems with triangular and tetrahedral meshes. In this paper, a coupled ES-/FS-FEM model is extended to solve the structural-acoustic problems consisting of a plate structure interacting with the fluid medium. Three-node triangular elements and four-node tetrahedral elements are used to discretize the two-dimensional (2D) plate and three-dimensional (3D) fluid, respectively, as they can be generated easily and even automatically for complicated geometries. The field variable in each element is approximated using the linear shape functions, which is exactly the same as that in the standard FEM. The gradient field of the problem is obtained particularly using the gradient smoothing operation over the edge-based and face-based smoothing domains in 2D and 3D, respectively. The gradient smoothing technique can provide a proper softening effect to the model, effectively solve the problems caused by the well-known “overly-stiff” phenomenon existing in the standard FEM, and hence significantly improve the accuracy of the solution for the coupled systems. Intensive numerical studies have been conducted to verify the effectiveness of the coupled ES-/FS-FEM for structural-acoustic problems.     

A multi-moment finite volume formulation for shallow water equations on unstructured mesh

A novel and accurate finite volume method has been presented to solve the shallow water equations on unstructured grid in plane geometry. In addition to the volume integrated average (VIA moment) for each mesh cell, the point values (PV moment) defined on cell boundary are also treated as the model variables. The volume integrated average is updated via a finite volume formulation, and thus is numerically conserved, while the point value is computed by a point-wise Riemann solver. The cell-wise local interpolation reconstruction is built based on both the VIA and the PV moments, which results in a scheme of almost third order accuracy. Efforts have also been made to formulate the source term of the bottom topography in a way to balance the numerical flux function to satisfy the so-called C-property. The proposed numerical model is validated by numerical tests in comparison with other methods reported in the literature.