Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “level set”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Impacts of High Variable Renewable Energy Futures on Electric-Sector Decision Making: Demand-Side Effects

Previous work by the Berkeley Lab describes how high shares of variable renewable energy (VRE) such as wind and solar power could change wholesale electricity price dynamics. These include the timing of when electricity is cheap or expensive, locational differences in the cost of electricity, and the degree of regularity or predictability in those costs. Many decentralized decision-makers on the demand-side may not yet have considered the implications of these possible future changes. In this report, we evaluate the sensitivity of a set of demand-side decisions to different levels of VRE penetration ranging from a low of 5-20% to a high of 40-50%. The analysis builds on hourly wholesale energy and capacity prices in different VRE scenarios for four wholesale markets in the United States for the year 2030 (CAISO, ERCOT, NYISO, and SPP). The principal question for this exploration is whether private and public electric-sector decisions that are made based on assumptions reflecting low VRE levels still achieve their intended objective in a high VRE scenario with 40-50% wind and solar? This scoping report evaluates the impacts of changing patterns of peak system needs on the benefits of demand reductions by examining the altered value of different energy efficiency (EE) measures. Similarly, we investigate new opportunities for large energy consumers that may arise from periods with very low wholesale electricity prices. We calculate the value of new process investments (e.g., hydrogen production and other generalized electro-commodities), showcase the varying value of new product storage investments (such as reservoir extensions at a desalination plant), and estimate the benefits of increased process flexibility that uses electricity as a process-input in addition to traditional fossil fuels (e.g., district energy systems). Finally, many decentralized decision-makers and end-use customers are not directly exposed to wholesale electricity prices but instead receive price signals from their retail electricity rates. As wind and solar shares increase, we compare the economic efficiency of flat retail rates relative to more dynamic time-of-use tariffs with and without critical peak-pricing events.

14 SOLAR ENERGY↗
DOE OSTI2022

Background Information on the Protection Requirements in IEEE Std 1547-2018

As the DER penetration level has risen significantly in recent times, parameter settings and configuration of installed DERs are having direct impact on local electric distribution utilities as well as bulk power systems during normal and abnormal grid conditions. On the other hand, all grid-connected DERs in the United States of America must conform to the interconnection requirements prescribed in the IEEE Std 1547-2018. Because of these reasons, it is crucial to understand the expected DER performance during abnormal grid conditions considering the requirements of the IEEE Std 1547-2018. The aim of this document is twofold: (1) to summarize requirements in IEEE Std 1547-2018 having direct implications in distribution network protection, and (2) to provide the reader with knowledge and information that will be required for users to apply the requirements specified. Since the focus of this document is on distribution network protection in the presence of DERs, it is assumed that the readers have basic understanding of distribution network protection as well as working principles of DERs. This document starts with a brief introduction on protection systems for distribution network followed by a discussion on the impact of DER on the protection systems. Current practices to provide network protection in the presence of DER is briefly discussed. The document then explores the DER performance requirements in IEEE Std. 1547-2018; especially on those related to distribution network protection. Section 6.1 in IEEE Std. 1547-2018 provides an overview of capabilities and control requirements for DER under abnormal operating conditions. This section also introduces abnormal operating performance categories I, II, III. The DER response to various types of faults and grid conditions such as short-circuit faults and open phase conditions are discussed in section 6.2 in the standard while requirements for coordination with the Area EPS reclosing scheme are provided in section 6.3. Section 6.4 in the IEEE Std 1547-2018 specifies requirements for mandatory voltage tripping and ride-through requirements during low and high voltage disturbances and section 6.5 specifies similar requirements for low and high frequency disturbances. IEEE Std. 1547 requires that the conformance of the DERs to IEEE Std 1547-2018 requirements should be verified in accordance with IEEE Std. 1547. 1. There are several parameters and settings of DERs that need to be properly selected for reliable operation during abnormal grid condition while applying the requirements of IEEE Std 1547-2018. Key decisions for proper selections of parameter and settings are: 1. Determination of required DER abnormal operating performance category, 2. Determination of DER response (shall trip) to abnormal voltages and 3. Determination of DER response (shall trip) to abnormal frequency.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

The Effect of Luminance Pattern on Nighttime Discomfort Glare Response (CRADA Final Report)

