一种自由界面追踪的模板化VOF方法

发展了一种模板化的volume-of-fluid (VOF)方法．该方法根据自由界面的法向建立一个模板,然后由已知的网格单元上的流体体积比值确定出自由界面的准确位置,使得在二维情形下一个网格单元被自由界面切割的形式只有3种．另一方面,引入了单元边流体占有长度的概念,在此基础上建立了一个统一的流体占有面积模型,可以使得自由界面输运方程的求解有统一的算法．该方法不受网格单元形式的限制,并且容易推广到三维情形．算例表明,该方法能保证自由界面的跟踪精度．     

多介质大变形流动的MOF-MMALE数值模拟研究

多介质大变形流动数值模拟的关键和难点是在精确追踪物质界面的同时又能够处理好流体的大变形运动.将MOF(moment-of-fluid)界面重构算法与多介质任意Lagrange-Euler方法(MMALE)相耦合,形成MOF-MMALE方法,并应用于多介质大变形流动问题的数值模拟研究.MOF-MMALE方法在传统的ALE方法基础上,允许计算网格边界跨过物质界面,允许存在混合网格,即一个网格内可以存在两种或两种以上物质;在混合网格内,利用MOF界面重构算法来确定物质界面的位置和方向.数值算例表明,MOF-MMALE方法是模拟多介质大变形流动的有效手段,并且具有较好的数值精度和界面分辨率.     

垂直于双材料非完美界面Ⅱ型裂纹的断裂分析

研究了垂直于双材料非完美界面的Ⅱ型裂纹问题,采用线性弹簧模型模拟非完美界面．然后用Fourier积分变换方法把边值问题转化为求解具有Cauchy核的奇异积分方程,获得了裂纹两端应力强度因子的数值解．详细研究了问题的几种特例,并用数值实例分析了界面的非完美性对应力强度因子的影响．结果表明应力强度因子与界面参量有关并在完美界面和分离界面所对应的结果中变化．     

双材料接合半无限体三维矩形界面裂纹应力强度因子分析

基于体积力法,研究了双材料接合半无限体三维矩形界面裂纹的应力强度因子问题．在数值计算中,未知的体积力密度采用基本密度函数和多项式乘积的形式来近似,其中基本密度函数是根据界面裂纹应力的振荡奇异性来选取的．计算结果表明,基于本算法得到的数值结果其收敛精度和计算误差都是令人满意的．算例中,给出了应力强度因子随矩形形状及双材料参数的变化规律．     

位于两不同正交各向异性半平面间张开型界面裂纹的性能分析

利用Schmidt方法分析了位于正交各向异性材料中的张开型界面裂纹问题．经富立叶变换使问题的求解转换为求解两对对偶积分方程,其中对偶积分方程的变量为裂纹面张开位移．最终获得了应力强度因子的数值解．与以前有关界面裂纹问题的解相比,没遇到数学上难以处理的应力振荡奇异性,裂纹尖端应力场的奇异性与均匀材料中裂纹尖端应力场的奇异性相同．同时当上下半平面材料相同时,可以得到其精确解．     

求解带有移动界面的线性对流方程的浸入界面有限体积法（英文）

本文为一类带有移动界面的守恒律方程提出了耦合高分辨率格式的数值算法.这种算法是在一致大小的笛卡尔网格上导出而满足标准的双曲型稳定条件.文末列举数值算例研究这种算法的收敛性和数值精度.     

非一致性界面热流固耦合作用的整体求解

提出了非一致性界面热流固耦合作用整体求解的一种方法．热流体求解基于Boussinesq假设和不可压缩的Navier-Stokes方程．流体区域的运动采用任意Lagrange-Euler（ALE）方法．拟固体元方法实现流体区域的变形．使用几何非线性的热弹性动力学描述固体运动．为了保证界面处应力和传热的平衡,采用了基于Gauss积分点的数据交换方法,对热流固耦合最终形成的强非线性方程实现整体求解．数值实例分析表明该方法的健壮性和有效性．     

基于AUSM分裂的二维通量分裂格式

