Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PORE STRUCTURE”

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

Synthesis and characteristics of porous ceramic tubes: a comparison of centrifugal casting and cold isostatic pressing methods

Short porous ceramic tubes in small quantities are needed for membrane-related research and development efforts and can be synthesized by centrifugal casting (CC) or cold isostatic pressing (CIP) methods. Here, this work evaluates the characteristics of porous ceramic tubes prepared by both methods and provides a strategy for selecting the proper method for ceramic tube fabrication, depending on its application. Samarium-doped ceria (SDC) with an average crystallite size of 10 nm was synthesized by a citrate method and used to make porous SDC tubes by the CC and CIP methods under different conditions. Porous SDC tubes were prepared from an aqueous suspension of the SDC powder with polyvinyl alcohol (PVA) using a cylindrical metallic mold, centrifugated at 4000–6000 rpm, followed by sintering in 1200–1500 °C. SDC tubes were also made by the CIP method from the SDC powder placed inside a cylindrical rubber bag followed by isostatic compression at 1379 bar. The pore structure of the SDC tubes was characterized by helium permeation and liquid nitrogen Archimedes method. Due to high compression pressure giving a low porosity green tube, the CIP method with lower sintering temperature produces SDC tubes with much lower porosity and average pore size and similar tortuosity compared to the CC methods. The CIP method should be used to prepare low porosity or dense ceramic tubes, and the CC method should be used for preparing high-porosity ceramic tubes or tubes with pore gradient tube wall.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural mechanism of intracellular autoregulation of zinc uptake in ZIP transporters

Abstract Zinc is an essential micronutrient that supports all living organisms through regulating numerous biological processes. However, the mechanism of uptake regulation by intracellular Zn 2+ status remains unclear. Here we report a cryo-electron microscopy structure of a ZIP-family transporter from Bordetella bronchiseptica at 3.05 Å resolution in an inward-facing, inhibited conformation. The transporter forms a homodimer, each protomer containing nine transmembrane helices and three metal ions. Two metal ions form a binuclear pore structure, and the third ion is located at an egress site facing the cytoplasm. The egress site is covered by a loop, and two histidine residues on the loop interact with the egress-site ion and regulate its release. Cell-based Zn 2+ uptake and cell growth viability assays reveal a negative regulation of Zn 2+ uptake through sensing intracellular Zn 2+ status using a built-in sensor. These structural and biochemical analyses provide mechanistic insight into the autoregulation of zinc uptake across membranes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Selenium infiltrated hierarchical hollow carbon spheres display rapid kinetics and extended cycling as lithium metal battery (LMB) cathodes

Lithium metal–selenium (Li–Se) batteries offer high volumetric energy but are limited in their cycling life and fast charge characteristics. Here a facile approach is demonstrated to synthesize hierarchically porous hollow carbon spheres that host Se (Se@HHCS) and allow for state-of-the-art electrochemical performance in a standard carbonate electrolyte (1 M LiPF 6 in 1:1 EC:DEC). The Se@HHCS electrodes display among the most favorable fast charge and cycling behavior reported. For example, they deliver specific capacities of 442 and 357 mA h g -1 after 1500 and 2000 cycles at 5C and 10C, respectively. At 2C, Se@HHCS delivers 558 mA h g -1 after 500 cycles, with cycling coulombic efficiency of 99.9%. Post-mortem microstructural analysis indicates that the structures remain intact during extended cycling. Per GITT analysis, Se@HHCS possesses significantly higher diffusion coefficients in both lithiation and delithiation processes as compared to the baseline. Furthermore, the superior performance of Se@HHCS is directly linked to its macroscopic and nanoscale pore structure: the hollow carbon sphere morphology as well as the remnant open nanoporosity accommodates the 69% volume expansion of the Li to Li 2 Se transformation, with the nanopores also providing a complementary fast ion diffusion path.

25 ENERGY STORAGE↗
DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method (MPM)

Access to fresh drinking water in the future requires effective management of industrial wastewater and water purification from available resources. In this study, we present the simulation methodology and the analysis of a Reverse Osmosis (RO) membrane under various pressure conditions using the Material Point Method (MPM). In contrast to other numerical methods, MPM solves the continuum governing equations on material points in a Lagrangian framework. The method does not require a grid connecting the material points hence making it suitable to simulate large deformations during membrane compaction. The time integration is carried out using the explicit Euler method, while the spatial discretization is performed using linear or cubic-spline shape functions. A series of planar images containing detailed pore structures obtained from X-ray tomography experiments is converted to a three-dimensional collection of material points to simulate the membrane. Compressive loads are applied to the top layer of the membrane to simulate the experiments. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compaction in detail.

ENGINEERING,MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Synthesis and Characterization of Carbon Microbeads

We report a microfluidic-based droplet generation platform for synthesizing micron-sized porous carbon microspheres. The setup employs carbon materials such as graphite, carbon nanotubes, graphene, fullerenes, and carbon black as starting materials. Custom composition, structure, and function are achieved through combinations of carbon materials, cross-linkers, and additives along with variations in process parameters. Carbon materials can be assembled into spheres with a mean diameter of units to hundreds of μm with relatively tight size distribution (<25% RSD). Pore structure and size (tens to hundreds of angstrom) can be modulated by incorporating porogen/coporogen dilutants during synthesis. The microbeads have excellent mechanical stability with an elastic modulus of hundreds of MPa. They can sustain high dynamic fluid flow pressures of up to 9000 psi. This work lays the foundation for synthesizing novel tailorable and customizable carbon microbeads. It opens avenues for applying these novel materials for composite and additive manufacturing, energy, life science, and biomedical applications.

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

The evolution of coal porosity during pyrolysis

Gasification of coal, municipal waste, or other organic materials is a potential hydrogen source that entails complex thermal decomposition and transport processes. This study provides a multiscale analysis of these processes for sub-bituminous (Usibelli, Healy, Alaska) and lignite (Center, North Dakota) coals and provides data useful for process design. The chemistry, mineralogy, and pore structures of pyrolyzed coal and their evolution with thermal decomposition are discussed. Samples pyrolyzed at 200–1000 °C were analyzed by small-angle neutron scattering; ultra-small, small-, and wide-angle X-ray scattering; and other complementary techniques. Scanning electron microscopy showed new pores in the high-temperature-pyrolyzed material. Upon heating, the coals became progressively denser, and the concentration of hydrogen decreased. Changes in pore volume fell into three temperature ranges: an initial, low-temperature range that, for the Usibelli coal, involved an increase in overall porosity; a mid-temperature range associated with pore volume loss; and a high-temperature range associated with significant porosity increase and char formation. This transformation was paralleled by changes in fractal dimension and correlation length. The higher the pyrolysis temperature the greater the small-pore-volume fraction and overall surface area became. Pyrolysis increased the lateral size of coal crystallites, decreased the amorphous fraction, and increased the aromatics fraction and overall coal rank. Comparisons of neutron and X-ray scattering data and subsequent water uptake studies showed that pre-dried coals can re-hydrate relatively rapidly upon exposure to air, which can significantly affect the porosity calculated from small-angle-scattering data. Fits to the cumulative porosity curves provide a method for modeling the physical and chemical transformation of hydrogen-containing feedstock during gasification.

