多洞室对矢量波散射引起半空间表面位移的边界元解

采用边界元方法研究了半空间中近表面多洞室对矢量波的散射问题,给出了以全空间格林函数为基本解且半空间表面离散的边界积分方程,在这一边界各分方程中,较好的消除了主值积分,在半空间表面进行离散时,采用无限单元与有限单元相结合的方法,大大减少了计算量,提出了精度。     

Simulation of non‐linear free surface motions in a cylindrical domain using a Chebyshev–Fourier spectral collocation method

When a liquid is perturbed, its free surface may experience highly non‐linear motions in response. This paper presents a numerical model of the three‐dimensional hydrodynamics of an inviscid liquid with a free surface. The mathematical model is based on potential theory in cylindrical co‐ordinates with a σ‐transformation applied between the bed and free surface in the vertical direction. Chebyshev spectral elements discretize space in the vertical and radial directions; Fourier spectral elements are used in the angular direction. Higher derivatives are approximated using a collocation (or pseudo‐spectral) matrix method. The numerical scheme is validated for non‐linear transient sloshing waves in a cylindrical tank containing a circular surface‐piercing cylinder at its centre. Excellent agreement is obtained with Ma and Wu's [Second order transient waves around a vertical cylinder in a tank. Journal of Hydrodynamics 1995; Ser. B4 : 72–81] second‐order potential theory. Further evidence for the capability of the scheme to predict complicated three‐dimensional, and highly non‐linear, free surface motions is given by the evolution of an impulse wave in a cylindrical tank and in an open domain. Copyright © 2001 John Wiley & Sons, Ltd.     

Refraction effects in the generation of helical surface waves on a cylindrical obstacle

We investigate the scattering of compressional waves from an infinite, circular-cylindrical obstacle, and the excitation during the scattering process of surface waves that propagate along helical paths over the cylinder surface. For the case of a rigid or soft obstacle, the surface waves are external, and are obtained via the use of a Watson transformation. For the case of a penetrable cylinder, additional internal, resonant surface waves are generated for which the phase and group velocity dispersion curves can be obtained from the Resonance Scattering Theory. We perform a detailed study of certain refraction effects which take place upon the generation of the surface waves by the incident plane wave.     

Inertial waves in a rotating annulus with inclined inner cylinder: comparing the spectrum of wave attractor frequency bands and the eigenspectrum in the limit of zero inclination

We investigate theoretically inertial waves inside a liquid confined between two co-rotating coaxial cylinders of finite length. We consider the case of small viscosity and high angular velocity (i.e., small Ekman numbers), a parameter range of interest for many geophysical applications. In this case, inertial waves propagating in the container show multiple reflections at the walls before the waves can be damped by weak diffusion. We allow for the inner cylinder wall to be parallel or inclined with respect to the annulus’ vector of rotation (truncated cone). For the limit of zero viscosity, the wave propagation is governed by a boundary value problem that is composed of a linear second-order hyperbolic partial differential equation and the impermeability boundary conditions. For the special case of vertical cylinder walls (no inclination of the inner cylinder), this boundary value problem is separable, the corresponding eigenmodes can analytically be found and they are regular. However, when the inner cylinder wall is inclined, the hyperbolicity of the governing equation leads to internal shear layers (corresponding to singularities for the inviscid case). The geometrical structure of the shear layers can be explained by inertial waves, trapped on limit cycles denoted as wave attractors. The shape of the limit cycles depends on the wave frequency. In fact, the spectrum of regular modes, existing for the case of vertical cylinder walls, vanishes almost completely when the inner wall is inclined. Instead of a spectrum of discrete frequencies and regular eigenmodes, a spectrum of wave attractor frequency bands and singular eigenmodes exist. The question addressed here is whether the spectrum of wave attractor intervals collapses to the discrete frequency spectrum when the inclination angle of the inner cylinder goes to zero. To answer this question, the attractor frequency intervals are evaluated numerically for a series of decreasing cylinder inclination angles and are compared with the analytically found eigenspectrum for the case of zero inclination. Goal is to better understand the asymptotic behavior of the problem for decreasing inclination angles. This understanding helps to interpret results from laboratory experiments with geometries that differ from the perfect annulus with parallel cylinder walls.     

