Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximation algorithm”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

The straightforward automatic-differentiation and the hand-differentiated incremental iterative methods are interwoven to produce a hybrid scheme that captures some of the strengths of each strategy. With this compromise, discrete aerodynamic sensitivity derivatives are calculated with the efficient incremental iterative solution algorithm of the original flow code. Moreover, the principal advantage of automatic differentiation is retained (i.e., all complicated source code for the derivative calculations is constructed quickly with accuracy). The basic equations for second-order sensitivity derivatives are presented; four methods are compared. Each scheme requires that large systems are solved first for the first-order derivatives and, in all but one method, for the first-order adjoint variables. Of these latter three schemes, two require no solutions of large systems thereafter. For the other two for which additional systems are solved, the equations and solution procedures are analogous to those for the first order derivatives. From a practical viewpoint, implementation of the second-order methods is feasible only with software tools such as automatic differentiation, because of the extreme complexity and large number of terms. First- and second-order sensitivities are calculated accurately for two airfoil problems, including a turbulent flow example; both geometric-shape and flow-condition design variables are considered. Several methods are tested; results are compared on the basis of accuracy, computational time, and computer memory. For first-order derivatives, the hybrid incremental iterative scheme obtained with automatic differentiation is competitive with the best hand-differentiated method; for six independent variables, it is at least two to four times faster than central finite differences and requires only 60 percent more memory than the original code; the performance is expected to improve further in the future.

Taylor, Arthur C., III↗
NASA NTRS1988

Optimal reentry guidance for aeroassisted orbit transfer vehicles

A three-state model is presented for analyzing the problem of optimal changes in heading with minimum energy loss for a hypersonic gliding vehicle. A further model order reduction to a single state model is examined using singular perturbation theory. The optimal solution for the reduced problem defines an optimal altitude profile dependent on the current energy of the vehicle. A separate boundary-layer analysis is used to account for altitude and flight path angle dynamics, and to obtain lift and bank angle control solutions. By considering alternative approximations to solve the boundary-layer problem, three guidance laws are obtained, each having a feedback form. The guidance laws are evaluated for a hypothetical vehicle, and compared to an optimal solution obtained using a multiple shooting algorithm.

Calise, Anthony J.↗
NASA NTRS1991

l1-optimal control of multivariable systems with output norm constraints

This paper considers the l1-optimal control problem for general rational plants. It is shown that, for plants with no poles or zeros on the unit circle, an optimal compensator exists and that the resulting closed loop transfer function is polynomial whenever there are at least as many controls as regulated outputs and at least as many measurements as exogeneous inputs. Exactly or approximately, optimal rational compensators can be obtained by solving a sequence of finite linear programs for the coefficients of a polynomial closed-loop transfer function. No assumptions on plant poles or zeros are required to obtain at least approximately optimal compensators. It is shown that constrained problems in which a set of outputs is regulated subject to l(infinity)-norm constraints on another set of outputs can be solved using a slight modification of the same algorithm.

Mcdonald, J. S.↗
NASA NTRS1992

Multigrid calculation of internal flows in complex geometries

The development, validation, and application of a general purpose multigrid solution algorithm and computer program for the computation of elliptic flows in complex geometries is presented. This computer program combines several desirable features including a curvilinear coordinate system, collocated arrangement of the variables, and Full Multi-Grid/Full Approximation Scheme (FMG/FAS). Provisions are made for the inclusion of embedded obstacles and baffles inside the flow domain. The momentum and continuity equations are solved in a decoupled manner and a pressure corrective equation is used to update the pressures such that the fluxes at the cell faces satisfy local mass continuity. Despite the computational overhead required in the restriction and prolongation phases of the multigrid cycling, the superior convergence results in reduced overall CPU time. The numerical scheme and selected results of several validation flows are presented. Finally, the procedure is applied to study the flowfield in a side-inlet dump combustor and twin jet impingement from a simulated aircraft fuselage.

Smith, K. M.↗
NASA NTRS1994

An IUE survey of interstellar H I Ly alpha absorption. 1: Column densities

