SEARCH · Engineering Papers

Results for “Process modelling”

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,009 records · Page 56

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Conservative force model performance for TOPEX/Poseidon precision orbit determination

The TOPEX/Poseidon spacecraft was launched on August 10, 1992 to study the Earth's oceans. To achieve maximum benefit from the altimetric data collected, mission requirements dictate that TOPEX/Poseidon's orbit must be computed at an unprecedented level of accuracy. In order to satisfy these requirements, a model which accounts for the satellite's complex geometry, attitude, and surface properties has been developed. This `box-wing' representation treats the spacecraft as the combination of flat plates arranged in the shape of a box and a connecetd solar array. The nonconservative forces acting on each of the eight surfaces are computed independently, yielding vector accelerations which are summed to compute the total aggregate effect on the satellite center-of-mass. Parameters associated with each flat plate were derived from a finite element analysis of the spacecraft. Certain parameters can be inferred from tracking data and have been adjusted to obtain a better representation of the satellite acceleration history. Changes in the nominal mission profile and the presence of an `anomalistic' force have complicated this tuning process. Model performance, parameter sensitivities, and the `anomalistic' force will be discussed.

Marshall, J. Andrew↗

DOE OSTI2026

Electrochemical behavior of SnCl 2 and influence of Cu and Ni ions in molten LiCl−KCl−CaCl 2 eutectic

Reliable transport and thermodynamic data for multivalent ions in complex molten salts are scarce, limiting model fidelity for electrorefining and impurity control. Here, we report a comprehensive electrochemical characterization of SnCl₂ in LiCl–KCl–CaCl₂ (50.5–44.2–5.3 mol%) at 685 K, including the effects of Ni 2+ and Cu + impurities. Using cyclic voltammetry (CV), chronoamperometry (CA), and chronopotentiometry (CP), we quantified Sn 2+ and Ni 2+ diffusion with exceptional agreement across methods: Sn 2+ averaged (1.03 ± 0.10) × 10 −5 cm 2 s −1 , and Ni 2+ averaged (0.75 ± 0.19) × 10 −5 cm 2 s −1 . The tight confidence-interval overlap across CV, CA, and CP strengthens confidence in these values and is uncommon in molten chloride studies. Open-circuit-potential measurements provided standard apparent reduction potentials that closely match LiCl–KCl literature, indicating minimal shift with CaCl₂ present. The Sn 2+ /Sn couple behaves as a reversible two-electron soluble–insoluble process at 685 K; the Sn 4+ /Sn 2+ couple transitions to soluble–soluble behavior near 788 K, which may correlate with the decomposition of surface bound chlorostannates, though direct characterization remains to be established. In mixed systems, Cu+/Cu overlaps Sn 2+ /Sn, limiting Cusingle bondSn electroseparation, whereas the larger potential gap between Ni 2+ /Ni and Sn 2+ /Sn supports selective Ni removal. These internally consistent transport and thermodynamic data establish a validated basis for process modeling and optimization of Sn electrorefining and impurity management in LiCl–KCl–CaCl₂.

Berzins-Delahay↗

DOE OSTI2024

Method for Real-Time Salt Aerosol Concentration and Size Measurements for Molten Salt Reactor Safety Assessments

Licensing a new nuclear reactor involves evaluating credible accident scenarios by developing models that simulate accident progression to predict possible outcomes and their impact on safety. The formation of radionuclide-bearing aerosols in the respirable size range can significantly affect dose consequence and threaten human health and is therefore a focus of nuclear reactor safety assessments. Recent reviews of the literature identified a lack of existing experimental data describing the mechanisms of formation and properties (size, concentration, and composition) of radionuclide-bearing aerosols that may be produced from molten salt reactor (MSR) facilities during postulated accident scenarios. Experiments on aerosol formation from molten salt systems are a high priority need that will support MSR licensing by indicating the radiological consequences of aerosol formation and providing the data required for model development and validation. The evolving conditions that occur during a MSR accident, such as a spill of molten fuel salt, may affect aerosol formation mechanisms and aerosol properties over time. There is a need to conduct experiments that simulate credible MSR accidents in a laboratory to generate aerosols with realistic characteristics and behaviors and a need to measure these aerosols in real time under accident-relevant conditions. This report describes the development of a method that can be used to quantify the size and concentration of salt aerosol particles that form from molten salt systems in real time. This method will be employed in future integral effects tests that are conducted at an engineering scale to simulate realistic MSR accidents and in future separate effects tests that will provide mechanistic insight into aerosol formation and properties to support process model development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Deep Search for Decaying Dark Matter with XMM-Newton Blank-Sky Observations

