SEARCH · Engineering Papers

Results for “variational fracture modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Analysis of Europan Cycloid Morphology and Implications for Formation Mechanisms

Europa's highly fractured crust has been shown to contain features with a range of differing morphologies. Most lineaments on Europa are believed to have initiated as cracks, although the type of cracking (e.g. tensile vs. shear) remains unclear and may vary for different morphologies. Arcuate lineaments, called cycloids or flexi, have been observed in nearly all imaged regions of Europa and have been modeled as tensile fractures that were initiated in response to diurnal variations in tides. Despite this hypothesis about the formation mechanism, there have been no detailed analyses of the variable morphologies of cycloids. We have examined Galileo images of numerous locations on Europa to develop a catalog of the different morphologies of cycloids. This study focuses on variations in morphology along individual cycloid segments and differences in cusp styles between segments, while illustrating how morphologic evidence can help unravel formation mechanisms. In so doing, we present evidence for cycloid cusps forming due to secondary fracturing during strike-slip sliding on pre-existing cycloid segments.

Marshall, S. T.↗

DOE OSTI2022

3D Simulations of gas injection on callovo-oxfordian claystone assuming spatial heterogeneity and anisotropy

A series of gas injection tests on Callovo-Oxfordian (COx) claystone from the Bure underground research laboratory (URL) in France were carried out at the British Geological Survey (BGS). The tests were performed using a triaxial apparatus specifically designed to capture small volumetric strains induced by the injected gas flow and consequent material dilatancy. The long-duration experiments were monitored throughout. Measurements also included pressure, stresses (axial and radial stresses prescribed for each test stage), rate of gas inflow, gas outflow volume as well as pore-pressures observed at various points of the sample. A coupled hydro-gas-mechanical 3D numerical model has been developed to simulate the tests. Initial permeability is assumed heterogeneous throughout the specimen and embedded fractures are incorporated in the formulation. Gas pressure-induced deformations during the test lead to variations of permeability due to changes in matrix porosity and, especially, fracture aperture as well as fracture orientation due to material anisotropy. A programme of sensitivity analyses involving the variation of different aspects and parameters of the model contributes to a better understanding of the phenomena and highlights its complexity. The model is able to reproduce the observed behaviour of the tests.

58 GEOSCIENCES↗

NASA NTRS1990

Thermal stresses, differential subsidence, and flexure at oceanic fracture zones

Geosat geoid undulations over four Pacific fracture zones have been analyzed. After correcting for the isostatic thermal edge effect, the amplitudes of the residuals are shown to be proportional to the age offset. The shape of the residuals seems to broaden with increasing age. Both geoid anomalies and available ship bathymetry data suggest that slip must sometimes occur on the main fracture zone or secondary faults. Existing models for flexure at fracture zones cannot explain the observed anomalies. A combination model accounting for slip and including flexure from thermal stresses and differential subsidence is presented. This model accounts for lateral variations in flexural rigidity from brittle and ductile yielding due to both thermal and flexural stresses and explains both the amplitudes and the shape of the anomalies along each fracture zone. The best fitting models have mechanical plate thicknesses that are described by the depth to the 600-700 C isotherms.

Wessel, Pal↗

NASA NTRS1990

Geoid anomalies over two South Atlantic fracture zones

Seasat altimetry profiles across the Falkland-Agulhas and Ascension fracture zones were examined for evidence of step-like geoid offsets predicted from thermal modeling of the lithosphere. Neither the plate nor the half-space model of lithospheric thermal evolution succeeds in predicting the variations in step height with age across the fracture zone branches considered. For ages less than about 30 Ma, the step offsets decrease in a manner suggesting a plate with a thickness of 50-75 km. At greater ages, the offsets show complex behavior that may be due to bathymetric features adjacent to the fracture zones. Similar geoid patterns on opposite branches of the Falkland-Agulhas fracture zone are indicative of processes that act symmetrically on both sides of the Mid-Atlantic Ridge. This is consistent both with small-scale convection occurring beneath the lithosphere and with bathymetric features originally produced along the ridge crest and now located symmetrically on opposite sides of the ridge. The west flank of the Ascension fracture zone displays a regrowth in step height at about 40 Ma consistent with small-scale convection and in agreement with other studies of Pacific and South Atlantic fracture zones.

Freedman, Adam P.↗

NASA NTRS1986

Convection beneath young oceanic lithosphere - Implications for thermal structure and gravity