Anovitz, Lawrence {Larry} [ORNL] (ORCID:0000000226↗
DOE OSTI2024

Soil pore characteristics and the fate of new switchgrass-derived carbon in switchgrass and prairie bioenergy cropping systems

Abstract Monoculture switchgrass and restored prairie are promising perennial feedstock sources for bioenergy production on the lands unsuitable for conventional agriculture. Such lands often display contrasting topography that influences soil characteristics and interactions between plant growth and soil C gains. This study aimed at elucidating the influences of topography and plant systems on the fate of C originated from switchgrass plants and on its relationships with soil pore characteristics. For that, switchgrass plants were grown in intact soil cores collected from two contrasting topographies, namely steep slopes and topographical depressions, in the fields in multi-year monoculture switchgrass and restored prairie vegetation. The 13 C pulse labeling allowed tracing the C of switchgrass origin, which X-ray computed micro-tomography enabled in-detail characterization of soil pore structure. In eroded slopes, the differences between the monoculture switchgrass and prairie in terms of total and microbial biomass C were greater than those in topographical depressions. While new switchgrass increased the CO 2 emission in depressions, it did not significantly affect the CO 2 emission in slopes. Pores of 18–90 µm Ø facilitated the accumulation of new C in soil, while > 150 µm Ø pores enhanced the mineralization of the new C. These findings suggest that polyculture prairie located in slopes can be particularly beneficial in facilitating soil C accrual and reduce C losses as CO 2 .

09 BIOMASS FUELS↗
DOE OSTI2021

Characterization of Natural Consolidated Halloysite Nanotube Structures

Halloysite is a unique 1:1 clay mineral frequently appearing with nanotubular morphology, and having surfaces of different polarity with interesting and important technological applications. HNTs can be consolidated naturally in the earth by pressure and thermal flows. In this study of natural consolidated HNTs, the strength and hardness of these materials were found to be dependent on the presence of impurities (gibbsite, alunite, quartz, and other silica minerals), which accounted for the increased stability of such samples. In the absence of impurities, the strength of consolidated HNTs was significantly lower. The first 3D mapping of the pore structure of natural consolidated HNT is provided. The contributions of the porosity within the nanotubes and between the nanotubes were delineated using a combination of non-invasive ultra-small and small-angle X-ray scattering (USAXS/SAXS) analyses, BET/BJH pore size analyses, and computed tomography studies. A total porosity of 40%, as determined by X-ray attenuation and He porosimetry, was found for the natural consolidated HNTs, of which about one-third was due to the inter-HNT porosity. Nano-X-ray computed tomography (nano-XCT) analyses also indicated that 76% of the inter-HNT pores were smaller than 150 nm in diameter. The intra-HNT pore size determined by combined USAXS/SAXS and BET/BJH was about 10 nm. This pore network information is essential for the utilization of natural consolidated HNTs as a model geomaterial to investigate the effects of surface characteristics on confined fluid flow.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interconnected nanoconfining pore networks enhance catalyst CO 2 interaction in electrified reactive capture

Systems that sequentially capture and upgrade CO 2 from air to fuels/fuel-intermediates, such as syngas and ethylene, rely on an energy-intensive CO 2 release process. Electrified reactive capture systems transform CO 2 obtained directly from carbonate capture liquid into products. Previous reactive capture systems show a decline in Faradaic efficiencies (FE) at current densities above 200 mA/cm 2 . Here we show the chemical origins of this problem, finding that prior electrocatalyst designs failed to arrest, activate, and reduce in situ-generated CO 2 (i-CO 2 ) before it traversed the catalyst layer and entered the tailgas stream. We develop a templated synthesis to define pore structures and the sites of Ni single atoms, and find that carbon-nitrogen-based nanopores are effective in accumulating i-CO 2 via short-range, non-electrostatic interactions between CO 2 molecules and the nanochannel walls. These interactions confine and enrich i-CO 2 within the pores, enhancing its binding and activation. We report as a result carbonate electrolysis at 300 mA/cm 2 with FE to CO of 50% ± 3%, and with <1% CO 2 in the tailgas outlet stream. This corresponds to a projected energy efficiency (EE) to 2:1 syngas of 46% at 300 mA/cm 2 when H 2 is added using a water electrolyzer.

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

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2021

Pore characterization of the Marcellus Shale by nitrogen adsorption and prediction of its gas storage capacity

In a shale gas reservoir, pore characterization is an important factor used to determine gas storage capacity. However, the nanometer (nm)-scale pore system in shale is difficult to explore by traditional optical, scanning electron microscopy, or even nuclear magnetic resonance well logging. We have investigated the pore structure and storage capacity of the Marcellus Shale through integration of petrophysical analysis from laboratory and well-logging data and nitrogen adsorption. The isotherm of Marcellus Shale is a composite isotherm, which has features of type I, type II, and type IV isotherms with type H4 of the hysteresis loop, suggesting slit-like pores developed in the Marcellus Shale. Quantitative analysis of pore volumes from the nitrogen adsorption indicates that density porosity may be more properly used to approximate the shale porosity and estimate the shale gas volume. In addition, the specific surface area, micropore, and mesopore volumes have a positive relationship with the kerogen volume and total organic content (TOC). By using the Langmuir and Brunauer-Emmet-Teller models, the simulated result indicates that the higher adsorbed quantity of the Marcellus Shale could be the result of the increase of micropore volume contributed, by the increase of kerogen or TOC content. The proposed equations rapidly compute TOC, a key parameter to predict gas storage capacity in overmature shale such as the Marcellus Shale.

Geochemistry & Geophysics↗
DOE OSTI2021

Modeling and experimental validation of a W f /W-fabrication by chemical vapor deposition and infiltration

Tungsten (W) has a unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall material in future fusion devices. However, its intrinsic brittleness and its further embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) with extrinsic toughening is being developed. A possible synthesis route is chemical vapor deposition (CVD) using heated W fabrics as substrate. The challenge is that the growing CVD-W can isolate domains from precursor access leading to strength-reducing pores. To deepen the process understanding and to optimize the CVD parameters, models were developed with COMSOL Multiphysics and validated experimentally. W deposition rate equations as function of the temperature and the partial pressures of the precursors H 2 and WF 6 were experimentally validated in previous work. In the present article, these equations are applied to obtain partial pressures within the CVD reactor. The results are taken as input for transient simulations in the microscale, in which W coatings, growing onto multiple adjacent W fibers, were simulated via mesh deformation and remeshing. The surface-to-surface contact of the W coatings and the corresponding potential pore formation were simulated by implementing sophisticated deposition rate stop conditions. Within the measuring uncertainties of ≃ ±1%, the models are validated successfully by experimental comparison regarding the deposition rate, pore structure, and relative densities ranging from 0.6 to 0.9.

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

