Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 649 records · Page 36

Geometry-aware training of factorized layers in tensor Tucker format

Reducing parameter redundancies in neural network architectures is crucial for achieving feasible computational and memory requirements during train and inference of large networks. Given its easy implementation and flexibility, one promising approach is layer factorization, which reshapes weight tensors into a matrix format and parameterizes it as the product of two rank-r matrices. However, this family of approaches often requires an initial full-model warm-up phase, prior knowledge of a feasible rank, and it is sensitive to parameter initialization.In this work, we introduce a novel approach to train the factors of a Tucker decomposition of the weight tensors. Our training proposal proves to be optimal in locally approximating the original unfactorized dynamics and stable for the initialization. Furthermore, the rank of each mode is dynamically updated during training.We provide a theoretical analysis of the algorithm, showing convergence, approximation and local descent guarantees. The method's performance is further illustrated through a variety of experiments, showing remarkable training compression rates and comparable or even better performance than the full baseline and alternative layer factorization strategies.

Zangrando, Emanuele [Gran Sasso Science Institute ↗

Asteroseismic modelling of solar-type stars: a deeper look at the treatment of initial helium abundance

ABSTRACT Detailed understanding of stellar physics is essential towards a robust determination of stellar properties (e.g. radius, mass, and age). Among the vital input physics used in the modelling of solar-type stars which remain poorly constrained, is the initial helium abundance. To this end, when constructing stellar model grids, the initial helium abundance is estimated either (i) by using the semi-empirical helium-to-heavy element enrichment ratio, (ΔY/ΔZ), anchored to the standard big bang nucleosynthesis value, or (ii) by setting the initial helium abundance as a free variable. Adopting 35 low-mass, solar-type stars with multiyear Kepler photometry from the asteroseismic ‘LEGACY’ sample, we explore the systematic uncertainties on the inferred stellar parameters (i.e. radius, mass, and age) arising from the treatment of the initial helium abundance in stellar model grids. The stellar masses and radii derived from grids with free initial helium abundance are lower compared to those from grids based on a fixed ΔY/ΔZ ratio. We find the systematic uncertainties on mean density, radius, mass, and age arising from grids which employ a fixed value of ΔY/ΔZ and those with free initial helium abundance to be ∼ 0.9 per cent, ∼ 2 per cent, ∼ 5 per cent, and ∼ 29 per cent, respectively. We report that the systematic uncertainties on the inferred masses and radii arising from the treatment of initial helium abundance in stellar grids lie within the expected accuracy limits of ESA’s PLATO, although this is not the case for the age.

Nsamba, Benard↗

Demonstrating Hierarchical System Development With the Common Community Physics Package Single‐Column Model: A Case Study Over the Southern Great Plains

This study demonstrates a specific application of the hierarchical system development (HSD) approach to investigate, analyze, and attribute model issues within the Unified Forecast System (UFS), with a focus on process isolation. By evaluating a non‐precipitating, shallow cumulus case at the Atmospheric Radiation Measurement Southern Great Plains site in the UFS global forecast against the observation, the investigation identifies a warmer and deeper daytime convective planetary boundary layer (PBL) and misrepresented nocturnal PBL transition. Hypothesis testing, which employs the Common Community Physics Package (CCPP) single‐column model (SCM) and uses the same physics as the UFS global model, confirms that these issues are attributed to the model physics and initialization. Specifically, misrepresented PBL processes are linked to problematic surface condition and a lack of cloud formation, which may stem from deficiencies in PBL and cloud microphysics parameterizations and their interactions. The UFS initial condition contributes to an earlier, excessively collapsed daytime convective boundary layer and a lack of decoupling between the stable boundary layer and residual layer late in the afternoon. This work introduces an avenue for the community to engage with the application of HSD, along with the CCPP and CCPP SCM, to understand the interplay of model physics, disentangle the roles of model components, as well as facilitate model and forecast improvement.

54 ENVIRONMENTAL SCIENCES↗

Evidence of Hexadecapole Deformation in Uranium-238 at the Relativistic Heavy Ion Collider

State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic 238 U + 238 U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding 238 U ions. Here, we show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like 238 U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of 238 U on high-energy collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Milestone 1.2.10: Steady-state H 2 “roll over” point data for aluminum alloys 1100 and 6061

Extended (> 50 years) dry storage is being evaluated by the U.S. Department of Energy (DOE) for the disposition of ~ 18 metric tons of aluminum-clad spent nuclear fuel (ASNF). Transition of the current ASNF inventory into dry storage—using the standard DOE canister—necessitates a rigorous, predictive understanding of the long-term physical and chemical factors that may influence the integrity of the proposed storage canister, including radiolytic molecular hydrogen (H 2 ) generation. Current model predictions employ initial radiolytic yields of H 2 , the values of which change as the cladding’s H 2 -precursor inventory is depleted and H 2 itself becomes progressively more involved in radiolytic and surface dissociation processes. Consequently, the absorbed radiation dose that this steady-state H 2 yield corresponds to is essential for the evaluation and improvement of model predictions. Here, we report our findings on the long-term generation of H 2 from the gamma irradiation (≤ 36 MGy) of corroded AA1100 and AA6061 coupons in helium environments at ambient temperature and ~ 50% RH. Our findings show that AA1100 systems reached steady-state by ~ 36 MGy, while higher doses were necessary for AA6061 systems. This discrepancy was attributed to the AA6061 coupons developing a thicker corrosion layer that led to the trapping of H 2 and its precursors, and potentially additional chemistries, ultimately delaying the depletion of H 2 precursors and the system’s “roll over” point. Further, current model predictions—based on previous AA1100 data—do not show steady-state attainment until above 120 MGy, which is not the case for the AA1100 data collected here. Consequently, the new alloy dependent data presented here are important for the continued improvement of predictive computer models for evaluating the feasibility of extended storage of ASNF in helium backfilled canisters.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

THERMAL STRESS ANALYSIS OF A SPENT NUCLEAR FUEL CANISTER

The potential for chloride induced stress corrosion cracking (CISCC) in spent nuclear fuel dry storage canisters is a current topic of research by the US Department of Energy Spent Fuel and Waste Science and Technology program. One of the important prerequisites for CISCC is a tensile residual stress state. This study utilizes computational models to provide an initial analysis of thermal stresses in a generic spent nuclear fuel canister. A STAR-CCM+ thermal fluid model provides a temperature profile analysis of the canister based on four different ambient temperature conditions. An ANSYS APDL finite element model incorporates the temperature profiles to analyze the thermal stresses in the canister. The calculated thermal stress magnitudes are in the range of 10 MPa to 80 MPa, which could be significant for crack propagation through the canister wall. The maximum thermal stress is not a concern for the structural failure of the canister, but it is high enough that it is a significant contribution to the total stress state of the canister wall, which also includes residual stress from fabrication, internal pressure from helium cover gas, and potential transient mechanical loads, such as earthquakes. This paper describes the initial finite element thermal stress analysis that was completed in 2020, reports the initial findings, and identifies areas where model refinement is still needed to complete this work in the future.

Jensen, Ben J.↗

Chromium Plume at LANL - From Conceptual Site Model to Development of Complex Remediation Strategies - 18682

A hexavalent chromium plume is present in deep groundwater beneath Los Alamos National Laboratory (LANL or the Laboratory). The plume originated from releases of potassium dichromate from cooling towers from 1956 - 1972. Estimates are that up to 72,000 kg of chromium was released during that period. Data from surface water-balance studies, core-hole and groundwater characterization, geomorphic investigations and geology and geophysics provide the basis for an evolved conceptual model that chromium initially traveled rapidly via an effluent-supported stream for several kilometers, then infiltrated through a stratigraphically complex 275-300 meters thick vadose zone and eventually into a deep (regional) aquifer that is the water supply for the Laboratory and surrounding municipality. The travel time from the release point to the regional aquifer is estimated at 5-10 years based on modeling and contaminant tracers. Monitoring data from two-screen wells and piezometers within the regional aquifer, and data derived from stratified sampling conducted during sonic coring and from analysis of sonic core indicate that the chromium contamination occurs within the upper 23 meters of the highly stratified regional aquifer. Different geologic units at the water table and structural dip of depositional bedding appear to have little effect on chromium migration and plume shape and thickness which instead appear to be a manifestation of multiple breakthrough locations, interconnectedness of preferential hydraulic strata, and a negligible vertical gradient. Recent monitoring data indicate that the down-gradient plume edge (as defined by the New Mexico groundwater standard of 50 μg/L) may be advancing near the Laboratory's boundary with a neighboring Indian reservation. The Laboratory is addressing the increasing concentrations with a regulatory-driven interim measure that uses hydraulic control to manage plume migration. The interim measure involves limited-duration (2-3 years) pumping, treatment, and injection to rapidly establish and maintain the 50 μg/L plume edge within the Laboratory boundary while a final remediation strategy is developed and approved by the State regulator. A final remediation strategy will depend heavily on the ability to cost-effectively access the deep aquifer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Spatially Explicit Individual-Based Modeling Approach to Evaluate the Cumulative Effects of Wind Energy Development on the Greater Sage-Grouse: Sensitivity Analysis and Validation of Key Parameters

To address the need for an approach to evaluate cumulative ecological impacts of wind energy development, particularly those on critical wildlife habitats and species, the authors developed a proof-of-concept, landscape-based, spatially explicit individual-based modeling (IBM) framework for populations of the greater sage-grouse (Centrocercus urophasianus) in Albany County, Wyoming, based on published and other available information on the life history of the species. This sage-grouse IBM estimates the spatial and temporal movements of sage-grouse based on modeled effects of habitat characteristics and developments on sage-grouse condition, survivorship, and reproduction. After completing the initial version of the model, we conducted a series of sensitivity analyses of key model parameters and validation of lek occurrences and spatial distributions to better understand performance of the model relative to observed values of important life history characteristics. The parameter sensitivity analysis was conducted as a series of one-factor-at-a-time tests in which the values of selected parameters were individually varied to evaluate the corresponding effects on model output. We examined the model’s response to changes in parameter values to gain insight into the overall robustness of the model, identify key driver parameters, and initiate qualitative validation. Five functions selected for sensitivity analyses were (1) the rate of change in body condition in relationship to habitat suitability index (HSI) of occupied habitat (condition change rate), (2) the effect of competition between individuals with overlapping home ranges expressed as a percentage reduction in HSI (competition factor), (3) value of condition at which survivorship probability is equal to zero (zero survivorship condition), (4) the relationship between female body condition and clutch size, and (5) the relationship between female body condition and nest success. We varied the parameter values to be tested for the sensitivity analysis around the original value used in the model (i.e., tested values were higher and lower than the model’s original value) to examine the effect of this variation on the predicted output of the model. Model outputs evaluated in the sensitivity analysis included (1) predicted age-class distribution (as measured by the percentage of immature individuals within the population); (2) predicted sex ratio; (3) predicted age-sex class distribution (percent of population in each of six age-sex classes, i.e., male and female juveniles, yearlings, and adults); (4) lifespan; (5) population size; and (6) spatial distribution. Of the five parameters analyzed, condition change rate, survivorship, clutch size, and nest success did not substantially affect any of the evaluated model outputs. For the predicted age-class distribution, the percentage of immature individuals predicted by the model for these four parameters ranged from 61.9% to 65.0%, which was reasonable with respect to the reference values of 51.4 to 57.8% in published studies. The ratio of females to males ranged from 1.09 to 1.14, which was comparable to the reference values of 1.2 to 3.0 presented in an existing study. The predicted mean lifespan ranged from 1.62 to 1.74 years, which appears to be reasonable with respect to the reference values of 0.9 to 1.1 years for sharp-tail grouse and greater prairie chickens presented by published studies. The predicted population size ranged from 1,877 to 3,140. This prediction would be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county. The predicted spatial distribution was not substantially affected by the parameter values examined. In our sensitivity analyses, the only tested parameter that had a noticeable effect on model output was competition factor. Over the range of values tested for this parameter, predicted population size exhibited the largest range (from 1,833 to 4,737) of all parameters tested. As the competition factor increased, the predicted population size decreased. This reflects the effect competition had on the number of individuals that could co occur in high-quality habitat patches. Despite this apparently significant effect, the range of values appears to be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county.

17 WIND ENERGY↗

Evaluating skill in predicting the Interdecadal Pacific Oscillation in initialized decadal climate prediction hindcasts in E3SMv1 and CESM1 using two different initialization methods and a small set of start years

Abstract It is a daunting challenge to conduct initialized hindcasts with enough ensemble members and associated start years to form a drifted climatology from which to compute the anomalies necessary to quantify the skill of the hindcasts when compared to observations. This limits the ability to experiment with case studies and other applications where only a few initial years are needed. Here we run a set of hindcasts with CESM1 and E3SMv1 using two different initialization methods for a limited set of start years and use the respective uninitialized free-running historical simulations to form the model climatologies. Since the drifts from the observed initial states in the hindcasts toward the uninitialized model state are large and rapid, after a few years the drifted initialized models approach the uninitialized model climatological errors. Therefore, hindcasts from the limited start years can use the uninitialized climatology to represent the drifted model states after about lead year 3, providing a means to compute forecast anomalies in the absence of a large hindcast sample. There is comparable skill for predicting spatial patterns of multi-year Pacific sea surface temperature anomalies in the domain of the Interdecadal Pacific Oscillation using this method compared to the conventional methodology with a large hindcast data set, though there is a model dependence to the drifts in the two initialization methods.

54 ENVIRONMENTAL SCIENCES↗

Automated Model Updating Effect on the Linear and Nonlinear Dynamic Responses of Historical Masonry Structures

Finite element models supply the simulation of the problems to obtain the structural behaviors under different loading conditions. The updating processes based on changing the uncertain parameters such as material properties, boundary conditions, meshing size etc. are made to create actual finite element models in order to represent the real conditions of the problems. This procedure can be made using two methods, namely manual model updating and global/local automatic model updating. The global/local automatic model updating procedures are very popular owing to the applicability to all kind of engineering structures, minimizing the differences nearly to zero, availability of damage localization and structural health monitoring. In this paper a historical masonry armory building selected as an application and its finite element analyses, experimental measurements and automated model updating processes are carried out. The finite element model of the structures was developed with ANSYS software and first three natural frequencies are attained. The experimental dynamic characteristics such as natural frequencies, damping ratios and mode shapes are obtained using ambient vibration method. The differences between experimental and numerical frequencies are minimized up to the range of 0–1%, by using local updating procedure. The linear and nonlinear seismic analyses of the building are conducted on initial and updated models. At the end of the analyses, maximum displacements, principal stresses and strains are given with detail. In addition, crack distributions diagrams are attained to validate the stresses accumulation points. It can be seen that there is an increasing trend in the displacement, stress and strain values after finite element model updating both linear and nonlinear analysis. In addition, linear and nonlinear analysis results are compared with each other to evaluate the nonlinear analysis effect. It is concluded that the displacement and strain values are increased after nonlinear analysis. But, there is not any agreement in stress values.

42 ENGINEERING↗

Monte Carlo event generator for initial conditions of conserved charges in nuclear geometry

At top collider energies where baryon stopping is negligible, the initial state of heavy-ion collisions is overall charge neutral and predominantly composed of gluons. Nevertheless, there can also be significant local fluctuations of the baryon number, strangeness, and electric charge densities about zero, perturbatively corresponding to the production of quark/antiquark pairs. These previously ignored local charge fluctuations can permit the study of charge diffusion in the quark-gluon plasma (QGP), even at top collider energies. In this paper we present a new model denoted ICCING (initial conserved charges in nuclear geometry) which can reconstruct the initial conditions of conserved charges in the QGP by sampling a (g→q¯q) splitting probability over the initial energy density. We find that the new charge distributions generally differ from the bulk energy density; in particular, the strangeness distribution is significantly more eccentric than standard bulk observables and appears to be associated with the geometry of hot spots in the initial state. The new information provided by these conserved charges opens the door to studying a wealth of new charge- and flavor-dependent correlations in the initial state and ultimately the charge transport parameters of the QGP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Combustion Characteristics of a Tier II Gasoline Certification Fuel and its Surrogate with Iso-butanol: experiments and detailed numerical modeling

This presentation will describe results of an experimental and numerical study of combustion of a Tier II gasoline certification fuel and its surrogate blended with iso-butanol as a potential biofuel additive. The experiments were carried out under conditions where gas transport originated entirely by fuel evaporation with external convection (both forced and natural) were minimized. The resulting near spherically symmetric droplet flames were amenable to detailed numerical modeling. The droplet initial diameters were in the narrow range of 0.55 mm to 0.62 mm while and burning is performed under standard atmosphere. The combustion characteristics such as evolution of droplet, flame and soot shell diameters for indolene/iso-butanol and surrogate/iso-butanol mixtures are compared. The presence of iso-butanol is experimentally found to reduce the formation of soot as the iso-butanol concentration was increased. Simulations are also reported using a detailed kinetic chemistry comprised of 324 species and 17608 reactions using the CRECK POLIMI mechanism. The detailed numerical model included unsteady gas and liquid transport, soot chemistry and radiative transport. Droplet burn rates were well predicted while flame diameters showed more variability depending on how the flame was defined in the numerical simulation (peak gas temperature and OH concentration). Simulations also show the influence of iso-butanol concentration on flame temperature, greenhouse gas emissions and particulate emissions.

Guo, Songtao↗

Constraints on cosmic strings using data from the first Advanced LIGO observing run

Cosmic strings are topological defects which can be formed in grand unified theory scale phase transitions in the early universe. They are also suggested to form in the context of string theory. The main mechanism for a network of Nambu-Goto cosmic strings to lose energy is through the production of loops and the subsequent emission of gravitational waves, thus offering an experimental signature for the existence of cosmic strings. Here we report on the analysis conducted to specifically search for gravitational-wave bursts from cosmic string loops in the data of Advanced LIGO 2015-2016 observing run (O1). No evidence of such signals was found in the data, and as a result we set upper limits on the cosmic string parameters for three recent loop distribution models. Here, we initially derive constraints on the string tension G μ and the intercommutation probability, using not only the burst analysis performed on the O1 data set but also results from the previously published LIGO stochastic O1 analysis, pulsar timing arrays, cosmic microwave background and big-bang nucleosynthesis experiments. We show that these data sets are complementary in that they probe gravitational waves produced by cosmic string loops during very different epochs. Lastly, we determine that the data sets exclude large parts of the parameter space of the three loop distribution models we consider.

79 ASTRONOMY AND ASTROPHYSICS↗

Similarity Criteria for the Study of Removal of Spherical Non-metallic Inclusions in Physical Models of Continuous Casting Tundishes: A More Fundamental Approach

Physical modeling is widely used for studying steelmaking processes. The accuracy of this technique depends on the use of appropriate similarity criteria during experimental design. The present work deduced rigorously the similarity criteria for the study of spherical non-metallic inclusions in isothermal physical models of tundishes. Initially, macroscopic quantities were used to simplify the multiphase problem. Then, the similarity criteria were obtained from the dimensionless form of the exact equation of motion for an inclusion in Stokes regime. The validity range of these criteria was discussed and they were compared with common literature practice. Some practical considerations were also presented. This work clarifies various implicit assumptions usually adopted in literature. Additionally, a restrictive criterion to choose the density of the inclusion simulating particles was obtained, which can be satisfied, for example, by hollow glass particles in water models. This criterion is necessary to adequately describe the inclusion motion in regions where fluid elements suffer strong acceleration (e.g., inside nozzles, near flow controllers or in the entering fluid jet). Moreover, it was established an upper limit for the concentration of particles injected in the inlet nozzle of physical models, which do not simulate the phenomenon of inclusion aggregation.

36 MATERIALS SCIENCE↗

Observations for Model Intercomparison Project (Obs4MIPs): status for CMIP6

The Observations for Model Intercomparison Project (Obs4MIPs) was initiated in 2010 to facilitate the use of observations in climate model evaluation and research, with a particular target being the Coupled Model Intercomparison Project (CMIP), a major initiative of the World Climate Research Programme (WCRP). To this end, Obs4MIPs (1) targets observed variables that can be compared to CMIP model variables; (2) utilizes dataset formatting specifications and metadata requirements closely aligned with CMIP model output; (3) provides brief technical documentation for each dataset, designed for nonexperts and tailored towards relevance for model evaluation, including information on uncertainty, dataset merits, and limitations; and (4) disseminates the data through the Earth System Grid Federation (ESGF) platforms, making the observations searchable and accessible via the same portals as the model output. Taken together, these characteristics of the organization and structure of obs4MIPs should entice a more diverse community of researchers to engage in the comparison of model output with observations and to contribute to a more comprehensive evaluation of the climate models. At present, the number of obs4MIPs datasets has grown to about 80; many are undergoing updates, with another 20 or so in preparation, and more than 100 are proposed and under consideration. A partial list of current global satellite-based datasets includes humidity and temperature profiles; a wide range of cloud and aerosol observations; ocean surface wind, temperature, height, and sea ice fraction; surface and top-of-atmosphere longwave and shortwave radiation; and ozone (O 3 ), methane (CH 4 ), and carbon dioxide (CO 2 ) products. A partial list of proposed products expected to be useful in analyzing CMIP6 results includes the following: alternative products for the above quantities, additional products for ocean surface flux and chlorophyll products, a number of vegetation products (e.g., FAPAR, LAI, burned area fraction), ice sheet mass and height, carbon monoxide (CO), and nitrogen dioxide (NO 2 ). While most existing obs4MIPs datasets consist of monthly-mean gridded data over the global domain, products with higher time resolution (e.g., daily) and/or regional products are now receiving more attention. Along with an increasing number of datasets, obs4MIPs has implemented a number of capability upgrades including (1) an updated obs4MIPs data specifications document that provides additional search facets and generally improves congruence with CMIP6 specifications for model datasets, (2) a set of six easily understood indicators that help guide users as to a dataset's maturity and suitability for application, and (3) an option to supply supplemental information about a dataset beyond what can be found in the standard metadata. With the maturation of the obs4MIPs framework, the dataset inclusion process, and the dataset formatting guidelines and resources, the scope of the observations being considered is expected to grow to include gridded in situ datasets as well as datasets with a regional focus, and the ultimate intent is to judiciously expand this scope to any observation dataset that has applicability for evaluation of the types of Earth system models used in CMIP.

54 ENVIRONMENTAL SCIENCES↗

A Synoptic System for Capturing Ecosystem Control Points Across Terrestrial‐Aquatic Interfaces

Interconnected landscape features such as terrestrial‐aquatic interfaces play an outsized role in biogeochemical cycles as ecosystem control points, but it is notoriously challenging to characterize these. Here, we document a synoptic sensor network design that is (a) flexible to accommodate diverse ecosystem interfaces and gradients, (b) adaptable to monitoring and modeling needs of small and large projects alike, (c) standardized for intercomparability across sites and field experiments, and (d) adequately replicated to capture heterogeneity of each parameter monitored. This real‐time monitoring of surface water, groundwater, soil, and vegetation supports configuration and evaluation of models that span upland, wetland, open water strata, and transitions between them. We established the network at seven sites along the Chesapeake Bay and Lake Erie coastlines, including large‐scale flood manipulation experiments in both regions. A central design element is “one data logger program to rule them all”—a collection of sensor‐specific modules deployed on 40 loggers controlling ∼2,000 sensors, with the goal of streamlining maintenance, debugging, and reproducible data processing. The network generates ∼6 M observations per month, capturing system dynamics at the broad spatial and fine temporal scales needed to initialize and benchmark models; measurement frequency can be modified remotely to capture events. This network design has also revealed behaviors not represented in Earth system models, such as transient groundwater oxygen pulses. Completely documented and open source, this standardized, flexible, and efficient sensor network design can reduce barriers to understanding environmental changes and ecosystem responses across systems and scales.

Ward, Nicholas D. [Pacific Northwest National Labo↗

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗

A knowledge-informed large language model framework for U.S. nuclear power plant shutdown initiating event classification for probabilistic risk assessment

Identifying and classifying shutdown initiating events (SDIEs) is critical for developing shutdown probabilistic risk assessment for nuclear power plants. Existing computational approaches cannot achieve satisfactory performance due to the challenges of unavailable large, labeled datasets, imbalanced event types, and label noise. To address these challenges, we propose a hybrid pipeline that integrates a knowledge-informed machine learning model to prescreen non-SDIEs and a large language model (LLM) to classify SDIEs into four types. In the prescreening stage, we proposed a set of 44 SDIE text patterns that consist of the most salient keywords and phrases from six SDIE types. Text vectorization based on the SDIE patterns generates feature vectors that are highly separable by using a simple binary classifier. The second stage builds Bidirectional Encoder Representations from Transformers (BERT)-based LLM, which learns generic English language representations from self-supervised pretraining on a large dataset and adapts to SDIE classification by fine-tuning it on an SDIE dataset. The proposed approaches are evaluated on a dataset with 10,928 events using precision, recall ratio, F 1 score, and average accuracy. In conclusion, the results demonstrate that the prescreening stage can exclude more than 97% non-SDIEs, and the LLM achieves an average accuracy of 95.1% for SDIE classification.

99 - GENERAL AND MISCELLANEOUS↗