Search NASA⌕ Search

SEARCH · Search NASA

Results for “map validation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
DOE OSTI2015

WE-EF-210-08: BEST IN PHYSICS (IMAGING): 3D Prostate Segmentation in Ultrasound Images Using Patch-Based Anatomical Feature

Purpose: Transrectal ultrasound (TRUS) is the standard imaging modality for the image-guided prostate-cancer interventions (e.g., biopsy and brachytherapy) due to its versatility and real-time capability. Accurate segmentation of the prostate plays a key role in biopsy needle placement, treatment planning, and motion monitoring. As ultrasound images have a relatively low signal-to-noise ratio (SNR), automatic segmentation of the prostate is difficult. However, manual segmentation during biopsy or radiation therapy can be time consuming. We are developing an automated method to address this technical challenge. Methods: The proposed segmentation method consists of two major stages: the training stage and the segmentation stage. During the training stage, patch-based anatomical features are extracted from the registered training images with patient-specific information, because these training images have been mapped to the new patient’ images, and the more informative anatomical features are selected to train the kernel support vector machine (KSVM). During the segmentation stage, the selected anatomical features are extracted from newly acquired image as the input of the well-trained KSVM and the output of this trained KSVM is the segmented prostate of this patient. Results: This segmentation technique was validated with a clinical study of 10 patients. The accuracy of our approach was assessed using the manual segmentation. The mean volume Dice Overlap Coefficient was 89.7±2.3%, and the average surface distance was 1.52 ± 0.57 mm between our and manual segmentation, which indicate that the automatic segmentation method works well and could be used for 3D ultrasound-guided prostate intervention. Conclusion: We have developed a new prostate segmentation approach based on the optimal feature learning framework, demonstrated its clinical feasibility, and validated its accuracy with manual segmentation (gold standard). This segmentation technique could be a useful tool for image-guided interventions in prostate-cancer diagnosis and treatment. This research is supported in part by DOD PCRP Award W81XWH-13-1-0269, and National Cancer Institute (NCI) Grant CA114313.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

FLS2-BAK1 Extracellular Domain Interaction Sites Required for Defense Signaling Activation

Signaling initiation by receptor-like kinases (RLKs) at the plasma membrane of plant cells often requires regulatory leucine-rich repeat (LRR) RLK proteins such as SERK or BIR proteins. The present work examined how the microbe-associated molecular pattern (MAMP) receptor FLS2 builds signaling complexes with BAK1 (SERK3). We first, using in vivo methods that validate separate findings by others, demonstrated that flg22 (flagellin epitope) ligand-initiated FLS2-BAK1 extracellular domain interactions can proceed independent of intracellular domain interactions. We then explored a candidate SERK protein interaction site in the extracellular domains (ectodomains; ECDs) of the significantly different receptors FLS2, EFR (MAMP receptors), PEPR1 (damage-associated molecular pattern (DAMP) receptor), and BRI1 (hormone receptor). Repeat conservation mapping revealed a cluster of conserved solvent-exposed residues near the C-terminus of models of the folded LRR domains. However, site-directed mutagenesis of this conserved site in FLS2 did not impair FLS2-BAK1 ECD interactions, and mutations in the analogous site of EFR caused receptor maturation defects. Hence this conserved LRR C-terminal region apparently has functions other than mediating interactions with BAK1. In vivo tests of the subsequently published FLS2-flg22-BAK1 ECD co-crystal structure were then performed to functionally evaluate some of the unexpected configurations predicted by that crystal structure. In support of the crystal structure data, FLS2-BAK1 ECD interactions were no longer detected in in vivo co-immunoprecipitation experiments after site-directed mutagenesis of the FLS2 BAK1-interaction residues S554, Q530, Q627 or N674. In contrast, in vivo FLS2-mediated signaling persisted and was only minimally reduced, suggesting residual FLS2-BAK1 interaction and the limited sensitivity of co-immunoprecipitation data relative to in vivo assays for signaling outputs. However, Arabidopsis plants expressing FLS2 with the Q530A+Q627A double mutation were impaired both in detectable interaction with BAK1 and in FLS2-mediated responses, lending overall support to current models of FLS2 structure and function.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2001

Geologic Validation of EO-1 Hyperion using AVIRIS

