双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分（PFD）格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过娄似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分（FPFD）格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。

Evolution of three-dimensional coherent structures in compressible axisymmetric jet

A high order difference scheme is used to simulate the spatially developing compressible axisymmetric jet. The results show that the Kelvin-Helmholtz instability appears first when the jet loses its stability, and then with development of jet the increase in nonlinear effects leads to the secondary instability and the formation of the streamwise vortices. The evolution of the three-dimensional coherent structure is presented. The computed results verify that in axisymmetric jet the secondary instability and formation of the streamwise vortices are the important physical mechanism of enhancing the flow mixing and transition occurring.

Liquefaction and Displacement of Saturated Sand Under Vertical Vibration Loading

In order to investigate the influence of the vertical vibration loading on the liquefaction of saturated sand, one dimensional model for the saturated sand with a vertical vibration is presented based on the two phase continuous media theory. The development of the liquefaction and the liquefaction region are analyzed. It is shown that the vertical vibration loading could induce liquefaction. The rate of the liquefaction increases with the increase of the initial limit strain or initial porosity or amplitude and frequency of loading, and increases with the decrease of the permeability or initial modulus. It is shown also that there is a phase lag in the sand column. When the sand permeability distribution is non-uniform, the pore pressure and the strain will rise sharply where the permeability is the smallest, and fracture might be induced. With the development of liquefaction, the strength of the soil foundation becomes smaller and smaller. In the limiting case, landslides or debris flows could occur.

Can frequency diversity provide performance gains for WSNs at 2.4GHz for the fire hydrant to above ground channel

Wireless Sensor Networks (WSNs) which utilise IEEE 802.15.4 technology offer the potential for low cost deployment and maintenance compared with conventional wired sensor networks, enabling effective and efficient condition monitoring of aged civil engineering infrastructure. We will address wireless propagation for a below to above ground scenario where one of the wireless nodes is located in a below ground fire hydrant chamber to permit monitoring of the local water distribution network. Frequency Diversity (FD) is one method that can be used to combat the damaging effects of multipath fading and so improve the reliability of radio links. However, no quantitative investigation concerning the potential performance gains from the use of FD at 2.4GHz is available for the outlined scenario. In this paper, we try to answer this question by performing accurate propagation measurements using modified and calibrated off-the-shelf 802.15.4 based sensor nodes. These measurement results are also compared with those obtained from simulations that employ our Modified 2D Finite-Difference Time-Domain (FDTD) approach. ©2009 IEEE.

Frequency diversity measurements at 2.4 GHz for wireless sensor networks deployed in tunnels

Wireless Sensor Networks (WSNs) which utilise IEEE 802.15.4 technology operate primarily in the 2.4 GHz globally compatible ISM band. However, the wireless propagation channel in this crowded band is notoriously variable and unpredictable, and it has a significant impact on the coverage range and quality of the radio links between the wireless nodes. Therefore, the use of Frequency Diversity (FD) has potential to ameliorate this situation. In this paper, the possible benefits of using FD in a tunnel environment have been quantified by performing accurate propagation measurements using modified and calibrated off-the-shelf 802.15.4 based sensor motes in the disused Aldwych underground railway tunnel. The objective of this investigation is to characterise the performance of FD in this confined environment. Cross correlation coefficients are calculated from samples of the received power on a number of frequency channels gathered during the field measurements. The low measured values of the cross correlation coefficients indicate that applying FD at 2.4 GHz will improve link performance in a WSN deployed in a tunnel. This finding closely matches results obtained by running a computational simulation of the tunnel radio propagation using a 2D Finite-Difference Time-Domain (FDTD) method. ©2009 IEEE.

非均布表面应力作用下薄板变形问题的求解

给出非均布表面应力作用下弹性薄板挠曲变形问题的控制方程及边界条件,通过与热应力问题进行物理比拟,对这一问题进行了求解,并采用这一方法对实验中观测到的局部弯曲现象进行了理论解释.

Effect of imperfections on the yielding of two-dimensional foams

