The enhanced volume source boundary point method for the calculation of acoustic radiation problem

The Volume Source Boundary Point Method(VSBPM) is greatly improved so that it will speed up the VSBPM‘s solution of the acoustic radiation problem caused by the vibrating body.The fundamental solution provided by Helmholtz equation is enforced in a weighted residual sense over a tetrahedron located on the normal line of the boundary node to replace the coefficient matrices of the system equation.Through the enhanced volume source boundary point analysis of various examples and the sound field of a vibrating rectangular box in a semi-anechoic chamber,it has revealed that the calculating speed of the EVSBPM is more than 10 times faster than that of the VSBPM while it workss on the aspects of its calculating precision and stability,adaptation to geometric shape of vibrating body as well as its ability to overcome the non-uniqueness problem.     

An application of radiation boundary condition to scattering problem of acoustic waves in fluids

A numerical method of solving acoustic wave scattering pnblemin fluids is described.Radiation boundary condition(RBC)obtained byfactorization method of Helmholtz equation is applied to transforming theexterior boundary value problem in unbounded region into one in a finiteregion.Combined with RBC and scatterer surface boundary condition,Helmholtz equation is solved numerically by the finite difference method.Computational results for sphere and prolate spheroidal scatterers are inexcellent agreement with eigenfunction solutions and much better than theresults of OSRC method.     

On the application of wavelet transform to the solution of integral equations for acoustic radiation and scattering

