Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model based definitions”

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

Hybrid Solar System (Final Scientific/Technical Report)

GTI Energy (GTI) teamed with the University of California at Merced (UCM) to scaleup the hybrid solar system (HSS) technology for demonstrating its performance at the US Gypsum (USG) plant in Plaster City, California. The technology integrates two-stage concentrating solar collector with matching particle thermal transport and storage (TSS) system to deliver cost-effective, and on-demand distributed high temperature industrial process heat up to 600°C with solar thermal, in this case to a gypsum kettle, to reduce its fuel use and carbon footprint. Current solar technologies, which reach these temperatures, are not distributable (towers) or cost-effective (dish). The research team developed a conceptual system design for host site retrofit, including preliminary heat balance, process flow diagram, particle to process heat exchanger and equipment placements at the site. Subsequently, parallel efforts were carried out at UCM to design, build and test a 12 m long commercial scale prototype concentrating thermal-only collector system and at GTI to design, build and test a matching 650°C capable particle TTS system. The nominal 50 kWth collector consists of a parabolic trough and three 4 m long two-stage receivers in series. Prior to on-sun testing, a 4 m long receiver was fabricated and successfully tested at 650 °C in a laboratory setting for 100 hrs of continuous operation showing less than 15% radiation loss. A 7 m wide x 17 m long parabolic trough was then installed at UCM for on-sun testing of the 12 m long receiver, and concurrently several 4 m long receivers were built. The optics of the parabolic trough were calibrated, and on-sun test were carried out on 12 m long receivers. During tests, the intense solar radiation (53x) caused the absorber tubes in the receivers to bend, reducing the overall optical efficiency. To address the bending issue, a self-consistent algorithm that includes ray tracing, thermal and deformation models was developed to perform thermal stress analysis on absorbers for parabolic solar collectors. Results obtained with this algorithm showed a dramatic rise in deformation as absorber tube length increases. A combined efficiency parameter that includes the occluded area for the mounts was developed to obtain an optimized tube length obtained. Based on the results, a length of 2.7 m for the absorber + 0.2 m for the coupler was chosen to minimize any bending and optimize optical efficiency while maintaining ease of mounting. The associated particle TTS system was designed, built and successfully tested at GTI. It includes storage, receiving and lock hoppers and piping that simulates the transfer of captured solar energy to an actual industrial furnace. Tests over 77 charge-discharge cycles demonstrated <2% particle degradation, with no problematic particle accumulations and no flow interruptions. The piping pressure drop was about 5 psi. The team also worked with Stanley Consultants (Stanley) to prepare conceptual and preliminary engineering packages to facilitate follow-on development and commercialization efforts. These include process and instrumentation diagram’s (P&ID’s), general arrangements, electrical one-line, project definitions document, equipment data sheets, schedule, and construction cost estimate for 2 MWth system. Updated HSS technology commercialization and customer engagement plans and detailed costs and evaluated market trade-offs and manufacturing.

03 NATURAL GAS↗
DOE OSTI2021

Extreme metrics from large ensembles: investigating the effects of ensemble size on their estimates

Abstract. We consider the problem of estimating the ensemble sizes required to characterize the forced component and the internal variability of a number of extreme metrics. While we exploit existing large ensembles, our perspective is that of a modeling center wanting to estimate a priori such sizes on the basis of an existing small ensemble (we assume the availability of only five members here). We therefore ask if such a small-size ensemble is sufficient to estimate accurately the population variance (i.e., the ensemble internal variability) and then apply a well-established formula that quantifies the expected error in the estimation of the population mean (i.e., the forced component) as a function of the sample size n, here taken to mean the ensemble size. We find that indeed we can anticipate errors in the estimation of the forced component for temperature and precipitation extremes as a function of n by plugging into the formula an estimate of the population variance derived on the basis of five members. For a range of spatial and temporal scales, forcing levels (we use simulations under Representative Concentration Pathway 8.5) and two models considered here as our proof of concept, it appears that an ensemble size of 20 or 25 members can provide estimates of the forced component for the extreme metrics considered that remain within small absolute and percentage errors. Additional members beyond 20 or 25 add only marginal precision to the estimate, and this remains true when statistical inference through extreme value analysis is used. We then ask about the ensemble size required to estimate the ensemble variance (a measure of internal variability) along the length of the simulation and – importantly – about the ensemble size required to detect significant changes in such variance along the simulation with increased external forcings. Using the F test, we find that estimates on the basis of only 5 or 10 ensemble members accurately represent the full ensemble variance even when the analysis is conducted at the grid-point scale. The detection of changes in the variance when comparing different times along the simulation, especially for the precipitation-based metrics, requires larger sizes but not larger than 15 or 20 members. While we recognize that there will always exist applications and metric definitions requiring larger statistical power and therefore ensemble sizes, our results suggest that for a wide range of analysis targets and scales an effective estimate of both forced component and internal variability can be achieved with sizes below 30 members. This invites consideration of the possibility of exploring additional sources of uncertainty, such as physics parameter settings, when designing ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Magnetic field-temperature phase diagram of the spin-$\frac{1}{2}$ triangular lattice antiferromagnet KYbSe 2

