Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Discretization error”

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 415 records · Page 23

Time maintenance of user clocks via the tracking and data relay satellite system

A system is described which uses the Tracking and Data Relay Satellite System (TDRSS) itself to compare the user satellite clock with a clock at the White Sands station that is referenced to Universal Time Coordinated (UTC). No command of the spacecraft by the system is required, and actual on-board clock corrections are made by the spacecraft control center at its discretion. Computer models were constructed using basic orbital parameters for user and TDRS satellites. With only first-order corrections and simple averaging techniques for constant clock rates, error measurement precision of better than one microsecond was obtained. More sophisticated computations should allow considerable improvement over this.

Whitworth, G.↗

Parametric conditions for stability of reduced-order linear time-varying control systems

Using a single framework, parametric conditions are derived which encompass those for both local and global BIBO stability of a linear multivariable discrete-time reduced-order time-varying control system. These conditions indicate that the system will be BIBO stable if the norm of the system-parameter error matrix is bounded by an l exp 1 function superimposed on an l exp infinity function.

Ma, C. C. H.↗

Establishing Best Practices for X-57 Maxwell CFD Database Generation

The X-57 Maxwell is NASA’s latest electric airplane concept that has been simulated for aerodynamic performance using the structured overset and unstructured grid solvers within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework as well as the unstructured polyhedral grid solver in Star-CCM+ for code-to-code comparison. In order to validate the predictions, comparisons were made between the CFD solutions and experimental data collected in the 12-foot Low-Speed Wind Tunnel at NASA Langley Research Center. The simulations are in preparation for the development of a comprehensive aerodynamic database which will assess aircraft performance at a variety of conditions. The findings from these simulations will establish the best practices for mesh resolution, numerical discretization, and turbulence modeling to be used for this database. Preliminary database results have shown that best-practices learned from the initial validation simulations will potentially reduce error in X-57 aerodynamic loads and moments relative to experiment by up to 14%.

Duensing, Jared C.↗

Numerical comparison of discrete Kalman filter algorithms - Orbit determination case study

Numerical characteristics of various Kalman filter algorithms are illustrated with a realistic orbit determination study. The case study of this paper highlights the numerical deficiencies of the conventional and stabilized Kalman algorithms. Computational errors associated with these algorithms are found to be so large as to obscure important mismodeling effects and thus cause misleading estimates of filter accuracy. The positive result of this study is that the U-D covariance factorization algorithm has excellent numerical properties and is computationally efficient, having CPU costs that differ negligibly from the conventional Kalman costs. Accuracies of the U-D filter using single precision arithmetic consistently match the double precision reference results. Numerical stability of the U-D filter is further demonstrated by its insensitivity to variations in the a priori statistics.

Bierman, G. J.↗

Domain decomposition methods for systems of conservation laws: Spectral collocation approximations

Hyperbolic systems of conversation laws are considered which are discretized in space by spectral collocation methods and advanced in time by finite difference schemes. At any time-level a domain deposition method based on an iteration by subdomain procedure was introduced yielding at each step a sequence of independent subproblems (one for each subdomain) that can be solved simultaneously. The method is set for a general nonlinear problem in several space variables. The convergence analysis, however, is carried out only for a linear one-dimensional system with continuous solutions. A precise form of the error reduction factor at each iteration is derived. Although the method is applied here to the case of spectral collocation approximation only, the idea is fairly general and can be used in a different context as well. For instance, its application to space discretization by finite differences is straight forward.

Quarteroni, Alfio↗

Optimal filtering in the presence of unmodeled time correlated driving disturbances.

In many realistic data filtering problems, the cross correlation of the state estimation error and the state forcing function is unknown due to the poor knowledge of the time history of the forcing function. In this paper, the conservative and minimal approximation to the cross correlation terms is presented. It requires only the knowledge of the estimation error covariance and the forcing function covariance, with the choice of an associated free parameter left to the user. If the estimation error covariance and/or the forcing function covariance are bounded from above but not known exactly, the cross correlation approximation using those upper bounds remains conservative. This cross correlation approximation leads to a conservative approximation to the estimation error covariance matrix differential equation between measurement times. The free parameter is determined as the analytic solution to an associated optimal control problem. The procedure is expanded to include discrete linear measurement incorporation.

