Computation of complex fluid flows using an adaptive grid method

Recently the concept of adaptive grid computation has received much attention in the computational fluid dynamics research community. This paper continues the previous efforts of multiple one-dimensional procedures in developing and asessing the ideas of adaptive grid computation. The focus points here are the issue of numerical stability induced by the grid distribution and the accuracy comparison with previously reported work. Two two-dimensional problems with complicated characteristics—namely, flow in a channel with a sudden expansion and natural convection in an enclosed square cavity—are used to demonstrate some salient features of the adaptive grid method. For the channel flow, by appropriate distribution of the grid points the numerical algorithm can more effectively dampen out the instabilities, especially those related to artificial boundary treatments, and hence can converge to a steady-state solution more rapidly. For a more accurate finite difference operator, which contains less undesirable numerical diffusion, the present adaptive grid method can yield a steady-state and convergent solution, while uniform grids produce non-convergent and numerically oscillating solutions. Furthermore, the grid distribution resulting from the adaptive procedure is very responsive to the different characteristics of laminar and turbulent flows. For the problem of natural convection, a combination of a multiple one-dimensional adaptive procedure and a variational formulation is found very useful. Comparisons of the solutions on uniform and adaptive grids with the reported benchmark calculations demonstrate the important role that the adaptive grid computation can play in resolving complicated flow characteristics.     

Dynamic stability of a curved rail under a moving load

The stability of a curved rail under a constant moving load has been investigated using a linear theory; critical speeds of the moving load, and the dynamic rail deflections and rotation were calculated. The effect of the foundation was included through distributed linear springs. It was assumed that the moving load remains in constant contact with the rail and travels along a fixed path on the rail head.     

Conformational analysis of methylphenidate: comparison of molecular orbital and molecular mechanics methods

Summary Methylphenidate (MP) binds to the cocaine binding site on the dopamine transporter and inhibits reuptake of dopamine, but does not appear to have the same abuse potential as cocaine. This study, part of a comprehensive effort to identify a drug treatment for cocaine abuse, investigates the effect of choice of calculation technique and of solvent model on the conformational potential energy surface (PES) of MP and a rigid methylphenidate (RMP) analogue which exhibits the same dopamine transporter binding affinity as MP. Conformational analysis was carried out by the AM1 and AM1/SM5.4 semiempirical molecular orbital methods, a molecular mechanics method (Tripos force field with the dielectric set equal to that of vacuum or water) and the HF/6-31G* molecular orbital method in vacuum phase. Although all three methods differ somewhat in the local details of the PES, the general trends are the same for neutral and protonated MP. In vacuum phase, protonation has a distinctive effect in decreasing the regions of space available to the local conformational minima. Solvent has little effect on the PES of the neutral molecule and tends to stabilize the protonated species. The random search (RS) conformational analysis technique using the Tripos force field was found to be capable of locating the minima found by the molecular orbital methods using systematic grid search. This suggests that the RS/Tripos force field/vacuum phase protocol is a reasonable choice for locating the local minima of MP. However, the Tripos force field gave significantly larger phenyl ring rotational barriers than the molecular orbital methods for MP and RMP. For both the neutral and protonated cases, all three methods found the phenyl ring rotational barriers for the RMP conformers/invertamers (denoted as cte, tte, and cta) to be: cte, tte> MP > cta. Solvation has negligible effect on the phenyl ring rotational barrier of RMP. The B3LYP/6-31G* density functional method was used to calculate the phenyl ring rotational barrier for neutral MP and gave results very similar to those of the HF/6-31G* method.     

Tail Calculus with Remainder,Applications to Tail Expansions for Infinite Order Moving Averages,Randomly Stopped Sums,and Related Topics

We derive asymptotic expansions for tails of infinite weighted convolutions of some heavy-tailed distributions. Applications are given to tail expansion of the marginal distribution of ARMA processes, randomly stopped sums, as well as limiting waiting time distribution. AMS 2000 Subject Classifications. Primary—62E99, Secondary—41A60, 44A35, 60G50, 60K25     

Smith-Purcell Radiation in the Sub-mm Wavelength Region

Spontaneous Smith-Purcell radiation at wavelengths of 100 and 220m has been observed. The angles of observation were 17° and 28° and the electron beam energies 1.9 and 1.56MeV, respectively. The radiated power was about 320nW at 17° and 92nW at 28°. These values are in reasonable agreement with what is theoretically expected if the emission process is due to the acceleration of image charges induced on the grating surface by the electrons of the beam.     

Mixed finite volume methods on nonstaggered quadrilateral grids for elliptic problems

We construct and analyze a mixed finite volume method on quadrilateral grids for elliptic problems written as a system of two first order PDEs in the state variable (e.g., pressure) and its flux (e.g., Darcy velocity). An important point is that no staggered grids or covolumes are used to stabilize the system. Only a single primary grid system is adopted, and the degrees of freedom are imposed on the interfaces. The approximate flux is sought in the lowest-order Raviart-Thomas space and the pressure field in the rotated- nonconforming space. Furthermore, we demonstrate that the present finite volume method can be interpreted as a rotated- nonconforming finite element method for the pressure with a simple local recovery of flux. Numerical results are presented for a variety of problems which confirm the usefulness and effectiveness of the method.

     

网络条件下列车追踪模型及延迟传播的研究

描述了一种移动闭塞下考虑不同类型列车混跑时基于元胞自动机的列车追踪模型.将该模型应用到铁路网，模拟了网络条件下列车运行情况.应用该模型模拟了移动闭塞下列车延迟传播的现象.分析了发车间隔松弛时间、初始延迟时间等因素对列车延迟的影响.将该模型的模拟结果和理论公式的计算结果进行了比较，验证了模型的可靠性和有效性. 关键词： 元胞自动机 铁路网络 列车追踪模型 移动闭塞     

一种基于正负差图像的运动目标检测新方法

运动目标检测领域中现有的差图像法是利用绝对值差图像检测差图像上运动目标区域，用现有方法检测时易受噪声干扰，而且当摄像机有自运动时需要进行背景运动补偿。因此，提出一种新算法，即首先分别计算正差图像与负差图像，然后利用运动目标区域在正差图像与负差图像中的幅值、形状以及运动等信息的对称性对其进行检测，最后给出针对飞机尾焰序列图像进行检测的结果。实验结果表明：该方法可提高运动目标检测的可靠性与效率。     

非结构网格上温度扩散方程的能流计算方法

讨论非结构网格上温度扩散方程的能流计算方法.应用有限点方法(Finite Point Method,简称FPM)导出基于有限点两点公式和三点公式的能流计算公式,该公式适用于任意多边形及非匹配网格等非结构网格;给出网格角点温度新的计算公式.数值试验表明:基于两点公式的离散解和基于三点公式的离散解均具有平方阶的收敛速度;基于三点公式的离散解的精度总优于基于两点公式的离散解.     

A local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.