A quantum spin liquid (QSL) is a state of matter characterized by fractionalized quasiparticle excitations, quantum entanglement, and a lack of long-range magnetic order. However, QSLs have evaded definitive experimental observation. Several Yb 3+ -based triangular lattice antiferromagnets with effective 𝑆 = $\frac{1}{2}$ have been suggested to stabilize the QSL state as the ground state. Here, in this work, we build a comprehensive magnetic temperature phase diagram of a high-quality single crystalline KYbSe 2 via heat capacity and magnetocaloric effect down to 30 mK with magnetic field applied along the 𝑎 axis. At zero magnetic field, we observe the magnetic long-range order at 𝑇 N =0.29 K entering 120 degrees ordered state in heat capacity, consistent with neutron scattering studies. Analysis of the low-temperature (𝑇) specific heat (𝐶) at zero magnetic field indicates linear 𝑇 dependence of 𝐶/𝑇 and a broad hump of 𝐶/𝑇 in the proximate QSL region above 𝑇 N . By applying magnetic field, we observe the up-up-down phase with 1/3 magnetization plateau and oblique phases, in addition to two new phases. These observations strongly indicate that while KYbSe 2 closely exhibits characteristics resembling an ideal triangular lattice, deviations may exist, such as the effect of the next-nearest-neighbor exchange interaction, calling for careful consideration for spin Hamiltonian modeling. Further investigations into tuning parameters, such as chemical pressure, could potentially induce an intriguing QSL phase in the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Completion design improvement using a deep convolutional network

Maximizing stimulated natural and hydraulic fracture network is one of the primary hydraulic fracturing concerns for economic production from a horizontal shale gas well. Geomechanical facies and preexisting fractures in each stage are identified based on similarities in formation characteristics to optimize the locations of perforation clusters. This often requires analyzing large volumes of drilling, Logging While Drilling (LWD) and Measurement While Drilling (MWD) data. In this paper, we develop a methodology that calculates the mechanical specific energy (MSE) using real-time drill string acceleration signals directly from its definition. High resolution vibration signals have been collected using a tri-axial accerlometer, which was an auxiliary tool included in acoustic borehole imager. This technique provides a cost-efficient solution for engineered completion design. Furthermore, we adopt deep Convolutional Neural Network (CNN) with signal processing to build a data pipeline that effectively extracts patterns from dynamic acceleration signals for rock lateral MSE classification. First, we apply discrete wavelet transform and Short-Time Fourier Transform (STFT) for signal denoising and pattern recognition. Then we construct an image dataset using multi-scale image fusion at pixel level from 3 sensor channels, including axial, lateral acceleration spectrograms and zero-padded revolutions per minute (RPM). The resulted RGB image dataset includes 4,000 images of 5 MSE ranges with various rock strength conditions. Our results demonstrate that the proposed deep learning model can achieve more than 90% classification accuracy. The deep learning results, as a reference source, were applied in selected Marcellus Shale Energy and Environmental Lab (MSEEL) wells engineered completion located in the Marcellus shale gas site.

03 NATURAL GAS↗
DOE OSTI2026

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

U.S. Offshore Pipeline and Reported Incident Datasets

