MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/348314 348314 Research Article SPH Simulation of Acoustic Waves: Effects of Frequency, Sound Pressure, and Particle Spacing Zhang Y. O. 1, 2 Zhang T. 1, 3 http://orcid.org/0000-0003-0312-0326 Ouyang H. 4 Li T. Y. 1, 3 Liew Kim M. 1 School of Naval Architecture and Ocean Engineering Huazhong University of Science and Technology Wuhan 430074 China hust.edu.cn 2 Department of Mechanical and Aerospace Engineering University of California San Diego, La Jolla, CA 92093 USA ucsd.edu 3 Hubei Key Laboratory of Naval Architecture and Ocean Engineering Hydrodynamics Huazhong University of Science and Technology Wuhan 430074 China hust.edu.cn 4 School of Engineering University of Liverpool The Quadrangle, L69 3GH Liverpool UK liv.ac.uk 2015 12102015 2015 30 08 2014 15 11 2014 15 11 2014 12102015 2015 Copyright © 2015 Y. O. Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Acoustic problems consisting of multiphase systems or with deformable boundaries are difficult to describe using mesh-based methods, while the meshfree, Lagrangian smoothed particle hydrodynamics (SPH) method can handle such complicated problems. In this paper, after solving linearized acoustic equations with the standard SPH theory, the feasibility of the SPH method in simulating sound propagation in the time domain is validated. The effects of sound frequency, maximum sound pressure amplitude, and particle spacing on numerical error and time cost are then subsequently discussed based on the sound propagation simulation. The discussion based on a limited range of frequency and sound pressure demonstrates that the rising of sound frequency increases simulation error, and the increase is nonlinear, whereas the rising sound pressure has limited effects on the error. In addition, decreasing the particle spacing reduces the numerical error, while simultaneously increasing the CPU time. The trend of both changes is close to linear on a logarithmic scale.

1. Introduction

Some classic numerical methods such as the finite element method (FEM) [1, 2], the boundary element method (BEM) , and other modified or coupled methods  are widely used for acoustic simulations. However, these mesh-based methods are not ideal for solving acoustic problems consisting of a variety of media or with deformable boundaries.

Meshfree methods can handle such complicated problems. The method of fundamental solutions (MFS) , the multiple-scale reproducing kernel particle method (RKPM) , the element-free Galerkin method (EFGM) , and other meshfree methods [10, 11] have been applied to these acoustic problems.

The smoothed particle hydrodynamics (SPH) method, as a meshfree, Lagrangian method, was first independently pioneered by Lucy  and Gingold and Monaghan  to solve astrophysical problems in 1977. In addition, the SPH method has been used in many different fields . It not only has most advantages of a meshfree method, but also is suitable for solving problems with material separation or large ranges of density as illustrated in recent reviews by Li and Liu , Springel , M. B. Liu and G. R. Liu , and Monaghan  due to its Lagrangian property. Introducing the SPH method to acoustic computation also brings its advantages to some fields like bubble acoustic, combustion noise, sound propagation in multiphase flows, and so on.

With the advance of the SPH method in acoustic simulation, some research literatures [21, 22] discussed solving fluid dynamic equations to simulate sound waves. In addition, we published a conference paper  that used the SPH method to solve linearized acoustic equations for modeling sound propagation and interference. Numerical results showed that the SPH method was capable of accurately modeling sound propagation, but the effects of frequency and sound pressure on the SPH simulation need further discussion. Therefore, the present paper focuses on discussing the effect of frequency and sound pressure of the acoustic waves on the numerical error caused by the SPH simulation.

The present paper is organized as follows. In Section 2, the standard SPH theory is used to solve the linearized acoustic equations. In Section 3, a one-dimensional sound propagation model is built. In Section 4, a numerical experiment is given based on standard SPH algorithms, and the effect of frequency, sound pressure, and particle spacing on the simulation is analyzed with considering the changes of particle spacing and Courant number. Section 5 summarizes the results of this work.

2. SPH Formulations of Sound Waves 2.1. Basic Formulations of SPH

As a meshfree, Lagrangian particle method, the SPH method is an important method widely used in recent years. Formulations in the SPH theory are represented in a particle approximation form. The properties of each particle are computed using an interpolation process over its neighboring particles . In this way, the integral of a field function f(r) can be represented as(1)fr=ΩfrWr-r,hdr,where · is the kernel approximation operator, f is a function of the vector r,Ω is the volume of the integral, W is the smoothing kernel, and h is the smoothing length.

