Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Task Performance and Analysis”

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,225 records · Page 68

Sequence2Self: Self-supervised image sequence denoising of pixel-level spray breakup morphology

Optical imaging of fast and transient phenomena such as the turbulent breakup of liquid sprays exhibit low signal-to-noise ratios due to the limited illumination intensity relative to the short exposure time. Image denoising is required to facilitate physical studies over these data but is challenging due to the absence of clean ground-truths and the stringency of the denoising task (e.g., strong and complex noise, limited resolution, preserving physical fidelity), preventing supervised and existing un-/self-supervised deep learning methods. To this end, Sequence2Self (Seq2S) is proposed, an extension of Self2Self (S2S) to image sequences that leverages both the signal’s spatial and temporal correlation. Seq2S is demonstrated on time-resolved x-ray phase contrast imaging of liquid jet fuel sprays in a gas turbine combustor, which possesses all of challenges detailed above. Experiments are conducted across four fuels with different breakup morphology using various state-of-the-art methods. Overall, many of the methods failed and Seq2S was most successful: (1) Accurate spray structures were reconstructed with consistent evolution across frames void of artifacts. (2) The performance was robust, invariant to the hyperparameter choice. (3) Computational time is short and can be made eligible for real-time denoising. In particular, the images denoised by Seq2S showed spray droplet diameter distributions with near-zero Kullback–Leibler divergence (0.01 ± 0.01) to a cleaner reference, whereas the second best method yielded 0.06 ± 0.03. In conclusion, this suggests that Seq2S can be reliably used prior to subsequent quantitative spray analyses as it retains (if not, improves) the statistical physical properties of the data.

97 MATHEMATICS AND COMPUTING↗

Applicability study of Bayesian optimization in core neutronic design using a toy model

At the Japan Atomic Energy Agency (JAEA), an innovative design approach named ARKADIA (Advanced Reactor Knowledge- and AI-aided Design Integration Approach through the whole plant life cycle) for advanced nuclear reactors is currently under development. One task in ARKADIA is to build a system that automatically optimizes core and fuel designs by conducting core neutronic and thermal-hydraulic calculations, fuel integrity evaluations, and plant dynamic analyses. This system will be implemented to automatically find an optimal design that minimizes (or maximizes) objective function defined by core performance while varying the core and fuel design parameters such as fuel pin diameter, core height and diameter. In this study, as the first step of system development, we focused only on core neutronic design and conducted a study of automatic optimization. As the optimization algorithm, Bayesian optimization (BO), an effective method for optimization problems with expensive computational cost of objective function, was utilized. The applicability of BO was studied based on single- and two-objective optimization examples of core neutronic design in a toy model. As a result, in the former, it was shown that BO can give the optimal solution, which matches the reference solution calculated by a brute force calculation well, with a small number of required calculations. Usability on core neutronic designs, where the computational cost per case is high, was confirmed. In the latter, it was found that BO can give a Pareto solutions-set that shows good agreement with the reference solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Laboratory Instrument Software Controlled Spread Spectrum Time Domain Reflectometry for Electrical Cable Testing

This research discusses development of a software-controlled laboratory instrument based spread spectrum time domain reflectometry system (SSTDR). This constitutes one task within PNNL’s Light Water Sustainability Program (LWRS) whose mission includes advancing nondestructive examination (NDE) techniques for off-line and on-line in-situ cable condition monitoring. In 2022, PNNL evaluated SSTDR for detection and characterization of a number of cable anomalies (Glass et al. 2022). The review included comparison of SSTDR to Frequency Domain Reflectometry (FDR) techniques which have enjoyed encouraging feedback and are starting to be used in nuclear power plants for periodic cable condition monitoring of cable systems as part of the plant’s overall cable aging management program. The FDR test introduces a broad-band chirp onto the cable at the cable end then listens for any reflection from a change of impedance along the cable caused by a damaged conductor or insulation, splices, contact with moisture, or other cable anomalies. The signal is captured in the frequency domain then transformed back to the time domain using an inverse Fourier transform (IFT). Based on the velocity of propagation, the impedance response signal is plotted against distance along the cable. Peak locations along the X-axis indicate the distance along the cable where a portion of the signal has been reflected back to the instrument as a result of a cable anomaly. The FDR test is considered the gold standard of reflectometry however it does require the cable to be de-energized to perform the test. The LIVEWIRE commercial SSTDR produces a similar plot to the FDR however all processing is in the time domain. A pseudo-random noise code (PN code) is input onto the cable conductor and the instrument listens for any reflected response from cable anomalies. The SSTDR processes the signal as an autocorrelation comparing the input PN code to any reflected signal detected. The autocorrelation analysis for thermal aging, water and water ingress detection, ground fault and phase-to-phase fault detection at various locations along the cable and with the cable attached and detached from a motor load, and on both energized and un-energized conditions were performed. These results were contrasted to Frequency Domain Reflectometry (FDR) measurements of the un-energized cable. Results were encouraging but indicated more work was warranted – particularly with the SSTDR, it seemed that the insulation damage would likely be better evaluated with multiple bandwidth cable tests particularly including larger bandwidths than were possible with the current commercial instrument. The commercial instrument’s bandwidth was set at 6, 12, 24, and 48MHz but note that SSTDR and FDR definitions of bandwidth trend similarly but are not the same. The FDR response could be more broadly adjusted, and the bandwidth of 100 to 500 MHz produced the best responses. FDR responses to anomalies were clearer than SSTDR responses and indications were that a broader bandwidth SSTDR may lead to improved SSTDR detection capability. This project used a laboratory instrument based SSTDR (primarily using an Arbitrary Waveform Generator (AWG) and a digital oscilloscope plus Python in-house software) that allowed software adjustment of the SSTDR bandwidth, window functions applied to the exciting Pseudo-random Noise (PN) code plus and other aspects of the SSTDR signal processing. Hereafter, this will be referred to as the PNNL SSTDR. Evaluating specific performance of the PNNL SSTDR is left to a separate report. This report documents hardware and software development to produce the SSTDR cable test system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Understanding the Effect of Na in Improving the Performance of CuInSe 2 Based Photovoltaics

Cu(In,Ga)Se 2 (CIGS) thin film photovoltaic technology is in the early stages of commercialization with an annual manufacturing capacity over 1 GW and has demonstrated the highest module efficiency of any of the thin film technologies. However there still is a lack of fundamental understanding of the relationship between the material properties and solar cell device operation. It is well known that the incorporation of a small amount of Na into the CIGS film during processing is essential for high efficiency devices. However, there are conflicting explanations for how Na behaves at the atomic scale. This report investigates how Na is incorporated into the CIGS device structure and evaluates the diffusion of Na into CIGS grain boundaries (GBs) and bulk crystallites. Participants: This project was carried out at the Institute of Energy Conversion at the University of Delaware, collaborating with the Rockett group at the University of Illinois Urbana-Champagne. Significant Findings: The significant outcomes of this project for each task include; Task 1.0: Effect of Na in Devices Fabricated on PVD Deposited CIGS; Na diffusion occurs through the Mo back contact via GBs driven by the presence of oxygen; Na reversibly compensates donor defects in CIGS GBs,Task 2.0: Na Incorporation in Single Crystal CIGS; and bulk Na diffusion proceeds rapidly such that grains are Na-saturated immediately following CIGS thin film manufacture. Industry Guidance: The presented results offer interesting concepts for modification of manufacturing processes of CIGS-based PV modules. Possible approaches to improve control of Na uptake and uniformly increase levels in CIGS films are highlighted for processes that employ either soda-lime glass or NaF as the Na source. Concepts include the potential of O 2 or oxidative based treatments of Mo back contacts to improve Na diffusion through the metal film and increase Na uptake into the growing CIGS. This project has also offered fundamental understanding of the behavior of Na in CIGS grains and GBs, particularly the confirmation that CIGS grains will be saturated with Na immediately following manufacture Summary of Results: Most commercially available CIGS modules are fabricated on soda-lime glass coated with Mo as the back electric contact, and Na in the glass diffuses through the Mo layer into the CIGS during film growth. In Task 1 the transport of Na through Mo was evaluated using x-ray photoelectron surface spectroscopy along with diffusion modeling to obtain diffusion coefficients at several temperatures. It was determined that Na diffusion in Mo only occurs along GBs and that oxygen provides an additional driving force to enhance Na transport. Device data revealed that older Mo substrates with a greater amount of surface oxide resulted in slightly higher efficiencies due to enhanced Na incorporation caused by the oxide. This finding shows that Mo substrates could potentially undergo an oxidation treatment prior to CIGS deposition to further improve and control the incorporation of Na. To determine if in-grain Na affects device performance, in Task 1 Na was selectively removed from GBs using heat/rinse cycles. Due to the low temperature of this treatment, Na at GBS remained mobile while diffusion within the bulk was too slow for Na removal from the grain interiors. Changes in electrical properties were evaluated using conductivity and Seebeck coefficient measurements, with both decreasing as Na was removed to reach values similar to Na-free controls samples. This can be explained by the compensation of donor defects by Na, causing an increase in the free carrier concentration. Devices showed a decrease in open-circuit voltage after Na removal confirming that the beneficial effects of GB Na. The findings of this project will provide guidance for rational optimization of Na incorporation procedures in the manufacturing of CIGS solar cells. While it is known that Na segregates at CIGS GBs, the nature and role of Na diffusion into grain interiors was less clear. In Task 2, single crystal CuInSe 2 was used as a model system to represent the grain interiors of CIGS. Crystals processed by two different methods of different compositions and dislocation densities, were evaluated. Diffusion coefficients were obtained at two temperatures after Na diffusion, giving near identical values, ~2x10 11 cm 2 /s and ~6x10 11 cm2/s at 420°C and 480°C, respectively, for each crystal. Characterization confirmed that dislocation densities were too low to significantly impact the effective diffusion coefficient. The Cu-poor crystal had a higher solubility suggesting that Na diffusion is mediated by Cu-vacancies, but was not accompanied by an expected increase in diffusion coefficient. The activation energy for diffusion was similar to values expected for interstitial diffusion, but the large size of Na + ions should result in a solubility that is much lower than what was experimentally measured. A hybrid interstitial-substitutional mechanism is proposed that combines the fast diffusion of interstitial atoms with the high solubility common for substitutional impurities. Lattice diffusion of Na proceeds fast enough that CIGS grain interiors will have Na concentrations near the solubility limit of 1018 cm -3 when manufactured under standard conditions. Na and K treated epitaxial CIS films showed a significant increase in cathodoluminescence emission intensity, indicating a reduction of non-radiative recombination pathways, which is consistent with improvements in CIGS device performance, though the mechanism is not clear. Pathways forward: Despite the success of this project, there are a number of questions remaining related to further the understanding of the chemistry of Na in CIGS films and devices. These include further elucidation of the mechanisms of Na passivation in CIGS GBs, with identification of which defects are involved and confirmation of the possible effects of in-grain Na on device performance. To complete analysis of the cell structure, conformation of the presence and possible chemistries of Na at the CIGS/CdS junction and/or in the front transparent contacts, and its effects on device performance, is needed.

14 SOLAR ENERGY↗

Concentrating Rare Earth Elements in Coal mine drainage Using Coal Combustion Products through Abandoned Mine Land Reclamation

