Engineering topics

Das, Santanu

Publications and source records attributed to Das, Santanu.

12 records

DOE OSTI2022

The Tianlai dish array low-z surveys forecasts

Here, we present the science case for surveys with the Tianlai dish array interferometer tuned to the [1300, 1400] MHz frequency range. Starting from a realistic generation of mock visibility data according to the survey strategy, we reconstruct maps of the sky and perform foreground subtraction. We estimate the level of residuals from imperfect subtraction, mostly due to mode mixing, i.e. distortions in the reconstructed 3D maps due to frequency-dependent instrument response. We show that a survey of the North Celestial Polar cap during a year of observations, covering an area of 150 deg 2 , would reach a sensitivity of $1.5-2 mK per 1 MHz x 0.25 2 deg 2 voxel and be marginally impacted by mode mixing. Tianlai would be able to detect ~10 nearby massive H i clumps as well as a very strong cross-correlation signal of 21 cm intensity maps with the North Celestial Cap Survey optical galaxies. We also studied the performance of a mid-latitude survey, covering ~ 1500 deg 2 overlapping the SDSS footprint. Despite a higher noise level for the mid-latitude survey, as well as significant distortions due to mode mixing, Tianlai would be able to detect a highly significant cross-correlation between the 21 cm signal and the Sloan spectroscopic galaxy sample. Using the extragalactic signals measured from either or both of these surveys, and comparing them with simulations such as those presented here will make it possible to assess the impact of various instrumental imperfections on the Tianlai dish array performance. This would pave the way for future intensity mapping surveys with higher sensitivity.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

AlgoSCR: an algorithm for solar contamination removal from radio interferometric data

Hydrogen intensity mapping is a new field in astronomy that promises to make three-dimensional maps of the matter distribution of the Universe using the redshifted 21cm line of neutral hydrogen gas (HI). Several ongoing and upcoming radio interferometers, such as Tianlai, CHIME, HERA, HIRAX, etc., are using this technique. These instruments are designed to map large swaths of the sky by drift scanning over periods of many months. One of the challenges of the observations is that the daytime data are contaminated by strong radio signals from the Sun. In the case of Tianlai, this results in almost half of the measured data being unusable. We try to address this issue by developing an algorithm for solar contamination removal (AlgoSCR) from the radio data. The algorithm is based on an eigenvalue analysis of the visibility matrix and hence is applicable only to interferometers. We apply AlgoSCR to simulated visibilities, as well as real daytime data from the Tianlai dish array. The algorithm can reduce strong solar contamination by about 95 per cent without seriously affecting other weaker sky signals and thus makes the data usable for certain applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Beam Measurements of the Tianlai Dish Radio Telescope using an Unmanned Aerial Vehicle

Precision measurement of the beam pattern of an antenna is very important for many applications. While traditionally such measurement is often made in a microwave anechoic chamber or at a test range, measurement using an unmanned aerial vehicle offers a number of advantages: the measurement can be made for the assembled antenna on site, thus reflecting the actual characteristics of the antenna of interest, and more importantly, it can be performed for larger antennas which cannot be steered or easily measured using the anechoic chamber and test range. Here we report our beam measurement experiment with UAV for a 6 meter dish used in the Tianlai array, which is a radio astronomy experiment. Due to the dish's small collecting area, calibration with an astronomical source only allows for determining the antenna beam pattern over a very limited angular range. We describe in detail the setup of the experiment, the components of the signal transmitting system, the design of the flight path and the procedure for data processing. We find the UAV measurement of the beam pattern agrees very well with the astronomical source measurement in the main lobe, but the UAV measurement can be extended to the fourth side lobe. The measured position and width of each lobe also shows good agreement with electromagnetic field simulation. This UAV-based approach of beam pattern measurement is flexible and inexpensive, and the technique may also be applied to other experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

The Tianlai dish pathfinder array: design, operation, and performance of a prototype transit radio interferometer

ABSTRACT The Tianlai Dish Pathfinder Array is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization universe as a means for measuring large-scale cosmic structure. It performs drift scans of the sky at constant declination. We describe the design, calibration, noise level, and stability of this instrument based on the analysis of about 5% of 6200 h of on-sky observations through 2019 October. Beam pattern determinations using drones and the transit of bright sources are in good agreement, and compatible with electromagnetic simulations. Combining all the baselines, we make maps around bright sources and show that the array behaves as expected. A few hundred hours of observations at different declinations have been used to study the array geometry and pointing imperfections, as well as the instrument noise behaviour. We show that the system temperature is below 80 K for most feed antennas and that noise fluctuations decrease as expected with integration time, at least up to a few hundred seconds. Analysis of long integrations, from 10 nights of observations of the North Celestial Pole (NCP), yielded visibilities with amplitudes of 20–30 mK, consistent with the expected signal from the NCP radio sky with ${\lt}10\,$ mK precision for 1 MHz × 1 min binning. Hi-pass filtering the spectra to remove smooth spectrum signal yields a residual consistent with zero signal at the $0.5\,$ mK level.