Over the last few years, the field of imaging spectrometry has grown rapidly as new instruments and analysis techniques have been developed. The launch of Hyperion as part of NASA's EO-1 program represents a significant landmark in the progression of the technology; the first spaceborne imaging spectrometry system. AIG, in cooperation with CSIRO, is evaluating, validating, and demonstrating use of EO-1 Hyperion hyperspectral data for geologic applications. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) plays a pivotal role in this effort. This manuscript describes the geologic test sites being used and progress on preparing 'ground truth' for anticipated Hyperion data collects using AVIRIS data. Preliminary findings are extrapolated to predict the effect of instrument characteristics and performance on geologic mapping using Hyperion.

Kruse, F. A.↗
NASA NTRS2012

Global Infrared Observations of Roughness Induced Transition on the Space Shuttle Orbiter

High resolution infrared observations made from a mobile ground based optical system captured the laminar-to-turbulent boundary layer transition process as it occurred during Space Shuttle Endeavour's return to earth following its final mission in 2011. The STS-134 imagery was part of a larger effort to demonstrate an emerging and reliable non-intrusive global thermal measurement capability and to complement a series of boundary layer transition flight experiments that were flown on the Shuttle. The STS-134 observations are believed to be the first time that the development and movement of a hypersonic boundary layer transition front has been witnessed in flight over the entire vehicle surface and in particular, at unprecedented spatial resolution. Additionally, benchmark surface temperature maps of the Orbiter lower surface collected over multiple flights and spanning a Mach range of 18 to 6 are now available and represent an opportunity for collaborative comparison with computational techniques focused on hypersonic transition and turbulence modeling. The synergy of the global temperature maps with the companion in-situ thermocouple measurements serve as an example of the effective leveraging of resources to achieve a common goal of advancing our understanding of the complex nature of high Mach number transition. It is shown that quantitative imaging can open the door to a multitude of national and international opportunities for partnership associated with flight-testing and subsequent validation of numerical simulation techniques. The quantitative imaging applications highlighted in this paper offer unique and complementary flight measurement alternatives and suggest collaborative instrumentation opportunities to advance the state of the art in transition prediction and maximize the return on investment in terms of developmental flight tests for future vehicle designs.

Horvath, Thomas J.↗
NASA NTRS2012

Advances in Spectral-Spatial Classification of Hyperspectral Images

Recent advances in spectral-spatial classification of hyperspectral images are presented in this paper. Several techniques are investigated for combining both spatial and spectral information. Spatial information is extracted at the object (set of pixels) level rather than at the conventional pixel level. Mathematical morphology is first used to derive the morphological profile of the image, which includes characteristics about the size, orientation and contrast of the spatial structures present in the image. Then the morphological neighborhood is defined and used to derive additional features for classification. Classification is performed with support vector machines using the available spectral information and the extracted spatial information. Spatial post-processing is next investigated to build more homogeneous and spatially consistent thematic maps. To that end, three presegmentation techniques are applied to define regions that are used to regularize the preliminary pixel-wise thematic map. Finally, a multiple classifier system is defined to produce relevant markers that are exploited to segment the hyperspectral image with the minimum spanning forest algorithm. Experimental results conducted on three real hyperspectral images with different spatial and spectral resolutions and corresponding to various contexts are presented. They highlight the importance of spectral-spatial strategies for the accurate classification of hyperspectral images and validate the proposed methods.

Fauvel, Mathieu↗
NASA NTRS2015

The Rapid Response Radiation Survey (R3S) Mission Using the HiSat Conformal Satellite Architecture

