A mimetic discretization of Westervelt's equation

William Barham,Philip J. Morrison
2024-07-20
Abstract:A broad class of nonlinear acoustic wave models possess a Hamiltonian structure in their dissipation-free limit and a gradient flow structure for their dissipative dynamics. This structure may be exploited to design numerical methods which preserve the Hamiltonian structure in the dissipation-free limit, and which achieve the correct dissipation rate in the spatially-discrete dissipative dynamics. Moreover, by using spatial discretizations which preserve the de Rham cohomology, the non-evolving involution constraint for the vorticity may be exactly satisfied for all of time. Numerical examples are given using a mimetic finite difference spatial discretization.
Numerical Analysis,Mathematical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to design a numerical method to preserve the Hamiltonian structure of the nonlinear acoustic wave model in the non - dissipative limit and achieve the correct dissipation rate in the spatially - discrete dissipative dynamics. Specifically, the paper focuses on the structure - preserving discretization method of the Westervelt equation. By using the spatially - discrete method that preserves de Rham cohomology, it can be ensured that the non - evolving constraint of vorticity is exactly satisfied throughout time. ### Main contributions of the paper 1. **Preservation of Hamiltonian structure**: - The paper proposes a framework based on mimetic discretization, which can preserve the Hamiltonian structure in the non - dissipative limit. - Through this method, it can be ensured that the numerical method preserves the known properties of the continuous model when dealing with the nonlinear acoustic wave model. 2. **Correct implementation of dissipation rate**: - By designing symmetric positive - definite brackets, the paper achieves the correct dissipation rate of the dissipative dynamics. - This method ensures that energy dissipates monotonically in numerical simulations, which is in line with physical expectations. 3. **Vorticity conservation**: - By using the spatially - discrete method that preserves de Rham cohomology, the paper ensures that the non - evolving constraint of vorticity is exactly satisfied throughout time. ### Specific content of the paper 1. **Continuous Hamiltonian theory**: - Section 2 details the continuous Hamiltonian structure of the nonlinear acoustic wave model, especially the Hamiltonian structure of the Westervelt equation. - The paper discusses how to represent these models as Hamiltonian systems and introduces the gradient - flow structure of dissipative dynamics. 2. **Discrete Hamiltonian theory**: - Section 3 presents the general framework for structure - preserving discretization of the nonlinear acoustic wave model. - The paper discusses how to use the mimetic discretization method to preserve the Hamiltonian structure and gives a specific discretization scheme. 3. **Structure - preserving discretization of the Westervelt equation**: - Section 4 specifically discusses the structure - preserving discretization method of the Westervelt equation. - The paper details how to use the mimetic finite - difference method for spatial discretization and gives the method of time integration. 4. **Numerical examples**: - Section 5 provides multiple numerical examples to verify the convergence of the method and the preservation of the properties of the continuous model. - It includes one - dimensional and two - dimensional examples, showing the performance of the method in different situations. ### Conclusion The paper successfully designs a structure - preserving numerical method for the numerical simulation of the nonlinear acoustic wave model (especially the Westervelt equation). This method not only preserves the Hamiltonian structure but also achieves the correct dissipation rate and ensures vorticity conservation. These characteristics make this method have important application value in the study of nonlinear acoustic waves.