Water wave radiation by a sphere submerged in water with an ice-cover

Using the multipoles method, we formulate the problems of radiation (both heave and sway) of water waves by a submerged sphere in deep as well as in uniform finite depth water with an ice-cover, with the ice-cover being modelled as an elastic plate of very small thickness. In each case this leads to an infinite system of linear equations which are solved numerically by standard techniques. The added-mass and damping coefficients for a heaving and swaying sphere are obtained and depicted graphically against the wave number for various values of the radius of the submerged sphere and flexural rigidity of the ice-cover to show the effect of the presence of ice-cover on these quantities. When the flexural rigidity is taken to be zero, the numerical results for the added-mass and damping coefficient for water with a free surface are recovered.     

Performance-based design and analysis of flexible composite propulsors

Advanced composite propellers, turbines, and jet engines have become increasingly popular in part because of their ability to provide improved performance over traditional metallic rotors through exploitation of the intrinsic bend–twist coupling characteristics of anisotropic composite materials. While these performance improvements can be significant from a conceptual perspective, the load-dependent deformation responses of adaptive blades make the design of these structures highly non-trivial. Hence, it is necessary to understand and predict the dependence of the deformations on the geometry, material constitution, and fluid–structure interaction responses across the entire range of expected loading conditions.The objective of this work is to develop a probabilistic performance-based design and analysis methodology for flexible composite propulsors. To demonstrate the method, it is applied for the design and analysis of two (rigid) metallic and (flexible) composite propellers for a twin-shafted naval combatant craft. The probabilistic operational space is developed by considering the variation of vessel thrust requirements as a function of the vessel speed and wave conditions along with the probabilistic speed profiles. The performance of the metallic and composite propellers are compared and discussed. The implications of load-dependent deformations of the flexible composite propeller on the operating conditions and the resulting performance with respect to propeller efficiency, power demand, and fluid cavitation are presented for both spatially uniform and varying flows. While the proposed framework is demonstrated for marine propellers, the methodology can be generally applied for any marine, aerospace, or wind energy structure that must operate in a wide range of loading conditions over its expected life.     

Phase velocity effects of the wave interaction with exponentially sheared current

The nonlinear interaction between the unidirectional bichromatic wave-train and exponentially sheared current in water of an infinite depth is investigated. The model is based on the vorticity transport equation and the exact free surface conditions, without any assumptions for the existence of small physical parameters. Earlier works of the wave–current interaction were mainly restricted to either current acted on the monochromatic wave or irregular waves limited to irrotational current. Different from these previous works, no constraint is made in our model for amplitudes of the primary wave, and the current owns an exponential type profile along the vertical line. To ensure that the effect of vorticity on the phase velocity is consistent with earlier derivation, the case of a small amplitude wave traveling on the exponentially sheared current is examined firstly. Then the effect of nonlinearity on the phase velocity of primary waves in a bichromatic wave-train is considered. Accurate high-order approximations of the phase velocity are obtained under consideration of both the nonlinear wave self–self and mutual interactions. Finally, the combined effect of vorticity and nonlinearity on the phase velocity is investigated through the case of a bichromatic wave-train propagating on an exponentially sheared current. It is found that the characteristic current slope determines the effect of vorticity on the phase velocity caused by nonlinear wave self–self and mutual interactions, and the surface current strength may amplify/reduce this effect.     

Numerical simulation of a solitary wave interaction with submerged multi-bodies

The problems of a solitary wave passing over rectangular cylinders have been analysed. The numerical simulation is based on the full nonlinear two-dimensional Navier-Stokes equations which are solved by the finite difference method. The free surface is dealt with by the Volume of Fluid method (VOF). Results for a solitary wave passing over a single cylinder are compared with the experimental data of Seabra-Santos, Penouard and Temperville^[2] and better agreement is obtained than those obtained from the long wave equation based on the potential flow theory. Results are also given for two cylinders with different gaps. The project supported by the National Natural Science Foundation of China and the Development Foundation of Science and Technology of Shanghai Education Committee and the Royal Society.     

Vortex shedding induced by a solitary wave propagating over a submerged vertical plate