LED luminaires with apertures containing visually resolvable “bright” (i.e., high luminance) spots are commonly used to illuminate outdoor environments. However, the impact of the distance between LEDs, which determines spatial frequency at a given viewing distance, on perceived discomfort caused by these luminaires remains unclear. This study involved 29 participants who were each shown 68 stimuli varying in spatial frequency, intensity, diffusion level, and ambient lighting in a dark laboratory setting. Participants reported their experiences of discomfort from glare, ability to resolve individual LEDs, and afterimages. The findings revealed that increased spatial frequency heightened both discomfort and the likelihood of experiencing afterimages. It is hypothesized that the point spread function of the eye contributed to these effects, where sources began to be perceived as a single larger source with equal or greater intensity. These results suggest that LED luminaire designs should favor configurations with lower spatial frequencies to minimize discomfort. A simple quantity, direct illuminance at the eye, is recommended for use to predict discomfort from glare because it performed like other more complex models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Nondetection of Helium in the Upper Atmospheres of Three Sub-Neptune Exoplanets

We present a search for helium in the upper atmospheres of three sub-Neptune-sized planets to investigate the origins of these ubiquitous objects. The detection of helium for a low-density planet would be strong evidence for the presence of a primary atmosphere accreted from the protoplanetary nebula because large amounts of helium are not expected in the secondary atmospheres of rocky planets. We used Keck+NIRSPEC to obtain high-resolution transit spectroscopy of the planets GJ 1214b, GJ 9827d, and HD 97658b around the 10833 Å He triplet feature. We did not detect helium absorption for any of the planets despite achieving a high level of sensitivity. We used the nondetections to set limits on the planets’ thermosphere temperatures and atmospheric loss rates by comparing grids of 1D models to the data. We also performed coupled interior structure and atmospheric loss calculations, which suggest that the bulk atmospheres (winds) of the planets would be at most modestly enhanced (depleted) in helium relative to their primordial composition. Our lack of detections of the helium triplet for GJ 1214b and GJ 9827d is highly inconsistent with the predictions of models for the present-day mass loss on these planets. Higher signal-to-noise data would be needed to detect the helium feature predicted for HD 97658b. We identify uncertainties in the extreme-ultraviolet fluxes of the host stars and the lack of detailed mass-loss models specifically for cool and metal-enhanced atmospheres as the main limitations to the interpretation of our results. Ultimately, our results suggest that the upper atmospheres of sub-Neptune planets are fundamentally different from those of gas giant planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Impact of High PV Penetration on Voltage Stability

Photovoltaics (PV) are becoming increasingly relevant in modern power systems. With this increase also comes reliability concerns as photovoltaics behave differently than conventional generators. One reliability concern is voltage stability. In this paper, voltage stability of the Austin area in the Electric Reliability Council of Texas (ERCOT) system is studied using dynamic models with varying levels of photovoltaic penetration. The base case is set at 0% renewable penetration. Additional cases include 15% wind penetration and up to 65% photovoltaic penetration. The study results show that voltage/var control capacity is critical to voltage stability, which PV lacks. Voltage regulation of photovoltaics may cause over-voltage and voltage collapse may be more abrupt under high regional photovoltaic penetration.

Till, Jeremy↗
DOE OSTI2025

Universal mass equation for equal-quantum excited-states sets II

We extend our recent study of the Universal Mass Equation for equal-quantum excited-states sets reported by Roper and Strakovsky (Eur Phys J A 61:102, 2025). The masses of twelve baryon sets and sixteen meson sets, with only two equal-quantum excited states in each set, using Breit-Wigner PDG2024 masses and their uncertainties at fixed J P for baryons and J PC for mesons, are fitted by a simple one-parameter logarithmic function, M n = α Ln(n) + M 1 , where n is the level of radial excitation. Two accurate masses that start a set are used to calculate four higher masses in the set accurately. It is noted that α values for $b\bar{b}$ equal-quantum excited-states sets accurately lie on a straight line, whose line parameters can be used to calculate α and predict higher mass states for $b\bar{b}$ sets that have only one known member.

