Numerical study on onset of oscillatory thermocapillary flow in rectangular liquid pool

Thermocapillary flow in a rectangular liquid pool of large Prandtl fluid (Pr = 105.6) is numerically studied in microgravity. Oscillatory thermocapillary flow arises when the imposed temperature difference between the sidewalls exceeds a critical value. The fluctuations of the oscillatory flow, accompanied by the propagation of the hydrothermal wave from the cold sidewall to the hot one, are much smaller than the time-averaged velocity and temperature fields. The corresponding disturbance cells arise in the centre of the liquid pool initially, and extend to the whole region with the increasing imposed temperature difference. The present study reveals the different characteristics of the oscillatory themocapillary flow in the rectangular liquid pool as compared to the cases in other configurations.

Fracture mechanics of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in an infinite piezoelectric or on the interface of piezoelectric bimaterials. For homogeneous materials it is found that the normal electric displacement D-2, induced by the crack, is constant along the crack faces which depends only on the remote applied stress fields. Within the crack slit, the perturbed electric fields induced by the crack are also constant and not affected by the applied electric displacement fields. For bimaterials, generally speaking, an interface crack exhibits oscillatory behavior and the normal electric displacement D-2 is a complex function along the crack faces. However, for bimaterials, having certain symmetry, in which an interface crack displays no oscillatory behavior, it is observed that the normal electric displacement D-2 is also constant along the crack faces and the electric field E-2 has the singularity ahead of the crack tip and has a jump across the interface. Energy release rates are established for homogeneous materials and bimaterials having certain symmetry. Both the crack front parallel to the poling axis and perpendicular to the poling axis are discussed. It is revealed that the energy release rates are always positive for stable materials and the applied electric displacements have no contribution to the energy release rates.

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.

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.

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.

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.

一个基于非均匀网格上的二阶基本无振荡差分格式

在基本无振荡格式的构造中，将通常的对流通量f的逼近方式推广到对通量导数的逼近，这一构造方法可以有效地应用到非均匀或非结构网格。直接基于非均匀网格上，构造了一个二阶的基本无振荡(ENO)差分格式。该格式具有形式简单，对网格的划分灵活，与传统格式相比不增加计算量等优点，儿个数值算例证明了格式的有效性。

Two Bifurcation Transition Processes in Floating Half Zone Convection of Larger Prandtl Number Fluid

Processes of the onset oscillation in the thermocapillaxy convection under the Earth's gravity are investigated by the numerical simulation and experiments in a floating half zone of large Prandtl number with different volume ratio. Both computational and experimental results show that the steady and axisymmetric convection turns to the oscillatory convection of m=1 for the slender liquid bridge, and to the oscillatory convection before a steady and 3D asymmetric state for the case of a fat liquid bridge. It implies that, there are two critical Marangoni numbers related, respectively, to these two bifurcation transitions for the fat liquid bridge. The computational results agree with the results of ground-based experiments.