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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Modeling of composite beams and plates for static and dynamic analysis

The main purpose of this research was to develop a rigorous theory and corresponding computational algorithms for through-the-thickness analysis of composite plates. This type of analysis is needed in order to find the elastic stiffness constants for a plate and to post-process the resulting plate solution in order to find approximate three-dimensional displacement, strain, and stress distributions throughout the plate. This also requires the development of finite deformation plate equations which are compatible with the through-the-thickness analyses. After about one year's work, we settled on the variational-asymptotical method (VAM) as a suitable framework in which to solve these types of problems. VAM was applied to laminated plates with constant thickness in the work of Atilgan and Hodges. The corresponding geometrically nonlinear global deformation analysis of plates was developed by Hodges, Atilgan, and Danielson. A different application of VAM, along with numerical results, was obtained by Hodges, Lee, and Atilgan. An expanded version of this last paper was submitted for publication in the AIAA Journal.

Hodges, Dewey H.↗
NASA NTRS1993

Analysis of the distribution of pitch angles in model galactic disks - Numerical methods and algorithms

An automated mathematical method capable of successfully isolating the many different features in prototype and observed spiral galaxies and of accurately measuring the pitch angles and lengths of these individual features is developed. The method is applied to analyze the evolution of specific features in a prototype galaxy exhibiting flocculent spiral structure. The mathematical-computational method was separated into two components. Initially, the galaxy was partitioned into dense regions constituting features using two different methods. The results obtained using these two partitioning algorithms were very similar, from which it is inferred that no numerical biasing was evident and that capturing of the features was consistent. Standard least-squares methods underestimated the true slope of the cloud distribution and were incapable of approximating an orientation of 45 deg. The problems were overcome by introducing a superior fit least-squares method, developed with the intention of calculating true orientation rather than a regression line.

Russell, William S.↗
NASA NTRS1994

Multigrid one shot methods for optimal control problems: Infinite dimensional control

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

Arian, Eyal↗
NASA NTRS2003

Transition-Independent Decentralized Markov Decision Processes

There has been substantial progress with formal models for sequential decision making by individual agents using the Markov decision process (MDP). However, similar treatment of multi-agent systems is lacking. A recent complexity result, showing that solving decentralized MDPs is NEXP-hard, provides a partial explanation. To overcome this complexity barrier, we identify a general class of transition-independent decentralized MDPs that is widely applicable. The class consists of independent collaborating agents that are tied up by a global reward function that depends on both of their histories. We present a novel algorithm for solving this class of problems and examine its properties. The result is the first effective technique to solve optimally a class of decentralized MDPs. This lays the foundation for further work in this area on both exact and approximate solutions.

Becker, Raphen↗
NASA NTRS2008

Development of arrays of position-sensitive microcalorimeters for Constellation-X

We are developing arrays of position-sensitive transition-edge sensor (POST) X-ray detectors for future astronomy missions such as NASA's Constellation-X. The POST consists of multiple absorbers thermally coupled to one or more transition-edge sensor (TES). Each absorber element has a different thermal coupling to the TES. This results in a distribution of different pulse shapes and enables position discrimination between the absorber elements. POST'S are motivated by the desire to achieve the largest possible focal plane area with the fewest number of readout channels and are ideally suited to increasing the Constellation-X focal plane area, without comprising on spatial sampling. Optimizing the performance of POST'S requires careful design of key parameters such as the thermal conductances between the absorbers, TES and the heat sink. as well as the absorber heat capacities. Using recently developed signal processing algorithms we have investigated the trade-off between position-sensitivity, energy resolution and pulse decay time. based on different device design parameters for PoST's. Our new generation of PoST's utilize technology successfully developed on high resolution (approximately 2.5eV) single pixels arrays of Mo/Au TESs. also under development for Constellation-X. This includes noise mitigation features on the TES and low resistivity electroplated absorbers. We report on the first experimental results from these new one and two-channel PoST"s, consisting of all Au and composite Au/Bi absorbers, which are designed to achieve an energy resolution of < 10 eV. coupled with count-rates of 100's per pixel per second and position sensitivity over the energy range 0.3-10 keV.

