Search

We study different approaches to implementing sparse-in-time observations into the the Azouani-Olson-Titi data assimilation algorithm. We propose a new method which introduces a "data assimilation window" separate from the observational time interval. We show that by making this window as small as possible, we can drastically increase the strength of the nudging parameter without losing stability. Previous methods used old data to nudge the solution until a new observation was made. In contrast, our method stops nudging the system almost immediately after an observation is made, allowing the system relax to the correct physics. We show that this leads to an order-of-magnitude improvement in the time to convergence in our 3D Navier-Stokes simulations. Moreover, our simulations indicate that our approach converges at nearly the same rate as the idealized method of direct replacement of low Fourier modes proposed by Hayden, Olson, and Titi (HOT). However, our approach can be readily adapted to non-idealized settings, such as finite element methods, finite difference methods, etc., since there is no need to access Fourier modes as our method works for general interpolants. It is in this sense that we think of our approach as ``second best;'' that is, the ``best'' method would be the direct replacement of Fourier modes as in HOT, but this idealized approach is typically not feasible in physically realistic settings. While our method has a convergence rate that is slightly sub-optimal compared to the idealized method, it is directly compatible with real-world applications. Moreover, we prove analytically that these new algorithms are globally well-posed, and converge to the true solution exponentially fast in time. In addition, we provide the first 3D computational validation of HOT algorithm., Comment: 25 pages, 2 figures consisting of 13 subfigures

We study a continuous data assimilation algorithm proposed by Azouani, Olson, and Titi (AOT) in the context of an unknown Reynolds number. We determine the large-time error between the true solution of the 2D Navier-Stokes equations and the assimilated solution due to discrepancy between an approximate Reynolds number and the physical Reynolds number. Additionally, we develop an algorithm that can be run in tandem with the AOT algorithm to recover both the true solution and the Reynolds number (or equivalently the true viscosity) using only spatially discrete velocity measurements. The algorithm we propose involves changing the viscosity mid-simulation. Therefore, we also examine the sensitivity of the equations with respect to the Reynolds number. We prove that a sequence of difference quotients with respect to the Reynolds number converges to the unique solution of the sensitivity equations for both the 2D Navier-Stokes equations and the assimilated equations. We also note that this appears to be the first such rigorous proof of existence and uniqueness to strong or weak solutions to the sensitivity equations for the 2D Navier-Stokes equations (in the natural case of zero initial data), and that they can be obtained as a limit of difference quotients with respect to the Reynolds number., Comment: 39 pages, 5 figures, 1 table

Data assimilation is a technique for increasing the accuracy of simulations of solutions to partial differential equations by incorporating observable data into the solution as time evolves. Recently, a promising new algorithm for data assimilation based on feedback-control at the PDE level has been proposed in the pioneering work of Azouani, Olson, and Titi (2014). The standard version of this algorithm is based on measurement from data points that are fixed in space. In this work, we consider the scenario in which the data collection points move in space over time. We demonstrate computationally that, at least in the setting of the 1D Allen-Cahn reaction diffusion equations, the algorithm converges with significantly fewer measurement points, up to an order or magnitude in some cases. We also provide an application of the algorithm to an inverse problem in the case of a uniform static grid.

We introduce three new nonlinear continuous data assimilation algorithms. These models are compared with the linear continuous data assimilation algorithm introduced by Azouani, Olson, and Titi (AOT). As a proof-of-concept for these models, we computationally investigate these algorithms in the context of the 1D Kuramoto-Sivashinsky equation. We observe that the nonlinear models experience super-exponential convergence in time, and converge to machine precision significantly faster than the linear AOT algorithm in our tests., Comment: 15 pages, 19 figures

We establish global existence and uniqueness theorems for the two-dimensional non-diffusive Boussinesq system with viscosity only in the horizontal direction, which arises in Ocean dynamics. This work improves the global well-posedness results established recently by R. Danchin and M. Paicu for the Boussinesq system with anisotropic viscosity and zero diffusion. Although we follow some of their ideas, in proving the uniqueness result, we have used an alternative approach by writing the transported temperature (density) as $��= ����$ and adapting the techniques of V. Yudovich for the 2D incompressible Euler equations. This new idea allows us to establish uniqueness results with fewer assumptions on the initial data for the transported quantity $��$. Furthermore, this new technique allows us to establish uniqueness results without having to resort to the paraproduct calculus of J. Bony. We also propose an inviscid $��$-regularization for the two-dimensional inviscid, non-diffusive Boussinesq system of equations, which we call the Boussinesq-Voigt equations. Global regularity of this system is established. Moreover, we establish the convergence of solutions of the Boussinesq-Voigt model to the corresponding solutions of the two-dimensional Boussinesq system of equations for inviscid flow without heat (density) diffusion on the interval of existence of the latter. Furthermore, we derive a criterion for finite-time blow-up of the solutions to the inviscid, non-diffusive 2D Boussinesq system based on this inviscid Voigt regularization. Finally, we propose a Voigt-$��$ regularization for the inviscid 3D Boussinesq equations with diffusion, and prove its global well-posedness. It is worth mentioning that our results are also valid in the presence of the $��$-plane approximation of the Coriolis force.