Hydrodynamic interaction between two vertical cylinders in water waves

The hydrodynamic interaction between two vertical cylinders in water waves is investigated based on the linearized potential flow theory. One of the two cylinders is fixed at the bottom while the other is articulated at the bottom and oscillates with small amplitudes in the direction of the incident wave. Both the diffracted wave and the radiation wave are studied in the present paper. A simple analytical expression for the velocity potential on the surface of each cylinder is obtained by means of Graf's addition theorem. The wave-excited forces and moments on the cylinders, the added masses and the radiation damping coefficients of the oscillating cylinder are all expressed explicitly in series form. The coefficients of the series are determined by solving algebraic equations. Several numerical examples are given to illustrate the effects of various parameters, such as the separation distance, the relative size of the cylinders, and the incident angle, on the first-order and steady second-order forces, the added masses and radiation-damping coefficients as well as the response of the oscillating cylinder.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

On the velocity of buckle propagation in a beam on a nonlinear elastic-foundation

The paper revisits a simple beam model used by Chater et al. (1983, Proc. IUTAM Symp. Collapse, Cambridge University Press) to examine the dynamics of propagating buckles on it. It was found that, if a buckle is initiated at a constant pressure higher than the propagation pressure of the model (P-p), the buckle accelerates and gradually reaches a constant velocity which depends upon the pressure, while if it is initiated at P-p, the buckle propagates at a velocity which depends upon the initial imperfection. The causes for the difference are also investigated.

The tangential discontinuous surface of velocity in solar chromosphere

The data of velocity and magnetic fields in the solar photosphere (5324 angstrom) and the chromosphere (4861 angstrom) clearly show the features of tangential discontinuity of velocity in the chromosphere. The velocity fields in and near the solar active region named No. 88029 by the Huairou Station have been analyzed in detail. A lot of magnetohydrodynamic discontinuous surfaces, especially the tangential discontinuities, are shown from the observations. The calculations of the thickness of discontinuous layer and the evolution time of instability agree with the observational results. The variations of the flow field will directly influence the evolutions and changes of the active region as the magnetic field are coupled closely with the plasma motion.

Basic equations for suspension flows with interphase mass-transfer

In the case of suspension flows, the rate of interphase momentum transfer M(k) and that of interphase energy transfer E(k), which were expressed as a sum of infinite discontinuities by Ishii, have been reduced to the sum of several terms which have concise physical significance. M(k) is composed of the following terms: (i) the momentum carried by the interphase mass transfer; (ii) the interphase drag force due to the relative motion between phases; (iii) the interphase force produced by the concentration gradient of the dispersed phase in a pressure field. And E(k) is composed of the following four terms, that is, the energy carried by the interphase mass transfer, the work produced by the interphase forces of the second and third parts above, and the heat transfer between phases. It is concluded from the results that (i) the term, (-alpha-k-nabla-p), which is related to the pressure gradient in the momentum equation, can be derived from the basic conservation laws without introducing the "shared-pressure presumption"; (ii) the mean velocity of the action point of the interphase drag is the mean velocity of the interface displacement, upsilonBAR-i. It is approximately equal to the mean velocity of the dispersed phase, upsilonBAR-d. Hence the work terms produced by the drag forces are f(dc) . upsilonBAR-d, and f(cd) . upsilonBAR-d, respectively, with upsilonBAR-i not being replaced by the mean velocity of the continuous phase, upsilonBAR-c; (iii) by analogy, the terms of the momentum transfer due to phase change are upsilonBAR-d-GAMMA-c, and upsilonBAR-d-GAMMA-d, respectively; (iv) since the transformation between explicit heat and latent heat occurs in the process of phase change, the algebraic sum of the heat transfer between phases is not equal to zero. Q(ic) and Q(id) are composed of the explicit heat and latent heat, so that the sum Q(ic) + Q(id)) is equal to zero.

Ballistic impact of ceramics

The impact behaviour of a range of glass and ceramic materials has been studied using high-speed photography. A gas gun was used to project hardened spheres at plate specimens in the velocity range 30 to 1000m s-1. The target materials included soda-lime glass, boron carbide and various glass ceramics and aluminas. The performance of a particular ceramic was found to depend on a combination of parameters but of key importance was the relative hardness of the projectile and target materials. The fracture toughness, K(IC), had only a secondary effect.

Simultaneous measurement of temperature, density and velocity in gas-flows by modulated photoluminescence

A novel possibility to determine the temperature, density and velocity simultaneously in gas flows by measuring the average value, amplitude of modulation and phase shift of the photoluminescence excited by a temporally or spatially modulated light source is investigated. Time-dependent equations taking the flow, diffusion, excitation and decay into account are solved analytically. Different experimental arrangements are proposed. Measurements of velocity with two components, and temporal and spatial resolutions in the measurements are investigated. Numerical examples are given for N z with biacetyl as the seed gas. Practical considerations for the measurements and the relation between this method and some existing methods of lifetime measurement are discussed.

Coflow jets with large velocity difference - a new aerodynamic principle of flame stabilization and combustion intensification

A new aerodynamic principle of flame stabilization and combustion intensification, the coflow jets with large velocity difference, is described. One or more small high-velocity jets of air or steam, injected off the axis and in the same direction as the low-velocity main fuel-air flow into the combustor, create a large recirculation zone of high turbulence intensity in which the combustibles and high temperature gases are effectively mixed, so that stable and intensive combustion can be maintained even for fuels with poor ignition. A pulverized coal combustor based on the principle mentioned above is shown to be characteristic of excellent combustoom and a simple structure. A number of precombustors of this type are in operation at some power stations and industrial boilers of China. Using such precombustor, successtul startups and part-load operation of the boilers have become available under conditions of unpreheated air and low-grade coal with volatiles as low as 15% and ash content as high as 30%. This principle shows good promise as an attractive new technology of combustion.

Distributions of temperature, velocity and viscosity at symmetrical axis of cylindrical plume in earth's mantle with variable viscosity

The parameters at the symmetrical axis of a cylindrical plume characterize the strength of this plume and provide a boundary condition which must be given to investigate the structure of a plume. For Newtonian fluid with a temperature-and pressure-dependence viscosity, an asymptotical solution of hydrodynamic equations at the symmetrical axis of the plume is found in the present paper. The temperature, upward velocity and viscosity at the symmetrical axis have been obtained as functions of depth, The calculated results have been given for two typical sets of Newtonian rheological parameters. The results obtained show that the temperature distribution along the symmetrical axis is nearly independent of the theological parameters. The upward velocity at the symmetrical axis, however, is strongly dependent on the rheological parameters.