Sterile neutrinos with masses in the keV range are well-motivated extensions to the Standard Model that could explain the observed neutrino masses while also making up the dark matter (DM) of the universe. If sterile neutrinos are DM then they may slowly decay into active neutrinos and photons, giving rise to the possibility of their detection through narrow spectral features in astrophysical x-ray data sets. Here, we perform the most sensitive search to date for this and other decaying DM scenarios across the mass range from 5 to 16 keV using archival XMM-Newton data. We reduce 547 Ms of data from both the MOS and PN instruments using observations taken across the full sky and then use this data to search for evidence of DM decay in the ambient halo of the Milky Way. We determine the instrumental and astrophysical baselines with data taken far away from the Galactic Center, and use Gaussian process modeling to capture additional continuum background contributions. No evidence is found for unassociated x-ray lines, leading us to produce the strongest constraints to date on decaying DM in this mass range.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

IM3 Projected US Data Center Locations

IM3 Projected US Data Center Locations This dataset contains model projections of new data center facilities in the contiguous United States (CONUS) through 2035 using the CERF – Data Centers model. Data center locations are modeled across four data center electricity demand growth scenarios (low, moderate, high, higher) and five market gravity scenarios (0%, 25%, 50%, 75%, 100%). Projected locations are intended to be regional representations of feasible siting locations in the future to assess potential grid and water stress impacts. The data center load growth scenarios correspond with the rates outlined in EPRI (2024) and include 3.71%, 5%, 10%, and 15% annual growth of electricity demand for data centers from 2023 values in 37 states across the CONUS. Market gravity scenarios correspond to the relative importance of proximity to data center markets or high population areas compared to locational cost in the siting algorithm. 0% market gravity means that siting decisions were entirely determined by the locational cost in each feasible location. 100% market gravity means that only market proximity was considered when siting. Other scenarios have weight placed on both components where total weight always equals 100%. Locational cost is dependent on facility cooling type and corresponding electricity cost, taxes, and other factors. Facility cooling type is spatially determined where high water stress and/or areas with high summer wet bulb temperatures are assumed to operate with mechanical cooling for a higher fraction of the year rather than evaporative cooling. Feasible data center siting areas are based on geospatial suitability raster data developed with open-source information. The following areas are excluded from siting: Areas within 300 m of a federal airport runway Waterbodies Areas with slope >16% Areas susceptible to sinkholes High coastal or inland flood risk areas Local, state, and federal parks, leisure areas, and cemeteries Areas >2 km away from electric substations Areas >5 km away from a municipal water supplier service area Areas >2 km away from high-speed fiber provider service territory Protected Areas Database of the United States (PAD-US) areas Railroads, major roadways, and minor roadways Military areas and training grounds NLCD developed lands Areas >0.8 km (0.5 miles) from NLCD developed lands Because we use open-source information, proprietary information that can influence siting decisions such as individual tax agreements with cities, detailed fiber line connectivity, electric grid power capacity agreements, and others, are not currently accounted for in the modeling process. Using specific building locations and footprints in the dataset for local planning purposes is not advised. Technical Information Geospatial data is provided in geojson format using the Albers Equal Area Conic (ESRI:102003) coordinate reference system. The datasets contain the following parameters: id - unique identification number within given scenario file growth_scenario – data center demand growth scenario market_gravity_weight – market gravity weight scenario (%) region – name of region (i.e., US State) total_cost_million_usd – locational siting cost ($million) campus_size_square_ft – total land acquired for data center facility (square ft) data_center_it_power_mw – IT power of data center facility (MW) mechanical_cooling_frac – fraction of year when data center uses mechanical cooling system water_cooling_frac– fraction of year when data center uses evaporative cooling system cooling_energy_demand_mwh – total annual facility energy demand for cooling (MWh) cooling_water_demand_mgy – total annual facility water demand for cooling (MG) cooling_water_consumption_mgy – total annual facility water consumed (MG) normalized_locational_cost – normalized total locational cost score for location normalized_gravity_score – normalized market gravity score for location weighted_siting_score – total weighted siting score of locational cost and gravity score geometry – polygon geometry of facility Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4.0 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall (ORCID:0000000328077088)↗