Rare earth elements (REEs) (including scandium, yttrium and a group of 15 lanthanides) are often considered to be critical components in the productions of renewable energy hardware, electric vehicles, health care and military equipment, and consumer electronic products. The demand of REEs has been projected to be growing at an annual rate of 5-9% in the next 25 years. In 2011, the global demand of total rare earth oxides (REOs) was estimated to be approximately 105,000 tons, which is expected to grow to 210,000 tons by 2025. China overwhelmingly dominates the current worldwide rare earth productions but has strategically restricted its exports, causing significant instability for the global market. In response to the increasing demand for REEs and the supply dominance of China, identifying alternative sources of REEs has become a critical issue for the United States and other countries. Coal, coal ash, and coal mine drainage (CMD) are considered to be the alternative sources of REEs. In the U.S., high REE concentrations have been reported to be closely associated with coal deposits, including the Appalachian Basins. When surface and/or groundwater come in contact with geologic strata containing sulfide minerals exposed by coal mining, the accelerated oxidation of sulfide minerals in the presence of ferric iron and/or oxygen can produce sulfuric acid. The process promotes the weathering of REE-bearing rocks and minerals in the host geologic strata. Compared to average river water and seawater, the concentrations of REEs can be orders of magnitude higher in CMD. In this study, we demonstrated a trap-extract-precipitate (TEP) process that can effectively recover REEs from CMD. The three-stage TEP process uses alkaline industrial by-products to capture REEs from CMD and then applies an extraction/precipitation procedure to produce a feedstock that can be economically processed to produce marketable rare earth oxides. The alkaline industrial by-products tested in this study include the residual from a water softening process (DRWP sludge) and two types of stabilized flue gas desulfurization materials (sFGDs). sFGD material is a mixture of lime (CaO) and two coal combustion by-products, calcium sulfite FGD by-product and fly ash. The objectives of this study are to (1) validate the effectiveness and feasibility; (2) determine mechanisms controlling the rare earth recovery, (3) quantify the associated economic and environmental benefits, and (4) evaluate the full-scale application. To achieve these objectives, tasks to be carried out in this proposed project are organized into three phases. In the first phase, the research team collaborated with Ohio Department of Natural Resources, American Electric Power, The Wilds (a nonprofit wildlife conservation organization), and a private landowner to carry out field investigations aimed to screen and evaluate the seasonal changes of rare earths in the CMD discharges that have high recovery potentials. Next, the recovery of REEs from CMD was tested using a series of lab-scale column and batch tests under, respectively, percolation and completely mixed conditions. Results obtained from these lab-scale studies show that all three tested solids are very effective in retaining REEs. Over 98% of the CMD REEs that contacted the solids were captured before the solids exhausted their neutralization capacities. We also determined an extraction process using a non-acid, organic ligand extraction solution that can effectively remobilize the retained REEs from the spent solids (over 90%). The REE concentrate (>7.5 wt. % of total REEs) is then formed in an aeration process. The TEP process uses environmentally benign industrial by-products and a naturally-occurring organic ligand to mitigate CMD and recover REEs. Techno-economic analysis (TEA) and life-cycle assessment (LCA) was carried out in the third phase. The engineering-economic costs and net energy, net CO2 emissions, and water and other requirements were investigated to understand the economic and environmental implications of this process. This work uses mass and energy balances from laboratory-scale experiments to estimate the economic costs and environmental impacts. The results suggest that passive treatment systems that use DRWP sludge are preferred over those that use sFGD material, because of lower economic costs ($89,300/yr with a unit cost of $86/gT-REE vs. $89,800/yr, or $278/gT-REE) and improved environmental performance across all indicators from two different impact assessment methods. These differences are largely attributable to the larger capacity of DRWP sludge in the passive treatment application. We envision this TEP process can be integrated with abandoned mine land (AML) reclamation to create an approach that can add economic incentives for AML reclamation, remediate CMD discharge, and eliminate public safety hazards and threats to local environment and ecological systems posed by AMLs. It can restore lands and communities that are adversely impacted by legacy mining.

01 COAL, LIGNITE, AND PEAT↗

Statistical model based iterative reconstruction (MBIR) in clinical CT systems. Part II. Experimental assessment of spatial resolution performance

Purpose: Statistical model based iterative reconstruction (MBIR) methods have been introduced to clinical CT systems and are being used in some clinical diagnostic applications. The purpose of this paper is to experimentally assess the unique spatial resolution characteristics of this nonlinear reconstruction method and identify its potential impact on the detectabilities and the associated radiation dose levels for specific imaging tasks. Methods: The thoracic section of a pediatric phantom was repeatedly scanned 50 or 100 times using a 64-slice clinical CT scanner at four different dose levels [CTDI{sub vol} =4, 8, 12, 16 (mGy)]. Both filtered backprojection (FBP) and MBIR (Veo{sup ®}, GE Healthcare, Waukesha, WI) were used for image reconstruction and results were compared with one another. Eight test objects in the phantom with contrast levels ranging from 13 to 1710 HU were used to assess spatial resolution. The axial spatial resolution was quantified with the point spread function (PSF), while the z resolution was quantified with the slice sensitivity profile. Both were measured locally on the test objects and in the image domain. The dependence of spatial resolution on contrast and dose levels was studied. The study also features a systematic investigation of the potential trade-off between spatial resolution and locally defined noise and their joint impact on the overall image quality, which was quantified by the image domain-based channelized Hotelling observer (CHO) detectability index d′. Results: (1) The axial spatial resolution of MBIR depends on both radiation dose level and image contrast level, whereas it is supposedly independent of these two factors in FBP. The axial spatial resolution of MBIR always improved with an increasing radiation dose level and/or contrast level. (2) The axial spatial resolution of MBIR became equivalent to that of FBP at some transitional contrast level, above which MBIR demonstrated superior spatial resolution than FBP (and vice versa); the value of this transitional contrast highly depended on the dose level. (3) The PSFs of MBIR could be approximated as Gaussian functions with reasonably good accuracy. (4) Thez resolution of MBIR showed similar contrast and dose dependence. (5) Noise standard deviation assessed on the edges of objects demonstrated a trade-off with spatial resolution in MBIR. (5) When both spatial resolution and image noise were considered using the CHO analysis, MBIR led to significant improvement in the overall CT image quality for both high and low contrast detection tasks at both standard and low dose levels. Conclusions: Due to the intrinsic nonlinearity of the MBIR method, many well-known CT spatial resolution and noise properties have been modified. In particular, dose dependence and contrast dependence have been introduced to the spatial resolution of CT images by MBIR. The method has also introduced some novel noise-resolution trade-off not seen in traditional CT images. While the benefits of MBIR regarding the overall image quality, as demonstrated in this work, are significant, the optimal use of this method in clinical practice demands a thorough understanding of its unique physical characteristics.

60 APPLIED LIFE SCIENCES↗

Diagnostics: Chapter 8 of the special issue: on the path to tokamak burning plasma operation

