A weakly compressible MPS method for modeling of open‐boundary free‐surface flow

A mesh‐free particle method, based on the moving particle semi‐implicit (MPS) interaction model, has been developed for the simulation of two‐dimensional open‐boundary free‐surface flows. The incompressibility model in the original MPS has been replaced with a weakly incompressible model. The effect of this replacement on the efficiency and accuracy of the model has been investigated. The new inflow–outflow boundary conditions along with the particle recycling strategy proposed in this study extend the application of the model to open‐boundary problems. The final model is able to simulate open‐boundary free surface flow in cases of large deformation and fragmentation of free surface. The models and proposed algorithms have been validated and applied to sample problems. The results confirm the model's efficiency and accuracy. Copyright © 2009 John Wiley & Sons, Ltd.     

超短超强激光脉冲在空气中长距离传输的研究进展

 介绍了超短超强激光脉冲在空气中长距离传输的最新研究进展,并分析了目前该项研究的几种理论模型,重点介绍了动态空间补给模型所给出的耦合方程及其研究结果,分析了不同情况下长距离传输模型的适用性。     

On the Infinitesimal Rigidity of Homogeneous Varieties

Let XP be a variety (respectively an open subset of an analytic submanifold) and let xX be a point where all integer valued differential invariants are locally constant. We show that if the projective second fundamental form of X at x is isomorphic to the second fundamental form of a point of a Segre P× P, n,m2, a Grassmaniann G(2,n+2), n4, or the Cayley plane OP², then X is the corresponding homogeneous variety (resp. an open subset of the corresponding homogeneous variety). The case of the Segre P²×P² had been conjectured by Griffiths and Harris in [GH]. If the projective second fundamental form of X at x is isomorphic to the second fundamental form of a point of a Veronese v₂(P) and the Fubini cubic form of X at x is zero, then X=v₂ (P) (resp. an open subset of v₂(P)). All these results are valid in the real or complex analytic categories and locally in the C category if one assumes the hypotheses hold in a neighborhood of any point x. As a byproduct, we show that the systems of quadrics I₂(P P) S²C, I₂(P¹× P) S²C and I₂(S₅) S²C¹⁶ are stable in the sense that if A S^* is an analytic family such that for t0,AA, then A₀A. We also make some observations related to the Fulton–:Hansen connectedness theorem.     

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.     

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     

钝体后湍流预混燃烧的PDF模拟

本文采用PDF方法对矩形燃烧室内钝体后的湍流预混火焰进行了数值模拟。脉动速度-频率-标量联合的PDF输运方程用Monte Carlo方法求解,质量、动量和能量的平均值由基于无结构网格的有限体积法求解,压力通过状态方程获得。PDF方程中所需的平均密度、平均速度和压力由有限体积法提供,并将用Monte Carlo方法求出的雷诺应力、化学反应源项和比热比传递给有限体积法。本文对丙烷和空气燃烧的不同简化化学反应机理进行了研究,并与实验结果进行比较,获得满意的结果。     

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.     

A non-conforming finite element method with anisotropic mesh grading for the Stokes problem in domains with edges

The solution of the Stokes problem in three-dimensional domainswith edges has anisotropic singular behaviour which is treatednumerically by using anisotropic finite element meshes. Thevelocity is approximated by Crouzeix–Raviart (nonconformingP₁ ) elements and the pressure by piecewise constants. Thismethod is stable for general meshes (without minimal or maximalangle condition). Denoting by N_e the number of elements in themesh, the interpolation and consistency errors are of the optimalorder h N_e^–1/3 which is proved for tensor product meshes.As a by-product, we analyse also nonconforming prismatic elementswith P₁ [oplus ] span {x₃²} as the local space for the velocitywhere x₃ is the direction of the edge.     

A posteriori L2 error estimation on anisotropic tetrahedral finite element meshes

A new a posteriori L₂ norm error estimator is proposed for thePoisson equation. The error estimator can be applied to anisotropictetrahedral or triangular finite element meshes. The estimatoris rigorously analysed for Dirichlet and Neumann boundary conditions. The lower error bound relies on specifically designed anisotropicbubble functions and the corresponding inverse inequalities.The upper error bound utilizes non-standard anisotropic interpolationestimates. Its proof requires H² regularity of the Poisson problem,and its quality depends on how good the anisotropic mesh resolvesthe anisotropy of the problem. This is measured by a so-called‘matching function’. A numerical example supports the anisotropic error analysis.     

NUMERICAL SOLUTION OF THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON THE CURVILINEAR HALF-STAGGERED MESH NUMERICAL SOLUTION OF THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON THE CURVILINEAR HALF-STAGGERED MESH

1. IntroductionLet us consider the unsteady incompressible Navier--Stokes equatinns (qSE)on a twedimensional reginn n with boundary On. Here w is the velocity vector; interms of its Caxtesian,components w = (u, v)"; p is the Pressure. The initial conditionis given assatisheng (1.2). We are concerned mainly with the solid wall boundary conditiollbut we will also briefly discuss the outflow boundary conditinn. Note the convectionterm can ajBo be written in conservative form (w' glad)w = div(…