Attractors for Fully Discrete Finite Difference Scheme of Dissipative Zakharov Equations

A fully discrete finite difference scheme for dissipative Zakharov equations is analyzed. On the basis of a series of the time-uniform priori estimates of the difference solutions, the stability of the difference scheme and the error bounds of optimal order of the difference solutions are obtained in L²×H¹×H² over a finite time interval (0, T]. Finally, the existence of a global attractor is proved for a discrete dynamical system associated with the fully discrete finite difference scheme.     

A Convergent Lagrangian Discretization for a Nonlinear Fourth-Order Equation

A fully discrete Lagrangian scheme for numerical solution of the nonlinear fourth-order DLSS equation in one space dimension is analyzed. The discretization is based on the equation’s gradient flow structure in the \(L^2\)-Wasserstein metric. By construction, the discrete solutions are strictly positive and mass conserving. A further key property is that they dissipate both the Fisher information and the logarithmic entropy. Our main result is a proof of convergence of fully discrete to weak solutions in the limit of vanishing mesh size. Convergence is obtained for arbitrary nonnegative, possibly discontinuous initial data with finite entropy, without any CFL-type condition. The key estimates in the proof are derived from the dissipations of the two Lyapunov functionals. Numerical experiments illustrate the practicability of the scheme.     

Variational formulations for Vlasov–Poisson–Fokker–Planck systems

Time‐discrete variational schemes are introduced for both the Vlasov–Poisson–Fokker–Planck (VPFP) system and a natural regularization of the VPFP system. The time step in these variational schemes is governed by a certain Kantorovich functional (or scaled Wasserstein metric). The discrete variational schemes may be regarded as discretized versions of a gradient flow, or steepest descent, of the underlying free energy functionals for these systems. For the regularized VPFP system, convergence of the variational scheme is rigorously established. Copyright © 2000 John Wiley & Sons, Ltd.     

H1‐Galerkin expanded mixed finite element methods for nonlinear pseudo‐parabolic integro‐differential equations

H¹‐Galerkin mixed finite element method combined with expanded mixed element method is discussed for nonlinear pseudo‐parabolic integro‐differential equations. We conduct theoretical analysis to study the existence and uniqueness of numerical solutions to the discrete scheme. A priori error estimates are derived for the unknown function, gradient function, and flux. Numerical example is presented to illustrate the effectiveness of the proposed scheme. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A convergent difference scheme for the infinity Laplacian: Construction of absolutely minimizing Lipschitz extensions

This article considers the problem of building absolutely minimizing Lipschitz extensions to a given function. These extensions can be characterized as being the solution of a degenerate elliptic partial differential equation, the ``infinity Laplacian', for which there exist unique viscosity solutions.

A convergent difference scheme for the infinity Laplacian equation is introduced, which arises by minimizing the discrete Lipschitz constant of the solution at every grid point. Existence and uniqueness of solutions to the scheme is shown directly. Solutions are also shown to satisfy a discrete comparison principle.

Solutions are computed using an explicit iterative scheme which is equivalent to solving the parabolic version of the equation.

     

Study of a Generalized Nonlinear Euler-Poisson-Darboux System: Numerical and Bessel Based Solutions

In this paper a nonlinear Euler-Poisson-Darboux system is considered. In a first part, we proved the genericity of the hypergeometric functions in the development of exact solutions for such a systemin some special cases leading to Bessel type differential equations. Next, a finite difference scheme in two-dimensional case has been developed. The continuous system is transformed into an algebraic quasi linear discrete one leading to generalized Lyapunov-Sylvester operators. The discrete algebraic system is proved to be uniquely solvable, stable and convergent based on Lyapunov criterion of stability and Lax-Richtmyer equivalence theorem for the convergence. A numerical example has been provided at the end to illustrate the efficiency of the numerical scheme developed in section 3. The present method is thus proved to be more accurate than existing ones and lead to faster algorithms.     

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     

On the relation between the continuous and discrete Painlevé equations

A method for deriving difference equations (the discrete Painlevé equations in particular) from the Bäcklund transformations of the continuous Painlevé equations is discussed. This technique can be used to derive several of the known discrete painlevé equations (in particular, the first and second discrete Painlevé equations and some of their alternative versions). The Painlevé equations possess hierarchies of rational solutions and one-parameter families of solutions expressible in terms of the classical special functions for special values of the parameters. Hence, the aforementioned relations can be used to generate hierarchies of exact solutions for the associated discrete Painlevé equations. Exact solutions of the Painlevé equations simultaneously satisfy both a differential equation and a difference equation, analogously to the special functions.     

Scheduling arrivals to a stochastic service mechanism

We consider the problem of scheduling the arrivals of a fixed number of customers to a stochastic service mechanism to minimize an expected cost associated with operating the system. We consider the special case of exponentially distributed service times and the problems in general associated with obtaining exact analytic solutions. For general service time distributions we obtain approximate numerical solutions using a stochastic version of gradient search employing Infinitesimal Perturbation Analysis estimates of the objective function gradient obtained via simulation.     

Difference schemes on uniform grids performed by general discrete operators

Our aim is to set the foundations of a discrete vectorial calculus on uniform n-dimensional grids, that can be easily reformulated on general irregular grids. As a key tool we first introduce the notion of tangent space to any grid node. Then we define the concepts of vector field, field of matrices and inner products on the space of grid functions and on the space of vector fields, mimicking the continuous setting. This allows us to obtain the discrete analogous of the basic first order differential operators, gradient and divergence, whose composition define the fundamental second order difference operator. As an application, we show that all difference schemes, with constant coefficients, for first and second order differential operators with constant coefficients can be seen as difference operators of the form for suitable choices of q, and  . In addition, we characterize special properties of the difference scheme, such as consistency, symmetry and positivity in terms of q, and  .     

An existence result for a free boundary shallow water model using a Lagrangian scheme

We propose a free boundary shallow water model for which we give an existence theorem. The proof uses an original Lagrangian discrete scheme in order to build a sequence of approximate solutions. The properties of this scheme allow to treat the difficulties linked to the boundary motion. These approximate solutions verify some compactness results which allow us to pass to the limit in the discrete problem.     

Stability and convergence of two‐grid Crank‐Nicolson extrapolation scheme for the time‐dependent natural convection equations

In this paper, we consider the Crank‐Nicolson extrapolation scheme for the 2D/3D unsteady natural convection problem. Our numerical scheme includes the implicit Crank‐Nicolson scheme for linear terms and the recursive linear method for nonlinear terms. Standard Galerkin finite element method is used to approximate the spatial discretization. Stability and optimal error estimates are provided for the numerical solutions. Furthermore, a fully discrete two‐grid Crank‐Nicolson extrapolation scheme is developed, the corresponding stability and convergence results are derived for the approximate solutions. Comparison from aspects of the theoretical results and computational efficiency, the two‐grid Crank‐Nicolson extrapolation scheme has the same order as the one grid method for velocity and temperature in H¹‐norm and for pressure in L²‐norm. However, the two‐grid scheme involves much less work than one grid method. Finally, some numerical examples are provided to verify the established theoretical results and illustrate the performances of the developed numerical schemes.     

非线性耦合Schrdinger-KdV方程组的一个局部能量守恒格式

本文利用平均值离散梯度给出了一个构造哈密尔顿偏微分方程的局部能量守恒格式的系统方法.并用非线性耦合Schrdinger-KdV方程组加以说明.证明了格式满足离散的局部能量守恒律,在周期边界条件下,格式也保持离散整体能量及系统的其它两个不变量.最后数值实验验证了理论结果的正确性.     

On Discrete Painlevé Equations Associated with the Lattice KdV Systems and the Painlevé VI Equation

A new integrable nonautonomous nonlinear ordinary difference equation is presented that can be considered to be a discrete analogue of the Painlevé V equation. Its derivation is based on the similarity reduction on the two-dimensional lattice of integrable partial differential equations of Korteweg–de Vries (KdV) type. The new equation, which is referred to as generalized discrete Painlevé equation (GDP), contains various "discrete Painlevé equations" as subcases for special values/limits of the parameters, some of which have already been given in the literature. The general solution of the GDP can be expressed in terms of Painlevé VI (PVI) transcendents. In fact, continuous PVI emerges as the equation obeyed by the solutions of the discrete equation in terms of the lattice parameters rather than the lattice variables that label the lattice sites. We show that the bilinear form of PVI is embedded naturally in the lattice systems leading to the GDP. Further results include the establishment of Bäcklund and Schlesinger transformations for the GDP, the corresponding isomonodromic deformation problem, and the self-duality of its bilinear scheme.     

Some results on the stability and dynamics of finite difference approximations to nonlinear partial differential equations

A miscellany of results on the nonlinear instability and dynamics of finite difference discretizations of the Burgers and Kortweg de Vries equations is obtained using a variety of phase-plane, functional analytic, and regularity methods. For the semidiscrete (space-discrete, time-continuous) schemes, large-wave-numer instabilities occurring in special exact solutions are investigated, and parameter values for which the semidiscrete scheme is monotone are considered. For fully discrete schemes (space and time discrete), large-wave-number instabilities introduced by various time-stepping schemes such as forward Euler, leapfrog, and Runge–Kutta schemes are analyzed. Also, a time step restriction for the monotonicity of the forward-Euler time-stepping scheme, and regularity of a 4-stage monotone/conservative Runge–Kutta time stepping are investigated. The techniques used here may be employed, in conjunction with bifurcation-theoretic and weakly nonlinear analyses, to analyze the stability of numerical schemes for other nonlinear partial differential equations of both dissipative and dispersive varieties. © 1993 John Wiley & Sons, Inc.     

基于Maxwell模型的线性粘弹性波传导问题的半离散/全离散杂交应力有限元法

1引言设ΩСR²为一个凸的有界开集,边界为ЭΩ;T为一个正常数.我们考虑如下基于Maxwell模型的二维粘弹性固体介质波传导问题。     

Morrey global bounds and quasilinear Riccati type equations below the natural exponent

We obtain global bounds in Lorentz–Morrey spaces for gradients of solutions to a class of quasilinear elliptic equations with low integrability data. The results are then applied to obtain sharp existence results in the framework of Morrey spaces for Riccati type equations with a gradient source term having growths below the natural exponent of the operator involved. A special feature of our results is that they hold under a very general assumption on the nonlinear structure, and under a mild natural restriction on the boundary of the ground domain.     

Analysis of new conservative difference scheme for two-dimensional Rosenau-RLW equation

In this article, numerical solution for the Rosenau-RLW equation in 2D is considered and a conservative Crank–Nicolson finite difference scheme is proposed. Existence of the numerical solutions for the difference scheme has been shown by Browder fixed point theorem. A priori bound and uniqueness as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability and a second-order accuracy on both space and time of the difference scheme are proved. Numerical experiments are given to support our theoretical results.