Smith, S. J.↗
NASA NTRS2015

Mars Double-Flyby Free Returns

A subset of Earth-originating Mars double-flyby ballistic trajectories is documented. The subset consists of those trajectories that, after the first Mars flyby, perform a half-revolution transfer with Mars before returning to Earth. This class of free returns is useful for both human and robotic Mars missions because of its low geocentric energy at departure and arrival, and because of its extended stay time in the vicinity of Mars. Ballistic opportunities are documented over Earth departure dates ranging from 2015 through 2100. The mission is viable over three or four consecutive Mars synodic periods and unavailable for the next four, with the pattern repeating approximately every 15 years. Over the remainder of the century, a minimum Earth departure hyperbolic excess speed of 3.16 km/s, a minimum Earth atmospheric entry speed of 11.47 km/s, and a minimum flight time of 904 days are observed. The algorithm used to construct these trajectories is presented along with several examples.

Jesick, Mark↗
NASA NTRSDate not supplied

Topographic Slant Range Modeling and Fault Detection for Precision Planetary Landing

This work presents a novel landing site relative topographic measurement model for aslant range sensor being utilized for precision planetary landing operations. The measurement model accounts for the local terrain the slant range sensor captures and leverages knowledge of the estimated landing site provided by the navigation filter. Notably, in contrast to previous works, the new model does not rely on surface normal approximation, reducing the model’s sensitivity to noisy digital elevation maps which represent the local topography. In addition to the measurement model, this work introduces a novel fault detection method, denoted the probabilistic inspection of topographic filter altitude likelihood (PITFAL) algorithm, that implements a statistical outlier rejection algorithm. PITFAL is designed for multi-beam slant range sensors, such as the Navigation Doppler LIDAR (NDL), and identifies statistically inconsistent range estimates through a consensus check on the set of apparent altitudes computed for each individual beam. These models are numerically validated by the Safe and Precise Landing Capability Evolution (SPLICE) project’s high-fidelity terrestrial and lunar lander simulations.

Davis W Adams↗
NASA NTRS2025

Topographic Slant Range Modeling and Fault Detection for Precision Planetary Landing

This work presents a novel landing site relative topographic measurement model for aslant range sensor being utilized for precision planetary landing operations. The measurement model accounts for the local terrain the slant range sensor captures and leverages knowledge of the estimated landing site provided by the navigation filter. Notably, in contrast to previous works, the new model does not rely on surface normal approximation, reducing the model’s sensitivity to noisy digital elevation maps which represent the local topography. In addition to the measurement model, this work introduces a novel fault detection method, denoted the probabilistic inspection of topographic filter altitude likelihood (PITFAL) algorithm, that implements a statistical outlier rejection algorithm. PITFAL is designed for multi-beam slant range sensors, such as the Navigation Doppler LIDAR (NDL), and identifies statistically inconsistent range estimates through a consensus check on the set of apparent altitudes computed for each individual beam. These models are numerically validated by the Safe and Precise Landing Capability Evolution (SPLICE) project’s high-fidelity terrestrial and lunar lander simulations.

Davis W Adams↗
NASA NTRS2013

Absorption Properties of Mediterranean Aerosols Obtained from Multi-year Ground-based and Satellite Remote Sensing Observations

