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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

UDU/T/ covariance factorization for Kalman filtering

There has been strong motivation to produce numerically stable formulations of the Kalman filter algorithms because it has long been known that the original discrete-time Kalman formulas are numerically unreliable. Numerical instability can be avoided by propagating certain factors of the estimate error covariance matrix rather than the covariance matrix itself. This paper documents filter algorithms that correspond to the covariance factorization P = UDU(T), where U is a unit upper triangular matrix and D is diagonal. Emphasis is on computational efficiency and numerical stability, since these properties are of key importance in real-time filter applications. The history of square-root and U-D covariance filters is reviewed. Simple examples are given to illustrate the numerical inadequacy of the Kalman covariance filter algorithms; these examples show how factorization techniques can give improved computational reliability.

Thornton, C. L.↗
NASA NTRS1985

Optimal generalized multistep integration formulae for real-time digital simulation

The problem of discretizing a dynamical system for real-time digital simulation is considered. Treating the system and its simulation as stochastic processes leads to a statistical characterization of simulator fidelity. A plant discretization procedure based on an efficient matrix generalization of explicit linear multistep discrete integration formulae is introduced, which minimizes a weighted sum of the mean squared steady-state and transient error between the system and simulator outputs.

Moerder, D. D.↗
NASA NTRS1987

A modern control theory based algorithm for control of the NASA/JPL 70-meter antenna axis servos

A digital computer-based state variable controller was designed and applied to the 70-m antenna axis servos. The general equations and structure of the algorithm and provisions for alternate position error feedback modes to accommodate intertarget slew, encoder referenced tracking, and precision tracking modes are descibed. Development of the discrete time domain control model and computation of estimator and control gain parameters based on closed loop pole placement criteria are discussed. The new algorithm was successfully implemented and tested in the 70-m antenna at Deep Space Network station 63 in Spain.

Hill, R. E.↗
NASA NTRS1990

Predicting thunderstorm evolution using ground-based lightning detection networks

Lightning measurements acquired principally by a ground-based network of magnetic direction finders are used to diagnose and predict the existence, temporal evolution, and decay of thunderstorms over a wide range of space and time scales extending over four orders of magnitude. The non-linear growth and decay of thunderstorms and their accompanying cloud-to-ground lightning activity is described by the three parameter logistic growth model. The growth rate is shown to be a function of the storm size and duration, and the limiting value of the total lightning activity is related to the available energy in the environment. A new technique is described for removing systematic bearing errors from direction finder data where radar echoes are used to constrain site error correction and optimization (best point estimate) algorithms. A nearest neighbor pattern recognition algorithm is employed to cluster the discrete lightning discharges into storm cells and the advantages and limitations of different clustering strategies for storm identification and tracking are examined.

Goodman, Steven J.↗
NASA NTRS1999

A Robust Locally Preconditioned Semi-Coarsening Multigrid Algorithm for the 2-D Navier-Stokes Equations

The goal of this thesis is to develop an efficient and robust locally preconditioned semi-coarsening multigrid algorithm for the two-dimensional Navier-Stokes equations. This thesis examines the performance of the multigrid algorithm with local preconditioning for an upwind-discretization of the Navier-Stokes equations. A block Jacobi iterative scheme is used because of its high frequency error mode damping ability. At low Mach numbers, the performance of a flux preconditioner is investigated. The flux preconditioner utilizes a new limiting technique based on local information that was developed by Siu. Full-coarsening and-semi-coarsening are examined as well as the multigrid V-cycle and full multigrid. The numerical tests were performed on a NACA 0012 airfoil at a range of Mach numbers. The tests show that semi-coarsening with flux preconditioning is the most efficient and robust combination of coarsening strategy, and iterative scheme - especially at low Mach numbers.

Cain, Michael D.↗
NASA NTRS1994

Eulerian-Lagrangian Simulations of Transonic Flutter Instabilities