1 IlltroductionIt is well Anown that Boundary Element Method (BEM), based on boundary integralequation, leads the reduction of the dimensionallty of the problem by one because the prob1emis formulated in terms of the flelds on botmdary ouly BEM, hOWver, generates algebraicequations with full matrices, whose solutions are moe eapensive than that of the bandedmatrices of FEM[1'2]. On the other hand, in boundary integral forInulations, an integral operatorhas the global behavior, that can b…     

Calculation of acoustic scattering of a nonrigid surface using physical acoustic method

In this paper the physical acoustic method or the Kirchhoff approxima-tion is extended to treat the scattering of a nonrigid surface in order to estimatethe target strength of targets with absorbing coatings.By using the locally planewave approximation,the relationship between the sound pressure and its normalderivative on the surface can be represented by the plane wave reflectioncoefficient and the acoustic impedance of the surface.The resulting modifiedKirchhoff approximation involves the plane wave reflection coefficient.For aimpedance sphere,a comparison between the physical acoustic method and theexact solution shows that the physical acoustic method still is a good approxima-tion at higher κα values.     

A coordinate–spectral method for rough surface scattering

Abstract

In this paper, we deal with the time-harmonic scattering by one-dimensional rough surfaces separating two homogeneous and isotropic media. The method is based on a rigorous integral formalism. The unknown of the integral equation is projected onto a Fourier basis while the equation itself is sampled as in a classical method of moments. The accuracy is tested against both other methods and experimental results. One of the main interests in choosing a Fourier basis lies in the ability to solve rigorously the scattering of a p polarized incident beam by a shallow metallic rough surface. The role of the surface waves is accurately taken into account and phenomena such as enhanced backscattering are well described. With this method, one can consider that the gap between the domain of validity of perturbation theories and the domain of practical use of rigorous methods is filled.     

Axial acoustic radiation force on an elastic spherical shell near an impedance boundary for zero-order quasi-Bessel–Gauss beam

Shell structures have increasingly widespread applications in biomedical ultrasound fields such as contrast agents and drug delivery,which requires the precise prediction of the acoustic radiation force under various circumstances to improve the system efficiency.The acoustic radiation force exerted by a zero-order quasi-Bessel-Gauss beam on an elastic spherical shell near an impedance boundary is theoretically and numerically studied in this study.By means of the finite series method and the image theory,a zero-order quasi-Bessel-Gauss beam is expanded in terms of spherical harmonic functions,and the exact solution of the acoustic radiation force is derived based on the acoustic scattering theory.The acoustic radiation force function,which represents the radiation force per unit energy density and per unit cross-sectional surface,is especially investigated.Some simulated results for a polymethyl methacrylate shell and an aluminum shell are provided to illustrate the behavior of acoustic radiation force in this case.The simulated results show the oscillatory property and the negative radiation force caused by the impedance boundary.An appropriate relative thickness of the shell can generate sharp peaks for a polymethyl methacrylate shell.Strong radiation force can be obtained at small half-cone angles and the beam waist only affects the results at high frequencies.Considering that the quasi-Bessel-Gauss beam possesses both the energy focusing property and the non-diffracting advantage,this study is expected to be useful in the development of acoustic tweezers,contrast agent micro-shells,and drug delivery applications.     

An efficient high-order boundary element method for nonlinear wave–wave and wave-body interactions

A highly efficient high-order boundary element method is developed for the numerical simulation of nonlinear wave–wave and wave-body interactions in the context of potential flow. The method is based on the framework of the quadratic boundary element method (QBEM) for the boundary integral equation and uses the pre-corrected fast Fourier transform (PFFT) algorithm to accelerate the evaluation of far-field influences of source and/or normal dipole distributions on boundary elements. The resulting PFFT–QBEM reduces the computational effort of solving the associated boundary-value problem from O(N^2～3) (with the traditional QBEM) to O(N ln N) where N represents the total number of boundary unknowns. Significantly, it allows for reliable computations of nonlinear hydrodynamics useful in ship design and marine applications, which are forbidden with the traditional methods on the presently available computing platforms. The formulation and numerical issues in the development and implementation of the PFFT–QBEM are described in detail. The characteristics of accuracy and efficiency of the PFFT–QBEM for various boundary-value problems are studied and compared to those of the existing accelerated (lower- and higher-order) boundary element methods. To illustrate the usefulness of the PFFT–QBEM, it is applied to solve the initial boundary-value problem in the generation of three-dimensional nonlinear waves by a moving ship hull. The predicted wave profile and resistance on the ship are compared to available experimental measurements with satisfactory agreements.     

Research on modal analysis of structural acoustic radiation using structural vibration modes and acoustic radiation modes

Modal analysis of structural acoustic radiation from a vibrating structure is discussed using structural vibration modes and acoustic radiation modes based on the quadratic form of acoustic power. The finite element method is employed for discretisizing the structure. The boundary element method and Rayleigh integral are used for modeling the acoustic fluid. It is shown that the power radiated by a single vibration mode is to increase the radiated power and the effect of modal interaction can lead to an increase or a decrease or no change in the radiated power, moreover, control of vibration modes is a good way to reduce both vibration and radiated sound as long as the influence of interaction of vibration modes on sound radiation is insignificant. Stiffeners may change mode shapes of a plate and thus change radiation efficiency of the plate‘s modes. The CHIEF method is adopted to obtain an acoustic radiation mode formulation without the nonuniqueness difficulty at critical frequencies for three-dimensional structures by using Moore-Penrose inverse. A pulsating cube is involved to verify the formulation. Good agreement is obtained between the numerical and analytical solutions. The shapes and radiation efficiencies of acoustic radiation modes of the cube are discussed. The structural acoustic control using structural vibration modes and acoustic radiation modes are compared and studied.     

An implicit immersed boundary method for three-dimensional fluid–membrane interactions

We present an implicit immersed boundary method for the incompressible Navier–Stokes equations capable of handling three-dimensional membrane–fluid flow interactions. The goal of our approach is to greatly improve the time step by using the Jacobian-free Newton–Krylov method (JFNK) to advance the location of the elastic membrane implicitly. The most attractive feature of this Jacobian-free approach is Newton-like nonlinear convergence without the cost of forming and storing the true Jacobian. The Generalized Minimal Residual method (GMRES), which is a widely used Krylov-subspace iterative method, is used to update the search direction required for each Newton iteration. Each GMRES iteration only requires the action of the Jacobian in the form of matrix–vector products and therefore avoids the need of forming and storing the Jacobian matrix explicitly. Once the location of the boundary is obtained, the elastic forces acting at the discrete nodes of the membrane are computed using a finite element model. We then use the immersed boundary method to calculate the hydrodynamic effects and fluid–structure interaction effects such as membrane deformation. The present scheme has been validated by several examples including an oscillatory membrane initially placed in a still fluid, capsule membranes in shear flows and large deformation of red blood cells subjected to stretching force.     

The immersed boundary method for advection–electrodiffusion with implicit timestepping and local mesh refinement

We describe an immersed boundary method for problems of fluid–solute-structure interaction. The numerical scheme employs linearly implicit timestepping, allowing for the stable use of timesteps that are substantially larger than those permitted by an explicit method, and local mesh refinement, making it feasible to resolve the steep gradients associated with the space charge layers as well as the chemical potential, which is used in our formulation to control the permeability of the membrane to the (possibly charged) solute. Low Reynolds number fluid dynamics are described by the time-dependent incompressible Stokes equations, which are solved by a cell-centered approximate projection method. The dynamics of the chemical species are governed by the advection–electrodiffusion equations, and our semi-implicit treatment of these equations results in a linear system which we solve by GMRES preconditioned via a fast adaptive composite-grid (FAC) solver. Numerical examples demonstrate the capabilities of this methodology, as well as its convergence properties.     

Time domain finite volume method for three-dimensional structural–acoustic coupling analysis

This work extends the application of finite volume method (FVM) to structural–acoustic problems. A three-dimensional time domain FVM (TDFVM) is proposed to predict the transient response and natural characteristics of structural–acoustic coupling systems. Acoustic wave equation in heterogeneous medium and structural dynamic equation are solved in fluid and solid sub-domains respectively. The structural–acoustic coupling is implemented according to normal components of particle acceleration continuity condition and normal traction equilibrium condition at the interface. The computational domain is discretized with four-node tetrahedral grid which is generated easily and has strong adaptability to complicated geometries. Numerical experiments are carried out to examine the accuracy of the method in both time domain and frequency domain. The results show good agreement with analytical solutions and numerical results. For structural–acoustic problem, TDFVM has the capability to consider the heterogeneity of both fluid and solid.     

The high frequency acoustic radiation from the boundary layer of an axisymmetric body

1 IntroductionWhen an edsynunetric body moves under water, a boundary layer is formed at the surfaceof the body. There are three regions in the boundary layer: lallilnar, transition from the laminarto turbulent, and turbulent region. The boundary layyer near the stagnation region remainslaminar and subsequently it goes through transition to turbulence. The acoustic radiation fromthe transition and turbulellt boundary 18y6r is the main components of self-noise of marines,and so its mechanical s…     

An improved penalty immersed boundary method for fluid–flexible body interaction

An improved penalty immersed boundary (pIB) method has been proposed for simulation of fluid–flexible body interaction problems. In the proposed method, the fluid motion is defined on the Eulerian domain, while the solid motion is described by the Lagrangian variables. To account for the interaction, the flexible body is assumed to be composed of two parts: massive material points and massless material points, which are assumed to be linked closely by a stiff spring with damping. The massive material points are subjected to the elastic force of solid deformation but do not interact with the fluid directly, while the massless material points interact with the fluid by moving with the local fluid velocity. The flow solver and the solid solver are coupled in this framework and are developed separately by different methods. The fractional step method is adopted to solve the incompressible fluid motion on a staggered Cartesian grid, while the finite element method is developed to simulate the solid motion using an unstructured triangular mesh. The interaction force is just the restoring force of the stiff spring with damping, and is spread from the Lagrangian coordinates to the Eulerian grids by a smoothed approximation of the Dirac delta function. In the numerical simulations, we first validate the solid solver by using a vibrating circular ring in vacuum, and a second-order spatial accuracy is observed. Then both two- and three-dimensional simulations of fluid–flexible body interaction are carried out, including a circular disk in a linear shear flow, an elastic circular disk moving through a constricted channel, a spherical capsule in a linear shear flow, and a windsock in a uniform flow. The spatial accuracy is shown to be between first-order and second-order for both the fluid velocities and the solid positions. Comparisons between the numerical results and the theoretical solutions are also presented.     

Absorbing boundary conditions for the Euler and Navier–Stokes equations with the spectral difference method

Two absorbing boundary conditions, the absorbing sponge zone and the perfectly matched layer, are developed and implemented for the spectral difference method discretizing the Euler and Navier–Stokes equations on unstructured grids. The performance of both boundary conditions is evaluated and compared with the characteristic boundary condition for a variety of benchmark problems including vortex and acoustic wave propagations. The applications of the perfectly matched layer technique in the numerical simulations of unsteady problems with complex geometries are also presented to demonstrate its capability.     

A matched interface and boundary method for solving multi-flow Navier–Stokes equations with applications to geodynamics

We have developed a second-order numerical method, based on the matched interface and boundary (MIB) approach, to solve the Navier–Stokes equations with discontinuous viscosity and density on non-staggered Cartesian grids. We have derived for the first time the interface conditions for the intermediate velocity field and the pressure potential function that are introduced in the projection method. Differentiation of the velocity components on stencils across the interface is aided by the coupled fictitious velocity values, whose representations are solved by using the coupled velocity interface conditions. These fictitious values and the non-staggered grid allow a convenient and accurate approximation of the pressure and potential jump conditions. A compact finite difference method was adopted to explicitly compute the pressure derivatives at regular nodes to avoid the pressure–velocity decoupling. Numerical experiments verified the desired accuracy of the numerical method. Applications to geophysical problems demonstrated that the sharp pressure jumps on the clast-Newtonian matrix are accurately captured for various shear conditions, moderate viscosity contrasts and a wide range of density contrasts. We showed that large transfer errors will be introduced to the jumps of the pressure and the potential function in case of a large absolute difference of the viscosity across the interface; these errors will cause simulations to become unstable.     

A temperature-related boundary Cauchy–Born method for multi-scale modeling of silicon nano-structures

The surface, edge and corner effects have significant influences in the electrical and optical properties of silicon nano-structures. In this paper, a novel hierarchical temperature-related multi-scale model is presented based on the boundary Cauchy–Born method to investigate not only the surface but also the edge and corner effects in thermal properties of diamond-like structures such as silicon nano-structures at finite temperature. A combined finite element method and molecular dynamics are respectively employed in macro- and micro-scale levels. The temperature-related Cauchy–Born rule is applied using the Helmholtz free energy, as the energy density of equivalent continua relating to the Tersoff inter-atomic potential. The model employs radial quadratures at the surface, edge and corner elements as an indicator of material behavior. The capability of computational algorithm is illustrated by numerical simulation of a nano-scale cube at finite temperature and the results are compared with the atomistic model.     

A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic–acoustic media. A velocity–strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic–acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic–acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.     

A new time–space domain high-order finite-difference method for the acoustic wave equation

A new unified methodology was proposed in Finkelstein and Kastner (2007) [39] to derive spatial finite-difference (FD) coefficients in the joint time–space domain to reduce numerical dispersion. The key idea of this method is that the dispersion relation is completely satisfied at several designated frequencies. We develop this new time–space domain FD method further for 1D, 2D and 3D acoustic wave modeling using a plane wave theory and the Taylor series expansion. New spatial FD coefficients are frequency independent though they lead to a frequency dependent numerical solution. We prove that the modeling accuracy is 2nd-order when the conventional (2M)$(2M)$th-order space domain FD and the 2nd-order time domain FD stencils are directly used to solve the acoustic wave equation. However, under the same discretization, the new 1D method can reach (2M)$(2M)$th-order accuracy and is always stable. The 2D method can reach (2M)$(2M)$th-order accuracy along eight directions and has better stability. Similarly, the 3D method can reach (2M)$(2M)$th-order accuracy along 48 directions and also has better stability than the conventional FD method. The advantages of the new method are also demonstrated by the results of dispersion analysis and numerical modeling of acoustic wave equation for homogeneous and inhomogeneous acoustic models. In addition, we study the influence of the FD stencil length on numerical modeling for 1D inhomogeneous media, and derive an optimal FD stencil length required to balance the accuracy and efficiency of modeling. A new time–space domain high-order staggered-grid FD method for the 1D acoustic wave equation with variable densities is also developed, which has similar advantages demonstrated by dispersion analysis, stability analysis and modeling experiments. The methodology presented in this paper can be easily extended to solve similar partial difference equations arising in other fields of science and engineering.     

Non—equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method

In this paper, we propose a new approach to implementing boundary conditions in the lattice Boltzmann method (LBM). The basic idea is to decompose the distribution function at the boundary node into its equilibrium and non-equilibrium parts, and then to approximate the non-equilibrium part with a first-order extrapolation of the non-equilibrium part of the distribution at the neighbouring fluid node. Schemes for velocity and pressure boundary conditions are constructed based on this method. The resulting schemes are of second-order accuracy. Numerical tests show that the numerical solutions of the LBM together with the present boundary schemes are in excellent agreement with the analytical solutions. Second-order convergence is also verified from the results. It is also found that the numerical stability of the present schemes is much better than that of the original extrapolation schemes proposed by Chen et al. (1996 Phys. Fluids 8 2527).