The influence of each of the six different types of morphological imperfection - waviness, non-uniform cell wall thickness, cell-size variations, fractured cell walls, cell-wall misalignments, and missing cells - on the yielding of 2D cellular solids has been studied systematically for biaxial loading. Emphasis is placed on quantifying the knock-down effect of these defects on the hydrostatic yield strength and upon understanding the associated deformation mechanisms. The simulations in the present study indicate that the high hydrostatic strength, characteristic of ideal honeycombs, is reduced to a level comparable with the deviatoric strength by several types of defect. The common source of this large knock-down is a switch in deformation mode from cell wall stretching to cell wall bending under hydrostatic loading. Fractured cell edges produce the largest knock-down effect on the yield strength of 2D foams, followed in order by missing cells, wavy cell edges, cell edge misalignments, Γ Voronoi cells, δ Voronoi cells, and non-uniform wall thickness. A simple elliptical yield function with two adjustable material parameters successfully fits the numerically predicted yield surfaces for the imperfect 2D foams, and shows potential as a phenomenological constitutive law to guide the design of structural components made from metallic foams.

Low-Dimensional Dynamical Systems for a Wake-Type Shear Flow with Vortex Dislocations

Low-dimensional systems are constructed to investigate dynamics of vortex dislocations in a wake-type shear flow. High-resolution direct numerical simulations are employed to obtain flow snapshots from which the most energetic modes are extracted using proper orthogonal decomposition (POD). The first 10 modes are classified into two groups. One represents the general characteristics of two-dimensional wake-type shear flow, and the other is related to the three-dimensional properties or non-uniform characteristics along the span. Vortex dislocations are generated by these two kinds of coherent structures. The results from the first 20 three-dimensional POD modes show that the low- dimensional systems have captured the basic properties of the wake-type shear flow with vortex dislocation, such as two incommensurable frequencies and their beat frequency.

非均质泥流的输移特性

非均质泥流是一种含有大量悬浮细颗粒的高浓度固液两相流体，也是泥石流的一种特例。作者依据水槽输沙试验资料建立的泥流不淤流速公式，通过实测资料检验了其可靠性。在此基础上考虑了推移部分影响，推导得出非均质泥流的不淤流速、临界坡降、沟道断面形态指标等重要的输移临界参数计算公式。最后列举计算实例，提出应用计算方法，并据此进一步揭示了非均质泥流的输移特性。

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Vertically forced surface wave in weakly viscous fluids bounded in a circular cylindrical vessel

Two-time scale perturbation expansions were developed in weakly viscous fluids to investigate surface wave motions by linearizing the Navier-Stokes equation in a circular cylindrical vessel which is subject to a vertical oscillation. The fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates a damping term and external excitation, was derived for the weakly viscid fluids. The condition for the appearance of stable surface waves was obtained and the critical curve was determined. In addition, an analytical expression for the damping coefficient was determined and the relationship between damping and other related parameters (such as viscosity, forced amplitude, forced frequency and the depth of fluid, etc.) was presented. Finally, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent.

A coupled model for multiphase fluid flow and sedimentation deformation in oil reservoir and its numerical simulation

A mathematical model for coupled multiphase fluid flow and sedimentation deformation is developed based on fluid-solid interaction mechanism. A finite difference-finite element numerical approach is presented. The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances, and the coupled model has practical significance for oilfield development.

Periodic vortex breakdown in wide spherical gaps

The flow field with vortex breakdown in wide spherical gaps was studied numerically by a finite difference method under the axisymmetric condition. The result shows that the flow bifurcates to periodic motion as the Reynolds number or the eccentricity of the spheres increases. (C) 1997 American Institute of Physics.

Numerical investigation on a simulation model of floating zone convection

A simulation model with adiabatic condition at the upper rod and constant temperature at the lower rod is studied numerically in this paper. The temperature distribution in a simulation model is closer to the one in the half part of a floating full zone in comparison with the one in a usual floating half zone model with constant temperature at both rods, because the temperature distribution of a floating full zone is symmetric for the middle plane in a microgravity environment. The results of the simulation model show that the temperature profiles and the how patterns are different from those of the usual floating half zone model. Another type of half zone model, with a special non-uniform temperature distribution at the upper rod and constant temperature at the lower rod, has been suggested by recent experiments. The temperature boundary condition of the upper rod has a maximum value in the center and a lower value near the free surface. This modified simulation model is also simulated numerically in the present paper. Copyright (C)1996 Elsevier Science Ltd.