近边界三维水下爆炸气泡动态特性研究

模拟了近壁、近自由面的水下爆炸气泡的非线性动态特性,假定水下爆炸气泡脉动阶段的流场是无旋、不可压缩的,采用高阶曲面三角形单元离散三维气泡表面,用边界积分法求解气泡的运动,在计算奇异积分时通过重新构造双层位势的主值积分消除双层奇异积分的奇异性,得到更精确的结果,并通过合理的加权方法精确的求解边界面上各节点的真实速度,结合弹性网格技术（elastic mesh technique, EMT）得到优化速度,在整个模拟过程中不需要采用数值光顺。将本文的三维模型与轴对称模型进行的对比分析表明,两种模型计算结果吻合很好,并用三维模型模拟了气泡与自由表面及圆筒的相互作用,水下爆炸气泡在自由表面及圆筒的联合作用下呈现出强非线性。     

Dynamic stress intensity factors around two rectangular cracks in an infinite elastic plate under impact load

A dynamic problem for two equal rectangular cracks in an infinite elastic plate is considered. The two cracks are placed perpendicular to the plane surfaces of the plate. An incoming shock tensile stress is returned by the cracks. In the Laplace transform domain, the boundary conditions at the two sides of the plate are satisfied using the Fourier transform technique. The mixed boundary conditions are reduced to dual integral equations. Crack displacement is expanded in a series of functions which are zero outside of the cracks. The unknown coefficients in the series are determined by the Schmidt method. The stress intensity factors are defined in the Laplace transform domain and these are inverted using a numerical method.     

One-dimensional linear elastic waves at moving property interface

Smart materials exhibit time-varying properties while time-varying external field is applied. To investigate the one-dimensional (1-D) homogeneous time-varying properties, a moving property interface (MPI) model is proposed, and the propagation of linear elastic waves at 1-D MPI is studied in this paper. Based on the idea of weak solutions and an infinity approximation, a novel method to deal with the difficulties in using the continuities to study the waves at MPI is also proposed. Some interesting phenomena are revealed: (i) besides wave impedance, the property interface motion and wave velocity are also very important factors that influence the wave propagation; (ii) at MPI, there may exist shock waves; (iii) the property interface motion has a significant impact on the wave frequency and energy. This research provides a theoretical viewpoint in the study of smart materials with a time-dependent mechanical properties at different loading conditions.     

A rising bubble in a polymer solution

We propose a boundary integral method to study the shape of a bubble rising under gravity in a dilute polymer solution. Constitutive properties are modelled using a FENE model [M.D. Chilcott, J.M. Rallison, J. Non-Newtonian Fluid Mech. 29 (1988) 381] with a pure surface tension interface. We employ a birefringent strand representation [O.G. Harlen, J.M. Rallison, M.D. Chilcott, High-Deborah-number flows of dilute polymer, J. Non-Newtonian Fluid Mech.34 (1990) 319–349] of the wake to simulate the shape and the time-dependent motion of the bubble. Steady and non-steady solutions reproduce qualitatively the bubble deformation seen in experiment with a small region of very high curvature near the rear stagnation point of the bubble. We find a limit point for steady axisymmetric solutions if the polymer concentration is increased or the surface tension is decreased. Rise speed jump discontinuities were not found.     

General solution for dynamical problem of infinite beam on an elastic foundation

Based on the theory of Euler-Bernoulli beam and Winkler assumption for elasticfoundation,a mathematical model is presented.By using Fourier transformation for spacevariable,Laplace transformation for time variable and convolution theorem for theirinverse transformations,a general solution for dynamical problem of infinite beam on anelastic foundation is obtained.Finally,the cases of free vibration,impulsive response andmoving load are also discussed.     

Benchmark computations of diffuse interface models for two‐dimensional bubble dynamics

Diffuse interface models for incompressible two‐phase flow with large density ratios are tested on benchmark configurations for a two‐dimensional bubble rising in liquid columns. The benchmark quantities circularity, center of mass, and mean rise velocity are compared with reference solutions from Hysing et al. Copyright © 2011 John Wiley & Sons, Ltd.     

