Piecewise Acoustic Source Imaging with Unknown Speed of Sound Using a Level-Set Method

We investigate the following inverse problem:starting from the acoustic wave equation,reconstruct a piecewise constant passive acoustic source from a single boundary temporal measurement without knowing the speed of sound.When the amplitudes of the source are known a priori,we prove a unique determination result of the shape and propose a level set algorithm to reconstruct the singularities.When the singularities of the source are known a priori,we show unique determination of the source amplitudes and propose a least-squares fitting algorithm to recover the source amplitudes.The analysis bridges the low-frequency source inversion problem and the inverse problem of gravimetry.The proposed algorithms are validated and quantitatively evaluated with numerical experiments in 2D and 3D.

A Fast Local Level Set Method for Inverse Gravimetry

We propose a fast local level set method for the inverse problem of gravimetry.The theoretical foundation for our approach is based on the following uniqueness result:if an open set D is star-shaped or x3-convex with respect to its center of gravity,then its exterior potential uniquely determines the open set D.To achieve this purpose constructively,the first challenge is how to parametrize this open set D as its boundary may have a variety of possible shapes.To describe those different shapes we propose to use a level-set function to parametrize the unknown boundary of this open set.The second challenge is how to deal with the issue of partial data as gravimetric measurements are only made on a part of a given reference domainΩ.To overcome this difficulty,we propose a linear numerical continuation approach based on the single layer representation to find potentials on the boundary of some artificial domain containing the unknown set D.The third challenge is how to speed up the level set inversion process.Based on some features of the underlying inverse gravimetry problem such as the potential density being constant inside the unknown domain,we propose a novel numerical approach which is able to take advantage of these features so that the computational speed is accelerated by an order of magnitude.We carry out numerical experiments for both two-and three-dimensional cases to demonstrate the effectiveness of the new algorithm.

A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS

We propose an effective stopping criterion for higher-order fast sweeping schemes for static Hamilton-Jacobi equations based on ratios of three consecutive iterations. To design the new stopping criterion we analyze the convergence of the first-order Lax-Friedrichs sweeping scheme by using the theory of nonlinear iteration. In addition, we propose a fifth-order Weighted PowerENO sweeping scheme for static Hamilton-Jacobi equations with convex Hamiltonians and present numerical examples that validate the effectiveness of the new stopping criterion.

An Adjoint StateMethod for Numerical Approximation of Continuous Traffic Congestion Equilibria

The equilibriummetric forminimizing a continuous congested trafficmodel is the solution of a variational problem involving geodesic distances.The continuous equilibrium metric and its associated variational problem are closely related to the classical discrete Wardrop’s equilibrium.We propose an adjoint state method to numerically approximate continuous traffic congestion equilibria through the continuous formulation.The method formally derives an adjoint state equation to compute the gradient descent direction so as to minimize a nonlinear functional involving the equilibrium metric and the resulting geodesic distances.The geodesic distance needed for the state equation is computed by solving a factored eikonal equation,and the adjoint state equation is solved by a fast sweeping method.Numerical examples demonstrate that the proposed adjoint state method produces desired equilibrium metrics and outperforms the subgradient marching method for computing such equilibrium metrics.