基于对流迎风分裂思想构造的AUSM类格式具有简单、高效、分辨率高等优点,在计算流体力学中得到了广泛的应用.传统的AUSM类格式在计算界面数值通量时只考虑网格界面法向的波系,忽略了网格界面横向波系的影响.使用Liou-Steffen通量分裂方法将二维Euler方程的通量分裂成对流通量和压力通量,采用AUSM格式来分别计算对流数值通量和压力数值通量.通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的AUSM通量分裂格式.在计算一维算例时,该格式保留了精确捕捉激波和接触间断的优点.在计算二维算例时,该格式不仅具有更高的分辨率而且表现出更好的鲁棒性,可以消除强激波波后的不稳定现象.此外,在多维问题的数值模拟中,该格式大大地提高了稳定性CFL数,具有更高的计算效率.因此,它是一种精确、高效并且强鲁棒性的数值方法.     

一维多介质可压缩流动的守恒型界面追踪方法

本文针对一维问题的ProntTracking方法,提出了一种较易实现的守恒型界面追踪方法.利用双波近似求解Riemann问题来确定界面处的数值通量,在固定的网格上采用统一的有限体积格式进行内点和交界面点的计算,通过守恒插值以及守恒量的重新分配,保证数值解在全场实现一致守恒,将该方法应用于一维多介质可压缩流动的模拟,给出了满意的数值模拟结果.     

预测纳米纤维复合材料有效弹性性能的界面模型和界面相模型北大核心CSCD

基于广义自洽法,同时采用Gurtin-Murdoch界面模型和界面相模型研究了纳米纤维复合材料的有效弹性性能,获得了两种模型下有效体积模量的封闭解析解和计算有效面内剪切模量数值解的全部公式.基于界面模型的解答,讨论了有效体积模量和有效面内剪切模量的界面效应.证明了界面模型的解答可由界面相模型的解答退化得到,其中有效体积模量可以实现解析退化,有效面内剪切模量则可以数值退化.以含纳米孔洞的金属铝为例,比较了两种模型计算结果的差异.结果表明,当纳米孔洞半径较小时,两个模型的结果存在很大差异,而当半径较大时两个模型的结果差别不大.     

A Characteristic Finite Volume Element Method for the Air Pollution Model

In this article, a characteristic finite volume element method is presented for solving air pollution models. The convection term is discretized using the characteristic method and diffusion term is approximated by finite volume element method. Compared with standard finite volume element method, our proposed method is more accurate and efficient, especially suitable to solve convection-dominated problems. The proposed numerical schemes are analyzed for convergence in L ₂ norm. Some numerical results are presented to demonstrate the efficiency and accuracy of the method.     

粘弹性方程全离散化有限体积元格式及数值模拟

本文利用有限体积元方法研究二维粘弹性方程, 给出一种时间二阶精度的全离散化有限体积元格式, 并给出这种全离散化有限体积元解的误差估计, 最后用数值例子验证数值结果与理论结果是相吻合的. 通过与有限元方法和有限差分方法相比较, 进一步说明了全离散化有限体积元格式是求解二维粘弹性方程数值解的最有效方法之一.     

A New Reduced-Order fve Algorithm Based on POD Method for Viscoelastic Equations

A proper orthogonal decomposition (POD) technique is used to reduce the finite volume element (FVE) method for two-dimensional (2D) viscoelastic equations. A reduced-order fully discrete FVE algorithm with fewer degrees of freedom and sufficiently high accuracy based on POD method is established. The error estimates of the reduced-order fully discrete FVE solutions and the implementation for solving the reduced-order fully discrete FVE algorithm are provided. Some numerical examples are used to illustrate that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced-order fully discrete FVE algorithm is one of the most effective numerical methods by comparing with corresponding numerical results of finite element formulation and finite difference scheme and that the reduced-order fully discrete FVE algorithm based on POD method is feasible and efficient for solving 2D viscoelastic equations.     

状态转换下美式跳扩散期权的修正Crank-Nicolson拟合有限体积法

