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 propose an approximate model for the 2D Kuramoto-Sivashinsky equations (KSE) of flame fronts and crystal growth. We prove that this new ``calmed'' version of the KSE is globally well-posed, and moreover, its solutions converge to solutions of the KSE on the time interval of existence and uniqueness of the KSE at an algebraic rate. In addition, we provide simulations of the calmed KSE, illuminating its dynamics. These simulations also indicate that our analytical predictions of the convergence rates are sharp. We also discuss analogies with the 3D Navier-Stokes equations of fluid dynamics., Comment: 35 pages, 8 figures

We adapt a continuous data assimilation scheme, known as the Azouani-Olson-Titi (AOT) algorithm, to the case of moving observers for the 2D incompressible Navier-Stokes equations. We propose and test computationally several movement patterns (which we refer to as "the bleeps, the sweeps and the creeps"), as well as Lagrangian motion and combinations of these patterns, in comparison with static (i.e. non-moving) observers. In several cases, order-of-magnitude improvements in terms of the time-to-convergence are observed. We end with a discussion of possible applications to real-world data collection strategies that may lead to substantial improvements in predictive capabilities.

We introduce a general and fast convolution-based method (FCBM) for peridynamics (PD). Expressing the PD integrals in terms of convolutions and computing them by fast Fourier transform (FFT), we reduce the computational complexity of PD models from O(N^2) to O(Nlog_2 N), with N being the total number of discretization nodes. Initial neighbor identification and storing neighbor information is not required, and, as a consequence, memory allocation scales with O(N) instead of O(N^2), common for existing methods. The method is applicable to bounded domains with arbitrary shapes and boundary conditions via an embedded constraint (EC) approach. We explain the FCBM-EC formulation for certain bond-based and state-based, linear and nonlinear PD models of elasticity and dynamic brittle fracture, as applications. We solve a 3D elastostatic problem and show that the FCBM reduces the computational time from days to hours and from years to days, compared with the original meshfree discretization for PD models. Large-scale computations of PD models are feasible with the new method, and we demonstrate its versatility by simulating, with ease, the difficult problem of multiple crack branching in a brittle plate.

We study a continuous data assimilation (CDA) algorithm for a velocity-vorticity formulation of the 2D Navier-Stokes equations in two cases: nudging applied to the velocity and vorticity, and nudging applied to the velocity only. We prove that under a typical finite element spatial discretization and backward Euler temporal discretization, application of CDA preserves the unconditional long-time stability property of the velocity-vorticity method and provides optimal long-time accuracy. These properties hold if nudging is applied only to the velocity, and if nudging is also applied to the vorticity then the optimal long-time accuracy is achieved more rapidly in time. Numerical tests illustrate the theory, and show its effectiveness on an application problem of channel flow past a flat plate., Comment: 28 pages, 8 figures, 2 tables

We introduce an efficient boundary-adapted spectral method for peridynamic diffusion problems with arbitrary boundary conditions. The spectral approach transforms the convolution integral in the peridynamic formulation into a multiplication in the Fourier space, resulting in computations that scale as O(NlogN). The limitation of regular spectral methods to periodic problems is eliminated using the volume penalization method. We show that arbitrary boundary conditions or volume constraints can be enforced in this way to achieve high levels of accuracy. To test the performance of our approach we compare the computational results with analytical solutions of the nonlocal problem. The performance is tested with convergence studies in terms of nodal discretization and the size of the penalization parameter in problems with Dirichlet and Neumann boundary conditions., Comment: 26 pages, 12 figures

We propose a new blow-up criterion for the 3D Euler equations of incompressible fluid flows, based on the 3D Euler-Voigt inviscid regularization. This criterion is similar in character to a criterion proposed in a previous work by the authors, but it is stronger, and better adapted for computational tests. The 3D Euler-Voigt equations enjoy global well-posedness, and moreover are more tractable to simulate than the 3D Euler equations. A major advantage of these new criteria is that one only needs to simulate the 3D Euler-Voigt, and not the 3D Euler equations, to test the blow-up criteria, for the 3D Euler equations, computationally.

We propose a new blow-up criterion for the 3D Euler equations of incompressible fluid flows, based on the 3D Euler-Voigt inviscid regularization. This criterion is similar in character to a criterion proposed in a previous work by the authors, but it is stronger, and better adapted for computational tests. The 3D Euler-Voigt equations enjoy global well-posedness, and moreover are more tractable to simulate than the 3D Euler equations. A major advantage of these new criteria is that one only needs to simulate the 3D Euler-Voigt, and not the 3D Euler equations, to test the blow-up criteria, for the 3D Euler equations, computationally.

In light of the question of finite-time blow-up vs. global well-posedness of solutions to problems involving nonlinear partial differential equations, we provide several cautionary examples which indicate that modifications to the boundary conditions or to the nonlinearity of the equations can effect whether the equations develop finite-time singularities. In particular, we aim to underscore the idea that in analytical and computational investigations of the blow-up of three-dimensional Euler and Navier-Stokes equations, the boundary conditions may need to be taken into greater account. We also examine a perturbation of the nonlinearity by dropping the advection term in the evolution of the derivative of the solutions to the viscous Burgers equation, which leads to the development of singularities not present in the original equation, and indicates that there is a regularizing mechanism in part of the nonlinearity. This simple analytical example corroborates recent computational observations in the singularity formation of fluid equations.

In light of the question of finite-time blow-up vs. global well-posedness of solutions to problems involving nonlinear partial differential equations, we provide several cautionary examples which indicate that modifications to the boundary conditions or to the nonlinearity of the equations can effect whether the equations develop finite-time singularities. In particular, we aim to underscore the idea that in analytical and computational investigations of the blow-up of three-dimensional Euler and Navier-Stokes equations, the boundary conditions may need to be taken into greater account. We also examine a perturbation of the nonlinearity by dropping the advection term in the evolution of the derivative of the solutions to the viscous Burgers equation, which leads to the development of singularities not present in the original equation, and indicates that there is a regularizing mechanism in part of the nonlinearity. This simple analytical example corroborates recent computational observations in the singularity formation of fluid equations.