DISCRETE GALERKIN METHOD FOR FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS

In this article, we develop a fully Discrete Galerkin（DG） method for solving ini- tial value fractional integro-differential equations（FIDEs）. We consider Generalized Jacobi polynomials（CJPs） with indexes corresponding to the number of homogeneous initial conditions as natural basis functions for the approximate solution. The fractional derivatives are used in the Caputo sense. The numerical solvability of algebraic system obtained from implementation of proposed method for a special case of FIDEs is investigated. We also provide a suitable convergence analysis to approximate solutions under a more general regularity assumption on the exact solution. Numerical results are presented to demonstrate the effectiveness of the proposed method.

基于高阶DG方法的非定常流场声辐射特性数值模拟研究

随着航空噪声越来越受到关注,计算声传播的算法成为研究热点。高阶间断伽辽金(Discontinuous Galerkin,DG)方法具有高精度、对网格质量要求低、适合自适应和并行计算等优点,可以以较高的效率对声场进行计算。文章运用高阶DG方法对线性化欧拉方程(Linearized Euler Equations,LEE)进行空间离散,并且基于离散后的线性化欧拉方程对带有背景流场的NACA0012翼型和30P30N多段翼型的声场进行数值计算。采用有限体积法计算得出流场信息后,通过插值将流场数据导入声场网格,并运用高阶DG方法进行声场计算。计算结果与参考文献中FW-H(Ffowcs Williams-Hawkings)算法对比一致性较好,验证了高阶DG算法的可行性。

Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

A Provable Positivity-Preserving Local Discontinuous Galerkin Method for the Viscous and Resistive MHD Equations

In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.

Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

A Local Macroscopic Conservative(LoMaC)Low Rank Tensor Method with the Discontinuous Galerkin Method for the Vlasov Dynamics

In this paper,we propose a novel Local Macroscopic Conservative(LoMaC)low rank tensor method with discontinuous Galerkin(DG)discretization for the physical and phase spaces for simulating the Vlasov-Poisson(VP)system.The LoMaC property refers to the exact local conservation of macroscopic mass,momentum,and energy at the discrete level.The recently developed LoMaC low rank tensor algorithm(arXiv:2207.00518)simultaneously evolves the macroscopic conservation laws of mass,momentum,and energy using the kinetic flux vector splitting;then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables.This paper is a generalization of our previous work,but with DG discretization to take advantage of its compactness and flexibility in handling boundary conditions and its superior accuracy in the long term.The algorithm is developed in a similar fashion as that for a finite difference scheme,by observing that the DG method can be viewed equivalently in a nodal fashion.With the nodal DG method,assuming a tensorized computational grid,one will be able to(i)derive differentiation matrices for different nodal points based on a DG upwind discretization of transport terms,and(ii)define a weighted inner product space based on the nodal DG grid points.The algorithm can be extended to the high dimensional problems by hierarchical Tucker(HT)decomposition of solution tensors and a corresponding conservative projection algorithm.In a similar spirit,the algorithm can be extended to DG methods on nodal points of an unstructured mesh,or to other types of discretization,e.g.,the spectral method in velocity direction.Extensive numerical results are performed to showcase the efficacy of the method.

Fully Discrete Nonlinear Galerkin Methods for Kuramoto-Sivashinsky Equation and Their Error Estimates

In this paper,the uniform error estimates with respect to t∈[0, ∞ ) of the nonlinear Galerkin method are given for the long time integration of the Kuramoto-Sivashinsky equation. The nonlinear Galerkin method is used to study the asymptotic behaviour of Kuramoto-Sivashinsky equation and to construct the bifurcation diagrams.

FULL DISCRETE NONLINEAR GALERKIN METHOD FOR THE NAVIER-STOKES EQUATIONS

This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.

ON GALERKIN DISCRETIZATION OF AXIALLY MOVING NONLINEAR STRINGS

A computational technique is proposed for the Galerkin discretization of axially moving strings with geometric nonlinearity. The Galerkin discretization is based on the eigenfunctions of stationary strings. The discretized equations are simplified by regrouping nonlinear terms to reduce the computation work. The scheme can be easily implemented in the practical programming. Numerical results show the effectiveness of the technique. The results also highlight the feature of Galerkin＇s discretization of gyroscopic continua that the term number in Galerkin＇s discretization should be even. The technique is generalized from elastic strings to viscoelastic strings.

On Stability of Full Discrete Nonlinear Gálerkin Method(Dedicated to Professor You Zhaoyong for his 60th brithday)

This paper discuss stability of the full discrete nonlinear Galerkin method based on the approximation inertial manifold method for some nonlinear evolution equation, for example, some nonlinear reactor equation and Navier-Stokes Equation. In the paper we provide some necessary and sufficient conditions of stability.

A multithreaded parallel upwind sweep algorithm for the S_(N) transport equations discretized with discontinuous finite elements

The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.

Kirchhoff型抛物方程的Galerkin有限元法的超收敛误差分析

文章研究了后向Euler全离散Galerkin格式下的Kirchhoff型抛物方程的超收敛误差分析。首先,讨论了数值解的先验误差估计,并证明了数值解的存在唯一性。其次,使用双线性元的高精度误差估计以及Ritz投影算子与插值算子相结合的技术,通过技巧性地处理非线性项得到了超逼近的误差估计结果。再次,通过插值后处理方法获得了整体的超收敛结果。最后,通过数值试验验证了理论分析的正确性。

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG／FV 混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度 DG 方法计算量、存储量小等优点，自提出以来已成功应用于一维、二维标量方程和 Euler／N-S 方程的求解。综述了 DG／FV 混合方法的研究进展，重点介绍了 DG／FV 混合方法的空间重构算法、针对 RANS 方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析，并给出了系列粘性流算例的计算结果，包括用于验证混合方法数值精度的库埃特流，以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明 DG／FV 混合方法达到了设计的精度阶，且相比同阶 DG 方法计算量减少约40％，而隐式方法能大幅提高定常流的收敛历程，较显式 Runge-Kutta 的收敛速度提高1～2个量级。

