Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “astronomical methods”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Defending Planet Earth: Near-Earth Object Surveys and Hazard Mitigation Strategies

The United States spends approximately four million dollars each year searching for near-Earth objects (NEOs). The objective is to detect those that may collide with Earth. The majority of this funding supports the operation of several observatories that scan the sky searching for NEOs. This, however, is insufficient in detecting the majority of NEOs that may present a tangible threat to humanity. A significantly smaller amount of funding supports ways to protect the Earth from such a potential collision or "mitigation." In 2005, a Congressional mandate called for NASA to detect 90 percent of NEOs with diameters of 140 meters of greater by 2020. Defending Planet Earth: Near-Earth Object Surveys and Hazard Mitigation Strategies identifies the need for detection of objects as small as 30 to 50 meters as these can be highly destructive. The book explores four main types of mitigation including civil defense, "slow push" or "pull" methods, kinetic impactors and nuclear explosions. It also asserts that responding effectively to hazards posed by NEOs requires national and international cooperation. Defending Planet Earth: Near-Earth Object Surveys and Hazard Mitigation Strategies is a useful guide for scientists, astronomers, policy makers and engineers.

Source record↗
NASA NTRS2018

Development of a Robust, Efficient Process to Produce Scalable, Superconducting Kilopixel Far-IR Detector Arrays

The far-IR band is uniquely suited to study the physical conditions in the interstellar medium from nearby sources out to the highest redshifts. FIR imaging and spectroscopy instrumentation using incoherent superconducting bolometers represents a high sensitivity technology for many future suborbital and space missions, including the Origins Space Telescope. Robust, high sensitivity detector arrays with several 104 pixels, large focal plane filling factors, and low cosmic ray cross sections that operate over the entire far-IR regime are required for such missions. These arrays could consist of smaller sub-arrays, in case they are tileable. The TES based Backshort Under Grid array architecture which our group has fielded in a number of FIR cameras, is a good candidate to meet these requirements: BUGs are tileable, and with the integration of the SQUID multiplexer scaleable beyond wafer sizes; they provide high filling factors, low cosmic cross section and have been demonstrated successfully in far-infrared astronomical instrumentation. However, the production of BUGs with integrated readout multiplexers has many time and resource consuming process steps. In order to meet the requirement of robustness and efficiency on the production of future arrays, we have developed a new method to provide the superconducting connection of BUG detectors to the readout multiplexers or general readout boards behind the detectors. This approach should allow us to reach the goal to produce reliable, very large detector arrays for future FIR missions.

Staguhn, Johannes↗
NASA NTRS2023

Unusually Long Path Length for a Nearly Scatter-Free Solar Particle Event Observed By Solar Orbiter at 0.43 au

Context. After their acceleration and release at the Sun, solar energetic particles (SEPs) are injected into the interplanetary medium and are bound to the interplanetary magnetic field (IMF) by the Lorentz force. The expansion of the IMF close to the Sun focuses the particle pitch-angle distribution, and scattering counteracts this focusing. Solar Orbiter observed an unusual solar particle event on 9 April 2022 when it was at 0.43 astronomical units (au) from the Sun. Aims. We show that the inferred IMF along which the SEPs traveled was about three times longer than the nominal length of the Parker spiral and provide an explanation for this apparently long path. Methods. We used velocity dispersion analysis (VDA) information to infer the spiral length along which the electrons and ions traveled and infer their solar release times and arrival direction. Results. The path length inferred from VDA is approximately three times longer than the nominal Parker spiral. Nevertheless, the pitch-angle distribution of the particles of this event is highly anisotropic, and the electrons and ions appear to be streaming along the same IMF structures. The angular width of the streaming population is estimated to be approximately 30 degrees. The highly anisotropic ion beam was observed for more than 12 h. This may be due to the low level of fluctuations in the IMF, which in turn is very probably due to this event being inside an interplanetary coronal mass ejection The slow and small rotation in the IMF suggests a flux-rope structure. Small flux dropouts are associated with very small changes in pitch angle, which may be explained by different flux tubes connecting to different locations in the flare region. Conclusions. The unusually long path length along which the electrons and ions have propagated virtually scatter-free together with the short-term flux dropouts offer excellent opportunities to study the transport of SEPs within interplanetary structures. The 9 April 2022 solar particle event offers an especially rich number of unique observations that can be used to limit SEP transport models.