Aerosol absorption properties are of high importance to assess aerosol impact on regional climate. This study presents an analysis of aerosol absorption products obtained over the Mediterranean Basin or land stations in the region from multi-year ground-based AERONET and satellite observations with a focus on the Absorbing Aerosol Optical Depth (AAOD), Single Scattering Albedo (SSA) and their spectral dependence. The AAOD and Absorption Angstrom Exponent (AAE) data set is composed of daily averaged AERONET level 2 data from a total of 22 Mediterranean stations having long time series, mainly under the influence of urban-industrial aerosols and/or soil dust. This data set covers the 17 yr period 1996-2012 with most data being from 2003-2011 (approximately 89 percent of level-2 AAOD data). Since AERONET level-2 absorption products require a high aerosol load (AOD at 440 nm greater than 0.4), which is most often related to the presence of desert dust, we also consider level-1.5 SSA data, despite their higher uncertainty, and filter out data with an Angstrom exponent less than 1.0 in order to study absorption by carbonaceous aerosols. The SSA data set includes both AERONET level-2 and satellite level-3 products. Satellite-derived SSA data considered are monthly level 3 products mapped at the regional scale for the spring and summer seasons that exhibit the largest aerosol loads. The satellite SSA dataset includes the following products: (i) Multi-angle Imaging SpectroRadiometer (MISR) over 2000-2011, (ii) Ozone Monitoring Instrument (OMI) near-UV algorithm over 2004-2010, and (iii) MODerate resolution Imaging Spectroradiometer (MODIS) Deep-Blue algorithm over 2005-2011, derived only over land in dusty conditions. Sun-photometer observations show that values of AAOD at 440 nm vary between 0.024 +/- 0.01 (resp. 0.040 +/- 0.01) and 0.050 +/- 0.01 (0.055 +/- 0.01) for urban (dusty) sites. Analysis shows that the Mediterranean urban-industrial aerosols appear "moderately" absorbing with values of SSA close to approximately 0.94-0.95 +/- 0.04 (at 440 nm) in most cases except over the large cities of Rome and Athens, where aerosol appears more absorbing (SSA approximately 0.89-0.90 +/- 0.04). The aerosol Absorption Angstrom Exponent (AAE, estimated using 440 and 870 nm) is found to be larger than 1 for most sites over the Mediterranean, a manifestation of mineral dust (iron) and/or brown carbon producing the observed absorption. AERONET level-2 sun-photometer data indicate the existence of a moderate East-West gradient, with higher values over the eastern basin (AAEEast. = 1.39/AAEWest. = 1.33) due to the influence of desert dust. The North-South AAE gradient is more pronounced, especially over the western basin. Our additional analysis of AERONET level-1.5 data also shows that organic absorbing aerosols significantly affect some Mediterranean sites. These results indicate that current climate models treating organics as nonabsorbing over the Mediterranean certainly underestimate the warming effect due to carbonaceous aerosols. A~comparative analysis of the regional SSA variability has been attempted using satellite data. OMI and MODIS data show an absorbing zone (SSA approximately 0.90 at 470-500 nm) over Northeastern Africa that does not appear in the MISR retrievals. In contrast, MISR seems able to observe the East-West SSA gradient during summer, as also detected by AERONET. Also, the analysis of SSA provided by satellites indicates that the aerosol over the Mediterranean Sea appears less absorbing during spring (MAM) than summer (JJA).

satellite remote sensing observations↗
NASA NTRS2001

Theoretical Basis for the Surface Spectral Reflectance Relationships Used in the MODIS Aerosol Algorithm

The analysis of data from the MODIS instrument on the Terra platform to derive global distribution of aerosols assumes a set of relationships between the blue, rho (sub blue), the red, rho (sub red), and 2.1 micrometers, rho (sub 2.1), spectral channels. These relations have been established from a series of measurements indicating that rho (sub blue) approximately 0.5 rho (sub red) approximately 0.25 rho (sub 2.1). Here we use a model to describe the transfer of radiation through a vegetation canopy composed of randomly oriented leaves to assess the theoretical foundations for these relationships. The influence of varying fractional vegetation coverage is simulated simply as a linear combination of pure soil and pure vegetation conditions, also known as Independent Pixel Approximation (IPA). Calculations for a wide range of leaf area indices and vegetation fractions show that rho (sub blue) is consistently about 1/4 of rho (sub 2.1) as used by MODIS for the whole range of analyzed cases, except for very dark soils, such as those found in burn scars. For its part, the ratio rho (sub red)/rho (sub 2.1) varies from less than the empirically derived value of 1/2 for dense and dark vegetation (rho (sub 2.1) less than 0.1), to more than 1/2 for bright mixture of soil and vegetation. This is in agreement with measurements over uniform dense vegetation, but not with measurements over mixed dark scenes. In the later case, the discrepancy is probably mitigated by shadows due to uneven canopy and terrain on a large scale. It is concluded that the value of this ratio should ideally be made dependent on the land cover type in the operational processing of MODIS data, especially over dense forests.

