Error Bounds for Uniform Asymptotic Expansions—Modified Bessel Function of Purely Imaginary Order

The authors modify a method of Olde Daalhuis and Temme for representing the remainder and coefficients in Airy-type expansions of integrals.By using a class of rational functions,they express these quantities in terms of Cauchy-type integrals;these expressions are natural generalizations of integral representations of the coe？cients and the remainders in the Taylor expansions of analytic functions.By using the new representation,a computable error bound for the remainder in the uniform asymptotic expansion of the modified Bessel function of purely imaginary order is derived.     

Material models in principal stretches for elastodynamics

In the present work we deal with the conserving integration of elastic bodies undergoing finite deformations. In particular, we make use of constitutive laws formulated in terms of principal stretches. Most material models for hyperelastic isotropic materials are described in terms of principal stretches (Simo and Taylor [1]), like the Neo–Hooke material which is a special case of the Ogden material, or in invariants. The main advantage of principal stretches is the fact that they can be measured directly, which means that the numerical results can be compared easily with experimental ones, see for example, Ogden [2]. Moreover, it is advantageous to describe viscoelastic material behaviour (e.g. rubberlike materials) in terms of principal stretches. Concerning the discretization in space we apply the enhanced assumed strain (EAS) method, see Simo and Armero [3]. For the discretization in time we aim at numerical integrators which inherit fundamental conservation laws from the underlying continuous system. In particular, we propose an energy and momentum conserving time–stepping scheme which relies on the notion of a discrete gradient (or derivative) in the sense of Gonzalez [4]. The proposed approach starts from our previous developments in [5]. Numerical examples demonstrate the advantageous properties of the present formulation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Discrete Hodge operators

Summary. Many linear boundary value problems arising in computational physics can be formulated in the calculus of differential forms. Discrete differential forms provide a natural and canonical approach to their discretization. However, much freedom remains concerning the choice of discrete Hodge operators, that is, discrete analogues of constitutive laws. A generic discrete Hodge operator is introduced and it turns out that most finite element and finite volume schemes emerge as its specializations. We reap the possibility of a unified convergence analysis in the framework of discrete exterior calculus. Received November 26, 1999 / Revised version received November 2, 2000 / Published online May 30, 2001     

非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅲ):时间二阶精度的全离散格式

In this paper, a fully discrete format of nonlinear Galerkin mixed element method with two-step discretization of time for the non stationary conduction-convection problems is presented. The existence and the convergence of the fully discrete mixed element solution are shown. On the basis of [9] and [10], we have proved that the schemes have second-order convergence accuracy for the time discretization.     

Discrete time,finite state space mean field games

In this paper we study a mean field model for discrete time, finite number of states, dynamic games. These models arise in situations that involve a very large number of agents moving from state to state according to certain optimality criteria. The mean field approach for optimal control and differential games was introduced by Lasry and Lions (2006, 2007) [3], [4], [5]. The discrete time, finite state space setting is motivated both by its independent interest as well as by numerical analysis questions which appear in the discretization of the problems introduced by Lasry and Lions. The main contribution of this paper is the exponential convergence to equilibrium of the initial-terminal value problem.     

A ROBUST DISCRETIZATION OF THE REISSNER-MINDLIN PLATE WITH ARBITRARY POLYNOMIAL DEGREE

A numerical scheme for the Reissner-Mindlin plate model is proposed.The method is based on a discrete Helmholtz decomposition and can be viewed as a generalization of the nonconforming finite element scheme of Arnold and Falk[SIAM J.Numer.Anal.,26(6):1276-1290,1989].The two unknowns in the discrete formulation are the in-plane rotations and the gradient of the vertical displacement.The decomposition of the discrete shear variable leads to equivalence with the usual Stokes system with penalty term plus two Poisson equations and the proposed method is equivalent to a stabilized discretization of the Stokes system that generalizes the Mini element.The method is proved to satisfy a best-approximation result which is robust with respect to the thickness parameter t.     

The Finite Difference Based Fast Adaptive Composite Grid Method

The fast adaptive composite grid (FAC) method is an iterative method for solving discrete boundary value problems on composite grids. McCormick introduced the method in [8] and considered the convergence behaviour for discrete problems resulting from finite volume element discretization on composite grids. In this paper we consider discrete problems resulting from finite difference discretization on composite grids. We distinguish between two obvious discretization approaches at the grid points on the interfaces between fine and coarse subgrids. The FAC method for solving such discrete problems is described. In the FAC method several intergrid transfer operators appear. We study how the convergence behaviour depends on these intergrid transfer operators. Based on theoretical insights, (quasi-)optimal intergrid transfer operators are derived. Numerical results illustrate the fast convergence of the FAC method using these intergrid transfer operators.     

Multiple grid and Osher's scheme for the efficient solution of the steady Euler equations

