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.

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.

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.

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.

Vertically forced surface wave in weakly viscous fluids bounded in a circular cylindrical vessel

Two-time scale perturbation expansions were developed in weakly viscous fluids to investigate surface wave motions by linearizing the Navier-Stokes equation in a circular cylindrical vessel which is subject to a vertical oscillation. 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 a damping term and external excitation, was derived for the weakly viscid fluids. The condition for the appearance of stable surface waves was obtained and the critical curve was determined. In addition, an analytical expression for the damping coefficient was determined and the relationship between damping and other related parameters (such as viscosity, forced amplitude, forced frequency and the depth of fluid, etc.) was presented. Finally, 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.

Wave generation by the falling rock in the two-dimensional wave tank

Wave generation by the falling rock in the two-dimensional wave tank is experimentally and numerically studied, where the numerical model utilizes the boundary element method to solve the fully nonlinear potential flow theory. The wave profiles at different times are measured in the laboratory, which are also used to test the numerical model. Comparisons show that the experimental and numerical results are in good agreement, and the numerical model can be used to simulate the wave generation due to the submarine rock falling. Further numerical tests on the influences of the rock size, density, initial position and the falling angle on the wave elevation of the generated waves are performed, respectively. The results show that the size and density of the rock have strong effects on the maximum elevation of the generated wave, while the effects of the initial position and the falling angle of the rock are also significant. When the size or the density of the rock increases, the maximum elevation of the generated wave increases. The same effect on the generated wave would be produced if the initial position of the rock becomes closer to the surface, or the falling angle between the falling route and the vertical direction turns larger. In addition, the present numerical tests reveal that the submarine rock falling provides a new generation method for the breaking wave in the wave tank.

Application of Electrical Resistance Tomography System to Monitor Gas/Liquid Two-Phase Flow in a Horizontal Pipe

The Electrical Resistance Tomography (ERT) technique possesses great potential in monitoring widely exiting industrial two/multi-phase flow. For vertical pipe flow and inclined pipe flow, some application studies with exciting results have been reported, but there is rarely a paper regarding the application of ERT to horizontal gas/liquid pipe flow. This paper addresses this issue and proposes a smart method, Liquid Level Detection method, to conventional ERT system. The enhanced ERT system using the new method can monitor horizontal pipe flow effectively and its application is no longer restricted by the flow conditions. Some experimental results from monitoring an air/water slug pipe flow are presented.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

Atomistic Origin Of Rate-Dependent Serrated Plastic Flow In Metallic Glasses

Nanoindentation simulations on a binary metallic glass were performed under various strain rates by using molecular dynamics. The rate-dependent serrated plastic flow was clearly observed, and the spatiotemporal behavior of its underlying irreversible atomic rearrangement was probed. Our findings clearly validate that the serration is a temporally inhomogeneous characteristic of such rearrangements and not directly dependent on the resultant shear-banding spatiality. The unique spatiotemporal distribution of shear banding during nanoindentation is highlighted in terms of the potential energy landscape (PEL) theory.

An optimal theory for an expansion of flow quantities to capture the flow structures

An optimal theory on how database analysis to capture the flow structures has been developed in this paper, which include the POD method as its special case. By means of the remainder minimization method in the Sobolev space, for more general optimal conditions the new theory has the potential to overcome an inherent limitation of the POD method, i.e., it cannot be used to the situations in which the optimal condition is other than the inner product global one. As an example, using the new theory, the database of a two-dimensional flow over a backward-facing step is analyzed in detail, with velocity and vorticity bases.

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.

The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

Air flow over a mountain and the convective boundary-layer

In this paper, the analytical model coupling the convective boundary layer (CBL) with the free atmosphere developed by Qi and Fu (1992) is improved. And by this improved model, the interaction between airflow over a mountain and the CBL is further discussed. The conclusions demonstrate: (1) The perturbation potential temperatures in the free atmosphere can counteract the effect of orographic thermal forcing through entraining and mixing in the CBL. If u(M)BAR > u(F)BAR, the feedback of the perturbation potential temperatures in the free atmosphere is more important than orographic thermal forcing, which promotes the effect of interfacial waves. If u(M)BAR < u(F)BAR, orographic thermal forcing is more important, which makes the interfacial height and the topographic height identical in phase, and the horizontal speeds are a maximum at the top of the mountain. (2) The internal gravity waves propagating vertically in the free atmosphere cause a strong downslope wind to become established above the lee slope in the CBL and result in the hydraulic jump at the top of the CBL. (3) With the CBL deepening, the interfacial gravity waves induced by the potential temperature jump at the top of the CBL cause the airflow in the CBL to be subcritical.

Computation of turbulent-flow in a square duct - aspects of the secondary flow and its origin

A numerical study of turbulent flow in a straight duct of square cross-section is made. An order-of-magnitude analysis of the 3-D, time-averaged Navier-Stokes equations resulted in a parabolic form of the Navier-Stokes equations. The governing equations, expressed in terms of a new vector-potential formulation, are expanded as a multi-deck structure with each deck characterized by its dominant physical forces. The resulting equations are solved using a finite-element approach with a bicubic element representation on each cross-sectional plane. The numerical integration along the streamwise direction is carried out with finite-difference approximations until a fully-developed state is reached. The computed results agree well with other numerical studies and compare very favorably with the available experimental data. One important outcome of the current investigation is the interpretation analytically that the driving force of the secondary flow in a square duct comes mainly from the second-order terms of the difference in the gradients of the normal and transverse Reynolds stresses in the axial vorticity equation.