Robert F. Wimmer-Schweingruber↗
DOE OSTI2024

Overview of the distributed image processing infrastructure to produce the Legacy Survey of Space and Time

The Vera C. Rubin Observatory is preparing to execute the most ambitious astronomical survey ever attempted, the Legacy Survey of Space and Time (LSST). Currently the final phase of construction is under way in the Chilean Andes, with the Observatory’s ten-year science mission scheduled to begin in 2025. Rubin’s 8.4-meter telescope will nightly scan the southern hemisphere collecting imagery in the wavelength range 320–1050 nm covering the entire observable sky every 4 nights using a 3.2 gigapixel camera, the largest imaging device ever built for astronomy. Automated detection and classification of celestial objects will be performed by sophisticated algorithms on high-resolution images to progressively produce an astronomical catalog eventually composed of 20 billion galaxies and 17 billion stars and their associated physical properties. In this article we present an overview of the system currently being constructed to perform data distribution as well as the annual campaigns which reprocess the entire image dataset collected since the beginning of the survey. These processing campaigns will utilize computing and storage resources provided by three Rubin data facilities (one in the US and two in Europe). Each year a Data Release will be produced and disseminated to science collaborations for use in studies comprising four main science pillars: probing dark matter and dark energy, taking inventory of solar system objects, exploring the transient optical sky and mapping the Milky Way. Also presented is the method by which we leverage some of the common tools and best practices used for management of large-scale distributed data processing projects in the high energy physics and astronomy communities. We also demonstrate how these tools and practices are utilized within the Rubin project in order to overcome the specific challenges faced by the Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

sOPTICS: a modified density-based algorithm for identifying galaxy groups/clusters and brightest cluster galaxies

A direct approach to studying the galaxy–halo connection is to analyse groups and clusters of galaxies that trace the underlying dark matter haloes, emphasizing the importance of identifying galaxy clusters and their associated brightest cluster galaxies (BCGs). In this work, we test and propose a robust density-based clustering algorithm that outperforms the traditional Friends-of-Friends (FoF) algorithm in the currently available galaxy group/cluster catalogues. Our new approach is a modified version of the Ordering Points To Identify the Clustering Structure (OPTICS) algorithm, which accounts for line-of-sight positional uncertainties due to redshift space distortions by incorporating a scaling factor, and is thereby referred to as sOPTICS. When tested on both a galaxy group catalogue based on semi-analytic galaxy formation simulations and observational data, our algorithm demonstrated robustness to outliers and relative insensitivity to hyperparameter choices. In total, we compared the results of eight clustering algorithms. The proposed density-based clustering method, sOPTICS, outperforms FoF in accurately identifying giant galaxy clusters and their associated BCGs in various environments with higher purity and recovery rate, also successfully recovering 115 BCGs out of 118 reliable BCGs from a large galaxy sample. Furthermore, when applied to an independent observational catalogue without extensive re-tuning, sOPTICS maintains high recovery efficiency, confirming its flexibility and effectiveness for large-scale astronomical surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Astrometry in two-photon interferometry using an Earth rotation fringe scan

Optical interferometers may not require a phase-stable optical link between the stations if instead sources of quantum-mechanically entangled pairs could be provided to them, enabling long baselines. We developed a new variation of this idea, proposing that photons from two different astronomical sources could be interfered at two decoupled stations. Interference products can then be calculated in post-processing or requiring only a slow, classical connection between stations. In this work, we investigated practical feasibility of this approach. We developed a Bayesian analysis method for the earth rotation fringe scanning technique and showed that in the limit of high signal-to-noise ratio it reproduced the results from a simple Fisher matrix analysis. Here we identify candidate stair pairs in the northern hemisphere, where this technique could be applied. With two telescopes with an effective collecting area of ~ 2 m 2 , we could detect fringing and measure the astrometric separation of the sources at ~100 μas precision in a few hours of observations, in agreement with previous estimates.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Red Noise–based False Alarm Thresholds for Astrophysical Periodograms via Whittle’s Approximation to the Likelihood

