Abstract:We propose a self-adaptive absorbing technique for quasilinear ultrasound waves in two- and three-dimensional computational domains. As a model for the nonlinear ultrasound propagation in thermoviscous fluids, we employ Westervelt's wave equation solved for the acoustic velocity potential. The angle of incidence of the wave is computed based on the information provided by the wave-field gradient which is readily available in the finite element framework. The absorbing boundary conditions are then updated with the angle values in real time. Numerical experiments illustrate the accuracy and efficiency of the proposed method.
What problem does this paper attempt to address?
This paper aims to solve the problem of boundary reflection in nonlinear ultrasonic propagation simulation. Specifically, it attempts to develop an adaptive absorbing boundary condition (ABC) technique to effectively reduce spurious reflections at the truncated boundaries of the computational domain, thereby improving the accuracy and efficiency of numerical simulations.
### Main Problems in the Paper
1. **Spurious Reflection Problem**:
- When performing numerical simulations of nonlinear ultrasonic propagation, in order to simplify the calculation, it is usually necessary to truncate the physical space.
- The boundaries of the truncated computational domain will introduce spurious reflections, which will affect the accuracy of the simulation results.
2. **Angle Sensitivity**:
- Traditional absorbing boundary conditions are very sensitive to the angle of incidence. When the range of the angle of incidence is large, their performance will decline significantly.
- This sensitivity limits the effectiveness of traditional methods in practical applications, especially in cases of complex geometries and multi - directional propagation.
### Solutions
To solve the above problems, the author proposes an adaptive absorbing boundary condition technique, whose main features include:
- **Real - Time Update Based on Gradient Information**:
- Calculate the wave - field gradient to determine the direction of the local wave vector and update the angle - of - incidence information in the absorbing boundary condition in real time.
- This method can dynamically adjust the boundary condition at each time step, so as to better adapt to different angles of incidence.
- **Applicable to Two - Dimensional and Three - Dimensional Computational Domains**:
- The proposed adaptive technique is applicable not only to two - dimensional computational domains but also to three - dimensional computational domains, ensuring the possibility of its wide application.
- **Combined with the Westervelt Equation**:
- Use the classic quasilinear acoustic model - the Westervelt equation as the model equation, considering the form of the acoustic velocity potential and the influence of non - zero angles of incidence.
- Through appropriate linearization of the Westervelt equation, derive the absorbing boundary condition applicable to a given angle of incidence.
### Mathematical Formula Representation
The key formulas involved in the paper are as follows:
- Westervelt Equation:
\[
\frac{1}{c^2}\psi_{tt}-\Delta\psi - \delta\Delta\psi_t=k\psi_t\psi_{tt}
\]
where:
\[
\delta = \frac{b}{c^2}, \quad k=\frac{1}{c^2}\left(\frac{B}{A}+2\right)
\]
- Adaptive Absorbing Boundary Condition (Considering Loss Terms):
\[
c^2\frac{\partial\psi}{\partial n}+b\frac{\partial\psi_t}{\partial n}=-c\sqrt{1-\frac{k}{2}\psi_t}\psi_t\cos\theta
\]
- Calculation of the Cosine of the Angle of Incidence:
\[
\cos\theta=\frac{|\nabla\psi\cdot\mathbf{n}|}{\sqrt{\psi_x^2+\psi_y^2+\psi_z^2}}
\]
Through these methods and techniques, the paper successfully improves the accuracy and efficiency of nonlinear ultrasonic propagation simulation and reduces the influence of boundary reflection on the results.