This paper presents an overview of recent applications of Eulerian-Lagrangian computational schemes in simulating transonic flutter instabilities. This approach, the fluid-structure system is treated as a single continuum dynamics problem, by switching from an Eulerian to a Lagrangian formulation at the fluid-structure boundary. This computational approach effectively eliminates the phase integration errors associated with previous methods, where the fluid and structure are integrated sequentially using different schemes. The formulation is based on Hamilton's Principle in mixed coordinates, and both finite volume and finite element discretization schemes are considered. Results from numerical simulations of transonic flutter instabilities are presented for isolated wings, thin panels, and turbomachinery blades. The results suggest that the method is capable of reproducing the energy exchange between the fluid and the structure with significantly less error than existing methods. Localized flutter modes and panel flutter modes involving traveling waves can also be simulated effectively with no a priori knowledge of the type of instability involved.

Bendiksen, Oddvar O.↗
NASA NTRS2019

Cloud Scattering Impact on Thermal Radiative Transfer and Global Longwave Radiation

The potential importance of longwave (LW) cloud scattering has been recognized but the actual estimate of this effect on thermal radiation varies greatly among different studies. General circulation models (GCMs) generally neglect or simplify the multiple scattering in the LW. In this study, we use a rigorous radiative transfer algorithm to explicitly consider LW multiple-scattering and apply the GCM to quantify the impact of cloud LW scattering on thermal radiation fluxes. Our study shows that the cloud scattering effect on downward thermal radiation at the surface is concentrated in the infrared atmospheric window spectrum (800–1250 cm9exp −1)). The scattering effect on the outgoing longwave radiation (OLR) is also present in the window region over low clouds but it is mainly in the far-infrared spectrum (300–600 cm(exp −1)) over high clouds. For clouds with small to moderate optical depth (τ < 10), the scattering effect on thermal fluxes shows large variation with the cloud τ and has a maximum at an optical depth of ∼3. For opaque clouds, the scattering effect approaches an asymptote and is smaller and less important. The 2-stream radiative transfer scheme could have an error over 10% with an RMS error around 3.5%–4.0% in the calculated LW flux. This algorithm error of the 2-stream approximation could readily exceed the no-scattering error in the LW, and thus it is worthless to include the time-consuming computation of multiple scattering in a 2-stream radiative transfer scheme. However, the calculation error rapidly decreases as stream number increases and the RMS error in LW flux using the 4-stream scheme is under 0.3%, an accuracy sufficient for most climate studies. We implement the 4-stream discrete-ordinate algorithm in the GISS GCM and run the GCM for 20 years with and without the LW scattering effect, respectively. When cloud LW scattering is included, we find that the global annual mean OLR is reduced by 2.7 W/m(exp 2), and the downward surface flux and the net atmospheric absorption are increased by 1.6 W/m2 and 1.8 W/m(exp 2), respectively. Using one year of ISCCP clouds and running the standalone radiative transfer offline, the global annual mean non-scattering errors in OLR, surface LW downward flux and net atmospheric absorption are 3.6 W/m(exp 2), −1.1 W/m(exp 2), and −2.5 W/m(exp 2), respectively. The global scattering impact of 2.7 W/m(exp 2) on the OLR is small when compared to the typical global OLR value of 240 W/m2, but it is significant when compared to cloud LW radiative forcing (30 W/m2) and net cloud forcing (−14 W/m(exp 2)). Overall, the effect of neglecting scattering on the thermal fluxes is comparable to the reported clear sky radiative effect of doubling CO2.

longwave cloud scattering↗
NASA NTRS1979

State and model error estimation for distributed parameter systems

In-flight estimation of large structure model errors in order to detect inevitable deficiencies in large structure controller/estimator models is discussed. Such an estimation process is particularly applicable in the area of shape control system design required to maintain a prescribed static structural shape and, in addition, suppress dynamic disturbances due to the vehicle vibrational modes. The paper outlines a solution to the problem of static shape estimation where the vehicle shape must be reconstructed from a set of measurements discretely located throughout the structure. The estimation process is based on the principle of least-squares that inherently contains the definition and explicit computation of model error estimates that are optimal in some sense. Consequently, a solution is provided for the problem of estimation of static model errors (e.g., external loads). A generalized formulation applicable to distributed parameters systems is first worked out and then applied to a one-dimensional beam-like structural configuration.