This chapter presents the activity conducted by the ITPA topical group (TG) on Diagnostics over about the last 15 years. Following a general introduction of the ITER Diagnostics led by their measurement roles, the document is organized in several subchapters detailing the design support, research and development activity conducted by each of the specialist working groups (WGs) of the TG. Please note that the magnetic diagnostics were supported at the TG without a specific WG. Their status is included in the general introduction. In the following some highlights of the subchapter’s contents are provided. Recent advances in ITER first wall (FW) diagnostics for the measurements of plasma-metallic wall interaction in support of the ITER research plan are reported. An InfraRed imaging Video Bolometer for ITER has been developed and tested on several tokamaks to measure the radiated power loss. A laser-induced breakdown spectroscopy (LIBS) technique which utilizes a pulsed laser beam to ablate locally by forming a crater, will measure local tritium inventory in the FW material. Real-time Residual Gas Analyzers will measure the neutral gas composition in a divertor port and an equatorial port during plasma operation. Due to the full metallic FW environment, the plasma-wall interaction in ITER will face several challenges such as the compromised radiated power and divertor heat flux measurements by reflection. Ray tracing and analysis codes have been developed to eliminate and correct the effects of reflection in the measurements. The characteristics of the reflecting surfaces depending on the roughness and angle of the incidence have been measured by dedicated experiments, and the results were applied to the reflection elimination. For the measurement of the metallic impurity radiation induced by eroded metallic atoms, a vacuum ultraviolet spectrometer has been developed and tested. An extensive thermonuclear diagnostic suite will be required to support the operation of ITER and the planned experimental program for future burning plasma experiments. Due to the harsh environmental conditions, the implementation of diagnostic systems in ITER is a major challenge. These conditions include high levels of neutron and gamma fluxes, neutron heating, particle bombardment. Therefore, the selection and design of diagnostic systems must take into account a number of phenomena previously unseen in diagnostic design. For this reason, the measurement of neutrons and confined or lost fast ions, with particular emphasis on alpha particles, is critical to ITER. The diagnostics associated with these measurements will be important for future plasma-burning experiments at ITER. The high neutron emission and very large plasma size in ITER make neutron diagnostics the main diagnostic method used to measure plasma parameters such as fusion power, fusion power density, ion temperature, energy of fast ions and their spatial distributions in the plasma core. Active spectroscopy techniques are methods where a neutral particle beam is injected into the plasma and information on plasma parameters is extracted from the measurement of line emission resulting from the beam-plasma interaction, either by plasma ions or by beam atoms. Spatial localization is achieved by crossing the beamline and multiple observation lines. The ITER plasma will be a high temperature, moderately dense, fully ionized collisional plasma. The plasma facing surfaces are principally metallic being fashioned from beryllium or tungsten but many other elements, arising from either structural or from operational needs, may enter this plasma. The energy range of the emitted photons range from meV (infra-red) to multi keV (x-rays) and originate from all areas of the plasma volume. The primary role of passive emission diagnostics is to identify what is in the plasma from spectral signatures. Extracting quantitative information from these measurements such as impurity content, ion temperature, rotation, degree of detachment and radiated power depends on calibrated instruments, a physics model of the atomic and molecular processes and plasma transport and an analysis workflow that takes into account environmental effects such as reflections. The particular needs for ITER have prompted a multi-machine, many-year effort to address all these aspects and this chapter reviews the work on diagnostic design, experiments and new analysis techniques. An overview of the laser diagnostics to be implemented on ITER is also provided in this paper. This includes descriptions of the Thomson scattering in the core, edge and divertor regions, polarimetry and interferometry diagnostics used for measuring plasma density and also measurements of helium density in the divertor using Laser Induced Flourescence. Techniques which can allow improvements on current measurements are also addressed in particular expanding poloidal polarimetry measurements to measure field fluctuations and proposed use of dispersion interferometery which has a number of advantages over existing methods. This paper identifies particular areas where further research and testing on existing tokamaks is useful even at this advanced stage to inform the design of diagnostics for ITER. Outstanding areas of concern for the implementation of laser diagnostics, in particular with a view to reliable operation are identified. An overview of the latest developments of microwave diagnostic systems and techniques is given. The primary focus is the contributions for ITER—the next step burning plasma experiment—which is supplemented by describing recent progress of techniques applicable for fusion experiments beyond ITER. The contributions are intentionally kept concise, and are being supplemented by a rich list of references for further studies. Radiation induced effects are receiving continuous and well-deserved attention of the ITER diagnostic community and they are in many cases one of the primary design drivers of the ITER diagnostic systems. The paper summarizes recent progress in this area focusing primarily on the ITER diagnostics but in some cases provides also outlook for the possible solutions for even more demanding radiation environment of fusion reactors beyond ITER. Despite advancements in the area of modeling and simulation of various radiation induced effects, experimental testing in a nuclear environment as close as possible to the target one is still seen as unavoidable for proper qualification of particular diagnostic functional elements. Recent advancement within three diagnostic areas: optical diagnostics, magnetics and bolometers is covered. Encouraging results on qualification of silica glass vacuum window assemblies are presented. In the area of magnetic sensors, progress of irradiation tests performed on ITER in-vessel LTCC inductive sensors is presented with outlook for novel technological approaches to inductive sensors utilizing thick printing and photolithography technologies being highlighted. Summary of advancements in the area of steady state magnetic field sensors based on Hall effect is given. New results of neutron irradiation test of the ITER borosilicate glass inserts for vacuum electrical feedthroughs are summarized finding negligible swelling at target level of neutron fluence. Off-line irradiation tests of fiber optic current sensors for plasma current measurement demonstrated that both for gamma doses up to 5 MGy and a total neutron fluence up to 10 15 cm −2 , radiation induced changes are still compatible with required measurement accuracy on ITER. The ITER bolometers are given as an example how considering radiation effects may influence the diagnostic design. Finally, outlook for future main R&D directions is outlined. All optical and laser-based diagnostics in ITER will be using mirrors to guide plasma radiation toward detectors, cameras and sensors. In the hostile plasma, radiation and particle environment the optical characteristics of diagnostic mirrors will degrade directly affecting the entire performance of involved diagnostic systems. An assessment of factors affecting mirror performance is provided. Among the prime adverse factors are deposition of plasma impurities, sputtering of mirror surface and steam ingress in the vicinity of mirrors. Within the International Tokamak Physics Activity with active support by ITER central team and domestic agencies, the structured research and development (R&D) program on mitigation of risks for diagnostic mirrors is underway. Within this program the mirror material development, the passive mitigation of mirror degradation by using diagnostic ducts and shutters along with an active mirror recovery program comprising the in-situ mirror cleaning and calibration is underway. Recent developments in diagnostic mirror R&D are described in this Chapter along with an example of their implementation of R&D solutions in ITER Infrared Thermography diagnostic. An assessment of still open engineering and physics questions, considerations on mirror risks during an early phase of ITER operation are given along with an overview of diagnostic mirror evolution in the late ITER operation stage toward the demonstration fusion power plant. Several crucial areas of diagnostic R&D outlined in ITER Research Plan are addressed. The basic control groups in a fusion reactor can be broken-down in five categories: (1) plasma position, magnetic configuration, and plasma current control, (2) profile control and confinement optimization, (3) MHD control and suppression, (4) edge dissipation control, radiation and plasma exhaust control and (5) break-down optimization. These categories are coupled via the physics (a control action in one domain will affect the other domains) and via shared actuators (e.g. ECRH for impurity accumulation avoidance, current density distribution control and MHD suppression). Consequently, a supervisory control system should determine the priority of the various control tasks, their couplings, and the interfaces with the safety and interlock system. For the systematic development of the various controllers taking the complexity of the plasma and the control system into account, a model-based approach is required. A short historical overview is given of the developments in systems and control theory and control engineering with special emphasis on those developments that are most relevant for Nuclear Fusion research and operation. An overview is given of the state of the field of fusion plasma control for the control categories. It will be shown how synthetic diagnostics are being developed in ITER and how they are used in diagnostic design and design validation and how they can be in model-based controller synthesis using relatively simple models. In modern control methods, multiple diagnostics are used to constrain relatively simple models. The constrained models provide an estimate for the state. This opens the route to state controllers, such as model predictive control. A major challenge in nuclear fusion research is the coherent combination of data from heterogeneous diagnostics and modeling codes for machine control and safety as well as physics studies. Measured data from different diagnostics often provide information about the same subset of physical parameters. Additionally, information provided by some diagnostics might be needed for the analysis of other diagnostics. A joint analysis of complementary and redundant data allows, e.g. to improve the reliability of parameter estimation, to increase the spatial and temporal resolution of profiles, to obtain synergistic effects, to consider diagnostics interdependencies and to find and resolve data inconsistencies. Physics-based modeling and parameter relationships provide additional information improving the treatment of ill-posed inversion problems. A coherent combination of all kind of available information within a probabilistic framework allows for improved data analysis results. The concept of integrated data analysis (IDA) in the framework of Bayesian probability theory is outlined and contrasted with conventional data analysis. Components of the probabilistic approach are summarized and specific ingredients beneficial for data analysis at fusion devices are discussed.