An iterative method is developed for the solution of the steady Euler equations for inviscid flow. The system of hyperbolic conservation laws is discretized by a finite-volume Osher-discretization. The iterative method is a multiple grid (FAS) iteration with symmetric Gauss-Seidel (SGS) as a relaxation method. Initial estimates are obtained by full multigrid (FMG). In the pointwise relaxation the equations are kept in block-coupled form and local linearization of the equations and the boundary conditions is considered. The efficient formulation of Osher's discretization of the 2-D non-isentropic steady Euler equations and its linearization is presented. The efficiency of FAS-SGS iteration is shown for a transonic model problem. It appears that, for the problem considered, the rate of convergence is almost independent of the gridsize and that for all meshsizes the discrete system is solved up to truncation error accuracy in only a few (2 or 3) iteration cycles.     

Numerical study of linear and nonlinear string vibrations by means of physical discretization

To solve partial differential equations numerically, discretization of the continuous model is required and may be achieved either mathematically or physically. This paper illustrates how physical discretization of a continuous string may be accomplished by employing discrete model theory which has as its essential substance Newtonian mechanics.Typical examples of wave motion in discretized ‘linear’ and ‘non-linear’ strings are discussed. They include the transverse vibrations of a string after having been subjected to a given initial displacement, reflection and superposition of wave pulses in the string, and resonance of the string when coupled to a harmonic vibrator. The equations that arise after application of discrete model theory to these problems, describe the subsequent motion of the string, and are solved numerically by computer. In all cases the results obtained for the discrete linear string agree remarkably well with those for the corresponding continuous physical string. The stability of the solutions obtained by discretization are also investigated.     

A fully discrete scheme for diffusive-dispersive conservation laws

Summary.   We introduce a fully discrete (in both space and time) scheme for the numerical approximation of diffusive-dispersive hyperbolic conservation laws in one-space dimension. This scheme extends an approach by LeFloch and Rohde [4]: it satisfies a cell entropy inequality and, as a consequence, the space integral of the entropy is a decreasing function of time. This is an important stability property, shared by the continuous model as well. Following Hayes and LeFloch [2], we show that the limiting solutions generated by the scheme need not coincide with the classical Oleinik-Kruzkov entropy solutions, but contain nonclassical undercompressive shock waves. Investigating the properties of the scheme, we stress various similarities and differences between the continuous model and the discrete scheme (dynamics of nonclassical shocks, nucleation, etc). Received November 15, 1999 / Revised version received May 27, 2000 / Published online March 20, 2001     

Entropic order quantity (EnOQ) model for deteriorating items

Jaber et al. [M.Y. Jaber, R.Y. Nuwayhid, M.A. Rosen, Price-driven economic order systems from a thermodynamic point of view, Int. J. Prod. Res. 42 (24) (2004) 5167–5184] suggested that it might be possible to improve production systems performance by applying the first and second laws of thermodynamics to reduce system entropy (or disorder). They then used these laws to modify the economic order quantity (EOQ) model to derive an equivalent entropic order quantity (EnOQ). The results suggested that larger quantities should be ordered than is suggested by the classical EOQ model.     

A FULLY DISCRETE IMPLICIT-EXPLICIT FINITE ELEMENT METHOD FOR SOLVING THE FITZHUGH-NAGUMO MODEL

This work develops a fully discrete implicit-explicit finite element scheme for a parabolic-ordinary system with a nonlinear reaction term which is known as the FitzHugh-Nagumo model from physiology. The first-order backward Euler discretization for the time derivative, and an implicit-explicit discretization for the nonlinear reaction term are employed for the model, with a simple linearization technique used to make the process of solving equations more efficient. The stability and convergence of the fully discrete implicit-explicit finite element method are proved, which shows that the FitzHugh-Nagumo model is accurately solved and the trajectory of potential transmission is obtained. The numerical results are also reported to verify the convergence results and the stability of the proposed method.     

Unbiased multi-index Monte Carlo

We introduce a new class of Monte Carlo-based approximations of expectations of random variables such that their laws are only available via certain discretizations. Sampling from the discretized versions of these laws can typically introduce a bias. In this paper, we show how to remove that bias, by introducing a new version of multi-index Monte Carlo (MIMC) that has the added advantage of reducing the computational effort, relative to i.i.d. sampling from the most precise discretization, for a given level of error. We cover extensions of results regarding variance and optimality criteria for the new approach. We apply the methodology to the problem of computing an unbiased mollified version of the solution of a partial differential equation with random coefficients. A second application concerns the Bayesian inference (the smoothing problem) of an infinite-dimensional signal modeled by the solution of a stochastic partial differential equation that is observed on a discrete space grid and at discrete times. Both applications are complemented by numerical simulations.     

Semi-Discrete and Fully Discrete Weak Galerkin Finite Element Methods for a Quasistatic Maxwell Viscoelastic Model