(cosmology:) large-scale structure of the Universe↗

DOE OSTI2019

Cosmic Microwave Background Anisotropy Numerical Solution (CMBAns) I: An Introduction to $C_l$ Calculation

Cosmological Boltzmann codes are often used by researchers for calculating the CMB angular power spectra from different theoretical models, for cosmological parameter estimation, etc. Therefore, the accuracy of a Boltzmann code is of utmost importance. Different Markov Chain Monte Carlo based parameter estimation algorithms typically require 10^3 - 10^4 iterations of Boltzmann code. This makes the time complexity of such codes another critical factor. In the last two decades, several Boltzmann packages, such as CMBFAST, CAMB, CMBEasy, CLASS etc., have been developed. In this paper, we present a new cosmological Boltzmann code, CMBAns, that can be used for accurate calculation of the CMB power spectrum. At present, CMBAns is developed for a flat background matrix. It is mostly written in the C language. However, we borrowed the concept of class from C++. This gives researchers the flexibility to develop their own independent package based on CMBAns, without an in-depth understanding of the source code. We also develop multiple stand-alone facilities which can be directly compiled and run on a given parameter set. In this paper, we discuss all the mathematical formulation, approximation schemes, integration methods etc., that are used in CMBAns. The package will be made available through github for public use in the near future.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2019

sitoolbox: A package for Bayesian estimation of the isotropy violation in the CMB sky

ABSTRACT The standard model of cosmology predicts a statistically isotropic (SI) cosmic microwave background (CMB) sky. However, the SI violation signals are always present in an observed sky-map. Different cosmological artefacts, measurement effects, and unavoidable effects during data analysis, etc. may lead to isotropy violation in an otherwise SI sky. Therefore, a proper data analysis technique should calculate all these SI violation signals, so that they can be matched with SI violation signals from the known sources, and then conclude if there is any intrinsic SI violation in the CMB sky. In one of our recent works, we presented a general Bayesian formalism for measuring the isotropy violation signals in the CMB sky in the presence of an idealized isotropic noise. In this paper, we have extended the mechanism and develop a software package, sitoolbox, for measuring SI violation in presence of anisotropic noise and masking.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2011

Sparse Solutions for Single Class SVMs: A Bi-Criterion Approach

In this paper we propose an innovative learning algorithm - a variation of One-class nu Support Vector Machines (SVMs) learning algorithm to produce sparser solutions with much reduced computational complexities. The proposed technique returns an approximate solution, nearly as good as the solution set obtained by the classical approach, by minimizing the original risk function along with a regularization term. We introduce a bi-criterion optimization that helps guide the search towards the optimal set in much reduced time. The outcome of the proposed learning technique was compared with the benchmark one-class Support Vector machines algorithm which more often leads to solutions with redundant support vectors. Through out the analysis, the problem size for both optimization routines was kept consistent. We have tested the proposed algorithm on a variety of data sources under different conditions to demonstrate the effectiveness. In all cases the proposed algorithm closely preserves the accuracy of standard one-class nu SVMs while reducing both training time and test time by several factors.

Das, Santanu↗

NASA NTRS2010

Multiple Kernel Learning for Heterogeneous Anomaly Detection: Algorithm and Aviation Safety Case Study

The world-wide aviation system is one of the most complex dynamical systems ever developed and is generating data at an extremely rapid rate. Most modern commercial aircraft record several hundred flight parameters including information from the guidance, navigation, and control systems, the avionics and propulsion systems, and the pilot inputs into the aircraft. These parameters may be continuous measurements or binary or categorical measurements recorded in one second intervals for the duration of the flight. Currently, most approaches to aviation safety are reactive, meaning that they are designed to react to an aviation safety incident or accident. In this paper, we discuss a novel approach based on the theory of multiple kernel learning to detect potential safety anomalies in very large data bases of discrete and continuous data from world-wide operations of commercial fleets. We pose a general anomaly detection problem which includes both discrete and continuous data streams, where we assume that the discrete streams have a causal influence on the continuous streams. We also assume that atypical sequence of events in the discrete streams can lead to off-nominal system performance. We discuss the application domain, novel algorithms, and also discuss results on real-world data sets. Our algorithm uncovers operationally significant events in high dimensional data streams in the aviation industry which are not detectable using state of the art methods

Das, Santanu↗

NASA NTRS2010

Near Real-Time Optimal Prediction of Adverse Events in Aviation Data

The prediction of anomalies or adverse events is a challenging task, and there are a variety of methods which can be used to address the problem. In this paper, we demonstrate how to recast the anomaly prediction problem into a form whose solution is accessible as a level-crossing prediction problem. The level-crossing prediction problem has an elegant, optimal, yet untested solution under certain technical constraints, and only when the appropriate modeling assumptions are made. As such, we will thoroughly investigate the resilience of these modeling assumptions, and show how they affect final performance. Finally, the predictive capability of this method will be assessed by quantitative means, using both validation and test data containing anomalies or adverse events from real aviation data sets that have previously been identified as operationally significant by domain experts. It will be shown that the formulation proposed yields a lower false alarm rate on average than competing methods based on similarly advanced concepts, and a higher correct detection rate than a standard method based upon exceedances that is commonly used for prediction.