Initial kinking of an interface crack between two elastic media

The initial kinking of a thin fracture process zone near the crack tip under plane strain is studied using the Wiener-Hopf method. The crack is located at the interface between dissimilar elastic media. The fracture process zone is modeled by a straight line of normal displacement discontinuity emerging from the crack tip at an angle to the interface. The angle between the process zone and the interface is determined from the condition of strain energy maximum in the process zone. The dependences of the length and angle of the process zone on the external load and other parameters of the problem are studied. The results are compared with theoretical and experimental data obtained by other researchers __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 10, pp. 28–41, October 2007.     

Scattering of harmonic elastic waves by a plane interface crack with linear adhesive tips in a layered half space

In this paper, the scattering of elastic waves by an interface crack with linear adhesive tips in a layered half space is considered. By use of integral transform and integral equation methods, the singular integral equations of this problem are derived, which are transformed into a set of algebraic equations by means of contour integration and Chebyshev polynomials expanding technique. The numerical results of the adhesive region and stress amplitudes are given in this paper.     

On certain bounds for the in-plane translational stiffness of a disc inclusion at a bi-material elastic interface

This paper presents a set of bounds that can be used to estimate the in-plane translational stiffness of a rigid circular disc inclusion that is embedded at the interface between two dissimilar elastic half-space regions.     

Dynamics of a cavitation bubble near a solid surface and the induced damage

Numerical and experimental studies of the dynamics of a cavitating bubble near a resilient metal surface were performed. To augment the experimental flow visualizations of a collapsing bubble, numerical simulations were conducted to more thoroughly identify the collapse dynamics and analyze the flow. A bubble collapse was captured using a high-speed camera and back illumination. The metal sample was made of pure aluminum placed near a collapsing cavitation bubble at various distances from the metal surface. Width, depth, and volume of the induced material deformations were measured using an optical microscope and a three-dimensional profilometer and then compared against existing experimental data from the literature. The cavitating bubble’s dynamics and the related flow were simulated numerically using the open source finite volume based flow solver CavitatingFOAM. This code solved the Navier–Stokes equations for compressible two-phase flows using an Euler–Euler approach, including the barotropic equations of state. Bubble shapes, collapse times, and obtained damage parameters were compared to experimental observations. Impact velocities, pressures, shear rates, and various flow phenomena were discussed, providing broad insight into bubble dynamics and the induced damage.     

Anti-plane waves near an interface between two piezoelectric half-spaces

We show the existence of certain new waves that can propagate near an interface between two half-spaces of different piezoelectric ceramics, where the interface is modeled by a membrane with the surface/interface elasticity [1]. The current configuration can be reduced into a number of well-known results as special cases, such as Love wave, Bleustein and Gulyaev wave. Together with our previous work for the imperfect interface [2], a full range of consideration of the interface affecting the anti-plane waves is now completed.     

Level set method for numerical simulation of a cavitation bubble,its growth,collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier–Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving interface is captured by the level set function, and the interface velocity is resolved by “one-side” velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the “shock wave”-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating. The project supported by the National Natural Science Foundation of China (10272032 and 10672043). The English text was polished by Keren Wang.     

自由场水中爆炸气泡的物理特性

将水中爆炸气泡运动阶段周围流场假设为无粘、无旋、不可压缩的理想流体,运用边界元法模拟自由场中气泡的运动,在气泡运动模拟过程中引入数值光顺技术及弹性网格技术,避免因网格扭曲而导致的数值发散,并开发计算程序。计算值与实验值吻合良好,误差小于10%。从自由场水中爆炸气泡的基本现象入手,基于本文中开发的程序系统地研究了自由场中气泡的动力学特性。对流场中不同方位的压力进行分析,得出气泡中心的迁移方向及射流的攻击方向压力载荷比其他方向均大,说明气泡射流的攻击方向压力载荷最大,对水中结构造成严重毁伤,表明了气泡载荷的不对称性。计算了流场中不同位置的速度变化曲线,结果表明随着距气泡中心距离的增大,气泡运动引起的滞后流的速度迅速减小,且随着气泡的膨胀和坍塌,滞后流的方向逆转,总结了滞后流的衰减及变化规律。     

