SEARCH · Engineering Papers

Results for “well data”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2022

$K^{*}(892)^0$ meson production in inelastic p+p interactions at 40 and 80 $\text{ GeV }\!/\!c$ beam momenta measured by NA61/SHINE at the CERN SPS

Measurements of $K^{*}(892)^0$ resonance production via its $K^{+}\pi ^{-}$ decay mode in inelastic p+p collisions at beam momenta 40 and 80 $\text{ GeV }\!/\!c$ ($\sqrt{s_{NN}}=8.8$ and 12.3 $\text{ GeV }$) are presented. The data were recorded by the NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The template method was used to extract the $K^{*}(892)^0$ signal. Transverse momentum and rapidity spectra were obtained. The mean multiplicities of $K^{*}(892)^0$ mesons were found to be $(35.1 \pm 1.3 \mathrm {(stat)} \pm 3.6 \mathrm {(sys))} \cdot 10^{-3}$ at 40 $\text{ GeV }\!/\!c$ and $(58.3 \pm 1.9 \mathrm {(stat)} \pm 4.9 \mathrm {(sys))} \cdot 10^{-3}$ at 80 $\text{ GeV }\!/\!c$. The NA61/SHINE results are compared with the Epos1.99 and Hadron Resonance Gas models as well as with world data. The transverse mass spectra of $K^{*}(892)^0$ mesons and other particles previously reported by NA61/SHINE were fitted within the Blast-Wave model. The transverse flow velocities are close to 0.1–0.2 of the speed of light and are significantly smaller than the ones determined in heavy nucleus-nucleus interactions at the same beam momenta.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Storage Field Development Plan: One Earth Energy

This Storage Field Development plan presents the Storage Complex characterization results, construction, monitoring, and operational plans, and costs associated with the proposed One Earth Sequestration Carbon Capture and Storage (OES-CCS) site in McLean County, Illinois, near Gibson City. The proposed storage complex, known as the Mt. Simon Storage Complex, comprises the Cambrian Mt. Simon Sandstone reservoir and the primary seal, the Cambrian Eau Claire Formation. The lowermost Underground Source of Drinking Water (USDW) identified for the site is the Ordovician St. Peter Sandstone. Geologic characterization of the Mt. Simon Storage Complex at the OES-CCS site was performed by the Illinois Storage Corridor CarbonSAFE Phase III project, which also prepared and submitted three UIC Class VI applications to construct three injection wells; the permit applications were submitted and are in the federal EPA review process. A characterization well, OEE #1, was drilled to collect site-specific data. These data were analyzed and used to develop the UIC Class VI applications. The OEE #1 well will be converted to an in-zone monitoring (IZM) well for the injection phase. The proposed buildout for the OES-CCS site includes (1) three injection wells (OES #1, OES #2, and OES #3), (2) two IZM wells, (3) two above confining zone (ACZ) monitoring wells, one of which will be used to monitor the lowermost USDW, (4) capture and compression facilities, and (5) transportation facilities, i.e., pipelines. A pre-operational testing program was proposed in the Class VI permit application and will be employed at the site pending approval. Additional pre-injection (baseline), syn-injection, and post-injection monitoring and site care procedures will be followed by OES to ensure that injection activities are protective of human health and the environment. Injection is scheduled to begin in 2025, distributed across the three injection wells in accordance with the permit operating conditions. One Earth Sequestration intends to inject up to 90 million tonnes of CO 2 over a period of approximately 20 years. Injection will begin at approximately 0.5 million tonnes of CO 2 annually and ramp up to a maximum of 4.5 million tonnes annually. Daily injection rates are expected to range from 1,400 to 1,500 tonnes per day initially and reach a maximum of approximately 4,225 tonnes per day, depending on site geology and injectivity at each injection well location, and CO 2 availability. The costs associated with the OES-CCS project include pre-operational costs (e. g. additional seismic data acquisition and well drilling), capture and transportation facility and equipment costs, predicted field operating expenditures (OpEx), and decommissioning and post-injection site care (PISC) costs. The risks associated with project activities, such as site construction, injection operations, and verification of secure storage were evaluated, and mitigation strategies proposed to alleviate those risks.