ITER↗

Energy metric prediction for double insertion mutants via the RoseNet deep learning framework

Studying the structural and functional implications of protein mutations is an important task in computational biology and bioinformatics. We leverage our previously proposed RoseNet neural network architecture to predict energy metrics of proteins with double amino acid insertions or deletions (InDels). We train models on previously generated benchmark datasets containing the exhaustive double InDel mutations for three proteins, as well as an additional three proteins for which ∼145k random mutants, each with two InDels, have been generated. We expand on our previous work by evaluating three additional proteins and analyzing domain features that impact the prediction capabilities of RoseNet. These features include InDels into secondary structures and the solvent accessible surface area (SASA) scores of the residues. We uncover further evidence to support that RoseNet has a higher proficiency of generalizing to unseen residue combinations than unseen insertion positions. We also observe that RoseNet produces higher-quality predictions when inserting into a β-sheet over an α-helix. Additionally, when the insertions fall in an area of high SASA, RoseNet often displays better performance than inserting into areas of low SASA.

59 BASIC BIOLOGICAL SCIENCES↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

Athena-I CUBIT Journal Files

The Monte Carlo N-Particle (MCNP)1 transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into con structive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is a time-consuming and error prone process as the complexities of geometries increase. The UM capability was originally designed to work with UM models created with the Abaqus software and ASCII input files that it generates. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with version 6.3, MCNP can also process HDF5 mesh input files. External codes must be used to generate Abaqus input files for MCNP UM calculations. CUBIT, the Sandia National Laboratory automated mesh generation toolkit, can generate a UM model formatted as an Abaqus input file. However, the Abaqus input files exported from CUBIT cannot be used for MCNP simulations because it lacks the proper syntax. A Python script was developed to convert an Abaqus file created by CUBIT to an Abaqus file format that MCNP can process. Creating UM models for complex geometries is not an easy task. The process of creating UM models in CUBIT for MCNP simulations is detailed in. CUBIT provides several user interface options including a graphical user interface (GUI) and a command line interface. A GUI provides an easy way to use CUBIT without learning the CUBIT command syntax. When using CUBIT with either interface option, command lines are written into an ASCII file known as a journal file; this journal file can be edited and archived so that it can be played back in CUBIT to automatically generate a UM model. This report describes the CUBIT journal files of the UM models developed for Athena-I. The Athena platform, an energy-tuning assembly, was developed to spectrally shape the National Ignition Facility (NIF) deuterium-tritium fusion neutron source to a thermonuclear (fusion) plus prompt fission neutron spectrum with capability to act as a short-pulse neutron source. MCNP6 was used for the Athena experiment design analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

ATR Spent Fuel Options Study

The Advanced Test Reactor (ATR) is a materials and fuels test nuclear reactor that performs irradiation services for the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Naval Reactors, the National Nuclear Security Administration (NNSA), and other research programs. ATR achieved initial criticality in 1967 and is expected to operate in support of needed missions until the year 2050 or beyond. It is anticipated that ATR will generate approximately 105 spent nuclear fuel (SNF) elements per year through the year 2050. Idaho National Laboratory (INL) currently stores 2,008 ATR SNF elements in dry storage, 976 in wet storage, and expects to have 1,000 elements in wet storage before January 2017. A capability gap exists at INL for long-term (greater than the year 2050) management, in compliance with the Idaho Settlement Agreement (ISA), of ATR SNF until a monitored retrievable geological repository is open. INL has significant wet and dry storage capabilities that are owned by the DOE Office of Environmental Management (EM) and operated and managed by Fluor Idaho, which include the Idaho Nuclear Technology and Engineering Center’s (INTEC’s) CPP-666, CPP-749, and CPP-603. In addition, INL has other capabilities owned by DOE-NE and operated and managed by Battelle Energy Alliance, LLC (BEA), which are located at the Materials and Fuel Complex (MFC). Additional storage capabilities are located on the INL Site at the Naval Reactors Facility (NRF). Current INL SNF management planning, as defined in the Fluor Idaho contract, shows INTEC dry fuel storage, which is currently used for ATR SNF, will be nearly full after transfer of an additional 1,000 ATR SNF from wet storage. DOE-NE tasked BEA with identifying and analyzing options that have the potential to fulfill this capability gap. BEA assembled a team comprised of SNF management experts from Fluor Idaho, Savannah River Site (SRS), INL/BEA, and the MITRE Corp with an objective of developing and analyzing options for fulfilling the capability gap. This management options analysis is not an alternatives analysis as defined by DOE Order 413.3B; rather, it is an evaluation of near-term, mid term and long-term actions needed to fulfill the capability gap. The actions are described in sufficient detail to inform stakeholders and DOE decision makers regarding a potential path forward. The recommended path forward will inform Fiscal Year 2019 budget formulation, support potential National Environmental Policy Act (NEPA) analyses, and may or may not include capital asset projects.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Unexpected topology of the temperature fluctuations in the cosmic microwave background