Rodriguez, G.↗
NASA NTRS1979

The discrete-time compensated Kalman filter

A suboptimal dynamic compensator to be used in conjunction with the ordinary discrete-time Kalman filter is derived. The resultant compensated Kalman filter has the property that steady-state bias estimation errors, resulting from modelling errors, are eliminated. The implementation of the compensated Kalman filter involves the use of accumulators in the residual channels in addition to the nominal dynamic model of the stochastic system.

Lee, W.-H.↗
NASA NTRS1980

Discrete Tchebycheff orthonormal polynomials and applications

Discrete Tchebycheff orthonormal polynomials offer a convenient way to make least squares polynomial fits of uniformly spaced discrete data. Computer programs to do so are simple and fast, and appear to be less affected by computer roundoff error, for the higher order fits, than conventional least squares programs. They are useful for any application of polynomial least squares fits: approximation of mathematical functions, noise analysis of radar data, and real time smoothing of noisy data, to name a few.

Lear, W. M.↗
NASA NTRS1983

A rate-controlled teleoperator task with simulated transport delays

A teleoperator-system simulation was used to examine the effects of two control modes (joint-by-joint and resolved-rate), a proximity-display method, and time delays (up to 2 sec) on the control of a five-degree-of-freedom manipulator performing a probe-in-hole alignment task. Four subjects used proportional rotational control and discrete (on-off) translation control with computer-generated visual displays. The proximity display enabled subjects to separate rotational errors from displacement (translation) errors; thus, when the proximity display was used with resolved-rate control, the simulated task was trivial. The time required to perform the simulated task increased linearly with time delay, but time delays had no effect on alignment accuracy. Based on the results of this simulation, several future studies are recommended.

Pennington, J. E.↗
NASA NTRS1996

Error and Complexity Analysis for a Collocation-Grid-Projection Plus Precorrected-FFT Algorithm for Solving Potential Integral Equations with LaPlace or Helmholtz Kernels

In this paper we derive error bounds for a collocation-grid-projection scheme tuned for use in multilevel methods for solving boundary-element discretizations of potential integral equations. The grid-projection scheme is then combined with a precorrected FFT style multilevel method for solving potential integral equations with 1/r and e(sup ikr)/r kernels. A complexity analysis of this combined method is given to show that for homogeneous problems, the method is order n natural log n nearly independent of the kernel. In addition, it is shown analytically and experimentally that for an inhomogeneity generated by a very finely discretized surface, the combined method slows to order n(sup 4/3). Finally, examples are given to show that the collocation-based grid-projection plus precorrected-FFT scheme is competitive with fast-multipole algorithms when considering realistic problems and 1/r kernels, but can be used over a range of spatial frequencies with only a small performance penalty.

Phillips, J. R.↗
NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗
NASA NTRS2017

Adjoint-Based Mesh Adaptation for the Sonic Boom Signature Loudness

The mesh adaptation functionality of FUN3D is utilized to obtain a mesh optimized to calculate sonic boom ground signature loudness. During this process, the coupling between the discrete-adjoints of the computational fluid dynamics tool FUN3D and the atmospheric propagation tool sBOOM is exploited to form the error estimate. This new mesh adaptation methodology will allow generation of suitable meshes adapted to reduce the estimated errors in the ground loudness, which is an optimization metric employed in supersonic aircraft design. This new output-based adaptation could allow new insights into meshing for sonic boom analysis and design, and complements existing output-based adaptation techniques such as adaptation to reduce estimated errors in off-body pressure functional. This effort could also have implications for other coupled multidisciplinary adjoint capabilities (e.g., aeroelasticity) as well as inclusion of propagation specific parameters such as prevailing winds or non-standard atmospheric conditions. Results are discussed in the context of existing methods and appropriate conclusions are drawn as to the efficacy and efficiency of the developed capability.

Rallabhandi, Sriram K.↗
NASA NTRS2018

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗
NASA NTRS1994

Dynamic modeling of orographically induced precipitation

