SEARCH · Search NASA

Results for “Source code instrumentation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2016

DEEP WIDEBAND SINGLE POINTINGS AND MOSAICS IN RADIO INTERFEROMETRY: HOW ACCURATELY DO WE RECONSTRUCT INTENSITIES AND SPECTRAL INDICES OF FAINT SOURCES?

Many deep wideband wide-field radio interferometric surveys are being designed to accurately measure intensities, spectral indices, and polarization properties of faint source populations. In this paper, we compare various wideband imaging methods to evaluate the accuracy to which intensities and spectral indices of sources close to the confusion limit can be reconstructed. We simulated a wideband single-pointing (C-array, L-Band (1–2 GHz)) and 46-pointing mosaic (D-array, C-Band (4–8 GHz)) JVLA observation using a realistic brightness distribution ranging from 1 μ Jy to 100 mJy and time-, frequency-, polarization-, and direction-dependent instrumental effects. The main results from these comparisons are (a) errors in the reconstructed intensities and spectral indices are larger for weaker sources even in the absence of simulated noise, (b) errors are systematically lower for joint reconstruction methods (such as Multi-Term Multi-Frequency-Synthesis (MT-MFS)) along with A-Projection for accurate primary beam correction, and (c) use of MT-MFS for image reconstruction eliminates Clean-bias (which is present otherwise). Auxiliary tests include solutions for deficiencies of data partitioning methods (e.g., the use of masks to remove clean bias and hybrid methods to remove sidelobes from sources left un-deconvolved), the effect of sources not at pixel centers, and the consequences of various other numerical approximations within software implementations. This paper also demonstrates the level of detail at which such simulations must be done in order to reflect reality, enable one to systematically identify specific reasons for every trend that is observed, and to estimate scientifically defensible imaging performance metrics and the associated computational complexity of the algorithms/analysis procedures.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Joint tests at INL and CEA of a transient hot wire needle probe for in-pile thermal conductivity measurement

Thermal conductivity is a key property that must be known for proper design, testing, and deployment of new fuels and structural materials in nuclear reactors. Thermal conductivity is highly dependent on the physical structure, chemical composition, and the state of the material. Typically, thermal conductivity changes that occur during irradiation are currently measured out-of-pile using a 'cook and look' approach. But repeatedly removing samples from a test reactor to make measurements is expensive, has the potential to disturb phenomena of interest, and only provides understanding of the sample's end state when each measurement is made. There are also limited thermo-physical property data available for advanced fuels; and such data are needed for simulation codes, the development of next generation reactors, and advanced fuels for existing nuclear plants. Being able to quickly characterize fuel thermal conductivity during irradiation can improve the fidelity of data, reduce costs of post-irradiation examinations, increase understanding of how fuels behave under irradiation, and confirm or improve existing thermal conductivity measurement techniques. This paper discusses efforts to develop and evaluate an innovative in-pile thermal conductivity sensor based on the transient hot wire thermal conductivity method (THWM), using a single needle probe (NP) containing a line heat source and thermocouple embedded in the fuel. The sensor that has been designed and manufactured by the Idaho National Laboratory (INL) includes a unique combination of materials, geometry, and fabrication techniques that make the hot wire method suitable for in-pile applications. In particular, efforts were made to minimize the influence of the sensor and maximize fuel hot-wire heating. The probe has a thermocouple-like construction with high temperature resistant materials that remain ductile while resisting transmutation and materials interactions. THWM-NP prototypes were fabricated for both room temperature proof-of-concept evaluations and high temperature testing. Evaluations have been performed jointly by the INL and the French Alternative Energies and Atomic Energy Commission (CEA), both in Idaho Falls (USA) and in Cadarache (France), in the framework of a collaborative program for instrumentation of Material Testing Reactors. Initial tests were conducted on samples with a large range of thermal conductivities and temperatures ranging from 20 deg. C to 600 deg. C. Particularly, tests were recently performed on a sample having thermal conductivity and dimensions similar to UO{sub 2} and MOX nuclear fuels, in order to validate the ability of this sensor to operate for in-pile characterization of Light Water Reactors fuels. The results of the tests already completed at INL and CEA indicate that the Transient Hot Wire Needle Probe offers an enhanced method for in-pile detection of thermal conductivity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

High current DC negative ion source for cyclotron

