Daniel G. Alfaro Vigo,Daniele C.R. Gomes,Bruno A. do Carmo,Mauro A. Rincon
Abstract:In this article, we present the mathematical analysis of the convergence of the linearized Crank-Nicolson Galerkin method for a nonlinear Schrodinger problem related to a domain with a moving boundary. The convergence analysis of the numerical method is carried out for both semi-discrete and fully discrete problems. An optimal error estimate in the $L^2$-norm with order ${O}(\tau^2+ h^s),~ 2\leq s\leq r$, where $h$ is the finite element mesh size parameter, $\tau$ is the time step, and $r-1$ represents the degree of the finite element polynomial basis. Numerical simulations are provided to confirm the consistency between theoretical and numerical results, validating the method and the order of convergence for different degrees $p\geq 1$ of the Lagrange polynomials and also for Hermite polynomials (degree $p=3$), which form the basis of the approximate solution.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of convergence analysis of numerical methods related to the nonlinear Schrödinger equation with moving boundaries. Specifically, the authors focus on the convergence of the linearized Crank - Nicolson Galerkin method in solving such equations, and provide theoretical analysis and numerical verification.
#### Main research contents
1. **Mathematical model**: Consider a nonlinear Schrödinger equation defined on a non - cylindrical domain:
\[
\begin{cases}
u'(x, t) - i\Delta u(x, t) + |u(x, t)|^\rho u(x, t) = \hat{f}(x, t) & \text{in } \hat{Q}, \\
u(x, t) = 0 & \text{on } \hat{\Sigma}, \\
u(x, 0) = u_0(x) & \text{in } \Omega_0,
\end{cases}
\]
where \( u' \) represents the partial derivative with respect to time \( t \), \( i \) is the imaginary unit, and \( \Delta \) is the Laplace operator.
2. **Numerical method**: The linearized Crank - Nicolson Galerkin method is used to discretize the above problem. This method combines finite - element space discretization and the Crank - Nicolson time - stepping scheme.
3. **Convergence analysis**: Convergence analysis is carried out for the semi - discrete and fully - discrete problems. For the semi - discrete problem, the variational form is mainly used for analysis; for the fully - discrete problem, a time - dependent Ritz projection is introduced, and the nonlinear term is dealt with by the error - splitting technique.
4. **Error estimate**: The optimal \( L^2 \)-norm error estimate, namely \( O(\tau^2 + h^s) \), is established, where \( \tau \) is the time step size, \( h \) is the finite - element mesh size parameter, and \( r - 1 \) represents the degree of the finite - element polynomial basis function.
5. **Numerical simulation**: The theoretical results are verified by numerical experiments, and the convergence orders under different - order Lagrange polynomials (\( p\geq1 \)) and Hermite polynomials (\( p = 3 \)) are shown.
#### Research significance
This research fills the gap in the previous theoretical analysis of numerical methods for the nonlinear Schrödinger equation with moving boundaries. Although the Schrödinger equation on cylindrical domains has been extensively studied in the existing literature, it is not easy to directly apply these methods to non - cylindrical domains. Therefore, through the adaptive adjustment of existing techniques, this paper provides a systematic theoretical framework and numerical verification, laying the foundation for further research in related fields.