Search NASA⌕ Search

SEARCH · Search NASA

Results for “arraying algorithms”

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 631 records · Page 35

On the VHF Source Retrieval Errors Associated with Lightning Mapping Arrays (LMAs)

This presentation examines in detail the standard retrieval method: that of retrieving the (x, y, z, t) parameters of a lightning VHF point source from multiple ground-based Lightning Mapping Array (LMA) time-of-arrival (TOA) observations. The solution is found by minimizing a chi-squared function via the Levenberg-Marquardt algorithm. The associated forward problem is examined to illustrate the importance of signal-to-noise ratio (SNR). Monte Carlo simulated retrievals are used to assess the benefits of changing various LMA network properties. A generalized retrieval method is also introduced that, in addition to TOA data, uses LMA electric field amplitude measurements to retrieve a transient VHF dipole moment source.

inversion↗

An introduction to financial option pricing on a qudit-based quantum computer

The financial sector is anticipated to be one of the first industries to benefit from the increased computational power of quantum computers, in areas such as portfolio optimisation and risk management to financial derivative pricing. Financial mathematics, and derivative pricing in particular, are not areas quantum physicists are traditionally trained in despite the fact that they often have the raw technical skills needed to understand such topics. On the other hand, most quantum algorithms have largely focused on qubits, which are comprised of two discrete states, as the information carriers. However, discrete higher-dimensional qudits, in addition to possibly possessing increased noise robustness and allowing for novel error correction protocols in certain hardware implementations, also have logarithmically greater information storage and processing capacity. In the current NISQ era of quantum computing, a wide array of hardware paradigms are still being studied and any potential advantage a platform offers is worth exploring. Here we introduce the basic concepts behind financial derivatives for the unfamiliar enthusiast as well as outline in great detail the quantum algorithm routines needed to price a European option, the simplest derivative. This is done within the context of a quantum computer comprised of qudits and employing the natural higher-dimensional analogue of a qubit-based pricing algorithm with its various subroutines. From these pieces, one should relatively easily be able to tailor the scheme to more complex, realistic financial derivatives. Finally, the entire stack is numerically simulated with the results demonstrating how the qudit-based scheme's payoff quickly approaches that of both a similarly-resourced classical computer as well as the true payoff, within error, for a modest increase in qudit dimension.

Bornman, Nicholas↗

Fpack and Funpack Utilities for FITS Image Compression and Uncompression

