Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIFFERENCE EQUATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Acceleration and Velocity Sensing from Measured Strain

A simple approach for computing acceleration and velocity of a structure from the strain is proposed in this study. First, deflection and slope of the structure are computed from the strain using a two-step theory. Frequencies of the structure are computed from the time histories of strain using a parameter estimation technique together with an autoregressive moving average model. From deflection, slope, and frequencies of the structure, acceleration and velocity of the structure can be obtained using the proposed approach. Simple harmonic motion is assumed for the acceleration computations, and the central difference equation with a linear autoregressive model is used for the computations of velocity. A cantilevered rectangular wing model is used to validate the simple approach. Quality of the computed deflection, acceleration, and velocity values are independent of the number of fibers. The central difference equation with a linear autoregressive model proposed in this study follows the target response with reasonable accuracy. Therefore, the handicap of the backward difference equation, phase shift, is successfully overcome.

shape sensing↗
NASA NTRS1985

Error Reduction Program

The details of a study to select, incorporate and evaluate the best available finite difference scheme to reduce numerical error in combustor performance evaluation codes are described. The combustor performance computer programs chosen were the two dimensional and three dimensional versions of Pratt & Whitney's TEACH code. The criteria used to select schemes required that the difference equations mirror the properties of the governing differential equation, be more accurate than the current hybrid difference scheme, be stable and economical, be compatible with TEACH codes, use only modest amounts of additional storage, and be relatively simple. The methods of assessment used in the selection process consisted of examination of the difference equation, evaluation of the properties of the coefficient matrix, Taylor series analysis, and performance on model problems. Five schemes from the literature and three schemes developed during the course of the study were evaluated. This effort resulted in the incorporation of a scheme in 3D-TEACH which is usuallly more accurate than the hybrid differencing method and never less accurate.

Syed, S. A.↗
NASA NTRS1989

Monte Carlo turbulence simulation for Shuttle reentry studies

A realistic model for the Space Shuttle reentry which will permit more rational selection of the reaction control fuel redlines than did the current turbulence model, is described, with special attention given to the turbulence simulation difference equations, which are revised in this paper. The difference equations for the third-order approximation are presented for both the longitudinal turbulence and the transverse turbulence equations. Each set of equations can be used to generate a probability distribution of fuel remaining at landing.

Campbell, Warren C.↗
NASA NTRS1977

The low-latitude circulation of Mars

A steady-state two-layer atmospheric model is used in an attempt to simulate the circulation of the equatorial regions of Mars, with consideration of both thermal and mechanical orographic effects. A steady-state solution is sought with the aim of producing mean seasonal circulations upon which tides can be superimposed, and the relative importance of orographic forcing at low latitudes is evaluated. It is assumed that the atmospheres over a topographic feature and over flat terrain possess the same radiative and convective properties, simple dissipation laws are introduced which represent a surface drag effect, and six first-order difference equations are obtained and then reduced to two linear complex second-order coupled difference equations. Analysis of the response of the equinoctial basic wind field to a Gaussian orographic feature situated at the equator and zero deg longitude reveals a number of similarities and differences between the Martian and terrestrial responses. Most of the similarities are explained in terms of the dynamics of a rotating atmosphere, while the major differences are at least partially explained by the different optical properties, lower boundaries, and smaller mass of the Martian atmosphere.

Webster, P. J.↗
NASA NTRS1984

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗
NASA NTRS1984

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗
DOE OSTI2026

Second Order System Study

During my education in mathematics, engineering and physics, I learned transform pairs and their usage mechanics but I never remember seeing the derivations of the solutions to second order ordinary differential equations (ODE) and difference equations. A solution to a question posed in a potential funder meeting put me on a path to solving second order systems using the five principal Fourier based methods: Fourier transform (FT), Z-transform (ZT), discrete time Fourier transform (DTFT), discrete Fourier transform (DFT) and Laplace transform (LT).

42 ENGINEERING↗
NASA NTRS2007

