Search NASA⌕ Search

SEARCH · Search NASA

Results for “Experimental evolution”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Phase-field modeling of the beta to omega phase transformation in Zr–Nb alloys

A three-dimensional elastoplastic phase-field model is developed, using the Finite Element Method (FEM), for modeling the athermal beta to omega phase transformation in Zr–Nb alloys by including plastic deformation and strain hardening of the material. The microstructure evolution during athermal transformation as well as under different stress states, e.g. uni-axial tensile and compressive, bi-axial tensile and compressive, shear and tri-axial loadings, is studied. The effects of plasticity, stress states and the stress loading direction on the microstructure evolution as well as on the mechanical properties are studied. The input data corresponding to a Zr – 8 at.% Nb alloy are acquired from experimental studies as well as by using the CALPHAD method. Our simulations show that the four different omega variants grow as ellipsoidal shaped particles. Our results show that due to stress relaxation, the athermal phase transformation occurs slightly more readily in the presence of plasticity compared to that in its absence. The evolution of omega phase is different under different stress states, which leads to the differences in the mechanical properties of the material. The variant selection mechanism, i.e. formation of different variants under different stress loading directions, is also nicely captured by our model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Predictive Tools for Customizing Heat Treatment of Additively Manufactured Aerospace Components

Laser-bed powder fusion (LBPF) additive manufacturing is increasingly being used to produce components of complex geometries using the Ni-base superalloy Inconel 718. The composition and the microstructure of the alloy are currently well optimized for wrought components made using conventional manufacturing processes such as rolling, forging, extrusion, etc. The attractive mechanical properties of the alloy result from the underlying austenitic matrix with fine equiaxed grains, and a high density and uniform distribution of the precipitation hardening phase, γ". Heat treatment steps such as homogenization, solutioning and aging are well documented for the wrought alloy. However, when the same wrought alloy compositions are used for the additive manufacturing (AM) processes, the asprocessed microstructure is significantly different, because of the different thermal history associated with LBPF, including rapid solidification and multiple temperature excursions that lead to multiple re-melting and reheating in the solid state. Rapid solidification introduces potential non-equilibrium effects at the moving solid-liquid interfaces that impact the extent of solute segregation, as well as the morphology of the dendritic grains that form. In order to recover the target mechanical properties, AM components have to undergo post-process heat treatments. However, such heat treatments have to be custom designed for the AM process and the component geometry because of the expected vast differences in the microstructure at various locations of a component with complex geometry. The homogenization and precipitation steps should be optimized for the component so that target mechanical properties can be obtained throughout the part. The objective of this research is to utilize High Performance Computing in phase field simulations of microstructure evolution during post-processing of AM components. The physics-based modeling will be beneficial in reducing the experimental effort required for heat treatment process selection, optimization, and certification, thus leading to a significant reduction in energy consumption for AM and post-processing heat treatment. The optimization study will help identify heat treatments steps that are critical for development of a final desired microstructure with the minimum energy input. This combined with shortening of the production cycle (time-to-market) by reducing the number of failed parts (property targets), and reduction in the number of iterations for process optimization, will enable 30-40% savings in the energy costs. Phase field simulations of the degree of homogenization and the effect of local matrix composition on the nucleation and growth of competing precipitating phases were performed using the Microstructure Evolution Using Massively Parallel Phase Field Simulations code developed in-house at the Oak Ridge National Laboratory. The simulations were able to successfully capture the kinetics of nucleation and growth, and morphologies of various precipitating phases as a function of local matrix compositions and composition gradients characteristic of local microstructures arising from location-dependent variations in the thermal conditions. Future work will involve extending the simulations to a length scale consisting of multiple dendrites, so that the effect of homogenization on the coarsening of the dendrites can be simulated and used as an additional input to the optimization of the heat treatment process.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Investigating ceramic breeder pebble beds thermomechanics evolution using novel pressure mapping technology

