Search NASA⌕ Search

SEARCH · Search NASA

Results for “exponential convergence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Research highlights: June 1990 - May 1991

Linear instability calculations at MSFC have suggested that the Geophysical Fluid Flow Cell (GFFC) should exhibit classic baroclinic instability at accessible parameter settings. Interest was in the mechanisms of transition to temporal chaos and the evolution of spatio-temporal chaos. In order to understand more about such transitions, high resolution numerical experiments for the physically simplest model of two layer baroclinic instability were conducted. This model has the advantage that the numerical code is exponentially convergent and can be efficiently run for very long times, enabling the study of chaotic attractors without the often devastating effects of low-order trunction found in many previous studies. Numerical algorithms for implementing an empirical orthogonal function (EOF) analysis of the high resolution numerical results were completed. Under conditions of rapid rotation and relatively low differential heating, convection in a spherical shell takes place as columnar banana cells wrapped around the annular gap, but with axes oriented along the axis of rotation; these were clearly evident in the GFFC experiments. The results of recent numerical simulations of columnar convection and future research plans are presented.

Source record↗

Theory and implementation of high-order adaptive hp methods for analysis of incompressible viscous flows

An account is given of 'smart' algorithms for CFD which change in structure and performance with during flow calculations to accommodate changing properties of the solution. Such algorithms prominently include adaptive FEM methods, which are designed to adjust mesh parameters for the control of numerical error; attention is presently given to those which change the mesh size h and the local spectral order p in order to achieve high accuracies with minimal numbers of degrees of freedom. These 'hp methods' produce exponentially convergent approximations through which flow features are resolved by automatically distributing element sizes and spectral orders. This leads to calculation of the local (elementwise) error in the approximation.

Oden, J. T.↗

On the Gibbs phenomenon 3: Recovering exponential accuracy in a sub-interval from a spectral partial sum of a piecewise analytic function

The investigation of overcoming Gibbs phenomenon was continued, i.e., obtaining exponential accuracy at all points including at the discontinuities themselves, from the knowledge of a spectral partial sum of a discontinuous but piecewise analytic function. It was shown that if we are given the first N expansion coefficients of an L(sub 2) function f(x) in terms of either the trigonometrical polynomials or the Chebyshev or Legendre polynomials, an exponentially convergent approximation to the point values of f(x) in any sub-interval in which it is analytic can be constructed.

Gottlieb, David↗

An HP Adaptive Discontinuous Galerkin Method for Hyperbolic Conservation Laws

This dissertation addresses various issues for model classes of hyperbolic conservation laws. The basic approach developed in this work employs a new family of adaptive, hp-version, finite element methods based on a special discontinuous Galerkin formulation for hyperbolic problems. The discontinuous Galerkin formulation admits high-order local approximations on domains of quite general geometry, while providing a natural framework for finite element approximations and for theoretical developments. The use of hp-versions of the finite element method makes possible exponentially convergent schemes with very high accuracies in certain cases; the use of adaptive hp-schemes allows h-refinement in regions of low regularity and p-enrichment to deliver high accuracy, while keeping problem sizes manageable and dramatically smaller than many conventional approaches. The use of discontinuous Galerkin methods is uncommon in applications, but the methods rest on a reasonable mathematical basis for low-order cases and has local approximation features that can be exploited to produce very efficient schemes, especially in a parallel, multiprocessor environment. The place of this work is to first and primarily focus on a model class of linear hyperbolic conservation laws for which concrete mathematical results, methodologies, error estimates, convergence criteria, and parallel adaptive strategies can be developed, and to then briefly explore some extensions to more general cases. Next, we provide preliminaries to the study and a review of some aspects of the theory of hyperbolic conservation laws. We also provide a review of relevant literature on this subject and on the numerical analysis of these types of problems.

Bey, Kim S.↗

On the Gibbs phenomenon 4: Recovering exponential accuracy in a sub-interval from a Gegenbauer partial sum of a piecewise analytic function

We continue our investigation of overcoming Gibbs phenomenon, i.e., to obtain exponential accuracy at all points (including at the discontinuities themselves), from the knowledge of a spectral partial sum of a discontinuous but piecewise analytic function. We show that if we are given the first N Gegenbauer expansion coefficients, based on the Gegenbauer polynomials C(sub k)(sup mu)(x) with the weight function (1 - x(exp 2))(exp mu - 1/2) for any constant mu is greater than or equal to 0, of an L(sub 1) function f(x), we can construct an exponentially convergent approximation to the point values of f(x) in any subinterval in which the function is analytic. The proof covers the cases of Chebyshev or Legendre partial sums, which are most common in applications.

Gottlieb, David↗

On the Gibbs phenomenon 5: Recovering exponential accuracy from collocation point values of a piecewise analytic function

The paper presents a method to recover exponential accuracy at all points (including at the discontinuities themselves), from the knowledge of an approximation to the interpolation polynomial (or trigonometrical polynomial). We show that if we are given the collocation point values (or a highly accurate approximation) at the Gauss or Gauss-Lobatto points, we can reconstruct a uniform exponentially convergent approximation to the function f(x) in any sub-interval of analyticity. The proof covers the cases of Fourier, Chebyshev, Legendre, and more general Gegenbauer collocation methods.

Gottlieb, David↗

Solution of acoustic workshop problems by a spectral multidomain method

We use a new staggered grid Chebyshev spectral multidomain method to solve three of the Workshop benchmark problems. The method defines solution unknowns at the nodes of the Chebyshev Gauss quadrature, and the fluxes at the nodes of the Chebyshev Gauss-Lobatto quadrature. The Chebyshev spectral method gives exponentially convergent phase and dissipation errors. The multidomain approximation gives the method flexibility. Using the method, we solve problems in Categories 1 and 5 of the benchmark problems.

Kopriva, Davis A.↗

Magnetic Field Suppression of Flow in Semiconductor Melt

One of the most promising approaches for the reduction of convection during the crystal growth of conductive melts (semiconductor crystals) is the application of magnetic fields. Current technology allows the experimentation with very intense static fields (up to 80 KGauss) for which nearly convection free results are expected from simple scaling analysis in stabilized systems (vertical Bridgman method with axial magnetic field). However, controversial experimental results were obtained. The computational methods are, therefore, a fundamental tool in the understanding of the phenomena accounting during the solidification of semiconductor materials. Moreover, effects like the bending of the isomagnetic lines, different aspect ratios and misalignments between the direction of the gravity and magnetic field vectors can not be analyzed with analytical methods. The earliest numerical results showed controversial conclusions and are not able to explain the experimental results. Although the generated flows are extremely low, the computational task is a complicated because of the thin boundary layers. That is one of the reasons for the discrepancy in the results that numerical studies reported. Modeling of these magnetically damped crystal growth experiments requires advanced numerical methods. We used, for comparison, three different approaches to obtain the solution of the problem of thermal convection flows: (1) Spectral method in spectral superelement implementation, (2) Finite element method with regularization for boundary layers, (3) Multiquadric method, a novel method with global radial basis functions, that is proven to have exponential convergence. The results obtained by these three methods are presented for a wide region of Rayleigh and Hartman numbers. Comparison and discussion of accuracy, efficiency, reliability and agreement with experimental results will be presented as well.

Fedoseyev, A. I.↗

A Coupled Nonlinear Spacecraft Attitude Controller/Observer With an Unknown Constant Gyro Bias

A nonlinear control scheme for attitude control of a spacecraft is combined with a nonlinear gyro bias observer for the case of constant gyro bias. The closed loop system is proven to be globally stable, with zero tracking error, thus proving a separation principle for the given system. The nonlinear observer incorporates persistency of excitation, resulting in exponential convergence of the gyro bias error.

Deutschmann, Julie↗

A Non-Linear Approach to Spacecraft Trajectory Control in the Vicinity of a Libration Point

An expanding interest in mission design strategies that exploit libration point regions demands the continued development of enhanced, efficient, control algorithms for station-keeping and formation maintenance. This paper discusses the development of a non-linear, station-keeping, control algorithm for trajectories in the vicinity of a libration point. The control law guarantees exponential convergence, based on a Lyaponov analysis. Controller performance is evaluated using FreeFlyer(R) and MATLAB(R) for a spacecraft stationed near the L2 libration point in the Earth-Moon system, tracking a pre-defined reference trajectory. Evaluation metrics are fuel usage and tracking accuracy. Simulation results are compared with a linear-based controller for a spacecraft tracking the same reference trajectory. Although the analysis is framed in the context of station-keeping, the control algorithm is equally applicable to a formation flying problem with an appropriate definition of the reference trajectory.

Luquette, Richard J.↗

A Nonlinear Spacecraft Attitude Controller and Observer with an Unknown Constant Gyro Bias and Gyro Noise