We measure Galactic interstellar neutral hydrogen column densities by analyzing archival interstellar Ly alpha absorption line data toward 554 B2 and hotter stars observed at high resolution with the IUE satellite. This study more than doubles the number of lines of sight with measures of N(H I) based on Ly alpha. We have included the scattered light background correction algorithm of Bianchi and Bohlin in our data reduction. We use the correlation between the Balmer discontinuity (c sub 1) index and the stellar Ly alpha absorption in order to assess the effects of stellar Ly alpha contamination. Approximately 40% of the B stars with measured (c sub 1) index, exhibit serious stellar Ly alpha contamination. One table contains the derived values of the interstellar N(H I) for 393 stars with at most small amounts of stellar contamination. Another lists the observed values of total N(H I) for 161 stars with suspected stellar Ly alpha contamination and/or uncertain stellar parameters.

Diplas, Athanassios↗
NASA NTRS1994

Rapid computation of spectrally integrated non-local thermodynamic equilibrium limb emission

The interpretation of infrared radiance measurements made by satellite-borne limb-scanning broadband radiometers requires accurate and computationally fast techniques with which to evaluate the equation of radiative transfer. This requirement is made even more stringent when analyzing measurements of non-local thermodynamic equilibrium (non-LTE) emission from the terrestrial mesosphere and lower thermosphere. In principle, line-by-line calculations which explicitly account for the departure from thermodynamic equilibrium in both the source functions and the transmittances are necessary. In this paper we extend the emissivity growth approximation (EGA) technique developed for local thermodynamic equilibrium (LTE) limb radiance for the molecular oxygen dayglow (1.27 micrometers and 762 nm), ozone and carbon dioxide in the 9- to 11-micrometer spectral interval, carbon monoxide (4.6 micrometers), nitric oxide (5.3 micrometers), and the carbon dioxide bands (15 micrometers) are presented. Using the non-LTE form of the EGA, the spectrally integrated limb emission is calculated for 35 tangent heights in the mesosphere and lower thermosphere (a total of 1200 atmospheric layers) with line-by-line accuracy in approximately 0.35 s of CPU time on readily available desktop computer hardware, while the corresponding line-by-line calculations may require several minutes. The non-LTE EGA technique will allow kinetic temperature and minor constituend retrieval algorithms to readily include non-LTE efects limited only by the a priori knowledge of the departure from LTE in the observed bands.

Mlynczak, Martin G.↗
NASA NTRS1999

Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions

We examine the nonlinear stability of the Wisdom-Holman (WH) symplectic mapping applied to the integration of perturbed, highly eccentric (e-0.9) two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories for this class of problems, unless the step size chosen 1s small enough that PeriaPse is always resolved, in which case the method is generically stable. This 'radial orbit instability' persists even for weakly perturbed systems. Using the Stark problem as a fiducial test case, we investigate the dynamical origin of this instability and argue that the numerical chaos results from the overlap of step-size resonances; interestingly, for the Stark-problem many of these resonances appear to be absolutely stable. We similarly examine the robustness of several alternative integration methods: a time-regularized version of the WH mapping suggested by Mikkola; the potential-splitting (PS) method of Duncan, Levison, Lee; and two original methods incorporating approximations based on Stark motion instead of Keplerian motion. The two fixed point problem and a related, more general problem are used to conduct a comparative test of the various methods for several types of motion. Among the algorithms tested, the time-transformed WH mapping is clearly the most efficient and stable method of integrating eccentric, nearly Keplerian orbits in the absence of close encounters. For test particles subject to both high eccentricities and very close encounters, we find an enhanced version of the PS method-incorporating time regularization, force-center switching, and an improved kernel function-to be both economical and highly versatile. We conclude that Stark-based methods are of marginal utility in N-body type integrations. Additional implications for the symplectic integration of N-body systems are discussed.

Rauch, Kevin P.↗
NASA NTRS2015

Modeling the Swift BAT Trigger Algorithm with Machine Learning

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien et al. (2014) is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of approximately greater than 97% (approximately less than 3% error), which is a significant improvement on a cut in GRB flux which has an accuracy of 89:6% (10:4% error). These models are then used to measure the detection efficiency of Swift as a function of redshift z, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of eta(sub 0) approximately 0.48(+0.41/-0.23) Gpc(exp -3) yr(exp -1) with power-law indices of eta(sub 1) approximately 1.7(+0.6/-0.5) and eta(sub 2) approximately -5.9(+5.7/-0.1) for GRBs above and below a break point of z(sub 1) approximately 6.8(+2.8/-3.2). This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting. The code used in this is analysis is publicly available online.

gamma rays: general↗
NASA NTRS2018

Formation Flight System Extremum-Seeking-Control Using Blended Performance Parameters

An extremum-seeking control system for formation flight that uses blended performance parameters in a conglomerate performance function that better approximates drag reduction than performance functions formed from individual measurements. Generally, a variety of different measurements are taken and fed to a control system, the measurements are weighted, and are then subjected to a peak-seeking control algorithm. As measurements are continually taken, the aircraft will be guided to a relative position which optimizes the drag reduction of the formation. Two embodiments are discussed. Two approaches are shown for determining relative weightings: "a priori" by which they are qualitatively determined (by minimizing the error between the conglomerate function and the drag reduction function), and by periodically updating the weightings as the formation evolves.

Ryan, John J.↗
NASA NTRS2022

Uncertainty Estimation for SMAP Level-1 Brightness Temperature Assimilation at Different Timescales

Soil Moisture Active Passive (SMAP) mission brightness temperature (T(b) ) observations are assimilated into NASA’s Catchment Land Surface Model using an Ensemble Kalman filter to update simulations of surface and root-zone soil moisture. Different time series components of the T(b) observations are assimilated including anomalies, inter-annual variations, and high frequency variations. To optimize the weights that the data assimilation (DA) puts on the observations, the ratio between the uncertainties of modeled and observed T(b) is approximated using modeled and observed soil moisture uncertainties estimated using triple collocation analysis. In a benchmark experiment, T(b) observations are assimilated using a spatially constant 4 Kelvin (K) observation uncertainty, as in the operational SMAP Level-4 algorithm. All DA experiments exhibit notable skill improvements in most regions. Improvements are largest for the inter-annual variations in the simulations of both surface and root-zone soil moisture (mean improvements in terms of Pearson correlation (-) are 0.08 and 0.06, respectively). Anomaly simulations improve similarly (0.07), and improvements in the high-frequency variations are only observed for surface soil moisture simulations (0.06). No notable difference in skill - neither improvement nor deterioration - is observed between the experiments that use optimized observation uncertainty parameters and the 4 K benchmark experiment. This may be explained by the presence of large observation operator errors, which are analytically shown to have the potential to render post-update uncertainty insensitive to inaccuracies in estimates of the Kalman gain. These results have important implications for the design of soil moisture DA systems, in particular for parameterizing model and observation uncertainties.

Hydrology↗
NASA NTRS2016

Management of Asymptomatic Renal Stones in Astronauts

Introduction: Management guidelines were created to screen and manage asymptomatic renal stones in U.S. astronauts. The risks for renal stone formation in astronauts due to bone loss and hypercalcuria are unknown. Astronauts have a stone risk which is about the same as commercial aviation pilots, which is about half that of the general population. However, proper management of this condition is still crucial to mitigate health and mission risks in the spaceflight environment. Methods: An extensive review of the literature and current aeromedical standards for the monitoring and management of renal stones was done. The NASA Flight Medicine Clinic's electronic medical record and Longitudinal Survey of Astronaut Health were also reviewed. Using this work, a screening and management algorithm was created that takes into consideration the unique operational environment of spaceflight. Results: Renal stone screening and management guidelines for astronauts were created based on accepted standards of care, with consideration to the environment of spaceflight. In the proposed algorithm, all astronauts will receive a yearly screening ultrasound for renal calcifications, or mineralized renal material (MRM). Any areas of MRM, 3 millimeters or larger, are considered a positive finding. Three millimeters approaches the detection limit of standard ultrasound, and several studies have shown that any stone that is 3 millimeters or less has an approximately 95 percent chance of spontaneous passage. For mission-assigned astronauts, any positive ultrasound study is followed by low-dose renal computed tomography (CT) scan, and flexible ureteroscopy if CT is positive. Other specific guidelines were also created. Discussion: The term "MRM" is used to account for small areas of calcification that may be outside the renal collecting system, and allows objectivity without otherwise constraining the diagnostic and treatment process for potentially very small calcifications of uncertain significance. However, a small asymptomatic MRM or stone within the renal collecting system may become symptomatic, and so affect launch and flight schedules, cause incapacitation during flight, and ultimately require medical evacuation. For exploration class missions, evacuation is unlikely. The new screening and management algorithm allows better management of mission risks, and will define the true incidence of renal stones in U.S. astronauts. This information will be used to refine future screening, countermeasures and treatment methods; and will also inform the needed capabilities to be flown on exploration-class missions.

Reyes, David↗
NASA NTRS2016

Management of Asymptomatic Renal Stones in Astronauts

Introduction: Management guidelines were created to screen and manage asymptomatic renal stones in U.S. astronauts. The risks for renal stone formation in astronauts due to bone loss and hypercalcuria are unknown. Astronauts have a stone risk which is about the same as commercial aviation pilots, which is about half that of the general population. However, proper management of this condition is still crucial to mitigate health and mission risks in the spaceflight environment. Methods: An extensive review of the literature and current aeromedical standards for the monitoring and management of renal stones was done. The NASA Flight Medicine Clinic's electronic medical record and Longitudinal Survey of Astronaut Health were also reviewed. Using this work, a screening and management algorithm was created that takes into consideration the unique operational environment of spaceflight. Results: Renal stone screening and management guidelines for astronauts were created based on accepted standards of care, with consideration to the environment of spaceflight. In the proposed algorithm, all astronauts will receive a yearly screening ultrasound for renal calcifications, or mineralized renal material (MRM). Any areas of MRM, 3 millimeters or larger, are considered a positive finding. Three millimeters approaches the detection limit of standard ultrasound, and several studies have shown that any stone that is 3 millimeters or less has an approximately 95 percent chance of spontaneous passage. For mission-assigned astronauts, any positive ultrasound study is followed by low-dose renal computed tomography (CT) scan, and flexible ureteroscopy if CT is positive. Other specific guidelines were also created. Discussion: The term "MRM" is used to account for small areas of calcification that may be outside the renal collecting system, and allows objectivity without otherwise constraining the diagnostic and treatment process for potentially very small calcifications of uncertain significance. However, a small asymptomatic MRM or stone within the renal collecting system may become symptomatic, and so affect launch and flight schedules, cause incapacitation during flight, and ultimately require medical evacuation. For exploration class missions, evacuation is unlikely. The new screening and management algorithm allows better management of mission risks, and will define the true incidence of renal stones in U.S. astronauts. This information will be used to refine future screening, countermeasures and treatment methods; and will also inform the needed capabilities to be flown on exploration-class missions.

Reyes, David↗
NASA NTRS2018

Gravitational Wave Sources as Timing References for LISA Data

In the megahertz gravitational-wave band, galactic ultra-compact binaries (UCBs) are continuous sources emitting at near-constant frequency. The signals from many of these galactic binaries will be sufficiently strong to be detectable by the Laser Interferometer Space Antenna (LISA) after approximately Omicron (1 week) of observing. In addition to their astrophysical value, these UCBs can be used to monitor the data quality of the observatory. This paper demonstrates the capabilities of galactic UCBs to be used as calibration sources for LISA by demanding signal coherence between adjacent week-long data segments separated by a gap in time of a priori unknown duration. A parameter for the gap duration is added to the UCB waveform model and used in a Markov-chain Monte Carlo algorithm simultaneously fitting for the astrophysical source parameters. Results from measurements of several UCBs are combined to produce a joint posterior on the gap duration. The measurement accuracy's dependence on how much is known about the UCBs through prior observing, and seasonal variations due to the LISA orbital motion, is quantified. The duration of data gaps in a two-week segment of data can be constrained to within approximately 0.2 s using Omicron (10) UCBs after one month of observing. The timing accuracy from UCBs improves to approximately or less than o.1 s after 1 year of mission operations. These results are robust to within a factor of approximately 2 when taking into account seasonal variations.