Experimental study was conducted on the vortex shedding process induced by the interaction between a solitary wave and a submerged vertical plate. Particle image velocimetry (PIV) was used for quantitative velocity measurement while a particle tracing technique was used for qualitative flow visualization. Vortices are generated at the tip of each side of the plate. The largest vortices at each side of the plate eventually grow to the size of the water depth. Although the fluid motion under the solitary wave is only translatory, vortices are shed in both the upstream and downstream directions due to the interaction of the generated vortices as well as the vortices with the plate and the bottom. The process can be divided into four phases: the formation of a separated shear layer, the generation and shedding of vortices, the formation of a vertical jet, and the impingement of the jet onto the free surface. Similarity velocity profiles were found both in the separated shear layer and in the vertical jet.     

Analytical study for oblique wave interaction with a submerged horizontal perforated plate near a partially reflecting vertical wall

Meccanica - A novel iterative analytical solution is developed to study the oblique wave interaction with a horizontally submerged perforated plate near a partially reflecting vertical wall. The...     

Shock wave interaction with a cloud of particles

The present paper is devoted to experimental and theoretical investigation of the shock wave (SW) propagation in a mixture of gas and solid particles in the presence of explicit boundaries of the two-phase region (cloud of particles). The effect of the qualitative change in the supersonic flow behind the SW in a cloud of particles within the range of the volume concentration of the disperse phase 0.1-3% is experimentally shown and theoretically grounded. Received 15 April 1996 / Accepted 3 June 1996     

Generic non-linear modelling of a bi-material composite beam with partial shear interaction

Composite members composed of two materials joined by shear connection find widespread use in engineering infrastructure, in both traditional practice and innovative applications. Studies in the literature dating back nearly 60 years have elucidated the mechanics of the behaviour of these composite structural members in which the solution for the slip at the interface between the materials was determined by solving a linear differential equation. However, these solutions are based on a linear formulation of the strain-displacement relationship, and in some applications this relationship must be represented in non-linear form, so that the second order effects in the member can be quantified correctly. This paper presents such a study for a composite member with two materials, being typical of a steel-concrete composite beam in structural engineering. It quantifies the restraint of the member ends by longitudinal and rotational elastic springs, so that the axial tension developed is a function of the transverse loading, material properties, cross-sectional properties and the restraint stiffness. The problem is treated using minimisation of the total potential stored in the two members, the elastic shear connection at their interface, the restraints at the ends and the work done by the transverse forces, for which the differential equations for the deformations can be determined from routine variational calculus. The non-linear equation of equilibrium relating the external loading to the internal actions is stated in closed form by invoking the static and kinematic boundary conditions for the member. The solution is compared with closed form treatments derived elsewhere, and a representative member is analysed so that the influences of the non-linearity, end restraint stiffness and degree of partial shear interaction on its behaviour can be examined.     

Elastic field of a composite cylinder with a spatially varying dynamic eigenstrain

The analytical solutions of displacements and stresses for an eigenstrain problem in a composite bi-layered coaxial cylinder are presented in this article. The inner cylinder is assumed to undergo a dynamic, spatially varying eigenstrain. The spatial distribution of the eigenstrian is taken to be a quadratic polynomial with arbitrary coefficients along the radial direction. Furthermore, the eigenstrain is assumed to be harmonically time-dependent. Elasticity equations are constructed to directly solve the problem. The effect of spatial distribution, as well as the angular frequency of the eigenstrain on the elastic response of the composite cylinder has been illustrated graphically.     

Numerical study of formation of solitary strain waves in a nonlinear elastic layered composite material

Propagation of nonlinear strain waves through a layered composite material is considered. The governing macroscopic wave equation for the long-wave case was obtained earlier by the higher-order asymptotic homogenization method (Andrianov et al., 2013). Non-stationary dynamic processes are investigated by a pseudo-spectral numerical procedure. The time integration is performed by the Runge–Kutta method; the approximation with respect to the spatial co-ordinate is provided by the Fourier series expansion. The convergence of the Fourier series is substantially improved and the Gibbs–Wilbraham phenomenon is reduced with the help of Padé approximants. As result, we explore how fast and under what conditions the solitary strain waves can be generated from an initial excitation. The numerical and analytical solutions (when the latter can be obtained) are in good agreement.     