The SMM Model as a Boundary Value Problem Using the Discrete Diffusion Equation

A generalized single step stepwise mutation model (SMM) is developed that takes into account an arbitrary initial state to a certain partial difference equation. This is solved in both the approximate continuum limit and the more exact discrete form. A time evolution model is developed for Y DNA or mtDNA that takes into account the reflective boundary modeling minimum microsatellite length and the original difference equation. A comparison is made between the more widely known continuum Gaussian model and a discrete model, which is based on modified Bessel functions of the first kind. A correction is made to the SMM model for the probability that two individuals are related that takes into account a reflecting boundary modeling minimum microsatellite length. This method is generalized to take into account the general n-step model and exact solutions are found. A new model is proposed for the step distribution.

Campbell, Joel↗
DOE OSTI2025

Ocpp 2.0.1. Interim Kpi Calculator

The project is split into four pieces. The first is a raw OCPP log parser. The second is a file splitter. The third is a message parser. The final piece is the Interim KPI calculator. The OCPP log parser was created from two different formats of raw OCPP 2.0.1 data. Its intended purpose is to extract device IDs and OCPP event messages from nontabular text logs. The parser looks for specific substrings in the logs to identify which of the two "standards" it should select from. The KPI generator does not perform any of its calculations in parallel. Instead, we opt for a naive batching approach. The splitter takes the file generated from the parser and creates many smaller files for each of the device IDs in the dataset. This allows the pandas queries in the log formatter to be iterate over a significantly smaller slice of data, increasing performance significantly. The message parser step takes messages from each of the files (containing distinct device IDs) and breaks the message out into pieces. The final result is a file with different columns specifying different attributes of the JSON message. The file is an aggregation of all different devices. This is the most complex portion of the code. The KPI calculator takes the parsed messages, as a single file, and calculates the KPI from that data. An excel file is produced with four sheets. These contain the metrics for Session Success, Charge Start Success, Charge End Success, and Charge Start Time. It includes the metrics for the different equations in the Interim KPI Implementation Guide as well as a weighted sum of the different equations for each KPI (excluding Charge End Success and Charge Start Time).

Quinn, Casey↗
NASA NTRS2000

Double-Plate Penetration Equations

This report compares seven double-plate penetration predictor equations for accuracy and effectiveness of a shield design. Three of the seven are the Johnson Space Center original, modified, and new Cour-Palais equations. The other four are the Nysmith, Lundeberg-Stern-Bristow, Burch, and Wilkinson equations. These equations, except the Wilkinson equation, were derived from test results, with the velocities ranging up to 8 km/sec. Spreadsheet software calculated the projectile diameters for various velocities for the different equations. The results were plotted on projectile diameter versus velocity graphs for the expected orbital debris impact velocities ranging from 2 to 15 km/sec. The new Cour-Palais double-plate penetration equation was compared to the modified Cour-Palais single-plate penetration equation. Then the predictions from each of the seven double-plate penetration equations were compared to each other for a chosen shield design. Finally, these results from the equations were compared with test results performed at the NASA Marshall Space Flight Center. Because the different equations predict a wide range of projectile diameters at any given velocity, it is very difficult to choose the "right" prediction equation for shield configurations other than those exactly used in the equations' development. Although developed for various materials, the penetration equations alone cannot be relied upon to accurately predict the effectiveness of a shield without using hypervelocity impact tests to verify the design.

Hayashida, K. B.↗
NASA NTRS1971

Bayesian recursive image estimation.

A procedure for recursively estimating images that are characterized statistically by the mean and correlation functions associated with the random process representing the brightness level is proposed for the case where the images are corrupted by additive noise. First, a dynamic model is developed with a response characteristic which matches that of the scanner output (the input of the estimator is the output of a horizontal line scanner) in a statistical sense. Such models have the form of an ordinary differential or difference equation with white noise input. An insignificant approximation is introduced by using a constant-coefficient model. The appropriate model is a vector valued difference equation with the solution representing a vector Markov process. The next step is to obtain the minimum mean square estimate of the image by using a Kalman filter. Since the image estimation is an interpolation problem, two successive runs over the observation are performed in opposite directions and the resultant estimates are averaged. Examples are included for illustration.