The particle approximation for the function f(r) at particle i within the support domain can be written as(2)fri=j=1Nmjρjfrj·Wij,where ri and rj are the position of particles i and j, N is the number of particles in the computational domain, mj is the mass of particle j, Wij=W(rij,h), and rij is the distance between particle i and particle j.

Similarly, the gradient of function f(r) at particle i is obtained as(3)·fri=j=1Nmjρjfrj·iWij,where iWij=ri-rj/rijWij/rij.

2.2. SPH Formulations of Sound Waves

In acoustic simulation, the governing equations for constructing SPH formulations are the laws of continuity, momentum, and state. The simplest and most common acoustical problem occurs under some assumptions. On one hand, the medium is lossless and at rest, so an energy equation is unnecessary; on the other hand, a small departure from quiet conditions occurs as follows:(4)ρ=ρ0+δρ,δρρ0(5)P=p0+p,pρ0c02(6)u=0+u,uc0,where ρ is the fluid density, ρ0 is the quiescent density which does not vary in time and space, δρ is the change of density, P is the instantaneous pressure at time t of the fluid, p0 is the quiescent pressure, p is the sound pressure, c0 is the speed of sound, and u is the flow velocity. δρ, p, and u are taken to be small quantities of first order.

By discarding second-order terms in the acoustic equations, the linearized continuity, momentum, and state equations (for ideal air) governing sound waves are obtained as(7)δρt=-ρ  ·u(8)ut=-1ρp(9)p=c02δρ.

Applying the SPH particle approximation (see (3)) to (7), the particle approximation equation of the continuity of acoustic waves is written as(10)δρit=ρ0+δρij=1Nmjρ0+δρjuijiWij,where ρ0 is the quiescent density which does not vary in time and space, p0 is the quiescent pressure, and uij=ui-uj.

The momentum equation in SPH method is obtained as(11)uit=j=1Nmjpiρ0+δρi2+pjρ0+δρj2iWij.

Particle approximation of the equation of state is(12)pi=c02δρi.

The cubic spline function given by Monaghan and Lattanzio  is used as the smoothing kernel in this paper, which is written as(13)Wr,h=αD1-32q2+34q30q1142-q31q20q2,where αD is 1/h in one-dimensional space, q=r/h, r is the distance between two particles, and h is the smoothing length which defines the influence area of smoothing function W.

The second-order leap-frog integration  is used in the present paper. All-pair search approach , as a direct and simple algorithm, is used to realize the neighbor particles searching in the acoustic wave simulation.

3. Sound Propagation Model

In order to evaluate the effect of sound pressure and frequency, a one-dimensional sound wave which propagates in a pipe with uniform cross section is used. The acoustic model is shown in Figure 1. The sound pressure at t=0 is plotted with a solid line, while the sound pressure at t=0.1 s is plotted with a dashed line.

One-dimensional sound propagation model.

Sound pressure of the source of sound wave in Figure 1 is(14)pt,x<0=pAmpsinwt-kx,where t represents time (propagation starts when t=0), x is the geometric position, pAmp is the maximum sound pressure of the acoustic wave, w is the circular frequency of the acoustic wave, k=w/c0 is the wave number, and the sound speed c0 is 340 m/s.

The sound propagates from x<0 to x0, with the SPH computational region going from −10 m to 70 m. In the computation, the region of the sound source is from −10 m to 0 m, which is determined by the theoretical solutions, while the region of sound propagation is from 0 m to 70 m, which is obtained by the SPH method. At the end of the sound propagation domain, the particles are set free in the computation. Since the SPH computational region is far longer than the sound propagation distance, any effects caused by ignoring this boundary will not propagate to the sound propagation domain and therefore can be neglected. The simulation results at the time t=0.1 s are used to compare with the theoretical solutions. A similar numerical model is used in our recent conference paper .

In the evaluation of numerical results, the nondimensional simulation error of the sound pressure is defined by(15)εpre=1N·pAmpj=1Npj-pxj2,where pj is the simulation pressure of the particle j and p(xj) is the theoretical pressure at the position of particle j.

4. Results and Discussion 4.1. Numerical Experiment