Ceramic breeder pebble beds undergo complex thermomechanical interactions during blanket operation due to stress build-up and relaxation under the effects of confined thermal expansion, thermal cycling, and creep. Understanding the evolution of such processes can aid in guiding blanket design/assembly, breeder materials developments, predicting performance and possible failure modes. This study introduces experimental techniques that enable us to create an environment in which the pebble bed stresses are organically generated as a result of relevant temperature gradients and magnitudes that allows us to observe the combined thermomechanical interaction effects. A novel non-intrusive in-situ tactile pressure sensing technology is used to generate real-time contact pressure maps revealing the spatial and temporal stress evolution with emphasis on understanding the roles that each of the thermo-mechanical forces play in dictating the pebble bed’s equilibrium conditions. Two types of bed-wall contact pressure drop were recorded: (1) within the subsequent cycles due to pebbles irreversible rearrangements, and (2) within the cycle itself as a result of creep/stress-relaxation. Consequently, by the fourth cycle, the contact pressure drops to negligible values and remains stable. This stability is attributed to two modes of self-regulation: (1) stress self-regulation as a result of pressure rise and fall due to thermal expansion and creep/thermal cycling, respectively, and (2) temperature self-regulation due to the locally enhanced thermal conductivity in the core region of the bed accompanied by deteriorated interface conductance. These two mechanisms are desirable as they lower the probability of the events of pebbles crushing and thermal runaways under high temperatures/stresses and poor heat extraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experimental and Theoretical Studies of Ultrafast Vibrational Energy Transfer Dynamics in Energetic Materials

Energy transfer through anharmonically-coupled vibrations influences the earliest chemical steps in shockwave-induced detonation in energetic materials. A mechanistic description of vibrational energy transfer is therefore necessary to develop predictive models of energetic material behavior. We performed transient broadband infrared spectroscopy on hundreds of femtoseconds to hundreds of picosecond timescales as well as density functional theory and molecular dynamics simulations to investigate the evolution of vibrational energy distribution in thin film samples of pentaerythritol tetranitrate (PETN) , 1,3,5 - trinitroperhydro - 1,3,5 - triazine (RDX) , and 2,4,6 - triamino 1,3,5 - trinitrobenzene (TATB). Experimental results show dynamics on multiple timescales, providing strong evidence for coupled vibrations in these systems, as well as material-dependent evolution on tens to hundreds of picosecond timescales. Theoretical results also reveal pathways and distinct timescales for energy transfer through coupled vibrations in the three investigated materials, providing further insight into the mechanistic underpinnings of energy transfer dynamics in energetic material sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Influence of surface nano-patterning on the placement of InAs quantum dots

We have examined the influence of spontaneous nano-patterning on the placement of InAs quantum dots (QDs) on (Al)GaAs surfaces using an experimental-computational approach. Both atomically flat and mounded surfaces, generated via a surface instability induced by the Ehrlich- Schwoebel barrier, are employed as templates for the subsequent deposition of InAs QDs. Using height profiles from atomic-force micrographs, we simulate QD deposition with a 2D phase field model, which describes the time evolution of the InAs layer driven by a chemical potential gradient. For flat surfaces, phase-field simulations result in QD densities comparable to experimental observations. For mounded surfaces, the simulations reveal QDs preferentially positioned in regions of positive curvature (substrate valleys), e.g., at the edge of surface mounds, consistent with the anisotropic QD placement observed experimentally. We discuss the role of curvaturedriven diffusion in the spontaneous ordering of QDs, demonstrating the applicability of this mechanism to AlGaAs mounds.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In situ study of microstructure evolution and α → ω phase transition in annealed and pre-deformed Zr under hydrostatic loading

