Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependence analysis”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Charge radii of the nucleon from its flavor dependent Dirac form factors

In this report we have determined the proton and the neutron charge radii from a global analysis of the proton and the neutron elastic form factors, after first performing a flavor decomposition of these form factors under charge symmetry in the light cone frame formulation. We then extracted the transverse mean-square radii of the flavor dependent quark distributions. In turn, these are related in a model-independent way to the proton and neutron charge radii but allow us to take into account motion effects of the recoiling nucleon for data at finite but high momentum transfer. In the proton case we find $\langle r_p \rangle$ = 0.852±0.002 (stat.) ± 0.009 (syst.) (fm), consistent with the proton charge radius obtained from muonic hydrogen spectroscopy. The current method improves on the precision of the $\langle r_p \rangle$; extraction based on the form factor measurements. Furthermore, we find no discrepancy in the $\langle r_p \rangle$ determination among the different electron scattering measurements, all of which, utilizing the current method of extraction, result in a value that is consistent with the smallest $\langle r_p \rangle$ extraction from the electron scattering measurements. Concerning the neutron case, past results relied solely on the neutron-electron scattering length measurements, which suffer from an underestimation of underlying systematic uncertainties inherent to the extraction technique. Utilizing the present method we have performed the first extraction of the neutron charge radius based on nucleon form factor data, and we find $\langle r^2_n \rangle$ = -0.122±0.004 (stat.) ± 0.010 (syst.) (fm 2 ).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A theoretical investigation of the effect of Ga alloying on thermodynamic stability, electronic-structure, and oxidation resistance of Ti 2 AlC MAX phase

We present a systematic investigation of thermodynamic stability, phase-reaction, and chemical activity of Al containing disordered Ti 2 (Al-Ga)C MAX phases using machine-learning driven high-throughput framework to understand the oxidation resistance behavior with increasing temperature and exposure to static oxygen. The A-site (at Al) disordering of Ti 2 AlC with Ga shows significant change in the chemical activity of Al with increasing temperature and exposure to static oxygen, which is expected to enable surface segregation of Al, thereby, the formation of Al 2 O 3 and improved oxidation resistance. We performed in-depth convex hull analysis of ternary Ti-Al-C, Ti-Ga-C, and Ti-Al-Ga-C based MAX phase, and provide detailed contribution arising from electronic, chemical and vibrational entropies. The thermodynamic analysis shows change in the Gibbs formation enthalpy (ΔG form ) at higher temperatures, which implies an interplay of temperature-dependent enthalpy and entropic contributions in oxidation resistance Ga doped Ti 2 AlC MAX phases. A detailed electronic structure and chemical bonding analysis using crystal orbital Hamilton population method reveal the origin of change in phases stability and in oxidation resistance in disorder Ti2(Al 1-x Ga x )C MAX phases. Our electronic structure analysis correlate well with the change in oxidation resistance of Ga doped MAX phases. We believe our study provides a useful guideline to understand to role of alloying on electronic, thermodynamic, and oxidation related mechanisms of bulk MAX phases, which can work as a precursor to understand oxidation behavior of twodimensional MAX phases, i.e., MXenes (transition metal carbides, carbonitrides and nitrides).

36 MATERIALS SCIENCE↗
DOE OSTI2024

3D multi-system Bayesian calibration with energy conservation to study rapidity-dependent dynamics of nuclear collisions

Considerable information about the early-stage dynamics of heavy ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems — large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late-time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Data and code for: "Soil respiration response to simulated precipitation change depends on ecosystem type and study duration"

Data, code, and figures from our meta-analysis of the effect of precipitation change on soil respiration, with an emphasis on ecosystem variability and the effect of study duration. We combined the results of 81 separate studies to determine the effect of altered rainfall on the amount of carbon that leaves soil in the form of carbon dioxide, a major component of global carbon cycling. In addition we looked at how long the changes lasted, and how different soil properties and the intensity of precipitation changes affected the study results. We found that more precipitation resulted in greater amounts of carbon dioxide leaving the soil, and less precipitation resulted in less. However, the changes weakened over time in ecosystems that typically receive plenty of rainfall (e.g., forests), in contrast to ecosystems that typically receive little rainfall (e.g., deserts) where changes strengthened over time. Changes in the amount of carbon dioxide that left the soil were also affected by the amount of carbon in the soil, which impacts how much water soil can hold.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Environmental Dependence of Type Ia Supernovae in Low-redshift Galaxy Clusters

