Optimal Neumann boundary and distributed control of the Westervelt equation with time-fractional attenuation
Vanja Nikolić, Belkacem Said-Houari
TL;DR
The paper addresses optimal control of nonlinear acoustic waves modeled by the Westervelt equation with time-fractional attenuation, encompassing both Neumann boundary and distributed controls. It develops a rigorous forward analysis, including two regularity regimes and a fixed-point argument to handle the quasilinear, fractionally damped state, and it extends well-posedness to inhomogeneous Neumann data via a data-extension technique. Existence of globally optimal controls is proven, with stability results under perturbations of the target pressure and vanishing regularization, and adjoint-based first-order optimality conditions are derived, despite the state-dependent adjoint coefficients. The results provide a solid theoretical foundation for adjoint-based optimization in fractionally damped Westervelt models, with potential guidance for numerical algorithms and applications in medical ultrasound where precise energy deposition is essential.
Abstract
Optimal control of nonlinear acoustic waves is relevant in many medical ultrasound technologies, ranging from cancer therapy to targeted drug delivery, where it can help guide the precise deposition of acoustic energy. In this work, we study Neumann boundary and distributed control problems for tracking a prescribed pressure field governed by the Westervelt equation with time-fractional dissipation. This model captures nonlinear ultrasonic wave propagation in biological media and accounts for the experimentally observed power-law attenuation. We begin by extending the existing well-posedness theory for time-fractional equations to include inhomogeneous Neumann boundary data used as control inputs, which requires constructing an appropriate data extension and regularization. Using these analytical results for the forward problem, we prove the existence of globally optimal controls and analyze the stability of the optimization problem with respect to perturbations in the target pressure field and to vanishing regularization parameters. Finally, we investigate the associated adjoint equation, which has state-dependent coefficients, and use it to derive first-order necessary optimality conditions.