The Rapid Response Radiation Survey (R3S) experiment, designed as a quick turnaround mission to make radiation measurements in Low Earth Orbit (LEO), will fly as a hosted payload in partnership with NovaWurks using their Hyper-integrated Satlet (HISat) architecture. The need for the mission arises as the Nowcast of Atmospheric Ionization Radiation for Aviation Safety (NAIRAS) model moves from a research effort into an operational radiation assessment tool. Currently, airline professionals are the second largest demographic of radiation workers and to date their radiation exposure is undocumented in the USA. The NAIRAS model seeks to fill this information gap. The data collected by R3S, in addition to the complementary data from a NASA Langley Research Center (LaRC) atmospheric balloon mission entitled Radiation Dosimetry Experiment (RaD-X), will validate exposure prediction capabilities of NAIRAS. The R3S mission collects total dose and radiation spectrum measurements using a Teledyne μDosimeter and a Liulin-6SA2 LED spectrometer. These two radiation sensors provide a cross correlated radiometric measurement in combination with the Honeywell HMR2300 Smart Digital Magnetometer. The magnetometer assesses the Earth's magnetic field in the LEO environment and allows radiation dose to be mapped as a function of the Earth's magnetic shielding. R3S is also unique in that the radiation sensors will be exposed on the outer surface of the spacecraft, possibly making this the first measurements of the LEO radiation environment with bare sensors. Viability of R3S as an extremely fast turnaround mission is due, in part, to the nature of the robust, well-defined interfaces of the conformal satellite HiSat Architecture. The HiSat architecture, which was developed with the support of the Defense Advanced Research Projects Agency's (DARPA's) Phoenix Program, enabled the R3S system to advance from the first concept to delivery of preliminary design review (PDR) level documents in 29 calendar days. The architecture allows for interface complexities between the specific devices and the satellite bus to be resolved in a standardized interface control document (ICD). The ICD provided a readymade framework to interface to the modular satellite bus. This modularity allowed for approximately 90% of the R3S system to be designed and fabricated in two months without constraint of the hosting satellite's development cycle. This paper discusses the development of the R3S experiment as made possible by use of the HiSat architecture. The system design and operational modes of the experiment are described, as well as the experiment interfaces to the HiSat satellite via the user defined adapter (UDA) provided by NovaWurks. This paper outlines the steps taken by the project to execute the R3S mission in the 4 months of design, build, and test. Additionally portrayed is the ground work done at LaRC to posture the organization for a fast response and the process by which the opportunity was identified as aligning with key strategic goals. Finally, a description of the engineering process is provided, including the use of facilitated rapid/concurrent engineering sessions, the associated documentation, and the review process employed.

Miller, Nathanael A.↗
NASA NTRS2018

Multibaseline TanDEM-X Mangrove Height Estimation: The Selection of the Vertical Wavenumber

We generated a large-scale mangrove forest height map using multiple TanDEM-X (TDX) interferometric synthetic aperture radar (InSAR) acquisitions with various spatial baselines in order to improve the height estimation accuracy across a wide range of forest heights. The forest height inversion using InSAR data is strongly dependent upon the vertical wavenumber (i.e., perpendicular baseline). First, we investigated the role of the vertical wavenumber in forest height inversion from InSAR data using the sensitivity of the interferometric (volume) coherence to forest height. We used corrected but lower resolution and accuracy Shuttle Radar Topography Mission (SRTM) mangrove height maps as a priori information over Akanda and Pongara National Parks in Gabon to estimate lower and upper boundaries of the vertical wavenumber over test sites from the measured coherence-to-height sensitivity. Only TDX acquisitions within the boundaries of the vertical wavenumber were selected and combined for multibaseline mangrove height inversion. Mangrove forest height was obtained with multibaseline TDX acquisitions and was validated against the reference height derived from field measurement data providing improvements in multibaseline inversion over existing height estimates (i.e., SRTM height) and single-baseline inversions (multibaseline inversion: r(exp 2) = 0.98, root mean square error (RMSE) of 2.73 m; SRTM height: r(exp 2) = 0.86, RMSE = 7.21 m; single-baseline inversions: r(exp 2) = 0.08-0.97, RMSE = 3.86-11.10 m). As a result, to accurately estimate forest heights over a wide range (3-60 m), multibaseline InSAR acquisitions (at least three different baselines) are needed to exclude biases associated with the vertical wavenumber in forest height inversion.

vertical wavenumber.↗
NASA NTRSDate not supplied

A Microphone Phased Array for Launch Acoustics Application

A new, portable, phased array of microphones is built at NASA Ames Research Center specifically for the harsh environment encountered in launch pads of rocket vehicles. It uses 70 piezoresistive, dynamic pressure sensors, optimally distributed on a 10.5ft diameter open frame dome structure. The open frame is light yet robust to sustain the high wind load of typical seaside launch pads, and the blast and acoustic loads from the launch. A 200-ft long cable bundle carries the microphone signals to a weather-protected cabinet containing the data systems, and allows for the placement of the array tall structures. The array is equipped with an infra-red camera and a visible wavelength camera for imaging the launch pad. The beamformed noise maps will be superimposed on the video footages collected by the cameras for correct identification of the noise sources. The array is tested with very loud speaker sources to determine the beamforming ability using different schemes, such as the conventional beamforming, orthogonal-functional beamforming, and spectral element methods. A comparative study was conducted to determine the minimum attainable spatial resolution. Data from a validation test in an outdoor setup, where the array was placed on a Telehandler at several different heights and distances from either a single or a pair of speakers, is presented in this paper. The final goal is to use the array during NASA’s Artemis-II launch to determine the distribution of noise sources on the Mobile Launcher.

