The application of B-P constitutive equations in finite element analysis of high velocity impact

We present in this paper the application of B-P constitutive equations in finite element analysis of high velocity impact. The impact process carries out in so quick time that the heat-conducting can be neglected and meanwhile, the functions of temperature in equations need to be replaced by functions of plastic work. The material constants in the revised equations can be determined by comparison of the one-dimensional calculations with the experiments of Hopkinson bar. It can be seen from the comparison of the calculation with the experiment of a tungsten alloy projectile impacting a three-layer plate that the B-P constitutive equations in that the functions of temperature were replaced by the functions of plastic work can be used to analysis of high velocity impact.

Diffusive Wave Solutions for Open Channel Flows with Uniform and Concentrated Lateral Inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow intoa river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form interms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomesS-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usuallyresulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral infloware analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively,hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infiniteand finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channeltypes.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

Propagation Properties of Electromagnetic Wave in a Plasma Slab

In this paper, the calculated results about the propagation properties of electromagnetic wave in a plasma slab are described. The relationship of the propagation properties with frequencies of electromagnetic wave, and parameters of plasma (electron temperature, electron density, dimensionless collision frequency and the size of the plasma slab) is analyzed.

Hydrothermal Wave in a Shallow Liquid Layer

The oscillatory thermocapillary convection and hydrothermal wave in a shallow liquid layer, where a temperature difference is applied between two parallel sidewalls, have been numerically investigated in a two-dimensional model. The oscillatory thermocapillary convection and hydrothermal wave appear if the Marangoni number is larger than a critical value. The critical phase speed and critical wave number of the hydrothermal wave agree with the ones given analytically by Smith and Davis in the microgravity environment, and it travels in the direction opposed to the surface flow. Another wave traveled downstream in addition to the hydrothermal wave traveled upstream was observed in the case of earth gravity condition.

A New Version of Smoothed Point Method to Simulate Linear Elastic Wave Propagation

By using the kernel function of the smoothed particle hydrodynamics (SPH) and modification of statistical volumes of the boundary points and their kernel functions, a new version of smoothed point method is established for simulating elastic waves in solid. With the simplicity of SPH kept, the method is easy to handle stress boundary conditions, especially for the transmitting boundary condition. A result improving by de-convolution is also proposed to achieve high accuracy under a relatively large smooth length. A numerical example is given and compared favorably with the analytical solution.

The Optical Methods on Measuring Surface Deformation and Surface Wave in The Thermal Capillary Convection

An optical diagnostic system consisting of the Michelson interferometer with the image processor has been developed for the study of the kinetics of the thermal capillary convection. The capillary convection, surface deformation, surface wave and the velocity field in a rectangular cavity with different temperature's sidewalls have been investigated by optical interference method and PIV technique. In order to calculate the surface deformation from the interference fringe, Fourier transformation is used to grating analysis. The quantitative results of the surface deformation and surface wave have been calculated from the interference fringe pattern.

A further discussion on basic equations for two-phase flow: On collision stress of solid phase

The following points are argued: (i) there are two independent kinds of interaction on interfaces, i.e. the interaction between phases and the collision interaction, and the jump relations on interfaces can accordingly be resolved; (ii) the stress in a particle can also be divided into background stress and collision stress corresponding to the two kinds of interaction on interfaces respectively; (iii) the collision stress, in fact, has no jump on interface, so the averaged value of its derivative is equal to the derivative of its averaged value; (iv) the stress of solid phase in the basic equations for two\|phase flow should include the collision stress, while the stress in the expression of the inter\|phase force contains the background one only. Based on the arguments, the strict method for deriving the equations for two\|phase flow developed by Drew, Ishii et al. is generalized to the dense two\|phase flow, which involves the effect of collision stress.

Roll Waves in Overland Flow

Roll waves are frequently observed in overland flow, especially in rill flow, which has an important effect on the development of soil erosion. Using one-dimensional St. Venant equations, this paper investigates the dynamics of periodic roll waves based on Dressler’s and Brock’s work. Under the assumption that the average flow depth equals the uniform flow depth, expressions of the roll-wave speed and roll-wave profile were obtained. Testing with the results observed by Brock (1970) for wave properties shows that these expressions can approximately describe the characteristics of periodic permanent roll waves. Numerical solutions of roll waves under specific conditions are found, which show that when a roll wave appears, the shear stress of flow increases, and the soil erosion accelerates.

Shock wave diffraction and reflection around a dusty square cavity

The diffraction and reflection of planar shock wave around a dusty square cavity is investigated numerically, which is embedded in the net bottom surface of a two-dimensional channel, and the induced gas-particle two-phase now. The wave patterns at different times are obtained for three different values of the particle diameter. The computational results show that the existence of particles affects appreciably the shock wave diffraction and cavity flow.

Variation of drag coefficient with wind wave development

Turbulent air flows over developing wind waves in the air-sea boundary layer are numerically simulated without considering wave breaking. Influences of wind waves on air flows are considered using a model of significant wave and surface roughness, with a formula proposed for calculating the surface roughness, k - epsilon model is adopted to simulate turbulent flows. The results of the drag coefficient and turbulence characteristics agree well with the observations.

Numerical Investigation of Toroidal Shock Wave Focusing in a Cylindrical Chamber

In this paper, focusing of a toroidal shock wave propagating from an annular shock tube into a cylindrical chamber was investigated numerically with the dispersion controlled dissipation (DCD) scheme. The first case for an incident Mach number of 1.5 was conducted and compared with experiments for validation. Then, several cases were calculated for higher incident Mach numbers varying from 2.0 to 5.0, and complicated flow structures were observed. The numerical study was mainly focused on two aspects: focusing process and flow structures. The process, including diffraction, focusing, and reflection, is displayed to reveal the focusing mechanism, and the flow structures at different incident. Mach numbers are used to demonstrate shock reflection styles and focusing characteristics.

Improved Analysis Method for Wave-Induced Pipeline Stability on Sandy Seabed

The existing Det Norske Veritas DNV Recommended Practice RP E305 for pipeline on-bottom stability is mainly based on the pipe–soil interaction model reported by Wagner et al. in 1987, and the wake model reported by Lambrakos et al. in 1987, to calculate the soil resistance and the hydrodynamic forces upon pipeline, respectively. Unlike the methods in the DNV Practice, in this paper, an improved analysis method is proposed for the on-bottom stability of a submarine pipeline, which is based on the relationships between Um/ gD 0.5 and Ws / D2 for various restraint conditions obtained by the hydrodynamic loading experiments, taking into account the coupling effects between wave, pipeline, and sandy seabed. The analysis procedure is illustrated with a detailed flow chart. A comparison is made between the submerged weights of pipeline predicted with the DNV Practice and those with the new method. The proposed analysis method may provide a helpful tool for the engineering practice of pipeline on-bottom stability design.

NS方程激波计算的摄动有限差分方法

摄动有限差分（PFD）方法从一阶迎风差分格式出发，将差分系数展开为网格步长的幂级数，通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式，因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导，将摄动有限差分格式引入时间相关法的计算，并构造了守恒形式的摄动有限差分格式，然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明：本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率，而且保持了一阶迎风格式的无振荡性质。

Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0.similar to 4.7) x 10(7)/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30 degrees for a sweepback angle of 67.6 degrees. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10 degrees and secondary separation were detected at deflection of theta greater than or equal to 20 degrees. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements at M greater than or equal to 6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers.