Search NASASearch

Engineering topics

Patera, A. T.

Publications and source records attributed to Patera, A. T..

Reliable Real-Time Solution of Parametrized Partial Differential Equations: Reduced-Basis Output Bound Methods

We present a technique for the rapid and reliable prediction of linear-functional outputs of elliptic (and parabolic) partial differential equations with affine parameter dependence. The essential components are (i) (provably) rapidly convergent global reduced basis approximations, Galerkin projection onto a space W(sub N) spanned by solutions of the governing partial differential equation at N selected points in parameter space; (ii) a posteriori error estimation, relaxations of the error-residual equation that provide inexpensive yet sharp and rigorous bounds for the error in the outputs of interest; and (iii) off-line/on-line computational procedures, methods which decouple the generation and projection stages of the approximation process. The operation count for the on-line stage, in which, given a new parameter value, we calculate the output of interest and associated error bound, depends only on N (typically very small) and the parametric complexity of the problem; the method is thus ideally suited for the repeated and rapid evaluations required in the context of parameter estimation, design, optimization, and real-time control.

Prudhomme, C.

Spectral collocation methods

This review covers the theory and application of spectral collocation methods. Section 1 describes the fundamentals, and summarizes results pertaining to spectral approximations of functions. Some stability and convergence results are presented for simple elliptic, parabolic, and hyperbolic equations. Applications of these methods to fluid dynamics problems are discussed in Section 2.

Hussaini, M. Y.

Spectral collocation methods

This review covers the theory and application of spectral collocation methods. Section 1 describes the fundamentals, and summarizes results pertaining to spectral approximations of functions. Some stability and convergence results are presented for simple elliptic, parabolic, and hyperbolic equations. Applications of these methods to fluid dynamics problems are discussed in Section 2.

Hussaini, M. Y.

A well-posed optimal spectral element approximation for the Stokes problem

A method is proposed for the spectral element simulation of incompressible flow. This method constitutes in a well-posed optimal approximation of the steady Stokes problem with no spurious modes in the pressure. The resulting method is analyzed, and numerical results are presented for a model problem.

Maday, Y.

An isoparametric spectral element method for solution of the Navier-Stokes equations in complex geometry

High-order (p-type) finite element methods combine the geometric flexibility of standard low-order finite element (or finite volume) techniques with the rapid convergence properties of spectral methods. Various p-type schemes have been proposed for elliptic problems. However, an application of p-type methods to more complex equations, in particular the passive scalar or Navier-Stokes equations, is complicated by the nature of the equations and by efficiency considerations as regards the solution procedure. Patera (1984) has proposed a p-type method for the Navier-Stokes equations, taking into account the arising additional complications. The present study is concerned with an isoparametric spectral element discretization for the solution of flow problems in fully general, curvy geometries. Attention is given to the isoparametric formulation for elliptic equations and a time-splitting scheme for the Navier-Stokes equations.

Korczak, K. Z.

Advances and future directions of research on spectral methods

Recent advances in spectral methods are briefly reviewed and characterized with respect to their convergence and computational complexity. Classical finite element and spectral approaches are then compared, and spectral element (or p-type finite element) approximations are introduced. The method is applied to the full Navier-Stokes equations, and examples are given of the application of the technique to several transitional flows. Future directions of research in the field are outlined.

Patera, A. T.

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.

Secondary instability of wall-bounded shear flows

The present analysis of a secondary instability in a wide class of wall-bounded parallel shear flows indicates that two-dimensional, finite amplitude waves are exponentially unstable to infinitessimal three-dimensional disturbances. The instability appears to be the prototype of transitional instability in such flows as Poiseuille flow, Couette flow, and flat plate boundary layers, in that it has the convective time scales observed in the typical transitions. The energetics and vorticity dynamics of the instability are discussed, and it is shown that the two-dimensional perturbation without directly providing energy to the disturbance. The three-dimensional instability requires that a threshold two-dimensional amplitude be achieved. It is found possible to identify experimental features of transitional spot structure with aspects of the nonlinear two-dimensional/linear three-dimensional instability.

Orszag, S. A.

Spectral methods for flows in complex geometries

Three methods for the application of spectral methods to flows in complex geometries are discussed. They are the spectral iteration method, the spectral embedding method, and the spectral element method. Applications are given to the flow over wavy walls, past a large-eddy breakup device (LEBU), and to flow over a step.

Orszag, S. A.

Three-dimensional instability of plane channel flows at subcritical Reynolds numbers

A three-dimensional instability that predicts subcritical transition in plane channel flows in good agreement with experiment is discussed. Both full simulations of the Navier-Stokes equations and a three-dimensional linear analysis about two-dimensional finite amplitude states confirm that two-dimensional secondary flows are strongly unstable to small three-dimensional perturbations.

Orszag, S. A.

Subcritical transition to turbulence in planar shear flows

The two-dimensional steady and time dependent properties of plane Poiseuille and plane Couette flows are analyzed using iterative techniques and full numerical simulation of the Navier-Stokes equations. It is shown that the finite-amplitude two-dimensional states investigated are strongly unstable to very small three-dimensional perturbations. It is also shown, through full numerical simulation, that this explosive secondary instability can explain the subcritical transitions that occur in real flows. Finally, it is shown that the three-dimensional instability can be analyzed by a linear stability analysis of a two-dimensional flow consisting of the basic parallel flow and a steady (or quasi-steady) finite-amplitude two-dimensional cellular motion.

Orszag, S. A.

Subcritical transition to turbulence in plane channel flows

A linear three dimensional mechanism for the transition of plane Poiseuille flows to turbulence is presented which provides good agreement with experimental observations. The mechanism is based on the evolution of states within a band of quasi-equilibria which slowly approach the stable upper branch solutions for the evolution of flow energy but which are strongly unstable to infinitesimal three-dimensional disturbances. Numerical simulation has shown that if two-dimensional flow persists long enough for the three-dimensional perturbations to attain finite amplitude, the resulting three dimensional flow quickly develops a turbulent character with nonperiodic behavior, and thus transition can be predicted from knowledge of the initial two- and three-dimensional energies and time scales. The mechanism predicts transition to turbulence at Reynolds numbers greater than 1000, as observed in experiments, and implies higher threshold three-dimensional energies in plane Couette flow.

Orszag, S. A.