Public Reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comment regarding this burden estimates or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and Finite-difference, time-domain ͑FDTD͒ calculations are typically performed with partial differential equations that are first order in time. Equation sets appropriate for FDTD calculations in a moving inhomogeneous medium ͑with an emphasis on the atmosphere͒ are derived and discussed in this paper. Two candidate equation sets, both derived from linearized equations of fluid dynamics, are proposed. The first, which contains three coupled equations for the sound pressure, vector acoustic velocity, and acoustic density, is obtained without any approximations. The second, which contains two coupled equations for the sound pressure and vector acoustic velocity, is derived by ignoring terms proportional to the divergence of the medium velocity and the gradient of the ambient pressure. It is shown that the second set has the same or a wider range of applicability than equations for the sound pressure that have been previously used for analytical and numerical studies of sound propagation in a moving atmosphere. Practical FDTD implementation of the second set of equations is discussed. Results show good agreement with theoretical predictions of the sound pressure due to a point monochromatic source in a uniform, high Mach number flow and with Fast Field Program calculations of sound propagation in a stratified moving atmosphere.

No abstract

Formulation and implementation of time-domain boundary conditions (TDBCs) at the surface of a reactive porous material are made challenging by the slow decay, complexity, or noncausal nature of many commonly used models of porous materials. In this paper, approaches are described that improve computational efficiency and enforce causality. One approach involves approximating the known TDBC for the modified Zwikker-Kosten impedance model as a summation of decaying exponential functions. A second approach, which can be applied to any impedance model, involves replacing the characteristic admittance with its Padé approximation. Then, approximating fractional derivatives with decaying exponentials, a causal and recursive TDBC is formulated.

No abstract

The classical wavefront method for interpreting seismic refraction arrival times is implemented on a digital computer, Modern finite-difference propagation algorithms are used to downward continue recorded refraction arrival times through a near-surface heterogeneousvelocity structure. Two such subsurface traveltime fields need to be reconstructed from the arrivals observed on a forward and reverse geophone spread.