Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Nondensely defined evolution impulsive differential inclusions with nonlocal conditions

In this paper, we shall establish sufficient conditions for the existence of integral solutions for some nondensely defined evolution impulsive differential inclusions in Banach spaces with nonlocal conditions.     

超重核266Hs形成截面的入射道依赖

计算和比较了²⁶Mg+²⁴⁴Cm, ²⁷Al+²⁴³Am和³²S+²³⁸U3个反应系统的俘获截面和复合核²⁷⁰Hs形成截面. 在俘获截面计算中, 考虑了靶核形变效应. 穿越库仑势垒后, 反应系统由熔合谷进入不对称裂变谷. 只有越过不对称裂变谷中的条件鞍点的事件才进入复合核组态. 我们用考虑中子流动和径向运动的二参量Smoluchowski扩散方程来处理中间阶段的动力学过程. 此外, 还计算了经4n蒸发形成超重核²⁶⁶Hs的截面. 研究表明, 入射道的势垒分布, 中间阶段的条件鞍点高度对俘获截面, 复合核形成几率, 以及最终的超重核形成截面有显著影响.     

多元线性回归应用于近代物理实验

本文介绍非线性方程转化直线性方程和多元线性回归法来分析近代物理实验中塞曼效应分裂干涉圆环多处选点测量的处理过程。     

A marker method for the solution of the damped Burgers' equation

A new method for the solution of the damped Burgers' equation is described. The marker method relies on the definition of a convective field associated with the underlying partial differential equation; the information about the approximate solution is associated with the response of an ensemble of markers to this convective field. Some key aspects of the method, such as the selection of the shape function and the initial loading, are discussed in some details. The marker method is applicable to a general class of nonlinear dispersive partial differential equations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Velocity spectrum for non-Markovian Brownian motion in a periodic potential

Non-Markovian Brownian motion in a periodic potential is studied by means of an electronic analogue simulator. Velocity spectra, the Fourier transforms of velocity autocorrelation functions, are obtained for three types of random force, that is, a white noise, an Ornstein—Uhlenbeck process, and a quasimonochromatic noise. The analogue results are in good agreement both with theoretical ones calculated with the use of a matrix-continued-fraction method, and with the results of digital simulations. An unexpected extra peak in the velocity spectrum is observed for Ornstein-Uhlenbeck noise with large correlation time. The peak is attributed to a slow oscillatory motion of the Brownian particle as it moves back and forth over several lattice spaces. Its relationship to an approximate Langevin equation is discussed.     

拟线性抛物型方程广义解的局部有界性

本文给出了一类较广泛的拟线性抛物型方程广义解的局部有界性,进而也得到了解的Holder连续性.     

动态模糊随机信息处理的数学方法

本文系统地概述了我们近年来提出的动态模糊随机信息处理的数学方法,内容包括模糊随机变量、模糊随机过程和模糊随机微分方程的基本解法等方面的基本概念、基本定义和某些重要的定理,以及动态系统的模糊随机响应与可靠性分析的方法等。这些方法是为我们研究工程实际问题的需要逐步发展起来的,对于处理某些类型的问题简便实用。