Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “algorithm development”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Demonstration of Pronghorn’s Subchannel Code Modeling of Liquid-Metal Reactors and Validation in Normal Operation Conditions and Blockage Scenarios

Pronghorn-SC is a subchannel code within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Initially designed to simulate flows in water-cooled, square lattice, subchannel assemblies, Pronghorn-SC has been expanded to simulate liquid-metal-cooled flows in triangular lattices, hexagonal subchannel assemblies. For this purpose, the algorithm of Pronghorn-SC was adapted to solve the subchannel equations as they are applicable to a hexagonal wire-wrapped sodium-cooled fast reactor. Cheng–Todreas models for pressure drop and cross-flow models were adopted and a coolant heat conduction term was added. To solve these equations, an improved implicit algorithm was developed robust enough to deal with the numerical issues, associated with low flow and recirculation phenomena. To confirm the prediction capability of Pronghorn-SC, calculations and comparisons with available experimental data of 19- and 37-pin assemblies were performed, as well as other subchannel codes. Finally, a flow blockage modeling feature was added. This capability was validated for both water-cooled square sub-assemblies and sodium-cooled hexagonal sub-assemblies, using experimental data of partially and fully blocked cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Series FACTS Devices for Increasing Resiliency in Severe Weather Conditions

Severe weather conditions are low-probability, high-impact events that affect grid operations. The majority of power outages are caused by severe weather conditions. Grid resiliency to weather events can be enhanced by decreasing the reliance on its affected sections. One way to do this is to reduce the power flow through lines vulnerable to severe weather. If a line is disconnected, its initial power flow is distributed through the neighbor lines, which may cause congestion in the grid. FACTS devices can be used to control the power flow of lines that have a higher chance of power outages. Most previous works do not consider weather events in power flow control. In this work, a linearized optimal power flow (OPF)–based algorithm is developed to minimize the real power flow of vulnerable lines considering the thermal limits of lines to prevent infeasible solutions; the simulation is fast, making it suitable for large-scale systems. The proposed optimization problem is presented as a mixed-integer linear program (MILP), making it capable of using short-term load forecasting due to its high solution speed. The proposed optimization problem considers multiple lines with different outage probabilities and the uncertainties of the weather forecast. Moreover, it estimates the power reduction in vulnerable lines due to changes in the series FACTS devices. The performance of the proposed optimization problem is tested on IEEE 14-, 30-, and 118-bus systems for several scenarios. The results are validated with the AC power flow results from MATPOWER.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Seeking Echoes of Circumstellar Disks in Kepler Light Curves

Light echoes of flares on active stars offer the opportunity for direct detection of circumstellar dust. We revisit the problem of identifying faint echoes in postflare light curves, focusing on debris disks from ongoing planet formation. Starting with simulations, we develop an algorithm for estimating the radial extent and total mass from disk echo profiles. We apply this algorithm to light curves from over 2100 stars observed by NASA’s Kepler mission, selected for multiple short-lived flares in either the long-cadence or short-cadence data sets. While flux uncertainties in light curves from individual stars preclude useful mass limits on circumstellar disks, catalog-averaged light curves yield constraints on disk mass that are comparable to estimates from known debris disks. The average mass in micron- to millimeter-sized dust around the Kepler stars cannot exceed 10% of an Earth mass in exo-Kuiper belts or 10% of a lunar mass in the terrestrial zone. We group stars according to IR excess, based on WISE W1–W3 color, as an indicator for the presence of circumstellar dust. The mass limits are greater for stars with strong IR excess, a hint that echoes are lurking not far beneath the noise in postflare light curves. With increased sensitivity, echo detection will let time-domain astronomy complement spectroscopic and direct-imaging studies in mapping how, when, and where planets form.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