DOE OSTI2026

IM3 Projected US Data Center Locations

IM3 Projected US Data Center Locations This dataset contains model projections of new data center facilities in the contiguous United States (CONUS) through 2035 using the CERF – Data Centers model. Data center locations are modeled across four data center electricity demand growth scenarios (low, moderate, high, higher) and five market gravity scenarios (0%, 25%, 50%, 75%, 100%). Projected locations are intended to be regional representations of feasible siting locations in the future to assess potential grid and water stress impacts. The data center load growth scenarios correspond with the rates outlined in EPRI (2024) and include 3.71%, 5%, 10%, and 15% annual growth of electricity demand for data centers from 2023 values in 37 states across the CONUS. Market gravity scenarios correspond to the relative importance of proximity to data center markets or high population areas compared to locational cost in the siting algorithm. 0% market gravity means that siting decisions were entirely determined by the locational cost in each feasible location. 100% market gravity means that only market proximity was considered when siting. Other scenarios have weight placed on both components where total weight always equals 100%. Locational cost is dependent on facility cooling type and corresponding electricity cost, taxes, and other factors. Facility cooling type is spatially determined where high water stress and/or areas with high summer wet bulb temperatures are assumed to operate with mechanical cooling for a higher fraction of the year rather than evaporative cooling. Feasible data center siting areas are based on geospatial suitability raster data developed with open-source information. The following areas are excluded from siting: Areas within 300 m of a federal airport runway or within an airport area boundary Waterbodies Areas with slope >16% Areas susceptible to sinkholes High coastal or inland flood risk areas Local, state, and federal parks, leisure areas, and cemeteries Areas >2 km away from electric substations Areas >5 km away from a municipal water supplier service area Areas >2 km away from high-speed fiber provider service territory USGS Protected Areas Database of the United States (PAD-US) GAP status 1, 2, or 3 areas US National Parks Wetlands USFWS critical habitats BIA land areas Railroads, major roadways, and minor roadways Military areas and training grounds NLCD developed lands Areas >0.8 km (0.5 miles) from NLCD developed lands Because we use open-source information, proprietary information that can influence siting decisions such as individual tax agreements with cities, detailed fiber line connectivity, electric grid power capacity agreements, and others, are not currently accounted for in the modeling process. Using specific building locations and footprints in the dataset for local planning purposes is not advised. Technical Information Geospatial data is provided in geojson format using the Albers Equal Area Conic (ESRI:102003) coordinate reference system. The datasets contain the following parameters: id - unique identification number within given scenario file growth_scenario – data center demand growth scenario market_gravity_weight – market gravity weight scenario (%) region – name of region (i.e., US State) total_cost_million_usd – locational siting cost ($million) campus_size_square_ft – total land acquired for data center facility (square ft) data_center_it_power_mw – IT power of data center facility (MW) mechanical_cooling_frac – fraction of year when data center uses mechanical cooling system water_cooling_frac– fraction of year when data center uses evaporative cooling system cooling_energy_demand_mwh – total annual facility energy demand for cooling (MWh) cooling_water_demand_mgy – total annual facility water demand for cooling (MG) cooling_water_consumption_mgy – total annual facility water consumed (MG) normalized_locational_cost – normalized total locational cost score for location normalized_gravity_score – normalized market gravity score for location weighted_siting_score – total weighted siting score of locational cost and gravity score geometry – polygon geometry of facility Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4.0 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall (ORCID:0000000328077088)↗

