skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Transmissibility in Interactive Nanocomposite Diffusion: The Nonlinear Double-Diffusion Model
Model analogies and exchange of ideas between physics or chemistry with biology or epidemiology have often involved intersectoral mapping of techniques. Material mechanics has benefited hugely from such interpolations from mathematical physics where dislocation patterning of plastically deformed metals and mass transport in nanocomposite materials with high diffusivity paths such as dislocation and grain boundaries, have been traditionally analyzed using the paradigmatic Walgraef-Aifantis (W-A) double-diffusivity (D-D) model. A long standing challenge in these studies has been the inherent nonlinear correlation between the diffusivity paths, making it extremely difficult to analyze their interdependence. Here, we present a novel method of approximating a closed form solution of the ensemble averaged density profiles and correlation statistics of coupled dynamical systems, drawing from a technique used in mathematical biology to calculate a quantity called the basic reproduction number R0, which is the average number of secondary infections generated from every infected. We show that the R0 formulation can be used to calculate the correlation between diffusivity paths, agreeing closely with the exact numerical solution of the D-D model. The method can be generically implemented to analyze other reaction-diffusion models.  more » « less
Award ID(s):
2015317
PAR ID:
10517579
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Frontiers
Date Published:
Journal Name:
Frontiers in Applied Mathematics and Statistics
Volume:
8
ISSN:
2297-4687
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A novel dislocation-density-based crystal plasticity model for nanocrystalline face-centered cubic metals is developed based on the thermally-activated mechanism of dislocations depinning from grain boundaries. Dislocations nucleated from grain boundary dislocation sources are assumed to be the primary carriers of plasticity in the nanocrystals. The evolution of the dislocation density thereby involves a competition between the nucleation of dislocations from grain boundary defect structures, such as ledges, and the absorption of dislocations into the grain boundary via diffusion processes. This model facilitates the simulation of plastic deformation in nanocrystalline metals, with consideration of the initial microstructure resulting from a particular processing method, to be computed as a direct result of dislocation-mediated plasticity only. The exclusion of grain boundary-mediated plasticity mechanisms in the formulation of the crystal plasticity model allows for the exploration of the fundamental role dislocations play in nanocrystalline plasticity. The combined effect of average grain size, grain size distribution shape, and initial dislocation density on the mechanical performance and strain-rate sensitivity are explored with the model. Further, the influence of the grain boundary diffusivity on post-yielding strain-hardening behavior is investigated to discern the impact that the choice of processing route has on the resulting deformation response of the material. 
    more » « less
  2. In this paper, we present concurrent atomistic-continuum (CAC) simulations of the hydrogen (H) diffusion along a grain boundary (GB), nearby which a large population of dislocations are piled up, in a plastically deformed bi-crystalline bcc iron sample. With the microscale dislocation slip and the atomic structure evolution at the GB being simultaneously retained, our main findings are: (i) the accumulation of tens of dislocations near the H-charged GB can induce a local internal stress as high as 3 GPa; (ii) the more dislocations piled up at the GB, the slower the H diffusion ahead of the slip–GB intersection; and (iii) H atoms diffuse fast behind the pileup tip, get trapped within the GB, and diffuse slowly ahead of the pileup tip. The CAC simulation-predicted local H diffusivity, Dpileup−tip, and local stresses, σ, are correlated with each other. We then consolidate such correlations into a mechanics model by considering the dislocation pileup as an Eshelby inclusion. These findings will provide researchers with opportunities to: (a) characterize the interplay between plasticity, H diffusion, and crack initiation underlying H-induced cracking (HIC); (b) develop mechanism-based constitutive rules to be used in diffusion–plasticity coupling models for understanding the interplay between mechanical and mass transport in materials at the continuum level; and (c) connect the atomistic deformation physics of polycrystalline materials with their performance in aqueous environments, which is currently difficult to achieve in experiments. 
    more » « less
  3. The plasma membranes of cells are thin viscous sheets in which some transmembrane proteins have two-dimensional mobility and some are immobilized. Previous studies have shown that immobile proteins retard the short-time diffusivity of mobile particles through hydrodynamic interactions and that steric effects of immobile proteins reduce the long-time diffusivity in a model that neglects hydrodynamic interactions. We present a rigorous derivation of the long-time diffusivity of a single mobile protein interacting hydrodynamically and thermodynamically with an array of immobile proteins subject to periodic boundary conditions. This method is based on a finite element method (FEM) solution of the probability density of the mobile protein diffusing with a position-dependent mobility determined through a multipole solution of Stokes equations. The simulated long-time diffusivity in square arrays decreases as the spacing in the array approaches the particle size in a manner consistent with a lubrication analysis. In random arrays, steric effects lead to a percolation threshold volume fraction above which long-time diffusion is arrested. The FEM/multipole approach is used to compute the long-time diffusivity far away from this threshold. An approximate analysis of mobile protein diffusion through a network of pores connected by bonds with resistances determined by the FEM/multipole calculations is then used to explore higher immobile area fractions and to evaluate the finite simulation cell size scaling behaviour of diffusion near the percolation threshold. Surprisingly, the ratio of the long-time diffusivity to the spatially averaged short-time diffusivity in these two-dimensional fixed arrays is higher in the presence of hydrodynamic interactions than in their absence. Finally, the implications of this work are discussed, including the possibility of using the methods developed here to investigate more complex diffusive phenomena observed in cell membranes. 
    more » « less
  4. Abstract. The multiple-diffusion domain (MDD) model empirically describes the diffusive behavior of noble gases in some terrestrial materials and has been commonly used to interpret 40Ar/39Ar stepwise degassing observations in K-feldspar. When applied in this manner, the MDD model can be used to test crustal exhumation scenarios by identifying the permissible thermal paths a rock sample could have undergone over geologic time, assuming the diffusive properties of Ar within the mineral are accurately understood. More generally, the MDD model provides a framework for quantifying the temperature-dependent diffusivity of noble gases in minerals. However, constraining MDD parameters that successfully predict the results of step-heating diffusion experiments is a complex task, and the assumptions made by existing numerical methods used to quantify model parameters can bias the absolute temperatures permitted by thermal modeling. For example, the most commonly used method assumes that no domains lose more than 60 % of their gas during early heating steps (Lovera et al., 1997). This assumption is unverifiable, and we show that the Lovera et al. (1997) procedure may bias predicted temperatures towards lower values when it is violated. To address this potential bias and to provide greater accessibility to the MDD model, we present a new open-source method for constraining MDD parameters from stepwise degassing experimental results, called the “MDD Tool Kit” (https://github.com/dgorin1/mddtoolkit, last access: 11 October 2024). This software optimizes all MDD parameters simultaneously and removes any need for user-defined Ea or regression fitting choices used by other tools. In doing so, this new method eliminates assumptions about the domain size distribution. To test the validity of our thermal predictions, we then use the MDD Tool Kit (https://github.com/dgorin1/mddtoolkit) to interpret 40Ar/39Ar results from the Grayback Fault, AZ, USA. Although the resulting thermal histories are consistently ∼ 60–75 °C higher than those found in previous studies, they agree with independent observations from apatite fission track, zircon fission track, and (U-Th)/He. 
    more » « less
  5. ABSTRACT Accounting for uncertainty in model selection is crucial for statistical inference and data‐driven decision‐making, particularly with high‐dimensional data. While multiple studies have focused on constructing model confidence sets, a practical and informative visualization tool to assist in decision‐making under such uncertainty has been lacking. This paper introduces an intuitive visualization tool, the graph of ranking from solution paths (GRASP), designed to provide instant insights into model selection uncertainty. Additionally, GRASP accounts for the uncertainty of variable importance, enabling decision‐makers to assess each variable under uncertainty. Based on an innovative selection procedure that utilizes the entire solution path, a feature importance score and bootstrap techniques, GRASP effectively visualizes the uncertainty of model selection, as demonstrated by our numerical examples. Furthermore, we propose a novel measure of uncertainty based on GRASP, providing a single‐number summary of selection uncertainty. This measure incorporates the concept of the flat norm, traditionally used in geometry and physics. Our simulation studies and numerical examples confirm that this measure accurately and robustly quantifies uncertainty. 
    more » « less