This paper considers weak Galerkin finite element approximations on polygonal/polyhedral meshes for a quasistatic Maxwell viscoelastic model. The spatial discretization uses piecewise polynomials of degree $k (k ≥ 1)$ for the stress approximation, degree $k+1$ for the velocity approximation, and degree $k$ for the numerical trace of velocity on the inter-element boundaries. The temporal discretization in the fully discrete method adopts a backward Euler difference scheme. We show the existence and uniqueness of the semi-discrete and fully discrete solutions, and derive optimal a priori error estimates. Numerical examples are provided to support the theoretical analysis.     

Discrete May–Leonard Competition Models III

A discretization method attributed to Kahan is used to approximate the May–Leonard (M–L) competition model for three species. The local dynamics of this discrete-time M–L model are analyzed. This model differs from the discrete M–L models being studied previously. This discrete model shows dynamical consistency with the continuous M–L model. Numerically, we showed that the discrete M–L model has a degenerate Hopf-bifurcation, which is consistent with the continuous M–L model.     

Well-balanced unstaggered central schemes for one and two-dimensional shallow water equation systems

We propose a new well-balanced unstaggered central finite volume scheme for hyperbolic balance laws with geometrical source terms. In particular we construct a new one and two-dimensional finite volume method for the numerical solution of shallow water equations on flat/variable bottom topographies. The proposed scheme evolves a non-oscillatory numerical solution on a single grid, avoids the time consuming process of solving Riemann problems arising at the cell interfaces, and is second-order accurate both in space and time. Furthermore, the numerical scheme follows a well-balanced discretization that first discretizes the geometrical source term according to the discretization of the flux terms, and then mimics the surface gradient method and discretizes the water height according to the discretization of the water level. The resulting scheme exactly satisfies the C-property at the discrete level. The proposed scheme is then applied and classical one and two-dimensional shallow water equation problems with flat or variable bottom topographies are successfully solved. The obtained numerical results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potential and efficiency of the proposed method.     

High order accurate semi-implicit WENO schemes for hyperbolic balance laws

In this paper we propose a family of well-balanced semi-implicit numerical schemes for hyperbolic conservation and balance laws. The basic idea of the proposed schemes lies in the combination of the finite volume WENO discretization with Roe’s solver and the strong stability preserving (SSP) time integration methods, which ensure the stability properties of the considered schemes [S. Gottlieb, C.-W. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM Rev. 43 (2001) 89-112]. While standard WENO schemes typically use explicit time integration methods, in this paper we are combining WENO spatial discretization with optimal SSP singly diagonally implicit (SDIRK) methods developed in [L. Ferracina, M.N. Spijker, Strong stability of singly diagonally implicit Runge-Kutta methods, Appl. Numer. Math. 58 (2008) 1675-1686]. In this way the implicit WENO numerical schemes are obtained. In order to reduce the computational effort, the implicit part of the numerical scheme is linearized in time by taking into account the complete WENO reconstruction procedure. With the proposed linearization the new semi-implicit finite volume WENO schemes are designed.A detailed numerical investigation of the proposed numerical schemes is presented in the paper. More precisely, schemes are tested on one-dimensional linear scalar equation and on non-linear conservation law systems. Furthermore, well-balanced semi-implicit WENO schemes for balance laws with geometrical source terms are defined. Such schemes are then applied to the open channel flow equations. We prove that the defined numerical schemes maintain steady state solution of still water. The application of the new schemes to different open channel flow examples is shown.     

Integrable Discrete Model for One‐Dimensional Soil Water Infiltration

We propose an integrable discrete model of one‐dimensional soil water infiltration. This model is based on the continuum model by Broadbridge and White, which takes the form of nonlinear convection–diffusion equation with a nonlinear flux boundary condition at the surface. It is transformed to the Burgers equation with a time‐dependent flux term by the hodograph transformation. We construct a discrete model preserving the underlying integrability, which is formulated as the self‐adaptive moving mesh scheme. The discretization is based on linearizability of the Burgers equation to the linear diffusion equation, but the naïve discretization based on the Euler scheme which is often used in the theory of discrete integrable systems does not necessarily give a good numerical scheme. Taking desirable properties of a numerical scheme into account, we propose an alternative discrete model that produces solutions with similar accuracy to direct computation on the original nonlinear equation, but with clear benefits regarding computational cost.     

An unconditionally stable uncoupled scheme for a triphasic Cahn–Hilliard/Navier–Stokes model

We propose an original scheme for the time discretization of a triphasic Cahn–Hilliard/Navier–Stokes model. This scheme allows an uncoupled resolution of the discrete Cahn–Hilliard and Navier‐Stokes system, which is unconditionally stable and preserves, at the discrete level, the main properties of the continuous model. The existence of discrete solutions is proved, and a convergence study is performed in the case where the densities of the three phases are the same. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq. 2013