Hakobyan, Hayk [Flatiron Institute, New York, NY (↗
DOE OSTI2023

Long-range transported continental aerosol in the eastern North Atlantic: three multiday event regimes influence cloud condensation nuclei

Abstract. The eastern North Atlantic (ENA) is a region dominated by pristine marine environment and subtropical marine boundary layer clouds. Under unperturbed atmospheric conditions, the regional aerosol regime in the ENA varies seasonally due to different seasonal surface-ocean biogenic emissions, removal processes, and meteorological regimes. However, during periods when the marine boundary layer aerosol in the ENA is impacted by particles transported from continental sources, aerosol properties within the marine boundary layer change significantly, affecting the concentration of cloud condensation nuclei (CCN). Here, we investigate the impact of long-range transported continental aerosol on the regional aerosol regime in the ENA using data collected at the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility on Graciosa Island in 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. We develop an algorithm that integrates number concentrations of particles with optical particle dry diameter (Dp) between 100 and 1000 nm, single scattering albedo, and black carbon concentration to identify multiday events (with duration >24 consecutive hours) of long-range continental aerosol transport in the ENA. In 2017, we detected nine multiday events of long-range transported particles that correspond to ∼ 7.5 % of the year. For each event, we perform HYSPLIT 10 d backward trajectories analysis, and we evaluate CALIPSO aerosol products to assess, respectively, the origins and compositions of aerosol particles arriving at the ENA site. Subsequently, we group the events into three categories, (1) mixture of dust and marine aerosols, (2) mixture of marine and polluted continental aerosols from industrialized areas, and (3) biomass burning aerosol from North America and Canada, and we evaluate their influence on aerosol population and cloud condensation nuclei in terms of potential activation fraction and concentrations at supersaturation of 0.1 % and 0.2 %. The arrival of plumes dominated by the mixture of dust and marine aerosol in the ENA in the winter caused significant increases in baseline Ntot. Simultaneously, the baseline particle size modes and CCN potential activation fraction remained almost unvaried, while cloud condensation nuclei concentrations increased proportionally to Ntot. Events dominated by a mixture of marine and polluted continental aerosols in spring, fall, and winter led to a statistically significant increase in baseline Ntot, a shift towards larger particular sizes, a higher CCN potential activation fractions, and cloud condensation nuclei concentrations of >170 % and up to 240 % higher than during baseline regime. Finally, the transported aerosol plumes characterized by elevated concentration of biomass burning aerosol from continental wildfires detected in the summertime did not statistically contribute to increase baseline aerosol particle concentrations in the ENA. However, particle diameters were larger than under baseline conditions, and CCN potential activation fractions were >75 % higher. Consequentially, cloud concentration nuclei concentrations increased by ∼ 115 % during the period affected by the biomass burning events. Our results suggest that, through the year, multiday events of long-range continental aerosol transport periodically affect the ENA and represent a significant source of CCN in the marine boundary layer. Based on our analysis, in 2017, the multiday aerosol plume transport dominated by a mixture of dust and marine aerosol, a mixture of marine and polluted continental aerosols, and biomass burning aerosols caused increases in the NCCN baseline regime of, respectively, 6.6 %, 8 %, and 7.4 % at SS 0.1 % (and, respectively, 6.5 %, 8.2 %, and 7.3 % at SS 0.2 %) in the ENA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Integration of Total-Sky Imager Data with a Physics-Based Smart Persistence Model for Intra-Hour Forecasting of Solar Radiation

Short-term solar forecasting models based solely on global horizontal irradiance (GHI) measurements are often unable to discriminate the forecasting of the factors affecting GHI from those that can be precisely computed by atmospheric models. Our previous study introduced a Physics-based Smart Persistence model for Intra-hour forecasting of solar radiation (PSPI) that decomposed the forecasting of GHI into the computation of extraterrestrial solar radiation and solar zenith angle and the forecasting of cloud albedo and cloud fraction. The extraterrestrial solar radiation and solar zenith angle were accurately computed by the Solar Position Algorithm (SPA) developed at the National Renewable Energy Laboratory (NREL). A cloud retrieval technique was used to estimate cloud albedo and cloud fraction from surface-based observations of GHI. With the assumption of persistent cloud structures, the cloud albedo and cloud fraction were predicted for future time steps using a two-stream approximation and a 5-minute exponential weighted moving average, respectively. The model evaluation indicated the estimation and forecast of cloud fraction mostly contributed to the uncertainty of the PSPI though it overcame the persistence and smart persistence models in all forecast time horizons between 5 and 60 minutes. This study aims to enhance the PSPI by ingesting surface-based observations of cloud fraction from a total sky imager (TSI). The estimation and forecast of cloud albedo is correspondingly improved by utilizing the cloud fraction observations and thus leads to more accurate GHI forecast. Various time-series analysis methods are also investigated on the forecasting of cloud fraction and cloud albedo for further improving the GHI forecast. These improvements are valuable for many applications, such as forecasting energy use for buildings, grid operations, and ultimately bringing down the cost of solar energy.

14 SOLAR ENERGY↗
DOE OSTI2020

Dynamic Restoration Strategy for Distribution System Resilience Enhancement

In electric power distribution systems, distributed energy resources (DERs) can act as controllable power sources and support utility operators to minimize power outages after extreme weather events (e.g., hurricane, earthquake, wildfire) and thus help enhance the grid's resilience. Meanwhile, the influences of extreme events and the capabilities of DERs are dynamic and difficult to predict. Hence, the desired distribution system restoration strategy should be able to evolve according to real-time fault/disturbance information and the availability of DERs. In this paper, we propose a new dynamic restoration strategy for distribution systems to enhance system resilience against potential hazards. An efficient reconfiguration algorithm is developed to eliminate the use of integer variables to relieve the computational burden. Model predictive control is implemented to adjust the system topology and DER operation set points based on the updated fault information and DER forecasts. The effectiveness of the proposed restoration model in enhancing distribution system resilience is validated through an IEEE 123-bus test system. Simulation results also validate that the proposed restoration model can mitigate the occurrence of unexpected events and the fluctuations of DERs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Model-Free Primal-Dual Methods for Network Optimization with Application to Real-Time Optimal Power Flow: Preprint

This paper examines the problem of real-time optimization of networked systems and develops online algorithms that steer the system towards the optimal trajectory without explicit knowledge of the system model. The problem is modeled as a dynamic optimization problem with time-varying performance objectives and engineering constraints. The design of the algorithms leverages the online zero-order primal-dual projected-gradient method. In particular, the primal step that involves the gradient of the objective function (and hence requires networked systems model) is replaced by its zero-order approximation with two function evaluations using a deterministic perturbation signal. The evaluations are performed using the measurements of the system output, hence giving rise to a feedback interconnection, with the optimization algorithm serving as a feedback controller. The paper provides some insights on the stability and tracking properties of this interconnection. Finally, the paper applies this methodology to a real-time optimal power flow problem in power systems, and shows its efficacy on the IEEE 37-node distribution test feeder for reference power tracking and voltage regulation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Progress on the National Solar Radiation Data Base (NSRDB): A new DNI computation

This study introduces a new technique to compute direct normal irradiance (DNI) for improving the National Solar Radiation Data Base (NSRDB). A finite-surface integration algorithm is developed to compute solar radiation in differential solid angles and efficiently infer its contribution to a surface perpendicular to the solar direction. A lookup table of cloud bi-directional transmittance distribution function (BTDF) is developed by use of the discrete ordinates radiative transfer (DISORT) model for possible solar and observing directions and various cloud optical and microphysical properties. In each solar incident direction, DNI is given by the cloud BTDFs from approximately 200 differential solid angles. The simulated DNI is calibrated and evaluated using surface observations at the National Renewable Energy Laboratory's (NREL's) Solar Energy Research Laboratory (SRRL) and the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) facility.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

High-Fidelity Heavy-Duty Vehicle Modeling Using Sparse Telematics Data

Heavy-duty commercial vehicles consume a significant amount of energy due to their large size and mass, directly leading to vehicle operators prioritizing energy efficiency to reduce operational costs and comply with environmental regulations. One tool that can be used for the evaluation of energy efficiency in heavy-duty vehicles is the evaluation of energy efficiency using vehicle modeling and simulation. Simulation provides a path for energy efficiency improvement by allowing rapid experimentation of different vehicle characteristics on fuel consumption without the need for costly physical prototyping. The research presented in this paper focuses on using real-world, sparsely sampled telematics data from a large fleet of heavy-duty vehicles to create high-fidelity models for simulation. Samples in the telematics dataset are collected sporadically, resulting in sparse data with an infrequent and irregular sampling rate. Captured in the dataset was geospatial information, time series measurements, and vehicle-specific metadata from a subset of 96 vehicles from varied geographic regions across North America. A series of custom algorithms was developed to process vehicle data and derive both vehicle model input parameters and representative drive cycles. Derived models provide a basis on which to simulate real-world vehicles and iterate on vehicle aerodynamics, auxiliary power loads, transmission shift schedules, and other parameters to achieve reduced fuel consumption and increase energy efficiency. Notably, these models were developed without the use of expensive field data collection, using only data collected through fleet telematics. Processed representative drive cycles are used to validate the fuel economy of derived models. The models developed through this research allow for more representative vehicle simulations with increased flexibility regarding vehicle-to-vehicle variations.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Towards a neutron multiplicity measurement with the Accelerator Neutrino Neutron Interaction Experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26 ton Gadolinium (Gd)-loaded water Cherenkov detector located on the Booster Neutrino Beam line at Fermilab. Its main goals are the measurement of the neutron multiplicity in neutrino-nucleus interactions as well as the cross-section of Charged Current Quasi-Elastic (CCQE) neutrino interactions on water. Besides the physics goals, the experiment also aims to be a testbed for new technologies such as Large Area Picosecond Photodetectors (LAPPDs) and Water-based Liquid Scintillators (WbLS). This thesis presents a preliminary measurement of the neutron multiplicity with {ANNIE}, using an analysis conducted on a fraction of the 2021 beam year. As preparatory measures, the efficiency of ANNIE's Front Muon Veto (FMV) was determined to be {$\bar{\varepsilon}_{\mathrm{FMV}} = (95.6 \pm 1.6)\%$} while the average efficiency for active scintillator paddles in the Muon Range Detector (MRD) was found to be {$\bar{\varepsilon}_ {\mathrm{MRD}} = (92.1 \pm 7.9)\%$}. Furthermore, the simulation framework used for ANNIE was validated and adapted to reproduce the experimental data by comparing the detector response for samples of Michel electrons, Americium Beryllium neutrons, and through-going muons. The analysis finds average neutron yields of {$\bar{n}_{\mathrm{data}} (\mathrm{beam}) = (0.272 \pm 0.010_{\mathrm{stat}})$} for an inclusive set of all identified muon neutrino candidates and {$\bar{n}_{\mathrm{data}} (\mathrm{beam,FV}) = (0.287 \pm 0.044_{\mathrm{stat}})$ for interactions which happened inside of the Fiducial Volume of ANNIE, which was optimized to increase the neutron detection acceptance. The presented neutron multiplicity values represent the number of detected neutrons after all event selection cuts and are not yet corrected for the neutron detection efficiency. An equivalent analysis on a simulated beam sample predicts neutron yields of $\bar{n}_{\mathrm{MC}}(\mathrm{beam}) = (0.515 \pm 0.0 07_{\mathrm{stat}})$ and $\bar{n}_{\mathrm{MC}}(\mathrm{beam,FV}) = (0.627 \pm 0.031_{\mathrm{stat}})$, indicating that the models tend to overpredict the number of neutrons produced in such interactions. Systematic errors have been briefly considered to contribute {$\sigma_{\mathrm{sys,FMV}} \sim 0.01\,$neutrons/$\nu$-interaction} due to the slight FMV inefficiency and {$\sigma_{\mathrm{sys,n}} \sim 0.05\,$neutrons/$\nu$-interaction} due to the neutron detection efficiency. Simulation studies further highlighted the importance of neutron detection in Diffuse Supernova Background (DSNB) searches. A combination of neutron tagging and Convolutional Neural Networks was found to reduce the most relevant Neutral Current Quasi-Elastic (NCQE) interaction background below the signal level, achieving a Signal-to-Background ratio of 4:1. In a further study, we investigated the positive impact of a deployment of a WbLS target on the energy reconstruction in ANNIE. WbLS provides a scintillation signal from hadronic recoils in addition to the charged lepton that can be included in neutrino energy reconstruction. It was found that a deployed WbLS volume in ANNIE improves the neutrino energy reconstruction from 14\% to 12\%, with the potential of going beyond this if more sophisticated reconstruction algorithms are developed in the future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Expanding Accurate Person Recognition to New Altitudes and Ranges: The BRIAR Dataset

Face recognition technology has advanced significantly in recent years due largely to the availability of large and increasingly complex training datasets for use in deep learning models. These datasets, however, typically comprise images scraped from news sites or social media plat-forms and, therefore, have limited utility in more advanced security, forensics, and military applications. These applications require lower resolution, longer ranges, and ele-vated viewpoints. To meet these critical needs, we collected and curated the first and second subsets of a large multi-modal biometric dataset designed for use in the research and development (R&D) of biometric recognition technolo-gies under extremely challenging conditions. Thus far, the dataset includes more than 350,000 still images and over 1,300 hours of video footage of approximately 1,000 sub-jects. To collect this data, we used Nikon DSLR cameras, a variety of commercial surveillance cameras, specialized long-rage R&D cameras, and Group 1 and Group 2 UAV platforms. The goal is to support the development of algorithms capable of accurately recognizing people at ranges up to 1,000 m and from high angles of elevation. These ad-vances will include improvements to the state of the art in face recognition and will support new research in the area of whole-body recognition using methods based on gait and anthropometry. This paper describes methods used to col-lect and curate the dataset, and the dataset's characteristics at the current stage.

Brogan, Joel↗
DOE OSTI2023

A Lookup Table Based Low Resource Usage Waveform Data Compression Block in FPGA

aA lookup table based lossless data compression scheme suitable for FPGA implementation, with very low silicon resource usage, is developed. Different algorithms can be implemented by simply changing the contents of the lookup tables in the functional block while keeping the logic circuits unchanged. It exhibits performance comparable to the zip compression software in Linux distributions. The compression block can fit into almost any existing FPGA used in the online data acquisition systems of high energy physics experiments. The decompression process can be performed in either software or in FPGA firmware, allowing users to process compressed data flexibly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Optimisation of the Search for CP-symmetry Violation at the Deep Underground Neutrino Experiment

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long baseline experiment, which will be situated in South Dakota. Its detectors will utilise liquid-argon time projection chamber technology, which is able to capture neutrino interactions with an incredible spatial and calorimetric resolution. With what will become the world’s most intense neutrino beam, a highly capable near detector, and four (10kt fducial mass) far detector modules, DUNE will be able to achieve an ambitious physics programme. Most notably, DUNE will determine whether charge-parity symmetry is broken in neutrino oscillations - a finding that would have significant implications for the understanding of the matter-antimatter asymmetry in our Universe.This thesis presents the optimisation of a CP-violation analysis at DUNE using thePandora pattern-recognition software. The analysis assumes a 3.5 year exposure (1.36 × 1023 protons on target) to a neutrino and an antineutrino beam (7 year total). Only the predicted data of the far detector modules is used; near detector samples are not included. The initial sensitivity to CP-violation is found to be 3.8σ+0.9σ−1.1σ in an estimate that includes oscillation parameter uncertainties, systematic uncertainties and statistical fluctuations, and assumes a normal-ordering of the neutrino mass hierarchy. The performance of the Pandora event reconstruction is linked to that of the analysis, which is found to be limited by the reconstruction of the initial track-like region of electrons and photons. The ShowerRefinement algorithm is developed in response to this and its implementation into the analysis workflow results in an improved sensitivity to CP-violation of 4.6σ+0.9σ−1.0σ. With a perfected neutrino interaction vertex placement, this is further increased to 5.1σ+1.0σ−1.1σ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Measurement of the pion charge exchange differential cross section on Argon with the ProtoDUNE-SP detector

