A Generalized Soft-Sphere for Monte Carlo Simulation

A new collision model, called the generalized soft-sphere (GSS) model, is introduced. It has the same total cross section as the generalized hard-sphere model [Phys. Fluids A 5, 738 (1993)], whereas the deflection angle is calculated by the soft-sphere scattering model [Phys. Fluids A 3, 2459 (1991)]. In virtue of a two-term formula given to fit the numerical solutions of the collision integrals for the Lennard-Jones (6-12) potential and for the Stockmayer potential, the parameters involved in the GSS model are determined explicitly that may fully reproduce the transport coefficients under these potentials. Coefficients of viscosity, self-diffusion and diffusion for both polar and nonpolar molecules given by the GSS model and experiment are in excellent agreement over a wide range of temperature from low to high.

Dynamic behavior of a cylindrical crack in a functionally graded interlayer under torsional loading

In this paper the problem of a cylindrical crack located in a functionally graded material (FGM) interlayer between two coaxial elastic dissimilar homogeneous cylinders and subjected to a torsional impact loading is considered. The shear modulus and the mass density of the FGM interlayer are assumed to vary continuously between those of the two coaxial cylinders. This mixed boundary value problem is first reduced to a singular integral equation with a Cauchy type kernel in the Laplace domain by applying Laplace and Fourier integral transforms. The singular integral equation is then solved numerically and the dynamic stress intensity factor (DSIF) is also obtained by a numerical Laplace inversion technique. The DSIF is found to rise rapidly to a peak and then reduce and tend to the static value almost without oscillation. The influences of the crack location, the FGM interlayer thickness and the relative magnitudes of the adjoining material properties are examined. It is found among others that, by increasing the FGM gradient, the DSIF can be greatly reduced.

Mixing process in estuaries

A theoretical expression for vertical profile-of horizontal velocity in terms of its depth-average is derived based on oscillatory boundary layer theory and estuarine flow characteristics. The derived theoretical profile is then incorporated into a vertical quasi-two-dimensional model, which is proved advantageous in more physical implications and less CPU time demand. To validate the proposed model, the calculated results are compared to the field data in the Yangtze River Estuary, exhibiting good agreement with observations. The proposed quasi-two-dimensional vertical model is used to study mixing process, especially dependence of salinity distribution and salt front strength on runoff and tides in estuaries.

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, 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 damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, 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. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

Scour of the seabed under a pipeline in oscillating flow

The scour of the seabed under a pipeline is studied experimentally in this paper. Tests are carried out in a U-shaped oscillatory water tunnel with a box imbedded in the bottom of the test section. By use of the standard sand, clay and plastic grain as the seabed material, the influence of the bed material on the scour is studied. The relationship between the critical initial gap-to-diameter ratio above which no scour occurs and the parameters of the oscillating flow is obtained. The self-burial phenomenon. which occurs for the pipeline not fixed to two sidewalls of the test section, is not observed for the Bred pipeline. The effect of the pipe on sand wave formation is discussed. The maximum equilibrium scour depths For different initial gap-to-diameter ratios, different Kc numbers and different bed sands are also given in this paper.

Wave-soil-pipe coupling effect on submarine pipeline on-bottom stability

Wave-soil-pipe coupling effect on the untrenched pipeline stability on sands is for the first time investigated experimentally. Tests are conducted in the U-shaped water tunnel, which generates an oscillatory how, simulating the water particle movements with periodically changing direction under the wave action. Characteristic times and phases during the instability process are revealed. Linear relationship between Froude number and non-dimensional pipe weight is obtained. Effects of initial embedment and loading history are observed. Test results between the wavesoil-pipe interaction and pipe-soil interaction under cyclic mechanical loading are compared. The mechanism is briefly discussed. For applying in the practical design, more extensive and systematic investigations are needed.

Experimental Study on Thermocapillary Convection

In the present paper, the experimental studies on thermocapillary convection are reviewed. The author's interest is mainly focused on the onset of oscillatory thermocapillary convection, the features of oscillatory flow pattern, and the critical Marangoni number related with temperature and free surface oscillation. The coordinated measurement in a microgravity environment of a drop shaft is also addressed.

