Energy Conservation with Non-Symplectic Methods: Examples and Counter-Examples

Energy conservation of numerical integrators is well understood for symplectic one-step methods. This article provides new insight into energy conservation with non-symplectic methods. Sufficient conditions and counter-examples are presented. AMS subject classification (2000) 65L06, 65P10, 37J99.Submitted June 2004. Accepted October 2004. Communicated by Syvert Nørsett.     

Phase dynamics in the complex Ginzburg-Landau equation

For αβ>−1, stable time periodic solutions A(X,T)=A_qe^iqX+iω_qT are the locally preferred planform for the complex Ginzburg-Landau equation

     

Riemann problem for kinematical conservation laws and geometrical features of nonlinear wavefronts

A pair of kinematical conservation laws (KCL) in a ray coordinatesystem (,t) are the basic equations governing the evolutionof a moving curve in two space dimensions. We first study elementarywave solutions and then the Riemann problem for KCL when themetric g, associated with the coordinate designating differentrays, is an arbitrary function of the velocity of propagationm of the moving curve. We assume that m>1 (m is appropriatelynormalized), for which the system of KCL becomes hyperbolic.We interpret the images of the elementary wave solutions inthe (,t)-plane to the (x,y)-plane as elementary shapes of themoving curve (or a nonlinear wavefront when interpreted in aphysical system) and then describe their geometrical properties.Solutions of the Riemann problem with different initial datagive the shapes of the nonlinear wavefront with different combinationsof elementary shapes. Finally, we study all possible interactionsof elementary shapes.     

BV estimates of Lax-Friedrichs' scheme for a class of nonlinear hyperbolic conservation laws

We give uniform BV estimates and -stability of Lax-Friedrichs' scheme for a class of systems of strictly hyperbolic conservation laws whose integral curves of the eigenvector fields are straight lines, i.e., Temple class, under the assumption of small total variation. This implies that the approximate solutions generated via the Lax-Friedrichs' scheme converge to the solution given by the method of vanishing viscosity or the Godunov scheme, and then the Glimm scheme or the wave front tracking method.

     

Entropy Production and Admissibility of Shocks

In shock wave theory there are two considerations in selecting the physically relevant shock waves.There is the admissibility criterion for the well-posedness of hyperbolic conservation laws.Another consideraztion concerns the entropy production across the shochs.The latter is natural from the physical point of view,but is not sufficient in its straightforward formulation,if the system is not genuinely nonlinear.In this paper we propose the principles of increasing entropy production and that of the superposition of shocks.These principles arc shown to be equivalent to the admissibility criterion.     

Equilibrium schemes for scalar conservation laws with stiff sources

We consider a simple model case of stiff source terms in hyperbolic conservation laws, namely, the case of scalar conservation laws with a zeroth order source with low regularity. It is well known that a direct treatment of the source term by finite volume schemes gives unsatisfactory results for both the reduced CFL condition and refined meshes required because of the lack of accuracy on equilibrium states. The source term should be taken into account in the upwinding and discretized at the nodes of the grid. In order to solve numerically the problem, we introduce a so-called equilibrium schemes with the properties that (i) the maximum principle holds true; (ii) discrete entropy inequalities are satisfied; (iii) steady state solutions of the problem are maintained. One of the difficulties in studying the convergence is that there are no estimates for this problem. We therefore introduce a kinetic interpretation of upwinding taking into account the source terms. Based on the kinetic formulation we give a new convergence proof that only uses property (ii) in order to ensure desired compactness framework for a family of approximate solutions and that relies on minimal assumptions. The computational efficiency of our equilibrium schemes is demonstrated by numerical tests that show that, in comparison with an usual upwind scheme, the corresponding equilibrium version is far more accurate. Furthermore, numerical computations show that equilibrium schemes enable us to treat efficiently the sources with singularities and oscillating coefficients.

     

Second class particles as microscopic characteristics in totally asymmetric nearest-neighbor -exclusion processes

We prove laws of large numbers for a second class particle in one-dimensional totally asymmetric -exclusion processes, under hydrodynamic Euler scaling. The assumption required is that initially the ambient particle configuration converges to a limiting profile. The macroscopic trajectories of second class particles are characteristics and shocks of the conservation law of the particle density. The proof uses a variational representation of a second class particle, to overcome the problem of lack of information about invariant distributions. But we cannot rule out the possibility that the flux function of the conservation law may be neither differentiable nor strictly concave. To give a complete picture we discuss the construction, uniqueness, and other properties of the weak solution that the particle density obeys.

     

非齐次守恒律方程分片光滑解的粘性方法

1　引　　言考虑非齐次守恒律方程ut+f(u) x =g(u) ,　　 -∞ 0 ,(1 .1 )u(x,0 ) =u0 (x) ,　　 -∞ 0 , (1 .5)g∈ C3且 g是 Lipschitz连续的 ,Lipschitz系数为 L . (1 .6 )对于一般守恒律齐次方程 ,粘性解逼近熵解的收敛阶为 O(ε ) [1 ] .在 f严格凸的条件下 ,其收敛速度可以提高到 O(ε|lnε|+ε) [2 ] ,[3] .本文考虑具有条件 (1 .5) (1 .6 )的非齐次方程(1 .1 ) ,在较广泛的一类初值条件下…