二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法（英文）

构造矩形网格下求解Lagrangian坐标系下气动方程组的单元中心型格式.空间离散采用控制体积间断Petrov-Galerkin方法,时间离散采用二阶TVD Runge-Kutta方法.利用限制器来抑制非物理震荡并保证RKCV算法的稳定性.构造的算法可以保证物理量的局部守恒.与Runge-Kutta间断Galerkin(RKDG)方法相比较,RKCV方法的计算公式少一项积分项使得计算较简单.给出一些数值算例验证了算法的可靠性及效率.     

基于DCS5格式的LDDRK算法

基于五阶线性耗散紧致格式(DCS5)和七级龙格原库塔时间积分算法,根据数值增长因子对精确增长因子的最佳逼近原则,提出与DCS5格式耗散性相适应的优化方法,并得到相应的七级五阶低耗散低色散龙格原库塔(LDDRK)算法.求解标量线性对流方程和线化Euler方程得到的一维波传播问题的数值结果显示,七级五阶LDDRK算法的精度优于七级七阶精度的标准龙格原库塔算法.     

二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法

构造矩形网格下求解Lagrangian坐标系下气动方程组的单元中心型格式. 空间离散采用控制体积间断Petrov-Galerkin方法,时间离散采用二阶TVD Runge-Kutta方法. 利用限制器来抑制非物理震荡并保证RKCV算法的稳定性. 构造的算法可以保证物理量的局部守恒. 与Runge-Kutta间断Galerkin（RKDG）方法相比较,RKCV方法的计算公式少一项积分项使得计算较简单. 给出一些数值算例验证了算法的可靠性及效率.     

耦合广义非线性薛定谔方程的相互作用表象龙格库塔算法及其误差分析

本文通过表象变换, 将耦合广义非线性薛定谔方程 (C-GNLSE) 变换成相互作用表象中的向量方程, 再利用向量形式的4阶龙格-库塔迭代格式, 建立了一种在频域内求解C-GNLSE的同步更新迭代算法. 通过将该向量形式的相互作用表象中的4阶龙格-库塔 (V-JH-RK4IP) 算法应用于高双折射光子晶体光纤中超连续谱产生的数值模拟, 验证了算法的有效性, 通过与现有其他典型算法的比较, 表明以V-JH-RK4IP算法求解C-GNLSE具有最高的计算精度和计算效率. 关键词： 耦合广义非线性薛定谔方程(C-GNLSE) 相互作用表象 4阶龙格-库塔算法 超连续谱产生     

聚变－裂变混合反应堆堆芯等离子体动力学方程组数值求解

用四阶龙格-库塔方法求解聚变-裂变混合反应堆堆芯等离子体动力学方程组,求解过程采用分段步长法,既可保证求解精度又可节省机时,所得出的解的曲线能充分反映堆芯等离子体的燃烧过程。求解过程还引用数值反馈控制方法用以对解曲线的优化和控制。这种数值反馈控制方法可以在工程技术上用来探索堆芯等离子体燃烧过程的控制。     

基于近似Riemann解的有限体积ALE方法

研究二维平面坐标系和二维轴对称坐标系中四边形网格上可压缩流体力学的有限体积ALE(Arbitrary Lagrangian Eulerian)方法.数值方法采用节点中心有限体积法,数值通量采用适用于任意状态方程的HLLC(Harten-Lax-Van Leer-Collela)通量.空间二阶精度通过用WENO(weighted essentially non-oscillatory)方法对原始变量进行重构获得,时间离散采用两步显式Runge-Kutta格式.数值例子显示,方法具有良好的激波分辨能力和高精度的数值逼近能力.     

隐-显积分因子间断Galerkin方法求解二维辐射扩散方程

提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是一种非常有效的数值算法.     

自适应间断有限元方法求解三维欧拉方程

将龙格库塔间断有限元方法(RDDG)与自适应方法相结合,求解三维欧拉方程.区域剖分采用非结构四面体网格,依据数值解的变化采用自适应技术对网格进行局部加密或粗化,减少总体网格数目,提高计算效率.给出四种自适应策略并分析不同自适应策略的优缺点.数值算例表明方法的有效性.     

粒子输运方程的线性间断有限元方法

将空间线性间断有限元方法应用于动态粒子输运方程的求解.数值算例表明,空间线性间断有限元方法在网格边界的数值精度方面明显高于指数格式和菱形格式,并且通量在时间上的微分曲线相对光滑,避免了指数格式、菱形格式数值解的非物理振荡现象.     

微重力环境下无限长柱体内角毛细流动解析近似解研究

应用同伦分析法研究无限长柱体内角毛细流动解析近似解问题,给出了级数解的递推公式.不同于其他解析近似方法,该方法从根本上克服了摄动理论对小参数的过分依赖,其有效性与所研究的非线性问题是否含有小参数无关,适用范围广.同伦分析法提供了选取基函数的自由,可以选取较好的基函数,更有效地逼近问题的解,通过引入辅助参数和辅助函数来调节和控制级数解的收敛区域和收敛速度,同伦分析法为内角毛细流动问题的解析近似求解开辟了一个全新的途径.通过具体算例,将同伦分析法与四阶龙格库塔方法数值解做了比较,结果表明,该方法具有很高的计算精度.     