主要研究了一类状态转换下美式跳扩散期权定价模型的修正Crank-Nicolson拟合有限体积法并且给出收敛性分析.文章所构造的新方法是对[Gan X T,Yin J F,Li R,Fitted finite volume method for pricing American options under regime-switching.jump-diffusion models based on penalty method.Adv.Appl.Math.Mech.,2020,12(3):748-773]中时间方向上Crank-Nicolson格式的改进.同时,还对求解非线性系统迭代方法的收敛性证明进行了补充.最后,数值实验验证了新方法的有效性.     

有限体积法与多重网格法用于提高激波捕捉质量及加速收敛

多重网格技术是一种非常有效的数值计算方法，本文采用多重网格的ＦＡＳ格式进行数值实验，计算加速效果十分明显，同时，结合矢通量分裂用有限体积法，大大提高了主激波的质量。     

An unstructured finite volume time domain method for structural dynamics

An unstructured finite volume time domain method (UFVTDM) is proposed to simulate stress wave propagation, in which the original variables of displacement and stress are solved based on the dynamic equilibrium equations. An Euler explicit and unstructured finite volume method is used for time dependent and spacial terms respectively. The displacements are stored on the cell vertex and a vertex based finite volume method is formed with that integral surface and the stresses are as assumed to be uniform in the cell. The present UFVTDM has several features. (1) The governing equations are discretized with the finite volume method which naturally follows conservation laws. (2) It can handle complex engineering problem. (3) This method is also able to analyze the natural characteristics and the numerical experiment shows that it is very efficient. Several cases are used to show the capability of the algorithm.     

非饱和土壤水流问题基于POD 方法的降阶有限体积元格式及外推算法实现

利用特征投影分解(proper orthogonal decomposition, 简记为POD) 方法对非饱和土壤水流问题的经典有限体积元格式做降阶处理, 建立一种具有足够高精度维数较低的降阶有限体积元格式, 并给出这种降阶有限体积元解的误差估计和外推算法的实现, 最后用数值例子说明数值结果与理论结果是相吻合的. 进一步表明了基于POD 方法的降阶有限体积元格式对求解非饱和土壤水流问题数值解是可靠和有效的.     

A posteriori error estimate for finite volume element method of the parabolic equations

In this article, residual‐type a posteriori error estimates are studied for finite volume element (FVE) method of parabolic equations. Residual‐type a posteriori error estimator is constructed and the reliable and efficient bounds for the error estimator are established. Residual‐type a posteriori error estimator can be used to assess the accuracy of the FVE solutions in practical applications. Some numerical examples are provided to confirm the theoretical results. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 259–275, 2017     

Stability and consistency of the semi-implicit co-volume scheme for regularized mean curvature flow equation in level set formulation

We show stability and consistency of the linear semi-implicit complementary volume numerical scheme for solving the regularized, in the sense of Evans and Spruck, mean curvature flow equation in the level set formulation. The numerical method is based on the finite volume methodology using the so-called complementary volumes to a finite element triangulation. The scheme gives the solution in an efficient and unconditionally stable way.     

A parallel finite volume scheme preserving positivity for diffusion equation on distorted meshes

Parallel domain decomposition methods are natural and efficient for solving the implicity schemes of diffusion equations on massive parallel computer systems. A finite volume scheme preserving positivity is essential for getting accurate numerical solutions of diffusion equations and ensuring the numerical solutions with physical meaning. We call their combination as a parallel finite volume scheme preserving positivity, and construct such a scheme for diffusion equation on distorted meshes. The basic procedure of constructing the parallel finite volume scheme is based on the domain decomposition method with the prediction‐correction technique at the interface of subdomains: First, we predict the values on each inner interface of subdomains partitioned by the domain decomposition. Second, we compute the values in each subdomain using a finite volume scheme preserving positivity. Third, we correct the values on each inner interface using the finite volume scheme preserving positivity. The resulting scheme has intrinsic parallelism, and needs only local communication among neighboring processors. Numerical results are presented to show the performance of our schemes, such as accuracy, stability, positivity, and parallel speedup.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2159–2178, 2017