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

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分（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.

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

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

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.

Spall damage in aluminum-alloy

A void growth relations for ductile porous materials under intense dynamic general loading condition is presented. The mathematical model includes the influence of inertial effects, material rate sensitivity, as well as the contribution of void surface energy and material work-hardening. Numerical analysis shows that inertia appears to resist the growth of voids. The inertial effects increase quickly with the loading rates. The theoretical analysis suggests that the inertial effects cannot be neglected at high loading rates. Plate-impact tests of aluminum alloy are performed with light gas gun. The processes of dynamic damage in aluminum alloy are successfully simulated with a finite-difference dynamic code in which the theoretical model presented in this paper is incorporated.

Dynamic failure in ductile porous materials

In this paper, a mathematical model of dynamic fracture in porous ductile materials under intense dynamic general loading is developed. The mathematical model includes the influence of inertial effects and material rate sensitivity, as well as the contribution of surface energy of a void and material work-hardening. In addition, the condition of the void compaction is considered as well. The threshold stresses for the void growth and compaction are obtained. A simple criterion for ductile fracture which is associated with material distention and plastic deformation is adopted. As an application of the theoretical model, the processes of two-dimensional spallation in LY12 aluminum alloy are successfully simulated by means of two-dimensional finite-difference Lagrangian code.

Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.