The Deep Underground Neutrino Experiment (DUNE) utilizes Liquid Argon Time Projection Chamber (LArTPC) technology in its pursuit of uncovering the origins of matter. As a crucial prototype for the DUNE far detector (FD) module, ProtoDUNE-SP at CERN Neutrino Platform Facility employs a hadron test beam with a momentum range of 0.3 to 7 GeV/c to study the detector response to particles that could be produced in the neutrino interactions at the DUNE FD. Furthermore, ProtoDUNE-SP is dedicated to testing event reconstruction and performing detector calibration under controlled conditions. Neutral pion reconstruction performance is an important benchmark of the calorimetric and tracking capability of LArTPC, and understanding neutral pion in LArTPC is crucial to characterize backgrounds to oscillation measurements and rare searches. To improve neutral pion reconstruction in LArTPC, a kinematic fitting algorithm is developed, achieving a 12% resolution in neutral pion energy in ProtoDUNE-SP. This enhancement enables the first measurements of the pion-argon charge exchange differential cross-section. By comparing these results with the current model, valuable constraints can be made to improve generator simulations in the future DUNE experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Big Data and AI at DoE's Legacy Sites - 20546

More than 30 years have passed since DOE started the decommissioning of nuclear weapon complexes and the clean-up of soil and groundwater. All the sites have been collecting and archiving soil and groundwater monitoring datasets; particularly contaminant concentration time-series. These datasets provide unparalleled opportunities to understand the system behavior (including more fundamental hydrological and geochemical processes, the response to various perturbations, the long-term trend and environmental decay rate towards the regulatory limit). This understanding is critical for providing multiple lines of evidences that can support site closure. In this study, we explore the machine learning (ML) and artificial intelligence (AI) applications to the long-term soil and groundwater management at DoE's legacy sites. ML can improve our understanding of the subsurface systems, which is critical for long-term monitoring and management of the sites, while AI can automate or support some of decision-making processes (e.g., anomaly detection, monitoring well placements). The particular focuses are to develop general algorithms to: (1) to identify distinct spatiotemporal patterns and to identify several groups that have similar temporal behaviors, using unsupervised clustering methods, (2) identify the different temporal scales of hydrological responses to climate perturbations by time-series analysis, and (3) reduce the number of monitoring wells by identifying the minimum sufficient number of wells to capture the heterogeneity of the groundwater contaminant plume and concentration distribution, using the Gaussian Process model. We demonstrate our methodology at the Savannah River Site F-Area. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Scaling and performance portability of the particle-in-cell scheme for plasma physics applications through mini-apps targeting exascale architectures

