Evaluation of the drag force by integrating the energy dissipation rate in stokes flow for 2D domains using the FEM

The total drag force on the surface of a body, which is the sum of the form drag and the skin friction drag in a 2D domain, is numerically evaluated by integrating the energy dissipation rate in the whole domain for an incompressible Stokes fluid. The finite element method is used to calculate both the energy dissipation rate in the whole domain as well as the drag on the boundary of the body. The evaluation of the drag and the energy dissipation rate are post-processing operations which are carried out after the velocity field and the pressure field for the flow over a particular profile have been obtained. The results obtained for the flow over three different but constant area profiles—a circle, an ellipse and a cross-section of a prolate spheroid—with uniform inlet velocity are presented and it is shown that the total drag force times the velocity is equal to the total energy dissipation rate in the entire finite flow domain. Hence, by calculating the energy dissipation rate in the domain with unit velocity specified at the far-field boundary enclosing the domain, the drag force on the boundary of the body can be obtained.     

The Möbius distribution on the disc

A simple new family of distributions is proposed which has support the unit disc in two dimensions. The density functions of the family are unimodal, monotonic or uniantimodal. The bivariate symmetric beta distributions, which include the uniform distribution, are special cases, but many members of the family are skew. The distributions have three parameters, one controlling orientation, one controlling degree of concentration and the third controlling skewness, or more precisely off-centredness. Importantly, these parameters are globally orthogonal. An illustrative example of fitting the model to data is given. Conditional and marginal distributions are considered. The new distributions are compared favourably with an earlier suggestion of the same author.     

光子晶体的原理与应用

光子晶体是一种在微米、亚微米等光波长量级上折射率呈现周期性变化的介质材料，它使某些频率范围内的光子态密度大大降低．甚至完全形成光子禁带．本文介绍了光子晶体的原理、制备及应用．     

Predicting the onset of DNA supercoiling using a non-linear hemitropic elastic rod

Elastic rod models provide a means to interpret single molecule DNA experiments as well as predict DNA behavior under physiological conditions. Here we use an elastic rod model to predict the stability boundary (critical torque vs. applied tension) for single molecule DNA experiments in which the molecule is subjected to applied tension and twist. We discuss the shortcomings of the usual isotropic rod model. We then derive a consistent non-linear material law from the general representation for a hemitropic (chiral) rod. Finally, we present results of a standard bifurcation analysis predicting the stability boundary. We find results from the non-linear hemitropic rod to match the data closely.     

Space–time integrated least squares: a time‐marching approach

A time‐marching formulation is derived from the space–time integrated least squares (STILS) method for solving a pure hyperbolic convection equation and is numerically compared to various known methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Analytical formulas for the velocity field induced by an infinitely thin vortex ring

Three different analytical solutions are presented for a potential vortex ring using three different streamfunctions. Verification studies confirm that all three approaches are valid. It is found that the solution obtained using the Biot–Savart law is the most efficient method due to its simplicity. It is shown that all analytical results are accurate to within machine accuracy and sample calculations are included. Copyright © 2004 John Wiley & Sons, Ltd.     

高剥离类镁等电子序列的磁偶极M1和电四极E2光谱跃迁理论研究

利用全相对论性多组态Dirac-Fock广义平均能级方法，系统地计算了类镁离子3s3p磁偶极Ml^3P2--^3P1和电四极E2 ^2P2--^3P0(Z=20-103)光谱跃迁的能级间隔、跃迁几率和振子强度。计算中考虑了原子核的有限体积效应，进行了高阶Breit修正和QED修正，所得到的能级间隔和最近的实验数据及理论计算值进行了比较。计算结果表明：高原子序数的高荷电离子的磁偶极矩M1和电四极矩E2跃迁几率和中性原子的电偶极E1的相当，在ICF和MCF高温激光等离子体中，磁偶极矩M1和电四极矩E2跃迁过程不容被忽视。     

CO（ a）+ NO（X）传能体系中的新传能通道

在相对碰撞平动能为 0 .0 5eV的分子束实验条件下 ,研究了亚稳态CO(a) +NO(X)的E E传能通道 .通过测量和分析交碰区的传能发射光谱 ,在 780和 860nm处观测到了NO(b -a)跃迁Ogawa带的△v =+ 4和△v =+ 3序的发射光谱 .从而首次在实验上直接证实了传能过程中第四通道的存在 (CO (a) +NO(X)→CO (X) +NO(b) ) .这一通道的发现解释了前人测量到的在CO(a)与NO(X)碰撞传能过程中CO(a)的猝灭速率远大于NO(A ,B)生成速率的实验结果 ,并进一步证实了这一“经典”E E传能体系为电子交换机理的传能观点     

Measure of Ordering in Statistical Systems

In this paper, we introduce a phenomenological measure of ordering in statistical systems. Using this characteristic, we construct a formal model of a system evolving according to the law of steepest ascent as the control parameter increases.