Surface wave scattering by a material filled rectangular impedance groove in a reactive plane

The problem of TE-polarized surface wave scattering from a rectangular impedance groove located on an infinite reactive plane which is filled with dielectric material is considered for a rather general case where the impedances of the horizontal and vertical sides of the groove have different values. The multiple interactions up to the second order between the edges of the groove are obtained to yield diffracted field. The diffraction problem is first reduced into a modified Wiener-Hopf equation and then solved approximately. The solution contains branch-cut integrals and two infinite sets of constants satisfying two infinite systems of linear algebraic equations. The approximate analytical evaluations of the corresponding integrals as well as the numerical solutions of the linear algebraic equation systems are obtained for various values of the parameters such as the surface reactance of the guiding plane, the vertical and horizontal wall impedances of the groove, the permittivity of the material loading, the width and the height of the groove which permit one to study the effect of these parameters on the diffraction phenomenon.     

Nonlinear azimuthal waves in a centrifuge

Azimuthal wave motions in a liquid which partially fills a cylinder (centrifuge) rapidly rotating about a horizontal axis are discussed in this paper. Under the action of centrifugal force the liquid is pressed to the wall of the cylinder and moves together with it about the central air core. The vibrations of the free surface which arise are called centrifugal waves [1]. The difficulties of their theoretical investigation are related to the nonlinearity both of the basic equations and also of the boundary condition for the pressure on the free surface; therefore they have previously been studied only by linear methods [1, 2]. Nonlinear azimuthal waves in a centrifuge with an infinite radius of the rotating cylinder are analytically described below. The waves found are an analog of Gerstner trochoidal waves on a cylindrical surface. An approximate solution for a centrifuge with a finite outer radius is constructed by matching the waves obtained to the known linear ones.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 86–89, May–June, 1984.In conclusion the author expresses his gratitude to E. I. Yakubovich for useful discussion.     

Fully non‐linear free‐surface simulations by a 3D viscous numerical wave tank

A finite difference scheme using a modified marker‐and‐cell (MAC) method is applied to investigate the characteristics of non‐linear wave motions and their interactions with a stationary three‐dimensional body inside a numerical wave tank (NWT). The Navier–Stokes (NS) equation is solved for two fluid layers, and the boundary values are updated at each time step by a finite difference time marching scheme in the frame of a rectangular co‐ordinate system. The viscous stresses and surface tension are neglected in the dynamic free‐surface condition, and the fully non‐linear kinematic free‐surface condition is satisfied by the density function method developed for two fluid layers. The incident waves are generated from the inflow boundary by prescribing a velocity profile resembling flexible flap wavemaker motions, and the outgoing waves are numerically dissipated inside an artificial damping zone located at the end of the tank. The present NS–MAC NWT simulations for a vertical truncated circular cylinder inside a rectangular wave tank are compared with the experimental results of Mercier and Niedzwecki, an independently developed potential‐based fully non‐linear NWT, and the second‐order diffraction computation. Copyright © 1999 John Wiley & Sons, Ltd.     

Axisymmetric longitudinal wave propagation in a finite prestrained circular cylinder embedded in a finite prestrained compressible medium

Axisymmetric longitudinal wave propagation in a finite prestrained circular cylinder contained in a finite prestrained infinite body is investigated within the scope of the piecewise-homogeneous body model by employing the three-dimensional linearized theory of elastic waves in a prestressed body. It is assumed that the materials of the cylinder and infinite body are compressible and that their elastic relations are described by a harmonic potential. Numerical results are presented and discussed for the case where the elastic constants of the cylinder are greater than those of the surrounding infinite body     

Water wave interaction with a sphere in a two-layer fluid flowing through a channel of finite depth

Using linear water wave theory, we consider a three-dimensional problem involving the interaction of waves with a sphere in a fluid consisting of two layers with the upper layer and lower layer bounded above and below, respectively, by rigid horizontal walls, which are approximations of the free surface and the bottom surface; these walls can be assumed to constitute a channel. The effects of surface tension at the surface of separation is neglected. For such a situation time-harmonic waves propagate with one wave number only, unlike the case when one of the layers is of infinite depth with the waves propagating with two wave numbers. Method of multipole expansions is used to find the particular solutions for the problems of wave radiation and scattering by a submerged sphere placed in either of the upper or lower layer. The added-mass and damping coefficients for heave and sway motions are derived and plotted against various values of the wave number. Similarly the exciting forces due to heave and sway motions are evaluated and presented graphically. The features of the results find good agreement with previously available results from the point of view of physical interpretation.     