Roper, L. David [Virginia Polytechnic Inst. and St↗
DOE OSTI2021

A causal data fusion method for the general exposure and outcome

Abstract With the advent of the big data era, the need to combine multiple individual data sets to draw causal effects arises naturally in many medical and biological applications. Especially each data set cannot measure enough confounders to infer the causal effect of an exposure on an outcome. In this article, we extend the method proposed by a previous study to causal data fusion of more than two data sets without external validation and to a more general (continuous or discrete) exposure and outcome. Theoretically, we obtain the condition for identifiability of exposure effects using multiple individual data sources for the continuous or discrete exposure and outcome. The simulation results show that our proposed causal data fusion method has unbiased causal effect estimate and higher precision than traditional regression, meta‐analysis and statistical matching methods. We further apply our method to study the causal effect of BMI on glucose level in individuals with diabetes by combining two data sets. Our method is essential for causal data fusion and provides important insights into the ongoing discourse on the empirical analysis of merging multiple individual data sources.

Li, Hongkai↗
DOE OSTI2024

Power Profile Monitoring and Tracking Evolution of System-Wide HPC Workloads

The power & energy demands of HPC machines have grown significantly. Modern exascale HPC systems require tens of megawatts of combined power for computing resources and cooling facilities at full capacity. The current energy trend is not sustainable for future HPC systems, and there is a need to work toward the energy efficiency aspect of HPC performance. Energy awareness of the HPC applications at the job level is essential for running an efficient HPC system. This work aims to develop a pipeline to provide a production-level system-wide overview of the HPC workloads' power profile while handling evolving workloads exhibiting new power trends. We developed an open-set classification model for HPC jobs based on the properties of power profiles to continuously provide a system-wide holistic view of recently completed jobs. The pipeline helps continuously monitor the job-level power usage pattern of HPC and enables us to capture the new trends in applications' power behavior. We employed a comprehensive set of techniques to generate job-level data, custom-designed feature extraction methods to extract critical features from jobs' power profiles, clustering techniques powered by generative modeling, and open-set classification for identifying job profiles into known classes or an unknown set. With extensive evaluations, we demonstrate the effectiveness of each component in our pipeline. We provide an analysis of the resulting clusters that characterize the power profile landscape of the Summit supercomputer from more than 60K jobs executed in a year. The open-set classification classifies the known data sets into known classes with high accuracy and identifies unknown data noints with over 85% accuracy.

Karimi, Ahmad Maroof↗
DOE OSTI2023

Relativistic coupled-cluster calculations of RaOH pertinent to spectroscopic detection and laser cooling

Here a relativistic coupled-cluster study of the low-lying electronic states in the radium monohydroxide molecule (RaOH), a radioactive polyatomic molecule of interest to laser cooling and to the search of new physics beyond the Standard Model, is reported. The level positions of the A 2 Π 1/2 and C 2 Σ states have been computed with an accuracy of around 200 cm –1 to facilitate spectroscopic observation of RaOH using laser induced fluorescence spectroscopy, thereby exploiting the systematic convergence of electron-correlation and basis-set effects in relativistic coupled-cluster calculations. The energy level for the B 2 Δ 3/2 state has also been calculated accurately to conclude that the B 2 Δ 3/2 state lies above the A 2 Π 1/2 state. This confirms X 2 Σ ↔ A 2 Π 1/2 as a promising optical cycling transition for laser cooling RaOH.

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

A Joint Roman Space Telescope and Rubin Observatory Synthetic Wide-Field Imaging Survey

We present and validate 20 deg(2) of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and 5 yr of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg(2) of LSST simulated imaging produced as part of Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2). Both synthetic surveys observe the same simulated DESC DC2 universe. For the synthetic Roman survey, we simulate for the first time fully chromatic images along with the detailed physics of the Sensor Chip Assemblies derived from lab measurements using the flight detectors. The simulated imaging and resulting pixel-level measurements of photometric properties of objects span a wavelength range of similar to 0.3-2.0 mu m. We also describe updates to the Roman simulation pipeline, changes in how astrophysical objects are simulated relative to the original DC2 simulations, and the resulting simulated Roman data products. We use these simulations to explore the relative fraction of unrecognized blends in LSST images, finding that 20-30 per cent of objects identified in LSST images with i-band magnitudes brighter than 25 can be identified as multiple objects in Roman images. These simulations provide a unique testing ground for the development and validation of joint pixel-level analysis techniques of ground- and space-based imaging data sets in the second half of the 2020s - in particular the case of joint Roman-LSST analyses.