We perform a scaling and performance portability study of the particle-in-cell scheme for plasma physics applications through a set of mini-apps we name "Alpine", which can make use of exascale computing capabilities. The mini-apps are based on Independent Parallel Particle Layer, a framework that is designed around performance portable and dimension independent particles and fields. We benchmark the simulations with varying parameters such as grid resolutions (5123 to 20483) and number of simulation particles (109 to 1011) with the following mini-apps: weak and strong Landau damping, bump-on-tail and two-stream instabilities, and the dynamics of an electron bunch in a charge-neutral Penning trap. We show strong and weak scaling and analyze the performance of different components on several pre-exascale architectures such as Piz-Daint, Cori, Summit and Perlmutter. While the scaling and portability study helps identify the performance critical components of the particle-in-cell scheme in the current state-of-the-art computing architectures, the mini-apps by themselves can be used to develop new algorithms and optimize their high performance implementations targeting exascale architectures.

Muralikrishnan, Sriramkrishnan↗
DOE OSTI2022

Fully Distributed Multi-Material Magnetic Sensing Structures for Multiparameter DAS Applications

This dissertation demonstrates the first of its kind distributed magnetic field sensor based on a fiber optic distributed acoustic sensing (DAS) scheme. Ferromagnetic nickel and Metglas® were dispersed internally within a fiber optic preform and then drawn on an in-house fiber optic draw tower to lengths in the kilometers. Due to the close proximity of the ferromagnetic metals and fiber optic core, the magnetostrictive strain response of the ferromagnetic materials when exposed to a magnetic field would perturbate within the fiber cladding and transfer that strain, internally, to the fiber optic core. Strain resulting from the magnetostrictive effect allows the DAS based sensor to accurately translate strain into readable magnetic field data. Due to the high sensitivity seen in this sensor design, multiparameter sources, acoustic and magnetic fields, were tested and validated and a three dimensional magnetic-field vector sensor was proposed. Numerical analysis of the novel sensor design was first implemented using COMSOL Multiphysics, where inputs such as magnetostrictive element shape, size, distance, and number were first investigated. Upon optimizing system constraints, the sensor design was further modified such that single mode operation was consistent across multiple fiber draws while retaining high strain transfer from the ferromagnetic elements to the fiber optic core. Ferromagnetic material selection was evaluated as a function of the saturation magnetostriction constants and a total of 4 modules were used to fully characterize the complex physics involved in this sensor design. All fabrication and testing were performed in-house using a full scale 3-story fiber draw tower and custom environmental testing stations to imitate naturally occurring events such as magnetic or acoustic point sources. A unique stacking method was used to embed ferromagnetic nickel and Metglas® into a fiber optic preform which when combined with a custom fiber draw process resulted in consistent multi-material fibers drawn to lengths of 1-km. In-house testing facilities included different types of electromagnetic generators, in addition to a soil test bed, and an outdoor test bed which allowed 100 meters of fiber to be tested simultaneously. All tested sensors demonstrated high strain transfer capabilities on the order of 0.01-10 μϵ depending on the materials used, ferromagnetic rod number, and core to metal spacing. Due to the sensitivity of the system the difference between AC and DC was distinct, and directional magnetostriction was studied. Transverse and longitudinal magnetic wave propagation was controlled through a solenoid and rectangular Helmholtz coil, both built in-house. A three-dimensional magnetic field vector sensor was proposed due to the success of the magnetic field sensor, and a design was proposed and initially tested to validate direction as a function of field strength and distance. To summarize, this dissertation explores the first fully distributed magnetic field sensor using DAS based techniques and one of the first multi-material fiber draw processes which can produce consistent single mode fiber up to 1-km. Due to extensive FEA modeling, multiple iterations of the magnetic sensor were fully characterized and an equation describing the relationship between sensor design and strain transfer has been created and validated experimentally. Multi-parameter tests including acoustic and magnetic fields were implemented and an algorithm was developed to separate the mixed signals. Finally, a test was performed to demonstrate the feasibility of sensing magnetic fields directionally. Cumulative results demonstrate a high-quality sensor alternative to current designs which may surpass other magnetic sensors due to innate multi-parameter capabilities, in addition to the inexpensive production cost and extremely long operating lengths.

02 PETROLEUM↗