Diffusion at solid-solid interfaces
Updated 14 Jan 2022
Finite volume methods implemented in Matlab to model 1D solute diffusion at solid-solid interfaces, accounting for energy barriers and interfacial segregation. The derivation of the models has been published in: F.D. León-Cázares and E.I. Galindo-Nava. Phys. Rev. Mat., 5 (2021) 123802. Please cite this publication if you benefit from this repository.
Coded in Matlab R2021a. There is no need for additional packages.
The Examples.m file showcases the usage of all the functions in this repository. The user can open the example of their preference and change the input parameters. If the simulation diverges, a smaller convergence criterion ~D*dt/dx^2<0.5 (printed by the functions) should be selected. Relevant plots of the concentrations and fluxes are also added for each diffusion case.
These functions solve the 1D diffusion equations for a compendium of model bicrystal interfaces. The geometry can be set to cartesian, cylindrical or spherical coordinates. Boundary conditions 'perm' of constant concentrations at both surfaces were implemented for all these cases:
- traps1D_perm.m - Interfacial energy barrier.
- traps1D_perm_t.m - Monolayer interfacial trap.
- traps1D_perm_t_diff.m - Diffuse interfacial trap.
- traps1D_perm_t_inh.m - Monolayer inhomogeneous interfacial trap.
Additionally, the monolayer interfacial trap case is implemented for other two sets of boundary conditions: 'open system' with zero flux (left) and constant concentration (right) surfaces, and 'closed system' with two no flux surfaces.
- traps1D_open_t.m - Monolayer interfacial trap.
- traps1D_closed_t.m - Monolayer interfacial trap.
- energy_landscape_plot.m - Plotting an energy landscape of a monolayer or no trap.
- energy_landscape_plot_diff.m - Plotting the energy landscape of a diffuse trap.
León-Cázares, Fernando D., and Enrique I. Galindo-Nava. “General Model for the Kinetics of Solute Diffusion at Solid-Solid Interfaces.” Physical Review Materials, vol. 5, no. 12, American Physical Society (APS), Dec. 2021, doi:10.1103/physrevmaterials.5.123802.
MATLAB Release Compatibility
Platform CompatibilityWindows macOS Linux
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!