A study of shock wave interaction with a rotating cylinder

Shock wave reflection over a rotating circular cylinder is numerically and experimentally investigated. It is shown that the transition from the regular reflection to the Mach reflection is promoted on the cylinder surface which rotates in the same direction of the incident shock motion, whereas it is retarded on the surface that rotates to the reverse direction. Numerical calculations solving the Navier-Stokes equations using extremely fine grids also reveal that the reflected shock transition from RR MR is either advanced or retarded depending on whether or not the surface motion favors the incident shock wave. The interpretation of viscous effects on the reflected shock transition is given by the dimensional analysis and from the viewpoint of signal propagation.Received: 24 April 2002, Accepted: 16 August 2002, Published online: 25 March 2003     

Experimental study of the response of a flexible plate to a harmonic forcing in a flow

Most aquatic animals propel themselves by flapping flexible appendages. To gain insight into the effect of flexibility on the swimming performance, we have studied experimentally an idealized system. It consists of a flexible plate whose leading edge is forced into a harmonic heave motion, and which is immersed in a uniform flow. As the forcing frequency is gradually increased, resonance peaks are evidenced on the plate's response. In addition to the forcing frequency, the Reynolds number, the plate rigidity and the forcing amplitude have also been varied. In the range of parameters studied, the main effect on the resonance is due to the forcing amplitude, which reveals that non-linearities are essential in this problem.     

Time-resolved PIV measurements of the interaction of polystyrene beads with near-wall-coherent structures in a turbulent channel flow

Time-resolved PIV measurements were performed in a dilute particle-laden flow tracking near-neutrally buoyant polystyrene beads and the velocity field of a near wall turbulent boundary layer. Data were taken in a vertical light sheet aligned in the streamwise direction at the center of a horizontal, closed loop, transparent square water channel facility. In addition, low speed measurements were performed characterizing the effects of the dispersed phase on mean and turbulence flow quantities. Reynolds shear stress slightly differed from clear water conditions whereas fluid mean and rms values were not affected. A case study for several beads revealed a clear relation between their movement and near-wall coherent structures. Several structures having 2D vorticity signatures of near-wall hairpin vortices and hairpin packets, directly affected bead movement. A statistical analysis showed that the mean streamwise velocity of ascending beads lagged behind the mean fluid velocity and bead rms values were higher than fluid ones. Particle Reynolds numbers based on the magnitude of the instantaneous relative velocity vector peaked near the wall; values not exceeding 100, too low for vortex shedding to occur. Quadrant analysis showed a clear preference for ascending beads to reside in ejections while for descending beads the preference for sweeps was less.     

A study of sloshing absorber geometry for structural control with SPH

A liquid sloshing absorber consists of a container, partially filled with liquid. The absorber is attached to the structure to be controlled, and relies on the structure's motion to excite the liquid. Consequently, a sloshing wave is produced at the liquid free-surface within the absorber, possessing energy dissipative qualities. The primary objective of this work is to numerically demonstrate the effect of a sloshing absorber's shape on its control performance. Smoothed Particle Hydrodynamics (SPH) is used to model fluid–structure interaction of the structure/sloshing absorber system in two dimensions. The structure to be controlled is a lightly damped single degree-of-freedom structure. The structure is subjected to a transient excitation and then allowed to respond dynamically, coming to rest either due to its own damping alone or with the added control of the sloshing absorber. It is identified that the control performance of the conventionally used rectangular container geometry can be improved by having inward-angled walls. This new arrangement is robust, and of significant advantage in situations when the external disturbance is of uncertain magnitude.     

Numerical investigation of the interaction of a vortex dipole with a deformable plate