Acoustic radiation force on an elastic cylinder in a Gaussian beam near an impedance boundary

Acoustic radiation force (ARF) is studied by considering an infinite elastic cylinder near an impedance boundary when the cylinder is illuminated by a Gaussian beam. The surrounding fluid is an ideal fluid. Using the method of images and the translation-addition theorem for the cylindrical Bessel function, the resulting sound field including the incident wave, its reflection from the boundary, the scattered wave from the elastic cylinder, and its image are expressed in terms of the cylindrical wave function. Then, we deduce the exact equations of the axial and transverse ARFs. The solutions depend on the cylinder position, cylinder material, beam waist, reflection coefficient, distance from the impedance boundary, and absorption in the cylinder. To analyze the effects of the various factors intuitively, we simulate the radiation force for non-absorbing elastic cylinders made of stainless steel, gold, and beryllium as well as for an absorbing elastic cylinder made of polyethylene, which is a well-known biomedical polymer. The results show that the impedance boundary, cylinder material, absorption in the cylinder, and cylinder position in the Gaussian beam significantly affect the magnitude and direction of the force. Both stable and unstable equilibrium regions are found. Moreover, a larger beam waist broadens the beam domain, corresponding to non-zero axial and transverse ARFs. More importantly, negative ARFs are produced depending on the choice of the various factors. These results are particularly important for designing acoustic manipulation devices operating with Gaussian beams.     

Static and dynamic snap-through behaviour of an elastic spherical shell

The deformation and snap-through behaviour of athin-walled elastic spherical shell statically compressed on aflat surface or impacted against a flat surface are studied theoretically and numerically in order to estimate the influenceof the dynamic effects on the response.A table tennis ballis considered as an example of a thin-walled elastic shell.Itis shown that the increase of the impact velocity leads to avariation of the deformed shape thus resulting in larger deformation energy.The increase of the contact force is causedby both the increased contribution of the inertia forces andcontribution of the increased deformation energy.The contact force resulted from deformation/inertia ofthe ball and the shape of the deformed region are calculated by the proposed theoretical models and compared withthe results from both the finite element analysis and somepreviously obtained experimental data.Good agreement isdemonstrated.     

Dynamics of a small bubble in a compressible fluid

An efficient (simplified) method for solving problems of spherically symmetric dynamics of a small gas bubble in a compressible fluid is proposed. The method is based on the joint use of the full problem statement (the gas dynamics equations for the gas and the fluid) and its relevant simplifications. Some approximate statements are discussed. In the proposed method, the rarefaction and compression of the gas during the slow motion of the bubble surface is assumed to be uniform over the bubble volume. At the same time the fluid in the thin zone adjacent to the bubble is considered to be slightly compressible. Otherwise the gas dynamics equations are used for the gas and the fluid. The dynamics of the fluid in the thick external zone are described by the linear acoustics only. The proposed simplified method and two others used in literature are estimated by comparison of their numerical results with those obtained in full statement. Copyright © 2000 John Wiley & Sons, Ltd.     

The asymptotic limit of an eigenvalue problem related to the buckling of rolled elastic strips

We examine the structure of the marginal stability curves of an eigenvalue problem related to the buckling deformations observed during cold rolling of sheet metal. The instability in question is characterised by a centre “wave” pattern and arises as the interplay between the self-equilibrating residual stresses associated with the rolling process, on the one hand, and the traction force acting on the strip, on the other. When the latter effect dominates, we show that singular perturbation methods can be used to unravel a number of novel mathematical features of the linear bifurcation equation. We also provide simple quantitative formulae that facilitate an easy interpretation of the corresponding physical phenomena.