Astronomers who search for periodic signals using Lomb–Scargle periodograms rely on false alarm level (FAL) estimates to identify statistically significant peaks. Although FALs are often calculated from white noise models, many astronomical time series suffer from red noise. Prewhitening is a statistical technique in which a continuum model is subtracted from the log power spectrum estimate, after which the observer can proceed with a white-noise treatment. Here we present a prewhitening-based method of calculating frequency-dependent FALs. We fit power laws and autoregressive models of order 1 to each Lomb–Scargle periodogram by minimizing the Whittle approximation to the negative log-likelihood (NLL), then calculate FALs based on the best-fit model power spectrum. Our technique is a novel extension of the Whittle NLL to datasets with uneven time sampling. We demonstrate FAL calculations using observations of α Cen B, GJ 581, HD 192310, synthetic data from the radial velocity (RV) fitting challenge, and Kepler observations of a differential rotator. The Kepler data analysis shows that only true rotation signals are detected by red noise FALs, while white noise FALs suggest all spurious peaks in the low-frequency range are significant. A high-frequency sinusoid injected into α Cen B logR$'$ HK observations exceeds the 1% red noise FAL despite having only 8.9% of the power of the dominant rotation signal. In a periodogram of HD 192310 RVs, peaks associated with differential rotation and planets are detected against the 5% red noise FAL without iterative model fitting or subtraction. The software for calculating red noise–based FALs is available on GitHub.

Astrostatistics (1882)↗
NASA NTRS1972

Solar wind physics

A double-chambered gas proportional counter was constructed to detect and identify solar wind ions after acceleration by a high voltage power supply. It was determined that the best method of detecting deuterium in the solar wind is to use a tritium target as proposed for IMP H and J. The feasibility of detecting H(+) and He(+) ions of interstellar origin is considered. A program is described to carry out ground-based astronomical observations of faint, diffuse optical emission lines from interstellar gas. Hydrogen and oxygen emission lines from galactic sources were detected and the galactic and geocoronal H alpha and beta lines were clearly resolved.

Source record↗
NASA NTRS1973

Approach to rapid mission design and planning

Methods and techniques are described for implementation in automated computer systems to assess parametric data, capabilities, requirements and constraints for planning earth orbit missions. Mission planning and design procedures are defined using two types of typical missions as examples. These missions were the high energy Astronomical Observatory Satellite missions, and Small Applications Technology Satellite missions.

Green, W. G.↗
NASA NTRS1981

Quasi-periodic climatic changes on Mars and earth

Evidence of climatic changes on Mars and the earth due to geologic and astronomical variations is discussed. Finely striped ice-free bands in the Martian polar caps have been taken to indicate that long term variations in the orbit and axial tilt of Mars have precipitated these features at the rate of a mm/yr. Photogrammetric and photometric methods have contributed to measurements of the composition and depth of the Martian caps (14-46 m), and observations of higher solar energy absorption in the northern ice cap implies greater dust deposition in that region than on the south cap; however, the transport mechanisms are not well understood. Comparisons of earth and Martian climatic variations data are made, noting a lack of information on the age intervals of marine and nonmarine sediments on the earth. The possibilities of using quantitative data other than layer thickness to constrain climate models are discussed, and the slope or albedo of layers, or the spacing of polar undulations are suggested.

Cutts, J. A.↗
DOE OSTI2024

Implementing Superresolution of Nonstationary Tides with Wavelets: An Introduction to CWT_Multi

Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.

Lobo, Matthew↗
NASA NTRS1982

Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data

Detection of a periodic signal hidden in noise is frequently a goal in astronomical data analysis. This paper does not introduce a new detection technique, but instead studies the reliability and efficiency of detection with the most commonly used technique, the periodogram, in the case where the observation times are unevenly spaced. This choice was made because, of the methods in current use, it appears to have the simplest statistical behavior. A modification of the classical definition of the periodogram is necessary in order to retain the simple statistical behavior of the evenly spaced case. With this modification, periodogram analysis and least-squares fitting of sine waves to the data are exactly equivalent. Certain difficulties with the use of the periodogram are less important than commonly believed in the case of detection of strictly periodic signals. In addition, the standard method for mitigating these difficulties (tapering) can be used just as well if the sampling is uneven. An analysis of the statistical significance of signal detections is presented, with examples

Scargle, J. D.↗
NASA NTRS1985

The search for extraterrestrial life: Recent developments; Proceedings of the Symposium, Boston University, MA, June 18-21, 1984

