Design and analysis of nonlinear-transformation-based broadband cloaking for acoustic wave propagation

A coordinate-transformation method can be used to design invisibility cloaks for many types of waves, including acoustic waves. The traditional method for designing a cloak depends on a transformation from a virtual space to a physical space. Previous acoustic cloaks that are mainly designed with linear-transformation-based acoustics have drawbacks that acoustic wave trajectories in the cloaks cannot be controlled and tuned. This work uses a nonlinear mapping from a ray trajectory perspective to construct acoustic cloaks with tunable non-singular material properties. Use of a ray trajectory equation is a straightforward and alternate way to study propagation characteristics of different types of waves, which allows more flexibility in controlling the waves. A broadband cylindrical cloak for acoustic waves in an inviscid fluid is realized with layered non-singular, homogeneous, and isotropic materials based on a nonlinear transformation. Some advantages and improvements of the invisibility nonlinear-transformation cloak over a traditional linear-transformation cloak are analyzed. The invisibility capability of the nonlinear-transformation cloak can be tuned by adjusting a design parameter that is shown to have influence on the acoustic wave energy flowing into the region inside the cloak. Numerical examples show that the nonlinear-transformation cloak is more effective for making a domain undetectable by acoustic waves in an inviscid fluid and shielding acoustic waves from outside the cloak than the linear-transformation cloak in a broad frequency range. The methodology developed here can be used to design nonlinear-transformation cloaks for other types of waves.     

Applications of active acoustics in particle technology

Earlier work by the authors in which active sound pressure signals and impulsive pressure disturbances were used to measure flow rates in gas solid systems was briefly reviewed. Work in progress with an emerging technology in which Helmholtz resonance is applied to the measurement of volume is outlined.     

High order methods for acoustic scattering: Coupling farfield expansions ABC with deferred-correction methods

Arbitrary high order numerical methods for time-harmonic acoustic scattering problems originally defined on unbounded domains are constructed. This is done by coupling recently developed high order local absorbing boundary conditions (ABCs) with finite difference methods for the Helmholtz equation. These ABCs are based on exact representations of the outgoing waves by means of farfield expansions. The finite difference methods, which are constructed from a deferred-correction (DC) technique, approximate the Helmholtz equation and the ABCs, with the appropriate number of terms, to any desired order. As a result, high order numerical methods with an overall order of convergence equal to the order of the DC schemes are obtained. A detailed construction of these DC finite difference schemes is presented. Additionally, a rigorous proof of the consistency of the DC schemes with the Helmholtz equation and the ABCs in polar coordinates is also given. The results of several numerical experiments corroborate the high order convergence of the novel method.     

Numerical interaction of boundary waves with perfectly matched layers in two space dimensional elastic waveguides

Perfectly matched layers (PMLs) are now a standard approach to simulate the absorption of waves in open domains. Wave propagation in elastic waveguides has the possibility to support back-propagating modes (propagating modes with oppositely directed group and phase velocities) with long wavelengths. Back-propagating modes can lead to temporally growing solutions in the PML. In this paper, we demonstrate that back-propagating modes in a two space dimensional isotropic elastic waveguide are not harmful to a discrete and finite width PML. Analysis and numerical experiments confirm the accuracy and stability of the PML.     

A boundary element method for calculating the SH waves scattering from arbitrarily shaped holes in an anisotropic medium

The scattering problem of elastic wave by arbitrarily shaped cavities in an infinite anisotropic medium is investigated by the boundary integral equation (BIE) method. The formulations of BIE are derived with the help of generalized Green's formula. The discretization of BIE is based upon constant elements. After confirmation of the accuracy of the present method, some numerical examples are given for various cavities in a full space, in which an isotropic body with a circular cylinder hole is used for comparison and good agreement is observed. It has been proved that the method developed in this paper is effective.     

Parametric interaction of acoustic waves in micro-inhomogeneous media with hysteretic nonlinearity and relaxation

A theoretical investigation of parametric processes that arise as a result of the interaction of powerful and weak longitudinal acoustic waves in micro-inhomogeneous media with hysteretic nonlinearity and relaxation was carried out. The case of degenerate interaction between a powerful high-frequency wave and a weak low-frequency one was considered. The nonlinear damping coefficient and the carrier frequency phase delay of the weak wave propagating under the action of the powerful wave were determined.     

求解二维污染扩散方程的时空任意阶的三点紧致显格式

通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间二阶导数任意精度的三点紧致的表达式,并在半高散方程中通过二维扩散方程本身把时间导数转换为空间导数,从而推导出了时空任意阶的三点紧致显格式.数值实验表明,本文格式的精度很高,而且具有使用简单,易于编程的优点,对求解二维污染扩散方程具有很好的应用前景.     

Three dimensional simulations of penetrative convection in a porous medium with internal heat sources