Energy harvesting from coherent fluid structures is a current research topic due to its application in the design of small self-powered sensors for underwater applications. The impact of a vortex dipole with a deformable cantilevered plate at the plate tip is herein studied numerically using a strongly coupled staggered fluid–structure interaction algorithm. Three dipole Reynolds numbers, Re=500, 1500, and 3000, are investigated for constant plate properties. As the dipole approaches the plate, positive vorticity is produced on the impact face, while negative vorticity is generated at the tip of the plate. Upon impact, the dipole splits into two, and two secondary dipoles are formed. The circulation and, therefore, the trajectories of these dipoles depend on both the Reynolds number and the elasticity of the plate, and these secondary dipoles may return for subsequent impacts. While the maximum deflection of the plate does not depend significantly on Reynolds number, the plate response due to subsequent impacts of secondary dipoles does vary with Reynolds number. These results elucidate the strong interdependency between plate deformation and vortex dynamics, as well as the effect of Reynolds number on both.     

Experimental study of extreme thrust on a tidal stream rotor due to turbulent flow and with opposing waves

Time-varying thrust has been measured on a rotor in shallow turbulent flow at laboratory scale. The onset flow has a turbulence intensity of 12% at mid depth and a longitudinal turbulence length scale of half the depth, about 5 times the vertical scale, typical of shallow flows. The rotor is designed to have thrust and power coefficient variations with tip speed ratio close to that of a full-scale turbine. Three extreme probability distributions give similar thrust exceedance values with the Type 1 Pareto in mid range which gives 1:100, 1:1000 and 1:10 000 exceedance thrust forces of 1.38, 1.5 and 1.59 times the mean value. With opposing waves superimposed the extreme thrust distribution has a very similar distribution to the turbulent flow only. Exceedance forces are predicted by superposition of a drag force with drag coefficient of 2.0 based on the wave particle velocity only and with an unchanged mean thrust coefficient of 0.89. These values are relevant for the design of support structures for marine turbines.     

Experimental and numerical study on a laminar fluid-structure interaction reference test case

With the rapid development of numerical codes for fluid-structure interaction computations, the demand for validation test cases increases. In this paper we present a comparison between numerical and experimental results for such a fluid-structure interaction reference test case. The investigated structural model consists of an aluminum front cylinder with an attached thin metal plate and a rear mass at the trailing edge. All the structure is free to rotate around the axle mounted in the center of the front cylinder. The model's geometry and mechanical properties are chosen in such a way as to attain a self-exciting periodical swiveling movement when exposed to a uniform laminar flow. Reproducibility of the coupled fluid-structure motion is the key criterion for the selection of the model in order to permit an accurate reconstruction of the results in the time-phase space. The Reynolds number of the tests varies up to 270 and within that range the structure undergoes large deformations and shows a strong nonlinear behavior. It also presents two different self-excitation mechanisms depending on the flow velocity. Hence, challenging tasks arise for both the numerical solution algorithm and the experimental measurements. To account for the two different excitation mechanisms observed on increasing the speed of the flow, results for two different velocities are considered: the first at 1.07 m/s (Re=140) and the second at 1.45 m/s (Re=195). The comparisons presented in this paper are carried out on the basis of the time trace of the front body angle, trailing edge coordinates, structure deformation and the time-phase resolved flow velocity field. They reveal very good agreement in some of the fluid-structure interaction modes whereas in others deficiencies are observed that need to be analyzed in more detail.     

Finite elements based on a first-order shear deformation moderate rotation shell theory with applications to the analysis of composite structures

The first-order shear deformation moderate rotation shell theory of Schmidt and Reddy [R. Schmidt and J. N. Reddy, J. Appl. Mech. 55, 611–617 (1988)] is used as a basis for the development of finite element models for the analysis of the static, geometrically non-linear response of anisotropic and laminated structures. The incremental, total Lagrangian formulation of the theory is developed, and numerical solutions are obtained by using the isoparametric Lagrangian 9-node and Serendipity 8-node shell finite elements. Various integration schemes (full, selective reduced, and uniformly reduced integration) are applied in order to detect and to overcome the effects of shear and membrane locking on the predicted structural response. A number of sample problems of isotropic, orthotropic, and multi-layered structures are presented to show the accuracy of the present theory. The von Kármán-type first-order shear deformation shell theory and continuum 2D theory are used for comparative analyses.