The Homogeneous Hamilton–Jacobi and Bernoulli Equations Revisited, II

It is shown that the admissible solutions of the continuity and Bernoulli or Burgers' equations of a perfect one-dimensional liquid are conditioned by a relation established in 1949–1950 by Pauli, Morette, and Van Hove, apparently, overlooked so far, which, in our case, stipulates that the mass density is proportional to the second derivative of the velocity potential. Positivity of the density implies convexity of the potential, i.e., smooth solutions, no shock. Non-elementary and symmetric solutions of the above equations are given in analytical and numerical form. Analytically, these solutions are derived from the original Ansatz proposed in Ref. 1 and from the ensuing operations which show that they represent a particular case of the general implicit solutions of Burgers' equation. Numerically and with the help of an ad hoc computer program, these solutions are simulated for a variety of initial conditions called compatible if they satisfy the Morette–Van Hove formula and anti-compatible if the sign of the initial velocity field is reversed. In the latter case, singular behaviour is observed. Part of the theoretical development presented here is rephrased in the context of the Hopf–Lax formula whose domain of applicability for the solution of the Cauchy problem of the homogeneous Hamilton–Jacobi equation has recently been enlarged.     

A composition law for Runge-Kutta methods applied to index-2 differential-algebraic equations

This paper presents a new composition law for Runge-Kutta methods when applied to index-2 differential-algebraic systems. Applications of this result to the study of the order of composite methods and of symmetric methods are given.     

Concretely Locally Presentable and Locally Generated Categories

We give concrete versions of the characterizations of locally -presentable (resp. -generated) categories as categories of models of (resp. weakly) -ary limit-theories and of (resp. weakly) -ary essentially algebraic theories. By concrete we mean that starting with a category of -structures, the theories obtained are extensions of the original ones, and the equivalences of categories are concrete isomorphisms. -presentable and -generated objects are also investigated from this viewpoint.     

On the iterative solution of the algebraic equations in fully implicit Runge-Kutta methods

This paper deals with the iterative solution of stage equations which arise when some fully implicit Runge-Kutta methods, in particular those based on Gauss, Radau and Lobatto points, are applied to stiff ordinary differential equations. The error behaviour in the iterates generated by Newton-type and, particularly, by single-Newton schemes which are proposed for the solution of stage equations is studied. We consider stiff systems y'(t) = f(t,y(t)) which are dissipative with respect to a scalar product and satisfy a condition on the relative variation of the Jacobian of f(t,y) with respect to y, similar to the condition considered by van Dorsselaer and Spijker in [7] and [17]. We prove new convergence results for the single-Newton iteration and derive estimates of the iteration error that are independent of the stiffness. Finally, some numerical experiments which confirm the theoretical results are presented. This revised version was published online in June 2006 with corrections to the Cover Date.     

Unconditionally Stable Complex Envelope Wave Equation PML Algorithm for Band Limited FDTD Simulations

Unconditionally stable complex envelope (CE) perfectly matched layer (PML) absorbing boundary conditions (ABCs) are presented for truncating the scalar wave-equation finite difference time domain (WE-FDTD) grids. The formulations are based on incorporating the alternating direction implicit (ADI) scheme into the CE FDTD implementations of the scalar wave-equation derived in the PML region at the domain boundaries. Numerical example carried out in two dimensional domain shows that the proposed formulations are more accurate than the classical ADI scalar wave equation PML formulations when it is used for modelling band limited electromagnetic applications.     

A second-order accurate numerical method for the two-dimensional fractional diffusion equation

Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are used in modeling practical superdiffusive problems in fluid flow, finance and others. In this paper, we present an accurate and efficient numerical method to solve a fractional superdiffusive differential equation. This numerical method combines the alternating directions implicit (ADI) approach with a Crank–Nicolson discretization and a Richardson extrapolation to obtain an unconditionally stable second-order accurate finite difference method. The stability and the consistency of the method are established. Numerical solutions for an example super-diffusion equation with a known analytic solution are obtained and the behavior of the errors are analyzed to demonstrate the order of convergence of the method.     

A conservative unstructured scheme for rapidly varied flows

This article introduces a new semi‐implicit, staggered finite volume scheme on unstructured meshes for the modelling of rapidly varied shallow water flows. Rapidly varied flows occur in the inundation of dry land during flooding situations. They typically involve bores and hydraulic jumps after obstacles such as road banks. Near such sudden flow transitions, the grid resolution is often low compared with the gradients of the bathymetry. Locally the hydrostatic pressure assumption may become invalid. In these situations, it is crucial to apply the correct conservation properties to obtain accurate results. An important feature of this scheme is therefore its ability to conserve momentum locally or, by choice, preserve constant energy head along a streamline. This is achieved using a special interpolation method and control volumes for momentum. The efficiency of inundation calculations with locally very high velocities, and in the case of unstructured meshes locally very small grid distances, is severely hampered by the Courant condition. This article provides a solution in the form of a locally implicit time integration for the advective terms that allows for an explicit calculation in most of the domain, while maintaining unconditional stability by implicit calculations only where necessary. The complex geometry of flooded urban areas asks for the flexibility of unstructured meshes. The efficient calculation of the pressure gradient in this, and other semi‐implicit staggered schemes, requires, however, an orthogonality condition to be put on the grid. In this article a simple method is introduced to generate unstructured hybrid meshes that fulfil this requirement. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical experiments with several variant WENO schemes for the Euler equations

Numerical experiments with several variants of the original weighted essentially non‐oscillatory (WENO) schemes (J. Comput. Phys. 1996; 126 :202–228) including anti‐diffusive flux corrections, the mapped WENO scheme, and modified smoothness indicator are tested for the Euler equations. The TVD Runge–Kutta explicit time‐integrating scheme is adopted for unsteady flow computations and lower–upper symmetric‐Gauss–Seidel (LU‐SGS) implicit method is employed for the computation of steady‐state solutions. A numerical flux of the variant WENO scheme in flux limiter form is presented, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of low‐order schemes. Computations of unsteady oblique shock wave diffraction over a wedge and steady transonic flows over NACA 0012 and RAE 2822 airfoils are presented to test and compare the methods. Various aspects of the variant WENO methods including contact discontinuity sharpening and steady‐state convergence rate are examined. By using the WENO scheme with anti‐diffusive flux corrections, the present solutions indicate that good convergence rate can be achieved and high‐order accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Copyright © 2008 John Wiley & Sons, Ltd.     

计算最小奇异组的一个精化调和 Lanczos双对角化方法

在很多实际应用中需要计算大规模矩阵的若干个最小奇异组.调和投影方法是计算内部特征对的常用方法,其原理可用于求解大规模奇异值分解问题.本文证明了,当投影空间足够好时,该方法得到的近似奇异值收敛,但近似奇异向量可能收敛很慢甚至不收敛.根据第二作者近年来提出的精化投影方法的原理,本文提出一种精化的调和Lanczos双对角化方法,证明了它的收敛性.然后将该方法与Sorensen提出的隐式重新启动技术相结合,开发出隐式重新启动的调和Lanczos双对角化算法(IRHLB)和隐式重新启动的精化调和Lanczos双对角化算法(IRRHLB).位移的合理选取是算法成功的关键之一,本文对精化算法提出了一种新的位移策略,称之为"精化调和位移".理论分析表明,精化调和位移比IRHLB中所用的调和位移要好,且可以廉价可靠地计算出来.数值实验表明,IRRHLB比IRHLB要显著优越,而且比目前常用的隐式重新启动的Lanczos双对角化方法(IRLB)和精化算法IRRLB更有效.