Search

This paper considers continuous data assimilation (CDA) in partial differential equation (PDE) discretizations where nudging parameters can be taken arbitrarily large. We prove that long-time optimally accurate solutions are obtained for such parameters for the heat and Navier-Stokes equations (using implicit time stepping methods), with error bounds that do not grow as the nudging parameter gets large. Existing theoretical results either prove optimal accuracy but with the error scaled by the nudging parameter, or suboptimal accuracy that is independent of it. The key idea to the improved analysis is to decompose the error based on a weighted inner product that incorporates the (symmetric by construction) nudging term, and prove that the projection error from this weighted inner product is optimal and independent of the nudging parameter. We apply the idea to BDF2 - finite element discretizations of the heat equation and Navier-Stokes equations to show that with CDA, they will admit optimal long-time accurate solutions independent of the nudging parameter, for nudging parameters large enough. Several numerical tests are given for the heat equation, fluid transport equation, Navier-Stokes, and Cahn-Hilliard that illustrate the theory.

This paper extends a low-rank tensor decomposition (LRTD) reduced order model (ROM) methodology to simulate viscous flows and in particular to predict a smooth branch of solutions for the incompressible Navier-Stokes equations. Additionally, it enhances the LRTD-ROM methodology by introducing a non-interpolatory variant, which demonstrates improved accuracy compared to the interpolatory method utilized in previous LRTD-ROM studies. After presenting both the interpolatory and non-interpolatory LRTD-ROM, we demonstrate that with snapshots from a few different viscosities, the proposed method is able to accurately predict flow statistics in the Reynolds number range $[25,400]$. This is a significantly wider and higher range than state of the art (and similar size) ROMs built for use on varying Reynolds number have been successful on. The paper also discusses how LRTD may offer new insights into the properties of parametric solutions.

The purpose of this paper is to develop a practical strategy to accelerate Newton's method in the vicinity of singular points. We present an adaptive safeguarding scheme with a tunable parameter, which we call adaptive gamma-safeguarding, that one can use in tandem with Anderson acceleration to improve the performance of Newton's method when solving problems at or near singular points. The key features of adaptive gamma-safeguarding are that it converges locally for singular problems, and it can detect nonsingular problems automatically, in which case the Newton-Anderson iterates are scaled towards a standard Newton step. The result is a flexible algorithm that performs well for singular and nonsingular problems, and can recover convergence from both standard Newton and Newton-Anderson with the right parameter choice. This leads to faster local convergence compared to both Newton's method, and Newton- Anderson without safeguarding, with effectively no additional computational cost. We demonstrate three strategies one can use when implementing Newton-Anderson and gamma-safeguarded Newton-Anderson to solve parameter-dependent problems near singular points. For our benchmark problems, we take two parameter-dependent incompressible flow systems: flow in a channel and Rayleigh-Benard convection., Comment: 22 pages, 14 figures

We consider local balances of momentum and angular momentum for the incompressible Navier-Stokes equations. First, we formulate new weak forms of the physical balances (conservation laws) of these quantities, and prove they are equivalent to the usual conservation law formulations. We then show that continuous Galerkin discretizations of the Navier-Stokes equations using the EMAC form of the nonlinearity preserve discrete analogues of the weak form conservation laws, both in the Eulerian formulation and the Lagrangian formulation (which are not equivalent after discretizations). Numerical tests illustrate the new theory.

This paper considers improving the Picard and Newton iterative solvers for the Navier-Stokes equations in the setting where data measurements or solution observations are available. We construct adapted iterations that use continuous data assimilation (CDA) style nudging to incorporate the known solution data into the solvers. For CDA-Picard, we prove the method has an improved convergence rate compared to usual Picard, and the rate improves as more measurement data is incorporated. We also prove that CDA-Picard is contractive for larger Reynolds numbers than usual Picard, and the more measurement data that is incorporated the larger the Reynolds number can be with CDA-Picard still being contractive. For CDA-Newton, we prove that the domain of convergence, with respect to both the initial guess and the Reynolds number, increases as the amount of measurement data is increased. Additionally, for both methods we show that CDA can be implemented as direct enforcement of measurement data into the solution. Numerical results for common benchmark Navier-Stokes tests illustrate the theory., Comment: 22 pages