Kaufman, Yoram J.↗
NASA NTRS1998

A Parallel Non-Overlapping Domain-Decomposition Algorithm for Compressible Fluid Flow Problems on Triangulated Domains

This paper considers an algebraic preconditioning algorithm for hyperbolic-elliptic fluid flow problems. The algorithm is based on a parallel non-overlapping Schur complement domain-decomposition technique for triangulated domains. In the Schur complement technique, the triangulation is first partitioned into a number of non-overlapping subdomains and interfaces. This suggests a reordering of triangulation vertices which separates subdomain and interface solution unknowns. The reordering induces a natural 2 x 2 block partitioning of the discretization matrix. Exact LU factorization of this block system yields a Schur complement matrix which couples subdomains and the interface together. The remaining sections of this paper present a family of approximate techniques for both constructing and applying the Schur complement as a domain-decomposition preconditioner. The approximate Schur complement serves as an algebraic coarse space operator, thus avoiding the known difficulties associated with the direct formation of a coarse space discretization. In developing Schur complement approximations, particular attention has been given to improving sequential and parallel efficiency of implementations without significantly degrading the quality of the preconditioner. A computer code based on these developments has been tested on the IBM SP2 using MPI message passing protocol. A number of 2-D calculations are presented for both scalar advection-diffusion equations as well as the Euler equations governing compressible fluid flow to demonstrate performance of the preconditioning algorithm.

Barth, Timothy J.↗
NASA NTRS1986

Modeling turbulent, reacting flow

Several of the approximations or models involved in the development of a numerical combustor flow code are examined. In the first section, the importance of numerical accuracy is illustrated, and the impact that improved-accuracy schemes have on slowing convergence is demonstrated. Solution algorithms that can speed convergence are discussed and some performance features of these algorithms are illustrated. A sample calculation displaying the importance of boundary conditions on a three-dimensional numerical prediction is presented. The inaccuracy of a current turbulence model in highly turbulent (nonequilibrium) regions is described. Finally, the surprisingly good performance of a six-flux model in describing radiation heat transfer is displayed. In all the areas examined, continued research is still needed, but valuable engineering tools are available today.

Claus, Russell W.↗
NASA NTRS1988

Computations of vortical supersonic flows over a blunt-nose-cylinder at high angles of attack

Three-dimensional, viscous, and separated flows over a blunt-nose-cylinder at 20, 32, and 44-deg attack angles were computed. The approaching freestream was at a Mach number of 1.6 and a unit Reynolds number of 2 x 10 to the 6th/ft with a total temperature of 585 R. The cylinder used for the computations had a length-to-diameter ratio of 6.67 with a base diameter of 3 in. The flowfield was dominated by large-scale and multiple vortices generated by crossflow separation. The effect of turbulence on the flow structure of one case was modeled algebraically with modifications to correct the length and velocity scales in the regions of separation. The mass averaged Navier-Stokes equations were solved by an approximately factored, upwind-biased, implicit, finite volume scheme. The initialization of the flows was enhanced by a mesh sequencing strategy applied to the diagonalized form of the discretized equations. The convergence to steady-state was accelerated by a multigrid algorithm and using the block inversions for the discretized equations. Calculations were compared with experimental results.