高阶精度 DG/FV 混合方法在二维粘性流动模拟中的推广

DG／FV 混合方法因其具有紧致性、易于推广至高阶及相比同阶 DGM 计算量、存储量小等优点，已成功应用于一维／二维标量方程和 Euler 方程的求解。在此基础上，将该方法推广于二维三角形／矩形混合网格上的 Navier-Stokes 方程数值模拟，将格式形式精度提高至4～5阶。物理量的空间重构及离散使用 DG／FV 混合重构方法；无粘通量计算采用 Roe 格式；粘性通量计算采用 BR2格式；时间方向离散采用高阶显式 R-K 方法或隐式方法。利用该方法计算了有解析解的 Couette 流动问题以验证几种格式的数值精度阶，并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明 DG／FV 混合方法达到了设计的精度阶，在较粗的网格上亦能得到高精度的计算结果；定性分析和数值结果表明相比同阶 DG 方法单步计算量减少约40％。

r型网格自适应在间断Galerkin有限元激波捕捉中的应用

间断Galerkin(DG)有限元方法因计算精度高、适用于非结构网格等特点得到广泛研究和应用,其在数值模拟包含强间断流场时存在残差收敛性和计算鲁棒性差问题,均匀分布的网格加剧这一问题并影响激波分辨率。针对该问题,发展了r型网格自适应方法,实现间断Galerkin有限元数值模拟过程中网格自适应加密。基于网格点归一化的压力值作为r型网格自适应中网格点移动驱动力的重要权值,并将网格自适应后的网格点位移变化量与网格点之间的初始位移之比作为驱动力的另一重要权值,实现网格沿激波方向各向异性自适应加密,并且激波附近网格点的相邻网格点同步向激波方向移动。发展了适合间断Galerkin有限元方法的Venkatakrishnan限制器。并列NACA0012翼型超声速算例及三维并列圆柱相互干扰算例结果表明:基于r型网格自适应的间断Galerkin有限元方法能够清晰锐利捕捉激波,提高模拟精度,具有良好的收敛性和鲁棒性。

离散系统运动方程的Galerkin有限元EEP法自适应求解

对于结构动力分析中的离散系统运动方程,现有算法的计算精度和效率均依赖于时间步长的选取,这是时间域问题求解的难点.基于EEP(element energy projection)超收敛计算的自适应有限元法,以EEP超收敛解代替未知真解,估计常规有限元解的误差,并自动细分网格,目前已对诸类以空间坐标为自变量的边值问题取得成功.对离散系统运动方程建立弱型Galerkin有限元解,引入基于EEP法的自适应求解策略,在时间域上自动划分网格,最终得到所求时域内任一时刻均满足给定误差限的动位移解,进而建立了一种时间域上的新型自适应求解算法.

Sobolev方程一个新的H^1-Galerkin混合有限元分析

研究了Sobolev方程的H^1-Galerkin混合有限元方法.利用不完全双二次元Q_2^-和一阶BDFM元,建立了一个新的混合元模式,通过Bramble-Hilbert引理,证明了单元对应的插值算子具有的高精度结果.进一步,对于半离散和向后欧拉全离散格式,分别导出了原始变量u在H^1-模和中间变量p在H(div)-模意义下的超逼近性质.

伪双曲型方程的一个H^1-Galerkin非协调混合元格式(英文)

讨论了一类伪双曲型方程的一个H1-Galerkin非协调混合有限元方法.利用插值算子的特殊性质,在半离散和全离散格式下,得到了与传统混合有限元相同的误差估计且不需要满足LBB条件.

拟线性粘弹性方程新H^1-Galerkin最低阶混合元格式的高精度分析

利用双线性元和零阶Raviart-Thomas元,针对拟线性粘弹性方程建立新的H^1-Galerkin混合元逼近格式.在半离散格式下,给出原始变量u的H^1模和应力=?ut的H(div;?)模的超逼近性和超收敛结果.同时,导出向后欧拉格式和Crank-Nicolson-Galerkin格式的最优误差估计.最后,通过数值算例表明逼近格式是有效的.

基于多区域面积分方程的电大非均匀等离子体电磁特性高效计算方法

针对临近空间等离子体包覆目标的电磁特性快速准确评估这一迫切需求,本文运用基于多区域面积分方程(multi-region surface integral equation, MR-SIE)的全波数值计算方法,展现电大尺度、非均匀、高负介电常数等离子体的电磁散射与辐射问题的仿真能力。首先,推导了非均匀等离子体的MR-SIE,使用间断伽辽金(discontinuous Galerkin, DG)方法提高对多区域复杂等离子体目标的建模与计算效率,研究了不同方程的数值性能;随后针对电子密度较大的高负介电常数区等离子体,运用一种简洁高效的截断策略,进一步提升了多层快速多极子算法的鲁棒性,避免了SIE计算高损耗介质内问题面临的数值不稳定性。数值实验表明,该方法在计算多种分割方式的多区域等离子体鞘套模型时具有良好的精度和效率,可用于大尺度复杂等离子体目标电磁辐射与散射特性的快速精确评估。