Abstract We present an analysis of 102 Type Ia supernovae (SNe Ia) in nearby ( z < 0.1), x-ray-selected galaxy clusters. This is the largest such sample to date and is based on archival data primarily from ZTF and ATLAS. We divide our SNe Ia into an inner cluster sample projected within r 500 of the cluster center and an outer cluster sample projected between r 500 and 2 r 500 . We compare these to field samples of SNe Ia at similar redshifts in both quiescent and star-forming host galaxies. Based on SALT3 fits to the light curves, we find that the inner cluster SNe Ia have a higher fraction of fast-evolving objects (SALT3 x 1 < −1) than the outer cluster or field quiescent samples. This implies an intrinsically different population of SNe Ia occurs in inner cluster environments, beyond known correlations based on host galaxy alone. Our cluster samples show a strongly bimodal x 1 distribution with a fast-evolving component that dominates the inner cluster objects (≳75%) but is just a small fraction of SNe Ia in field star-forming galaxies (≲10%). We do not see strong evidence for variations in the color (SALT3 c ) distributions among the samples and find only minor differences in SN Ia standardization parameters and Hubble residuals. We suggest that the age of the stellar population drives the observed distributions, with the oldest populations nearly exclusively producing fast-evolving SNe Ia.

Larison, Conor (ORCID:0000000320374619)↗
DOE OSTI2026

Nonlinear Information from DESI Luminous Red Galaxies: An Emulator-Based Analysis of Pre- and Post-Reconstruction Power Spectra

We present joint measurements of the pre- and post-reconstruction power spectra, $P_{\rm pre}$ and $P_{\rm post}$, together with their cross-power spectrum, $P_{\rm cross}$, for the Luminous Red Galaxies (LRGs) in the DESI Data Release 1 (DR1). We jointly analyse these observables with an emulator-based full-shape modeling framework, thereby, for the first time, we extract complementary nonlinear information from the galaxy density field before and after reconstruction in real survey data. Specifically, including $P_{\rm post}$ and $P_{\rm cross}$ in addition to $P_{\rm pre}$ (hereafter $P_{\rm all}$) yields an improvement of approximately $18$-$27\%$ in the $σ_8$ constraint in both $Λ$CDM and $w$CDM, depending on the redshift bin, relative to the $P_{\rm pre}$-only analysis with the cosmic microwave background distance priors (hereafter CMB). In $w$CDM, the joint CMB+$P_{\rm all}$ analysis can tighten the constraints on $w$ by approximately $5$-$15\%$ across the two LRG redshift bins, compared to the CMB+$P_{\rm pre}$ combination. Further incorporating the Type Ia supernova dataset and comparing the cosmological constraints in $w$CDM from each individual power-spectrum component with those from the full combination, we find that $P_{\rm all}$ consistently provides the tightest constraints. From the joint CMB+$P_{\rm all}$+DES-Dovekie dataset, we obtain $Ω_m = 0.314 \pm 0.0048$ and $w = -0.988 \pm 0.023$ for the \texttt{LRG1} sample, and $Ω_m = 0.318 \pm 0.0046$ and $w = -0.988 \pm 0.025$ for \texttt{LRG2}. These results demonstrate that combining pre- and post-reconstruction power spectra with their cross-correlation enables DESI to harvest additional nonlinear information, leading to tighter constraints on cosmological parameters.