Nahi, N. E.↗
NASA NTRS1978

Steady state spurious errors in shock-capturing numerical schemes

The behavior of the steady state spurious error modes of the MacCormack scheme and the upwind scheme of Warming and Beam was obtained from a linearized difference equation for the steady state error. It was shown that the spurious errors can exist either as an eigensolution of the homogeneous part of this difference equation or because of excitation from large discretization errors near oblique shocks. It was found that the upwind scheme does not permit spurious oscillations on the upstream side of shocks. Examples are given for the inviscid Burgers' equation and for one and two dimensional gasdynamic flows.

Karlsen, L. K.↗
NASA NTRS2006

On the equivalence of the dual-wavelength and polarimetric equations for estimation of the raindrop size distribution

In writing the integral equations for the median mass diameter and particle concentration, or comparable parameters of the raindrop size distribution, it is apparent that when attenuation effects are included, the forms of the equations for polarimetric and dual wavelength radars are identical. In both sets of equations, differences in the backscattering and extinction cross sections appear: in the polarimetric equations, the differences are taken with respect polarization at a fixed frequency while for the dual wavelength equations, the differences are taken with respect to wavelength at a fixed polarization. Because the forms of the equations are the same, the ways in which they can be solved are similar as well. To avoid instabilities in the forward recursion procedure, the equations can be expressed in the form of a final-value. Solving the equations in this way traditionally has required estimates of the path attenuations to the final gate: either the attenuations at horizontal and vertical polarizations at the same frequency or attenuations at two frequencies with the same polarization. This has been done for dual-frequency (air/spaceborne case) and polarimetric radars by the respective use of the surface reference technique and the differential phase shift. An alternative to solving the constrained version of the equations is an iterative procedure recently proposed in which independent estimates of path attenuation are not required. Although the procedure has limitations, it appears to be quite useful. Simulations of the retrievals help clarify the relationship between the constrained and unconstrained approaches and their application to the polarimetric and dual-wavelength equations.

Meneghini, Robert↗
NASA NTRS1977

Finite-difference theory for sound propagation in a lined duct with uniform flow using the wave envelope concept

Finite difference equations are derived for sound propagation in a two dimensional, straight, soft wall duct with a uniform flow by using the wave envelope concept. This concept reduces the required number of finite difference grid points by one to two orders of magnitude depending on the length of the duct and the frequency of the sound. The governing acoustic difference equations in complex notation are derived. An exit condition is developed that allows a duct of finite length to simulate the wave propagation in an infinitely long duct. Sample calculations presented for a plane wave incident upon the acoustic liner show the numerical theory to be in good agreement with closed form analytical theory. Complete pressure and velocity printouts are given to some sample problems and can be used to debug and check future computer programs.

Baumeister, K. J.↗
NASA NTRS1979

On reliable control system designs with and without feedback reconfigurations

This paper contains an overview of a theoretical framework for the design of reliable multivariable control systems, with special emphasis on actuator failures and necessary actuator redundancy levels. Using a linear model of the system, with Markovian failure probabilities and quadratic performance index, an optimal stochastic control problem is posed and solved. The solution requires the iteration of a set of highly coupled Riccati-like matrix difference equations; if these converge one has a reliable design; if they diverge, the design is unreliable, and the system design cannot be stabilized. In addition, it is shown that the existence of a stabilizing constant feedback gain and the reliability of its implementation is equivalent to the convergence properties of a set of coupled Riccati-like matrix difference equations. In summary, these results can be used for offline studies relating the open loop dynamics, required performance, actuator mean time to failure, and functional or identical actuator redundancy, with and without feedback gain reconfiguration strategies.

Birdwell, J. D.↗