The U.S. Offshore Pipeline and Reported Incident Datasets provide a compilation of data from a variety of credible resources, spatially-temporally integrated into multivariate resources. This spatial resource includes more than 80,000 points along existing and abandoned pipelines in the Gulf with matched incidents based on similar lease blocks and temporal timelines (e.g., the incident date occurs within reported pipeline lifespan), structural characteristics, geologic and seafloor data, and meteorological, oceanographic, and biochemical statistics spatially and temporally matched to each point. This is provided as both a feature class in a file geodatabase, as well as a CSV file for ease of use. The pipeline incidents table is a CSV file containing more than 900 reported incidents from 1986 to 2021, including incident date, area (Outer Continental Shelf (OCS) lease block and area code), reported causes, reported incident information, and results (i.e., cost, repairs, inspections), along with quantitative severity metrics. Field dictionaries are included for both the pipeline locations and incidents datasets, which detail field definitions. The pipeline locations field dictionary includes original resource reference information.

Advanced Infrastructure Integrity Model↗
DOE OSTI2025

LandScan HD: a high-resolution gridded ambient population methodology for the world

Unwarned population distributions accounting for routine human activities are needed to address many global human security challenges, including disasters, conflict, and infrastructure demand. LandScan High Definition (LSHD) supports this need through gridded ambient population estimates that measure average human presence between daytime and nighttime at a high spatial resolution of 3 arcseconds (approximately 90 m). Although LSHD has traditionally been produced on a country-specific basis, advances in global foundational data and computational resources now enable scaling its methodology to the world. Combining aspects of top-down and bottom-up gridded population methods, LSHD allocates subnational population totals from authoritative statistics to built-up areas based on occupancy estimates for multiple facility types (e.g., residential, commercial) and then reaggregates these estimates to a global population grid. We scale this approach by organizing the LSHD data stack into a 1° resolution tileset of vector analytic features, enabling an efficient and repeatable workflow for all countries worldwide. Examining the Philippines as an output of the global LSHD baseline dataset, we contrast unwarned and residential (WorldPop) population distributions by (1) exploring a practical application of flood risk assessment and (2) evaluating their congruence with outcomes of collective human activities (subnational CO 2 emissions). Finally, we discuss plans to address current LSHD limitations through data/modeling and uncertainty quantification improvements and provide outlook for workflow automation and extending the model to social, demographic and economic population characteristics.

Building morphology↗
DOE OSTI2021

Spectral properties and enhanced superconductivity in renormalized Migdal-Eliashberg theory

Migdal-Eliashberg theory describes the properties of the normal and superconducting states of electron-phonon-mediated superconductors based on a perturbative treatment of the electron-phonon interactions. It is necessary to include both electron and phonon self-energies self-consistently in Migdal-Eliashberg theory in order to match numerically exact results from determinantal quantum Monte Carlo in the adiabatic limit. Here in this work we provide a method to obtain the real-axis solutions of the Migdal-Eliashberg equations with electron and phonon self-energies calculated self-consistently. Our method avoids the typical challenge of computing cumbersome singular integrals on the real axis and is numerically stable and exhibits fast convergence. Analyzing the resulting real-frequency spectra and self-energies of the two-dimensional Holstein model, we find that self-consistently including the lowest-order correction to the phonon self-energy significantly affects the solution of the Migdal-Eliashberg equations. The calculation captures the broadness of the spectral function, renormalization of the phonon dispersion, enhanced effective electron-phonon coupling strength, minimal increase in the electron effective mass, and the enhancement of superconductivity which manifests as a superconducting ground state despite strong competition with charge-density-wave order. We discuss surprising differences in two common definitions of the electron-phonon coupling strength derived from the electron mass and the density of states, quantities which are accessible through experiments such as angle-resolved photoemission spectroscopy and electron tunneling. An approximate upper bound on 2Δ/T c for conventional superconductors mediated by retarded electron-phonon interactions is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Compactly‐Supported Nonstationary Kernels for Computing Exact Gaussian Processes on Big Data