求解辐射传递的非结构混合有限体积／有限元法

本文给了一种适用于任意非结构网格的有限体积/有限元法的混合算法用于求解多维半透明吸收、发射、散射性灰矩形介质内的辐射传递.该方法使用有限元法进行角度离散,有限体积法进行空间离散.与基于辐射传递离散坐标方程的方法不同的是,该方法在迭代求解的过程中,针对每一个空间体元,所有角度方向的辐射强度同时耦合求出.通过两个算例验证了该解法的正确性.     

多块对接网格技术在电磁场散射问题中的应用

采用多块对接网格技术和时域有限体积法(FVTD)研究了典型目标和多体的电磁场散射问题.控制方程采用三维一般曲线坐标系下的时变麦克斯韦方程组.时间方向采用四步Runge-Kutta方法,空间离散采用基于通量雅可比矩阵特征结构的通矢量差分分裂,依赖变量采用MUSCL插值.时间方向的计算精度为2阶,空间方向的计算精度可达3阶.典型算例的雷达散射截面(RCS)的计算结果与理论解吻合很好.对于多体问题计算与文献结果相当一致,说明该算法具有对复杂拓扑结构外形(包括多体问题)进行数值模拟的能力.     

Prediction of radiative heat transfer in 3D complex geometries using the unstructured control volume finite element method

In this paper, a 3D algorithm for the treatment of radiative heat transfer in emitting, absorbing, and scattering media is developed. The numerical approach is based on the utilization of the unstructured control volume finite element method (CVFEM) which, to the knowledge of the authors, is applied for the first time to simulate radiative heat transfer in participated media confined in 3D complex geometries. This simulation makes simultaneously the use of the merits of both the finite element method and the control volume method. Unstructured 3D triangular element grids are employed in the spatial discretization and azimuthal discretization strategy is employed in the angular discretization. The general discretization equation is presented and solved by the conditioned conjugate gradient squared method (CCGS). In order to test the efficiency of the developed method, several 3D complex geometries including a hexahedral enclosure, a 3D equilateral triangular enclosure, a 3D L-shaped enclosure and 3D elliptical enclosure are examined. The results are compared with the exact solutions or published references and the accuracy obtained in each case is shown to be highly satisfactory. Moreover, this approach required a less CPU time and iterations compared with those of even parity formulation of the discrete ordinates method.     

各向同性紊流能量衰减的大涡模拟

在多重网格驱动的,高效且得到充分验证的有限体积方法的基础上发展了可压缩流大涡模拟的方法.空间离散采用Jameson的中心格式附加二阶和四阶耗散的方法,时间推进则采用了双时间步长的方法.亚格子剪切应力张量和亚格子热通量用Smagorinsky模型进行模拟.通过对各向同性紊流能量衰减的模拟来验证本方法的准确性和高效性,模拟得到的能量谱和紊流动能衰减历程与过滤后的CBC实验数据吻合良好.     

On maximum-principle-satisfying high order schemes for scalar conservation laws

We construct uniformly high order accurate schemes satisfying a strict maximum principle for scalar conservation laws. A general framework (for arbitrary order of accuracy) is established to construct a limiter for finite volume schemes (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO) schemes) or discontinuous Galerkin (DG) method with first order Euler forward time discretization solving one-dimensional scalar conservation laws. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. It is straightforward to extend the method to two and higher dimensions on rectangular meshes. We also show that the same limiter can preserve the maximum principle for DG or finite volume schemes solving two-dimensional incompressible Euler equations in the vorticity stream-function formulation, or any passive convection equation with an incompressible velocity field. Numerical tests for both the WENO finite volume scheme and the DG method are reported.     

A RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in Lagrangian coordinate

In this paper, Runge-Kutta Discontinuous Galerkin (RKDG) finite element method is presented to solve the one-dimensional inviscid compressible gas dynamic equations in Lagrangian coordinate. The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method. A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method. For multi-medium fluid simulation, the two cells adjacent to the interface are treated differently from other cells. At first, a linear Riemann solver is applied to calculate the numerical flux at the interface. Numerical examples show that there is some oscillation in the vicinity of the interface. Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical flux at the interface, which suppress the oscillation successfully. Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.     

An RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in a Lagrangian coordinate

In this paper,Runge-Kutta Discontinuous Galerkin(RKDG) finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ?ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ?ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.     

一种高精度的修正Hermite-ENO格式

设计一种基于三单元具有六阶精度的修正Hermite-ENO格式（CHENO）,求解一维双曲守恒律问题.CHENO格式利用有限体积法进行空间离散,在空间层上,使用ENO格式中的Newton差商法自适应选择模板.在重构半节点处的函数值及其一阶导数值时,利用Taylor展开给出修正Hermite插值使其提高到六阶精度,并设计了间断识别法与相应的处理方法以抑制间断处的虚假振荡;在时间层上采用三阶TVD Runge-Kutta法进行函数值及一阶导数值的推进.其主要优点是在达到高阶精度的同时具有紧致性.数值实验表明对一维双曲守恒律问题的求解达到了理论分析结果,是有效可行的.