gravitational lensing↗
DOE OSTI2025

Day-Ahead Probabilistic Forecasting of Net-Load and Demand Response Potentials with High Penetration of Behind-the-Meter Solar-plus-Storage

The goal of this project is to develop advanced methods for day-ahead net-load forecasting, by leveraging the state-of-the-art machine learning techniques. The developed models produce both point and probabilistic forecasts for a variety of use cases, and are versatile to work with different types of data sets. The innovation lies in the novel design of the architectures, leveraging the most recent advances in machine learning that have not been explored in power systems, accompanied by techniques in the broader artificial intelligence fields such as fuzzy systems. This project has achieved the following accomplishments: (1) preprocessing of over 10 data sets covering varying geographical regions, time horizons, and system levels, which form a robust foundation for training and evaluating forecasting models across a wide range of realistic grid scenarios; (2) development of an interactive web app that enables exploratory analysis of load and generation data, and supports better understanding of data trends, anomalies, and correlations, facilitating model development and stakeholder engagement; (3) implementation of over 10 benchmark models for point and probabilistic forecasting, which include a mix of conventional machine learning methods and state-of-the-art deep learning approaches, providing a comprehensive baseline for performance comparison and validation of the proposed models; (4) development of a fuzzy system based gradient boosting model, tailored for small (less than 3 years) data sets, which achieves a mean absolute percentage error (MAPE) of 4% for point forecasting and a 20% improvement in average pinball loss for probabilistic forecasting; (5) development of a Transformer (a state-of-the-art deep learning architecture) based neural network model, tailored for large (3 years or more) data sets, which achieves a MAPE of 2% for point forecasting and a 20% improvement in average pinball loss for probabilistic forecasting; (6) development of a methodology for quantifying DR potential, and extensions of the previous models for multi-target forecasting of net load and DR potential, which achieve a MAPE of 10% for DR potential.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Lifting of level-1 states in the D1D5 CFT

The D1D5 CFT has a large set of states that are supersymmetric at the ‘free’ orbifold point in moduli space. When we perturb away from this point, some of these states join into long multiplets and lift in energy, while others remain supersymmetric. The count of unlifted states can be bounded below by an index, but the index does not yield the pattern of lifting; i.e., which states join into a long multiplet and how much this multiplet lifts. In this paper we consider the simple case of the D1D5 CFT where the orbifold CFT is a sigma model with targets space (T 4 ) 2 /S 2 and consider states at energy level 1. There are 2688 states at this level. The lifted states form a triplet of long multiplets, and we compute their lift at second order in perturbation theory. Half the members of the long multiplet are in the untwisted sector and half are in the twisted sector. This and other similar studies should help in the understanding of fuzzball states that describe extremal holes, since CFT sectors with low twist describe shallow throats in the dual gravity solution while sectors with high twist describe deep throats.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Molecular Mechanisms of Nanomaterial-Bacterial Interactions Revealed by Omics—The Role of Nanomaterial Effect Level

Nanotechnology is employed across a wide range of antibacterial applications in clinical settings, food, pharmaceutical and textile industries, water treatment and consumer goods. Depending on type and concentration, engineered nanomaterials (ENMs) can also benefit bacteria in myriad contexts including within the human body, in biotechnology, environmental bioremediation, wastewater treatment, and agriculture. However, to realize the full potential of nanotechnology across broad applications, it is necessary to understand conditions and mechanisms of detrimental or beneficial effects of ENMs to bacteria. To study ENM effects, bacterial population growth or viability are commonly assessed. However, such endpoints alone may be insufficiently sensitive to fully probe ENM effects on bacterial physiology. To reveal more thoroughly how bacteria respond to ENMs, molecular-level omics methods such as transcriptomics, proteomics, and metabolomics are required. Because omics methods are increasingly utilized, a body of literature exists from which to synthesize state-of-the-art knowledge. Here we review relevant literature regarding ENM impacts on bacterial cellular pathways obtained by transcriptomic, proteomic, and metabolomic analyses across three growth and viability effect levels: inhibitory, sub-inhibitory or stimulatory. As indicated by our analysis, a wider range of pathways are affected in bacteria at sub-inhibitory vs. inhibitory ENM effect levels, underscoring the importance of ENM exposure concentration in elucidating ENM mechanisms of action and interpreting omics results. In addition, challenges and future research directions of applying omics approaches in studying bacterial-ENM interactions are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