The conference presents papers on the history of the search for extraterrestrial life, the scientific rationale and methods used in the search for other planetary systems, the detection of distant planets with the Space Telescope, planetary searches using optical astrometric interferometers, and infrared spectral identification of complex organic molecules in interstellar grains. Also considered are universal protein ancestors from hydrogen cyanide and water, astronomical sources of polarized light and their role in determining molecular chirality on earth, the universal diagrams and life in the universe, the precambrian evolution of terrestrial life and a thermodynamic approach to the occurrance and appearance of galactic life forms. Papers are also presented on the Ohio Seti program, lunar reflections of terrestrial radio leakage, the multichannel spectrum analyzer, software implementation of detection algorithms for the MCSA, the Serendip II design, galactic colonization and competition in a young galactic disk, implications of ancient and future migrations, extraterrestrial intelligence, the ineviability and the possible structures of supercivilizations, planetary, interplanetary and interstellar organic matter, and universal aspects of biological evolution.

Papagiannis, M. D.↗
NASA NTRS1994

Multiple Core Galaxies: Implications for M31

It is generally perceived that two cores cannot survive very long within the nuclear regions of a galaxy. The recent HST discovery of a double nucleus in M31 brings this question into prominence. Physical conditions in the nuclear regions of a typical galaxy help a second core survive so it can orbit for a long time, possibly for thousands of orbits. Given the nearly uniform mass density in a core, tidal forces within a core radius are compressive in all directions and help the core survive the buffeting it takes as it orbits near the center of the galaxy. We use numerical experiments to illustrate these physical principles. Modifications to the experimental method allow the full power of the experiments to be concentrated on the nuclear regions. Spatial resolution of about 0.2 parsec comfortably resolves detail within the 1.4 parsec core radius of the second, but brighter, core (P1) in M31. The same physical principles apply in other astronomical situations, such as dumbbell galaxies, galaxies orbiting near the center of a galaxy cluster, and subclustering in galaxy clusters. The experiments also illustrate that galaxy encounters and merging are quite sensitive to external tidal forces, such as those produced by the gravitational potential in a group or cluster of galaxies.

Smith, B. F.↗
NASA NTRS2021

What it Takes to Compute Highly Accurate Rovibrational Line Lists for Use in Astrochemistry

We review the Best Theory + Reliable High-resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is two-fold: a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH3, or have a high density of states such as SO2; and b) characterization of the atmospheres of exoplanets which will be one of the main areas of research in astronomy in the coming decades. The first motivation arises due to the need to eliminate lines due to common molecules in an astronomical observation in order to identify lines from new molecules, while the second motivation arises due to the need to obtain accurate molecular opacities in order to characterize the atmosphere of an exoplanet. The BTRHE strategy first consists of using high quality ab initio quantum chemical methods to obtain a global potential energy surface (PES) and dipole moment surface (DMS) that contains the proper physics. The global PES is then refined using a subset of the reliable high-resolution experimental data. The refined PES then gives energy level predictions to an accuracy similar to the reproduction accuracy of the experimental data used in the refinement step in the interpolation region (i.e., within the range of the experimental data used in the refinement step). The accuracy of the energy levels will slowly degrade as they are extrapolated to spectral regions beyond the high-resolution experimental data used in the refinement step. However, because the degradation is slow, the predicted energy levels can be used to assign new high-resolution experiments, and the data from these can then be used in a subsequent refinement step. In this way, the global PES eventually can yield highly accurate energy levels for all desired spectral regions including to very high energies and high J values. We show that IR intensities computed with the BTRHE rovibrational wavefunctions and the DMS can be very accurate provided one has minimized the fitting error of the DMS and tested the completeness of the DMS. Some examples of our work on NH3, CO2, and SO2 are given to highlight the usefulness of the BTRHE strategy and to provide ideas on how to further improve its predictive power in the future. In particular, it is shown how successive refinement steps, once new high-resolution data is available, can lead to PESs that yield highly accurate transition energies to larger spectral regions. The importance of including non-adiabatic corrections to reduce the J-dependence of errors for H-containing molecules is shown with work on NH3. Another very important aspect of the BTRHE approach is the consistency across isotopologues, which allows for highly accurate line lists for any isotopologue once one is obtained for the main isotopologue (which has more high-resolution data available for refinement).

Xinchuan Huang↗