09 BIOMASS FUELS↗

DOE OSTI2024

Detecting fractures and monitoring hydraulic fracturing processes at the first EGS Collab testbed using borehole DAS ambient noise

Enhanced geothermal systems (EGS) require cost-effective monitoring of fracture networks. We validate the capability of using borehole distributed acoustic sensing (DAS) ambient noise for fracture monitoring using core photos and core logs. The EGS Collab project has conducted 10 m scale field experiments of hydraulic fracture stimulation using 50–60 m deep experimental wells at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. The first EGS Collab testbed is located at 1616.67 m (4850 ft) depth at SURF and consists of one injection well, one production well, and six monitoring wells. All wells are drilled subhorizontally from an access tunnel called a drift. The project uses a single continuous fiber-optic cable installed sequentially in the six monitoring wells to record DAS data for monitoring hydraulic fracturing during stimulation. We analyze 60 s time records of the borehole DAS ambient noise data and compute the noise root-mean-square (rms) amplitude on each channel (points along the fiber cable) to obtain DAS ambient noise rms amplitude depth profiles along the monitoring wellbore. Our noise rms amplitude profiles indicate amplitude peaks at distinct depths. We compare the DAS noise rms amplitude profiles with borehole core photos and core logs and find that the DAS noise rms amplitude peaks correspond to the locations of fractures or lithologic changes indicated in the core photos or core logs. We then compute the hourly DAS noise rms amplitude profiles in two monitoring wells during three stimulation cycles in 72 h and find that the DAS noise rms amplitude profiles vary with time, indicating the fracture opening/growth or closing during the hydraulic stimulation. Our results demonstrate that borehole DAS passive ambient noise can be used to detect fractures and monitor fracturing processes in EGS reservoirs.

58 GEOSCIENCES↗

DOE OSTI2020

Historic and modern approaches for discovery of abandoned wells for methane emissions mitigation in Oil Creek State Park, Pennsylvania

Hundreds of oil wells were drilled along Oil Creek in Pennsylvania in the mid-16 1800s, birthing the modern oil industry. No longer in operation, many wells are now classified as abandoned, and, due to their age, their locations are either unknown or inaccurately recorded. These historic well sites present environmental, safety, and economic concerns in the form of possible methane leaks and physical hazards. Airborne magnetic and LiDAR surveys were conducted in the Pioneer Run watershed in Oil Creek State Park to find abandoned wells in a historically significant but physically challenging location. Wells were drilled in this area prior to modern geolocation and legal documentation. Although a large number of old wells were abandoned summarily without remediation of the site, much of the land area within Oil Creek State Park is now covered in trees and dense underbrush, which can obscure wellheads. The thick vegetation and steep terrain limited the possibility of ground-based surveys to easily find well sites for methane emissions studies. The data from remote sensing surveys were used to corroborate potential well locations from historic maps and photographs. Potential well sites were verified in a ground-based field survey and monitored for methane emissions. Two historic photographs documenting oil activity in the late 1800s were georeferenced using a combination of magnetic and LiDAR data. LiDAR data, which were more useful in georeferencing and in field verification, identified 290 field locations in the Pioneer Run watershed, 86% of which were possible well sites. Sixty-two percent of the ground-verified wells remained unplugged and comprised the majority of leaking wells. The mean methane emissions factor for unplugged wells was 0.027 ± 0.099 kg/day, lower than other Appalachian Basin methane emissions estimates. LiDAR was used for the first time, in combination with an airborne magnetic survey, to reveal underground oil industry features and inform well identification and remediation efforts in difficult-to-navigate regions. In the oldest oil fields, where well casing has been removed or wood conductor casing was installed, historic photographs provide additional lines of evidence for oil wells where ground disturbances have concealed surface features. Identification of well sites is necessary for mitigation efforts, as unplugged wells emit methane, a potent greenhouse gas.

02 PETROLEUM↗

DOE OSTI2021

Machine learning enables national assessment of wind plant controls with implications for land use