Fpack is a utility program for optimally compressing images in the FITS (Flexible Image Transport System) data format (see http://fits.gsfc.nasa.gov). The associated funpack program restores the compressed image file back to its original state (as long as a lossless compression algorithm is used). These programs may be run from the host operating system command line and are analogous to the gzip and gunzip utility programs except that they are optimized for FITS format images and offer a wider choice of compression algorithms. Fpack stores the compressed image using the FITS tiled image compression convention (see http://fits.gsfc.nasa.gov/fits_registry.html). Under this convention, the image is first divided into a user-configurable grid of rectangular tiles, and then each tile is individually compressed and stored in a variable-length array column in a FITS binary table. By default, fpack usually adopts a row-by-row tiling pattern. The FITS image header keywords remain uncompressed for fast access by FITS reading and writing software. The tiled image compression convention can in principle support any number of different compression algorithms. The fpack and funpack utilities call on routines in the CFITSIO library (http://hesarc.gsfc.nasa.gov/fitsio) to perform the actual compression and uncompression of the FITS images, which currently supports the GZIP, Rice, H-compress, and PLIO IRAF pixel list compression algorithms.

Pence, W.↗

Implementation of Complex Signal Processing Algorithms for Position-Sensitive Microcalorimeters

We have recently reported on a theoretical digital signal-processing algorithm for improved energy and position resolution in position-sensitive, transition-edge sensor (POST) X-ray detectors [Smith et al., Nucl, lnstr and Meth. A 556 (2006) 2371. PoST's consists of one or more transition-edge sensors (TES's) on a large continuous or pixellated X-ray absorber and are under development as an alternative to arrays of single pixel TES's. PoST's provide a means to increase the field-of-view for the fewest number of read-out channels. In this contribution we extend the theoretical correlated energy position optimal filter (CEPOF) algorithm (originally developed for 2-TES continuous absorber PoST's) to investigate the practical implementation on multi-pixel single TES PoST's or Hydras. We use numerically simulated data for a nine absorber device, which includes realistic detector noise, to demonstrate an iterative scheme that enables convergence on the correct photon absorption position and energy without any a priori assumptions. The position sensitivity of the CEPOF implemented on simulated data agrees very well with the theoretically predicted resolution. We discuss practical issues such as the impact of random arrival phase of the measured data on the performance of the CEPOF. The CEPOF algorithm demonstrates that full-width-at- half-maximum energy resolution of < 8 eV coupled with position-sensitivity down to a few 100 eV should be achievable for a fully optimized device.

Smith, Stephen J.↗

Machine-Learning Identification of the Sensing Descriptors Relevant in Molecular Interactions with Metal Nanoparticle-Decorated Nanotube Field-Effect Transistors

We report that carbon nanotube-based field-effect transistors (NTFETs) are ideal sensor devices as they provide rich information regarding carbon nanotube interactions with target analytes and have potential for miniaturization in diverse applications in medical, safety, environmental, and energy sectors. Herein, we investigate chemical detection with cross-sensitive NTFETs sensor arrays comprised of metal nanoparticle-decorated single-walled carbon nanotubes (SWCNTs). By combining analysis of NTFET device characteristics with supervised machine-learning algorithms, we have successfully discriminated among five selected purine compounds, adenine, guanine, xanthine, uric acid, and caffeine. Interactions of purine compounds with metal nanoparticle-decorated SWCNTs were corroborated by density functional theory calculations. Furthermore, by testing a variety of prepared as well as commercial solutions with and without caffeine, our approach accurately discerns the presence of caffeine in 95% of the samples with 48 features using a linear discriminant analysis and in 93.4% of the samples with only 11 features when using a support vector machine analysis. Lastly, we also performed recursive feature elimination and identified three NTFET parameters, transconductance, threshold voltage, and minimum conductance, as the most crucial features to analyte prediction accuracy.

36 MATERIALS SCIENCE↗

PyDA: A Hands-On Introduction to Dynamical Data Assimilation with Python

Dynamic data assimilation offers a suite of algorithms that merge measurement data with numerical simulations to predict accurate state trajectories. Meteorological centers rely heavily on data assimilation to achieve trustworthy weather forecast. With the advance in measurement systems, as well as the reduction in sensor prices, data assimilation (DA) techniques are applicable to various fields, other than meteorology. However, beginners usually face hardships digesting the core ideas from the available sophisticated resources requiring a steep learning curve. In this tutorial, we lay out the mathematical principles behind DA with easy-to-follow Python module implementations so that this group of newcomers can quickly feel the essence of DA algorithms. We explore a series of common variational, and sequential techniques, and highlight major differences and potential extensions. We demonstrate the presented approaches using an array of fluid flow applications with varying levels of complexity.

97 MATHEMATICS AND COMPUTING↗

A Search for Pulsars around Sgr A* in the First Event Horizon Telescope Data Set

In 2017 the Event Horizon Telescope (EHT) observed the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz (λ = 1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT data sets. The high observing frequency means that pulsars—which typically exhibit steep emission spectra—are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope, and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the fast folding algorithm, and single-pulse searches targeting both pulsars and burst-like transient emission. We use the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction (≲2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Sequence information signal processor for local and global string comparisons

A sequence information signal processing integrated circuit chip designed to perform high speed calculation of a dynamic programming algorithm based upon the algorithm defined by Waterman and Smith. The signal processing chip of the present invention is designed to be a building block of a linear systolic array, the performance of which can be increased by connecting additional sequence information signal processing chips to the array. The chip provides a high speed, low cost linear array processor that can locate highly similar global sequences or segments thereof such as contiguous subsequences from two different DNA or protein sequences. The chip is implemented in a preferred embodiment using CMOS VLSI technology to provide the equivalent of about 400,000 transistors or 100,000 gates. Each chip provides 16 processing elements, and is designed to provide 16 bit, two's compliment operation for maximum score precision of between -32,768 and +32,767. It is designed to provide a comparison between sequences as long as 4,194,304 elements without external software and between sequences of unlimited numbers of elements with the aid of external software. Each sequence can be assigned different deletion and insertion weight functions. Each processor is provided with a similarity measure device which is independently variable. Thus, each processor can contribute to maximum value score calculation using a different similarity measure.

Peterson, John C.↗

I’ve Looked at Clouds from Both Sides, Now: Viewpoints from Surface and Spaceborne Lidar Systems

Clouds play a critical role in the Earth’s climate system because they are inextricably linked to the hydrological cycle and radiation budget. Information about cloud height, thickness, occurrence, and amount are critical inputs for a host of numerical applications involving climate research. Therefore, it is important to have highly accurate and quantitative data records of cloud properties that span several years and geographic regions. Verification of even the most basic modeling processes demands long term and continuous observations of global cloud occurrence, if there is to be any confidence in their fidelity. For a number of reasons, however, it is impossible to gauge the complexities of clouds from a single source. Fundamentally, an array of remote sensing methods is needed in order to provide a complete picture. Within the Micropulse Lidar Network (MPLNET), we have developed a new algorithm to improve the quality of our cloud products. The largest impact of the changes to the cloud detection algorithm is evident with high clouds (those with cloud base > 5 km). Furthermore, polarized measurements enable us to determine the cloud thermodynamic phase. Naturally, the next step is to investigate how these improvements compare with observations from spaceborne lidars (e.g. CALIOP) which have a better (unobstructed) view of high clouds. Recent advances in remote sensing have revealed that cirrus clouds are the most common cloud genus observed in the atmosphere. Furthermore, cirrus skew highly towards relatively low cloud optical depths, as observed from both surface and spaceborne viewpoints. The radiative impacts of these findings are quite significant, considering the cumulative effect cirrus exhibit when compared to low clouds.

Jasper R Lewis↗

The OMPS Limb Profiler Instrument: An Alternative Data Analysis and Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth fs limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. The focal plane of the LP spectrometer is a two ]dimensional CCD array comprised of 340 x 740 pixels. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. The paper will describe an alternative algorithm which will retrieve ozone density and aerosol extinction directly from radiance data collected on individual CCD pixels. This alternative method uses an optimal estimation approach to retrieve ozone and aerosol in the 10-60 km range from the information contained within an ensemble of about 50000 down-linked pixels. Tangent height registration is performed using the Rayleigh Scattering Attitude Sensor (RSAS) technique applied to columns of pixels in the 340-360 nm range. Cloud height is determined by analyzing the radiance first derivative along pixel columns at longer wavelengths. Wavelength registration is performed using rows of pixels and identifying Fraunhofer solar lines within the measured spectra. Special attention is given to stray-light decontamination and modeling of the measured finite spectral/spatial line shape functions.

Rault, Didier F.↗

Quantum Adiabatic Optimization with Rydberg Arrays: Localization Phenomena and Encoding Strategies

Quantum adiabatic optimization seeks to solve combinatorial problems using quantum dynamics, requiring the Hamiltonian of the system to align with the problem of interest. However, these Hamiltonians are often incompatible with the native constraints of quantum hardware, necessitating encoding strategies to map the original problem into a hardware-conformant form. While the classical overhead associated with such mappings is easily quantifiable and typically polynomial in problem size, it is much harder to quantify their overhead on the quantum algorithm, e.g., the transformation of the adiabatic timescale. In this work, we address this challenge on the concrete example of the encoding scheme proposed in [Nguyen , PRX Quantum , 010316 (2023)], which is designed to map optimization problems on arbitrarily connected graphs into Rydberg atom arrays. We consider the fundamental building blocks underlying this encoding scheme and determine the scaling of the minimum gap with system size along adiabatic protocols. Even when the original problem is trivially solvable, we find that the encoded problem can exhibit an exponentially closing minimum gap. We show that this originates from a quantum coherent effect, which gives rise to an unfavorable localization of the ground-state wave function. On the QuEra Aquila neutral atom machine, we observe such localization and its effect on the success probability of finding the correct solution to the encoded optimization problem. Finally, we propose quantum-aware modifications of the encoding scheme that avoid this quantum bottleneck and lead to an exponential improvement in the adiabatic performance. This highlights the crucial importance of accounting for quantum effects when designing strategies to encode classical problems onto quantum platforms. Published by the American Physical Society 2025

Bombieri, Lisa (ORCID:0009000950422897)↗

Muon track reconstruction in a segmented bolometric array using multi-objective optimization

Recent advances in segmented solid-state detector arrays for rare-event searches have allowed the technology to approach the ton-scale in detector mass and the scale of meters in size. Often focused around searches for neutrinoless double-beta decay or direct dark matter detection, such experiments also have the capability to search for exotic particles that leave track-like signatures across their volume. However, the segmented nature of such detector arrays often sets the spatial resolution and makes the problem of reconstructing track-like paths non-trivial. Here, in this paper, we present an algorithm that improves reconstruction of track-like events in segmented detectors using multi-objective optimization — a computational technique that optimizes more than one cost function at a time without specifying a quantitative weighting between them. Such a technique allows the reconstruction of tracks through a detector and the determination of path-lengths through individual elements. When combined with the reconstructed energy depositions in each element this allows for a calculation of the stopping power of track-like particles and opens the door to searches for particles with abnormal stopping power like monopoles or lightly-ionizing particles (LIPs). Results are presented which evaluate the precision of the reconstruction tools as they currently stand against Monte Carlo generated data. The algorithm is presented in the context of the CUORE experiment, but has applications to other segmented calorimeter detectors.

47 OTHER INSTRUMENTATION↗

Active isolation of vibrations with adaptive structures

Vibration transmission in structures is controlled by means of a technique which employs distributed arrays of piezoelectric transducers bonded to the supporting structure. Distributed PVDF piezoelectric strips are employed as error sensors, and a two-channel feedforward adaptive LMS algorithm is used for minimizing error signals and thereby controlling the structure. A harmonic force input excites a thick plate, and a receiving plate is configured with three pairs of piezoelectric actuators. Modal analyses are performed to determine the resonant frequencies of the system, and a scanning laser vibrometer is used to study the shape of the response of the receiving plate during excitation with and without the control algorithm. Efficient active isolation of the vibrations is achieved with modal suppression, and good control is noted in the on-resonance cases in which increased numbers of PVDF sensors and piezoelectric actuators are employed.

Guigou, C.↗

MOLA II Laser Transmitter Calibration and Performance

The goal of the document is to explain the algorithm for determining the laser output energy from the telemetry data within the return packets from MOLA II. A simple algorithm is developed to convert the raw start detector data into laser energy, measured in millijoules. This conversion is dependent on three variables, start detector counts, array heat sink temperature and start detector temperature. All these values are contained within the return packets. The conversion is applied to the GSFC Thermal Vacuum data as well as the in-space data to date and shows good correlation.

Afzal, Robert S.↗

Development of an Aerosol Opacity Retrieval Algorithm for Use with Multi-Angle Land Surface Images

In 1998, the Multi-angle Imaging SpectroRadiometer (MISR) will fly aboard the EOS-AM1 spacecraft. MISR will enable unique methods for retrieving the properties of atmospheric aerosols, by providing global imagery of the Earth at nine viewing angles in four visible and near-IR spectral bands. As part of the MISR algorithm development, theoretical methods of analyzing multi-angle, multi-spectral data are being tested using images acquired by the airborne Advanced Solid-State Array Spectroradiometer (ASAS). In this paper we derive a method to be used over land surfaces for retrieving the change in opacity between spectral bands, which can then be used in conjunction with an aerosol model to derive a bound on absolute opacity.

Advanced Solid-State Array Spectroradiometer ASAS ↗

Shear Stress Sensing with Elastic Microfence Structures

In this work, elastic microfences were generated for the purpose of measuring shear forces acting on a wind tunnel model. The microfences were fabricated in a two part process involving laser ablation patterning to generate a template in a polymer film followed by soft lithography with a two-part silicone. Incorporation of a fluorescent dye was demonstrated as a method to enhance contrast between the sensing elements and the substrate. Sensing elements consisted of multiple microfences prepared at different orientations to enable determination of both shear force and directionality. Microfence arrays were integrated into an optical microscope with sub-micrometer resolution. Initial experiments were conducted on a flat plate wind tunnel model. Both image stabilization algorithms and digital image correlation were utilized to determine the amount of fence deflection as a result of airflow. Initial free jet experiments indicated that the microfences could be readily displaced and this displacement was recorded through the microscope.

Cisotto, Alexxandra↗

Genetic algorithm optimization of tin–copper graded shielding for improved plutonium safeguards measurements

Plutonium nuclear material safeguards is essential to ensure that states do not divert plutonium to clandestine nuclear weapons programs. Traditional safeguards methods rely on correlated neutron measurements. Detection systems based on organic scintillators, which are sensitive to fast neutrons and photons, are gaining support for use in safeguards applications due to increased capabilities and shorter measurement times compared to the current state-of-the-art. A drawback of organic scintillators in plutonium measurements is data throughput limitations; 241 Am builds up in plutonium samples and emits 60-keV gamma rays with a high specific activity. In this work, a genetic algorithm is used to design graded shielding made of tin and copper to attenuate the 60-keV gamma rays. The shield is designed to limit the attenuation of the higher-energy, fission photons, because fission photons can be used as a signature. The genetic algorithm optimization takes days to run, while a complete search space exploration would take years. The genetic algorithm was run for different levels of 60-keV gamma-ray attenuation, and the thickness of the required shielding layers was determined at each level. Tin–copper shielding was placed in front of an array of organic scintillation detectors for a measurement of approximately 4.5 kg of alpha-phase, weapons-grade plutonium, an AmLi source, and a 252 Cf source. We show that the shield successfully attenuated low-energy, non-fission gamma rays, which significantly improved the organic scintillator performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Digital Control Algorithm for Magnetic Suspension Systems

An ongoing program exists to investigate and develop magnetic suspension technologies and modelling techniques at NASA Langley Research Center. Presently, there is a laboratory-scale large air-gap suspension system capable of five degree-of-freedom (DOF) control that is operational and a six DOF system that is under development. Those systems levitate a cylindrical element containing a permanent magnet core above a planar array of electromagnets, which are used for levitation and control purposes. In order to evaluate various control approaches with those systems, the Generic Real-Time State-Space Controller (GRTSSC) software package was developed. That control software package allows the user to implement multiple control methods and allows for varied input/output commands. The development of the control algorithm is presented. The desired functionality of the software is discussed, including the ability to inject noise on sensor inputs and/or actuator outputs. Various limitations, common issues, and trade-offs are discussed including data format precision; the drawbacks of using either Direct Memory Access (DMA), interrupts, or program control techniques for data acquisition; and platform dependent concerns related to the portability of the software, such as memory addressing formats. Efforts to minimize overall controller loop-rate and a comparison of achievable controller sample rates are discussed. The implementation of a modular code structure is presented. The format for the controller input data file and the noise information file is presented. Controller input vector information is available for post-processing by mathematical analysis software such as MATLAB1.

Britton, Thomas C.↗