Small-scale convection under the oceanic lithosphere which begins in the first 5 m.y. of cooling can produce a gravity signal with the amplitude and wavelength observed for large areas of the central Pacific and southern Indian oceans using Seasat altimeter data. The trend of the observed anomalies is parallel to the direction of plate motion as might be expected if they were produced by small-scale convection. Models predict that the wavelengths of gravity anomalies increase more rapidly with age than is observed. The persistence of short relatively uniform wavelength anomalies (less than 200 km) to crustal ages of 50 Ma may indicate that they were produced when the lithosphere was very young and thin and were 'frozen in' as cooling thickened the elastic lithosphere. Small-scale convection which begins under very young lithosphere does not violate other geophysical data such as the rate of seafloor subsidence and variations of geoid height with age. After convection has begun, the subsidence due to thermal contraction within the lithosphere varies linearly with age, in the absence of mantle heat sources, although the rate of change of these quantities is affected by convection. Much of the variation of the geoid height across fracture zones can be fit by a model which includes small-scale convection.

Buck, W. R.↗

NASA NTRS2002

Prediction of Long-Term Strength of Thermoplastic Composites Using Time-Temperature Superposition

Accelerated tests for composite failure were investigated. Constant ramp transverse strength tests on thermoplastic composite specimens were conducted at four temperatures from 300 F to 450 F and five duration times from 0.5 sec to 24 hrs. Up to 400 F, the time-temperature-superposition method produces a master curve allowing strength at longer times to be estimated from strength tests conducted over shorter times but at higher temperatures. The shift factors derived from compliance tests applied well to the strength data. To explain why strength behaved similar to compliance, a viscoelastic fracture model was investigated based on the hypothesis that the work of fracture for crack initiation at some critical flaw remains constant with time and temperature. The model, which used compliance as input, was found to fit the strength data only if the critical fracture energy was allowed to vary with stress rate. Fracture tests using double cantilever beam specimens were conducted from 300 F to 450 F over time scales similar to the strength study. The toughness data showed a significant change with loading rate, less variation with temperature, did not form a master curve, and could not be correlated with the fracture model. Since the fracture model did not fit the fracture data, an alternative explanation based on the dilatational strain energy density was proposed. However the usefulness of this model is severely limited because it relies on a critical parameter which varies with loading rate.

Reeder, James R.↗

DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1989

Thermoplastic matrix composites - Finite-element analysis of mode I and mode II failure specimens

A finite-element analysis was conducted to evaluate the stress distributions within mode I and mode II failure specimens, assuming both isotropic and orthotropic elastic material properties. The effects of anisotropy on both the magnitude and the location of the highest stress concentration at the vicinity of the crack tip are significant. The results from modeling realistic blunt crack tip geometry and resin-rich zones imply that local variations in the microstructure strongly influence the stress state near the crack tip and therefore the measured fracture properties. In addition, the features of a viscoelastic model for thermoplastic matrices are described. This model will be used in future investigations of matrix-dominated failure phenomena.

Bankert, Ray J.↗

DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1998

Accelerated Strength Testing of Thermoplastic Composites

Constant ramp strength tests on unidirectional thermoplastic composite specimens oriented in the 90 deg. direction were conducted at constant temperatures ranging from 149 C to 232 C. Ramp rates spanning 5 orders of magnitude were tested so that failures occurred in the range from 0.5 sec. to 24 hrs. (0.5 to 100,000 MPa/sec). Below 204 C, time-temperature superposition held allowing strength at longer times to be estimated from strength tests at shorter times but higher temperatures. The data indicated that a 50% drop in strength might be expected for this material when the test time is increased by 9 orders of magnitude. The shift factors derived from compliance data applied well to the strength results. To explain the link between compliance and strength, a viscoelastic fracture model was investigated. The model, which used compliance as input, was found to fit the strength data only if the critical fracture energy was allowed to vary with temperature reduced stress rate. This variation in the critical parameter severely limits its use in developing a robust time-dependent strength model. The significance of this research is therefore seen as providing both the indication that a more versatile acceleration method for strength can be developed and the evidence that such a method is needed.

Reeder, J. R.↗

DOE OSTI2020

Effects of secondary void-initiating particles on the steady-state crack growth resistance of high-strength steel

In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.

42 ENGINEERING↗

DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗

DOE OSTI2024

DECOVALEX-2023: Task A Final Report