Gravitation Wave Sources↗
NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗
NASA NTRS1983

A Kalman filter algorithm for terminal-area navigation with sensors of moderate accuracy

Translational state estimation in terminal area operations, using a set of commonly available position, air data, nd acceleration sensors, is described. Kalman filtering is applied to obtain maximum estimation occuracy from the sensors but feasibility in real-time computations requires a variety of approximations and devices aimed at minimizing the required computation time with only negligible loss of accuracy. Accuracy behavior throughout the terminal area, its relation to sensor accuracy, its effect on trajectory tracking errors and control activity in an automatic flight control system, and its adequacy in terms of existing criteria for various terminal area operations are examined. The principal investigative tool is a simulation of the system. Previously announced in STAR as N83-29193

Cicolani, L. S.↗
NASA NTRS1990

Optical neural network system for pose determination of spinning satellites

An optical neural network architecture and algorithm based on a Hopfield optimization network are presented for multitarget tracking. This tracker utilizes a neuron for every possible target track, and a quadratic energy function of neural activities which is minimized using gradient descent neural evolution. The neural net tracker is demonstrated as part of a system for determining position and orientation (pose) of spinning satellites with respect to a robotic spacecraft. The input to the system is time sequence video from a single camera. Novelty detection and filtering are utilized to locate and segment novel regions from the input images. The neural net multitarget tracker determines the correspondences (or tracks) of the novel regions as a function of time, and hence the paths of object (satellite) parts. The path traced out by a given part or region is approximately elliptical in image space, and the position, shape and orientation of the ellipse are functions of the satellite geometry and its pose. Having a geometric model of the satellite, and the elliptical path of a part in image space, the three-dimensional pose of the satellite is determined. Digital simulation results using this algorithm are presented for various satellite poses and lighting conditions.

Lee, Andrew↗
NASA NTRS2009

Combining Satellite Microwave Radiometer and Radar Observations to Estimate Atmospheric Latent Heating Profiles

In this study, satellite passive microwave sensor observations from the TRMM Microwave Imager (TMI) are utilized to make estimates of latent + eddy sensible heating rates (Q1-QR) in regions of precipitation. The TMI heating algorithm (TRAIN) is calibrated, or "trained" using relatively accurate estimates of heating based upon spaceborne Precipitation Radar (PR) observations collocated with the TMI observations over a one-month period. The heating estimation technique is based upon a previously described Bayesian methodology, but with improvements in supporting cloud-resolving model simulations, an adjustment of precipitation echo tops to compensate for model biases, and a separate scaling of convective and stratiform heating components that leads to an approximate balance between estimated vertically-integrated condensation and surface precipitation. Estimates of Q1-QR from TMI compare favorably with the PR training estimates and show only modest sensitivity to the cloud-resolving model simulations of heating used to construct the training data. Moreover, the net condensation in the corresponding annual mean satellite latent heating profile is within a few percent of the annual mean surface precipitation rate over the tropical and subtropical oceans where the algorithm is applied. Comparisons of Q1 produced by combining TMI Q1-QR with independently derived estimates of QR show reasonable agreement with rawinsonde-based analyses of Q1 from two field campaigns, although the satellite estimates exhibit heating profile structure with sharper and more intense heating peaks than the rawinsonde estimates. 2

Grecu, Mircea↗
NASA NTRS1973

Optimal control of a variable spin speed CMG system for space vehicles

Many future NASA programs require very high accurate pointing stability. These pointing requirements are well beyond anything attempted to date. This paper suggests a control system which has the capability of meeting these requirements. An optimal control law for the suggested system is specified. However, since no direct method of solution is known for this complicated system, a computation technique using successive approximations is used to develop the required solution. The method of calculus of variations is applied for estimating the changes of index of performance as well as those constraints of inequality of state variables and terminal conditions. Thus, an algorithm is obtained by the steepest descent method and/or conjugate gradient method. Numerical examples are given to show the optimal controls.

Liu, T. C.↗