The Gaussian process (GP) is a widely used method for analyzing large-scale data sets, including spatio-temporal measurements of nonlinear processes that are now commonplace in the environmental sciences. Traditional implementations of GPs involve stationary kernels (also termed covariance functions) that limit their flexibility, and exact methods for inference that prevent application to data sets with more than about 10,000 points. Modern approaches to address stationarity assumptions generally fail to accommodate large data sets, while all attempts to address scalability focus on approximating the Gaussian likelihood, which can involve subjectivity and lead to inaccuracies. In this work, we explicitly derive an alternative kernel that can discover and encode both sparsity and nonstationarity. We embed the kernel within a fully Bayesian GP model and leverage high-performance computing resources to enable the analysis of massive data sets. We demonstrate the favorable performance of our novel kernel relative to existing exact and approximate GP methods across a variety of synthetic data examples. Furthermore, we conduct space–time prediction based on more than 1 million measurements of daily maximum temperature and verify that our results outperform state-of-the-art methods in the Earth sciences. More broadly, having access to exact GPs that use ultra-scalable, sparsity-discovering, nonstationary kernels allows GP methods to truly compete with a wide variety of machine learning methods.

Gaussian processes↗
DOE OSTI2022

Main and Satellite Features in the Ni 2p XPS of NiO

The origin and assignment of the complex main and satellite XPS features of the cations in ionic compounds has been the subject of extensive theoretical studies using different methods. There is agreement that within a molecular orbital model one needs to take into account different types of configurations. Specifically, those where a core electron is removed but no other configuration changes are made and those where in addition to ionization there are also shake or charge transfer changes to the ionic configuration. However, there are strong disagreements about the assignment of XPS features to these configurations. The present work is directed toward resolving the origin of main and satellite features for the Ni 2p XPS of NiO based on ab initio molecular orbital wavefunctions for a cluster model of NiO. A major problem in earlier ab initio XPS studies of ionic compounds has been the use of a common set of orbitals that was not able to properly describe all the ionic configurations that contribute to the full XPS spectra. This is resolved in the present work by using orbitals that are optimized for averages of the occupations of the different configurations that contribute to the XPS. The approach of using State Averaged orbitals is validated through comparisons between different averages and through use of higher order excitations in the wavefunctions for the ionic states. It represents a major extension of our earlier work on the main and satellite features of the Fe 2p XPS of Fe 2 O 3 and proves the reliability and the generality of the assignments of the character and origin of the different features of the XPS obtained with orbitals optimized for State Averages. These molecular orbital methods permit the characterization of the ionic states in terms of the importance of shake excitations and of the coupling of ionization of 2p 1/2 and 2p 3/2 spin-orbit split sub shells. Further, the work lays the foundation for definitive assignments of the character of main and satellite XPS features and points to their origin in the electronic structure of the material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

On de Sitter spacetime and string theory

We review various aspects of de Sitter spacetime in string theory: its status as an Effective Field Theory spacetime solution, its relation to the vacuum energy problem in string theory, its (global) holographic definition in terms of two entangled and noncanonical conformal field theories as well as a realization of a realistic de Sitter universe endowed with the observed visible matter and the necessary dark sector in order to reproduce the realistic cosmological structure. In particular, based on the new insight regarding the cosmological constant problem in string theory, we argue that in a doubled, [Formula: see text]-duality-symmetric, phase-space-like and noncommutative generalized-geometric formulation, string theory can naturally lead to a small and positive cosmological constant that is radiatively stable and technically natural. Such a formulation is fundamentally based on a quantum spacetime, but in an effective spacetime description of this general formulation of string theory, the curvature of the dual spacetime is the cosmological constant of the observed spacetime, while the size of the dual spacetime is the gravitational constant of the same observed spacetime. Also, the three scales associated with intrinsic noncommutativity of string theory, the cosmological constant scale, the Planck scale as well as the Higgs scale, can be arranged to satisfy various seesaw-like formulae. Along the way, we show that these new features of string theory can be implemented in a particular deformation of cosmic-string-like models.

Astronomy & Astrophysics↗
DOE OSTI2021

Investigation of the background in coherent J/ψ production at the EIC