A filament driven multi-cusp negative ion source has been developed for proton cyclotrons in medical applications. In Cs-free operation, continuous H{sup −} beam of 10 mA and D{sup −} beam of 3.3 mA were obtained stably at an arc-discharge power of 3 kW and 2.4 kW, respectively. In Cs-seeded operation, H{sup −} beam current reached 22 mA at a lower arc power of 2.6 kW with less co-extracted electron current. The optimum gas flow rate, which gives the highest H{sup −} current, was 15 sccm in the Cs-free operation, while it decreased to 4 sccm in the Cs-seeded operation. The relationship between H{sup −} production and the design/operating parameters has been also investigated by a numerical study with KEIO-MARC code, which gives a reasonable explanation to the experimental results of the H{sup −} current dependence on the arc power.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2001

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth Observing System) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the "Independent Pixel Approximation" or IPA. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

DOE OSTI2017

Multilevel acceleration of scattering-source iterations with application to electron transport

Acceleration/preconditioning strategies available in the SCEPTRE radiation transport code are described. A flexible transport synthetic acceleration (TSA) algorithm that uses a low-order discrete-ordinates (S N ) or spherical-harmonics (P N ) solve to accelerate convergence of a high-order S N source-iteration (SI) solve is described. Convergence of the low-order solves can be further accelerated by applying off-the-shelf incomplete-factorization or algebraic-multigrid methods. Also available is an algorithm that uses a generalized minimum residual (GMRES) iterative method rather than SI for convergence, using a parallel sweep-based solver to build up a Krylov subspace. TSA has been applied as a preconditioner to accelerate the convergence of the GMRES iterations. The methods are applied to several problems involving electron transport and problems with artificial cross sections with large scattering ratios. These methods were compared and evaluated by considering material discontinuities and scattering anisotropy. Observed accelerations obtained are highly problem dependent, but speedup factors around 10 have been observed in typical applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2019

Divergence study and emittance measurements for the electron beam emitted from a diamond pyramid

This paper reports the results of a divergence study and emittance measurements for an electron beam emitted from a single diamond pyramid in a sparse diamond field-emitter array (DFEA) cathode. DFEA cathodes are arrays of micrometer scale diamond pyramids with nanometer scale tips. A single diamond pyramid produces currents as high as 10 μA from a small surface area. DFEA cathodes are expected to be a good beam source for compact dielectric laser accelerators. For the electron beam divergence experiment, we have designed and assembled a test stand consisting of a DFEA cathode, a mesh anode, and a screen. We measured current and the size of the beam on the screen as functions of the distance between the cathode and the anode, the size of a pyramid’s base, and applied voltage on the cathode. In order to perform beam emittance measurements, we designed and fabricated a magnetic lens to be placed between the mesh and the screen. By measuring the size of the beam for different distances between the lens and the screen, we estimate the beam’s emittance. Here, we also simulated the electron beam dynamics with Computer Simulation Technology Studio and General Particle Tracer codes. Experimentally measured divergence angles and normalized transverse emittance show good agreement with simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2001

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM (atmospheric radiation measurement) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth observing system) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of 'radiative smoothing' of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. 'Offbeam' Lidar returns are now being used to estimate the cloud 'radiative Green's function' (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the 'Independent Pixel Approximation' or IPA. This and other measurements of 3D (three dimensional) transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international 'Intercomparison of 3D Radiation Codes' (I3RC) program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

DOE OSTI2019

The Effects of Radiation and Emitted Light Transport on the Positional Response of 11 cm × 42.5 cm × 5.5 cm NaI(Tl) Detectors

Experiments were performed with 30 11 cm x 42.5 x cm 5.5 cm NaI(Tl) detectors to better understand their positional response. Spectra were collected using 0.02 to 0.15 MBq point sources of 241 Am, 137 Cs, 60 Co, and 133 Ba positioned on lines parallel and perpendicular to the long axis of the crystal along both the narrow and wide detector faces as well as at different distances from them. A greater density of positions was sampled at the ends of the detector, and repeated measurements were made to examine potential gain drifts during the experiment. Spectroscopic peak counts, spectroscopic pulse heights, and net counts were analyzed. Empirical equations were fit to the aforementioned data for each specific source energy as a function of source position. In addition, a Monte Carlo radiation transport code was used to simulate the expected positionally variable response based solely upon radiation absorption. The simulated radiation transport efficiency functions were compared to the experimental data. The effects of the geometric radiation efficiency, the attenuation and scattering of emitted light within the scintillation crystal, and combined effects such as non-uniformity of the photomultiplier tube photocathode response and crystal irregularities were then distinguished. Functions describing each effect were derived. The results suggest potential new corrections to data obtained with large scintillation detectors as well as a novel approach to partial positional gamma-ray detection with minimal collimation, given that the energy resolution is within reason for particular photopeaks.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

An open source digital servo for atomic, molecular, and optical physics experiments