A sound propagation model is built to confirm agreement between SPH simulation results and the theoretical solution. The parameters used in the numerical experiment of the sound waves are listed in Table 1. The physical time required for the modeling was 141.6 seconds.

Model parameters and computational parameters for the numerical experiment.

 Model parameters p Amp w c 0 t 50 Pa 50 rad/s 340 m/s 0.1 s Computational parameters Particle spacing Smoothing length Courant number Particle mass 0.04 m 0.13 m 0.05 0.04 kg

A comparison between the SPH results and theoretical solutions of sound pressure at t=0.1 s is shown in Figure 2. The theoretical solution is plotted with a solid line while the SPH results are plotted with red points. In order to clearly identify the numerical results, the points are plotted at intervals of 12.

Sound pressure comparison between SPH results and theoretical solutions. (pAmp=  50 Pa, w=50 rad/s, and t=0.1 s).

As shown in Figure 2, the sound pressure of the sound wave changes along with the position, and a peak and a valley appear in the final wave at t=0.1 s. Comparing with the analytical data, it can be seen that the SPH method can accurately simulate the sound pressure of an acoustic wave as the values of sound pressure can be obtained precisely, and the simulation error is 1.6 × 10−4. The trend of both the SPH results and the theoretical solutions also has a similar tendency.

4.2. Effect of Frequency

Different sound propagation models with the circular frequency of sound changing from 10 rad/s to 200 rad/s are simulated, and the effect of circular frequency is discussed. In the discussion, different computational parameters, namely, the particle spacing and the Courant number, are also considered.

Figure 3 gives the numerical error for different computational cases with different w and particle spacing. Figure 4 shows the sound pressure error for different cases with different values of circular frequency and Courant number. Each point in both figures stands for a computational case. The Courant number used in the computational cases in Figure 3 is 0.05, and the particle spacing used in all cases in Figure 4 is 0.06 m.

Sound pressure error versus circular frequency for different particle spacings.

Sound pressure error versus circular frequency for different Courant numbers.

As can be seen from Figure 3, εpre increases along with the increase of w, and the increase of εpre is nonlinear. Along with the increase of w, the increase of εpre is larger. On the other hand, we can see from the line graph that εpre is much larger than other positions when w is 200 rad/s. The numerical error grows faster when the particle spacing is larger.

As shown in Figure 4, εpre increases along with the increase of w. What is more, three lines are represent three different Courant numbers, namely, three different choices of time step when the particle spacing is the same. εpre obtained from three different Courant numbers has almost the same value when w is under 100 rad/s. After that, some difference is shown among three lines when w becomes larger, but the Courant number affects the effect of w less in general.

4.3. Effect of Sound Pressure

In the present section, different sound propagation models with the maximum sound pressure amplitude changing from 10 Pa to 200 Pa are simulated, and the effects of sound pressure are discussed. At the same time, different computational parameters, namely, the particle spacing and the Courant number, are also used in the simulation. The w of the sound in all computational cases in this section is 50 rad/s.

Figure 5 illustrates the changes of sound pressure error when the numerical model uses different pAmp and particle spacing. Figure 6 describes the changes of error with different pAmp and different Courant numbers. Each point in both figures stands for a computational case. The Courant number used in all computational cases in Figure 5 is 0.05, and the particle spacing used in all cases in Figure 6 is 0.06 m.

Sound pressure error versus maximum sound pressure amplitude for different particle spacings.

Sound pressure error versus maximum sound pressure amplitude for different Courant numbers.

It can be seen from the line graph that εpre has limited change, which is a little increase, along with the increase of pAmp. The trend of the line is approximately linear when pAmp changes from 10 Pa to 200 Pa. The value of three lines in the figure is around 2.7 × 10−4, 3.7 × 10−4, and 5.2 × 10−4 corresponding to the particle spacing changes from 0.06 m to 0.1 m. In other words, when the particle spacing increases, the error shows an upward trend.

Similarly, εpre rises a little along with the increase of the amplitude of sound wave as shown in Figure 6. Specifically, εpre stays at about 2.6 × 10−4 when the amplitude of sound changes from 10 Pa to 200 Pa. Although with different Courant number, the error of numerical results has similar values. The largest difference of εpre among three lines is 2.0 × 10−5, which happens when pAmp reaches 200 Pa.

4.4. Effect of Particle Spacing

