Discontinuous Galerkin method for hyperbolic equations involving \delta -singularities: negative-order norm error estimates and applications

In this paper, we develop and analyze discontinuous Galerkin (DG) methods to solve hyperbolic equations involving $\delta $ -singularities. Negative-order norm error estimates for the accuracy of DG approximations to $\delta $ -singularities are investigated. We first consider linear hyperbolic conservation laws in one space dimension with singular initial data. We prove that, by using piecewise $k$ th degree polynomials, at time $t$ , the error in the $H^{-(k+2)}$ norm over the whole domain is $(k+1/2)$ th order, and the error in the $H^{-(k+1)}(\mathbb R \backslash \mathcal R _t)$ norm is $(2k+1)$ th order, where $\mathcal R _t$ is the pollution region due to the initial singularity with the width of order $\mathcal O (h^{1/2} \log (1/h))$ and $h$ is the maximum cell length. As an application of the negative-order norm error estimates, we convolve the numerical solution with a suitable kernel which is a linear combination of B-splines, to obtain $L^2$ error estimate of $(2k+1)$ th order for the post-processed solution. Moreover, we also obtain high order superconvergence error estimates for linear hyperbolic conservation laws with singular source terms by applying Duhamel’s principle. Numerical examples including an acoustic equation and the nonlinear rendez-vous algorithms are given to demonstrate the good performance of DG methods for solving hyperbolic equations involving $\delta $ -singularities.     

Enhanced accuracy by post-processing for finite element methods for hyperbolic equations

We consider the enhancement of accuracy, by means of a simple post-processing technique, for finite element approximations to transient hyperbolic equations. The post-processing is a convolution with a kernel whose support has measure of order one in the case of arbitrary unstructured meshes; if the mesh is locally translation invariant, the support of the kernel is a cube whose edges are of size of the order of only. For example, when polynomials of degree are used in the discontinuous Galerkin (DG) method, and the exact solution is globally smooth, the DG method is of order in the -norm, whereas the post-processed approximation is of order ; if the exact solution is in only, in which case no order of convergence is available for the DG method, the post-processed approximation converges with order in , where is a subdomain over which the exact solution is smooth. Numerical results displaying the sharpness of the estimates are presented.

     

NONLINEAR INITIAL-BOUNDARY VALUE PROBLEM FOR QUASILINEAR HYPERBOLIC SYSTEM

Consider the nonlinear inltial-boundary value problem for quasilinear hyperbolicsystem:Let k≥2[n/2] 6,(F,g)∈ H~k(R_ ;Ω)×H~k(R_ ;Ω),and their traces at t=0 are zeroup to the order k-1.If for u=0,the problem(*)at t=0 is a Kreiss hyperbolic system,and the boundaryconditions satisfy the uniformly Lopatinsky criteria,then there exists a T>0 such that(*)has a unique H~k soluton in(0,T).In the Appendix,for symmetric hyperbolic systems,a comparison between theuniformly Lopatinsky condition and the stable admissible condition is given.     

一阶双曲问题间断有限元的后验误差分析

一阶双曲问题的有限元后验误差估计至今没有得到很好的解决.本文对d维区域上一阶双曲问题的k次间断有限元逼近提出了一种新的后验误差分析方法, 进而建立了间断有限元解在DG范数下(强于L₂范数)基于误差余量型的后验误差估计. 数值计算验证了本文理论分析的有效性. 本文方法也适用于其他变分问题有限元逼近的后验误差分析.     

Constructing Order Two Superconvergent WG Finite Elements on Rectangular Meshes

In this paper, we introduce a stabilizer free weak Galerkin (SFWG) finite element method for second order elliptic problems on rectangular meshes. With a special weak Gradient space, an order two superconvergence for the SFWG finite element solution is obtained, in both $L^2$ and $H^1$ norms. A local post-process lifts such a $P_k$ weak Galerkin solution to an optimal order $P_{k+2}$ solution. The numerical results confirm the theory.     

Convergence of Weak Galerkin Finite Element Method for Second Order Linear Wave Equation in Heterogeneous Media

Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space $(\mathcal{P}_k(K), \mathcal{P}_{k−1}(∂K), [\mathcal{P}_{k−1}(K)]^2).$ Optimal order a priori error estimates for both space-discrete scheme and implicit fully discrete scheme are derived in $L^∞(L^2)$ norm. This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes. Finite element algorithm presented here can contribute to a variety of hyperbolic problems where physical domain consists of heterogeneous media.     

Finite Difference Methods for the Heat Equation with a Nonlocal Boundary Condition

