Time Domain Acoustic Propagation Analysis Using 2-D Pseudo-spectral Modeling for Ocean Environment

해양환경에서 2차원 유사 스펙트럴 모델링을 이용한 시간 영역 음 전달 해석

  • 김기산 (서울대학교 조선해양공학과) ;
  • 이근화 (서울대학교 조선해양공학과) ;
  • 성우제 (서울대학교 조선해양공학과)
  • Published : 2004.11.01

Abstract

A computer code that is based on the Pseudo-spectral finite difference algorithm using staggered grid is developed for the wave propagation modeling in the time domain. The advantage of a finite difference approximation is that any geometrically complicated media can be modeled. Staggered grids are advantageous as it provides much more accuracy than using a regular grid. Pseudo-spectral methods are those that evaluate spatial derivatives by multiplying a wavenumber by the Fourier transform of a pressure wave-field and performing the inverse Fourier transform. This method is very stable and reduces memory and the number of computations. The synthetic results by this algorithm agree with the analytic solution in the infinite and half space. The time domain modeling was implemented in various models. such as half-space. Pekeris waveguide, and range dependent environment. The snapshots showing the total wave-field reveals the Propagation characteristic or the acoustic waves through the complex ocean environment.

본 논문에서는 시간영역에서 음파 전달 모델링을 위해 엇갈림 격자에서 유사 스펙트럴 유한 차분 알고리듬을 기반으로 한 전산조직을 개발하였다. 유한 차분 근사는 기하학적으로 복잡한 매질에서 모델링을 가능하게 하고, 엇갈림 격자는 정규 격자에 비해 훨씬 정확한 해를 제공한다. 유사 스펙트럴 방법은 파수 영역에서 파수에 음압의 푸리에 변환을 곱한 후 이를 역푸리에 변환하므로서 공간 도함수를 구하는 방법이다. 이 방법은 매우 안정적이며, 메모리와 계산시간을 감소시키는 장점을 지니고 있다. 무한 및 반무한 영역에서 이 알고리듬에 의한 결과가 해석해와 잘 일치함을 확인하였고, 무한영역과 Pekeris도파관, 거리종속 해양환경에서 시간영역 모델링을 수행하여 스냅사진을 얻어내었으며 이 스냅사진을 통해 복잡한 해양환경에서 신호의 전파 양상을 파악할 수 있었다.

Keywords

References

  1. F. B. Jensen. W. A. Kuperman, M. B. Porter. and H. Schmidt, Computational Ocean Acoustics, (AlP Press, New York, 1994). Chap.8, 477-514
  2. J. E. Murphy, 'Finite-difference treatment of a time-domain parabolic equation: Theory,' .J. Acousl. Soc. Am. 77, 1958-1960. 1985 https://doi.org/10.1121/1.391844
  3. B. E. McDonald and W. A. Kuperman, 'Time domain formulation for pulse propagation including nonlinear behavior at a caustic,' J. Acousl. Soc. Am. 81. 1406-1417, 1987 https://doi.org/10.1121/1.394546
  4. M. D. Collins, 'The time-domain solution of the wide angle parabolic equation including the effects of sediment dispersion,' J. Acousl. Soc. Am. 84. 2114-2125, 1988 https://doi.org/10.1121/1.397057
  5. M. B. Porter, 'The time-marcred fast-field program(FFP) for modeling acoustic pulse propagation,' J. Acoust. Soc. Am. 87, 2013-2083, 1990 https://doi.org/10.1121/1.399329
  6. D. Kosloff and E. Baysal, 'Forward modeling by a Fourier method,' Geophysics 47, 1402-1412, 1980 https://doi.org/10.1190/1.1441288
  7. M. Reshef, D. Kosi off, M. Edwards, and C. Hsiung, 'Threedimensional acoustic modeling by the Fourier method,' Geophysics 53, 1184-1193, 1988 https://doi.org/10.1190/1.1442558
  8. T. $\ddot{O}zdenvar$ and G. A. McMechan, 'Causes and reduction of numerical artefacts in pseudo--spectral wavefield xtrapolation.' Geophys. J. Int. 126, 819--828, 199 https://doi.org/10.1111/j.1365-246X.1996.tb04705.x
  9. T. $\ddot{O}zdenvar$ and G. A. McMechan, 'Algorithms for staggered grid computations for poroelastic, elastic and acoustic, and scalar wave equations,' Geophysical Prospecting 45, 403-420, 1997 https://doi.org/10.1046/j.1365-2478.1997.390275.x
  10. Aki, K. and P. G. Richards, Quantitative seismology; the second edition, (Univ. Science Books, Sansalito, CA, 2002), Chap.6, 218-235
  11. P. Gerstoft, CABRILLO 1.0 : Acoustic, elastic and poroelastic finite difference modeling. User's manual. U. C. S. D. Marine Physical Laboratory, La Jolla. USA, 2002
  12. C. Cerjan, D. Koloff, R. Kosloff, and M. Reshef, 'A nonreflecting boundary condition for discrete acoustic and elastic wave equations,' Geophysics 50, 705-708, 1985 https://doi.org/10.1190/1.1441945