DOE OSTI2026

Tabulated data for Lignin-Derived Phenolic Compounds and Water Are Effective Cosolvents for Reductive Catalytic Fractionation

Reductive catalytic fractionation (RCF) is an effective lignin-first biorefining method to extract lignin as a stabilized oil from lignocellulosic biomass. To realize RCF at scale, process modeling has shown that minimizing the use of exogenous organic solvents is critical. To this end, here we investigate the ability of lignin-derived monomers to act as either solvents or cosolvents for RCF. We begin by examining the influence of lignin-derived aromatic compounds (4-propylguaiacol, 4-propylphenol, and propylbenzene) on RCF monomer yields and subsequently extend our analysis to mixtures of 4-propylguaiacol and either methanol or water. We demonstrate that 4-propylguaiacol is an effective solvent for lignin extraction and depolymerization during RCF, especially when used in combination with water as a cosolvent. Cosolvent mixtures of 4-propylguaiacol and water enable up to 81% lignin extraction, monomer yields up to 25 wt %, and postreaction phase separation. However, unlike methanol, water as a cosolvent fails to inhibit aromatic ring hydrogenation when conducted over Ru/C as a catalyst, potentially leading to excess hydrogen consumption in a process utilizing this approach. Nonetheless, these results suggest a promising strategy for eliminating external organic solvents from RCF by utilizing mixtures of lignin-derived compounds and water as alternative extraction solvents. This is the tabulated data for the paper

biorefining↗

DOE OSTI2026

Calculation of monomer yields, verification of PG concentration approximation, catalyst comparison for neat aromatic RCF, SEM micrographs for postreaction pulps, yields of aqueous-phase sugar derivatives after PG-H2O RCF, reactor pressure transients, GPC comparing pre- and postreaction organic solvent composition, PG-H2O miscibility testing, GPC probing dimer/oligomer fractionation in PG-H2O mixtures for

biorefining↗

DOE OSTI2026

Lignin-Derived Phenolic Compounds and Water Are Effective Cosolvents for Reductive Catalytic Fractionation

biorefining↗

DOE OSTI2025

Lignin-Derived Phenolic Compounds and Water Are Effective Cosolvents for Reductive Catalytic Fractionation

Reductive catalytic fractionation (RCF) is an effective lignin-first biorefining method to extract lignin as a stabilized oil from lignocellulosic biomass. To realize RCF at scale, process modeling has shown that minimizing the use of exogenous organic solvents is critical. To this end, here we investigate the ability of lignin-derived monomers to act as either solvents or cosolvents for RCF. We begin by examining the influence of lignin-derived aromatic compounds (4-propylguaiacol, 4-propylphenol, and propylbenzene) on RCF monomer yields and subsequently extend our analysis to mixtures of 4-propylguaiacol and either methanol or water. Here, we demonstrate that 4-propylguaiacol is an effective solvent for lignin extraction and depolymerization during RCF, especially when used in combination with water as a cosolvent. Cosolvent mixtures of 4-propylguaiacol and water enable up to 81% lignin extraction, monomer yields up to 25 wt %, and postreaction phase separation. However, unlike methanol, water as a cosolvent fails to inhibit aromatic ring hydrogenation when conducted over Ru/C as a catalyst, potentially leading to excess hydrogen consumption in a process utilizing this approach. Nonetheless, these results suggest a promising strategy for eliminating external organic solvents from RCF by utilizing mixtures of lignin-derived compounds and water as alternative extraction solvents.

09 BIOMASS FUELS↗

NASA NTRS2003

