Remarks on the links between low‐order DG methods and some finite‐difference schemes for the Stokes problem

In this paper we demonstrate that some well‐known finite‐difference schemes can be interpreted within the framework of the local discontinuous Galerkin (LDG) methods using the low‐order piecewise solenoidal discrete spaces introduced in (SIAM J. Numer. Anal. 1990; 27 (6): 1466–1485). In particular, it appears that it is possible to derive the well‐known MAC scheme using a first‐order Nédélec approximation on rectangular cells. It has been recently interpreted within the framework of the Raviart–Thomas approximation by Kanschat (Int. J. Numer. Meth. Fluids 2007; published online). The two approximations are algebraically equivalent to the MAC scheme, however, they have to be applied on grids that are staggered on a distance h/2 in each direction. This paper also demonstrates that both discretizations allow for the construction of a divergence‐free basis, which yields a linear system with a ‘biharmonic’ conditioning. Both this paper and Kanschat (Int. J. Numer. Meth. Fluids 2007; published online) demonstrate that the LDG framework can be used to generalize some popular finite‐difference schemes to grids that are not parallel to the coordinate axes or that are unstructured. Copyright © 2008 John Wiley & Sons, Ltd.     

A high‐order element‐based Galerkin method for the barotropic vorticity equation

A high‐order element‐based Galerkin method is developed to solve the non‐divergent barotropic vorticity equation (BVE). The solution process involves solving a conservative transport equation for the vorticity fields and a Poisson equation for the stream function fields. The discontinuous Galerkin method is employed for solving the transport equation and a spectral element method (continuous Galerkin) is used for the Poisson equation. A third‐order strong stability preserving explicit Runge–Kutta scheme is used for time integration. A series of tests have been performed to validate the model, which include the evolution of an idealized tropical cyclone and interaction of dual vortices in close proximity. The numerical convergence study is performed by solving the BVE on the sphere where the analytic solution is known. The test results are consistent with physical observations, and the model exhibits exponential convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.     

On the scarcity of contravariant left adjunctions

It is shown that any pair of contravariant mutually left adjoint functors between varieties of algebrasC andD factors through the retractions of these categories onto the complete lattices of epimorphs of their respective initial objects, and is thus fairly degenerate. However, some interesting contravariant left adjoint functors do exist between categories of finitely generated (or otherwise restricted) algebras.This work was done while the author was partly supported by NSF contract DMS 85-02330.Presented by F. E. J. Linton.     

TVD scheme for computing open channel wave flows

For the shallow water equations in the first approximation (Saint-Venant equations), a TVD scheme is developed for shock-capturing computations of open channel flows with discontinuous waves. The scheme is based on a special nondivergence approximation of the total momentum equation that does not involve integrals related to the cross-section pressure force and the channel wall reaction. In standard divergence difference schemes, most of the CPU time is spent on the computation of these integrals. Test computations demonstrate that the discontinuity relations reproduced by the scheme are accurate enough for actual discontinuous wave propagation to be numerically simulated. All the qualitatively distinct solutions for a dam collapsing in a trapezoidal channel with a contraction in the tailwater area are constructed as an example.     

Nonlinear stability of large k-value shallow spherical shells with a symmetrically distributed line load

In this paper, we treat the nonlinear stability problem of shallow spherical shells with large values ofk(k=12(1–v) · 2f/h,f = shell rise,h = shell thickness) under the action of uniformly distributed line load along a circle concentric with the shell boundary. Load-deflection curves are computed at successive increments of uniformly distributed line loads by using both cubic B-spline approximations and iterative techniques. Our algorithm yields fairly good convergent results for values ofk as large as 400. The limiting case in which shells are loaded along a circle of small radius has been specially investigated and the computed critical loads are compared with those obtained with central point loads by other authors.     

具线性退化特征场对角型方程组Goursat问题的奇性形成机制

For an inhomogeneous quasilinear hyperbolic system of diagonal form, under the assumptions that the system is linearly degenerate and the C^1 norm of the boundary data is bounded, we show that the mechanism of the formation of singularities of C^1 classical solution to the Goursat problem with C^1 compatibility conditions at the origin must be an ODE type. The similar result is also obtained for the weakly discontinuous solution with C^0 compatibility conditions at the origin.     

Optimal Runge–Kutta smoothers for the p-multigrid discontinuous Galerkin solution of the 1D Euler equations

This work presents a family of original Runge–Kutta methods specifically designed to be effective relaxation schemes in the numerical solution of the steady state solution of purely advective problems with a high-order accurate discontinuous Galerkin space discretization and a p-multigrid solution algorithm. The design criterion for the construction of the Runge–Kutta methods here developed is different form the one traditionally used to derive optimal Runge–Kutta smoothers for the h-multigrid algorithm, which are designed in order to provide a uniform damping of the error modes in the high-frequency range only. The method here proposed is instead designed in order to provide a variable amount of damping of the error modes over the entire frequency spectrum. The performance of the proposed schemes is assessed in the solution of the steady state quasi one-dimensional Euler equations for two test cases of increasing difficulty. Some preliminary results showing the performance for multidimensional applications are also presented.     

Artificial boundary conditions for the numerical solution of the Euler equations by the discontinuous galerkin method

We present artificial boundary conditions for the numerical simulation of compressible flows using high-order accurate discretizations with the discontinuous Galerkin (DG) finite element method. The construction of the proposed boundary conditions is based on characteristic analysis and applied for boundaries with arbitrary shape and orientation. Numerical experiments demonstrate that the proposed boundary treatment enables to convect out of the computational domain complex flow features with little distortion. In addition, it is shown that small-amplitude acoustic disturbances could be convected out of the computational domain, with no significant deterioration of the overall accuracy of the method. Furthermore, it was found that application of the proposed boundary treatment for viscous flow over a cylinder yields superior performance compared to simple extrapolation methods.     

水下物体视觉形象的分析与模拟研究

水下物体的视觉形象是一个看似简单但极易出错的问题.本文对此问题进行了较为系统的分析和模拟计算,纠正了一些流行的错误看法,指出在不同情况下水下物体的表观位置、表观大小及表观形状均会发生通常直观想象难以预见的复杂变化.该问题的讨论提供了一个培养学习者创新思维和利用基本知识分析及解决实际问题的能力的极好例证.