We consider the numerical solution by finite difference methods of the heat equation in one space dimension, with a nonlocal integral boundary condition, resulting from the truncation to a finite interval of the problem on a semi-infinite interval. We first analyze the forward Euler method, and then the $θ$-method for $0 < θ ≤ 1$, in both cases in maximum-norm, showing $O(h^2 + k)$ error bounds, where $h$ is the mesh-width and $k$ the time step. We then give an alternative analysis for the case $θ = 1/2$, the Crank-Nicolson method, using energy arguments, yielding a $O(h^2$ + $k^{3/2}$) error bound. Special attention is given the approximation of the boundary integral operator. Our results are illustrated by numerical examples.     

Uniqueness of Quasi-Regular Solutions for a Bi-Parabolic Elliptic Bi-Hyperbolic Tricomi Problem

The Tricomi equation $ yu_{xx} + u_{yy} = 0 $ was established in 1923 by Tricomi who is the pioneer of parabolic elliptic and hyperbolic boundary value problems and related problems of variable type. In 1945 Frankl established a generalization of these problems for the well-known Chaplygin equation $ K(\,y)u_{xx} + u_{yy} = 0 $ subject to the Frankl condition 1 + 2( K / K ')' > 0, y <0. In 1953 and 1955 Protter generalized these problems even further by improving the above Frankl condition. In 1977 we generalized these results in R n ( n > 2). In 1986 Kracht and Kreyszig discussed the Tricomi equation and transition problems. In 1993 Semerdjieva considered the hyperbolic equation $ K_1 (\,y)u_{xx} + (K_2 {\rm (\,}y{\rm )}u_y )_y + ru = f $ for y<0. In this paper we establish uniqueness of quasi-regular solutions for the Tricomi problem concerning the more general mixed type partial differential equation $ K_1 (\,y)(M_2 {\rm (}x{\rm )}u_x )_x + M_1 (x)(K_2 {\rm (\,}y{\rm )}u_y )_y + ru = f $ which is parabolic on both lines x = 0; y = 0, elliptic in the first quadrant x > 0, y > 0 and hyperbolic in both quadrants x< 0, y > 0; x > 0, y< 0. In 1999 we proved existence of weak solutions for a particular Tricomi problem. These results are interesting in fluid mechanics.     

一类高阶线性微分方程解的复振荡

In this paper,we shall use Nevanlinna theory of meromorphic functions to investigate the complex oscillation theory of solutions of some higher order linear differential equation.Suppose that A is a transcendental entire function with ρ(A)<1/2.Suppose that k≥2 and f(k)+A(z)f=0 has a solution f with λ(f)<ρ(A),and suppose that A1=A+h,where h≡0 is an entire function with ρ(h)<ρ(A).Then g(k)+A1(z)g=0 does not have a solution g with λ(g)<∞.     

Projektive Skalen,Tschebyscheff-Polynome und Fibonacci-Zahlen

Summary. Let $\widehat{\widehat T}_n$ and $\overline U_n$ denote the modified Chebyshev polynomials defined by $\widehat{\widehat T}_n (x) = {T_{2n + 1} \left(\sqrt{x + 3 \over 4} \right) \over \sqrt{x + 3 \over 4}}, \quad \overline U_{n}(x) = U_{n} \left({x + 1 \over 2}\right) \qquad (n \in \mathbb{N}_{0},\ x \in \mathbb{R}).$ For all $n \in \mathbb{N}_{0}$ define $\widehat{\widehat T}_{-(n + 1)} = \widehat{\widehat T}_n$ and $\overline U_{-(n + 2)} = - \overline U_n$, furthermore $\overline U_{-1} = 0$. In this paper, summation formulae for sums of type $\sum\limits^{+\infty}_{k = -\infty} \mathbf a_{\mathbf k}(\nu; x)$ are given, where $\bigl(\mathbf a_{\mathbf k}(\nu; x)\bigr)^{-1} = (-1)^k \cdot \Bigl( x \cdot \widehat{\widehat T}_{\left[k + 1 \over 2\right] - 1} (\nu) +\widehat{\widehat T}_{\left[k + 1 \over 2\right]}(\nu)\Bigr) \cdot \Bigl(x \cdot \overline U_{\left[k \over 2\right] - 1} (\nu) + \overline U_{\left[k \over 2\right]} (\nu)\Bigr)$ with real constants $ x, \nu $. The above sums will turn out to be telescope sums. They appear in connection with projective geometry. The directed euclidean measures of the line segments of a projective scale form a sequence of type $(\mathbf a_{\mathbf k} (\nu;x))_{k \in \mathbb{Z}}$ where $ \nu $ is the cross-ratio of the scale, and x is the ratio of two consecutive line segments once chosen. In case of hyperbolic $(\nu \in \mathbb{R} \setminus] - 3,1[)$ and parabolic $\nu = -3$ scales, the formula $\sum\limits^{+\infty}_{k = -\infty} \mathbf a_{\mathbf k} (\nu; x) = {\frac{1}{x - q_{{+}\atop(-)}}} - {\frac{1}{x - q_{{-}\atop(+)}}} \eqno (1)$ holds for $\nu > 1$ (resp. $\nu \leq - 3$), unless the scale is geometric, that is unless $x = q_+$ or $x = q_-$. By $q_{\pm} = {-(\nu + 1) \pm \sqrt{(\nu - 1)(\nu + 3)} \over 2}$ we denote the quotient of the associated geometric sequence.

     

Delta-Wave in a Simple 1-D 2x2 Hyperbolic Problem

A simple one-dimensional $2\times 2$ hyperbolic system is considered in the paper. The model contains a linear hyperbolic equation, as well as a hyperbolic equation of which the coefficients are about the solution of the linear one. The exact solution is presented and discussed, then numerical experiments are given by TVD (or MmB) type schemes for Riemann problems. From the results, we know that the solutions do have $\delta-$waves for some suitable initial data.     

On Poincaré-Friedrichs Type Inequalities for the Broken Sobolev Space ${\rm W}^{2,1}$

We are concerned with the derivation of Poincaré-Friedrichs type inequalities in the broken Sobolev space $W^{2,1}$($Ω$; $\mathcal{T}_h$) with respect to a geometrically conforming, simplicial triagulation $\mathcal{T}_h$ of a bounded Lipschitz domain $Ω$ in $\mathbb{R}^d$ , $d$ $∈$ $\mathbb{N}$. Such inequalities are of interest in the numerical analysis of nonconforming finite element discretizations such as ${\rm C}^0$ Discontinuous Galerkin (${\rm C}^0$${\rm DG}$) approximations of minimization problems in the Sobolev space $W^{2,1}$($Ω$), or more generally, in the Banach space $BV^2$($Ω$) of functions of bounded second order total variation. As an application, we consider a ${\rm C}^0$${\rm DG}$ approximation of a minimization problem in$BV^2$($Ω$) which is useful for texture analysis and management in image restoration.     

The Initial Value Problems and Nonlinear Boundary Value Problems for a Class of the Second order Quasilinear Parabolio Systems

In this paper initial value problems and nonlinear mixed boundary value problems for the quasilinear parabolic systems below $\[\frac{{\partial {u_k}}}{{\partial t}} - \sum\limits_{i,j = 1}^n {a_{ij}^{(k)}} (x,t)\frac{{{\partial ^2}{u_k}}}{{\partial {x_i}\partial {x_j}}} = {f_k}(x,t,u,{u_x}),k = 1, \cdots ,N\]$ are discussed.The boundary value conditions are $\[{u_k}{|_{\partial \Omega }} = {g_k}(x,t),k = 1, \cdots ,s,\]$ $\[\sum\limits_{i = 1}^n {b_i^{(k)}} (x,t)\frac{{\partial {u_k}}}{{\partial {x_i}}}{|_{\partial \Omega }} = {h_k}(x,t,u),k = s + 1, \cdots N.\]$ Under some "basically natural" assumptions it is shown by means of the Schauder type estimates of the linear parabolic equations and the embedding inequalities in Nikol'skii spaces,these problems have solutions in the spaces $\[{H^{2 + \alpha ,1 + \frac{\alpha }{2}}}(0 < \alpha < 1)\]$.For the boundary value problem with $\[b_i^{(k)}(x,t) = \sum\limits_{j = 1}^n {a_{ij}^{(k)}} (x,t)\cos (n,{x_j})\]$ uniqueness theorem is proved.     

A Study on CFL Conditions for the DG Solution of Conservation Laws on Adaptive Moving Meshes

The selection of time step plays a crucial role in improving stability and efficiency in the Discontinuous Galerkin (DG) solution of hyperbolic conservation laws on adaptive moving meshes that typically employs explicit stepping. A commonly used selection of time step is a direct extension based on Courant-Friedrichs-Levy (CFL) conditions established for fixed and uniform meshes. In this work, we provide a mathematical justification for those time step selection strategies used in practical adaptive DG computations. A stability analysis is presented for a moving mesh DG method for linear scalar conservation laws. Based on the analysis, a new selection strategy of the time step is proposed, which takes into consideration the coupling of the $α$-function (that is related to the eigenvalues of the Jacobian matrix of the flux and the mesh movement velocity) and the heights of the mesh elements. The analysis also suggests several stable combinations of the choices of the $α$-function in the numerical scheme and in the time step selection. Numerical results obtained with a moving mesh DG method for Burgers' and Euler equations are presented. For comparison purpose, numerical results obtained with an error-based time step-size selection strategy are also given.     

Tverberg-Type Theorems for Separoids

Let $S$ be a $d$-dimensional separoid of $(k-1)(d+1)+1$ convex sets in some "large-dimensional" Euclidean space $\E^N$. We prove a theorem that can be interpreted as follows: if the separoid $S$ can be mapped with a monomorphism to a $d$-dimensional separoid of points $P$ in general position, then there exists a $k$-colouring $\varsigma\colon \ S\to K_k$ such that, for each pair of colours $i,j\in K_k$, the convex hulls of their preimages do intersect---they are not separated. Here, by a monomorphism we mean an injective function such that the preimage of separated sets are separated. In a sense, this result is "dual" to the Hadwiger-type theorems proved by Goodman and Pollack (1988) and Arocha et al. (2002). We also introduce $\T(k,d)$, the minimum number $n$ such that all $d$-dimensional separoids of order at least $n$ can be $k$-coloured as before. By means of examples and explicit colourings, we show that for all $k>2$ and $d>0$, \[(k-1)(d+1)+1<\T(k,d)<{k\choose2}(d+1)+1.\] Furthermore, by means of a probabilistic argument, we show that for each $d$ there exists a constant $C=C(d)$ such that for all $k$, $\T(k,d)\leq Ck\log k$.     

An Extension of the Win Theorem: Counting the Number of Maximum Independent Sets

Win proved a well-known result that the graph G of connectivity κ(G) withα(G) ≤κ(G) + k-1(k ≥ 2) has a spanning k-ended tree, i.e., a spanning tree with at most k leaves. In this paper, the authors extended the Win theorem in case when κ(G) = 1 to the following: Let G be a simple connected graph of order large enough such that α(G) ≤ k + 1(k ≥ 3) and such that the number of maximum independent sets of cardinality k + 1 is at most n-2k-2. Then G has a spanning k-ended tree.     

Stabilized approximation to degenerate transport equations via filtering

We analyze a stabilization technique for degenerate transport equations. Of particular interest are coupled parabolic/hyperbolic problems, when the diffusion coefficient is zero in part of the domain. The unstabilized, computed approximations of these problems are highly oscillatory, and several techniques have been proposed and analyzed to mitigate the effects of the sub-grid errors that contribute to the oscillatory behavior. In this paper, we modify a time-relaxation algorithm proposed in [1] and further studied in [10]. Our modification introduces the relaxation operator as a post-processing step. The operator is not time-dependent, so the discrete (relaxation) system need only be factored once. We provide convergence analysis for our algorithm along with numerical results for several model problems.     

一类高阶周期微分方程解的性质

本文研究了高阶线性微分方程$$f^{(k)}(z)+A_{k-2}(z)f^{(k-2)}(z)+\cdots+A_0(z)f(z)=0,\eqno(*)$$解的线性相关性,其中$A_j(z)(j=0,2,\ldots,k-2)$是常数, $A_1$为非常数的的整周期函数,周期为$2\pi i$,且是$e^z$的有理函数.在一定条件下,我们给出了方程(*)解的表示.     

On the Positivity of Scattering Operators for Poincaré-Einstein Manifolds

In this paper, we mainly study the scattering operators for a Poincaré-Einstein manifold $(X^{n+1}, g_+)$, which define the fractional GJMS operators $P_{2\gamma}$ of order $2\gamma$ for $0<\gamma<\frac{n}{2}$ for the conformal infinity $(M, [g])$. We generalise Guillarmou-Qing's positivity results in [8] to the higher order case. Namely, if $(X^{n+1}, g_+)$ $(n\geq 5)$ is a hyperbolic Poincaré-Einstein manifold and there exists a smooth representative $g$ for the conformal infinity such that the scalar curvature $R_g$ is a positive constant and $Q_4$ is semi-positive on $(M, g)$, then $P_{2\gamma}$ is positive for $\gamma\in [1,2]$ and the first real scattering pole is less than $\frac{n}{2}-2$.     

Series Representation of Daubechies' Wavelets

1,Iotroduction.InthispaPerwe8tudytherepresentationofDaubechies'wavelets.DaubechiesI1]constructedaf4milyofcompartlysupportedregularscallngfUnctionsrk.(x)andtheassoci4tedregularwpeletsop.(x)(N32):where4.eL'(R)definedbythep0lyn0mia:withZq.(k)=1'q.(k)ER,k=0,1,')N-1.Itisknownthat[1]f0reachN32,k=Osuppgh.=[0,2N-l],suppop.=[-(N-1),N]andthewaveletop.generatesbyitsdilatiOnsandtranslati0nsan0rth0rn0rmalbasis{m.(2ix-k)}i,k6Z0fL'(R).Thefunctionsrk.andop.havebeenprovedtobeveryusefulinnumericalanal…