Environmental Snapshots for Satellite Multi-Angle Aerosol Retrieval Validation During the ACE-Asia Field Campaign

On five occasions spanning the ACE-Asia field experiment in spring 2001, the multi-angle imaging MISR instrument, flying aboard the NASA Earth Observing System s Terra satellite, took quarter-kilometer data over a 400-km-wide swath, coincident with high-quality observations by multiple instruments on two or more participating surface and airborne platforms. The cases capture a range of clean, polluted, and dusty aerosol conditions. They represent some of the best opportunities during ACE- Asia for comparative studies among intensive and extensive aerosol observations in their environmental context. We inter-compare related measurements and discuss the implications of apparent discrepancies for each case, at a level of detail appropriate to the analysis of satellite observations. With a three-stage optical modeling process, we synthesize data from multiple sources into layer-by-layer snapshots that summarize what we know about the state of the atmosphere and surface at key locations during each event, to be used for satellite vicarious calibration and aerosol retrieval validation. Aerosols within a few kilometers of the surface were composed primarily of pollution and Asian dust mixtures, as expected. Accumulation and coarse-mode particle size distributions varied little among the events studied, but column aerosol optical depth changed by more than a factor of four, and the near-surface proportion of dust ranged from about 25% to 50%. The amount of absorbing material in the sub-micron fraction was highest when near-surface winds crossed Beijing and the Korean Peninsula, and was considerably lower for all other cases. Ambiguities remain in segregating size distributions by composition; having simultaneous single scattering albedo measurements at more than a single wavelength would significantly reduce the resulting optical model uncertainties, as would integral constraints from surface and atmospheric radiative flux observations. The consistency of component particle micro-physical properties among the five events, even in this relatively complex aerosol environment, suggests that global, satellite-derived maps of aerosol-air-mass-type extent, combined with targeted in situ measurements, can provide a detailed global picture of aerosol behavior. Further joint satellite and in situ analysis is needed to assess the spatial variability of both intensive and extensive aerosol properties within aerosol air masses in two spatial dimensions.

Kahn, Ralph↗

DOE OSTI2024

Analysis of differential scanning calorimetry data for aged plutonium

Differential scanning calorimetry data for samples of a 52 year old plutonium alloy with 3.3 at. % Ga that were heated beyond the melting point is analyzed using transition state theory to find activation energies for the δ to ε and ε to liquid phase transitions. A Bayesian statistical method involving a Gaussian process model is used to find mean values and confidence intervals for the activation energies. The activation energy for the δ to ε phase transition increases by 3.3 ± 3.8% per decade, relative to the case when all age related plutonium lattice point defects have been removed through annealing. The corresponding increase in activation energy for the ε to liquid transition is shown to be 7.1 ± 1.8% per decade. It is postulated that the change in activation energy with age for both phase transitions is caused, in part, by the accumulation of the same type of lattice point defects associated with the observed increase in elastic bulk modulus over time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Decoupling of redox processes from soil saturation in Arctic tundra

Permafrost thaw in warming Arctic landscapes alters hydrology and saturation-driven biogeochemical processes. Models assume that aerobic respiration occurs in drained soils while saturated soils support methanogenesis; however, saturated soils maintain redox gradients that host a range of anaerobic metabolisms. We evaluated how redox potential and redox-active solutes vary with soil moisture in the active layer of permafrost-affected acidic and non-acidic tundra hillslopes. Oxidizing conditions persisted in highly permeable organic horizons of both unsaturated tussock tundra and saturated wet sedge meadows. Redox potential decreased with depth in all soils as increasing soil bulk density restricted groundwater flow and oxygen diffusion. High concentrations of dissolved iron, phosphate, and organic carbon coincided with redox boundaries below the soil surface in acidic tundra, indicating active iron redox cycling and potential release of adsorbed phosphate during iron (oxyhydr)oxide dissolution. In non-acidic tundra, weatherable minerals affected nutrient dynamics more than redox-driven iron cycling, especially in low-lying, saturated areas where thaw reached mineral soils. The role of thaw depth and the ability of saturated soils to maintain oxidizing conditions in organic surface layers highlight the importance of soil physical properties and hydrology in predicting biogeochemical processes and greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Factors influencing mode choice of adults with travel-limiting disability

