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     

Harmonic Analysis on the Compact Homogeneous Spaces──about Fourier Multiplier ofH~P（M） Space

ＨａｒｍｏｎｉｃＡｎａｌｙｓｉｓｏｎｔｈｅＣｏｍｐａｃｔＨｏｍｏｇｅｎｅｏｕｓＳｐａｃｅｓ──ａｂｏｕｔＦｏｕｒｉｅｒＭｕｌｔｉｐｌｉｅｒｏｆＨ~Ｐ（Ｍ）Ｓｐａｃｅ￥ＤｏｎｇＤａｏｚｈｅｎ；ＺｈｅｎｇＸｕｅａｎ（ＨｅｎａｎＵｎｉｖｅｒｓｉｔｙ，Ｋａｉ?..     

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

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

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.     

Carathéodory's method for a class of second order differential equations on the half-line

The method of equivalent variational methods, originally due to Carathéodory for free problems in the calculus of variations is extended to investigate boundary value problems for a class of second order differential equations on the half-line. Some applications are presented to illustrate the potential of this method.     

ON DIRECT METHOD OF SOLUTION FOR A CLASS OF SINGULAR INTEGRAL EQUATIONS

In this article, by introducing characteristic singular integral operator and associate singular integral equations (SIEs), the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.     

Light Scattering by Cylindrical Fibers with High Aspect Ratio Using the Null‐Field Method with Discrete Sources

The problem of computing light scattering by cylindrical fibers with high aspect ratio in the framework of the Null‐Field method with discrete sources is treated. Numerical experiments for investigating the scattering properties of two fiber geometries are performed using distributed spherical vector wave functions as discrete sources.     

关于有限Abel群直接因子的平均数

对于有限Ａｂｅｌ群Ｇ，ｔ（Ｇ）定义为其直接因子的个数，在［５］中，Ｅ．Ｋｒａｔｚｅｌ对于和式得到一渐近公式，对应余项为△＿１（ｘ）．本文利用Ｅ．Ｆｏｕｖｒｙ和Ｈ　．Ｉｗａｎｉｅｃ的方法，得到了关于△＿１（ｘ）的一个新的估计，改进了Ｈ．Ｍｅｎｚｅｒ等人的结果．