Wang, Yuting [NAOC, Beijing; Beijing, GUCAS] (ORCI↗
DOE OSTI2024

Skeletal reaction models for methane combustion

A local-sensitivity-analysis technique is employed to generate new skeletal reaction models for methane combustion from the foundational fuel chemistry model (FFCM-1). Here, the sensitivities of the thermo-chemical variables with respect to the reaction rates are computed via the forced-optimally time dependent (f-OTD) methodology. In this methodology, the large sensitivity matrix containing all local sensitivities is modeled as a product of two low-rank time-dependent matrices. The evolution equations of these matrices are derived from the governing equations of the system. The modeled sensitivities are computed for the auto-ignition of methane at atmospheric and high pressures with different sets of initial temperatures, and equivalence ratios. These sensitivities are then analyzed to rank the most important (sensitive) species. A series of skeletal models with different number of species and levels of accuracy in reproducing the FFCM-1 results are suggested. The performances of the generated models are compared against FFCM-1 in predicting the ignition delay, the laminar flame speed, and the flame extinction. The results of this comparative assessment suggest the skeletal models with 24 and more species generate the FFCM-1 results with an excellent accuracy.

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

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Efficient analysis of magnetic field line behavior in toroidal plasmas

The confinement of plasmas in tokamaks and stellarators depends on magnetic field lines lying in nested toroidal surfaces. The transition near the plasma edge away from the lines lying in magnetic surfaces defines properties of divertors. The transition in time defines properties of disruptions. Divertor design and disruption analyses require a detailed understanding of these transitions. The use of a Fourier transform coupled with a Gaussian window function allows far more information to be extracted about these transitions using far shorter field line integrations than can be obtained using traditional methods based on Poincaré plots. The physics of divertors and disruptions is reviewed to clarify why the type of information that can be gained from more efficient methods of analysis is of central importance to the fusion program based on magnetic confinement.

Boozer, Allen H. [Columbia Univ., New York, NY (Un↗
DOE OSTI2025

Geometry and size dependent microstructure and crack formation in Rene 41 superalloy fabricated by laser powder bed fusion

This study provides a systematic investigation into the size and geometry-dependent microstructural evolution and cracking susceptibility of LPBF-fabricated Rene 41. By coupling experimental microstructural analysis with thermal modeling, this research uniquely identifies the relationship between geometry-specific thermal histories, carbide coarsening, and liquation cracking. Four different geometries with varying thickness were fabricated with the identical process parameters. It was found that the grain size and morphology are not affected by the part size. However, the thinner parts exhibited coarser sub-grain structures compared to the thicker ones. The crack formation was observed for the parts with cross-sections smaller than 1 mm, whereas thicker parts had high density without any defects. Further, the cracks were observed in the interdendritic regions, suggesting that liquation cracking was the active micro-crack formation mechanism. The detailed microstructural analysis combined with a thermal finite element analysis proved that the heat extraction efficiency was lower for thinner parts causing a lower cooling rate and coarser carbides, making them more susceptible to constitutional liquation. Microhardness measurements were conducted for each geometry and correlated with the observed microstructural variations. The findings highlight the critical need for geometry-specific optimization of LPBF process parameters to mitigate cracking and achieve microstructural uniformity, offering valuable insights into the fabrication of complex, high-performance aerospace components.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Analysis methodology of coherent oscillations in time- and angle-resolved photoemission spectroscopy

Oscillatory signals from coherently excited phonons are regularly observed in ultrafast pump–probe experiments on condensed matter samples. Electron–phonon coupling implies that coherent phonons also modulate the electronic band structure. These oscillations can be probed with energy and momentum resolution using time- and angle-resolved photoemission spectroscopy (trARPES), which reveals the orbital dependence of the electron–phonon coupling for a specific phonon mode. However, a comprehensive analysis remains challenging when multiple coherent phonon modes couple to multiple electronic bands. Complex spectral line shapes due to strong correlations in quantum materials add to this challenge. In this work, we examine how the frequency domain representation of trARPES data facilitates a quantitative analysis of coherent oscillations of the electronic bands. Here, we investigate the frequency domain representation of the photoemission intensity and the first moment of the energy distribution curves. Both quantities provide complementary information and are able to distinguish oscillations of binding energy, linewidth, and intensity. We analyze a representative trARPES dataset of the transition metal dichalcogenide WTe 2 and construct composite spectra, which intuitively illustrate how much each electronic band is affected by a specific phonon mode. We also show how a linearly chirped probe pulse can generate extrinsic artifacts that are distinct from the intrinsic coherent phonon signal.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Deep learning of structural morphology imaged by scanning X-ray diffraction microscopy

Scanning X-ray nanodiffraction microscopy is a powerful technique for spatially resolving nanoscale structural morphologies by diffraction contrast. One of the critical challenges in experimental nanodiffraction data analysis is posed by the convergence angle of nanoscale focusing optics which creates simultaneous dependency of the far-field scattering data on three independent components of the local strain tensor-corresponding to dilation and two potential rigid body rotations of the unit cell. All three components are in principle resolvable through a spatially mapped sample tilt series; however, traditional data analysis is computationally expensive and prone to artifacts. In this study, we implement NanobeamNN, a convolutional neural network specifically tailored to the analysis of scanning probe X-ray microscopy data. NanobeamNN learns lattice strain and rotation angles from simulated diffraction of a focused X-ray nanobeam by an epitaxial thin film and can directly make reasonable predictions on experimental data without the need for additional fine-tuning. We demonstrate that this approach represents a significant advancement in computational speed over conventional methods, as well as a potential improvement in accuracy over the current standard.

Luo, Aileen [Cornell Univ., Ithaca, NY (United Sta↗
DOE OSTI2023

Market analysis for the integration of new power technologies: A case study of the deployment of hybrid fossil-based generator plus energy storage (ES-FE)

This study examines the national landscape of hybridized fossil energy (FE) power plants with energy storage (ES) technologies (“ES-FE”) and presents the compilation of an ES-FE dataset, which includes over 65 ES-FE projects and concepts in the United States, comprising approximately 500 MWh of co-located ES capacity with FE power plants. This study also estimates the economic feasibility of adding ES to existing FE power plants by characterizing the potential revenues that can be generated by the ES component through flexibility and capacity value. The analysis focuses on ES technologies with 2- to 10-h. durations located in four U.S. independent system operators (ISOs): Midcontinent ISO (MISO), Electric Reliability Council of Texas (ERCOT), PJM Interconnection (PJM), and California ISO (CAISO), which have +70,000 MW of combined FE power capacity that could add ES. Annual revenues are estimated for the ES component using a what-if-analysis approach, for capacity value, price arbitrage, or ancillary services provision. The results show that annual revenues depend on the end-use storage service, wholesale electricity and capacity market prices, and ES technology operation parameters such as discharging duration and cycling frequency. When performing a sensitivity analysis, ES accrues $7–178/kW-yr. via price arbitrage and ancillary services provision in the four ISOs, and $13–92/kW-yr. when providing capacity value only in MISO and PJM. A cash flow analysis is performed to estimate the net present value (NPV) of the ES addition using a range of ES costs. The study finds that for most ES technologies considered, these revenues alone are insufficient to achieve economic feasibility. In conclusion, of the 1645 total runs analyzed, 115 had positive NPVs (7%). Therefore, other revenue streams or monetizable benefits are necessary to achieve the break-even point.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Orientational analysis of atomic pair correlations in nanocrystalline indium oxide thin films

The application of grazing-incidence total X-ray scattering (GITXS) for pair distribution function (PDF) analysis using >50 keV X-rays from synchrotron light sources has created new opportunities for structural characterization of supported thin films with high resolution. Compared with grazing-incidence wide-angle X-ray scattering, which is only useful for highly ordered materials, GITXS/PDFs expand such analysis to largely disordered or nanostructured materials by examining the atomic pair correlations dependent on the direction relative to the surface of the supporting substrate. A characterization of nanocrystalline In 2 O 3 -derived thin films is presented here with in-plane-isotropic and out-of-plane-anisotropic orientational ordering of the atomic structure, each synthesized using different techniques. The atomic orientations of such films are known to vary based on the synthetic conditions. Here, an azimuthal orientational analysis of these films using GITXS with a single incident angle is shown to resolve the markedly different orientations of the atomic structures with respect to the planar support and the different degrees of long-range order, and hence, the terminal surface chemistries. It is anticipated that orientational analysis of GITXS/PDF data will offer opportunities to extend structural analyses of thin films by providing a means to qualitatively determine the major atomic orientation within nanocrystalline and, eventually, non-crystalline films.

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

Strain-Driven Evolution of Structural and Carrier Recombination Dynamics in Germanium

Strain engineering in group-IV semiconductors, in particular germanium (Ge), offers a promising route for advancing next-generation optoelectronic device performance. High minority carrier lifetime serves as a direct indicator of superior material quality and device efficiency. Several factors govern the minority carrier lifetime in both strained Ge and strained GeSn, such as strain level, epilayer thickness, defects, surface roughness, and Sn alloying. In this study, the impact of biaxial strain on effective minority carrier lifetimes of tensile-strained Ge (0% ..epsilon..-Ge/Al?As, 0.78% ..epsilon..-Ge/In0.115?Ga0.885?As, and 0.93% ..epsilon..-Ge/In0.16?Ga0.84?As) and compressively strained GeSn (-0.53% c-Ge0.97?Sn0.03/Al?As) epilayers, grown by molecular beam epitaxy, was experimentally investigated via the microwave-reflection photoconductive decay (..mu..-PCD) technique as a simultaneous function of thickness and strain. The structural and strain-relaxation properties were analyzed by high-resolution x-ray diffraction and Raman spectroscopy, while surface morphology analysis using atomic force microscopy provided insights into the surface roughness of strained Ge-based epilayers. Defect characterization by cross-section transmission electron microscopy revealed superior crystalline quality of ..epsilon..-Ge in one epilayer, while twin boundary formation was observed in the other. We investigated the carrier recombination dynamics in epitaxially strained Ge by extracting the bulk recombination lifetime and surface recombination velocity using the thickness-dependent minority carrier lifetime. High effective minority carrier lifetimes of approximately 60-340 ns were achieved at high strain levels with minimal dependence on the thickness of the strained Ge-based epilayers by means of the ..mu..-PCD technique. Minority carrier lifetimes in Ge-based epilayers exhibit a quadratic dependence on the biaxial strain, ranging from -0.88% (compressive) to +1.6% (tensile). The analysis of the carrier recombination dynamics in strained Ge-based heterostructures offers a pathway to designing group-IV based optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Transverse Λ polarization in e + e – processes within a TMD factorization approach and the polarizing fragmentation function

We perform a re-analysis of Belle data for the transverse Λ and Λ¯ polarization in e + e – annihilation processes within a transverse momentum dependent (TMD) factorization approach. We consider two data sets, one referring to the associated production of Λ’s with a light unpolarized hadron in an almost back-to-back configuration, and one for the inclusive Λ production, with the reconstruction of the thrust axis. We adopt the Collins-Soper-Sterman framework and employ the recent formulation on the factorization of single-inclusive hadron production. This extends a previous phenomenological study carried out in a more simplified TMD approach, leading to a new extraction of the polarizing fragmentation function (FF). While confirming several features of the previous analysis, here we include the proper QCD scale dependence of this TMD FF and carefully exploit the role of its nonperturbative component. Here, the compatibility, within a unique TMD factorization scheme, of double and single-inclusive hadron production is discussed in detail. Moreover, we consider another data set for inclusive Λ production, at much larger energies, from the OPAL Collaboration, in order to test the consistency of the entire approach. Finally, we elaborate on some fundamental issues like the role of SU(2) isospin symmetry and the heavy quark contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

X-ray Emission Spectroscopy of Single Protein Crystals Yields Insights into Heme Enzyme Intermediates

Enzyme reactivity is often enhanced by changes in oxidation state, spin state, and metal-ligand covalency of associated metallocofactors. The development of spectroscopic methods for studying these processes coincidentally with structural rearrangements is essential for elucidating metalloenzyme mechanisms. Herein, we demonstrate the feasibility of collecting X-ray emission spectra of metalloenzyme crystals at a third-generation synchrotron source. In particular, we report the development of a von Hamos spectrometer for the collection of Fe Kβ emission optimized for analysis of dilute biological samples. Here, we further showcase its application in crystals of the immunosuppressive heme-dependent enzyme indoleamine 2,3-dioxygenase. Spectra from protein crystals in different states were compared with relevant reference compounds. Complementary density functional calculations assessing covalency support our spectroscopic analysis and identify active site conformations that correlate to high- and low-spin states. These experiments validate the suitability of an X-ray emission approach for determining spin states of previously uncharacterized metalloenzyme reaction intermediates.

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

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