自由电子激光振荡器的二维数值模拟

在柱二维坐标系中,用数值计算方法求解了摇摆器内电子模拟方程组、光场方程组和无源区(不含摇摆器区)经非自适应方法变换后的无源旁轴波动方程并设计了相应的程序R_2D。用该程序计算的数值结果表明,在有源条件下二维数值模拟结果与高斯光束近似下的数值模拟结果符合较好;取完全相同的计算条件R-2D程序与洛斯阿拉莫斯实验室(LANL)的FELEX程序的计算结果符合较好;初步认为LANL实验数据的物理图象合理,数值结果基本可信;最后给出北京自由电子激光器(BFEL)的一些物理数据。     

A parallel inexact newton method for stochastic programs with recourse

A parallel inexact Newton method with a line search is proposed for two-stage quadratic stochastic programs with recourse. A lattice rule is used for the numerical evaluation of multi-dimensional integrals, and a parallel iterative method is used to solve the quadratic programming subproblems. Although the objective only has a locally Lipschitz gradient, global convergence and local superlinear convergence of the method are established. Furthermore, the method provides an error estimate which does not require much extra computation. The performance of the method is illustrated on a CM5 parallel computer.This work was supported by the Australian Research Council and the numerical experiments were done on the Sydney Regional Centre for Parallel Computing CM5.     

The many facets of linear programming

We examine the history of linear programming from computational, geometric, and complexity points of view, looking at simplex, ellipsoid, interior-point, and other methods. Received: June 22, 2000 / Accepted: April 4, 2001?Published online October 2, 2001     

曙光一号自由电子激光器的理论计算

系统总结曙光一号自由电子激光器理论计算的主要结果:包括曙光一号装置主要参数的选取和理解;磁场失谐曲线的计算;常参数摇摆器和变参数摇摆器的主要结果;高阶波导模的贡献;电子束参数扰动对激光性能的影响;空间电荷效应等。计算结果表明,常参数摇摆器激光输出功率可达80MW,效率约50％;变参数摇摆器激光输出功率可达250MW左右,效率约16％。     

Mixed methods on quadrilaterals and hexahedra

We describe a new family of discrete spaces suitable for use with mixed methods on certain quadrilateral and hexahedral meshes. The new spaces are natural in the sense of differential geometry, so all the usual mixed method theory, including the hybrid formulation, carries over to these new elements with proofs unchanged. Because transforming general quadrilaterals into squares introduces nonlinearity and because mixed methods involve the divergence operator, the new spaces are more complicated than either the corresponding Raviart-Thomas spaces for rectangles or corresponding finite element spaces for quadrilaterals. The new spaces are also limited to meshes obtained from a rectangular mesh through the application of a single global bilinear transformation. Despite this limitation, the new elements may be useful in certain topologically regular problems, where initially rectangular grids are deformed to match features of the physical region. They also illustrate the difficulties introduced into the theory of mixed methods by nonlinear transformations.     

Numerical solution of eigenvalue problems for Hamiltonian systems

This paper discusses the numerical solution of eigenvalue problems for Hamiltonian systems of ordinary differential equations. Two new codes are presented which incorporate the algorithms described here; to the best of the author’s knowledge, these are the first codes capable of solving numerically such general eigenvalue problems. One of these implements a new new method of solving a differential equation whose solution is a unitary matrix. Both codes are fully documented and are written inPfort-verifiedFortran 77, and will be available in netlib/aicm/sl11f and netlib/aicm/sl12f.     

Ranking Functions Induced by Weighted Average of Fuzzy Numbers

In this paper we present two definitions of possibilistic weighted average of fuzzy numbers, and by them we introduce two different rankings on the set of real fuzzy numbers. The two methods are dependent on several parameters. In the first case, the parameter is constant and the results generalize what Carlsson and Fuller have obtained in (2001). In the second case, the parameter is a function, not fixed a priori by the decision maker, but it depends on the position of the interval on the real axe. In all the two cases we call the parameter degree of risk, which takes into account of a risk-tendency or aversion of the decision maker.     

Numerical propagation of partially coherent laser beams through optical systems

The propagation of partially coherent beams through optical systems is computed numerically in one transverse dimension. The optical system is divided into different elementary segments, through which the propagation of light can be calculated by appropriate operators, working on the coherence function or the Wigner distribution function respectively. For the necessary changes between these two functions describing the partially coherent beams, the use of the remarkable z-transform seems to be an advantage. With this algorithm the grid and the resolution in the spatial frequency domain can be arbitrarily chosen in contrast to the usual Fourier transform, the influence of phase aberrations on the focusability of Gauss-Schell model beams is discussed as an application example of the numerical model. With the help of this tool, practical beam guiding systems can be simulated for use with multimode laser radiation.     

Cyclic reduction and FACR methods for piecewise hermite bicubic orthogonal spline collocation

Cyclic reduction and Fourier analysis-cyclic reduction (FACR) methods are presented for the solution of the linear systems which arise when orthogonal spline collocation with piecewise Hermite bicubics is applied to boundary value problems for certain separable partial differential equations on a rectangle. On anN×N uniform partition, the cyclic reduction and Fourier analysis-cyclic reduction methods requireO(N ²log₂ N) andO(N ²log₂log₂ N) arithmetic operations, respectively.