Introduction: Despite the plethora of research devoted to analyzing the impact of disability on travel behavior, not enough studies have investigated the varying impact of social and environmental factors on the mode choice of people with disabilities that restrict their ability to use transportation modes efficiently. This research gap can be addressed by investigating the factors influencing the mode choice behavior of people with travel-limiting disabilities, which can inform the development of accessible and sustainable transportation systems. Additionally, such studies can provide insights into the social and economic barriers faced by this population group, which can help policymakers to promote social inclusion and equity. Method: This study utilized a Random Parameters Logit model to identify the individual, trip, and environmental factors that influence mode selection among people with travel-limiting disabilities. Here, using the 2017 National Household Travel Survey data for New York State, which included information on respondents with travel-limiting disabilities, the analysis focused on a sample of 8,016 people. In addition, climate data from the National Oceanic and Atmospheric Administration were integrated as additional explanatory variables in the modeling process. Results: The results revealed that people with disabilities may be inclined to travel longer distances walking in the absence of suitable accommodation facilities for other transportation modes. Furthermore, people were less inclined to walk during summer and winter, indicating a need to consider weather conditions as a significant determinant of mode choice. Moreover, low-income people with disabilities were more likely to rely on public transport or walking. Conclusion: Based on this study’s findings, transportation agencies could design infrastructure and plan for future expansions that is more inclusive and accessible, thus catering to the mobility needs of people with travel-limiting disabilities.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2021

DPC Cement Filler Development (Progress Report)

The fillers research and development (R&D) program, mostly experimental, is part of a broader R&D program that includes new process modeling and performance assessment of criticality effects and the overall importance of criticality to repository performance (consequence screening). A literature research and consultation effort with experts by Hardin and Brady (2018) identified several potentially effective and workable filler materials including cements (primarily phosphate based), moltenmetal alloys, and low-temperature glasses. Filler attributes were defined, and the preliminary lists were compared qualitatively. Further comparative analysis will be done (e.g., cost estimates) after experimental screening has narrowed the list of alternatives. The following cement filler compositions were selected for experimental development work and accelerated testing in FY20: Aluminum phosphate cements (APCs); more specifically aluminum oxide / aluminum phosphate (Al 2 O 3 / AlPO 4 ) cements in which Al 2 O 3 serves as the filler material bound by an AlPO 4 binder formed by the reaction of Al 2 O 3 with H 3 PO 4 ; Calcium phosphate cements (CPCs); more specifically composed of pure or nearly pure hydroxyapatite or HAP (Ca 5 (PO 4 ) 3 (OH)); Wollastonite phosphate cements (WPC), specifically wollastonite and aluminum or calcium aluminum phosphates in which CaSiO 3 serves as the filler material and the phosphate serves as the binder. The FY20 effort focused on the optimization of compositions and subsequent processing of these three materials to achieve dense and well-consolidated monolithic samples with 30 to 40% porosity and permeabilities of 1 millidarcy. At the close of this progress report the aluminum phosphate cements (APCs) and the wollastonite phosphate cements (WPCs) appear to show the most promise for continued development. Less progress has been made with the calcium phosphate cements (CPCs); their slurry viscosities are high (and difficult to measure) and they exhibit relatively short cure times of 2 to 3 hours with concomitant and excessive volatile (e.g. CO 2 ) generation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Differentiable modelling to unify machine learning and physical models for geosciences