The detailed study of the effect of the initial microstructure on its evolution under hydrostatic compression before, during, and after the irreversible α → ω phase transformation and during pressure release in Zr using in situ x-ray diffraction is presented. Two samples were studied: one is plastically pre-deformed Zr with saturated hardness and the other is annealed. Phase transformation α → ω initiates at lower pressure for a pre-deformed sample but for a volume fraction of ω Zr, c > 0.7, a larger volume fraction is observed for the annealed sample. This implies that the proportionality between the athermal resistance to the transformation and the yield strength in the continuum phase transformation theory is invalid; an advanced version of the theory is outlined. Phenomenological plasticity theory under hydrostatic loading is outlined in terms of microstructural parameters, and plastic strain is estimated. During transformation, the first rule is suggested, i.e., the average domain size, microstrain, and dislocation density in ω Zr for c < 0.8 are functions of the volume fraction, c of ω Zr only, which are independent of the plastic strain tensor prior to transformation and pressure. The microstructure is not inherited during phase transformation. Surprisingly, for the annealed sample, the final dislocation density and the average microstrain after pressure release in the ω phase are larger than for the severely pre-deformed sample. The results suggest that an extended experimental basis is required for the predictive models for the combined pressure-induced phase transformations and microstructure evolutions.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Experimental Alteration of Basalt to Support Interpretation of Remote Sensing and In Situ Meausrements from Mars

Major occurrences of hydrous alteration minerals on Mars have been found in Noachian impact craters formed in basaltic targets and detected using visible/near infrared (VNIR) spectroscopy. Until recently phyllosilicates were detected only in craters in the southern hemisphere [1, 2]. However, it has been reported that at least nine craters in the northern plains apparently excavated thick layers of lava and sediment to expose phyllosilicates [3] as well. The MER (Mars Exploration Rovers) rovers previously reported results of in situ measurement indicating the presence of alteration minerals on Mars [4,5] and it was recently reported that the Mars Curiosity rover has detected alteration phases in situ at Yellowknife Bay in Gale crater as well [6,7]. An important discovery for Mars geochronology is that the Chemistry and Mineralogy (CheMin) x-ray diffraction (XRD) instrument on Curiosity detected phyllosilicates indicating that phyllosilicate formation on Mars extended beyond the Noachian Epoch [8]. These discoveries indicate that Mars was globally altered by water in the past but does not constrain formation conditions for alteration phase occurrences, which have important implications for the evolution of the surface and the biological potential on Mars. Understanding the alteration assemblages produced by a range of conditions is vital for the interpretation of phyllosilicate spectral signatures as well as in situ measurements and to decipher the environment and evolution of early Mars. The martian surface has been intensely altered by meteorite impacts whose effects include brecciation and melting of target materials as well as the initiation of hydrothermal circulation in a hydrous target [9,10,11,12]. Impact effects may facilitate aqueous alteration of a basaltic target because the rate of silicate dissolution is a function of the degree of crystallinity, surface area, and temperature. The resultant alteration mineralogies from shocked basaltic target material are a function of the original mineral assemblage in the parent rocks, the chemistry of fluids that interacted with the rocks, and physico-chemical conditions (pH, temperatures, and pressure) during the time of mineral formation. Understanding the alteration assemblages produced by a range of conditions is vital for the interpretation of phyllosilicate spectral signatures and to decipher the environment and evolution of early Mars, and especially for identifying habitable niches in which life could be initiated and sustained. No experimentally controlled and well characterized analog materials that simulate martian shock metamorphism and alteration conditions currently exist for calibrating either remote sensing or in situ measurements of Mars. A series of experiments was initiated to assess the effects of systematic changes in the physico-chemical conditions on Mars analog materials thereby providing samples to ground-truth Mars remote sensing observations from CRISM (Compact Reconnaissance Imaging Spectrometer for Mars) and in situ measurements from Opportunity's Mössbauer and Curiosity's CHEMIN (Chemistry and Mineralogy X-Ray Diffraction/XRay Fluorescence) instruments. Results of initial experimental runs as analysed by SEM-EDS (Secondary Electron Microscopy -Energy Dispersive Spectroscopy) and X-ray Diffraction (XRD) analysis are reported here and lay the foundation for comparison with shocked and altered samples that will be characterized in the next phase of this work.

Bell, M. S.↗
NASA NTRS1974

Trapping of hydrogen atoms in X-irradiated salts at room temperature and the decay kinetics

