An official website of the United States government
The formation of desiccation cracks in unsaturated soils as a discontinuity phenomenon can compromise the integrity of civil infrastructure on unsaturated soils. Because of the singularity at such discontinuities, the mathematical modeling of desiccation cracking is challenging. In this study, we apply a coupled nonlocal peridynamic poroelastic framework to model desiccation cracking in unsaturated soils. The soil skeleton is modeled by a nonlocal peridynamic elastic solid. A peridynamic equivalence of the generalized Darcy’s law is utilized to model unsaturated fluid flow. Cracking is determined by a critical stretch criterion between material points as well as an energy criterion. We present numerical simulations of desiccation cracking in soil bars and thin soil discs for one-dimensional cracking and two-dimensional cracking networks, respectively. The numerical results have demonstrated that the proposed nonlocal mathematical framework is a promising and robust method for modeling desiccation cracking in unsaturated soils.
more »
« less
This content will become publicly available on March 1, 2026
Nonlocal Nernst-Planck-Poisson System for Modeling Electrochemical Corrosion in Biodegradable Magnesium Implants
Abstract This paper provides a comprehensive derivation and application of the nonlocal Nernst-Planck-Poisson (NNPP) system for accurate modeling of electrochemical corrosion with a focus on the biodegradation of magnesium-based implant materials under physiological conditions. The NNPP system extends and generalizes the peridynamic bi-material corrosion model by considering the transport of multiple ionic species due to electromigration. As in the peridynamic corrosion model, the NNPP system naturally accounts for moving boundaries due to the electrochemical dissolution of solid metallic materials in a liquid electrolyte as part of the dissolution process. In addition, we use the concept of a diffusive corrosion layer, which serves as an interface for constitutive corrosion modeling and provides an accurate representation of the kinetics with respect to the corrosion system under consideration. Through the NNPP model, we propose a corrosion modeling approach that incorporates diffusion, electromigration and reaction conditions in a single nonlocal framework. The validity of the NNPP-based corrosion model is illustrated by numerical simulations, including a one-dimensional example of pencil electrode corrosion and a three-dimensional simulation of a Mg-10Gd alloy bone implant screw decomposing in simulated body fluid. The numerical simulations correctly reproduce the corrosion patterns in agreement with macroscopic experimental corrosion data. Using numerical models of corrosion based on the NNPP system, a nonlocal approach to corrosion analysis is proposed, which reduces the gap between experimental observations and computational predictions, particularly in the development of biodegradable implant materials.
more »
« less
- Award ID(s):
- 1953346
- PAR ID:
- 10597882
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Journal of Peridynamics and Nonlocal Modeling
- Volume:
- 7
- Issue:
- 1
- ISSN:
- 2522-896X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Summary Geomaterials such as sand and clay are highly heterogeneous multiphase materials. Nonlocality (or a characteristic length scale) in modeling geomaterials based on the continuum theory can be associated with several factors, for instance, the physical interactions of material points within finite distance, the homogenization or smoothing process of material heterogeneity, and the particle or problem size‐dependent mechanical behavior (eg, the thickness of shear bands) of geomaterials. In this article, we formulate a nonlocal elastoplastic constitutive model for geomaterials by adapting a local elastoplastic model for geomaterials at a constant suction through the constitutive correspondence principle of the state‐based peridynamics theory. We numerically implement this nonlocal constitutive model via the classical return‐mapping algorithm of computational plasticity. We first conduct a one‐dimensional compression test of a soil sample at a constant suction through the numerical model with three different values of the nonlocal variable (horizon)δ. We then present a strain localization analysis of a soil sample under the constant suction and plane strain conditions with different nonlocal variables. The numerical results show that the proposed nonlocal model can be used to simulate the inception and propagation of shear banding as well as to capture the thickness of shear bands in geomaterials at a constant suction.more » « less
-
The original two-dimensional bond-based peridynamic (BBPD) framework, which only considers the pairwise forces (compression and tension) between two material points, is extended by incorporating the effect of shear deformation in the calculations and its influence on the failure of the bonds. To this end, each bond is considered as a short Timoshenko beam, and by doing so, the traditional BBPD is enhanced into a more comprehensive model known as multi-polar peridynamic (MPPD). The proposed novel approach explicitly considers the shear influence factor used in Timoshenko beams and introduces a strain-based shear deformation failure criterion. The model is then validated against two benchmark experimental tests (i.e., a standard pure mode I edge crack, and a Kalthoff-Winkler configuration) reported in the literature under in-plane dynamic loading and plane stress conditions. In most cases, the developed model is shown to be more accurate in predicting the crack paths obtained from the experimental results when compared to other theoretical methods delineated in the literature. Furthermore, a noticeable change in crack branching and crack path is observed in a study on the effects of Poisson’s ratio and the loading rate. This investigation also demonstrated that the proposed MPPD model can accommodate materials with Poisson’s ratios up to 1/3, expanding the range beyond the traditional BBPD limitations.more » « less
-
Abstract A multiphysics phase field framework for coupled electrochemical and elastoplastic behaviors is presented, where the evolution of complex solid-electrolyte is described by the variation of the phase field variable with time. The solid-electrolyte interface kinetics nonlinearly depends on the thermodynamic driving force and can be accelerated by mechanical straining according to the film rupture-dissolution mechanism. A number of examples in two- and three- dimensions are demonstrated based on the finite element-based MOOSE framework. The model successfully captures the pit-to-crack transition under simultaneous electrochemical and mechanical effects. The crack initiation and growth has been demonstrated to depend on a variety of materials properties. The coupled corrosion and crystal plasticity framework also predict the crack initiation away from the perpendicular to the loading direction.more » « less
-
Ion transport within saturated porous media is an intricate process in which efficient ion delivery is desired in many engineering problems. However, controlling the behavior of ion transport proves challenging, as ion transport is influenced by a variety of driving mechanisms, which requires a systematic understanding. Herein, we study a coupled advection–diffusion–electromigration system for controlled ion transport within porous media using the scaling analysis. Using the Lattice–Boltzmann–Poisson method, we establish a transport regime classification based on an Advection Diffusion Index (ADI) and a novel Electrodiffusivity Index (EDI) for a two-dimensional (2D) microchannel model under various electric potentials, pressure gradients, and concentration conditions. The resulting transport regimes can be well controlled by changing the applied electric potential, the pressure field, and the injected ions concentration. Furthermore, we conduct numerical simulations in a synthetic 2D porous media and an x-ray microcomputed tomography sandstone image to validate the prevailing transport regime. The simulation results highlight that the defined transport regime observed in our simple micromodel domain is also observed in the synthetic two- and three-dimensional domains, but the boundary between each transport regime differs depending on the variation of the pore size within a given domain. Consequently, the proposed ADI and EDI emerge as dimensionless indicators for controlled ion transport. Overall, our proof-of-concept for ion transport control in porous media is demonstrated under advection–diffusion–electromigration transport, demonstrating the richness of transport regimes that can develop and provide future research directions for subsurface engineering applications.more » « less
