Parameter Error Analysis for the 3D Modified Leray-alpha Model: Analytical and Numerical Approaches

Débora A.F. Albanez,Maicon J. Benvenutti,Samuel Little,Jing Tian

2024-11-25

Abstract:In this study, we conduct a parameter error analysis for the 3D modified Leray-$\alpha$ model using both analytical and numerical approaches. We first prove the global well-posedness and continuous dependence of initial data for the assimilated system. Furthermore, given sufficient conditions on the physical parameters and norms of the true solution, we demonstrate that the true solution can be recovered from the approximation solution, with an error determined by the discrepancy between the true and approximating parameters. Numerical simulations are provided to validate the convergence criteria.

Numerical Analysis,Mathematical Physics,Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to analyze the parameter errors in the three - dimensional modified Leray - α model (3D modified Leray - α model). Specifically, the authors explore, through analytical and numerical methods, how to use the Continuous Data Assimilation (CDA) technique to estimate the model parameter errors when the length - scale parameter \(\alpha>0\) is uncertain, and verify the error between the approximate solution and the true solution.

### Main problem decomposition:

1. **Global well - posedness and continuous dependence**:
- The author first proves the global well - posedness of the modified Leray - α model under given initial data and its continuous dependence on the initial data. This ensures the mathematical rationality of the model.
- The formula is expressed as follows:
\[
\begin{cases}
\frac{\partial v}{\partial t}-\nu\Delta v+(v\cdot\nabla)u + \nabla p=f,\\
v=(1-\alpha^{2}\Delta)u,\\
\nabla\cdot u=\nabla\cdot v = 0,\\
u(0)=u_{0},
\end{cases}
\]
where \(v\) and \(u\) are the non - filtered and filtered velocity fields respectively, \(\nu>0\) is the kinematic viscosity, \(f = f(x,t)\) is the external force, and \(p=p(x,t)\) is the scalar pressure field.

2. **Parameter error analysis**:
- Given sufficient conditions for physical parameters and the norm of the true solution, the author shows that the true solution can be recovered by the approximate solution, and the error depends on the difference between the true parameter \(\alpha\) and the approximate parameter \(\beta\).
- Using the CDA technique, a "guessed" parameter \(\beta>0\) is introduced to replace the unknown \(\alpha\), and the following CDA system is constructed:
\[
\begin{cases}
\frac{\partial z}{\partial t}-\nu\Delta z+(z\cdot\nabla)w+\nabla p=f-\eta(I-\beta^{2}\Delta)(I_{h}(w)-I_{h}(u)),\\
z=(1-\beta^{2}\Delta)w,\\
\nabla\cdot w=\nabla\cdot z = 0.
\end{cases}
\]

3. **Numerical simulation verification**:
- To verify the results of the theoretical analysis, the authors provide numerical simulation experiments. These simulations are carried out using a new Python package Dedalus to ensure flexibility and accuracy.
- The simulation results show the convergence and divergence between the solutions of the original system and the assimilation system under different parameter settings.

### Summary:
The main purpose of this paper is to analyze the parameter errors in the 3D modified Leray - α model through analytical and numerical methods, and prove that under appropriate conditions, the true solution can be recovered by the approximate solution, and the error depends on the difference between the true parameter and the approximate parameter. In addition, the numerical simulation verifies the effectiveness of the theoretical analysis.