Challenges in Pulsed-Field Gradient Nuclear Magnetic Resonance on Magnetically Heterogeneous Interfaces: Sequence and Field-Dependent Apparent Diffusion Coefficients

It is well known that the internal gradient (gi) that exists within pores haunts the diffusion coefficient (D) as measured by the pulsed-field gradient (PFG) nuclear magnetic resonance (NMR). Several PFG-NMR methods developed to determine an accurate D were not successful. Then, the steady-state diffusion coefficient (Dapp,8) for the cation [C4mim]+ of [C4mim][Tf2N]; [1-butyl-3-methylimidazolium][bis(trifluoromethylsulfonyl)imde] ionic liquid confined in ordered mesoporous carbon (OMC) were determined by comparing Dapp,8 obtained from 1H PFG-NMR performed with three different stimulated echo sequences: STE, APFG, and MPFG under the two external magnetic field strength, B0 = 9.4 and 14.1 Tesla. The measured Dapp,8 which is an order of magnitude smaller than D of bulk [C4mim][Tf2N], is in good agreement between APFG and MPFG both in B0 = 9.4 and 14.1 Tesla. However, the strong gi artifact, which caused apparent diffusion coefficient (Dapp) depending strongly and weakly on B0 and temperature, respectively, in diffusion-time dependent Dapp, Dapp(?) obtained from a sequence with monopolar gradients (STE) was suppressed by using sequences employing bipolar gradients (APFG and MPFG) in the region of steady-state diffusion. But incompletely suppressed gi artifact resulting in the different behaviors of the early part of Dapp(?) between the sequences leads a ˜ 0.6 and 0.9 in MPFG and APFG, respectively, in the relationship between mean squared displacement and diffusion time: = 2Dta, where a = 0.5 and 1 for 1-dimensional single file diffusion and 3-dimensional bulk diffusion, respectively. The above observations clearly show that the diffusion behavior of ions/molecules within the pores and pore structure, such as the surface-to-volume ratio? (D?_app (?)=D_0 [1-4/(9vp) S/V v(D_0 ?)]) and tortuosity (T = D0/Dapp,8), are possible to be misunderstood, especially in the systems with a non-negligible gi. This work demonstrates that it may be necessary to test several PFG sequences under multiple external magnetic fields for the correct determination of the diffusion behavior of ions/molecules in the pores with a larger internal gradient, gi.

Han, Kee Sung↗
DOE OSTI2025

Phase-field modeling of aging-induced microstructure evolution in pentaerythritol tetranitrate thin films and ramifications for shock initiation

Aging of energetic materials may change performance and affect their safety and reliability, but the relationship between microstructure changes induced by aging and consequent performance changes has not been fully established. This work presents results of phase-field method simulations used to model microstructure evolution of vapor-deposited pentaerythritol tetranitrate (PETN) thin films. Simulated aging is shown to induce grain coarsening and substantial changes of the configuration of porosity in the film: Specifically, we show that porosity tends to concentrate in large pores to a greater degree in aged films, a state that is arrived at by closure or consolidation of small pores. To evaluate the performance of the as-deposited and aged films, we perform two-dimensional hydrocode flyer-film impact simulations that incorporate the phase-field output microstructures directly, permitting us to connect features therein to changes in reactivity, a key metric of energy output for shock initiation. The results demonstrate that declining sensitivity obtained for the simulated aged films can be correlated with the loss of fine-structured pores relatively early in the aging process, while long-term microstructure evolution that gradually alters the shape of large, branching pores is less impactful. Finally, we discuss commonalities and discrepancies between our simulation results and high-throughput initiation experiments on shock initiation of aged PETN thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Off-spec fly ash-based lightweight aggregate properties and their influence on the fresh, mechanical, and hydration properties of lightweight concrete: A comparative study