Simulating shock to detonation transition: Algorithm and results

An algorithm based on flux-corrected transport and the Lagrangian finite element method is presented for solving the problem of shock dynamics. It is verified through the model problem of one-dimensional strain elastoplastic shock wave propagation that the algorithm leads to stable, non-oscillatory results. Shock initiation and detonation wave propagation is simulated using the algorithm, and some interesting results are obtained. (C) 1999 Academic Press.

General Solution to Two-Dimensional Nonslipping JKR Model with a Pulling Force in an Arbitrary Direction

In this paper, a generalized JKR model is investigated, in which an elastic cylinder adhesively contacts with an elastic half space and the contact region is assumed to be perfect bonding. An external pulling force is acted on the cylinder in an arbitrary direction. The contact area changes during the pull-off process, which can be predicted using the dynamic Griffith energy balance criterion as the contact edge shifts. Full coupled solution with an oscillatory singularity is obtained and analyzed by numerical calculations. The effect of Dundurs' parameter on the pull-off process is analyzed, which shows that a nonoscillatory solution can approximate the general one under some conditions, i.e., larger pulling angle (pi/2 is the maximum value), smaller a/R or larger nondimensional parameter value of Delta gamma/E*R. Relations among the contact half width, the external pulling force and the pulling angle are used to determine the pull-off force and pull-off contact half width explicitly. All the results in the present paper as basic solutions are helpful and applicable for experimenters and engineers.

Dynamic SIF of interface crack between two dissimilar viscoelastic bodies under impact loading

In this paper, the transient dynamic stress intensity factor (SIF) is determined for an interface crack between two dissimilar half-infinite isotropic viscoelastic bodies under impact loading. An anti-plane step loading is assumed to act suddenly on the surface of interface crack of finite length. The stress field incurred near the crack tip is analyzed. The integral transformation method and singular integral equation approach are used to get the solution. By virtue of the integral transformation method, the viscoelastic mixed boundary problem is reduced to a set of dual integral equations of crack open displacement function in the transformation domain. The dual integral equations can be further transformed into the first kind of Cauchy-type singular integral equation (SIE) by introduction of crack dislocation density function. A piecewise continuous function approach is adopted to get the numerical solution of SIE. Finally, numerical inverse integral transformation is performed and the dynamic SIF in transformation domain is recovered to that in time domain. The dynamic SIF during a small time-interval is evaluated, and the effects of the viscoelastic material parameters on dynamic SIF are analyzed.

A simulation model of a floating half zone

A simulation model of a floating half zone was suggested by steady numerical simulation and experiment respectively, in the previous papers [Q.S. Chen, W.R. Hu, Int. J. Mass Heat Transfer 40 (1997) 757; J.H. Han, Y. Ar, R. Zhou, W.R. Hue, Int. J. Mass Heat Transfer 40 (1997) 2671]. In the present paper, the simulation model is studied by using the method of unsteady and three-dimensional numerical simulation, and the transient process from steady convection to oscillatory convection is especially analyzed. Comparison of onsets of oscillation for both simulation model and the usual model were obtained, and the results show that the critical Marangoni number of the simulation model is obviously smaller than that of the usual model for the same slender liquid bridge. This implies that the usual model of a floating half zone gives a lower estimation on the onset of oscillation for floating zone convection.

Damping of Vertically Excited Surface Wave in Weakly Viscous Fluid

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress,intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored.

Two bifurcation transitions of the floating half zone convection in a fat liquid bridge of larger Pr

The transient process of the thermocapillary convection was obtained for the large Pu floating half zone by using the method of three-dimensional and unsteady numerical simulation. The convection transits directly from steady and axisymmetric state to oscillatory flow for slender liquid bridge, and transits first from steady and axisymmetric convection to the steady and non-axisymmetric convection, then, secondly to the oscillatory convection for the fatter liquid bridge. This result implies that the volume of liquid bridge is not only a sensitive critical parameter for the onset of oscillation, but also relates to the new mechanism for the onset of instability in the floating half zone convection even in case of large Prandtl number fluid.

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.