Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle fracture”

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 OSTI2024

A high-resolution pseudo-polygon discrete element model for regional sea ice

Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A new non-diffusional gas bubble production route in used nuclear fuel: Implications for fission gas release, cladding corrosion, and next generation fuel design

The relationship between noble metal phase particles and fission gas bubble production in used nuclear fuel is described. The majority of Te atoms within noble metal phase undergo radioactive decay to form stable Xe within a few hours after particle formation. This results in the production of clusters of Xe atoms contained within the solid metal matrix exhibiting an equivalent gas bubble pressure approaching 1 GPa. These high pressure bubbles are stabilized by the UO2 within the bulk of the fuel. However, when these bubbles form near the fuel / cladding interface, in combination with local and temporal damage caused by fission recoil, they are capable of overcoming the fracture strength of the UO2 and rupturing catastrophically. The force of the resulting bubble rupture is sufficient to eject noble metal phase particles several microns into the cladding. This proposed mechanism explains the observance of noble metal phase in cladding and is consistent with a host of morphological features found near the fuel / cladding interface.

Schwantes, Jon M.↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
DOE OSTI2023

An Innovative Insight into Performance Degradation of NCM111 Cathode Induced by Suspension of Operation

The lifespan of lithium-ion batteries varies enormously from fundamental study to practical applications. This big difference has been typically ascribed to the high degree of uncertainty in unpredictable and complicated operation conditions in real-life applications. Here, in this paper, we report that the pause of the charging–discharging process, which is frequently operated in practice but rarely studied in academics, is an important reason for the performance degradation of the NCM111 cathode. It is found that the pause during cycling could trigger a remarkable drop in capacity, giving rise to ~30% more capacity decay compared with the continuously cycled sample. In situ synchrotron X-ray diffraction analysis reveals that the harmful H1–H2 phase transition, which typically appears in the initial cycle but disappears in subsequent cycles, is reactivated by the pausing process. The anisotropic lattice strains that occur during the H1–H2 transition result in mechanical fractures that terminate with an inert NiO-type rock-salt phase on the surface of particles. The present study indicates that the discontinuous usage of rechargeable batteries is also a key factor for cycle life, which might provide a distinct perspective on the performance decay in practical applications.

25 ENERGY STORAGE↗
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↗
DOE OSTI2020

Mobility of Radionuclides in Fractured Carbonate Rocks: Lessons from a Field-Scale Transport Experiment

Current research on radionuclide disposal is mostly conducted in granite, clay, saltstone, or volcanic tuff formations. These rock types are not always available to host a geological repository in every nuclear waste-generating country, but carbonate rocks may serve as a potential alternative. To assess their feasibility, a forced gradient cross-borehole tracer experiment was conducted in a saturated fractured chalk formation. The mobility of stable Sr and Cs (as analogs for their radioactive counterparts), Ce (an actinide analog), Re (a Tc analog), bentonite particles, and fluorescent dye tracers through the flow path was analyzed. The migration of each of these radionuclide analogs (RAs) was shown to be dependent upon their chemical speciation in solution, their interactions with bentonite, and their sorption potential to the chalk rock matrix. The brackish groundwater resulted in flocculation and immobilization of most particulate RAs. Nevertheless, the high permeability of the fracture system allowed for fast overall transport times of all aqueous RAs investigated. This study suggests that the geochemical properties of carbonate rocks may provide suitable conditions for certain types of radionuclide storage (in particular, brackish, high-porosity, and low-permeability chalks). Nevertheless, careful consideration should be given to high-permeability fracture networks that may result in high radionuclide mobility.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Evaluating sealability of blended smart polymer and fiber additive for geothermal drilling with the effect of fracture opening size

Geothermal formations often contain extensive fracture networks. These fracture networks contribute to the significant loss of drilling fluids during geothermal drilling. Multiple loss circulation materials (LCM) such as fiber, granules, and pills have been proposed to tackle this problem but with only limited success. Recent advances in materials science have led to the development of thermoset shape memory polymers (SMP) to address the lost circulation problem. In this paper, we evaluate a thermoset SMP performance in sealing near wellbore fractures of different sizes in geothermal wells. The SMP performance was assessed using granite disks and cylindrical granite cores having fracture sizes of 1000 μm and 3000 μm. A static filtration test was performed using cedar fiber, CaCO 3 , and SMP. Results showed cedar fiber performed better than the CaCO 3 ., reducing fluid loss by 89% and improving sealing pressure by 200 psi. A novel dynamic testing unit that allows for high-temperature testing under flowing conditions was used in this study. The analysis showed that 3% by weight SMP and fiber blends could bridge and plug the 1000 μm fracture. For a larger fracture of 3000 μm width, there was a need to increase the weight concentration of the SMP to 6% to plug the fracture opening effectively. We showed the influence of key parameters such as the type of LCM, concentration, and particle size distribution in optimizing the performance of drilling fluid loss treatment.

02 PETROLEUM↗
DOE OSTI2026

An Anisotropic Yield and Damage Material Model to Improve the Contact Pressure Analysis in a Biomass Shredding System

Size reduction systems used in biomass processing break biomass into smaller pieces by utilizing the kinetic energy from the sharp rotating blades. Abrasive and/or erosive wear caused by biomass comminution results in blade wear of the sharp edged cutters, deteriorating the process efficiency. Here, this study aims to optimize the blade design and improve the system efficiency by attempting to understand the interactions between the blades and biomass particles. Since real-time monitoring of these interactions is impractical during operation, mechanical simulations offer a viable alternative for investigating the shredding process. Yet, the irregular geometry and complex mechanical properties of biomass—such as the anisotropic nature of woodchips and their nonlinear fracture behavior—pose significant challenges for accurately simulating contact pressure. In this work an anisotropic yield material model, along with a damage initiation and evolution function, is applied to the woodchip particle to study the contact pressure on shredder blade, offering a scientific basis for improved blade design and process efficiency. This approach can be extended to other biomass processing systems with similar anisotropic feedstocks, making it a valuable tool for advancing sustainable biomass utilization.

09 - BIOMASS FUELS↗
DOE OSTI2023

The effect of powder feedstock and heat treatment on the thermal and mechanical properties of binder jet printed ZrC

In this study, zirconium carbide (ZrC) disks were fabricated using binder jet printing to study the effect of powder feedstock, print parameters, and heat treatment on flowability and final materials properties. A median volumetric particle size smaller than 10 μm was shown to cause the powder to stop flowing during printing. Disks were printed using ZrC with suitable flowability and then heat-treated at temperatures between 1800°C and 2200°C for 1 or 5h. The density, part shrinkage, thermal diffusivity, and fracture strength all increased with increasing temperature and time. The heat-treated disks were then heated to 2200°C for 5h and the properties converged for disks of the same particle size, indicating the hottest temperature and longest time of exposure dictates the final properties. Lastly, it was shown that larger particles produce lower density materials with worse thermal diffusivity, most likely because of poor connectivity between particles after heat treatment.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of deposition rate on microstructure and mechanical properties of 17-4 PH stainless steel fabricated by laser engineered net shaping

Laser engineered net shaping (LENS) is an additive manufacturing technique for fabricating and repairing metallic components. However, the relationships between its process parameters and the resulting mechanical properties of alloys such as 17-4 precipitation-hardening (PH) stainless steel require further investigation to enable reliable application. This study examines the specific effect of LENS deposition rate on the microstructure and mechanical properties of 17-4 PH stainless steel. Specimens were fabricated at two different deposition rates (8.47 and 9.31 mm s−1), subjected to subsequent solution and H900 aging treatment, and then evaluated via tensile testing, hardness measurements, and microscopy. A higher deposition rate results in a finer grain structure but increased porosity, leading to greater ultimate tensile strength (∼1294 MPa) yet lower ductility (strain at failure ∼5.8%) compared to the slower deposition rate (∼1266 MPa, ∼9.0% strain). Hardness follows the same trend as tensile strength. Tensile fracture surfaces for both conditions exhibited a mixed mode of ductile dimples and brittle quasi-cleavage regions, with chromium/silicon oxide particles identified within dimples. Complementary finite element modeling indicates that small void fractions primarily reduce ductility by enhancing localized plasticity, with only a marginal decrease in tensile strength. These integrated experimental and numerical results elucidate the mechanical properties linked to deposition rate, providing insight into tailoring LENS processes to achieve desired properties of 17-4 PH stainless steel.

17-4 PH stainlesssteel↗
DOE OSTI2024

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗
DOE OSTI2021

Supramolecular dynamic binary complexes with pH and salt-responsive properties for use in unconventional reservoirs

Hydraulic fracturing of unconventional reservoirs has seen a boom in the last century, as a means to fulfill the growing energy demand in the world. The fracturing fluid used in the process plays a substantial role in determining the results. Hence, several research and development efforts have been geared towards developing more sustainable, efficient, and improved fracturing fluids. Herein, we present a dynamic binary complex (DBC) solution, with potential to be useful in the hydraulic fracturing domain. It has a supramolecular structure formed by the self-assembly of low molecular weight viscosifiers (LMWVs) oleic acid and diethylenetriamine into an elongated entangled network under alkaline conditions. With less than 2 wt% constituents dispersed in aqueous solution, a viscous gel that exhibits high viscosities even under shear was formed. Key features include responsiveness to pH and salinity, and a zero-shear viscosity that could be tuned by a factor of ~280 by changing the pH. Furthermore, its viscous properties were more pronounced in the presence of salt. Sand settling tests revealed its potential to hold up sand particles for extended periods of time. In conclusion, this DBC solution system has potential to be utilized as a smart salt-responsive, pH-switchable hydraulic fracturing fluid.

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

pH-SWITCHABLE DYNAMIC BINARY COMPLEXES AS VISCOSITY MODIFYING AGENTS

The study of supramolecular structures has been gaining a lot of momentum in recent years due to their potential application in several interesting areas such as medicine, catalysis, energy, oil and gas, etc. Hydraulic fracturing is one area where the ongoing demand for novel viscosity modifiers has prompted researchers to explore various facets of chemistry. In this study, a dynamic binary complex (DBC) solution is presented, with the potential to be useful in the hydraulic fracturing domain. It has a supramolecular structure that is formed by the self-assembly of oleic acid and diethylenetriamine into an elongated network under alkaline conditions. With less than 2 wt% of constituents in aqueous solution, a viscous gel that exhibits high viscosities even under shear was formed. It was sensitive to pH and salinity, and its zero-shear viscosity could be tuned by a factor of ~280 by changing the pH. Furthermore, its viscous properties were pronounced in the presence of salt. Sand settling tests revealed its potential to hold up sand particles for extended periods of time, longer than 24 hours. In conclusion, this DBC solution system has potential to be utilized as a smart pH-switchable hydraulic fracturing fluid that can be prepared using seawater.

02 PETROLEUM↗
DOE OSTI2022

Novel Proppant Logging Technique for Infill Drilling of Unconventional Shale Wells

Summary During the development of an unconventional play, wells are drilled and completed in batches, and depending on the development plans, current and expected energy market trends, as well as other developmental considerations, new wells are drilled and hydraulically fractured later near existing producing laterals. This creates challenges in terms of optimizing resource recovery and reducing interwell communication. A novel approach is proposed that utilizes systematic composite sampling and analysis of drilling mud returns to look for and quantitatively identify sand particles. The workflow involves cleaning, drying, and segregation of samples into sizes of interest to us (size distribution of pumped proppant in offset parent wells). These samples are imaged at a very high resolution and analyzed for grains using characteristic optical imaging properties to classify proppant sand particles using computer vision algorithms. Further analysis, such as elemental compositional analysis, is used to validate the results from the imaging workflow. We present a case study from the Permian Basin, where a new child well was used as a test case to prove this technology at the Hydraulic Fracturing Test Site (HFTS-2) in Delaware Basin. We introduce new proppant parameters that help identify sustained proppant zones vs. localized propped fractures. We have used additional diagnostics and data collected at the test site to validate observations from the proppant log and have successfully interpreted significantly propped vs. unpropped zones. A key finding from this test has been the significant proppant transport distances observed away from parent wells. Observable proppant was found at a lateral distance of approximately 425 m for one set of parent wells and more than 915 m for another set of parent wells. While a major limitation of this technique is the sampling rate, given adequate sampling, the proposed technology represents a systematic and one-of-a-kind interpretation of spatial proppant distribution while drilling infill wells. It provides us with unique opportunities to better understand the current state of the reservoir being targeted, including zones that are likely highly drained relative to others, and how the planned hydraulic fracturing of child wells can be improved.

Energy & Fuels↗
DOE OSTI2026

The Effects of Realistic Impact-processed Asteroid Topology on Atmospheric Airbursts

When small asteroids (∼10–200 m) enter Earth’s atmosphere they can break up before hitting the ground, releasing a considerable amount of energy in what is known as an atmospheric airburst. To date, numerous studies have examined how an asteroid’s strength, size, shape, and entry characteristics affect the breakup process; but few have assessed the effects of heterogenous internal topologies. We present a two-step high-fidelity simulation approach, which we use to model the atmospheric breakup of asteroids with realistic impact-processed internal topologies and their associated ground effects. In the first step, we model the asteroid as it is impacted to generate an “impact-processed” topology. We then map the resulting fractures (i.e., damage field) from that simulation to the initial state of an atmospheric breakup simulation. For both phases, we use the smoothed particle hydrodynamics code Spheral++. As a demonstrative example, we apply our new simulation pipeline to a Tunguska-scale impactor (50 m diameter, 20 km s −1 entry velocity, and 45° angle) in a case study comparing impact-processed and fully intact asteroids in both planar 2D and full 3D. For sufficiently strong asteroids, we find that impact processing generally causes the asteroid to break up at higher altitudes. We also observe that impact-processed asteroids break up into smaller fragments than objects with equivalent properties that have not been impact-processed.

Astronomy and AstroPhysics↗