Summary As deployment of wind energy continues to expand, computationally efficient tools for predicting wind plant performance over a wide range of layout designs, technology innovations, and spatial locations are increasingly important for policy and investment decisions. We demonstrate two approaches to training a surrogate model to predict annual energy production (AEP) of parameterized wind plant layouts: one using a Gaussian process (GP) and the other using a fully convolutional neural network (FCNN). We leverage the powerful FCNN architecture by encoding wind plant design parameters and output response surface as an image. The FCNN produces more accurate results than the GP with mean absolute errors equivalent to 1% and 1.9% of plant rated power, respectively, although the GP performs well under limited training data and provides useful uncertainty information. We also evaluate a surrogate model for wake steering, enabling a nationwide assessment of the impact of plant control strategies and plant layout decisions. Across two million locations, we find that wake steering strategies boost AEP with relative gains upwards of 3%. Gains are most pronounced at sites without a dominant wind direction and where layout optimization is less fruitful. Additionally, we perform a nationwide sensitivity analysis showing that wake steering can mitigate wake losses from higher density plant layouts. Our results suggest that regions which have not been previously viable for wind deployment due to moderate wind resources are especially well enhanced by wake steering strategies that could help overcome land constraints and inflexible layout options, potentially identifying new deployment opportunities.

17 WIND ENERGY↗

DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Structure of Offshore Low-Level Jet Turbulence and Implications to Mesoscale-to-Microscale Coupling

This paper explores realistic nonstationary atmospheric boundary layer (ABL) turbulence arising from nonstationarity at the mesoscale, particularly within offshore low-level jets with implications to offshore wind farms, using high-fidelity multiscale large-eddy simulations (LES). To this end, we analyzed the single-point turbulence statistical structure of a North-Atlantic offshore LLJ event simulated using high-resolution LES (AMR-Wind). The nonstationary LLJ is simulated using a mesoscale-to-microscale coupled (MMC) simulation procedure involving data assimilation of mesoscale velocity and temperature data from the Weather Research and Forecasting (WRF) model. Unlike the assimilation of mesoscale velocity data into the LES, the direct assimilation of temperature profiles had a strong impact on turbulence stratification, thereby causing erroneous predictions of turbulence both above and within the jet layer. Various approaches to mitigate this effect have resulted in multiple (four) variants of this MMC strategy. Outcomes from this work clearly show that the turbulence within the low-level jet is a strong function of the MMC approach as the turbulence structure within the low-level jet is dependent on the flux of residual turbulence from outside the jet, which in turn depends on the temperature forcing history. Additionally, the turbulence predicted by all these different methods (as well as the observation data) show similar deviations from equilibrium as evidenced by comparisons with idealized atmospheric turbulence structure obtained using the same numerical method. In general, we observe that the predicted LLJ turbulence tends to differ from canonical ABL turbulence with comparable shear. Particularly, the combination of shear and turbulence observed in such nonstationary low-level turbulence cannot be matched using equilibrium settings and therefore, represents a critical use-case for both testing and leveraging meso–micro coupling strategies.

17 WIND ENERGY↗

DOE OSTI2023

A Novel Technique of Extracting UCN Decay Lifetime from Storage Chamber Measurements Dominated by Scattering Losses

The neutron’s lifetime is a critical parameter in the standard model. Its measurements, particularly measurements using both beamline and ultracold neutron storage techniques, have revealed significant tension. In this work, we review the status of the tension between various measurements, especially in light of the insights provided by the β-decay correlation measurements. We revisit the lifetime measurement in a material storage chamber, dominated by losses from scattering off the walls of the storage chamber. The neutron energy spectra and associated uncertainties were, for the first time, well-characterized using storage data alone. Such models have applications in the extraction of the mean time between wall bounces, which is a key parameter for neutron storage disappearance experiments in search of neutron oscillation. A comparison between the loss model and the number of neutrons stored in a single chamber allowed us to extract a neutron lifetime of τn*=880(+158/−78)stat.(+230/−114)sys.s (68.3% C.I.). Though the uncertainty of this lifetime is not competitive with currently available measurements, the highlight of this work is that we precisely identified the systematic sources of uncertainty that contribute to the neutron lifetime measurements in material storage bottles, namely from the uncertainty in the energy spectra, as well as from the storage chamber surface parameters of the Fermi potential and loss per bounce. In doing so, we highlight the underestimation of the uncertainties in the previous Monte Carlo simulations of experiments using the technique of ultracold neutron storage in material bottles.