We study the topology generated by the temperature fluctuations of the cosmic microwave background (CMB) radiation, as quantified by the number of components and holes, formally given by the Betti numbers, in the growing excursion sets. We compare CMB maps observed by the Planck satellite with a thousand simulated maps generated according to the $Λ$CDM paradigm with Gaussian distributed fluctuations. The comparison is multi-scale, being performed on a sequence of degraded maps with mean pixel separation ranging from 0.05 to 7.33°. The survey of the CMB over S 2 is incomplete due to obfuscation effects by bright point sources and other extended foreground objects like our own galaxy. To deal with such situations, where analysis in the presence of “masks” is of importance, we introduce the concept of relative homology. The parametric χ 2 -test shows differences between observations and simulations, yielding p -values at percent to less than permil levels roughly between 2 and 7°, with the difference in the number of components and holes peaking at more than 3 σ sporadically at these scales. The highest observed deviation between the observations and simulations for b 0 and b 1 is approximately between 3 σ and 4 σ at scales of 3–7°. There are reports of mildly unusual behaviour of the Euler characteristic at 3.66° in the literature, computed from independent measurements of the CMB temperature fluctuations by Planck ’s predecessor, the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. The mildly anomalous behaviour of the Euler characteristic is phenomenologically related to the strongly anomalous behaviour of components and holes, or the zeroth and first Betti numbers, respectively. Further, since these topological descriptors show consistent anomalous behaviour over independent measurements of Planck and WMAP, instrumental and systematic errors may be an unlikely source. These are also the scales at which the observed maps exhibit low variance compared to the simulations, and approximately the range of scales at which the power spectrum exhibits a dip with respect to the theoretical model. Non-parametric tests show even stronger differences at almost all scales. Crucially, Gaussian simulations based on power-spectrum matching the characteristics of the observed dipped power spectrum are not able to resolve the anomaly. Understanding the origin of the anomalies in the CMB, whether cosmological in nature or arising due to late-time effects, is an extremely challenging task. Regardless, beyond the trivial possibility that this may still be a manifestation of an extreme Gaussian case, these observations, along with the super-horizon scales involved, may motivate the study of primordial non-Gaussianity. Alternative scenarios worth exploring may be models with non-trivial topology, including topological defect models.

79 ASTRONOMY AND ASTROPHYSICS↗

Mid-Atlantic U.S. Offshore Carbon Storage Resource Assessment Project (Final Technical Report)

The greatest potential for carbon storage in the northeastern United States lies in the offshore geologic formations comprising the outer continental shelf (Monteverde et al., 2011). Offshore storage can be linked to large point-sources of carbon dioxide (CO 2 ) while avoiding many of the logistical difficulties and potential risks encountered when siting onshore projects, especially in densely populated areas of the East Coast. Offshore CO 2 storage assessment and research in the United States is underway to address uncertainty in potential storage resources, particularly in the mid- and north-Atlantic offshore area. Given the current knowledge base and access to publicly available data, the objectives of the Mid-Atlantic U.S. Offshore Carbon Storage Resource Assessment Project (MAOCSRAP) were fourfold: (1) complete a systematic carbon storage resource assessment of the mid-Atlantic offshore coastal region from the Georges Bank Basin (GBB) through the Long Island Platform (LIP) to the southern Baltimore Canyon Trough (BCT); (2) define key input parameters to reduce uncertainty for offshore storage resource and efficiency estimates; (3) perform a preliminary assessment of risk factors, uncertainties, and data gaps; and (4) engage industry and regulatory stakeholders through development of a Road Map to assist future project planning and implementation. The MAOCSRAP project has compiled, inventoried, and assimilated various publicly available data sets to provide a strong technical basis on which future carbon storage studies and applications can be built. The knowledge infrastructure necessary to support the development of full-scale offshore carbon storage include: Defining the geologic characteristics of candidate storage sites (Chapter 2). Using existing seismic data to better define the continuity of the storage zones and seals (Chapter 3). Cataloguing the hydrogeologic properties of mid-Atlantic offshore storage sites (Chapter 4). Calculating prospective CO 2 storage resources using net effective pore volumes and fluid displacement properties specific to offshore lithologies (Chapter 5). Examining risk factors related to offshore storage (Chapter 6). Communicating with industry and other stakeholders about the future prospects for offshore storage (Chapter 7). Ensuring technology transfer to industry and other stakeholders (Chapter 8). Led by Battelle, this project was conducted by public and private entities with expertise in offshore geology and resources for the study region, including state geological surveys of Delaware, Maryland, and Pennsylvania; United States Geological Survey-Woods Hole Coastal and Marine Science Center and Haifa University; Rutgers University; Harvard University; and Lamont-Doherty Earth Observatory (LDEO) at Columbia University. The storage resource assessment was completed for a broad region offshore of the U.S. East Coast, from Massachusetts to Virginia. The results include high-level storage resource estimates for areas not previously characterized and improved storage resource estimates for geographically expansive portions of offshore geologic units. Results of the project provide a foundation for developing CO 2 storage along offshore areas of the mid-Atlantic United States. A combination of information was used to provide confidence in project results. Key findings of the major project tasks are summarized as follows. Subsurface Data Analysis: Legacy seismic, well-log, core, and biostratigraphic data were digitized, reprocessed, and analyzed using modern techniques, augmenting previous characterization efforts. The sequence stratigraphic framework for the Cretaceous and Jurassic strata was developed by correlating well-log depositional sequences to seismic reflectors identified in the BCT, LIP, and GBB. The integration of well-log stratigraphy with seismic stratigraphy and petrophysical analysis helped to identify and map potential storage zones within sequences of the Logan Canyon (LC), Missisauga, and Mohawk formations, as well as the regionally extensive caprock comprised of the Dawson Canyon shale. Approximately 4,000 kilometers (km) of seismic data was reprocessed to create a structural framework that was used to tie together all available well data, constrain sequence stratigraphic interpretations, and make regional assessments of offshore carbon storage resources. Prospective storage resource estimates were 150 to 1136 megatons (Mt) for the combined storage zones. This suggests mid-Atlantic U.S. offshore formations can store decades of CO 2 from industrial sources in the region. Preliminary local reservoir simulations suggest that injection rates of 1 Mt CO 2 /year may be sustained for 30 years in single injection wells in the BCT. Advanced geologic modeling and new data acquisition are needed to address data gaps and advance carbon capture and storage (CCS) in key offshore areas selected for further investigation. Offshore Risk Factors: Offshore geologic risk factors include soft sediment deformation, unit continuity, sedimentological and structural features, seismicity and hydrates. CO 2 storage risks include inadequate seals, migration/leakage, and chemical interactions leading to decreased storage. Sensitive habitats, environmental impacts, disturbance to seafloor, and other risks need to be identified in advance of project activities and integrated into detailed mitigation plans for all project phases. Stakeholder Engagement: Input and participation from government, industry, and environmental groups provided input into the Road Map and address next steps needed for project deployment. Early engagement and ongoing communication, as well as policy framework development, is key to large-scale deployment of CCS. Next Steps: This project represents an important first step by completing a high-level CO 2 storage resource assessment and building the knowledge infrastructure necessary to improve quantitative storage resource estimates. The data sets that have been curated under this project provide an opportunity to conduct research and development (R&D) needed to address data gaps and reduce risk and uncertainty. Offshore characterization and validation strategies that are systematically designed to provide data and infrastructure that can be upscaled to meet commercial requirements could be developed off of this foundation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Multiview Incomplete Knowledge Graph Integration with application to cross-institutional EHR data harmonization