Scattering of implusive sound waves by a rigid cylinder

Summary Analytic time solutions for the scattering of impulsive waves (the time-space Green function) by a rigid circular cylinder in an annular domain are obtained through a finite integral transform of the spatial variable. By a similar generalized procedure the scattering of impulsive waves by a rigid cylinder in an infinite medium is described in order to obtain the Green function of the reduced wave equation in terms of a series of propagation modes.

---

Sommario In questo articolo si presentano alcuni metodi di soluzione analitica di problemi di diffrazione di onde impulsive in un dominio anulare. Le soluzioni sono ottenute con l'impiego di una trasformata integrale finita della variabile spaziale. Attraverso un procedimento analogo generalizzato viene descritta la diffrazione di onde impulsive causate da un cilindro circolare rigido in un mezzo infinito.

---

This work was supported by C.N.R., Committee for Mathematical Sciences.     

船行波对水中直立圆柱作用的数值模拟

针对船行波对水中直立圆柱的作用进行了数值计算。主要包括船行波的计算以及波浪对水中直立圆柱的作用。对于船行波的计算应用薄船理论和Noblesse给出的格林函数的简化形式,计算出Wigley船型在静水中匀速直线运动产生的船行波的远场解。以该船行波为入射波,应用Rankine源方法对船行波对水中直立圆柱的作用力及力矩进行了数值计算,得出了船只经过圆柱周围时圆柱的受力情况,对计算结果进行了定性的分析。     

Internal Waves Excited by a Moving Source in a Medium of Variable Buoyancy

The problem of the far field of internal gravity waves generated by an oscillating point perturbation source moving in a vertically infinite layer of a stratified medium of variable buoyancy is considered. The analytical solution of the problem is obtained by two ways for a model quadratic buoyancy frequency distribution. In the first case the solution is expressed in terms of the eigenfunctions of the vertical spectral problem and the Hermite polynomials. In the second case the solution in the form of the Green’s characteristic function is represented in terms of the functions of parabolic cylinder. The analytical solutions obtained make it possible to describe the amplitudephase characteristics of the far fields of internal gravity waves in a stratified medium with variable Brunt-Väisäläfrequency.     

Linear dispersive shear waves in two-layer elastic medium

In linear isotropic elasticity, cases exist where pure shear waves are possible even in bounded media; these have a nondispersive mode with propagation speed equal to that in an infinite medium. For these cases, consideration of a two-layer medium shows the existence of dispersion which vanishes only with equality of the propagation speeds. The present study uses the method of asymptotic expansions; a uniformly valid approximation is obtained to describe the speed and dispersive nature of these waves. The elegance of this approach is brought out by derivation of the basic result of this study viz., Jeffrey equation to describe the farfield structure of these waves.This work was done towards partial fulfillment of the requirements for Ph. D at Iowa State University, Ames, IOWA, in 1972.     

LOCAL AND GLOBAL INSTABILITY OF FLUID-CONVEYING PIPES ON ELASTIC FOUNDATIONS

We investigate the relationship between the local and global bending motions of fluid-conveying pipes on an elastic foundation. The local approach refers to an infinite pipe without taking into account its finite ends, while in the global approach we consider a pipe of finite length with a given set of boundary conditions. Several kinds of propagating disturbances are identified from the dispersion relation, namely evanescent, neutral and unstable waves. As the length of the pipe is increased, the global criterion for instability is found to coincide with local neutrality, whereby a local harmonic forcing only generates neutral waves. For sets of boundary conditions that give rise only to static instabilities, the criterion for global instability of the long pipe is that static neutral waves exist. Conversely, for sets of boundary conditions that allow dynamic instabilities, the criterion for global instability of the long pipe corresponds to that for the existence of neutral waves of finite nonzero frequency. These results are discussed in relation with the work of Kulikovskii and other similar approaches in hydrodynamic stability theory.     