Martin, Rodney Alexander↗

NASA NTRS2009

nu-Anomica: A Fast Support Vector Based Novelty Detection Technique

In this paper we propose nu-Anomica, a novel anomaly detection technique that can be trained on huge data sets with much reduced running time compared to the benchmark one-class Support Vector Machines algorithm. In -Anomica, the idea is to train the machine such that it can provide a close approximation to the exact decision plane using fewer training points and without losing much of the generalization performance of the classical approach. We have tested the proposed algorithm on a variety of continuous data sets under different conditions. We show that under all test conditions the developed procedure closely preserves the accuracy of standard one-class Support Vector Machines while reducing both the training time and the test time by 5 - 20 times.

Das, Santanu↗

NASA NTRS2009

Highly Scalable Matching Pursuit Signal Decomposition Algorithm

Matching Pursuit Decomposition (MPD) is a powerful iterative algorithm for signal decomposition and feature extraction. MPD decomposes any signal into linear combinations of its dictionary elements or atoms . A best fit atom from an arbitrarily defined dictionary is determined through cross-correlation. The selected atom is subtracted from the signal and this procedure is repeated on the residual in the subsequent iterations until a stopping criterion is met. The reconstructed signal reveals the waveform structure of the original signal. However, a sufficiently large dictionary is required for an accurate reconstruction; this in return increases the computational burden of the algorithm, thus limiting its applicability and level of adoption. The purpose of this research is to improve the scalability and performance of the classical MPD algorithm. Correlation thresholds were defined to prune insignificant atoms from the dictionary. The Coarse-Fine Grids and Multiple Atom Extraction techniques were proposed to decrease the computational burden of the algorithm. The Coarse-Fine Grids method enabled the approximation and refinement of the parameters for the best fit atom. The ability to extract multiple atoms within a single iteration enhanced the effectiveness and efficiency of each iteration. These improvements were implemented to produce an improved Matching Pursuit Decomposition algorithm entitled MPD++. Disparate signal decomposition applications may require a particular emphasis of accuracy or computational efficiency. The prominence of the key signal features required for the proper signal classification dictates the level of accuracy necessary in the decomposition. The MPD++ algorithm may be easily adapted to accommodate the imposed requirements. Certain feature extraction applications may require rapid signal decomposition. The full potential of MPD++ may be utilized to produce incredible performance gains while extracting only slightly less energy than the standard algorithm. When the utmost accuracy must be achieved, the modified algorithm extracts atoms more conservatively but still exhibits computational gains over classical MPD. The MPD++ algorithm was demonstrated using an over-complete dictionary on real life data. Computational times were reduced by factors of 1.9 and 44 for the emphases of accuracy and performance, respectively. The modified algorithm extracted similar amounts of energy compared to classical MPD. The degree of the improvement in computational time depends on the complexity of the data, the initialization parameters, and the breadth of the dictionary. The results of the research confirm that the three modifications successfully improved the scalability and computational efficiency of the MPD algorithm. Correlation Thresholding decreased the time complexity by reducing the dictionary size. Multiple Atom Extraction also reduced the time complexity by decreasing the number of iterations required for a stopping criterion to be reached. The Course-Fine Grids technique enabled complicated atoms with numerous variable parameters to be effectively represented in the dictionary. Due to the nature of the three proposed modifications, they are capable of being stacked and have cumulative effects on the reduction of the time complexity.

Christensen, Daniel↗

NASA NTRS2008

Algorithms for Spectral Decomposition with Applications to Optical Plume Anomaly Detection

The analysis of spectral signals for features that represent physical phenomenon is ubiquitous in the science and engineering communities. There are two main approaches that can be taken to extract relevant features from these high-dimensional data streams. The first set of approaches relies on extracting features using a physics-based paradigm where the underlying physical mechanism that generates the spectra is used to infer the most important features in the data stream. We focus on a complementary methodology that uses a data-driven technique that is informed by the underlying physics but also has the ability to adapt to unmodeled system attributes and dynamics. We discuss the following four algorithms: Spectral Decomposition Algorithm (SDA), Non-Negative Matrix Factorization (NMF), Independent Component Analysis (ICA) and Principal Components Analysis (PCA) and compare their performance on a spectral emulator which we use to generate artificial data with known statistical properties. This spectral emulator mimics the real-world phenomena arising from the plume of the space shuttle main engine and can be used to validate the results that arise from various spectral decomposition algorithms and is very useful for situations where real-world systems have very low probabilities of fault or failure. Our results indicate that methods like SDA and NMF provide a straightforward way of incorporating prior physical knowledge while NMF with a tuning mechanism can give superior performance on some tests. We demonstrate these algorithms to detect potential system-health issues on data from a spectral emulator with tunable health parameters.

Srivastava, Askok N.↗