A possible counterexample to well posedness of entropy solutions and to Godunov scheme convergence

A particular case of initial data for the two-dimensional Euler equations is studied numerically. The results show that the Godunov method does not always converge to the physical solution, at least not on feasible grids. Moreover, they suggest that entropy solutions (in the weak entropy inequality sense) are not well posed.

     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Instabilities in the Chapman-Enskog Expansion and Hyperbolic Burnett Equations

It is well-known that the classical Chapman-Enskog procedure does not work at the level of Burnett equations (the next step after the Navier-Stokes equations). Roughly speaking, the reason is that the solutions of higher equations of hydrodynamics (Burnett's, etc.) become unstable with respect to short-wave perturbations. This problem was recently attacked by several authors who proposed different ways to deal with it. We present in this paper one of possible alternatives. First we deduce a criterion for hyperbolicity of Burnett equations for the general molecular model and show that this criterion is not fulfilled in most typical cases. Then we discuss in more detail the problem of truncation of the Chapman-Enskog expansion and show that the way of truncation is not unique. The general idea of changes of coordinates (based on analogy with the theory of dynamical systems) leads finally to nonlinear Hyperbolic Burnett Equations (HBEs) without using any information beyond the classical Burnett equations. It is proved that HBEs satisfy the linearized H-theorem. The linear version of the problem is studied in more detail, the complete Chapman-Enskog expansion is given for the linear case. A simplified proof of the Slemrod identity for Burnett coefficients is also given.     

Comparisons between model equations for long waves

Summary Considered here are model equations for weakly nonlinear and dispersive long waves, which feature general forms of dispersion and pure power nonlinearity. Two variants of such equations are introduced, one of Korteweg-de Vries type and one of regularized long-wave type. It is proven that solutions of the pure initial-value problem for these two types of model equations are the same, to within the order of accuracy attributable to either, on the long time scale during which nonlinear and dispersive effects may accumulate to make an order-one relative difference to the wave profiles.This research was supported in part by the National Science Foundation. A considerable portion of the project was completed while the first author was resident at the Institute for Mathematics and Its Applications, University of Minnesota.     

Parametric oscillation at 2.1 μm in periodically poled RbTiOAsO4

An Nd:YAG laser-pumped periodically poled RbTiOAsO₄ (PPRTA) optical parametric oscillator (OPO) near degeneracy is reported. At a pump power of 3.7 W, a net conversion efficiency of 43% was obtained. These crystals were found to be phase-matchable for the higher order quasi-phase-matched second-harmonic generation (SHG) in the visible range. Detailed analyses of the phase-matching properties for these processes were made by using the published Sellmeier and thermo-optic dispersion formulas.     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

解奇异非光滑方程组的牛顿法和不精确牛顿法

本文主要解决奇异非光滑方程组的解法。应用一种新的次微分的外逆，我们提出了牛顿法和不精确牛顿法，它们的收敛性同时也得到了证明。这种方法能更容易在一引起实际应用中实现。这种方法可以看作是已存在的解非光滑方程组的方法的延伸。     

从立方抛物线谈起(1)--余弦型振动,椭圆余弦型振动与双曲线型非周期运动

从立方抛物线的特性谈起 ,用较初浅的方法 ,借助于雅可比椭圆函数求椭圆方程的解来说明一大类一维保守系统的余弦振动、椭圆余弦型振动与双曲线型非周期运动 .并以杜芬振子、单摆、倒摆、转动圆环上的小环运动为典型实例作了详细讨论