We describe a general purpose digital servo optimized for feedback control of lasers in atomic, molecular, and optical physics experiments. The servo is capable of feedback bandwidths up to roughly 1 MHz (limited by the 320 ns total latency); loop filter shapes up to fifth order; multiple-input, multiple-output control; and automatic lock acquisition. The configuration of the servo is controlled via a graphical user interface, which also provides a rudimentary software oscilloscope and tools for measurement of system transfer functions. We illustrate the functionality of the digital servo by describing its use in two example scenarios: frequency control of the laser used to probe the narrow clock transition of {sup 27}Al{sup +} in an optical atomic clock, and length control of a cavity used for resonant frequency doubling of a laser.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Ptychography at all wavelengths

Ptychography is a computational imaging technique that operates across multiple wavelength regimes, from electron (picometres) to X-ray (~0.1 nm), extreme ultraviolet (~10 nm) and visible light (micrometres). By reconstructing both amplitude and phase from diffraction patterns, ptychography enables high-resolution, quantitative imaging without conventional limitations imposed by lens-based optics. Ptychography has enabled advances across a range of scales: achieving deep-sub-angstrom resolution with electron microscopy, becoming an indispensable tool at X-ray synchrotron facilities worldwide and overcoming the trade-offs between resolution and field-of-view in optical imaging. This Primer provides a unified treatment of ptychography across these wavelength regimes. First, we discuss theoretical foundations, reconstruction algorithms, experimental considerations and wavelength-specific challenges. We then give examples of raw and processed data from various configurations and wavelengths. Next, we highlight key applications of ptychography in life sciences, materials science and industry. We also discuss data standards, open-source software implementations and best practices for ensuring reproducibility across different wavelength regimes. Finally, we consider limitations and future opportunities for ptychography. Together with accompanying datasets and code implementations, this Primer aims to serve newcomers and experienced practitioners in the field, facilitating broader adoption of ptychography across different disciplines.

47 OTHER INSTRUMENTATION↗

DOE OSTI2015

WE-AB-BRB-07: Alanine and Monte Carlo Determined Beam Quality Corrections for Nonstandard Fields of the Varian TrueBeam Accelerator

Purpose: To determine beam quality correction factors for flattened and unflattened beams of the Varian TrueBeam™ accelerator using alanine measurements and Monte Carlo calculations. Methods: Measurements were performed using L-α-alanine pellets encased in Virtual Water™ paddles. These were irradiated in liquid water using various field sizes of the 6MV and 10MV beam energies of the TrueBeam™, with both flattening-filter-free (FFF) and flattened beam modes. Measurements were also performed in a 10×10 cm{sup 2} {sup 60}Co field. Pellets were read out by the National Physical Laboratory (NPL, Teddington, UK) using electron paramagnetic resonance (EPR) spectrometry, and results were corrected for energy dependence using graphite calorimetry. Beam quality correction factors (k{sub Q}, k{sub Q}msr,Q, and k{sub Q}clin,Q{sub msr}) were determined for three ionization chambers: an Exradin A12 Farmer-type chamber, an Exradin A1SL scanning chamber, and an Exradin A26 microchamber. All chambers were modeled in the egs-chamber user code of EGSnrc. Sources were created using Varian-supplied IAEA-compliant phase spaces and the BEAMnrc user code, and were validated by comparing measured and simulated dose profiles. Results: Measured and calculated k{sub Q} values agreed within uncertainties, showing the chamber models to be reliable. A comparison of measured and calculated k{sub Q}msr,Q results with TG-51-based values showed that TG-51 adequately accounts for variations in beam quality between flattened and unflattened 10×10 cm{sup 2} fields. Small field corrections, k{sub Q}clin,Q{sub msr}, were determined to be up to 5%, referenced to a 10×10 cm{sup 2} field of the same energy and mode. Conclusion: Beam quality corrections were determined for several beam energies of the TrueBeam™ accelerator using Monte Carlo calculations and were validated using measurements with alanine. Values of k{sub Q} determined using TG-51 were found to be adequate for FFF reference fields. However, small field corrections are necessary for both flattened and unflattened beams to ensure greater treatment planning accuracy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1994

Accommodating Heterogeneity in a Debugger for Distributed Computations