Fraser, D. C.↗

Aircraft landing response in a discrete multipath environment

This paper considers the problem of discrete multipath reflections upon an aircraft in the landing phase. A model is developed for the communication channel for a typical receiver. Simulation studies are presented showing the effects of discrete multipath upon the aircraft's ability to follow a specified flight path. A development is presented for the analytical determination of the probability density function of the angular errors.

Guarino, C. R.↗

Speech versus manual control of camera functions during a telerobotic task

Voice input for control of camera functions was investigated in this study. Objective were to (1) assess the feasibility of a voice-commanded camera control system, and (2) identify factors that differ between voice and manual control of camera functions. Subjects participated in a remote manipulation task that required extensive camera-aided viewing. Each subject was exposed to two conditions, voice and manual input, with a counterbalanced administration order. Voice input was found to be significantly slower than manual input for this task. However, in terms of remote manipulator performance errors and subject preference, there was no difference between modalities. Voice control of continuous camera functions is not recommended. It is believed that the use of voice input for discrete functions, such as multiplexing or camera switching, could aid performance. Hybrid mixes of voice and manual input may provide the best use of both modalities. This report contributes to a better understanding of the issues that affect the design of an efficient human/telerobot interface.

Bierschwale, John M.↗

Stiffness, thermal expansion, and thermal bending formulation of stiffened, fiber-reinforced composite panels

A method is presented for formulating stiffness terms and thermal coefficients of stiffened, fiber-reinforced composite panels. The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. Nonlinear, temperature and load dependent constitutive material data of each laminate are used to 'build-up' the stiffened panel membrane, bending, and membrane-bending coupling stiffness terms and thermal coefficients. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN analysis program (FEA) is revealed that provides accurate solutions of entire airframes or engines with coarsely meshed models. An example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. The erroneous results and the correct ones produced from this method are compared to analysis from discretely meshed three-dimensional FEA.

Collier, Craig S.↗

A Discontinuous Galerkin Method for Parabolic Problems with Modified hp-Finite Element Approximation Technique

A recent paper is generalized to a case where the spatial region is taken in R(sup 3). The region is assumed to be a thin body, such as a panel on the wing or fuselage of an aerospace vehicle. The traditional h- as well as hp-finite element methods are applied to the surface defined in the x - y variables, while, through the thickness, the technique of the p-element is employed. Time and spatial discretization scheme based upon an assumption of certain weak singularity of double vertical line u(sub t) double vertical line 2, is used to derive an optimal a priori error estimate for the current method.

Kaneko, Hideaki↗

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

Performance of a Space-Based Wavelet Compressor for Plasma Count Data on the MMS Fast Plasma Investigation

Data compression has been a staple of imaging instruments for years. Recently, plasma measurements have utilized compression with relatively low compression ratios. The Fast Plasma Investigation (FPI) on board the Magnetospheric Multiscale (MMS) mission generates data roughly 100 times faster than previous plasma instruments, requiring a higher compression ratio to fit within the telemetry allocation. This study investigates the performance of a space-based compression standard employing a Discrete Wavelet Transform and a Bit Plane Encoder (DWT/BPE) in compressing FPI plasma count data. Data from the first 6 months of FPI operation are analyzed to explore the error modes evident in the data and how to adapt to them. While approximately half of the Dual Electron Spectrometer (DES) maps had some level of loss, it was found that there is little effect on the plasma moments and that errors present in individual sky maps are typically minor. The majority of Dual Ion Spectrometer burst sky maps compressed in a lossless fashion, with no error introduced during compression. Because of induced compression error, the size limit for DES burst images has been increased for Phase 1B. Additionally, it was found that the floating point compression mode yielded better results when images have significant compression error, leading to floating point mode being used for the fast survey mode of operation for Phase 1B. Despite the suggested tweaks, it was found that wavelet-based compression, and a DWT/BPE algorithm in particular, is highly suitable to data compression for plasma measurement instruments and can be recommended for future missions.

