Existence, uniqueness, and numerical solutions of the nonlinear periodic Westervelt equation

In this paper, we study the nonlinear periodic Westervelt equation with excitations located within a bounded domain in \mathbb{R}^d, where d \in \{2,3\}, subject to Robin boundary conditions. This problem is of particular interest for advancing imaging techniques %in invasive ultrasound, such as intravascular ultrasound imaging. that exploit nonlinearity of the acoustic propagation. We establish the existence and uniqueness of solutions in both the linear and the nonlinear setting, thereby allowing for spatially varying coefficients as relevant in quantitative imaging. %, highlighting the crucial role of small excitations in ensuring existence and uniqueness in the nonlinear setting. Derivation of a multiharmonic formulation enables us to show the generation of higher harmonics (that is, responses at multiples of the fundamental frequency) due nonlinear wave propagation. An iterative scheme for solving the resulting system is proposed that relies on successive resolution of these higher harmonics, and its convergence under smallness conditions on the excitation is proven. Furthermore, we investigate the numerical solution of the resulting system of Helmholtz equations, employing a conforming finite element method for its discretization. Through an implementation of the proposed methodology, we illustrate how acoustic waves propagate in nonlinear media. This study aims to enhance our understanding of ultrasound propagation dynamics, which is essential for obtaining high-quality images from limited in vivo and boundary measurements.