Van der Waals状态方程可压缩多流体流动的PPM方法

基于流体体积分数的混合型多流体数值模型，将Piecewise Parabolic Method（PPM）方法应用于可压缩多流体流动的数值模拟，采用双波近似求解多流体van der Waals状态方程的Riemann问题．模拟高密度比且含有激波的可压缩多流体流动，典型的纯界面平移问题模拟结果表明，在接触间断的界面附近，压力和速度没有任何的振荡且界面数值耗散都被控制在2—3个网格之内；一维和二维算例表明，该数值方法可以有效地处理接触间断、激波和多维滑移线等物理问题，并能够比其它多流体数值方法更精细地模拟多流体交界面．     

High-order kinetic flux vector splitting schemes in general coordinates for ideal quantum gas dynamics

A class of high-order kinetic flux vector splitting schemes are presented for solving ideal quantum gas dynamics based on quantum statistical mechanics. The collisionless quantum Boltzmann equation approach is adopted and both Bose–Einstein and Fermi–Dirac gases are considered. The formulas for the split flux vectors are derived based on the general three-dimensional distribution function in velocity space and formulas for lower dimensions can be directly deduced. General curvilinear coordinates are introduced to treat practical problems with general geometry. High-order accurate schemes using weighted essentially non-oscillatory methods are implemented. The resulting high resolution kinetic flux splitting schemes are tested for 1D shock tube flows and shock wave diffraction by a 2D wedge and by a circular cylinder in ideal quantum gases. Excellent results have been obtained for all examples computed.     

Conceptual linearization of Euler governing equations to solve high speed compressible flow using a pressure‐based method

The main objective of the current work is to introduce a new conceptual linearization strategy to improve the performance of a primitive shock‐capturing pressure‐based finite‐volume method. To avoid a spurious oscillatory solution in the chosen collocated grids, both the primitive and extended methods utilize two convecting and convected momentum expressions at each cell face. The expressions are obtained via a physical‐based discretization of two inclusive statements, which are constructed via a novel incorporation of the continuity and momentum governing equations. These two expressions in turn provide a strong coupling among the Euler conservative statements. Contrary to the primitive work, the linearization in the current work respects the definitions and essence of physics behind deriving the Euler governing equations. The accuracy and efficiency of the new formulation are then investigated by solving the shock tube as a problem with moving normal and expansion waves and the converging‐diverging nozzle as a problem with strong stationary normal shock. The results show that there is good improvement in performance of the primitive pressure‐based shock‐capturing method while its superior accuracy is not deteriorated at all. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Proficiency testing in occupational and environmental health – Current practice and international initiatives for harmonisation

An increasing number of proficiency testing schemes (PT schemes) related to occupational and environmental health have been organised. Most schemes emphasise the importance of the validation of analytical results. With regards to harmonisation of the schemes, there are many differences between the schemes at many levels. These include factors concerning their relationship with legislation, national status, type and quality of proficiency testing material, analytical range and priorities for future development. Since differences between PT schemes have been recognised at the European level by organisers of PT schemes and external quality assessment schemes (EQASs) it seems appropriate to reinforce collaboration between scheme organisers in order to improve the quality of analytical performance in occupational and environmental health.     

Numerical experiments with several variant WENO schemes for the Euler equations

Numerical experiments with several variants of the original weighted essentially non‐oscillatory (WENO) schemes (J. Comput. Phys. 1996; 126 :202–228) including anti‐diffusive flux corrections, the mapped WENO scheme, and modified smoothness indicator are tested for the Euler equations. The TVD Runge–Kutta explicit time‐integrating scheme is adopted for unsteady flow computations and lower–upper symmetric‐Gauss–Seidel (LU‐SGS) implicit method is employed for the computation of steady‐state solutions. A numerical flux of the variant WENO scheme in flux limiter form is presented, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of low‐order schemes. Computations of unsteady oblique shock wave diffraction over a wedge and steady transonic flows over NACA 0012 and RAE 2822 airfoils are presented to test and compare the methods. Various aspects of the variant WENO methods including contact discontinuity sharpening and steady‐state convergence rate are examined. By using the WENO scheme with anti‐diffusive flux corrections, the present solutions indicate that good convergence rate can be achieved and high‐order accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Copyright © 2008 John Wiley & Sons, Ltd.     

Response time and vendor–assembler relationship in a supply chain

Relationships between an assembler and a vendor in a supply chain are investigated in two-period models when the assembler wants to reduce response time by incentive systems. The assembler may offer myopic or farsighted incentive contracts to the vendor, under short-term or long-term relationships. Incentive schemes, effort levels, and expected payoffs under different perspectives and relationships are examined. We find that a farsighted assembler provides the vendor with a higher incentive than a myopic assembler in the first period. A long (short)-term relationship is preferred if the value of farsightedness under a long-term relationship is greater (less) than the switching option value under a short-term relationship. We propose several sufficient conditions regarding which perspectives and relationships are preferred.     

Isotropic Subspaces of Polylinear Forms

The dimensions of maximal isotropic subspaces of symmetric and skew-symmetric forms of degree d in general position are found.     

A finite strain contact model for mixed lubricated surfaces in forming processes

For an accurate simulation of forming processes, it is of paramount importance to model the different lubrication regimes that can develop at the contact interface. These might vary from zone to zone of the forming piece, and from one regime to another, resulting in forces of different nature and magnitude. In these cases, the use of the classical Coulomb friction law will be clearly not sufficient to capture, in a suitable manner, the variety of forces applied on the forming piece.