Process-based modelling offers interpretability and physical consistency in many domains of geosciences but struggles to leverage large datasets efficiently. Machine-learning methods, especially deep networks, have strong predictive skills yet are unable to answer specific scientific questions. Here, in this Perspective, we explore differentiable modelling as a pathway to dissolve the perceived barrier between process-based modelling and machine learning in the geosciences and demonstrate its potential with examples from hydrological modelling. ‘Differentiable’ refers to accurately and efficiently calculating gradients with respect to model variables or parameters, enabling the discovery of high-dimensional unknown relationships. Differentiable modelling involves connecting (flexible amounts of) prior physical knowledge to neural networks, pushing the boundary of physics-informed machine learning. It offers better interpretability, generalizability, and extrapolation capabilities than purely data-driven machine learning, achieving a similar level of accuracy while requiring less training data. Additionally, the performance and efficiency of differentiable models scale well with increasing data volumes. Under data-scarce scenarios, differentiable models have outperformed machine-learning models in producing short-term dynamics and decadal-scale trends owing to the imposed physical constraints. Differentiable modelling approaches are primed to enable geoscientists to ask questions, test hypotheses, and discover unrecognized physical relationships. Future work should address computational challenges, reduce uncertainty, and verify the physical significance of outputs.

58 GEOSCIENCES↗

DOE OSTI2024

The ATLAS experiment software on ARM

With an increased dataset obtained during the Run 3 of the LHC at CERN and the even larger expected increase of the dataset by more than one order of magnitude for the HL-LHC, the ATLAS experiment is reaching the limits of the current data processing model in terms of traditional CPU resources based on x86_64 architectures and an extensive program for software upgrades towards the HL-LHC has been set up. The ARM architecture is becoming a competitive and energy efficient alternative. Some surveys indicate its increased presence in HPCs and commercial clouds, and some WLCG sites have expressed their interest. Chip makers are also developing their next generation solutions on ARM architectures, sometimes combining ARM and GPU processors in the same chip. Consequently it is important that the ATLAS software embraces the change and is able to successfully exploit this architecture. We report on the successful porting to ARM of the Athena software framework, which is used by ATLAS for both online and offline computing operations. Furthermore we report on the successful validation of simulation workflows running on ARM resources. For this we have set up an ATLAS Grid site using ARM compatible middleware and containers on Amazon Web Services (AWS) ARM resources. The ARM version of Athena is fully integrated in the regular software build system and distributed in the same way as other software releases. In addition, the workflows have been integrated into the HEPscore benchmark suite which is the planned WLCG wide replacement of the HepSpec06 benchmark used for Grid site pledges. In the overall porting process we have used resources on AWS, Google Cloud Platform (GCP) and CERN. A performance comparison of different architectures and resources will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects

Researchers at New Mexico Tech, Los Alamos National Laboratory (LANL), Sandia National Laboratory (SNL), and other collaborators are developing a methodology to improve characterization of stress in the subsurface by means of a model-based inversion of reservoir engineering data, time lapse seismic measurements, and microseismicity. Historically, various geophysical techniques have been used in efforts to understand the state of stress in the subsurface through direct and indirect imaging of stress sensitive features (faults and fractures) and observations of transient stress related observations (time variant elastic moduli and microseismicity). Standard direct and indirect techniques for seismic fault and fracture imaging suffer from detectability limits due to reliance on high data quality and multiplicity. Variations in effective elastic properties from inversion of high quality seismic (Vertical Seismic Profile (VSP)) data may be used to infer stress near the wellbore through integration with independent experimental characterization of the stress-velocity relationship. Given a sufficiently robust observation network, microseismic emissions may be inverted to characterize focal mechanisms which, together with supporting assumptions and constraints, inform estimates of in-situ stress. While each of these techniques contributes in part to the characterization of stress within limited spatial and temporal domains, no one method provides an unambiguous stress measurement or a predictive capability over a site scale spatial extent. Challenges associated with solution non-uniqueness, measurement ambiguity, and irregular sampling may be greatly minimized through combination of one or more of independent measurements within a common framework in which realistic geological, hydrodynamic, geomechanical, and seismological constituent process models may act as constraints.

02 PETROLEUM↗