In an ongoing project at NASA Ames Research Center, we are building debugger for distributed computations running on a heterogeneous set of machines. Historically, such debuggers have been built as front-ends to existing source-level debuggers on the target platforms. In effect, these back-end debuggers are providing a collection of debugger services to a client. The major drawback is that because of inconsistencies among the back-end debuggers, the front-end must use a different protocol when talking to each back-end debugger. This can make the front-end quite complex. We have avoided this complexity problem by defining the client-server debugger protocol. While it does require vendors to adapt their existing debugger code to meet the protocol, vendors are generally interested in doing so because the approach has several advantages. In addition to solving the heterogenous platform debugging problem, it will be pos~sible to write interesting debugger user interfaces that can be easily ported across a variety of machines. This will likely encourage investment in application-domain specific debuggers. In fact, the user interface of our debugger will be geared to scientists developing computational fluid dynamics codes. This paper describes some of the problems encountered in developing a portable debugger for heterogenous, distributed computing and how the architecture of our debugger avoids them. It then provides a detailed description of the debugger client-server protocol. Some of the more interesting attributes of the protocol are: (1) It is object-oriented; (2) It uses callback functions to capture the asynchronous nature of debugging in a procedural fashion; (3) It contains abstractions, such as in-line instrumentation, for the debugging of computationally intensive programs; (4) For remote debugging, it has operations that enable the implementor to optimize message passing traffic between client and server. The soundness of the protocol is being tested through the implementation of a prototype client and server pair.

Hood, Robert↗

DOE OSTI2017

(U) An Analytic Examination of Piezoelectric Ejecta Mass Measurements

Ongoing efforts to validate a Richtmyer-Meshkov instability (RMI) based ejecta source model [1, 2, 3] in LANL ASC codes use ejecta areal masses derived from piezoelectric sensor data [4, 5, 6]. However, the standard technique for inferring masses from sensor voltages implicitly assumes instantaneous ejecta creation [7], which is not a feature of the RMI source model. To investigate the impact of this discrepancy, we define separate “areal mass functions” (AMFs) at the source and sensor in terms of typically unknown distribution functions for the ejecta particles, and derive an analytic relationship between them. Then, for the case of single-shock ejection into vacuum, we use the AMFs to compare the analytic (or “true”) accumulated mass at the sensor with the value that would be inferred from piezoelectric voltage measurements. We confirm the inferred mass is correct when creation is instantaneous, and furthermore prove that when creation is not instantaneous, the inferred values will always overestimate the true mass. Finally, we derive an upper bound for the error imposed on a perfect system by the assumption of instantaneous ejecta creation. When applied to shots in the published literature, this bound is frequently less than several percent. Errors exceeding 15% may require velocities or timescales at odds with experimental observations.

47 OTHER INSTRUMENTATION↗

DOE OSTI2016

Photon Next-Event Estimator Implementation in MCATK

A photon next-event fluence estimator at a point has been implemented in the Monte Carlo Application Toolkit (MCATK). The next-event estimator provides an expected value estimator for the flux at a point due to all source and collision events. An advantage of the next-event estimator over track-length estimators, which are normally employed in MCATK, is that flux estimates can be made in locations that have no random walk particle tracks. The next-event estimator allows users to calculate radiographs and estimate response for detectors outside of the modeled geometry. The next-event estimator is not yet accessable through the MCATK FlatAPI for C and Fortran. The next-event estimator in MCATK has been tested against MCNP6 using 5 suites of test problems. No issues were found in the MCATK implementation. One issue was found in the exclusion radius approximation in MCNP6. The theory, implementation, and testing are described in this document.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2016

MAGNETO-FRICTIONAL MODELING OF CORONAL NONLINEAR FORCE-FREE FIELDS. I. TESTING WITH ANALYTIC SOLUTIONS

We report our implementation of the magneto-frictional method in the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The method aims at applications where local adaptive mesh refinement (AMR) is essential to make follow-up dynamical modeling affordable. We quantify its performance in both domain-decomposed uniform grids and block-adaptive AMR computations, using all frequently employed force-free, divergence-free, and other vector comparison metrics. As test cases, we revisit the semi-analytic solution of Low and Lou in both Cartesian and spherical geometries, along with the topologically challenging Titov–Démoulin model. We compare different combinations of spatial and temporal discretizations, and find that the fourth-order central difference with a local Lax–Friedrichs dissipation term in a single-step marching scheme is an optimal combination. The initial condition is provided by the potential field, which is the potential field source surface model in spherical geometry. Various boundary conditions are adopted, ranging from fully prescribed cases where all boundaries are assigned with the semi-analytic models, to solar-like cases where only the magnetic field at the bottom is known. Our results demonstrate that all the metrics compare favorably to previous works in both Cartesian and spherical coordinates. Cases with several AMR levels perform in accordance with their effective resolutions. The magneto-frictional method in MPI-AMRVAC allows us to model a region of interest with high spatial resolution and large field of view simultaneously, as required by observation-constrained extrapolations using vector data provided with modern instruments. The applications of the magneto-frictional method to observations are shown in an accompanying paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