Task A, also known as HGFrac, examines the fracturing processes that may occur in the Callovo-Oxfordian claystone (COx) in the context of the high-level (HLW) and intermediate-level long-lived (ILW-LL) radioactive waste repository in France. Understanding these processes and improving numerical models to reproduce them will aid in the design, optimization, and safety of the repository. Heat and gas fracturing are studied in two independent subtasks following a stepwise approach: laboratory tests/benchmark exercises, in-situ experiment, and finally, an application case. The in-situ heater experiment aimed to thermally induce a hydraulic fracture; temperature and pore pressure were monitored to detect any evidence of fracturing. The in-situ gas injection experiment aimed to study the effect of the stress orientation and gas injection kinetics on the gas fracturing process. The occurrence of fracturing was monitored by gas pressure measured in the injection interval. In both experiments, the excavation-induced fracture network around the heater/injection boreholes played an important role in the reduction of the compressive stress state, leading to both a tensile and shear failure response of the COx. During the first two years of the project, the research teams working on each task developed and/or proposed numerical approaches for reproducing the occurrence of fracturing in the in-situ experiments. The failure criteria were defined by reproducing the measurements from laboratory extension tests for the heat fracturing subtask. In the gas fracturing subtask, their approaches were used for simulating several benchmark exercises, and an inter-comparison between models was carried out. In both tasks, the developed approaches were compared with a simplified approach considering poro-elasticity for the mechanical behaviour of the COx. Most of the developed approaches are based on a continuous medium that takes into account variations in hydraulic properties due to mechanical degradation, such as plastic deformation or damage. Other approaches implicitly modelled weak planes or embedded discontinuities to reproduce fracture propagation. The potential for fracture initiation was also studied through of a discrete approach. In the second half of the project, the research teams mainly focused on interpretative modelling of two in-situ experiments and a blind prediction exercise to test their respective approaches. The models developed by the research teams were also applied at the repository scale to evaluate fracture initiation in a case study under vi unfavourable conditions, particularly in terms of spacing between High-Level Waste cells. The results showed that the poro-elasticity approach could be an efficient tool for understanding the main processes occurring in the COx. One example is the explicit representation of the excavation-induced fracture network around the boreholes, which yielded acceptable results compared to the measurement data. However, advanced approaches were needed to evaluate the potential increase of the excavation-induced fracture network extend and better understand fracture initiation. The stress analyses carried out by the teams revealed that hydraulic boundary conditions had a strong impact on fracture initiation in the heater experiment. Furthermore, in most cases, the results required higher pore pressure increments to reach fracturing than those measured in the experiment. This implies that the measurements may have been biased by the packer’s capacity to fully isolate the piezometric chambers, leading to lower pressures. On the contrary, there was no agreement on the fracturing mode, as some reported either shear or tensile fracturing, while others reported a combination of the two modes. In the case of gas fracturing, the research teams were limited to the comparison of a single point, which complicated their task. Nonetheless, the numerical results were able to reproduce the measurements and capture processes such as longitudinal gas flow through the excavation-induced fracture network, as suggested by some evidence in the observation piezometric chambers. The numerical models also agreed with the measurements in the sense of higher probability of developing along the injection borehole than radially towards the sound rock. The approaches developed by the research teams showed that they are capable of analysing and reproducing fracture initiation in the COx. However, areas of future work should focus on the fracture propagation and fracture aperture, which were out of the scope of this task. To this end, additional data must be gathered for the parameter characterisation and validation of the numerical models. Nonetheless, various approaches showed promising results as they were able to reproduce fracture development under certain conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗

NASA NTRS1976

On a 3-D singularity element for computation of combined mode stress intensities

A special three-dimensional singularity element is developed for the computation of combined modes 1, 2, and 3 stress intensity factors, which vary along an arbitrarily curved crack front in three dimensional linear elastic fracture problems. The finite element method is based on a displacement-hybrid finite element model, based on a modified variational principle of potential energy, with arbitrary element interior displacements, interelement boundary displacements, and element boundary tractions as variables. The special crack-front element used in this analysis contains the square root singularity in strains and stresses, where the stress-intensity factors K(1), K(2), and K(3) are quadratically variable along the crack front and are solved directly along with the unknown nodal displacements.

Atluri, S. N.↗

DOE OSTI2023

Radionuclide transport in fractured chalk under abrupt changes in salinity

Internationally, it has been agreed that geologic repositories for spent fuel and radioactive waste are considered the internationally agreed upon solution for intermediate and long-term disposal. In countries where traditional nuclear waste repository host rocks (e.g., clay, salt, granite) are not available, other low permeability lithologies must be studied. Here, chalk is considered to determine its viability for disposal. Despite chalk's low bulk permeability, it may contain fracture networks that can facilitate radionuclide transport. In arid areas, groundwater salinity may change seasonally due to the mixing between brackish groundwater and fresh meteoric water. Such salinity changes may impact the radionuclides' mobility. In this study, radioactive U(VI) and radionuclide simulant tracers (Sr, Ce and Re) were injected into a naturally fractured chalk core. The mobility of tracers was investigated under abrupt salinity variations. Two solutions were used: a low ionic strength (IS) artificial rainwater (ARW; IS ~0.002) and a high IS artificial groundwater (AGW; IS ~0.2). During the experiments, the tracers were added to ARW, then the carrier was changed to AGW, and vice versa. Ce was mobile only in colloidal form, while Re was transported as a conservative tracer. Both Re and Ce demonstrated no change in mobility due to salinity changes. In contrast, U and Sr showed increased mobility when AGW was introduced and decreased mobility when ARW was introduced into the core. These experimental results, supported by reactive transport modeling, suggest that saline groundwater solutions promote U and Sr release via ion-exchange and enhance their migration in fractured chalk. Here the study emphasizes the impact of salinity variations near spent fuel repositories and their possible impact on radionuclide mobility.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2021

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