Source Identification for the Time-Fractional Diffusion Equation With Robin Boundary Conditions北大核心CSCD

对Robin边界条件时间分数阶扩散方程的源项辨识问题进行了研究。这类问题是不适定的,因此提出了一种迭代型正则化方法,得到了源项辨识问题的正则近似解。给出了先验和后验参数选取规则下正则近似解和精确解之间的误差估计,数值算例验证了该方法的有效性。     

A High-Order Finite Difference Scheme for 3D Unsteady Convection Diffusion Reaction Equations北大核心CSCD

针对三维非稳态对流扩散反应方程,构造了一种高精度紧致有限差分格式,对空间的离散采用四阶紧致差分方法,对时间的离散采用Taylor级数展开和余项修正技术,所提格式在时间上的精度为二阶、在空间上的精度为四阶。利用Fourier稳定性分析法证明了该格式是无条件稳定的。最后给出数值算例验证了理论结果。     

Asymptotic preserving HLL schemes

This work concerns the derivation of HLL schemes to approximate the solutions of systems of conservation laws supplemented by source terms. Such a system contains many models such as the Euler equations with high friction or the M1 model for radiative transfer. The main difficulty arising from these models comes from a particular asymptotic behavior. Indeed, in the limit of some suitable parameter, the system tends to a diffusion equation. This article is devoted to derive HLL methods able to approximate the associated transport regime but also to restore the suitable asymptotic diffusive regime. To access such an issue, a free parameter is introduced into the source term. This free parameter will be a useful correction to satisfy the expected diffusion equation at the discrete level. The derivation of the HLL scheme for hyperbolic systems with source terms comes from a modification of the HLL scheme for the associated homogeneous hyperbolic system. The resulting numerical procedure is robust as the source term discretization preserves the physical admissible states. The scheme is applied to several models of physical interest. The numerical asymptotic behavior is analyzed and an asymptotic preserving property is systematically exhibited. The scheme is illustrated with numerical experiments. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1396–1422, 2011     

Weakly Nonlinear Geometric Optics for Hyperbolic Systems of Conservation Laws

We present a new approach to analyze the validation of weakly nonlinear geometric optics for entropy solutions of nonlinear hyperbolic systems of conservation laws whose eigenvalues are allowed to have constant multiplicity and corresponding characteristic fields to be linearly degenerate. The approach is based on our careful construction of more accurate auxiliary approximation to weakly nonlinear geometric optics, the properties of wave front-tracking approximate solutions, the behavior of solutions to the approximate asymptotic equations, and the standard semigroup estimates. To illustrate this approach more clearly, we focus first on the Cauchy problem for the hyperbolic systems with compact support initial data of small bounded variation and establish that the L ¹-estimate between the entropy solution and the geometric optics expansion function is bounded by O(?²), independent of the time variable. This implies that the simpler geometric optics expansion functions can be employed to study the behavior of general entropy solutions to hyperbolic systems of conservation laws. Finally, we extend the results to the case with non-compact support initial data of bounded variation.     

A Class of Cauchy Problem for Quasilinear Hyperbolic Systems with the Weakly Dissipative Terms

ＡＣｌａｓｓｏｆＣａｕｃｈｙＰｒｏｂｌｅｍｆｏｒＱｕａｓｉｌｉｎｅａｒＨｙｐｅｒｂｏｌｉｃＳｙｓｔｅｍｓｗｉｔｈｔｈｅＷｅａｋｌｙＤｉｓｓｉｐａｔｉｖｅＴｅｒｍｓ￥ＬｉｕＦａｇｕｉ；ＹａｎｇＱｉａｏ（ＮｏｒｔｈＣｈｉｎａＩｎｓｔｉｔｕｔｅｏｆＷａｔｅ...     

带有源项的广义Chaplygin气体磁流体Euler方程组Riemann解的极限

研究了带有源项的广义Chaplygin气体磁流体Euler方程组Riemann解的极限.由于非齐次项的影响,带有源项的广义Chaplygin气体磁流体Euler方程组Riemann解不再是自相似的.当压力和磁感强度同时消失时,它的解会收敛到零压流输运方程组的Riemann解,解中会出现δ-激波和真空现象.同时研究还得到了仅当磁感强度消失时,它的解会收敛到非齐次广义Chaplygin气体Euler方程组的Riemann解,并且解中只出现δ-激波.     

求解双曲守恒律方程的三阶修正模板WENO格式北大核心CSCD

为了降低经典的三阶加权本质无振荡(WENO)格式的数值耗散,提出了一种新的三阶WENO格式的修正模板近似方法.改进了经典WENO-JS格式中各候选模板上数值通量的一阶多项式逼近,通过加入二次项使模板逼近达到三阶精度.计算了相应的候选通量,并且通过引入可调函数φ(x),使得新的格式具有ENO性质.最后给出了一系列数值算例,证明了该方法的有效性.     

具耗散项拟线性双曲型方程组的边值问题

本文考虑具耗散项一阶拟线性双曲型方程组的具有自由边界的典型自由边值问题．在一些合理假设下,证明了其经典解的整体存在性定理．     

用WENO方法求解双曲型守恒律方程组的初(边)值问题

本文用WENO算法解决双曲型守恒律方程组初（边值）问题．给出一种满足熵条件、Sδ熵条件和边界熵条件的WENO算法．通过这个算法就能得到守恒律方程组的数值解，数值解和理论解是非常吻合的．     

An Efficient Method for Computing Hyperbolic Systems with Geometrical Source Terms Having Concentrations

We propose a simple numerical method for calculating both unsteady and steady state solution of hyperbolic system with geometrical source terms having concentrations. Physical problems under consideration include the shallow water equations with topography,and the quasi one-dimensional nozzle flows. We use the interface value, rather than the cell-averages, for the source terms, which results in a well-balanced scheme that can capture the steady state solution with a remarkable accuracy. This method approximates the source terms via the numerical fluxes produced by an (approximate) Riemann solver for the homogeneous hyperbolic systems with slight additional computation complexity using Newton‘s iterations and numerical integrations. This method solves well the subor super-critical flows, and with a transonic fix, also handles well the transonic flows over the concentration. Numerical examples provide strong evidence on the effectiveness of this new method for both unsteady and steady state calculations.     

ON THE CENTRAL RELAXING SCHEME Ⅱ: SYSTEMS OF HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are interested in construction of the central reltalng sChemes for system of noIilinearhyperbolic conservation lawswith initial data U(0, x) = Uo(x), x = (x1 ? ...! xd), based on the local relaJxation approkimationof Eq.(1.1) [2, 3, 6, 8, 9, 12].To i11ustrate the basic idea of the relaalng schemes, for the sake of simplicity in the presentation, we restrict our attention to onedimensional scalar conservaioll lawsFirst, introduce a linear hyperbollc system with a stiff sourc…     

一种新的求解双曲型守恒律的三阶中心差分格式

本文提出了一种求解双曲型守恒律新的三阶中心差分格式，主要是引入了一种推广的三阶重构，并证明了这种重构在网格边界无振荡．所提的格式保持了中心差分格式简单的优点，不需用Riemann解算器，避免了进行特征解耦．数值试验结果表明本文格式是高精度、高分辨率的。