The salts (hypophosphites, formates, a phosphite, a phosphate, and an oxalate) were X-irradiated, whereby hydrogen formed chemically by a radiolytic process becomes trapped in the solid. By room temperature vacuum extraction, the kinetics for the evolution of this trapped hydrogen was studied mass spectrometrically. All salts except two exhibited second-order kinetics. The two exceptions (NaH2PO2(H2O) and K2HPO4) showed first-order kinetics. Based on experimental results, the escape of hydrogen involves three steps: the diffusion of hydrogen atoms from the bulk to the surface, association of these atoms on the surface (rate controlling step for second-order hydrogen evolution), and the desorption of molecular hydrogen from the surface. The hydrogen does not escape if the irradiated salt is stored in air, apparently because adsorbed air molecules occupy surface sites required in the escape mechanism.

May, C. E.↗
NASA NTRSDate not supplied

Development of a Terrain Mapping/Crater Evolution Measurement using Diffractive Optical Elements

When landing on the moon, understanding the interaction of the engine exhaust plume with the lunar surface is critical for the success of the descent and landing flight phases. Two evaluation tools currently used are computational simulations and ground test measurements. Computational simulations require experimental measurements for comparison/validation, but ground test measurements cannot accurately emulate all aspects of an actual lunar landing; flight tests remain the only method of obtaining fully representative data. A terrain mapping/crater evolution measurement system was developed for potential inclusion on a future lander mission. This system uses two stereo cameras viewing a laser dot grid pattern projected on the ground, where the grid is created by shining a laser through one or two diffractive optical elements. CAD simulations of the stereo imaging system are first used to validate the proposed design. Laboratory testing of the system using both a large-scale fixed-geometry crater and a small-scale evolving-geometry crater validate the use of the system for terrain mapping measurements, and for its potential inclusion on a future lander mission.

Joshua M Weisberger↗
NASA NTRSDate not supplied

Development of a Terrain Mapping/Crater Evolution Measurement using Diffractive Optical Elements

When landing on the moon, understanding the interaction of the engine exhaust plume with the lunar surface is critical for the success of the descent and landing flight phases. Two evaluation tools currently used are computational simulations and ground test measurements. Computational simulations require experimental measurements for comparison/validation, but ground test measurements cannot accurately emulate all aspects of an actual lunar landing; flight tests remain the only method of obtaining fully representative data. A terrain mapping/crater evolution measurement system was developed for potential inclusion on a future lander mission. This system uses two stereo cameras viewing a laser dot grid pattern projected on the ground, where the grid is created by shining a laser through one or two diffractive optical elements. CAD simulations of the stereo imaging system are first used to validate the proposed design. Laboratory testing of the system using both a large-scale fixed-geometry crater and a small-scale evolving-geometry crater validate the use of the system for terrain mapping measurements. High-speed, front-illumination shadow particle tracking of particles ejected from the evolving geometry crater is also performed, demonstrating another diagnostic that can be used to further the understanding of plume-surface interactions.

Joshua M Weisberger↗
DOE OSTI2024

Correlated topological flat bands in rhombohedral graphite

Flat bands and nontrivial topological physics are two important topics of condensed matter physics. With a unique stacking configuration analogous to the Su–Schrieffer–Heeger model, rhombohedral graphite (RG) is a potential candidate for realizing both flat bands and nontrivial topological physics. Here, in this study, we report experimental evidence of topological flat bands (TFBs) on the surface of bulk RG, which are topologically protected by bulk helical Dirac nodal lines via the bulk-boundary correspondence. Moreover, upon in situ electron doping, the surface TFBs show a splitting with exotic doping evolution, with an order-of-magnitude increase in the bandwidth of the lower split band, and pinning of the upper band near the Fermi level. These experimental observations together with Hartree–Fock calculations suggest that correlation effects are important in this system. Our results demonstrate RG as a platform for investigating the rich interplay between nontrivial band topology, correlation effects, and interaction-driven symmetry-broken states.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Following atomistic kinetics on experimental timescales with the kinetic Activation–Relaxation Technique

