468 research outputs found
Computability of differential equations
In this chapter, we provide a survey of results concerning the computability and computational complexity of differential equations. In particular, we study the conditions which ensure computability of the solution to an initial value problem for an ordinary differential equation (ODE) and analyze the computational complexity of a computable solution. We also present computability results concerning the asymptotic behaviors of ODEs as well as several classically important partial differential equations.info:eu-repo/semantics/acceptedVersio
Three regularization models of the Navier-Stokes equations
We determine how the differences in the treatment of the subfilter-scale
physics affect the properties of the flow for three closely related
regularizations of Navier-Stokes. The consequences on the applicability of the
regularizations as SGS models are also shown by examining their effects on
superfilter-scale properties. Numerical solutions of the Clark-alpha model are
compared to two previously employed regularizations, LANS-alpha and Leray-alpha
(at Re ~ 3300, Taylor Re ~ 790) and to a DNS. We derive the Karman-Howarth
equation for both the Clark-alpha and Leray-alpha models. We confirm one of two
possible scalings resulting from this equation for Clark as well as its
associated k^(-1) energy spectrum. At sub-filter scales, Clark-alpha possesses
similar total dissipation and characteristic time to reach a statistical
turbulent steady-state as Navier-Stokes, but exhibits greater intermittency. As
a SGS model, Clark reproduces the energy spectrum and intermittency properties
of the DNS. For the Leray model, increasing the filter width decreases the
nonlinearity and the effective Re is substantially decreased. Even for the
smallest value of alpha studied, Leray-alpha was inadequate as a SGS model. The
LANS energy spectrum k^1, consistent with its so-called "rigid bodies,"
precludes a reproduction of the large-scale energy spectrum of the DNS at high
Re while achieving a large reduction in resolution. However, that this same
feature reduces its intermittency compared to Clark-alpha (which shares a
similar Karman-Howarth equation). Clark is found to be the best approximation
for reproducing the total dissipation rate and the energy spectrum at scales
larger than alpha, whereas high-order intermittency properties for larger
values of alpha are best reproduced by LANS-alpha.Comment: 21 pages, 8 figure
Increased Accuracy and Efficiency in Finite Element Computations of the Leray-Deconvolution Model of Turbulence
This thesis develops, analyzes and tests a finite element method for approximating solutions to the Lerayâdeconvolution regularization of the NavierâStokes equations. The scheme combines three ideas in order to create an accurate and effective algorithm: the use of an incompressible filter, a linearization that decouples the velocityâpressure system from the filtering and deconvolution operations, and a stabilization that works well with the linearization. A rigorous and complete numerical analysis of the scheme is given, and numerical experiments are presented that show clear advantages of the scheme
The GNAT method for nonlinear model reduction: effective implementation and application to computational fluid dynamics and turbulent flows
The Gauss--Newton with approximated tensors (GNAT) method is a nonlinear
model reduction method that operates on fully discretized computational models.
It achieves dimension reduction by a Petrov--Galerkin projection associated
with residual minimization; it delivers computational efficency by a
hyper-reduction procedure based on the `gappy POD' technique. Originally
presented in Ref. [1], where it was applied to implicit nonlinear
structural-dynamics models, this method is further developed here and applied
to the solution of a benchmark turbulent viscous flow problem. To begin, this
paper develops global state-space error bounds that justify the method's design
and highlight its advantages in terms of minimizing components of these error
bounds. Next, the paper introduces a `sample mesh' concept that enables a
distributed, computationally efficient implementation of the GNAT method in
finite-volume-based computational-fluid-dynamics (CFD) codes. The suitability
of GNAT for parameterized problems is highlighted with the solution of an
academic problem featuring moving discontinuities. Finally, the capability of
this method to reduce by orders of magnitude the core-hours required for
large-scale CFD computations, while preserving accuracy, is demonstrated with
the simulation of turbulent flow over the Ahmed body. For an instance of this
benchmark problem with over 17 million degrees of freedom, GNAT outperforms
several other nonlinear model-reduction methods, reduces the required
computational resources by more than two orders of magnitude, and delivers a
solution that differs by less than 1% from its high-dimensional counterpart
Output error estimation strategies for discontinuous Galerkin discretizations of unsteady convectionâdominated flows
We study practical strategies for estimating numerical errors in scalar outputs calculated from unsteady simulations of convectionâdominated flows, including those governed by the compressible NavierâStokes equations. The discretization is a discontinuous Galerkin finite element method in space and time on static spatial meshes. Timeâintegral quantities are considered for scalar outputs and these are shown to superconverge with temporal refinement. Output error estimates are calculated using the adjointâweighted residual method, where the unsteady adjoint solution is obtained using a discrete approach with an iterative solver. We investigate the accuracy versus computational cost tradeâoff for various approximations of the fineâspace adjoint and find that exact adjoint solutions are accurate but expensive. To reduce the cost, we propose a local temporal reconstruction that takes advantage of superconvergence properties at Radau points, and a spatial reconstruction based on nearestâneighbor elements. This inexact adjoint yields output error estimates at a computational cost of less than 2.5 times that of the forward problem for the cases tested. The calculated error estimates account for numerical error arising from both the spatial and temporal discretizations, and we present a method for identifying the percentage contributions of each discretization to the output error. Copyright © 2011 John Wiley & Sons, Ltd.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/88080/1/3224_ftp.pd
- âŠ