光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.

The accuracy of Kirchhoff's approximation in describing the far field speckles produced by random self-affine fractal surfaces

Based on the rigorous formulation of integral equations for the propagations of light waves at the medium interface, we carry out the numerical solutions of the random light field scattered from self-affine fractal surface samples. The light intensities produced by the same surface samples are also calculated in Kirchhoff's approximation, and their comparisons with the corresponding rigorous results show directly the degree of the accuracy of the approximation. It is indicated that Kirchhoff's approximation is of good accuracy for random surfaces with small roughness value w and large roughness exponent alpha. For random surfaces with larger w and smaller alpha, the approximation results in considerable errors, and detailed calculations show that the inaccuracy comes from the simplification that the transmitted light field is proportional to the incident field and from the neglect of light field derivative at the interface.

Approximate analytical description for fundamental-mode fields of graded-index fibers: Beyond the Gaussian approximation

An approximate analytical description for fundamental-mode fields of graded-index fibers is explicitly presented by use of the power-series expansion method, the maximum-value condition at the fiber axis, the decay properties of fundamental-mode fields at large distance from the fiber axis, and the approximate modal parameters U obtained from the Gaussian approximation. This analytical description is much more accurate than the Gaussian approximation and at the same time keep the simplicity of the latter. As two special examples, we present the approximate analytical formulas for the fundamental-mode fields of a step profile fiber and a Gaussian profile fiber, and we find that they are both highly accurate in the single-mode range by comparing them with the corresponding exact solutions.

Corrections to the paraxial approximation of an arbitrary free-propagation beam

A relatively simple transform from an arbitrary solution of the paraxial wave equation to the corresponding exact solution of the Helmholtz wave equation is derived in the condition that the evanescent waves are ignored and is used to study the corrections to the paraxial approximation of an arbitrary free-propagation beam. Specifically, the general lowest-order correction field is given in a very simple form and is proved to be exactly consistent with the perturbation method developed by Lax et nl. [Phys. Rev. A 11, 1365 (1975)]. Some special examples, such as the lowest-order correction to the paraxial approximation of a fundamental Gaussian beam whose waist plane has a parallel shin from the z = 0 plane, are presented. (C) 1998 Optical Society of America.

Corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media

The coupled differential recurrence equations for the corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media are established in terms of the perturbation method. All the corrections (including the longitudinal field corrections) to the paraxial approximation solutions are presented in the weak-guidance approximation. As a concrete application, the first-order longitudinal field correction and the second-order transverse field correction to the paraxial approximation of a Gaussian beam propagating in a transversely quadratic refractive index medium are analytically investigated. (C) 1999 Optical Society of America [S0740-3232(99)00310-5].