The properties of materials, even at the atomic level, evolve on macroscopic time scales. Following this evolution through simulation has been a challenge for many years. For lattice-based activated diffusion, kinetic Monte Carlo has turned out to be an almost perfect solution. Various accelerated molecular dynamical schemes, for their part, have allowed the study on long time scale of relatively simple systems. There is still a desire and need, however, for methods able to handle complex materials such as alloys and disordered systems. In this paper, we review the kinetic Activation–Relaxation Technique (k-ART), one of a handful of off-lattice kinetic Monte Carlo methods, with on-the-fly cataloging, that have been proposed in the last few years.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental Verification of Ir 5d Orbital States and Atomic Structures in Highly Active Amorphous Iridium Oxide Catalysts

In iridium oxide catalysts, the electronic states whose energies are in the range of energetics and charge transfer kinetics of the oxygen evolution reaction (OER) originate from the Ir 5d orbital states. However, the understanding of the atomic structures and orbital states underlying catalytic reactivity in amorphous iridium oxide oxygen evolving catalysts (Ir-OECs) is incomplete compared to that of crystalline oxides, owing to a lack of direct experimental verification. Here, we present experimental approaches using resonant inelastic X-ray scattering (RIXS) to directly access Ir 5d orbital excitations at the Ir L 3 edge and atomic pair distribution function (PDF) measurements to characterize electronic and coordination structures at the atomic scale. The so-called iridium blue layer (IrBL) and IrO x were formed from the organometallic precursor complex [Cp*Ir(H 2 O) 3 ]SO 4 and the inorganic precursor IrCl 3 , respectively. Ex situ IrBL and IrO x films for RIXS and PDF measurements were prepared by conditioning electrodeposited films at a low voltage. The incident energy RIXS profile of IrO x exhibited extra weak resonantly enhanced excitation below 2 eV energy loss. The feature was clearly different from a single high-energy excitation above 3 eV of IrBL related to the interband transition between π- and σ-antibonding states. The atomic structure refinement based on PDF measurements revealed the atomic structure domains to have edge- and corner-shared IrO 6 octahedra with trigonal-type distortion. Furthermore, density functional theory calculations guided by the refined atomic structures shed light on the electronic structure corresponding to experimental results, including insulating and metallic phases in ex situ IrBL and IrO x films, respectively. Our study establishes different Ir 5d orbital states and atomic structures in two amorphous Ir oxide OER catalysts in their reduction states.

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

Nontrivial strength of van der Waals epitaxial interaction in soft perovskites

The van der Waals (vdW) bond is traditionally believed to be orders of magnitude lower than the typical chemical bond (i.e., ionic, metallic, or covalent), and hence the effects caused by a vdW interface are thought to be trivial. In this work, by investigating the vdW epitaxial growth of a mechanically soft perovskite on the vdW substrate, we obtained the solid proof of strong non-negligible vdW interaction. The experimental results illustrate the formation of cracks and holes for the relaxation of the vdW strain as well as a lattice-constant-dependent epitaxial angle evolution and a pronounced band-structure change. The first-principles calculations indicate that the contribution of the vdW interaction energy at the epitaxial interface reaches up to more than a quarter of the overall interaction energy far exceeding the traditionally recognized vdW bonding strength. Both the experimental and the theoretical work show that given the appropriate material system, the vdW interaction could be strong enough to significantly manipulate material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Investigation of dynamics of ELM crashes and their mitigation techniques