Mohanmurthy, Prajwal (ORCID:0000000275737010)↗

DOE OSTI2025

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗

DOE OSTI2024

Forced Oscillation Source Location in Power Systems Using MVMD-assisted DEF in TF Plane

Fast and accurate forced oscillation source location (FOSL) is critical for mitigating forced oscillations in power systems. Here, to overcome the shortcomings of traditional dissipating energy flow (DEF) in decomposing the forced oscillation components from the measurement responses, a multivariate variational mode decomposition (MVMD)-assisted FOSL is proposed in this paper to locate the forced oscillation sources from measurements in time-frequency (TF) plane. The multi-channel measurement matrix of each generator is formed, then the multi-channel modes are decomposed from the formed multi-channel measurement matrix using the MVMD. Further, the intrinsic mode functions (IMFs) associated with the forced oscillation mode are distinguished from the decomposed modes according to the relative energy weights. With the distinguished IMFs, the MVMD-assisted DEF calculation method is developed to locate the forced oscillation source. Finally, the performance of the proposed method is evaluated using the simulation data of the WECC 179-bus test system and IEEE-NASPI Oscillation Source Location Contest as well as the actual PMU data of ISO New England. The results validated the effectiveness and accuracy of the proposed method in the FOSL.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Electroanalytical characterization of Np( VI )/Np( V ) redox in a pentadentate ligand environment and stabilization of [Np V O 2 ] + by hydrogen bonding

The redox chemistry of the actinyl cations (AnO 2 n+ ) heavily influences their reactivity and speciation in solution, but the redox properties of the actinyls in non-aqueous media have received far less attention than they deserve. Here, the non-aqueous electrochemistry of a chemically reversible Np(VI)/Np(V) redox manifold is reported in both protic (CH 3 OH) and aprotic (CH 3 CN) organic media. Using a neutral Np(VI) complex supported by a chelating and strongly donating pentadentate ligand that was fully characterized in prior work, a clean Np(VI)/Np(V) redox couple was found to be accessible under ambient conditions. Coupled electrochemical and spectroscopic studies, as well as simulations of cyclic voltammetry data, confirm the 1e− nature of this couple and establish it to be chemically reversible and nearly electrochemically reversible as well. Bulk electrolysis of a solution of the neutral Np(VI) complex facilitated isolation of the corresponding anionic and monomeric Np(V) species. Data from X-ray diffraction analysis as well as optical and vibrational spectroscopies provide strong evidence in support of metal-centered reduction and the Np(V) oxidation state, findings that are in accord with the measured reduction potentials. Distinctive hydrogen bonding interactions between the terminal (yl) oxo groups and water molecules appear to stabilize the isolated Np(V) species in the solid state, providing insight into the features that afford the uncommon chemically reversible redox encountered in this system.

Mikeska, Emily R.↗

DOE OSTI2021

A new hypothesis for the origin of Amazonian Dark Earths

Amazonian Dark Earths (ADEs) are unusually fertile soils characterised by elevated concentrations of microscopic charcoal particles, which confer their distinctive colouration. Frequent occurrences of pre-Columbian artefacts at ADE sites led to their ubiquitous classification as Anthrosols (soils of anthropic origin). However, it remains unclear how indigenous peoples created areas of high fertility in one of the most nutrient-impoverished environments on Earth. Here, we report new data from a well-studied ADE site in the Brazilian Amazon, which compel us to reconsider its anthropic origin. The amounts of phosphorus and calcium—two of the least abundant macronutrients in the region—are orders of magnitude higher in ADE profiles than in the surrounding soil. The elevated levels of phosphorus and calcium, which are often interpreted as evidence of human activity at other sites, correlate spatially with trace elements that indicate exogenous mineral sources rather than in situ deposition. Stable isotope ratios of neodymium, strontium, and radiocarbon activity of microcharcoal particles also indicate exogenous inputs from alluvial deposition of carbon and mineral elements to ADE profiles, beginning several thousands of years before the earliest evidence of soil management for plant cultivation in the region. Our data suggest that indigenous peoples harnessed natural processes of landscape formation, which led to the unique properties of ADEs, but were not responsible for their genesis. If corroborated elsewhere, this hypothesis would transform our understanding of human influence in Amazonia, opening new frontiers for the sustainable use of tropical landscapes going forward.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Perturbative and nonperturbative studies of CFTs with MN global symmetry

