2,177 research outputs found
An Ensemble-Proper Orthogonal Decomposition Method for the Nonstationary Navier-Stokes Equations
The definition of partial differential equation (PDE) models usually involves
a set of parameters whose values may vary over a wide range. The solution of
even a single set of parameter values may be quite expensive. In many cases,
e.g., optimization, control, uncertainty quantification, and other settings,
solutions are needed for many sets of parameter values. We consider the case of
the time-dependent Navier-Stokes equations for which a recently developed
ensemble-based method allows for the efficient determination of the multiple
solutions corresponding to many parameter sets. The method uses the average of
the multiple solutions at any time step to define a linear set of equations
that determines the solutions at the next time step. To significantly further
reduce the costs of determining multiple solutions of the Navier-Stokes
equations, we incorporate a proper orthogonal decomposition (POD) reduced-order
model into the ensemble-based method. The stability and convergence results for
the ensemble-based method are extended to the ensemble-POD approach. Numerical
experiments are provided that illustrate the accuracy and efficiency of
computations determined using the new approach
The instanton method and its numerical implementation in fluid mechanics
A precise characterization of structures occurring in turbulent fluid flows
at high Reynolds numbers is one of the last open problems of classical physics.
In this review we discuss recent developments related to the application of
instanton methods to turbulence. Instantons are saddle point configurations of
the underlying path integrals. They are equivalent to minimizers of the related
Freidlin-Wentzell action and known to be able to characterize rare events in
such systems. While there is an impressive body of work concerning their
analytical description, this review focuses on the question on how to compute
these minimizers numerically. In a short introduction we present the relevant
mathematical and physical background before we discuss the stochastic Burgers
equation in detail. We present algorithms to compute instantons numerically by
an efficient solution of the corresponding Euler-Lagrange equations. A second
focus is the discussion of a recently developed numerical filtering technique
that allows to extract instantons from direct numerical simulations. In the
following we present modifications of the algorithms to make them efficient
when applied to two- or three-dimensional fluid dynamical problems. We
illustrate these ideas using the two-dimensional Burgers equation and the
three-dimensional Navier-Stokes equations
Robust and efficient solution of the drum problem via Nystrom approximation of the Fredholm determinant
The drum problem-finding the eigenvalues and eigenfunctions of the Laplacian
with Dirichlet boundary condition-has many applications, yet remains
challenging for general domains when high accuracy or high frequency is needed.
Boundary integral equations are appealing for large-scale problems, yet certain
difficulties have limited their use. We introduce two ideas to remedy this: 1)
We solve the resulting nonlinear eigenvalue problem using Boyd's method for
analytic root-finding applied to the Fredholm determinant. We show that this is
many times faster than the usual iterative minimization of a singular value. 2)
We fix the problem of spurious exterior resonances via a combined field
representation. This also provides the first robust boundary integral
eigenvalue method for non-simply-connected domains. We implement the new method
in two dimensions using spectrally accurate Nystrom product quadrature. We
prove exponential convergence of the determinant at roots for domains with
analytic boundary. We demonstrate 13-digit accuracy, and improved efficiency,
in a variety of domain shapes including ones with strong exterior resonances.Comment: 21 pages, 7 figures, submitted to SIAM Journal of Numerical Analysis.
Updated a duplicated picture. All results unchange
- …