Search

We present algorithms for solving high-frequency acoustic scattering problems in complex domains. The eikonal and transport partial differential equations from the WKB/geometric optic approximation of the Helmholtz equation are solved recursively to generate boundary conditions for a tree of eikonal/transport equation pairs, describing the phase and amplitude of a geometric optic wave propagating in a complicated domain, including reflection and diffraction. Edge diffraction is modeled using the uniform theory of diffraction. For simplicity, we limit our attention to domains with piecewise linear boundaries and a constant speed of sound. The domain is discretized into a conforming tetrahedron mesh. For the eikonal equation, we extend the jet marching method to tetrahedron meshes. Hermite interpolation enables second order accuracy for the eikonal and its gradient and first order accuracy for its Hessian, computed using cell averaging. To march the eikonal on an unstructured mesh, we introduce a new method of rejecting unphysical updates by considering Lagrange multipliers and local visibility. To handle accuracy degradation near caustics, we introduce several fast Lagrangian initialization algorithms. We store the dynamic programming plan uncovered by the marcher in order to propagate auxiliary quantities along characteristics. We introduce an approximate origin function which is computed using the dynamic programming plan, and whose 1/2-level set approximates the geometric optic shadow and reflection boundaries. We also use it to propagate geometric spreading factors and unit tangent vector fields needed to compute the amplitude and evaluate the high-frequency edge diffraction coefficient. We conduct numerical tests on a semi-infinite planar wedge to evaluate the accuracy of our method. We also show an example with a more realistic building model with challenging architectural features. Finally, we demonstrate a simple approach to extending the method to handle nonconstant speeds of sound by modifying the semi-Lagrangian updates to account for a varying speed. • A second-order marcher for solving the eikonal equation on a tetrahedron mesh is introduced. • Algorithms for multiple scattering due to reflecting/diffracting acoustic waves in 3D. • Dynamic programming for UTD-based edge diffraction. • Dynamic programming is used to approximate geometric optic reflection and shadow boundaries. • Real-world numerical examples of sound wave propagation in complicated environments. [ABSTRACT FROM AUTHOR]

Abstract Nongradient SDEs with small white noise often arise when modeling biological and ecological time-irreversible processes. If the governing SDE were gradient, the maximum likelihood transition paths, transition rates, expected exit times, and the invariant probability distribution would be given in terms of its potential function. The quasipotential plays a similar role for nongradient SDEs. Unfortunately, the quasipotential is the solution of a functional minimization problem that can be obtained analytically only in some special cases. We propose a Dijkstra-like solver for computing the quasipotential on regular rectangular meshes in 3D. This solver results from a promotion and an upgrade of the previously introduced ordered line integral method with the midpoint quadrature rule for 2D SDEs. The key innovations that have allowed us to keep the CPU times reasonable while maintaining good accuracy are (i) a new hierarchical update strategy, (i i) the use of Karush–Kuhn–Tucker theory for rejecting unnecessary simplex updates, and (i i i) pruning the number of admissible simplexes and a fast search for them. An extensive numerical study is conducted on a series of linear and nonlinear examples where the quasipotential is analytically available or can be found at transition states by other methods. In particular, the proposed solver is applied to Tao's examples where the transition states are hyperbolic periodic orbits, and to a genetic switch model by Lv et al. (2014) [21]. The C source code implementing the proposed algorithm is available at M. Cameron's web page. Highlights • Ordered line integral method with midpoint rule for finding the quasipotential in 3D. • Using the KKT theory to reject unnecessary simplex updates. • An upgraded hierarchical update strategy to reduce the number of simplex updates. • Pruning the number of admissible simplexes and a fast search for them. • Testing the solver on a genetic switch model by Lv et al. [ABSTRACT FROM AUTHOR]