IGRINS RV: A Precision Radial Velocity Pipeline for IGRINS Using Modified Forward Modeling in the Near-infrared

Application of the radial velocity (RV) technique in the near-infrared is valuable because of the diminished impact of stellar activity at longer wavelengths, making it particularly advantageous for the study of late-type stars but also for solar-type objects. In this paper, we present the IGRINS RV open-source python pipeline for computing infrared RV measurements from reduced spectra taken with IGRINS, an R ≡ λ/Δλ ∼ 45,000 spectrograph with simultaneous coverage of the H band (1.49–1.80 μm) and K band (1.96–2.46 μm). Using a modified forward-modeling technique, we construct high-resolution telluric templates from A0 standard observations on a nightly basis to provide a source of common-path wavelength calibration while mitigating the need to mask or correct for telluric absorption. Telluric standard observations are also used to model the variations in instrumental resolution across the detector, including a yearlong period when the K band was defocused. Without any additional instrument hardware, such as a gas cell or laser frequency comb, we are able to achieve precisions of 26.8 m s{sup −1} in the K band and 31.1 m s{sup −1} in the H band for narrow-line hosts. These precisions are empirically determined by a monitoring campaign of two RV standard stars, as well as the successful retrieval of planet-induced RV signals for both HD 189733 and τ Boo A; furthermore, our results affirm the presence of the Rossiter–McLaughlin effect for HD 189733. The IGRINS RV pipeline extends another important science capability to IGRINS, with publicly available software designed for widespread use.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Correlated Signal Analysis for Nuclear Emergency Response

In the context of nuclear emergency response (NER) scenarios, the critical task is to quickly identify a "black box" as a potential threat, as failing to do so could have catastrophic consequences. Techniques for passive assay of a “black box” typically include gamma-ray spectroscopy and neutron coincidence/multiplicity counting. However, there are significant challenges associated with these type of measurements. First, the presence of intervening materials can obstruct the detection of relevant signatures. Second, the presence of strong non-fission neutron sources, like (α, n) emitters, can add uncertainties to the neutron multiplicity analysis. In our LDRD-MFR Phase II work, a portable neutron spectrometer, called the Compact Fast Neutron Spectrometer (CFNS), was developed for NER applications. The CFNS system leverages information-rich neutron energy spectra to derive actionable information. This work investigates the use of correlated signals in the CFNS from special nuclear material (SNM) to characterize physical properties, such as intervening shielding material and fission to non-fission neutron contributions. In this report, the Phase II results from bulk SNM measurements at the National Criticality Experiments Research Center (NCERC) are briefly discussed along with the motivation for this work. Afterwards, simulations using the MCNPX-PoliMi transport code are discussed, which were used to expand our correlated signal study. Signal triggered analysis for neutron multiplicity extraction will also be discussed. Lastly, the use neutron-photon correlations for intervening material identification are shown, along with the use of correlated neutron energy spectra for α-ratio extraction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2017

astroABC : An Approximate Bayesian Computation Sequential Monte Carlo sampler for cosmological parameter estimation

Given the complexity of modern cosmological parameter inference where we arefaced with non-Gaussian data and noise, correlated systematics and multi-probecorrelated data sets, the Approximate Bayesian Computation (ABC) method is apromising alternative to traditional Markov Chain Monte Carlo approaches in thecase where the Likelihood is intractable or unknown. The ABC method is called"Likelihood free" as it avoids explicit evaluation of the Likelihood by using aforward model simulation of the data which can include systematics. Weintroduce astroABC, an open source ABC Sequential Monte Carlo (SMC) sampler forparameter estimation. A key challenge in astrophysics is the efficient use oflarge multi-probe datasets to constrain high dimensional, possibly correlatedparameter spaces. With this in mind astroABC allows for massive parallelizationusing MPI, a framework that handles spawning of jobs across multiple nodes. Akey new feature of astroABC is the ability to create MPI groups with differentcommunicators, one for the sampler and several others for the forward modelsimulation, which speeds up sampling time considerably. For smaller jobs thePython multiprocessing option is also available. Other key features include: aSequential Monte Carlo sampler, a method for iteratively adapting tolerancelevels, local covariance estimate using scikit-learn's KDTree, modules forspecifying optimal covariance matrix for a component-wise or multivariatenormal perturbation kernel, output and restart files are backed up everyiteration, user defined metric and simulation methods, a module for specifyingheterogeneous parameter priors including non-standard prior PDFs, a module forspecifying a constant, linear, log or exponential tolerance level,well-documented examples and sample scripts. This code is hosted online athttps://github.com/EliseJ/astroABC

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