Fixed points in three dimensions described by conformal field theories with \ensuremath{M N}_{m,n} = O(m)^n\rtimes S_n M N m , n = O ( m ) n ⋊ S n global symmetry have extensive applications in critical phenomena. Associated experimental data for m=n=2 m = n = 2 suggest the existence of two non-trivial fixed points, while the \varepsilon ε expansion predicts only one, resulting in a puzzling state of affairs. A recent numerical conformal bootstrap study has found two kinks for small values of the parameters m m and n n , with critical exponents in good agreement with experimental determinations in the m=n=2 m = n = 2 case. In this paper we investigate the fate of the corresponding fixed points as we vary the parameters m m and n n . We find that one family of kinks approaches a perturbative limit as m m increases, and using large spin perturbation theory we construct a large m m expansion that fits well with the numerical data. This new expansion, akin to the large N N expansion of critical O(N) O ( N ) models, is compatible with the fixed point found in the \varepsilon ε expansion. For the other family of kinks, we find that it persists only for n=2 n = 2 , where for large m m it approaches a non-perturbative limit with \Delta_\phi\approx 0.75 Δ ϕ ≈ 0.75 . We investigate the spectrum in the case \ensuremath{M N}_{100,2} M N 100 , 2 and find consistency with expectations from the lightcone bootstrap.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Interferometric SAR modelling of near surface data to improve geological model in the Surat Basin, Australia

Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.

58 GEOSCIENCES↗

DOE OSTI2020

Recent evolutions in the LTDP CDF project

In the latest years, CNAF (the national center of the Italian Institute for Nuclear Physics INFN dedicated to Research and Development on Information and Communication Technologies) has been working on the Long Term Data Preservation (LTDP) project for the CDF experiment, active at Fermilab from 1990 to 2011. The main aims of the project are to protect the most relevant part of the CDF RUN-2 data collected between 2001 and 2011 and already stored on tape at CNAF (4 PB), as well as to ensure the availability and the access to the analysis facility to those data over time. Lately, the CDF database, hosting information about CDF datasets such as their structure, file locations and metadata, has been imported from Fermilab to CNAF. Also, the Sequential Access via Metadata (SAM) station data handling tool for CDF data management, that allows to manage data transfers and to retrieve information from the CDF database, has been properly installed and configured at CNAF. This was a fundamental step in the perspective of a complete decommissioning of CDF services on Fermilab side. An access system has been designed and tested to submit CDF analysis jobs, using CDF software distributed via CERN Virtual Machine File System (CVMFS) and requesting delivery of CDF files stored on CNAF tapes, as well as data present only on Fermilab storage archive. Moreover, the availability and the correctness of all CDF data stored on CNAF tapes has been verified. This paper describes all these recent evolutions in detail, presenting the future plans for the LTDP project at CNAF.

Fornari, Federico↗

DOE OSTI2023

Assessment of Sodium Thermal Stratification Models Utilizing the TSTF Benchmark

