Abstract:We design a Mixed Virtual Element Method for the approximated solution to the first-order form of the acoustic wave equation. In absence of external load, the semi-discrete method exactly conserves the system energy. To integrate in time the semi-discrete problem we consider a classical theta-method scheme. We carry out the stability and convergence analysis in the energy norm for the semi-discrete problem showing optimal rate of convergence with respect to the mesh size. We further study the property of energy conservation for the fully-discrete system. Finally, we present some verification tests as well as engineering application of the method.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a Mixed Virtual Element Method (MVEM) to solve the first - order form of the acoustic wave equation. Specifically, the main objectives of the paper include:
1. **Design of the discrete method**: Construct an MVEM framework for the spatial discretization of the acoustic wave equation on general polygonal meshes. This method can handle complex geometries and irregular meshing.
2. **Energy conservation**: Ensure that in the absence of external loads, the semi - discrete method can accurately maintain the energy conservation of the system. This is achieved by designing an appropriate discrete format.
3. **Time - integration scheme**: Adopt the classical θ - method scheme for time integration and analyze its stability and convergence. In particular, the energy conservation characteristics of the system under different θ values are studied.
4. **Stability and convergence analysis**: Conduct stability analysis under the energy norm and prove the optimal convergence rate. Specifically, the paper shows the optimal convergence order with respect to the mesh size.
5. **Numerical verification**: Verify the effectiveness and accuracy of the proposed method through a series of verification tests and engineering applications.
### Mathematical model
The acoustic wave equation considered in the paper can be represented as the following first - order hyperbolic partial differential equation system:
\[
\begin{cases}
\mathbf{u}_t-\nabla p = 0 & \text{in } \Omega\times(0,T], \\
c^{-2}p_t-\nabla\cdot\mathbf{u} = f & \text{in } \Omega\times(0,T],
\end{cases}
\]
where $\mathbf{u}$ is the particle velocity, $p$ is the acoustic pressure, $c$ is the characteristic wave speed in the medium, and $f$ is the source term.
The boundary conditions include:
- Dirichlet boundary condition: $p(x,t) = g_D(x,t)$ on $\Gamma_D\times(0,T]$,
- Neumann boundary condition: $\mathbf{u}(x,t)\cdot\mathbf{n} = g_N(x,t)$ on $\Gamma_N\times(0,T]$,
- Robin boundary condition: $\mathbf{u}(x,t)\cdot\mathbf{n}+\alpha^{-1}c^{-1}p(x,t) = g_R(x,t)$ on $\Gamma_R\times(0,T]$.
The initial conditions are:
\[
\begin{cases}
p(x,0) = p_0(x) & \text{in } \Omega, \\
\mathbf{u}(x,0) = \mathbf{u}_0(x) & \text{in } \Omega.
\end{cases}
\]
### Main contributions
- **Energy conservation**: The paper proves that when there is no external force and the Robin boundary is empty, the system is energy - conserving.
- **Optimal convergence**: By introducing the Fortin operator and performing interpolation estimates, it is proved that the method has the optimal convergence order.
- **Energy characteristics of the time - integration scheme**: Analyze the energy conservation characteristics under different θ values, especially the Crank - Nicolson method ($\theta = 1/2$) can maintain the discrete energy conservation.
In summary, this paper aims to provide an efficient and accurate numerical method for solving the acoustic wave equation and ensure its stability and energy conservation characteristics under complex geometries.