Understanding various fundamental properties of nucleons and nuclei is among the most important scientific goals at the upcoming Electron-Ion Collider (EIC). With the unprecedented opportunity provided by the next-generation machine, the EIC might provide definitive answers to many standing puzzles and open questions in modern nuclear physics. We investigate one of the golden measurements proposed at the EIC, which is to obtain the spatial gluon density distribution within a lead (Pb) nucleus. The proposed experimental process is the exclusive J/ψ vector-meson production off the Pb nucleus: e + Pb → e' + J/ψ + Pb'. The Fourier transformation of the momentum transfer |t| distribution of the coherent diffraction is the transverse gluon spatial distribution. In order to measure it, the experiment has to overcome an overwhelmingly large background arising from the incoherent diffractive production, where the nucleus Pb' mostly breaks up into fragments of particles in the far-forward direction close to the hadron-going beam rapidity. We systematically study the rejection of incoherent J/ψ production by vetoing products from these nuclear breakups—protons, neutrons, and photons—which is based on the BeAGLE event generator and the most up-to-date EIC Far-forward Interaction Region design. The achieved vetoing efficiency, the ratio between the numbers of vetoed events and total incoherent events, ranges from about 80% to 99% depending on |t|. Assuming a 5% smearing applied to the reconstructed |t| resolution in the Sar t re model, this vetoing efficiency can suppress the incoherent background to at least the first minimum of the coherent |t| distribution. Experimental and accelerator machine challenges as well as potential improvements are discussed herein.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Sensitivities of Large Eddy Simulations of Aerosol Plume Transport and Cloud Response

Abstract Cloud responses to surface‐based sources of aerosol perturbation partially depend on how turbulent transport of the aerosol to cloud base affects the spatial and temporal distribution of aerosol. Here, scenarios of plume injection below a marine stratocumulus cloud are modeled using large eddy simulations coupled to a prognostic bulk aerosol and cloud microphysics scheme. Both passive plumes, consisting of an inert tracer, and active plumes are investigated, where the latter are representative of saltwater droplet plumes such as have been proposed for marine cloud brightening. Passive plume scenarios show higher in‐plume cloud brightness (relative to out‐of‐plume) due to the predominant transport of the passive plume tracer from the near‐surface to the cloud layer within updrafts. These updrafts rise into brighter areas within the cloud deck, even in the absence of an aerosol perturbation associated with an active plume. Comparing albedo at in‐plume to out‐of‐plume locations associates the inert plume with the brightest cloud locations, without any causal effect of the plume on the cloud. Numerical sensitivities are first assessed to establish a suitable model configuration. Then sensitivity to particle injection rate is investigated. Trade‐offs are identified between the number of injected particles and the suppressive effect of droplet evaporation on plume loft and spread. Furthermore, as the near‐field in‐plume brightening effect does not depend significantly on injection rate given a suitable definition of perturbed versus unperturbed regions of the flow, plume area is a key controlling factor on the overall cloud brightening effect of an aerosol perturbation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Difference in load predictions obtained with effective turbulence vs. a dynamic wake meandering modeling approach

Abstract. According to the international standard for wind turbine design, the effects of wind turbine wakes on structural loads can be considered in two ways: (1) by augmenting the ambient turbulence levels with the effective turbulence model (EFF) and then calculating the resulting loads and (2) by performing dynamic wake meandering (DWM) simulations, which compute wake effects and loads for all turbines on a farm at once. There is no definitive answer in scientific literature as to the consequences of choosing one model over the other, but the two approaches are unarguably very different. The work presented here expounds on these differences and investigates to what extent they affect the simulated structural loads. We consider an idealized 4×4 rectangular array of National Renewable Energy Laboratory 5 MW wind turbines with a spacing of 5 by 8 rotor diameters and three wind speed scenarios at high ambient turbulence. Load simulations are performed in OpenFAST with EFF and in FAST.Farm with the DWM implementation. We compare ambient turbulence, wind farm turbulence, and loads between both approaches. When omnidirectional results are compared, EFF wind farm turbulence intensity is consistently higher by 0.2 % (above-rated wind speed) to 2.7 % (below-rated wind speed). However, for certain wind directions, the EFF turbulence can be lower than FAST.Farm by almost 9 %. Wind speeds within the farm were found to differ by up to 3 m s−1 due to the lack of wake deficits in the EFF approach, leading to longer tails toward low values in the FAST.Farm mean load distributions. Consistent with the turbulence results, the median EFF load standard deviations are also consistently higher, by a maximum of 20 % and 17 % for blade-root out-of-plane and tower-base fore-aft moments, respectively.

17 WIND ENERGY↗
DOE OSTI2021

A dynamical mechanism for the Page curve from quantum chaos

