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.     

TWO-POINT BOUNDARY VALUE PROBLEM FOR A CLASS OF SECOND-ORDER DIFFERENCE EQUATIONS

In this paper, we apply a critical point theorem and obtain the condition for the existence of three solutions to two-point boundary value problem of a second order nonlinear difference equation.     

Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

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.     

High-order accurate equations describing vibrations of thin bars

A method for deriving one-dimensional wave propagation equations in thin inhomogeneous anisotropic bars based on the mathematical homogenization theory for periodic media is used to obtain equations governing the longitudinal and transverse vibrations of a homogeneous circular bar. The equations are derived up to O(ε⁸) terms and take into account variable body forces and surface loads. Here, ε is the ratio of the bar’s typical thickness to the typical wavelength.     

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.     

两类方程亚纯解(n,1)级的上界估计

本文给出Riccati方程及另外一类具有代表性微分方程的亚纯解（n,1)级的上界估计，在一定条件下确立了文[2]中的猜测的正确性。     

Approximate Inverse Preconditioners for Some Large Dense Random Electrostatic Interaction Matrices

A sparse mesh-neighbour based approximate inverse preconditioner is proposed for a type of dense matrices whose entries come from the evaluation of a slowly decaying free space Green’s function at randomly placed points in a unit cell. By approximating distant potential fields originating at closely spaced sources in a certain way, the preconditioner is given properties similar to, or better than, those of a standard least squares approximate inverse preconditioner while its setup cost is only that of a diagonal block approximate inverse preconditioner. Numerical experiments on iterative solutions of linear systems with up to four million unknowns illustrate how the new preconditioner drastically outperforms standard approximate inverse preconditioners of otherwise similar construction, and especially so when the preconditioners are very sparse. AMS subject classification (2000) 65F10, 65R20, 65F35, 78A30