Capillary Effect on Vertically Excited Surface Wave in Circular Cylindrical Vessel

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions was developed in inviscid fluids to investigate the motion of single free surface standing wave including the effect of surface tension.

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.

Abel Inversion of Axially-Symmetric Shock Wave Flows

Finite-fringe interferograms produced for axisymmetric shock wave flows are analyzed by Fourier transform fringe analysis and an Abel inversion method to produce density field data for the validation of numerical models. For the Abel inversion process, we use basis functions to model phase data from axially-symmetric shock wave structure. Steady and unsteady flow problems are studied, and compared with numerical simulations. Good agreement between theoretical and experimental results is obtained when one set of basis functions is used during the inversion process, but the shock front is smeared when another is used. This is because each function in the second set of basis functions is infinitely differentiable, making them poorly-suited to the modelling of a step function as is required in the representation of a shock wave.

Study on current-random wave forces acting on a framework

The forces of random wave plus current acting on a simplified offshore platform (jacket) model have been studied numerically and experimentally. The numerical results are in good agreement with experiments. The mean force can be approximated as a function of equivalent velocity parameter and the root-mean-square force as a function of equivalent significant wave height parameter.

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.

Non-Linear Wave-Induced Transient Response of Soil Around a Trenched Pipeline

Submarine pipelines are always trenched within a seabed for reducing wave loads and thereby enhancing their stability. Based on Biot’s poroelastic theory, a two-dimensional finite element model is developed to investigate non-linear wave-induced responses of soil around a trenched pipeline, which is verified with the flume test results by Sudhan et al. [Sudhan, C.M., Sundar, V., Rao, S.N., 2002. Wave induced forces around buried pipeline. Ocean Engineering, 29, 533–544] and Turcotte et al. [Turcotte, B.R., Liu, P.L.F., Kulhawy, F.H., 1984. Laboratory evaluation of wave tank parameters for wave-sediment interaction. Joseph H. Defree Hydraulic Laboratory Report 84-1, School of Civil and Environmental Engineering, Cornell University]. Non-linear wave-induced transient pore pressure around pipeline at various phases of wave loading is examined firstly. Unlike most previous investigations, in which only a single sediment layer and linear wave loading were concerned, in this study, the influences of the non-linearity of wave loading, the physical properties of backfill materials and the geometry profile of trenches on the excess pore pressures within the soil around pipeline, respectively, were explored, taking into account the in situ conditions of buried pipeline in the shallow ocean zones. Based on the parametric study, it is concluded that the shear modulus and permeability of backfill soils significantly affect the wave-induced excess pore pressures around trenched pipeline, and that the effect of wave non-linearity becomes more pronounced and comparable with that of trench depth, especially at high wave steepness in shallow water.

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.

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.

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.

Pressure and thermal effects on elastic properties of single crystal alpha-Al2O3 from Monte-Carlo simulation

A rectangular structural unit cell of a-Al2O3 is generated from its hexagonal one. For the rectangular structural crystal with a simple interatomic potential [Matsui, Mineral Mag. 58A, 571 (1994)], the relations of lattice constants to homogeneous pressure and temperature are calculated by using Monte-Carlo method at temperature 298K and 0 GPa, respectively. Both numerical results agree with experimental ones fairly well. By comparing pair distribution function, the crystal structure of a-Al2O3 has no phase transition in the range of systematic parameters. Based on the potential model, pressure dependence of isothermal bulk moduli is predicted. Under variation of general strains, which include of external and internal strains, elastic constants of a-Al2O3 in the different homogeneous load are determined. Along with increase of pressure, axial elastic constants increase appreciably, but nonaxial elastic constants are slowly changed.

Modeling Ammonothermal Growth of GaN Single Crystals: The Role of Transport

Single crystal gallium nitride (GaN) is an important technological material used primarily for the manufacture of blue light lasers. An important area of contemporary research is developing a viable growth technique. The ammonothermal technique is an important candidate among many others with promise of commercially viable growth rates and material quality. The GaN growth rates are a complicated function of dissolution kinetics, transport by thermal convection and crystallization kinetics. A complete modeling effort for the growth would involve modeling each of these phenomena and also the coupling between these. As a first step, the crystallization and dissolution kinetics were idealized and the growth rates as determined purely by transport were investigated. The growth rates thus obtained were termed ‘transport determined growth rates’ and in principle are the maximum growth rates that can be obtained for a given configuration of the system. Using this concept, a parametric study was conducted primarily on the geometric and the thermal boundary conditions of the system to optimize the ‘transport determined growth rate’ and determine conditions when transport might be a bottleneck.

How important are shock waves to single-bubble sonoluminescence?

By solving numerically the full set of hydrodynamic equations governing the pulsation of a bubble,we show that shock waves are often absent in a stable sonoluminescing bubble. Nevertheless, for a wide range of physical parameters, a continuous compressional wave emerges and heats up the bubble, and the resulting black-body radiations have pulse heights and widths that agree with experimental data. Shock waves, being much less robust, are not essential for stable single-bubble sonoluminescence.

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.

Effect of surface tension on the mode selection of vertically excited surface waves in a circular cylindrical vessel

Singular perturbation theory of two-time-scale expansions was developed in inviscid fluids to investigate patternforming, structure of the single surface standing wave, and its evolution with time in a circular cylindrical vessel subject to a vertical oscillation. A nonlinear slowly varying complex amplitude equation, which involves a cubic nonlinear term, an external excitation and the influence of surface tension, was derived from the potential flow equation. Surface tension was introduced by the boundary condition of the free surface in an ideal and incompressible fluid. The results show that when forced frequency is low, the effect of surface tension on the mode selection of surface waves is not important. However, when the forced frequency is high, the surface tension cannot be neglected. This manifests that the function of surface tension is to cause the free surface to return to its equilibrium configuration. In addition, the effect of surface tension seems to make the theoretical results much closer to experimental results.

Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media

In this paper, we examine the characteristics of elastic wave propagation in viscoelastic porous media, which contain simultaneously both the Biot-flow and the squirt-flow mechanisms (BISQ). The frequency-domain Green's functions for viscoelastic BISQ media are then derived based on the classic potential function methods. Our numerical results show that S-waves are only affected by viscoelasticity, but not by squirt-flows. However, the phase velocity and attenuation of fast P-waves are seriously influenced by both viscoelasticity and squirt-flows; and there exist two peaks in the attenuation-frequency variations of fast P-waves. In the low-frequency range, the squirt-flow characteristic length, not viscoelasticity, affects the phase velocity of slow P-waves, whereas it is opposite in the high-frequency range. As to the contribution of potential functions of two types of compressional waves to the Green's function, the squirt-flow length has a small effect, and the effects of viscoelastic parameter are mainly in the higher frequency range. Crown Copyright (C) 2006 Published by Elsevier Ltd. All rights reserved.