Acoustics↗
DOE OSTI2015

Validating Laser-Induced Birefringence Theory with Plasma Interferometry

Intense laser beams crossing paths in plasma is theorized to induce birefringence in the medium, resulting from density and refractive index modulations that affect the polarization of incoming light. The goal of the associated experiment, conducted on Janus at Lawrence Livermore’s Jupiter Laser Facility, was to create a tunable laser-plasma waveplate to verify the relationship between dephasing angle and beam intensity, plasma density, plasma temperature, and interaction length. Interferometry analysis of the plasma channel was performed to obtain a density map and to constrain temperature measured from Thomson scattering. Various analysis techniques, including Fast Fourier transform (FFT) and two variations of fringe-counting, were tried because interferograms captured in this experiment contained unusual features such as fringe discontinuity at channel edges, saddle points, and islands. The chosen method is flexible, semi-automated, and uses a fringe tracking algorithm on a reduced image of pre-traced synthetic fringes. Ultimately, a maximum dephasing angle of 49.6° was achieved using a 1200 μm interaction length, and the experimental results appear to agree with predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

SU-F-I-79: Validation of a Monte Carlo Radiographic Simulation

Purpose: To validate a radiographic simulation in order to estimate patient dose due to clinically-used radiography protocols. Methods: A Monte Carlo simulation was created to simulate a radiographic x-ray beam using GEANT4. Initial validation was performed according to a portion of TG 195. Computational NURBS-based phantoms were used simulate patients of varying ages and sizes. The deposited energy in the phantom is output by the simulation. The exposure in air from a clinically used radiography unit was measured at 100 cm for various tube potentials. 10 million photons were simulated with 1 cubic centimeter of air located 100 cm from the source, and the total absorbed dose was noted. The normalization factor was determined by taking a ratio of the measured dose in air to the simulated dose in air. Dose to individual voxels is calculated using the energy deposition map along with the voxelized and segmented phantom and the normalization factor. Finally, the effective dose is calculated using the ICRP methodology and tissue weighting factors. Results: This radiography simulation allows for the calculation and visualization of the energy deposition map within a voxelized phantom. The ratio of exposure, measured using an ionization chamber, to air in the simulation was determined. Since the simulation output is calibrated to match the exposure of a given clinical radiographic x-ray tube, the dose map may be visualized. This will also allow for absorbed dose estimation in specific organs or tissues as well as a whole body effective dose estimation. Conclusion: This work indicates that our Monte Carlo simulation may be used to estimate the radiation dose from clinical radiographic protocols. This will allow for an estimate of radiographic dose from various examinations without the use of traditional methods such as thermoluminescent dosimeters and body phantoms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Effects of digital fringe projection operational parameters on detecting powder bed defects in additive manufacturing

Additive manufacturing is a technology transforming traditional production timelines. Specifically, metal additive manufacturing (MAM) has been increasingly adopted by a variety of industries, not only to prototype, but also to fulfill full production scale applications with much lower lead times. Like any maturing manufacturing technology, developments in verifying and validating processes are necessary to support continuous growth. Due to the complex nature of MAM, part quality and repeatability remain integral challenges that inhibit further adoption of MAM for critical component production. In this study, we present data taken from a developing in-process monitoring system designed to measure and detect powder bed defects (PBDs) in powder bed fusion MAM systems using surface height maps created with structured light illumination. We showcase the feasibility of the monitoring technique for in-process implementation by detecting streak PBDs with varying severities (height, width) created in a lab environment. We present results of powder bed measurements for varying experimental parameters of the structured light system such as illumination angle, illumination pattern, and number of illuminations. We also present an expression used to determine experimental height noise based on input parameters for PBD detection based on the instrument transfer function of the structured light monitoring system for arbitrary pixel intensity noise contributions. In conclusion, with the results of PBD detection across multiple experimental measurement parameters, we provide a best practices approach to in-process implementation of the monitoring system in powder bed fusion manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2013