The problem of penetrative convection in a fluid saturated porous medium heated internally is analysed. The linear instability theory and nonlinear energy theory are derived and then tested using three dimensions simulation.Critical Rayleigh numbers are obtained numerically for the case of a uniform heat source in a layer with two fixed surfaces. The validity of both the linear instability and global nonlinear energy stability thresholds are tested using a three dimensional simulation. Our results show that the linear threshold accurately predicts the onset of instability in the basic steady state. However, the required time to arrive at the basic steady state increases significantly as the Rayleigh number tends to the linear threshold.     

Panel flutter prediction in two dimensional flow with enhanced piston theory

Piston theory may be used in the high Mach number supersonic flow region and/or in very high frequency subsonic or supersonic flow. In this flow model, the pressure at a point on the fluid-solid interface only depends on the downwash at the same point. However the classical piston theory may not be sufficient for some phenomena in aeroelasticity and aeroacoustics (far field prediction). Dowell and Bliss have created an extension of piston theory that allows for higher order effects that take into account the effect the distribution of downwash on pressure at any point. For simple harmonic motion, expansions in reduced frequency, inverse reduced frequency and/or inverse (square of) Mach number have all been created; The effects of higher order terms in these several expansion in creating an enhanced piston theory was illustrated for plunge and pitch motion of an airfoil (discrete system) by Ganji and Dowell. In the present paper, flutter prediction for a flexible panel in two –dimensional flow is investigated using enhanced piston theory. The goal of the present paper is to demonstrate that an enhance version of piston theory can analyze single degree of freedom flutter of a panel as compared to the classical piston theory and quasi-steady aerodynamic models which can only treat coupled mode flutter.     

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.     

Stress analysis in two dimensional electrostrictive material with an elliptic rigid conductor

This paper presents the governing equations of electrostrictive materials. The stress and electric field solutions for an infinite plate with a rigid elliptic conductor under applied load at infinity are given. The asymptotic expansions of the solution for a narrow elliptic conductor show that the stresses and the electric fields near the end of a narrow elliptic conductor possess r⁻¹ and r^−1/2 forms respectively in a local coordinate system with the origin at its focus.     

Comparison and analysis of two Coulomb friction models on the dynamic behavior of slider-crank mechanism with a revolute clearance joint

The objective of this study is to investigate the effects of the Coulomb dry friction model versus the modified Coulomb friction model on the dynamic behavior of the slider-crank mechanism with a revolute clearance joint. The normal and tangential forces acting on the contact points between the journal and the bearing are described by using a Hertzian-based contact force model and the Coulomb friction models, respectively. The dynamic equations of the mechanism are derived based on the Lagrange equations of the first kind and the Baumgarte stabilization method. The frictional force is solved via the linear complementarity problem (LCP) algorithm and the trial-and-error algorithm. Finally, three numerical examples are given to show the influence of the two Coulomb friction models on the dynamic behavior of the mechanism. Numerical results show that due to the stick friction, the slider-crank mechanism may exhibit stick-slip motion and can balance at some special positions, while the mechanism with ideal joints cannot.     

On the propagation waves in the theory of thermoelasticity with microtemperatures

The present paper studies the propagation of plane time harmonic waves in an infinite space filled by a thermoelastic material with microtemperatures. It is found that there are seven basic waves traveling with distinct speeds: (a) two transverse elastic waves uncoupled, undamped in time and traveling independently with the speed that is unaffected by the thermal effects; (b) two transverse thermal standing waves decaying exponentially to zero when time tends to infinity and they are unaffected by the elastic deformations; (c) three dilatational waves that are coupled due to the presence of thermal properties of the material. The set of dilatational waves consists of a quasi-elastic longitudinal wave and two quasi-thermal standing waves. The two transverse elastic waves are not subjected to the dispersion, while the other two transverse thermal standing waves and the dilatational waves present the dispersive character. Explicit expressions for all these seven waves are presented. The Rayleigh surface wave propagation problem is addressed and the secular equation is obtained in an explicit form. Numerical computations are performed for a specific model, and the results obtained are depicted graphically.     

Rayleigh waves with impedance boundary conditions in anisotropic solids

The paper is concerned with the propagation of Rayleigh waves in an elastic half-space with impedance boundary conditions. The half-space is assumed to be orthotropic and monoclinic with the symmetry plane _x3=0$x_3=0$. The main aim of the paper is to derive explicit secular equations of the wave. For the orthotropic case, the secular equation is obtained by employing the traditional approach. It is an irrational equation. From this equation, a new version of the secular equation for isotropic materials is derived. For the monoclinic case, the method of polarization vector is used for deriving the secular equation and it is an algebraic equation of eighth-order. When the impedance parameters vanish, this equation coincides with the secular equation of Rayleigh waves with traction-free boundary conditions.     

Modeling acoustic waves in locally enhanced meshes with a staggered-grid finite difference approach

