skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

Abstract

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratifiedmore » wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.« less

Authors:
; ;  [1]
  1. Sorbonne Universités, UPMC Univ Paris 06, CNRS, UMR 7190, Institut Jean Le Rond d’Alembert, F-75005, Paris (France)
Publication Date:
OSTI Identifier:
22492650
Resource Type:
Journal Article
Journal Name:
AIP Conference Proceedings
Additional Journal Information:
Journal Volume: 1685; Journal Issue: 1; Conference: 20. international symposium on nonlinear acoustics, Ecully (France), 29 Jun - 3 Jul 2015, 2. international sonic boom forum, Ecully (France), 29 Jun - 3 Jul 2015; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 0094-243X
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ABSORPTION; AMPLITUDES; COMPUTERIZED SIMULATION; DIFFRACTION; NOISE; NONLINEAR PROBLEMS; SHOCK WAVES; SOUND WAVES; THREE-DIMENSIONAL CALCULATIONS; TURBINE BLADES; TURBOMACHINERY; TURBULENT FLOW; VELOCITY; WAVE EQUATIONS; WIND

Citation Formats

Luquet, David, Marchiano, Régis, and Coulouvrat, François. 3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium. United States: N. p., 2015. Web. doi:10.1063/1.4934447.
Luquet, David, Marchiano, Régis, & Coulouvrat, François. 3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium. United States. https://doi.org/10.1063/1.4934447
Luquet, David, Marchiano, Régis, and Coulouvrat, François. 2015. "3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium". United States. https://doi.org/10.1063/1.4934447.
@article{osti_22492650,
title = {3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium},
author = {Luquet, David and Marchiano, Régis and Coulouvrat, François},
abstractNote = {Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.},
doi = {10.1063/1.4934447},
url = {https://www.osti.gov/biblio/22492650}, journal = {AIP Conference Proceedings},
issn = {0094-243X},
number = 1,
volume = 1685,
place = {United States},
year = {Wed Oct 28 00:00:00 EDT 2015},
month = {Wed Oct 28 00:00:00 EDT 2015}
}