The accurate prediction of H-mode pedestal dynamics is critical for planning experiments in existing tokamaks and in the design of future tokamaks such as ITER and DEMO. The main objective of the proposed research is to advance the understanding of the physics of H-mode pedestal. Through advances in coupled kinetic-MHD simulations, a new model for H-mode pedestal and ELM crashes as well as an improved model for the bootstrap current will be developed. ELMmitigation techniques will also be investigated. The proposed research will help design efficient confinement scenarios and reduce transient heat loads on the divertor and plasma facing components. During the last two years, the principal investigator (PI) of this proposal actively participated in physics studies related to the DOE Joint Research Targets. These studies include the modeling of divertor heat load in the DIII-D, Alcator C-Mod, and NSTX tokamaks in 2010, and the modeling of H-mode pedestal structure in the DIII-D tokamak in 2011. It is proposed that this close collaboration with experimentalists from major US tokamaks continue during the next funding period. Verification and validation will be a strong component of the proposed research. During the course of the project, advances will be made in the following areas; Dynamics of the H-mode pedestal buildup and recovery after ELM crashes – The effects of neutral fueling, particle and thermal pinches will be explored; Dynamics of ELM crashes in realistic tokamak geometries – Heat loads associated with ELM crashes will be validated against experimental measurements. An improved model for ELM crashes will be developed; ELM mitigation – The effect of resonant magnetic perturbations on ELMs stability and their evolution will be investigated; Development of a new bootstrap current model – A reduced model for will be developed through careful verification of existing models for bootstrap current against first-principle kinetic neoclassical simulations. The primary focus here will be on plasma regimes where the existing reduced models deviate from each other, such as high-beta and low aspect ratio plasmas; and Results will be verified against analytical models and simulation results from other codes. Simulation results will be compared with diagnostic experimental data such as Thomson measurements (HRTC) or flux measurements if available. We have arrangements with experimentalists from DIII-D and JET to use diagnostic edge data for this validation. In summary, we propose to study the H-mode pedestal buildup, the stability and evolution of ELMs using a model that will account for dynamically evolving neoclassical sources.

42 ENGINEERING↗
DOE OSTI2024

DECOVALEX-2023: Task G, SAFENET Final Report

DECOVALEX Task G deals with fracture mechanics at several scales using a combined approach of experimental work, related model development, benchmarking and experimental analysis. Fig. 0.1 provides a graphical abstract for Task G. The experimental basis for the related mechanical (M), hydro-mechanical (HM), and thermo-mechanical (TM) processes comes from the rock mechanics laboratories of Universities of Freiberg (TUBAF) and Edinburgh, as well as the Korea Institute of Civil Engineering and Building Technology (KICT). The experimental work in the rock laboratories is closely linked to the underground research laboratories (URLs) Reiche Zeche (Germany), KURT (Korea), and Mont Terri (Switzerland). The scientific key questions are related to fracture permeability evolution under THM conditions, anisotropy effects on fracturing processes, and thermal fracture slip, which are being addressed in specific steps of the task. A large variety of numerical methods have been developed and applied for experimental analysis, ranging from continuum to discontinuum approaches. A detailed comparison of mechanical and hydro-mechanical processes is given in the section 3 via the benchmarking exercises. As a result of DECOVALEX-2023 Task G we further improved our understanding and predictability of fracturing processes under THM conditions. This was based on robust numerical simulation methods and in-depth experimental analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Interface Coating Design for Dynamic Voltage Stability of Solid‐State Batteries

Abstract Intrinsic or interface thermodynamic voltage windows of solid electrolytes are often narrower than the operational voltage range needed by a full battery, thus various interface decomposition reactions can happen in a practical solid‐state battery. Experimentally, it is found that a proper battery design utilizing the reactions can lead to a dynamic evolution from interface instability to stability, giving the so‐called dynamic voltage stability for advanced battery performance. Here, first the state‐of‐the‐art understanding is articulated about how the dynamic voltage stability should be interpreted in physical picture and treated in computation, emphasizing the potential importance of nonequilibrium reaction pathways. The constrained ensemble computational approach is further applied across most types of solid‐state electrolytes to systematically evaluate and compare their dynamic stability voltage windows in response to the mechanical constriction effect. High‐throughput calculations are used to search for coating materials for different interfaces between sulfide, halide, and oxide electrolytes and typical cathode materials with enhanced dynamic voltage stability. A comparison with experiment is given to highlight the value of these computational predictions.

Wang, Yichao↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