Abstract:In this paper we discuss energy conservation issues related to the numerical solution of the nonlinear wave equation, when a Fourier expansion is considered for the space discretization. The obtained semi-discrete problem is then solved in time by means of energy-conserving Runge-Kutta methods in the HBVMs class.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is related to the energy conservation in the numerical solutions of nonlinear wave equations. Specifically, the author focuses on how to ensure that the numerical method can accurately preserve the energy - conservation characteristics of the system when using Fourier expansion for spatial discretization. This is not limited to nonlinear wave equations and can also be extended to other types of Hamiltonian partial differential equations (Hamiltonian PDEs).
### Main Problem Description
The core of the paper's research is the energy conservation problem of nonlinear wave equations, especially how to preserve energy conservation during the time - integration process when using Fourier expansion for spatial discretization. The specific equation is as follows:
\[
u_{tt}(x,t)=\alpha^{2}u_{xx}(x,t)-f'(u(x,t)),\quad(x,t)\in(0,1)\times(0,\infty)
\]
The initial conditions are:
\[
u(x,0)=\psi_{0}(x),\quad u_{t}(x,0)=\psi_{1}(x),\quad x\in(0,1)
\]
And periodic boundary conditions, Dirichlet boundary conditions and Neumann boundary conditions are considered.
### Research Objectives
1. **Energy Conservation**: Ensure that the numerical method can accurately preserve the energy - conservation characteristics of the system.
2. **Selection of Numerical Method**: Use the energy - conserving Runge - Kutta method (such as HBVMs - type methods) to solve the semi - discretized Hamiltonian system.
3. **Handling of Different Boundary Conditions**: Research how to effectively apply these methods under different types of boundary conditions (periodic, Dirichlet, Neumann).
### Solutions
- **Spatial Discretization**: Use the Fourier - Galerkin spectral method for spatial discretization.
- **Time Integration**: Use the energy - conserving Runge - Kutta method (such as HBVMs - type methods) for time integration.
- **Boundary Condition Handling**: For different boundary conditions (periodic, Dirichlet, Neumann), corresponding handling methods are proposed to ensure the preservation of energy - conservation characteristics.
Through these methods, the author aims to develop an efficient algorithm that can numerically and accurately preserve energy - conservation characteristics, thereby better simulating the behavior of Hamiltonian partial differential equations.