Objective: The growing availability of electronic health records (EHR) data opens opportunities for integrative analysis of multi-institutional EHR to produce generalizable knowledge. A key barrier to such integrative analyses is the lack of semantic interoperability across different institutions due to coding differences. We propose a Multiview Incomplete Knowledge Graph Integration (MIKGI) algorithm to integrate information from multiple sources with partially overlapping EHR concept codes to enable translations between healthcare systems. Methods: The MIKGI algorithm combines knowledge graph information from (i) embeddings trained from the co-occurrence patterns of medical codes within each EHR system and (ii) semantic embeddings of the textual strings of all medical codes obtained from the Self-Aligning Pretrained BERT (SAPBERT) algorithm. Due to the heterogeneity in the coding across healthcare systems, each EHR source provides partial coverage of the available codes. MIKGI synthesizes the incomplete knowledge graphs derived from these multi-source embeddings by minimizing a spherical loss function that combines the pairwise directional similarities of embeddings computed from all available sources. MIKGI outputs harmonized semantic embedding vectors for all EHR codes, which improves the quality of the embeddings and enables direct assessment of both similarity and relatedness between any pair of codes from multiple healthcare systems. Results: With EHR co-occurrence data from Veteran Affairs (VA) healthcare and Mass General Brigham (MGB), MIKGI algorithm produces high quality embeddings for a variety of downstream tasks including detecting known similar or related entity pairs and mapping VA local codes to the relevant EHR codes used at MGB. Based on the cosine similarity of the MIKGI trained embeddings, the AUC was 0.918 for detecting similar entity pairs and 0.809 for detecting related pairs. For cross-institutional medical code mapping, the top 1 and top 5 accuracy were 91.0% and 97.5% when mapping medication codes at VA to RxNorm medication codes at MGB; 59.1% and 75.8% when mapping VA local laboratory codes to LOINC hierarchy. When trained with 500 labels, the lab code mapping attained top 1 and 5 accuracy at 77.7% and 87.9%. MIKGI also attained best performance in selecting VA local lab codes for desired laboratory tests and COVID-19 related features for COVID EHR studies. Compared to existing methods, MIKGI attained the most robust performance with accuracy the highest or near the highest across all tasks. Conclusions: The proposed MIKGI algorithm can effectively integrate incomplete summary data from biomedical text and EHR data to generate harmonized embeddings for EHR codes for knowledge graph modeling and cross-institutional translation of EHR codes.

Zhou, Doudou↗

Photovoltaic Cable Jackets: A Comparison of Representative Products Using Combined-Accelerated Stress Testing [Slides]

Photovoltaic (PV) cables facilitate the distribution of electricity collected from modules to the energy grid. Durable cabling enables continuous operation of PV installations, whereas cables with a lifetime less than the modules must be replaced - reducing electricity generation and adding to the operating expense. This study primarily focusses on the aging of the key cable types using the combined-accelerated stress testing (C-AST) protocol. Representative cables for utility, building, and control/auxiliary applications were examined. Cable jacket materials examined include: polyolefin, polyethylene, polyamide, poly(vinyl chloride), chlorinated polyethylene, thermoplastic elastomer, and ethylene propylene diene monomer rubber. Specimen characterizations applied include: optical microscopy, mechanical profilometry, instrumented indentation, scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDS), Fourier-transform infrared spectroscopy (FTIR), thermogravimetric analysis (TGA), and differential scanning calorimetry (DSC). A variety of performance and durability characteristics were observed, depending on the base material, polymer formulation, and jacket color. The results from C-AST are analyzed and discussed relative to a recent industry survey on electronic balance of system components in addition to a recent study where similar cables were aged using steady state ultraviolet weathering (International Electrotechnical Commission Technical Specification 62788-7-2). Recommendations are made for the screening, industry qualification, and service life prediction of PV cable jackets.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Deep data analysis via physically constrained linear unmixing: universal framework, domain examples, and a community-wide platform

Abstract Many spectral responses in materials science, physics, and chemistry experiments can be characterized as resulting from the superposition of a number of more basic individual spectra. In this context, unmixing is defined as the problem of determining the individual spectra, given measurements of multiple spectra that are spatially resolved across samples, as well as the determination of the corresponding abundance maps indicating the local weighting of each individual spectrum. Matrix factorization is a popular linear unmixing technique that considers that the mixture model between the individual spectra and the spatial maps is linear. Here, we present a tutorial paper targeted at domain scientists to introduce linear unmixing techniques, to facilitate greater understanding of spectroscopic imaging data. We detail a matrix factorization framework that can incorporate different domain information through various parameters of the matrix factorization method. We demonstrate many domain-specific examples to explain the expressivity of the matrix factorization framework and show how the appropriate use of domain-specific constraints such as non-negativity and sum-to-one abundance result in physically meaningful spectral decompositions that are more readily interpretable. Our aim is not only to explain the off-the-shelf available tools, but to add additional constraints when ready-made algorithms are unavailable for the task. All examples use the scalable open source implementation from https://github.com/ramkikannan/nmflibrary that can run from small laptops to supercomputers, creating a user-wide platform for rapid dissemination and adoption across scientific disciplines.

47 OTHER INSTRUMENTATION↗

Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (Final Report)

