差分流线扩散法的局域稳定性分析

本文讨论差分流线扩散法的局部稳定性分析，在理论上证明了该方法依旧保持着流线扩散法扰动局部影响的特性。     

标的股价服从混合过程的期权定价公式及有限元算法

本文将马尔科夫跳跃过程叠加于 Ito过程 ,形成混合过程 ,并用该过程来刻画股价走势情况。而后在标的股价服从混合过程的基础上 ,推导出了欧式看涨期权的定价公式 ,并对美式看跌期权定价给出了有限元算法。     

Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.     

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     

Approximation of Matrices and a Family of Gander Methods for Polar Decomposition

Two matrix approximation problems are considered: approximation of a rectangular complex matrix by subunitary matrices with respect to unitarily invariant norms and a minimal rank approximation with respect to the spectral norm. A characterization of a subunitary approximant of a square matrix with respect to the Schatten norms, given by Maher, is extended to the case of rectangular matrices and arbitrary unitarily invariant norms. Iterative methods, based on the family of Gander methods and on Higham’s scaled method for polar decomposition of a matrix, are proposed for computing subunitary and minimal rank approximants. Properties of Gander methods are investigated in details. AMS subject classification (2000) 65F30, 15A18     

用区域分解算法结合蒙特卡罗法求坦克温度场和红外辐射出射度

本文将区域分解算法和蒙特卡罗法相结合，求坦克温度场和红外辐射出射度．采用蒙特卡罗法计算辐射传递系数，可以考虑界面的复杂辐射特性，如：镜反射、各向异性发射、各向异性反射等；可直接考虑界面的复杂几何特性，如:相互遮挡、太阳入射方向上的投影面积等．引入辐射传递系数，分离了计算的难点，使得在时间域（计算步骤）上能把整个问题分解为若干个子问题并行处理、在空间域（计算区域）上将坦克分解为若干个子区域，缩小计算规模，并可使用多个处理器并行计算；同时减轻了蒙特卡罗法编程的难度，缩短了计算时间．     

Multi-cracks problem for piezoelectric materials strip subjected to dynamic loading

This paper analyses and models the dynamic interaction among permeable multi-cracks in a piezoelectric strip under anti-plane shear waves by the Schmidt method. The Fourier transform is applied and then two pairs of triple integral equations can be solved using the Schmidt method. The results show that the stress and the electric displacement intensity factors of cracks depend on not only the crack length and the piezoelectric coefficient, but also the thickness of the piezoelectric strip, the distance between multi-cracks and the frequency of incident wave.     

Optimized analysis of the critical behavior in polymer mixtures from Monte Carlo simulations

A complete outline is given for how to determine the critical properties of polymer mixtures with extrapolation methods similar to the Ferrenberg-Swendsen techniques recently devised for spin systems. By measuring not only averages but the whole distribution of the quantities of interest, it is possible to extrapolate the data obtained in only a few simulations nearT _c over the entire critical region, thereby saving at least 90% of the computer time normally needed to locate susceptibility peaks or cumulant intersections and still getting more precise results. A complete picture of the critical properties of polymer mixtures in the thermodynamic limit is then obtained with finite-size scaling functions. Since the amount of information extracted from a simulation in this way is drastically increased as compared to conventional methods, the investigation of mixtures with long chains or built-in asymmetries is now possible. As an example, the critical points, exponents, and amplitudes of dense, symmetric polymer mixtures with chain lengths ranging fromN=16 up toN=256 are determined within the framework of the 3D bond fluctuation model using grand canonical simulation techniques. As an example for an asymmetry, the generalization of the method to asymmetric monomer potentials is briefly discussed.     

水闸闸室结构三维样条边界元分析

用三维样边界元法分析水闸闸室结构。底板，闸墩和载水墙等为其子结构，交通桥，工作桥和胸墙等处理为内部支撑。地基和边载可以是任意的，只要能给定地表位移面力关系。在各种工况下，不论是设置平板门还是弧形门，是平底板是反拱底板，即使在稀疏剖分下也能给出高精度的位移场，应力场和地基反力场。     

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.