In addition to the accuracy, the convergence and the computational efficiency are also worth concern. Particle spacing, as an important parameter in the SPH simulation, affects these indexes a lot. In this section, the effects of the particle spacing are analyzed based on the model in Section 4.1 with the particle spacing (Δx) changes from 0.03 m to 0.3 m.

Different sound pressure errors corresponding to the SPH simulation results with different Δx are shown in Figure 7. Three different types of points are used in the figure to represent Courant numbers of 0.05, 0.10, and 0.15 in the computation. In order to clearly show the trend, logarithmic scales are used for both axes in the figure.

Sound pressure error versus particle spacing for different Courant numbers.

As can be seen from Figure 7, there is an approximately linear relation between εpre and Δx with double logarithmic coordinates. As Δx increases, εpre also increases. This result indicates a convergence of the simulation, as a smaller Δx leads to a more accurate result. Note that the errors generated with the three different Courant numbers all have similar values. A linear fitted curve obtained by using all εpre results produces a slope of approximately 1.4 for the fitted line.

The required CPU cost in terms of time to perform the SPH simulation corresponding to different Δx values is given in Figure 8. Three different point types were also used to represent different Courant numbers in the computation. The figure also shows the linear fit lines for each Courant number, as well as the calculated slopes.

CPU time versus particle spacing for different Courant numbers.

The graph shows that the CPU time logarithmically decreases from hundreds of seconds to about one second as Δx increases from 0.03 m to 0.3 m. This decreasing trend is roughly linear, with an absolute slope value of approximately 2.5. Although the CPU time is different for each Courant number, the trends are similar. Thus, an appropriate particle spacing value can be decided within a certain range by comparing the results in Figures 7 and 8 for a given application.

5. Conclusions

In this work, the meshfree, Lagrangian SPH method is used to solve linearized acoustic equations. The feasibility of the SPH method in simulating acoustic waves is evaluated using a numerical experiment of one-dimensional sound propagation model. The effect of frequency, sound pressure, and particle spacing is then subsequently discussed based on the numerical error with considering the Courant number.

The results of computational cases, with the acoustic wave frequency changes from 10 rad/s to 200 rad/s, show that the rising of sound frequency increases simulation error, and the increase is nonlinear.

As the maximum sound pressure amplitude changes from 10 Pa to 200 Pa, the rising sound pressure increases the simulation error. However, the resulting change in error is limited and marginal.

Decreasing the particle spacing in the computation reduces the numerical error, while simultaneously increasing the CPU time. The trend of both changes is close to linear on a logarithmic scale. The Courant number has little effect on the error generated with different sound frequency or sound pressure, whereas the particle spacing has a more significant influence on both properties.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was supported by the Independent Innovation Foundation for National Defence of Huazhong University of Science and Technology (no. 01-18-140019). Y. O. Zhang is supported by the China Scholarship Council (201406160032).

