skip to main content

Title: A time accurate, adaptive discretization for fluid flow problems
This report presents a low computational and cognitive complexity, stable, time accurate and adaptive method for the Navier-Stokes equations. The improved method requires a minimally intrusive modification to an existing program based on the fully implicit / backward Euler time discretization, does not add to the computational complexity, and is conceptually simple. The backward Euler approximation is simply post-processed with a two-step, linear time filter. The time filter additionally removes the overdamping of Backward Euler while remaining unconditionally energy stable, proven herein. Even for constant stepsizes, the method does not reduce to a standard / named time stepping method but is related to a known 2-parameter family of A-stable, two step, second order methods. Numerical tests confirm the predicted convergence rates and the improved predictions of flow quantities such as drag and lift.
; ;
Award ID(s):
Publication Date:
Journal Name:
International journal of numerical analysis and modeling
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Adaptive time stepping methods for metastable dynamics of the Allen–Cahn and Cahn–Hilliard equations are investigated in the spatially continuous, semi-discrete setting. We analyse the performance of a number of first and second order methods, formally predicting step sizes required to satisfy specified local truncation error $$\sigma $$ σ in the limit of small length scale parameter $$\epsilon \rightarrow 0$$ ϵ → 0 during meta-stable dynamics. The formal predictions are made under stability assumptions that include the preservation of the asymptotic structure of the diffuse interface, a concept we call profile fidelity. In this setting, definite statements about the relative behaviour of time stepping methods can be made. Some methods, including all so-called energy stable methods but also some fully implicit methods, require asymptotically more time steps than others. The formal analysis is confirmed in computational studies. We observe that some provably energy stable methods popular in the literature perform worse than some more standard schemes. We show further that when Backward Euler is applied to meta-stable Allen–Cahn dynamics, the energy decay and profile fidelity properties for these discretizations are preserved for much larger time steps than previous analysis would suggest. The results are established asymptotically for general interfaces, withmore »a rigorous proof for radial interfaces. It is shown analytically and computationally that for most reaction terms, Eyre type time stepping performs asymptotically worse due to loss of profile fidelity.« less
  2. Abstract We present a discontinuity aware quadrature (DAQ) rule and use it to develop implicit self-adaptive theta (SATh) schemes for the approximation of scalar hyperbolic conservation laws. Our SATh schemes require the solution of a system of two equations, one controlling the cell averages of the solution at the time levels, and the other controlling the space-time averages of the solution. These quantities are used within the DAQ rule to approximate the time integral of the hyperbolic flux function accurately, even when the solution may be discontinuous somewhere over the time interval. The result is a finite volume scheme using the theta time stepping method, with theta defined implicitly (or self-adaptively). Two schemes are developed, self-adaptive theta upstream weighted (SATh-up) for a monotone flux function using simple upstream stabilization, and self-adaptive theta Lax–Friedrichs (SATh-LF) using the Lax–Friedrichs numerical flux. We prove that DAQ is accurate to second order when there is a discontinuity in the solution and third order when it is smooth. We prove that SATh-up is unconditionally stable, provided that theta is set to be at least 1/2 (which means that SATh can be only first order accurate in general). We also prove that SATh-up satisfies the maximummore »principle, and is total variation diminishing under appropriate monotonicity and boundary conditions. General flux functions require the SATh-LF scheme, so we assess its accuracy through numerical examples in one and two space dimensions. These results suggest that SATh-LF is also stable and satisfies the maximum principle (at least at reasonable Courant-Friedrichs-Lewy numbers). Compared to the solutions of finite volume schemes using Crank–Nicolson and backward Euler time stepping, SATh-LF solutions often approach the accuracy of the former, but without oscillation, and they are numerically less diffuse than the latter.« less
  3. This paper develops, analyzes and tests a time-accurate partitioned method for the Stokes-Darcy equations. The method combines a time filter and Backward Euler scheme, is second order accurate and provide, at no extra complexity, an estimated the temporal error. This approach post-processes the solutions of Backward Euler scheme by adding three lines to original codes to increase the time accuracy from first order to second order. We prove long time stability and error estimates of Backward Euler plus time filter with constant time stepsize. Moreover, we extend the approach to variable time stepsize and construct adaptive algorithms. Numerical tests show convergence of our method and support the theoretical analysis.
  4. This paper is concerned with the optimized Schwarz waveform relaxation method and Ventcel transmission conditions for the linear advection-diffusion equation. A mixed formulation is considered in which the flux variable represents both diffusive and advective flux, and Lagrange multipliers are introduced on the interfaces between nonoverlapping subdomains to handle tangential derivatives in the Ventcel conditions. A space-time interface problem is formulated and is solved iteratively. Each iteration involves the solution of time-dependent problems with Ventcel boundary conditions in the subdomains. The subdomain problems are discretized in space by a mixed hybrid finite element method based on the lowest-order Raviart-Thomas space and in time by the backward Euler method. The proposed algorithm is fully implicit and enables different time steps in the subdomains. Numerical results with discontinuous coefficients and various Peclét numbers validate the accuracy of the method with nonconforming time grids and confirm the improved convergence properties of Ventcel conditions over Robin conditions.

  5. Abstract

    Our aim is to approximate a reference velocity field solving the two-dimensional Navier–Stokes equations (NSE) in the absence of its initial condition by utilizing spatially discrete measurements of that field, available at a coarse scale, and continuous in time. The approximation is obtained via numerically discretizing a downscaling data assimilation algorithm. Time discretization is based on semiimplicit and fully implicit Euler schemes, while spatial discretization (which can be done at an arbitrary scale regardless of the spatial resolution of the measurements) is based on a spectral Galerkin method. The two fully discrete algorithms are shown to be unconditionally stable, with respect to the size of the time step, the number of time steps and the number of Galerkin modes. Moreover, explicit, uniform-in-time error estimates between the approximation and the reference solution are obtained, in both the $L^2$ and $H^1$ norms. Notably, the two-dimensional NSE, subject to the no-slip Dirichlet or periodic boundary conditions, are used in this work as a paradigm. The complete analysis that is presented here can be extended to other two- and three-dimensional dissipative systems under the assumption of global existence and uniqueness.