Baysal, O.↗
NASA NTRS1989

A method for solution of the Euler-Bernoulli beam equation in flexible-link robotic systems

An efficient numerical method for solving the partial differential equation (PDE) governing the flexible manipulator control dynamics is presented. A finite-dimensional model of the equation is obtained through discretization in both time and space coordinates by using finite-difference approximations to the PDE. An expert program written in the Macsyma symbolic language is utilized in order to embed the boundary conditions into the program, accounting for a mass carried at the tip of the manipulator. The advantages of the proposed algorithm are many, including the ability to (1) include any distributed actuation term in the partial differential equation, (2) provide distributed sensing of the beam displacement, (3) easily modify the boundary conditions through an expert program, and (4) modify the structure for running under a multiprocessor environment.

Tzes, Anthony P.↗
NASA NTRS1991

Partitioning of unstructured problems for parallel processing

Many large-scale computational problems are based on unstructured computational domains. Primary examples are unstructured grid calculations based on finite volume methods in computational fluid dynamics, or structural analysis problems based on finite element approximations. The question of how to distribute such unstructured computational domains over a large number of processors in a MIMD machine with distributed memory is addressed. A graph theoretical framework for these problems is established. Based on this framework three decomposition algorithms are introduced. In particular a new decomposition algorithm is discussed, which is based on the computation of an eigenvector of the Laplacian matrix associated with the graph. Numerical comparisons on large-scale two- and three-dimensional problems demonstrate the superiority of the new spectral bisection algorithm.

Simon, H. D.↗
NASA NTRS1993

Damage detection using experimentally measured mass and stiffness matrices

A method is presented for locating physical damage or change in a structure using experimentally measured mass and stiffness matrices. The approach uses a recently developed algorithm for transforming a state-space realization into a second order structural model with physical displacements as the generalized coordinates. This is accomplished by first rotating a state-space model of the identified structural dynamics into modal coordinates and approximating the mass normalized modal vectors for the output measurement set. Next, the physical mass, damping and stiffness matrices are synthesized directly from the measured modal parameters. This yields experimental mass and stiffness matrices for the structure without the use of a finite element model or a numerical search. The computed mass and stiffness are asymptotically equivalent to a static condensation of the global physical coordinate model. Techniques for solving the inverse connectivity problem are then developed whereby it is possible to assess the stiffness in a region of the structure bounded by several sensors. Applications to both simulated data and experimental data are used to discuss the effectiveness of the approach.

Peterson, L. D.↗
NASA NTRS1995

Effects of Raman scattering on the water-leaving radiance

The contribution of Raman scattering to the water-leaving radiance is examined using Monte Carlo simulations. Exit angle information is retained, allowing a comparison of different satellite viewing directions. Chlorophyll values of 0.0, 0.01, 0.1, and 1.0 mg Chl/cu m are simulated. Little directional variability is found, with the exception of the direct solar backscatter direction. The wavelength variability is greatest for low chlorophyll concentrations and is negligible for 1.0 mg Chl/cu m. At 550 nm the Raman contribution ranges from approximately 18% of the total water-leaving radiance for pure water to 3% for 1.0 mg Chl/cu m. At 440 nm the range is from 6% to 2%, indicating that Raman scattering will impact radiance ratios for ocean color satellite algorithms.

Waters, Kirk J.↗
NASA NTRS2005

Elastic model for crimped collagen fibrils

A physiologic constitutive expression is presented in algorithmic format for the nonlinear elastic response of wavy collagen fibrils found in soft connective tissues. The model is based on the observation that crimped fibrils in a fascicle have a three-dimensional structure at the micron scale that we approximate as a helical spring. The symmetry of this wave form allows the force/displacement relationship derived from Castigliano's theorem to be solved in closed form: all integrals become analytic. Model predictions are in good agreement with experimental observations for mitral-valve chordae tendinece.

Evaluation Studies↗