Ihlenburg F. Finite Element Analysis of Acoustic Scattering 1998 Springer 10.1007/b98828 MR1639879 Harari I. A survey of finite element methods for time-harmonic acoustics Computer Methods in Applied Mechanics and Engineering 2006 195 13–16 1594 1607 10.1016/j.cma.2005.05.030 MR2203983 2-s2.0-30944443982 Kythe P. K. An Introduction to Boundary Element Methods 1995 CRC Press MR1396530 Warszawski A. Soares D. Jr. Mansur W. J. A FEM-BEM coupling procedure to model the propagation of interacting acoustic-acoustic/acoustic-elastic waves through axisymmetric media Computer Methods in Applied Mechanics and Engineering 2008 197 45–48 3828 3835 10.1016/j.cma.2008.03.005 2-s2.0-50049129895 Soares D. Jr. Coupled numerical methods to analyze interacting acoustic-dynamic models by multidomain decomposition techniques Mathematical Problems in Engineering 2011 2011 28 245170 10.1155/2011/245170 MR2826916 2-s2.0-81555199746 Feng J. Zheng X. Wang H. Zou Y. Liu Y. Yao Z. Low-frequency acoustic-structure analysis using coupled FEM-BEM method Mathematical Problems in Engineering 2013 2013 8 583079 MR3121582 10.1155/2013/583079 Tadeu A. Castro I. Coupling the BEM/TBEM and the MFS for the numerical simulation of wave propagation in heterogeneous fluid-solid media Mathematical Problems in Engineering 2011 2011 26 159389 10.1155/2011/159389 MR2853312 2-s2.0-84855596813 Uras R. A. Chang C. T. Chen Y. Liu W. K. Multiresolution reproducing kernel particle methods in acoustics Journal of Computational Acoustics 1997 5 1 71 94 10.1142/S0218396X9700006X 2-s2.0-0031506587 Bouillard P. Suleau S. Element-free Galerkin solutions for Helmholtz problems: formulation and numerical assessment of the pollution effect Computer Methods in Applied Mechanics and Engineering 1998 162 1–4 317 335 10.1016/S0045-7825(97)00350-2 2-s2.0-0032136822 Wang S. Li S. Huang Q. Li K. An improved collocation meshless method based on the variable shaped radial basis function for the solution of the interior acoustic problems Mathematical Problems in Engineering 2012 2012 20 632072 10.1155/2012/632072 MR2992423 Wang F. Z. Zheng K. H. Analysis of the boundary knot method for 3D Helmholtz-type equation Mathematical Problems in Engineering 2014 2014 9 853252 10.1155/2014/853252 MR3182439 2-s2.0-84897492320 Lucy L. B. A numerical approach to the testing of the fission hypothesis The Astronomical Journal 1977 82 1013 1024 10.1086/112164 Gingold R. A. Monaghan J. J. Smoothed Particle Hydrodynamics-theory and application to non-spherical stars Monthly Notices of the Royal Astronomical Society 1977 181 375 389 Caleyron F. Combescure A. Faucher V. Potapov S. SPH modeling of fluid-solid interaction for dynamic failure analysis of fluid-filled thin shells Journal of Fluids and Structures 2013 39 126 153 2-s2.0-84880169999 10.1016/j.jfluidstructs.2013.02.023 Hopkins P. F. A general class of Lagrangian smoothed particle hydrodynamics methods and implications for fluid mixing problems Monthly Notices of the Royal Astronomical Society 2013 428 4 2840 2856 10.1093/mnras/sts210 2-s2.0-84873737244 Ma G. W. Wang Q. S. Yi X. W. Wang X. J. A modified SPH method for dynamic failure simulation of heterogeneous material Mathematical Problems in Engineering 2014 2014 14 808359 10.1155/2014/808359 2-s2.0-84900993186 Li S. Liu W. K. Meshfree and particle methods and their applications Applied Mechanics Reviews 2002 55 1 1 34 10.1115/1.1431547 2-s2.0-0000969263 Springel V. Smoothed particle hydrodynamics in astrophysics Annual Review of Astronomy and Astrophysics 2010 48 391 430 10.1146/annurev-astro-081309-130914 2-s2.0-77956361551 Liu M. B. Liu G. R. Smoothed particle hydrodynamics (SPH): an overview and recent developments Archives of Computational Methods in Engineering. State of the Art Reviews 2010 17 1 25 76 10.1007/s11831-010-9040-7 MR2593940 2-s2.0-77952876665 Monaghan J. J. Smoothed particle hydrodynamics and its diverse applications Annual Review of Fluid Mechanics 2012 44 Palo Alto, Calif, USA Annual Reviews 323 346 Annual Review of Fluid Mechanics 10.1146/annurev-fluid-120710-101220 MR2882600 2-s2.0-84555202477 Wolfe C. T. Semwal S. K. Acoustic Modeling of Reverberation Using Smoothed Particle Hydrodynamics 2007 Springs, Colo, USA University of Colorado Hahn P. Negrut D. On the use of meshless methods in acoustic simulations Proceedings of the International Mechanical Engineering Congress and Exposition (IMECE '09) November 2009 185 199 10.1115/IMECE2009-11351 2-s2.0-77954287027 Zhang Y. O. Zhang T. Ouyang H. Li T. Y. SPH simulation of sound propagation and interference Proceedings of the 5th International Conference on Computational Method (ICCM '14) July 2014 Cambridge, UK Liu G. R. Liu M. B. Smoothed Particle Hydrodynamics: a Meshfree Particle Method 2003 World Scientific Monaghan J. J. Lattanzio J. C. A refined particle method for astrophysical problems Astronomy and Astrophysics 1985 149 135 143 Zbl0622.76054 Kelager M. Lagrangian Fluid Dynamics Using Smoothed Particle Hydrodynamics 2006 Copenhagen, Denmark University of Copenhagen