A mixed quantum-classical Liouville study of the population dynamics in a model photo-induced condensed phase electron transfer reaction

We apply two approximate solutions of the quantum-classical Liouville equation (QCLE) in the mapping representation to the simulation of the laser-induced response of a quantum subsystem coupled to a classical environment. These solutions, known as the Poisson Bracket Mapping Equation (PBME) and the Forward-Backward (FB) trajectory solutions, involve simple algorithms in which the dynamics of both the quantum and classical degrees of freedom are described in terms of continuous variables, as opposed to standard surface-hopping solutions in which the classical degrees of freedom hop between potential energy surfaces dictated by the discrete adiabatic state of the quantum subsystem. The validity of these QCLE-based solutions is tested on a non-trivial electron transfer model involving more than two quantum states, a time-dependent Hamiltonian, strong subsystem-bath coupling, and an initial energy shift between the donor and acceptor states that depends on the strength of the subsystem-bath coupling. In particular, we calculate the time-dependent population of the photoexcited donor state in response to an ultrafast, on-resonance pump pulse in a three-state model of an electron transfer complex that is coupled asymmetrically to a bath of harmonic oscillators through the optically dark acceptor state. Within this approach, the three-state electron transfer complex is treated quantum mechanically, while the bath oscillators are treated classically. When compared to the more accurate QCLE-based surface-hopping solution and to the numerically exact quantum results, we find that the PBME solution is not capable of qualitatively capturing the population dynamics, whereas the FB solution is. However, when the subsystem-bath coupling is decreased (which also decreases the initial energy shift between the donor and acceptor states) or the initial shift is removed altogether, both the PBME and FB results agree better with the QCLE-based surface-hopping results. These findings highlight the challenges posed by various conditions such as a time-dependent external field, the strength of the subsystem-bath coupling, and the degree of asymmetry on the accuracy of the PBME and FB algorithms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Using a pulsed array to determine the position of a defect

The frequency dependent directivity of periodic permanent magnet electromagnetic acoustic transducers has been demonstrated and utilized to create a new NDT technique. The transducer, when pulsed with an appropriate broadband signal creates an interesting wave-front. This generated wave-front is smooth and continuous over a large solid angle, and has a frequency that varies as a function of angle to the array. Analytic models created to simulate this phenomena have been validated with experimental results using ultrasonic waves in metal samples. This pulsed approach provides a rapid means of flooding a region of space with a wave-front, whereby any wave that scatters or reflects off a body to a detector will have a distinct arrival time and frequency. This was experimentally verified using a periodic permanent magnetic EMAT array to create a two dimensional map of defects in an aluminum bar. Steering of the ultrasonic beam, by driving the array with a narrow-band signal, was also demonstrated and used as a comparison for the pulsed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker–Planck equations). Here, we analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker–Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

LOCAL INTERSTELLAR MEDIUM: SIX YEARS OF DIRECT SAMPLING BY IBEX

The Interstellar Boundary Explorer (IBEX) has been directly observing neutral atoms from the local interstellar medium for the last six years (2009–2014). This paper ties together the 14 studies in this Astrophysical Journal Supplement Series Special Issue, which collectively describe the IBEX interstellar neutral results from this epoch and provide a number of other relevant theoretical and observational results. Interstellar neutrals interact with each other and with the ionized portion of the interstellar population in the “pristine” interstellar medium ahead of the heliosphere. Then, in the heliosphere's close vicinity, the interstellar medium begins to interact with escaping heliospheric neutrals. In this study, we compare the results from two major analysis approaches led by IBEX groups in New Hampshire and Warsaw. We also directly address the question of the distance upstream to the pristine interstellar medium and adjust both sets of results to a common distance of ∼1000 AU. The two analysis approaches are quite different, but yield fully consistent measurements of the interstellar He flow properties, further validating our findings. While detailed error bars are given for both approaches, we recommend that for most purposes, the community use “working values” of ∼25.4 km s{sup −1}, ∼75.°7 ecliptic inflow longitude, ∼ −5.°1 ecliptic inflow latitude, and ∼7500 K temperature at ∼1000 AU upstream. Finally, we briefly address future opportunities for even better interstellar neutral observations to be provided by the Interstellar Mapping and Acceleration Probe mission, which was recommended as the next major Heliophysics mission by the NRC's 2013 Decadal Survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