This document presents final technical findings for the project Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (FE0031782). The project is part of the U.S. Department of Energy Oil and Gas Program to develop and advance technologies that can significantly improve the recovery efficiencies of unconventional oil and gas resources. The overall objective of this project is to improve the ultimate recovery from unconventional oil and gas (UOG) resources in the United States by developing a method for using unrefined natural gas liquids (NGLs) as treatment fluids to improve hydrocarbon production. Horizontal extended-lateral drilling coupled with high volume hydraulic fracturing has significantly increased production from UOG resources in the U.S. However, the recovery efficiency is low compared to the estimated oil and gas in place. Recent data indicate that less than 10% of the oil in the liquid-rich UOG reservoirs is produced. Alternative completion methods using NGLs could increase production (Battelle, 2016; Wan et al., 2013; Wan, 2013; Downey et al, 2021); however, field validation tests are needed to develop an approach that is economical, efficient, and compatible in the UOG setting to advance towards commercial deployment. This project aims to develop and field test a method to improve recovery of oil resources in UOG shale plays by using Y-Grade NGLs, or a similar combination of NGLs, as treatment fluids. Refined NGLs have been used as a hydraulic treatment fluid in UOG plays for decades and are shown to be particularly effective because their miscibility with oil allows oil to flow more freely; however, the use of Y-Grade (unrefined) NGLs has not been studied. The use of Y-Grade NGLs would be advantageous over refined NGLs because Y-Grade NGLs do not require infrastructure or investment in refining and are already being produced from many UOG reservoirs. The concept was tested and monitored in the field at a commercial well site owned by project partner Hopco, Ltd. The project team, which consists of multiple oil and gas operators, Linde Gas North America LLC (Linde) and the Ohio Division of Geological Survey (ODGS), has extensive experience with oil and gas production in the Appalachian basin and the ability to work together quickly to solve technical issues and research needs. A key part of the proposed work was the use of existing wells for field testing and monitoring. A total of four wells (three vertical and one horizontal) were available for this project. One of the vertical wells was utilized as the test well for the NGL treatment test. A nearby vertical well was used for microseismic monitoring. The remaining vertical well and the horizontal well provided a baseline for typical UOG production in the oil window of the Utica/Point Pleasant (UPP). Major technical tasks of the project include characterization of the geotechnical properties of the UPP with an emphasis on the field site; design and planning for the NGLs testing; field testing and monitoring; analysis and integration of field data; and economic and resource/reserve assessment. The Shoman monitoring well and Doughty NGL treatment well were successfully plugged back during September-October 2020 in preparation for treatment and monitoring. A nitrogen diagnostic fracture injection test (DFIT) was completed on the Doughty well on July 22, 2021, consisting of 133,000 scf (91 Bbl.) of nitrogen. A nitrogen foam frac was completed in the Utica-Point Pleasant interval on August 17, 2021 with funding from outside sources. A microseismic monitoring array was installed in the Shoman well and monitored microseismic activity during the Doughty well frac job. Y-Grade NGL injection commenced on August 26, 2021. A total of 215 Bbl. was injected but the job was shut down due to a small leak on the suction hose on the pump truck. The Y-Grade treatment resumed on August 27, 2021 and an additional 726 Bbl. of Y-Grade NGL was injected at a well head pressure of 3850 psi. Total volume of injected Y-Grade over the two days of injection was 941 Bbl. The well was shut-in for 17 days following injection to allow the Y-Grade NGLs to soak on the formation. Y-Grade treatment flow back commenced on 9/13/2021 on a weekly basis. Production data, including surface pressures, oil, nitrogen, natural gas, and flow times was measured and recorded. Periodic gas samples were collected and analyzed to determine composition of flowback gas. As of July 2022, the treatment well had produced 726 Bbl of oil and 2,888 mcf gas. In August 2022, tubing and packer in the well were removed and a pump was installed to enhance oil recovery. Currently, the operator is producing the well about 2 days a week for a few hours. The performance of the NGL treatment test was evaluated based on reservoir simulations of the treatment process, processing of well testing data, analysis of micro seismic monitoring data, and production data analysis. This analysis suggested that oil production in the small test would continue through 2025. Upscaling the treatment to a horizontal Utica Point-Pleasant well would allow more oil recovery, but the process would involve more investment, services, and operational support. An economic analysis was conducted for scenarios aimed at upscaling the NGL treatment process for more typical horizontal UPP wells in the Appalachian Basin.

02 PETROLEUM↗

PRIMA-X - Performance Retargeting of Instrumentation, Measurement, and Analysis Technologies for Exascale Computing

The PRIMA-X (Performance Retargeting of Instrumentation, Measurement, and Analysis Technologies for Exascale Computing) project is the successor of the DOE PRIMA (Performance Refactoring of Instrumentation, Measurement, and Analysis Technologies for Petascale Computing) project, which addressed the challenge of creating a core measurement infrastructure that would serve as a common platform for both integrating leading parallel performance systems (notably TAU and Scalasca) and developing next-generation scalable performance tools. The PRIMA-X project shifts the focus away from refactorization of robust performance tools towards a re-targeting of the parallel performance measurement and analysis architecture for extreme scales. The massive concurrency, asynchronous execution dynamics, hardware heterogeneity, and multi-objective prerequisites (performance, power, resilience) that identify exascale systems introduce fundamental constraints on the ability to carry forward existing performance methodologies. In particular, there must be a deemphasis of per-thread observation techniques to significantly reduce the otherwise unsustainable flood of redundant performance data. Instead, it will be necessary to assimilate multi-level resource observations into macroscopic performance views, from which resilient performance metrics can be attributed to the computational features of the application. This requires a scalable framework for node-level and system-wide monitoring and runtime analyses of dynamic performance information. Also, the interest in optimizing parallelism parameters with respect to performance and energy drives the integration of tool capabilities in the exascale environment further. Initially, PRIMA-X was a collaborative project between the University of Oregon (lead institution) and the German Research School for Simulation Sciences (GRS). Because Prof. Wolf, the PI at GRS, accepted a position as full professor at Technische Universitaet Darmstadt (TU Darmstadt) starting February 1st, 2015, the work of GRS was continued at TU Darmstadt. This report covers the work at TU Darmstadt after the transition. The first main accomplishment of TU Darmstadt is the development of different techniques to aggregate performance data on the level of threads. TU Darmstadt evaluated different schemes designed at GRS during the first phase of this project and integrated them into Score-P, a widely used performance measurement framework for HPC applications. The second main accomplishment is a substantial increase of Score-P’s scalability, achieved by improving its internal representation of the underlying system. Third, we extended our performance-modeling tool Extra-P to automatically create models with more than one input parameter (e.g., process count and input size), which was previously impossible. These models allow the performance behavior of an application to be extrapolated to larger machines. Specific applications of the new multi-parameter capability realized in this project include modeling the isoefficiency function of task-based applications and a novel co-design methodology. Further extensions of Extra-P comprise online modeling functionality, automatic refinement of the model search space, and segmented performance modeling. Further contributions include new locality metrics for both shared- and distributed-memory applications, a new way of measuring the impact of inter-application inference on the execution time of an application, and finally an automatic method that selects appropriate OpenMP constructs for the parallelization of sequential code and classifies the affected variables based on their sharing semantics.

97 MATHEMATICS AND COMPUTING↗