As a result of certain transient scenarios, a thermally stratified layer of liquid sodium can develop in the bulk coolant volumes of a sodium-cooled fast reactor (SFR). In addition to the effects a stratification layer has on the temperature of the heat transport system, a stratification layer can also influence the transition to and establishment of natural circulation flow, which plays an important role in passive cooling and the inherent safety of a pool-type SFR. Therefore, the ability to accurately capture thermal stratification phenomena is important when demonstrating the safety basis of a pool-type SFR during transient sequences. The present work assesses various computational models with different fidelities in their ability to predict thermal stratification in the upper plenum of an SFR. Each computational model will be assessed using the data generated at the Thermal Stratification Test Facility (TSTF) located at the University of Wisconsin-Madison. Using measured flow rate and inlet temperature data, the measured temperature distributions of the tests are compared to the predictions of the lumped volume-based models in SAS4A/SASSYS-1, a 1D-based model in SAM, and a 3-D computational fluid dynamics (CFD) model using STAR-CCM+. The relative performance of the various computational methods is assessed with respect to key metrics such as bulk coolant temperature distribution and plenum exit temperature. A total of eight tests are analyzed, covering different combinations of flow rates (3 and 10 GPM) and upper internal structure (UIS) configurations (none, solid, porous, and open) The perfect mixing model of SAS4A/SASSYS-1 provides the highest accuracy when the flow rate is high and there is no UIS in the test vessel, as high flow rate injection promotes thermal mixing of the sodium in the test vessel. For most of the analyzed tests, the stratified volume model of SAS4A/SASSYS-1 is able to predict the delay in the outlet temperature drop and temperature distribution in the test vessel by a small number of layers to represent thermal stratification. However, the stratified volume model can only simulate a maximum of three temperature layers within a volume and when a layer approaches the elevation of the outlet, the predicted outlet temperature can demonstrate rapid, non-physical changes. The 1-D axial mixing model of SAM provides results that agree reasonably well with the measured data in the prediction of the temporal evolution of the outlet temperature with the exception of the case with a high flow rate and no UIS. The SAM 1-D model has a similar level of accuracy to CFD results when it comes to predicting the outlet temperature. CFD shows overall good agreement in predicting the temperature distribution in the test vessel and outlet temperature. As CFD can model the test vessel geometry in detail, it performs well in the cases of complex geometries such as tests that included a UIS and internal flow through the UIS resulting in active mixing of the coolant in the test vessel. Each of the models discussed in the present work has the potential to be useful during the various stages of reactor design, analysis, and licensing. The lumped-volume approach can be applied for fast turnaround safety calculations to obtain overall reactor behavior during transients. The 1-D models provide improved accuracy when stratification is expected for a relatively low increase in the computational cost. The CFD model can be utilized for confirmatory analysis of the 1-D model, when experimental measurements are not available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Methanol oxidation up to 100 atm in a supercritical pressure jet-stirred reactor

Methanol (CH 3 OH) has attracted considerable attention as a renewable fuel or fuel additive with low greenhouse gas emissions. Methanol oxidation was studied using a recently developed supercritical pressure jet-stirred reactor (SP-JSR) at pressures of 10 and 100 atm, at temperatures from 550 to 950 K, and at equivalence ratios of 0.1, 1.0, and 9.0 in experiments and simulations. Here, the experimental results show that the onset temperature of CH 3 OH oxidation at 100 atm is around 700 K, which is more than 100 K lower than the onset at 10 atm and this trend cannot be predicted by the existing kinetics models. Furthermore, a negative temperature coefficient (NTC) behavior was clearly observed at 100 atm at fuel rich conditions for methanol for the first time. To understand the observed temperature shift in the reactivity and the NTC effect, we updated some key elementary reaction rates of relevance to high pressure CH 3 OH oxidation from the literature and added some new low-temperature reaction pathways such as CH 2 O + HO 2 = HOCH 2 O 2 (RO 2 ), RO 2 + RO 2 = HOCH 2 O (RO) + HOCH 2 O (RO) + O 2 , and CH 3 OH + RO 2 = CH 2 OH + HOCH 2 O 2 H (ROOH). Although the model with these updates improves the prediction somewhat for the experimental data at 100 atm and reproduces well high-temperature ignition delay times and laminar flame speed data in the literature, discrepancies still exist for some aspects of the 100 atm low-temperature oxidation data. In addition, it was found that the pressure-dependent HO 2 chemistry shifts to lower temperature as the pressure increases such that the NTC effect at fuel-lean conditions is suppressed. Therefore, as shown in the experiments, the NTC phenomenon was only observed at the fuel-rich condition where fuel radicals are abundant and the HO 2 chemistry at high pressure is weakened by the lack of oxygen resulting in comparatively little HO 2 formation.

33 ADVANCED PROPULSION SYSTEMS↗