SIW21-95: Hybridizing Synchronous Condensers with Grid-Forming Battery Energy Storage Systems

One of the main challenges associated with deployment of high shares of inverter-based resources (IBRs) in power grid is not only reduced system inertia but also degrading system strength that may cause severe stability impacts. A minimum level of system strength is needed for the power system to remain stable under normal conditions and to return to a steady state condition following a system disturbance. Significant system strength reduction is expected in almost all planned areas for solar and wind generation deployment. Synchronous condensers (SC) have been considered as one main technology to address the system strength issues for the areas with high levels of IBRs. SCs can help improving reliability and resiliency of power system but they do not provide the full range of services needed by power systems for reliable and economic integration of high shares of inverter-coupled variable generation such as PV generation. Services related to active power controls cannot be provided by SCs due to lack of prime mover. Even for provision of reactive power SCs have certain constraints based on their thermal and stability limits. NREL has been conducting research on a hybridized concept that combines SCs with grid forming (GFM) battery energy storage systems (BESS). This super flexible AC transmission system (SuperFACTS) that combines these two technologies in a single plant under the same controller offers a unique scalable set of services to the power system at all levels (transmission, sub-transmission, distribution, islands and isolated microgrids). Depending on use cases, SuperFACTS can be controlled to provide fully dispatchable and flexible operation using energy storage component, provide a full range of existing and future ancillary and reliability services to the grid (similar or better than conventional sources), maintain adequate levels of grid strength and inertia, and provide fault current for proper operation of protection systems. Therefore, this economic and easy to commercialize solution has potential for significant impacts on certain segments of global energy sector. All types of gird services (market based, reliability and resiliency) can be provided by SuperFACTS plants. GFM BESS can act as a self-black start source for each SuperFACTS module, which in turn can act a black start resource for co-located PV and wind power plants, segments of transmission and distribution networks, for conventional power plants, etc. The issue of in-rush currents during black start is addressed by overcurrent capability of SCs. In this paper we describe the results of modeling for SuperFACTS concept.

grid forming↗
DOE OSTI2025

Effect on compression of lowering the design adiabat in the SQ-n campaign

A set of experiments in the SQ-n campaign attempts to increase the compression in layered implosions by reducing the design adiabat from 3.0 down to ∼ 2.6. The experiments utilize nanocrystalline high density carbon ablators driven by a ramped laser pulse drive. Contrary to expectations, the lower design adiabat layered implosion reaches the same level of compression, within the measurement uncertainty, as the previous set of SQ-n layered implosions at the higher design adiabat. Here, we examine potential explanations for this result, including incomplete melting or refreezing of the ablator after the first shock passes into the fuel layer, reproducibility, and dopant areal density. We also discuss additional potential paths to reaching lower adiabat including increased dopant areal density, ramped dopant profiles, and increased initial capsule size which are also based on experimental campaign results from the National Ignition Facility.

Physics - Plasma physics↗
DOE OSTI2026

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

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

Lifting at higher levels in the D1D 5 CFT

The D1D 5 P system has a large set of BPS states at its orbifold point. Perturbing away from this ‘free’ point leads to some states joining up into long supermultiplets and lifting, while other states remain BPS. We consider the simplest orbifold which exhibits this lift: that with N = 2 copies of the free c = 6 CFT. We write down the number of lifted and unlifted states implied by the index at all levels upto 6. We work to second order in the perturbation strength λ. For levels upto 4, we find the wavefunctions of the lifted states, their supermultiplet structure and the value of the lift. All states that are allowed to lift by the index are in fact lifted at order O(λ 2 ). We observe that the unlifted states in the untwisted sector have an antisymmetry between the copies in the right moving Ramond ground state sector, and extend this observation to find classes of states for arbitrary N that will remain unlifted to O(λ 2 )

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