We study continuous data assimilation (CDA) applied to projection and penalty methods for the Navier-Stokes (NS) equations. Penalty and projection methods are more efficient than consistent NS discretizations, however are less accurate due to modeling error (penalty) and splitting error (projection). We show analytically and numerically that with measurement data and properly chosen parameters, CDA can effectively remove these splitting and modeling errors and provide long time optimally accurate solutions., Comment: 21 pages, 15 figures

This work introduces, analyzes and demonstrates an efficient and theoretically sound filtering strategy to ensure the condition of the least-squares problem solved at each iteration of Anderson acceleration. The filtering strategy consists of two steps: the first controls the length disparity between columns of the least-squares matrix, and the second enforces a lower bound on the angles between subspaces spanned by the columns of that matrix. The combined strategy is shown to control the condition number of the least-squares matrix at each iteration. The method is shown to be effective on a range of problems based on discretizations of partial differential equations. It is shown particularly effective for problems where the initial iterate may lie far from the solution, and which progress through distinct preasymptotic and asymptotic phases., Comment: 22 pages, 11 figures, 2 tables

Notwithstanding its long history, the study of nonlinear optics from metal surfaces is still an active field of research. For instance, in view of the presence of absorption questions remain concerning the possibility of significantly enhancing harmonic conversion efficiencies in the visible and UV ranges. While to many it may seem that metals do not easily lend themselves to that purpose, they are nevertheless crucial materials in the development of nanophotonics, and more generally, to electromagnetism at the nanoscale. Here, we report our experimental observations and numerical simulations of second and third harmonic generation from a gold nanograting, which exhibits a plasmonic resonance whose spectral position depends on incident angle. All things being equal, the enhancement of nonlinear optical processes from the UV to the near IR range manifests itself in dramatic manner: second harmonic generation conversion efficiencies increase more than three orders of magnitude compared to a flat gold mirror, while third harmonic generation conversion efficiency increases by nearly four orders of magnitude, both in excellent agreement with predictions. The clear inferences one may draw from our results are that our model describes the dynamics with unprecedented accuracy, and that much remains to be revealed in the development of nonlinear optics of metals at the nanoscale.

We consider two modifications of the Arrow-Hurwicz (AH) iteration for solving the incompressible steady Navier-Stokes equations for the purpose of accelerating the algorithm: grad-div stabilization, and Anderson acceleration. AH is a classical iteration for general saddle point linear systems and it was later extended to Navier-Stokes iterations in the 1970's which has recently come under study again. We apply recently developed ideas for grad-div stabilization and divergence-free finite element methods along with Anderson acceleration of fixed point iterations to AH in order to improve its convergence. Analytical and numerical results show that each of these methods improves AH convergence, but the combination of them yields an efficient and effective method that is competitive with more commonly used solvers.

We investigate both theoretically and numerically the consistency between the nonlinear discretization in full order models (FOMs) and reduced order models (ROMs) for incompressible flows. To this end, we consider two cases: (i) FOM-ROM consistency, i.e., when we use the same nonlinearity discretization in the FOM and ROM; and (ii) FOM-ROM inconsistency, i.e., when we use different nonlinearity discretizations in the FOM and ROM. Analytically, we prove that while the FOM-ROM consistency yields optimal error bounds, FOM-ROM inconsistency yields additional terms dependent on the FOM divergence error, which prevent the ROM from recovering the FOM as the number of modes increases. Computationally, we consider channel flow around a cylinder and Kelvin-Helmholtz instability, and show that FOM-ROM consistency yields significantly more accurate results than FOM-ROM inconsistency.