A nonlinear control scheme for attitude control of a spacecraft is combined with a nonlinear gyro bias observer for the case of constant gyro bias, in the presence of gyro noise. The observer bias estimates converge exponentially to a mean square bound determined by the standard deviation of the gyro noise. The resulting coupled, closed loop dynamics are proven to be globally stable, with asymptotic tracking which is also mean square bounded. A simulation of the proposed observer-controller design is given for a rigid spacecraft tracking a specified, time-varying attitude sequence to illustrate the theoretical claims.

Deutschmann, Julie↗

A Non-Linear Approach to Spacecraft Formation Control in the Vicinity of a Collinear Libration Point

An expanding interest in mission design strategies that exploit libration point regions demands the continued development of enhanced, efficient, control algorithms for station-keeping and formation maintenance. This paper discusses the development of a non-linear, formation maintenance, control algorithm for trajectories in the vicinity of a libration point. However, the formulation holds for any trajectory governed by the equations of motion for the restricted three body problem. The control law guarantees exponential convergence, based on a Lyaponov analysis. FreeFlyer and MATLAB provide the simulation environment for controller performance evaluation. The simulation, modeled after the MAXIM Pathfinder mission, maintains the relative position of a "follower" spacecraft with respect to a "leader" spacecraft, stationed near the L2 libration point in the Sun-Earth system. Evaluation metrics are fuel usage and tracking accuracy.

Luquette, Richard J.↗

A Nonlinear Observer for Gyro Alignment Estimation

A nonlinear observer for gyro alignment estimation is presented. The observer is composed of two error terms, an attitude error and an alignment error. The observer is globally stable with exponential convergence of the attitude errors. The gyro alignment estimate converges to the true alignment when the system is completely observable.

Thienel, J.↗

Observers for Systems with Nonlinearities Satisfying an Incremental Quadratic Inequality

We consider the problem of state estimation for nonlinear time-varying systems whose nonlinearities satisfy an incremental quadratic inequality. These observer results unifies earlier results in the literature; and extend it to some additional classes of nonlinearities. Observers are presented which guarantee that the state estimation error exponentially converges to zero. Observer design involves solving linear matrix inequalities for the observer gain matrices. Results are illustrated by application to a simple model of an underwater.

observers↗

Joint Chance-Constrained Dynamic Programming

This paper presents a novel dynamic programming algorithm with a joint chance constraint, which explicitly bounds the risk of failure in order to maintain the state within a specified feasible region. A joint chance constraint cannot be handled by existing constrained dynamic programming approaches since their application is limited to constraints in the same form as the cost function, that is, an expectation over a sum of one-stage costs. We overcome this challenge by reformulating the joint chance constraint into a constraint on an expectation over a sum of indicator functions, which can be incorporated into the cost function by dualizing the optimization problem. As a result, the primal variables can be optimized by a standard dynamic programming, while the dual variable is optimized by a root-finding algorithm that converges exponentially. Error bounds on the primal and dual objective values are rigorously derived. We demonstrate the algorithm on a path planning problem, as well as an optimal control problem for Mars entry, descent and landing. The simulations are conducted using a real terrain data of Mars, with four million discrete states at each time step.

Ono, Masahiro↗

Simultaneous Trajectory Tracking Control and Online Mass Estimation for a Regolith Excavating Robot via Integral Concurrent Learning

An adaptive controller is developed for a regolith excavation robot to determine the mass of excavated material and to account for the effects of gravity and friction while on the surface of other celestial bodies. A data-based integral concurrent learning (ICL) parameter update law accounts for and estimates the unknown mass, gravity, and friction parameters. A Lyapunov-based analysis proves that the trajectory tracking error and the parameter estimate errors exponentially converge to zero. An estimation of the mass of regolith excavated by the robot is calculated from the estimated parameters. A simulation study is performed to show the performance of the developed technique. Simulation results show that for 3 kilograms of excavated material, the mass estimate has an error of 4.9 grams.

Max L. Greene↗

Linear Algorithms for Worst Case Identification in H infinity with Applications to Flexible Structures

This paper is concerned with linear algorithms for identification in which have been studied in [9]. It is shown that the two different linear algorithms in [9] can be unified into a single one which can be further extended to nonuniformly spaced frequency response samples with exponential convergence for the noise free case. Improved upper bounds for the corresponding identification errors are derived. Applications to the identification of lightly damped systems such as flexible structures are also considered.

linear↗