Local orography governs the triggering of cloud formation and the enhancement of processes such as condensation and hydrometeor nucleation and growth in mountainous regions. Intense, lengthy precipitation events are typical upwind of the topographic divide, with sharply decreasing magnitude and duration on the lee side. Differences in mean annual precipitation of several hundred percent between windward slopes of orographic barriers and adjacent valleys or lee side slopes are not unusual. Because much of the streamflow in areas such as the western United States is derived from mountainous areas that are remote and often poorly instrumented, modeling of orographic precipitation has important implications for water resources management. Models of orographically induced precipitation differ by their treatment of atmospheric dynamics and by the extent to which they rely on bulk parameterization of cloud and precipitation physics. Adiabatic ascent and a direct proportionality between efficiency and orographically magnified updrafts are the most frequent assumptions in orographic precipitation modeling. Space-time discretization (i.e., resolution) is a major issue because of the high spatial variability of orographic precipitation. For a specific storm, relative errors as large as 50 to 100% are common in the forecast/hindcast of precipitation intensity and can be even larger in the case of catastrophic storms. When monthly or seasonal timescales are used to evaluate model performance, the magnitude of such errors decreases dramatically, reaching values as low as 10 to 15%. Current research is focusing on the development of data assimilation techniques to incorporate radar and satellite observations, and on the development of aggregation and disaggregation methodologies to address the implications of modeling a multiscale problem at restricted spatial and temporal resolutions.

Barros, Ana Paula↗
NASA NTRS2002

Method and Apparatus for Evaluating the Visual Quality of Processed Digital Video Sequences

A Digital Video Quality (DVQ) apparatus and method that incorporate a model of human visual sensitivity to predict the visibility of artifacts. The DVQ method and apparatus are used for the evaluation of the visual quality of processed digital video sequences and for adaptively controlling the bit rate of the processed digital video sequences without compromising the visual quality. The DVQ apparatus minimizes the required amount of memory and computation. The input to the DVQ apparatus is a pair of color image sequences: an original (R) non-compressed sequence, and a processed (T) sequence. Both sequences (R) and (T) are sampled, cropped, and subjected to color transformations. The sequences are then subjected to blocking and discrete cosine transformation, and the results are transformed to local contrast. The next step is a time filtering operation which implements the human sensitivity to different time frequencies. The results are converted to threshold units by dividing each discrete cosine transform coefficient by its respective visual threshold. At the next stage the two sequences are subtracted to produce an error sequence. The error sequence is subjected to a contrast masking operation, which also depends upon the reference sequence (R). The masked errors can be pooled in various ways to illustrate the perceptual error over various dimensions, and the pooled error can be converted to a visual quality measure.

Watson, Andrew B.↗
NASA NTRS2004

Chandra Observations of the X-ray Environs of SN 1998bw/GRB 980425

We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the GRB. Eight X-ray point sources were localized, three and five each in the original error boxes - S1 and S2 - assigned for variable X-ray counterparts to the GRB by BeppoSAX. The sum of the discrete X-ray sources plus continuous emission in S2 observed by CXO on day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by BeppoSAX. We conclude that S2 is the sum of several variable sources that have not disappeared, and therefore is not associated with the GRB. Within S1, clear evidence is seen for a decline of approximately a factor of 12 between day 200 and day 1281. One of the sources in S1, Sla, is coincident with the well-determined radio location of SN 1998bw, and is certainly the remnant of that explosion. The nature of the other sources is also discussed. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning approximately 1300 days, is obtained by assuming the burst and supernova were coincident at 35.6 Mpc. When this X-ray iight curve is compared with those of the X-ray afterglows of ordinary GRBs, X-ray Flashes, and ordinary supernovae, evidence emerges for at least two classes of lightcurves, perhaps bounding a continuum. By three to ten years, all these phenomena seem to converge on a common X-ray luminosity, possibly indicative of the supernova underlying them all. This convergence strengthens the conclusion that SN 1998bw and GRB 980425 took place in the same object. One possible explanation for the two classes is a (nearly) standard GRB observed at different angles, in which case X-ray afterglows with intermediate luminosities should eventually be discovered. Finally, we comment on the contribution of GRBs to the ULX source population.

Kouveliotou, C.↗