This paper studies a finite element discretization of the regularized Bingham equations that describe viscoplastic flow. An efficient nonlinear solver for the discrete model is then proposed and analyzed. The solver is based on Anderson acceleration (AA) applied to a Picard iteration, and we show accelerated convergence of the method by applying AA theory (recently developed by the authors) to the iteration, after showing sufficient smoothness properties of the associated fixed point operator. Numerical tests of spatial convergence are provided, as are results of the model for 2D and 3D driven cavity simulations. For each numerical test, the proposed nonlinear solver is also tested and shown to be very effective and robust with respect to the regularization parameter as it goes to zero.

In this paper, we propose, analyze, and test an efficient algorithm for computing ensemble average of incompressible magnetohydrodynamics (MHD) flows, where instances/members correspond to varying kinematic viscosity, magnetic diffusivity, body forces, and initial conditions. The algorithm is decoupled in Els\"asser variables and permits a shared coefficient matrix for all members at each time-step. Thus, the algorithm is much more computationally efficient than separately computing simulations for each member using usual MHD algorithms. We prove the proposed algorithm is unconditionally stable and convergent. Several numerical tests are given to support the predicted convergence rates. Finally, we test the proposed scheme and observe how the physical behavior changes as the coupling number increases in a lid-driven cavity problem with mean Reynolds number $Re\approx 15000$, and as the deviation of uncertainties in the initial and boundary conditions increases in a channel flow past a step problem., Comment: 22 pages, 6 figures

We propose a numerical approximation method for the Cahn-Hilliard equations that incorporates continuous data assimilation in order to achieve long time accuracy. The method uses a C$^0$ interior penalty spatial discretization of the fourth order Cahn-Hilliard equations, together with a backward Euler temporal discretization. We prove the method is long time stable and long time accurate, for arbitrarily inaccurate initial conditions, provided enough data measurements are incorporated into the simulation. Numerical experiments illustrate the effectiveness of the method on a benchmark test problem.

This paper considers an enhancement of the classical iterated penalty Picard (IPP) method for the incompressible Navier-Stokes equations, where we restrict our attention to $O(1)$ penalty parameter, and Anderson acceleration (AA) is used to significantly improve its convergence properties. After showing the fixed point operator associated with the IPP iteration is Lipschitz continuous and Lipschitz continuously (Frechet) differentiable, we apply a recently developed general theory for AA to conclude that IPP enhanced with AA improves its linear convergence rate by the gain factor associated with the underlying AA optimization problem. Results for several challenging numerical tests are given and show that IPP with penalty parameter 1 and enhanced with AA is a very effective solver.

We introduce a mimetic dual-field discretization which conserves mass, kinetic energy and helicity for three-dimensional incompressible Navier-Stokes equations. The discretization makes use of a conservative dual-field mixed weak formulation where two evolution equations of velocity are employed and dual representations of the solution are sought for each variable. A temporal discretization, which staggers the evolution equations and handles the nonlinearity such that the resulting discrete algebraic systems are linear and decoupled, is constructed. The spatial discretization is mimetic in the sense that the finite dimensional function spaces form a discrete de Rham complex. Conservation of mass, kinetic energy and helicity in the absence of dissipative terms is proven at the discrete level. Proper dissipation rates of kinetic energy and helicity in the viscous case is also proven. Numerical tests supporting the method are provided., Comment: 26 pages, 11 figures

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

This paper develops an efficient and robust solution technique for the steady Boussinesq model of non-isothermal flow using Anderson acceleration applied to a Picard iteration. After analyzing the fixed point operator associated with the nonlinear iteration to prove that certain stability and regularity properties hold, we apply the authors' recently constructed theory for Anderson acceleration, which yields a convergence result for the Anderson accelerated Picard iteration for the Boussinesq system. The result shows that the leading term in the residual is improved by the gain in the optimization problem, but at the cost of additional higher order terms that can be significant when the residual is large. We perform numerical tests that illustrate the theory, and show that a 2-stage choice of Anderson depth can be advantageous. We also consider Anderson acceleration applied to the Newton iteration for the Boussinesq equations, and observe that the acceleration allows the Newton iteration to converge for significantly higher Rayleigh numbers that it could without acceleration, even with a standard line search.