Finite difference is a well-suited technique for modeling acoustic wave propagation in heterogeneous media as well as for imaging and inversion. Typically, the method aims at solving a set of partial differential equations for the unknown pressure field by using a regularly spaced grid. Although finite differences can be fast and cheap to implement, the accuracy of the solution is always restricted by the computational resources. This is a fundamental key point to treat when dealing with large-scale problems. In this work, we present and test a method that uses a non-uniform distribution of grid points to improve on accuracy or to reduce the required computational resources. The applied grid is generated through a coordinate transformation. Differential geometry and generalized coordinates are used to handle and analyze the effect of using a non-uniform grid. Results obtained with the presented method show that the applied transformation as well as the number of points-per-wavelength influences the stability and dispersion in the solution. We exploit this observation to locally improve the accuracy of our simulations. The work presented in this paper allows us to conclude that differential geometry for finite differences can be used to reduce dispersion and hence improve the accuracy when modeling acoustic wave propagation in heterogeneous media. In addition, it can be used to avoid oversampling through the optimization of the number of grid nodes required to have an accurate solution or just honor to the boundaries.     

Applications of the dual integral formulation in conjunction with fast multipole method to the oblique incident wave problem

In this paper, the dual integral formulation is derived for the modified Helmholtz equation in the propagation of oblique incident wave passing a thin barrier (zero thickness) by employing the concept of fast multipole method (FMM) to accelerate the construction of an influence matrix. By adopting the addition theorem, the four kernels in the dual formulation are expanded into degenerate kernels that separate the field point and the source point. The source point matrices decomposed in the four influence matrices are similar to each other or only to some combinations. There are many zeros or the same influence coefficients in the field point matrices decomposed in the four influence matrices, which can avoid calculating the same terms repeatedly. The separable technique reduces the number of floating‐point operations from O((N)²) to O(N log^a(N)), where N is the number of elements and a is a small constant independent of N. Finally, the FMM is shown to reduce the CPU time and memory requirement, thus enabling us to apply boundary element method (BEM) to solve water scattering problems efficiently. Two‐moment FMM formulation was found to be sufficient for convergence in the singular equation. The results are compared well with those of conventional BEM and analytical solutions and show the accuracy and efficiency of the FMM. Copyright © 2008 John Wiley & Sons, Ltd.     

理想导体中具有两个松弛时间的电磁热弹性波

研究处于均布磁场中的理想导体的二维电磁热弹性耦合问题，引入势函数使控制方程转化为3个偏微分方程．运用Laplace变换和Fourier变换得到该问题在变换域内的精确表达式，再通过级数展开和Laplace逆变换法求得在时间较短时的逆变换，得到时间-空间域内问题的解．运用此方法研究了表面受到热冲击的半无限空间问题．给出了电磁热弹性波、膨胀波和横向波传播的速度，并通过数值计算，给出了各个场量的分布图．所得结论与已有的结论一致．     

SIMULATION OF SEDIMENTATION OF TWO CIRCULAR PARTICLES WITH COLLISION CONSIDERED IN VERTICAL CHANNEL

D2Q9 model of lattice Boltzmann equation method was used to simulate the sedimentation of two circular particles in a bounded two dimension channel. The characteristics of the sedimentation shows some periodicity for the Reynolds number Re chosen, 0.1-20. The larger the Reynolds number, the stronger the interaction between the two particles and the larger the transversal displacements. For large Re, the two particles leading alternately; for small Re, the initially leading particle will keep its leading position and for moderate Re, the initially upper particle will get leading position and keep it. The influence of the initially relative position of the two particles on sedimentation is small. The width of the channel won＇t change the characteristics of the sedimentation as a whole, but will change the period of the sedimentation. The wider the channel, the longer the period will be.     

Deformation due to time harmonic sources in micropolar thermoelastic medium possessing cubic symmetry with two relaxation times

The response of a micropolar thermoelastic medium possessing cubic symmetry with two relaxation times due to time harmonic sources is investigated. Fourier transform is employed and the transform is inverted by using a numerical inversion technique. The components of displacement, stress, microrotation and temperature distribution in the physical domain are obtained numerically. The results of normal displacement, normal force stress, tangential couple stress and temperature distribution are compared for micropolar cubic crystal and micropolar isotropic solid. The numerical results are illustrated graphically for a particular material. Some special cases are also deduced.     

Features of how surface waves form in a cylinder made of a hard material and filled with a fluid

The properties of harmonic surface waves in an elastic cylinder made of a rigid material and filled with a fluid are studied. The problem is solved using the dynamic equations of elasticity and the equations of motion of a perfect compressible fluid. It is shown that two surface (Stoneley and Rayleigh) waves exist in this waveguide system. The first normal wave generates a Stoneley wave on the inner surface of the cylinder. If the material is rigid, no normal wave exists to transform into a Rayleigh wave. The Rayleigh wave on the outer surface forms on certain sections of different dispersion curves. The kinematic and energy characteristics of surface waves are analyzed. As the wave number increases, the phase velocities of all normal waves, except the first one, tend to the sonic velocity in the fluid from above __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 9, pp. 48–62, September 2007.