An off-spec fly ash-based spherical lightweight aggregate (LWA), designated as Spherical Porous Reactive Aggregate (SPoRA), was manufactured through a lab pilot-scale production and its engineering properties, including specific gravity, dry rodded unit weight, water absorption, mechanical performance, and pore structure, were evaluated. Using SPoRA, lightweight concrete (LWC) samples were made and their fresh, mechanical, and hydration properties were assessed and compared with LWC samples made using two commercial LWA available in the US market. The results indicated that fine and coarse SPoRA had 72 h water absorption capacities of 16.4 % and 20.9 %, respectively, which were higher than that of the two commercial LWAs. Higher saturated surface dry specific gravity of SPoRA compared to commercial LWAs led to a higher fresh density for the corresponding LWC. Using X-ray computed tomography, large spherical type pores were observed in SPoRA similar to those in the commercial slate-based LWA. Additionally, the pore size distribution of SPoRA, characterized by a dynamic vapor sorption analyzer, indicated that more than 97 % of the pores had diameters greater than 50 nm. SPoRA’s average bulk crushing strength was 6.8 MPa which was smaller than commercial LWA and was potentially attributed to the differences in manufacturing processes. Nonetheless, SPoRA LWC with 28 days compressive strength of 29 MPa passed the ASTM C330 requirement and had a comparable strength with LWC prepared with the commercial LWAs indicating lesser importance for the LWA strength. SPoRA LWC had a flowability of 5 % compared to LWC prepared with commercial LWAs having flowability of 18 %. This observation was potentially related to the fluxing agent used in the SPoRA production. SPoRA had a 70 % degree of hydration at 12 days which was comparable with that of LWC prepared with commercial LWA.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Surface and bulk characterization of reservoir and cap-rocks: Past, present, and future

Geologic formations provide potentially some of the largest volume capacities for CO 2 storage or sequestration. Potential storage sites can be deep saline aquifers, depleted oil reservoirs, and coal seams, surrounded by sealing layers to prevent CO 2 from leaking. It is therefore critical to understand mechanisms contributing to CO 2 trapping and CO 2 leaks. Both phenomena are governed by reactions at the interfaces of the reservoir and cap-rocks and are controlled by the complex chemistry and pore structures of rocks. Mechanisms at the macroscale are affected by the processes occurring at the nanoscale. This review highlights the necessity of multitechnique, multiscale characterization of rocks and points to the importance of surface analysis and surface science studies. Two shale rocks (seals) from Niobrara and Agardhfjellet formations with complex surface chemistry are used as examples throughout the paper. Typically, evaluation of rocks with x-ray diffraction, thermogravimetric analysis, Rock-Eval pyrolysis, gas adsorption, and electron microscopy combined with energy dispersive x-ray spectroscopy is conducted to provide valuable information about the bulk mineralogy, elemental composition, pore volume, and adsorbed species on the sample surface. These studies are necessary prior to designing surface sensitive experiments with x-ray photoelectron spectroscopy (XPS), guiding both sample preparation and sample analysis. XPS has been widely used to study the surface composition of rocks during the investigations of their fine-scale wettability, and the main findings are highlighted here. This paper also reviews the existing literature on ambient-pressure XPS, which provides new opportunities to study in situ chemical alteration due to interactions with CO 2 and offers recommendations for adapting this technique to study rock-fluid interactions, allowing for the identification of fundamental interactions during CO 2 sequestration and guide selection of formation sites for improved storage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multiple Degradation Mechanisms in Reinforced Concrete Structures, Modeling and Risk Analysis

The overarching goal of this project is to complement ongoing Department of Energy (DOE) Light Water Reactor Sustainability (LWRS)-funded Grizzly concrete modeling development efforts by providing improved multi-physics models for incorporation into the Grizzly and BlackBear codes. Objectives of this work identified at the outset of this project include: 1. Coupling between mechanical damage and transport processes. The mechanical damage will be characterized by nonlinear mechanical constitutive models appropriate for concrete used in nuclear power plants (NPPs). 2. Improved representation of the coupling among the transport process, building on the current moisture and thermal transport models in Grizzly. The parameters for coupling among heat and mass transport in the multi-physics model will be experimentally determined. 3. Coupling among different scale levels of concrete constituents, fine and coarse aggregates, cement paste, hydration products, and pore structure. 4. Probabilistic analysis of the random nature of the heterogeneous concrete structure and the environmental factors (ambient temperature and humidity) and their multiple effects on concrete deterioration. 5. Benchmark verification, validation and uncertainty quantification of the thermo-hygro- chemo-mechanical (THCM) concrete formulation. The end goal of this work was to provide a comprehensive and robust simulation capability that can be applied at the engineering scale for analysis of realistic deterioration scenarios in NPP concrete structures.

36 MATERIALS SCIENCE↗