plasma↗

Estimation With Range Depended Sensor Model

This paper focuses on the improvement of object detection accuracy taking into account the sensor’s reading degradation as the relative range increases. The approach is based on the assumption that range measurement error depends on the actual range. Specifically, we model the measurement error as a proportional to the actual range term plus a zero-mean, Gaussian distributed, and uncorrelated process. The tracking problem is considered in a mixed continuous-discrete time domain, where the target dynamics is in continuous-time and the measurements are in discrete-time, which is an optimal choice in many tracking and navigation applications. We adopt a commonly used continuous time coordinate-uncoupled white-noise acceleration model for a point object to describe the target motion, and use Extended Kalman Filter (EKF) framework to estimate the target's state and the unknown proportionality coefficient.

Vahram Stepanyan↗

Small Body GN&C Research Report: A G-Guidance Enhancement to Increase Mission Functionality

G-Guidance is a robust G&C (guidance and control) algorithm developed under the small-body GN&C task. The G-Guidance scheme utilizes a model predictive control approach, along with a convexification of the governing dynamics, control constraints, and trajectory/state constraints. This report details an enhancement to the FY2005 G-Guidance algorithm; the addition of a fire-second approach helps to nullify velocity errors and hit desired final velocities much more precisely than the original Fire-first scheme developed in FY2005. The enhancement preserves the primary benefit of the algorithm, which is to ensure required thruster silent times during trajectory maneuvers. The fire-second scheme increases the versatility of the algorithm for missions employing G-Guidance. For instance, a landing sequence could employ the fire-second scheme to ensure a null final velocity at landing, whereas an ascent sequence could utilize a fire-first scheme to ensure there is no drift of the spacecraft toward the ground at the beginning of the maneuver. Examples are provided within to demonstrate a fire-first versus fire-second guidance scheme. As in the existing G-Guidance algorithm, the examples and schemes incorporate gravity models and thruster firing times into discrete dynamics that are solved as a optimal control problem to minimize fuel consumption or thruster energy expenditure.

velocity errors↗

Boundary conditions for the numerical solution of elliptic equations in exterior regions

Elliptic equations in exterior regions frequently require a boundary condition at infinity to ensure the well-posedness of the problem. Examples of practical applications include the Helmholtz equation and Laplace's equation. Computational procedures based on a direct discretization of the elliptic problem require the replacement of the condition on a finite artificial surface. Direct imposition of the condition at infinity along the finite boundary results in large errors. A sequence of boundary conditions is developed which provides increasingly accurate approximations to the problem in the infinite domain. Estimates of the error due to the finite boundary are obtained for several cases. Computations are presented which demonstrate the increased accuracy that can be obtained by the use of the higher order boundary conditions. The examples are based on a finite element formulation but finite difference methods can also be used.

Bayliss, A.↗

Spectral methods for solution of the boundary-layer equations

The basic principles of spectral methods are reviewed, and their application to the two-dimensional incompressible boundary-layer equations is discussed. It is shown that spectral collocation methods provide an accuracy sufficient for engineering applications on extremely coarse grids, and the machine time requirements are small. Although for a given stream-wise discretization, the fully spectral scheme requires more computational effort than the marching scheme, the higher accuracy of the former yields a large net gain in efficiency when considered at equivalent error levels.

Streett, C. L.↗

Numerical solution of the unsteady Navier-Stokes equation

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws are discussed. These schemes share many desirable properties with total variation diminishing schemes, but TVD schemes have at most first-order accuracy, in the sense of truncation error, at extrema of the solution. In this paper a uniformly second-order approximation is constructed, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Osher, Stanley J.↗

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - 2

The global asymptotic nonlinear behavior of 1 1 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODES) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDES.

Yee, H. C.↗