各向异性介质的弹性波散射问题的边界元方法

本文引用加权残数法建立了各向异性介质内含任意形式异质夹杂时的散射问题的边界积分方程式,导出了相应的辐射条件,计算了内含圆柱体,椭圆柱体、界面裂纹情形下对SH 波的散射位移场、应力场以及散射横截面.数值结果表明本方法用于解答各向异性介质的弹性波散射问题具有良好的精度和应用前景.     

An explicit formulation for the evolution of nonlinear surface waves interacting with a submerged body

An explicit formulation to study nonlinear waves interacting with a submerged body in an ideal fluid of infinite depth is presented. The formulation allows one to decompose the nonlinear wave–body interaction problem into body and free‐surface problems. After the decomposition, the body problem satisfies a modified body boundary condition in an unbounded fluid domain, while the free‐surface problem satisfies modified nonlinear free‐surface boundary conditions. It is then shown that the nonlinear free‐surface problem can be further reduced to a closed system of two nonlinear evolution equations expanded in infinite series for the free‐surface elevation and the velocity potential at the free surface. For numerical experiments, the body problem is solved using a distribution of singularities along the body surface and the system of evolution equations, truncated at third order in wave steepness, is then solved using a pseudo‐spectral method based on the fast Fourier transform. A circular cylinder translating steadily near the free surface is considered and it is found that our numerical solutions show excellent agreement with the fully nonlinear solution using a boundary integral method. We further validate our solutions for a submerged circular cylinder oscillating vertically or fixed under incoming nonlinear waves with other analytical and numerical results. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of the theory of generalized rays to diffractions of transient waves by a cylinder

The theory of generalized rays was developed to analyze transient waves in layered media where incident circular or spherical waves are reflected and refracted by plane boundaries. The theory has been recently extended to analyze the diffraction of transient waves by a spherical or a cylindrical boundary. In this paper, the generalized ray integrals, which represent the fourier transformed diffracted waves, are formulated for the diffraction of an incident spherical pulse by a circular cylinder. The ray integral involves a double integration with respect to two variables of wave slowness. Through a simultaneous transformation of variables, the inverse Fourier transform of these double integrals are completed by applying the Cagniard method.     

On the asymptotical behaviour of solutions of a class of nonlinear control systems

Ｃｏｎｓｉｄｅｒｔｈｅｆｏｌｌｏｗｉｎｇｎｏｎｌｉｎｅａｒｃｏｎｔｒｏｌｓｙｓｔｅｍ：·ｘ＝Ａｘ＋Ｂ（σ）,σ＝ＣＴｘ,（１）ｗｈｅｒｅｘ＝（ｘ１,ｘ２,…,ｘｎ）Ｔ∈Ｒｎ,ｔｈｅｎｏｔｅＴｅｘｐｒｅｓｓｅｓｔｈｅｔｒａｎｓｐｏｓｅｏｆｍａｔｒｉｘａｎｄｖｅｃｔｏｒ,Ａ,ＢａｎｄＣａｒｅｎ×ｎ,ｎ×ｍａｎｄｎ×ｍｃｏｎｓｔａｎｔｍａｔｒｉｃｅｓｒｅｓｐｅｃｔｉｖｅｌｙ,σ＝（σ１,σ２,…,σｍ）Ｔ∈Ｒｍａｎｄ（σ）＝（１（σ１）,２（σ２）,…,ｍ（σｍ））Ｔ：Ｒｍ→Ｒｍｉｓｃｏｎｔ…     

DYNAMIC INTERACTION OF PLANE WAVES WITH A UNILATERALLY FRICTIONALLY CONSTRAINED INCLUSION-TIME DOMAIN BOUNDARY ELEMENT ANALYSIS

A 2D time domain boundary element method (BEM) is developed to solve the transient scattering of plane waves by a unilaterally frictionally constrained inclusion. Coulomb friction is assumed along the contact interface. The incident wave is assumed strong enough so that localized slip and separation take place along the interface. The present problem is in effect a nonlinear boundary value problem since the mixed boundary conditions involve unknown intervals (slip, separation and stick regions). In order to determine the unknown intervals, an iterative technique is developed. As an example, we consider the scattering of a circular cylinder embeddedin an infinite solid.