An algorithm for interpolation and numerical differentiation
Easily automated technique for computing continuous function to approximate discrete set of data
Engvall, J. L.↗
SEARCH · Search NASA
Results for “algorithmic differentiation”
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.
Simulation and analysis of differential global positioning system for civil helicopter operations
A Differential Global Positioning System (DGPS) computer simulation was developed, to provide a versatile tool for assessing DGPS referenced civil helicopter navigation. The civil helicopter community will probably be an early user of the GPS capability because of the unique mission requirements which include offshore exploration and low altitude transport into remote areas not currently served by ground based Navaids. The Monte Carlo simulation provided a sufficiently high fidelity dynamic motion and propagation environment to enable accurate comparisons of alternative differential GPS implementations and navigation filter tradeoffs. The analyst has provided the capability to adjust most aspects of the system, the helicopter flight profile, the receiver Kalman filter, and the signal propagation environment to assess differential GPS performance and parameter sensitivities. Preliminary analysis was conducted to evaluate alternative implementations of the differential navigation algorithm in both the position and measurement domain. Results are presented to show that significant performance gains are achieved when compared with conventional GPS but that differences due to DGPS implementation techniques were small. System performance was relatively insensitive to the update rates of the error correction information.
Denaro, R. P.↗
Accuracy study of the Upper Atmosphere Research Satellite (UARS) definitive attitude determination
The Upper Atmosphere Research Satellite (UARS) has two definitive attitude determination requirements: the definitive attitude of the Modular Attitude Control Subsystem (MACS) and the definitive attitude of the gimbaled Solar-Stellar Pointing Platform (SSPP). The onboard computer (OBC) will compute the MACS attitude using a Kalman filter and will transform this attitude solution through the SSPP gimbals to calculate the SSPP attitude. The attitude ground support system (AGSS) will compute the MACS attitude using a batch least-squares differential corrector algorithm and will also transform this solution through the gimbals to obtain the SSPP attitude. This paper reports the results of a prelaunch study to predict the accuracy of the OBC attitude solutions and the accuracy of the AGSS attitude solutions. The OBC and AGSS solution accuracies are then compared to establish the relative quality. The effects of star observability, sensor noise, and sensor misalignment uncertainties on attitude determination accuracy are analyzed for each case.
Snow, Frank↗
Automation of orbit determination functions for National Aeronautics and Space Administration (NASA)-supported satellite missions
The Flight Dynamics Facility (FDF) at Goddard Space Flight Center (GSFC) provides spacecraft trajectory determination for a wide variety of National Aeronautics and Space Administration (NASA)-supported satellite missions, using the Tracking Data Relay Satellite System (TDRSS) and Ground Spaceflight and Tracking Data Network (GSTDN). To take advantage of computerized decision making processes that can be used in spacecraft navigation, the Orbit Determination Automation System (ODAS) was designed, developed, and implemented as a prototype system to automate orbit determination (OD) and orbit quality assurance (QA) functions performed by orbit operations. Based on a machine-resident generic schedule and predetermined mission-dependent QA criteria, ODAS autonomously activates an interface with the existing trajectory determination system using a batch least-squares differential correction algorithm to perform the basic OD functions. The computational parameters determined during the OD are processed to make computerized decisions regarding QA, and a controlled recovery process isactivated when the criteria are not satisfied. The complete cycle is autonomous and continuous. ODAS was extensively tested for performance under conditions resembling actual operational conditions and found to be effective and reliable for extended autonomous OD. Details of the system structure and function are discussed, and test results are presented.
Mardirossian, H.↗
Automation of orbit determination functions for National Aeronautics and Space Administration (NASA)-supported satellite missions
The Flight Dynamics Facility (FDF) at Goddard Space Flight Center (GSFC) provides spacecraft trajectory determination for a wide variety of National Aeronautics and Space Administration (NASA)-supported satellite missions, using the Tracking Data Relay Satellite System (TDRSS) and Ground Spaceflight and Tracking Data Network (GSTDN). To take advantage of computerized decision making processes that can be used in spacecraft navigation, the Orbit Determination Automation System (ODAS) was designed, developed, and implemented as a prototype system to automate orbit determination (OD) and orbit quality assurance (QA) functions performed by orbit operations. Based on a machine-resident generic schedule and predetermined mission-dependent QA criteria, ODAS autonomously activates an interface with the existing trajectory determination system using a batch least-squares differential correction algorithm to perform the basic OD functions. The computational parameters determined during the OD are processed to make computerized decisions regarding QA, and a controlled recovery process is activated when the criteria are not satisfied. The complete cycle is autonomous and continuous. ODAS was extensively tested for performance under conditions resembling actual operational conditions and found to be effective and reliable for extended autonomous OD. Details of the system structure and function are discussed, and test results are presented.
Mardirossian, H.↗
Evolutionary optimization of a Genetically Refined Truss
Structural optimization is a field of research that has experienced noteworthy growth for many years. Researchers in this area have developed optimization tools to successfully design and model structures, typically minimizing mass while maintaining certain deflection and stress constraints. Numerous optimization studies have been performed to minimize mass, deflection and stress on a benchmark cantilever truss problem. Predominantly traditional optimization theory is applied to this problem. The cross-sectional area of each member is optimized to minimize the aforementioned objectives. This paper will present a structural optimization technique that has been previously applied to compliant mechanism design. This technique demonstrates a method that combines topology optimization, geometric refinement, finite element analysis, and two forms of evolutionary computation: Genetic Algorithms and Differential Evolution to successfully optimize a benchmark structural optimization problem. An non-traditional solution to the benchmark problem is presented in this paper, specifically a geometrically refined topological solution. The design process begins with an alternate control mesh formulation, multilevel geometric smoothing operation, and an elastostatic structural analysis. The design process is wrapped in an evolutionary computing optimization toolset.
Hull, Patrick V.↗
Evolutionary Optimization of a Geometrically Refined Truss
Structural optimization is a field of research that has experienced noteworthy growth for many years. Researchers in this area have developed optimization tools to successfully design and model structures, typically minimizing mass while maintaining certain deflection and stress constraints. Numerous optimization studies have been performed to minimize mass, deflection, and stress on a benchmark cantilever truss problem. Predominantly traditional optimization theory is applied to this problem. The cross-sectional area of each member is optimized to minimize the aforementioned objectives. This Technical Publication (TP) presents a structural optimization technique that has been previously applied to compliant mechanism design. This technique demonstrates a method that combines topology optimization, geometric refinement, finite element analysis, and two forms of evolutionary computation: genetic algorithms and differential evolution to successfully optimize a benchmark structural optimization problem. A nontraditional solution to the benchmark problem is presented in this TP, specifically a geometrically refined topological solution. The design process begins with an alternate control mesh formulation, multilevel geometric smoothing operation, and an elastostatic structural analysis. The design process is wrapped in an evolutionary computing optimization toolset.
Hull, P. V.↗
Evolutionary Computational Methods for Identifying Emergent Behavior in Autonomous Systems
A technique based on Evolutionary Computational Methods (ECMs) was developed that allows for the automated optimization of complex computationally modeled systems, such as autonomous systems. The primary technology, which enables the ECM to find optimal solutions in complex search spaces, derives from evolutionary algorithms such as the genetic algorithm and differential evolution. These methods are based on biological processes, particularly genetics, and define an iterative process that evolves parameter sets into an optimum. Evolutionary computation is a method that operates on a population of existing computational-based engineering models (or simulators) and competes them using biologically inspired genetic operators on large parallel cluster computers. The result is the ability to automatically find design optimizations and trades, and thereby greatly amplify the role of the system engineer.
Terrile, Richard J.↗
Differences in the Vertical and Microphysical Evolution of Volcanic and Pyrocumulonimbus Stratospheric Aerosol Plumes as Observed by CALIOP and CATS Satellite Lidar
For some time, volcanic eruptions have been thought to be the only significant direct injectors of aerosols in the stratosphere. However, recent fire seasons have featured fire events resulting in large volcanic-sized pyrocumulonimbus plumes of smoke aerosols reaching many kilometers into the lower stratosphere. To understand and model the effects of these pyrocumulonimbus events on stratospheric composition and climate, a natural analogy lies with better studied volcanic events; however, differences in plume composition may limit this comparison. Using satellite lidar from both CATS and CALIOP, we show that the stratospheric aerosol plumes from the record-setting Pacific Northwest pyrocumulonimbus event of 2017 and the Calbuco volcanic eruption of 2015 evolve differently both vertically and microphysically. Specifically, depolarization ratios indicate that this pyrocumulonimbus event’s aerosol particles became more irregularly shaped over time in contrast to volcanic aerosols which become more spherical over time. Accounting for these changes in aerosol properties may be significant in assessing the effects of pyrocumulonimbus events on the Earth’s radiative balance and aid in refining stratospheric aerosol typing algorithms to differentiate volcanic from pyrocumulonimbus plumes.
Christian, Kenneth↗
Evaluation of a Multizone Impedance Eduction Method
A computational study is used to evaluate the PyCHE impedance eduction method developed at the NASA Langley Research Center. This method combines an aeroacoustic duct propagation code based on numerical solution to the convected Helmholtz equation with a global optimizer that uses the Differential Evolution algorithm. The efficacy of this method is evaluated with acoustic pressure data simulated to represent that measured with one-zone, two-zone, and three-zone liners mounted in the NASA Langley Grazing Flow Impedance Tube. The PyCHE method has a normalized impedance error of approximately 0.2 for (uniform) one-zone liners with a length of at least 5”, and produces quite reasonable results for liners as short as 2”. Whereas the impedance of the liner has an effect on eduction accuracy, the amount of attenuation is shown to be the dominant parameter. Similar results are observed for two-zone liners, for which the impedance of each zone is unique. The two-zone results also indicate it is more difficult to accurately educe resistance than reactance, and a zone length of at least 6” (slightly longer than for uniform liners) is needed to limit the normalized error to 0.2. The PyCHE method is also demonstrated to successfully educe the impedances for each zone of a three-zone liner. These results are sufficiently encouraging to warrant the continued usage of the PyCHE impedance eduction method for single and multizone liners.
Jones, M. G.↗
Development and Analysis of a Thick Cloud Layers Database for Lightning Launch Commit Criteria Improvement
Lightning can pose a potential threat to space launch vehicles. In response to this, rules were created called the Lightning Launch Commit Criteria (LLCC) that help weather personnel evaluate the potential for natural and rocket-triggered lightning. One of the ten LLCC with the least research is called the Thick Cloud Layers rule. To further understand electrification of thick cloud layers and potentially improve the Thick Cloud Layers rule, a database of thick cloud layers that occurred over the Eastern Range was created. This database is then used to create an algorithm for identifying and differentiating thick cloud layers from other cloud types based on radar characteristics, temperature levels in reference to cloud height, and the surface electric field. By analyzing and identifying thick cloud events, this project could help narrow down when thick clouds are occurring and potentially minimize unnecessary launch delays. Events that caused LLCC violations involving the Thick Cloud Layers rule were analyzed by hand using Level-2 NEXRAD radar data from the National Weather Service WSR-88D radar in Melbourne with the program GR2Analyst. Cases that were found to be isolated and not involved with convection were recorded (date, start/end time, location) in a database. Radar data associated with these cases was collected and gridded using Python radar packages. Once gridded, I calculated and recorded for each radar scan the following radar reflectivity driven parameters within an 11x11 km bin centered on each 1 square km grid point: the mean reflectivity colder than 0 degrees Celsius, Maximum Radar Reflectivity (MRR) colder than 0 degrees Celsius, Volume Averaged Height Integrated Radar Reflectivity (VAHIRR), Hydrometeor Identification (HID), the difference between the maximum and mean reflectivity, the cloud depth colder than 0 degrees Celsius, the overall cloud depth, the cloud top, and the cloud bottom. Soundings for each event were used to determine cloud temperature levels, and where the cloud is in relation to the freezing level. Electric field mill data collected over the Eastern Range was used to determine surface electric fields below each cloud. All parameters were analyzed in depth for several thick cloud cases to gain an understanding of typical thick cloud characteristics. Cases of thick clouds and other isolated cloud types were also recorded for training purposes to see if enough differences exist between cloud types to differentiate them with an algorithm. Each case along with its corresponding characteristics was recorded in a database, and this database was used to compare differing cloud types, as well as train the algorithm to detect thick clouds.
Lightning↗
The coverage of elliptical orbits using Ergodic theory
One of the key performance metrics for satellite constellations is the statistics of the visibility periods between the satellites and poins on the ground. Associated with this are other desirable communications statistis such as data through-put, link qualities, etc. Typically, the computation of coverage statistics requires the propagation of the trajectories. In this paper a new algorithm using differential geometry enables us to extend the Ergodic theory to eliptical orbits.
Lo, Martin W.↗
Tropospheric Ozone Retrieval By A Combination of TROPOMI/S5P Measurements With BASCOE Assimilated Data
We present a new tropospheric ozone dataset based on TROPOspheric Monitoring Instrument (TROPOMI)/Sentinel-5 Precursor (S5P) total ozone measurements combined with stratospheric ozone data from the Belgian Assimilation System for Chemical ObsErvations (BASCOE) constrained by assimilating ozone observations from the Microwave Limb Sounder (MLS). The BASCOE stratospheric data are interpolated to the S5P observations and subtracted from the TROPOMI total ozone data. The difference is equal to the tropospheric ozone residual column from the surface up to the tropopause. The tropospheric ozone columns are retrieved at the full spatial resolution of the TROPOMI sensor (5.5×3.5 km 2 ) with daily global coverage. Compared to the Ozone Mapping and Profiler Suite Modern-Era Retrospective analysis for Research and Applications 2 (OMPS-MERRA-2) data, a global mean positive bias of 3.3 DU is found for the analysed period April 2018 to June 2020. A small negative bias of about −0.91 DU is observed in the tropics relative to the operational TROPOMI tropical tropospheric data based on the convective cloud differential (CCD) algorithm throughout the same period. The new tropospheric ozone data (S5P-BASCOE) are compared to a set of globally distributed ozonesonde data integrated up to the tropopause level. We found 2254 comparisons with cloud-free TROPOMI observations within 25 km of the stations. In the global mean, S5P-BASCOE deviates by 2.6 DU from the integrated ozonesondes. Depending on the latitude the S5P-BASCOE deviate from the sondes and between −4.8 and 7.9 DU, indicating a good agreement. However, some exceptional larger positive deviations up to 12 DU are found, especially in the northern polar regions (north of 70∘). The monthly mean tropospheric column and time series for selected areas showed the expected spatial and temporal pattern, such as the wave one structure in the tropics or the seasonal cycle, including a summer maximum, in the mid-latitudes.
Satellite Retrieval↗
Bandelier Ecological Conservation: Mapping Invasive Species Along the Rio Grande Corridor in Bandelier National Monument
The Southwest U.S. has experienced a growth of invasive riparian species, specifically Elaeagnus angustifolia (Russian olive), Tamarix ramosissima (saltcedar), and Ulmus pumila (Siberian elm), which alter local soil chemistry and outcompete native species. Locating these exotic species is critical for ecological conservation; however, field identification can be resource intensive. NASA DEVELOP partnered with the National Park Service (NPS) at Bandelier National Monument (BAND) to assess the feasibility of using Earth observation data to map invasive species along the Rio Grande corridor of the park. The team used Landsat 8 OLI, Sentinel-2 MSI, and ISS DESIS imagery to compute principal components based on spectral bands, vegetation indices, and terrain indices. Using the first five principal components, the team created classification maps using both a k-means classification algorithm and a random forest algorithm to differentiate between native and non-native species. The team derived maps for the three invasive riparian species in the region for the last five years. The team found that invasive species covered 33% of the park's river corridor in 2023, and the invasive species extent has increased by 5.7% from 2019 to 2023. The methods will serve as a guide for aiding historic and present invasive species identification in riparian regions, and the NPS staff at BAND will use the results to inform local mitigation practices and advocate for invasive species removal.
Evan Barrett↗
Active-Interrogation Measurements of Induced-Fission Neutrons from Low-Enriched Uranium
Protection and control of nuclear fuels is paramount for nuclear security and safeguards; therefore, it is important to develop fast and robust controlling mechanisms to ensure the safety of nuclear fuels. Through both passive- and active-interrogation methods we can use fast-neutron detection to perform real-time measurements of fission neutrons for process monitoring. Active interrogation allows us to use different ranges of incident neutron energy to probe for different isotopes of uranium. With fast-neutron detectors, such as organic liquid scintillation detectors, we can detect the induced-fission neutrons and photons and work towards quantifying a sample’s mass and enrichment. Using MCNPX-PoliMi, a system was designed to measure induced-fission neutrons from U-235 and U-238. Measurements were then performed in the summer of 2010 at the Joint Research Centre in Ispra, Italy. Fissions were induced with an associated particle D-T generator and an isotopic Am-Li source. The fission neutrons, as well as neutrons from (n, 2n) and (n, 3n) reactions, were measured with five 5” by 5” EJ-309 organic liquid scintillators. The D-T neutron generator was available as part of a measurement campaign in place by Padova University. The measurement and data-acquisition systems were developed at the University of Michigan utilizing a CAEN V1720 digitizer and pulse-shape discrimination algorithms to differentiate neutron and photon detections. Low-enriched uranium samples of varying mass and enrichment were interrogated. Acquired time-of-flight curves and cross-correlation curves are currently analyzed to draw relationships between detected neutrons and sample mass and enrichment. In the full paper, the promise of active-interrogation measurements and fast-neutron detection will be assessed through the example of this proof-of-concept measurement campaign. Additionally, MCNPX-PoliMi simulation results will be compared to the measured data to validate the MCNPX-PoliMi code when used for active-interrogation simulations.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Developing a Drilling Optimization System for Improved Overall Rate of Penetration in Geothermal Wells
Geothermal energy is renewable, reliable and environmentally friendly source of energy. The major cost in the development of geothermal wells is the actual drilling of the wells. The main objective of this paper is to introduce a new real-time drilling optimization system designed for granite formation to reduce the overall drilling cost. In this study, a drilling optimization system is verified using drilling data from Utah-Forge well 58-32. The drilling optimization system used the Utah-Forge well 58-32 data to achieve real-time unconfined compressive strength (UCS). Based on the UCS value from the previous feet, the system simulates the ROP for the next drilling feet. The drilling optimization system utilizes the Differential Evolution Algorithm (DEA), which is a metaheuristic method to search the space of solution, to find the best operating parameters (i.e. WOB and RPM) for the next drilling foot. The optimization algorithm takes a maximum cutter temperature into account as a constraint and avoids the accelerated wear. The developed drilling optimization system improves ROP responses and reduces the drilling cost of geothermal wells. The simulated ROP results from the system show a good agreement with the ROP from Utah-Forge well 58-32 drilling data. The drilling time before and after optimization for both intervals were presented.
15 GEOTHERMAL ENERGY↗
Application of Polynomial Chaos Expansion in One-Dimensional Inverse Transport Problems
In the problem of inverse radiation transport, measurements of particle leakages from radioactive source/shield systems are used to infer unknown parameters within the systems. This reconstruction can be accomplished by finding the physical parameters of the unknown system that minimize the difference between calculated detector responses and measured detector responses. The inverse transport solver should also propagate uncertainties from the detector responses to the reconstructed parameter values. Recently, the Differential Evolution Adaptive Metropolis (DREAM) method was shown to be a robust method for uncertainty quantification in inverse problems. The DREAM method uses the differential evolution algorithm to choose efficient steps in a Markov chain Monte Carlo (MCMC) simulation. Although the DREAM method has been shown to be more efficient than traditional MCMC approaches, it still requires thousands of transport computations to accurately quantify uncertainty. For the simple ray-tracing model, it was possible to perform these calculations in a reasonable amount of time (under an hour). However, inverse transport problems that require computationally-intensive transport solvers will need prohibitively long run times when using the DREAM method. A way to overcome the computational burden of the DREAM method is to perform surrogate modeling using the polynomial chaos expansion (PCE) approach. It was recently shown that the PCE method provides both accurate and computationally inexpensive uncertainty quantification in forward radiation source/shield problems. PCE was also recently applied to an inverse problem in nonlinear point reactor kinetics within a standard MCMC approach. PCE was first used with the DREAM method where it was applied to inverse problems in groundwater modeling.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Radial differential interior tomography and its image reconstruction with differentiated backprojection and projection onto convex sets
Purpose: Interior tomography has been recognized as one of the most effective approaches in computed tomography (CT) to reduce radiation dose rendered to patients. In this work, the authors propose and evaluate an imaging method of radial differential interior tomography.Methods: In interior tomography, an x-ray beam is collimated to only irradiate the region of interest (ROI) with suspected lesions while the surrounding area/volume of normal tissues/organs is spared. In the proposed imaging method of radial differential interior tomography, the outcome is a ROI image that has gone through a radial differential filtering. The image reconstruction algorithm for the radial differential interior tomography is kept in the fashion of differentiated backprojection and projection onto convex sets, but the required a priori knowledge in a small round area becomes zero and may be more readily available in practice.Results: Using the projection data simulated by computer and acquired by CT scanner, the authors evaluate and verify the performance of the proposed radial differential interior tomography method and its associated image reconstruction algorithm. The preliminary results show that the proposed imaging method can generate an image that is the radial differentiation of a conventional tomographic image and is robust over noise that inevitably exist in practice.Conclusions: It is believed that the proposed imaging method may find its utility in advanced clinical applications wherein a ROI-based image processing and analysis is required for lesion visualization, characterization, and diagnosis.
62 RADIOLOGY AND NUCLEAR MEDICINE↗