If the evaporation of a black hole formed from a pure state is unitary, the entanglement entropy of the Hawking radiation should follow the Page curve, increasing from zero until near the halfway point of the evaporation, and then decreasing back to zero. The general argument for the Page curve is based on the assumption that the quantum state of the black hole plus radiation during the evaporation process is typical. In this paper, we show that the Page curve can result from a simple dynamical input in the evolution of the black hole, based on a recently proposed signature of quantum chaos, without resorting to typicality. Our argument is based on what we refer to as the “operator gas” approach, which allows one to understand the evolution of the microstate of the black hole from generic features of the Heisenberg evolution of operators. One key feature which leads to the Page curve is the possibility of dynamical processes where operators in the “gas” can “jump” outside the black hole, which we refer to as void formation processes. Such processes are initially exponentially suppressed, but dominate after a certain time scale, which can be used as a dynamical definition of the Page time. In the Hayden-Preskill protocol for young and old black holes, we show that void formation is also responsible for the transfer of information from the black hole to the radiation. We conjecture that void formation may provide a microscopic explanation for the recent semi-classical prescription of including islands in the calculation of the entanglement entropy of the radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Geophysical and Environmental Monitoring Data, and Subsurface Flow Modelling Results for Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO

This dataset includes geoelectrical monitoring data acquired between October 2021 and November 2022, soil moisture and temperature data, groundwater data obtained from borehole SNIB covering the period from June 2021 to September 2022, and hydrological modelling results. The data were acquired to investigate how variations in bedrock type and topography, and vegetation cover control subsurface flow dynamics. To provide insights into the subsurface flow dynamics and their controls, a monitoring transect was installed at the Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO, measuring the spatio-temporal variations of soil moisture, soil and snow temperature, subsurface electrical resistivity variations, and groundwater dynamics. Field data are organized in a folder structure, with Electrical Resistivity Tomography (ERT) data being provided as one file per measurement, and data of the soil moisture and temperature sensors being provided as text files covering the entire monitoring period. The ‘Locations.csv’ file contains the location of all sensors, given in NAD83 – UTM Zone 13N. ERT monitoring data has been processed to filter data based on reciprocal errors (data with errors > 30% were removed), a linear error model was fitted to each survey, and to ensure a constant set of measurements for time-lapse inversion, filtered data were interpolated and assigned a 100% measurement error. Soil moisture and temperature data were acquired at 15 min intervals, and averaged to provide 1h data. Weather data and borehole data (groundwater depth, conductivity and temperature) were acquired at 30 min intervals, and are provided as daily measurements; all measurements are averaged, except of precipitation values, which are given as daily accumulation. The hydrological model was set up along the ERT monitoring transect, and net infiltration was used as surface boundary condition and derived from the weather data. Four different results are provided, (1) results for a parameterization using hydraulic permeability and porosity as derived from the ERT data through petrophysical relationships, and (2) three simplified model results, using 1 to 3 geological layers above the bedrock. Modelling was performed using PFLOTRAN, and for each model the PFLOTRAN input files are provided. The result files include weekly hydrological modelling results (e.g., saturation, velocities, pressures), as well as the model parameterization. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High dopant activation in arsenic doped single-crystal CdTe thin films: Insights from MBE growth and rapid thermal processing

Single-crystal model systems are valuable tools to investigate fundamental material properties. In this work, we use molecular beam epitaxy to deposit in situ arsenic (As) doped single-crystal CdTe films on large area Si substrates to better understand As doping for photovoltaic applications. We found that As incorporation is highly temperature dependent: a substrate temperature difference of 50 °C can lead to several orders of magnitude difference in As concentration. Cd overpressure during in situ doping may limit out-diffusion of As but decrease As incorporation, especially at lower growth temperatures. Carrier concentrations greater than 10 16 cm −3 can be achieved with or without Cd overpressure when annealed at temperatures above 500 °C. However, unlike the low (∼1% to 5%) dopant activation commonly observed in polycrystalline CdTe, our films achieve significantly higher activation ratios—exceeding 50%, and in some cases approaching 80%. These values are consistent with or exceed prior reports in single-crystal CdTe systems. In addition to as-deposited arsenic concentrations, we also consider arsenic distribution after different rapid thermal processing temperatures. We propose a detailed definition and description of how arsenic incorporation is considered and calculated. Due to carrier concentration saturation, As incorporation also needs to be controlled to average levels of 10 17 cm −3 to achieve high activation. These findings suggest that higher annealing temperature regimes may be beneficial to polycrystalline CdTe based PV devices.

36 MATERIALS SCIENCE↗