×
Home
Current Archive Editorial Board News Contact
Research paper

NUMERICAL HOMOGENIZATION AS A METHOD FOR MODELING CONSTRAINED DIFFUSION

By
Miloš Kojić ,
Miloš Kojić
Miljan Milošević ,
Miljan Milošević
Alessandro Grattoni ,
Alessandro Grattoni
Mauro Ferrari
Mauro Ferrari

Abstract

Material transport by diffusion is present in nature, as well as in living organisms as a vital process. The basic and generally used is a phenomenologically established law, named Fick’s law, It states that the material flux is proportional to the concentration gradient and that the material flows from the higher to lower concentration. However, if the diffusion occurs in nano space, the interaction of the transported molecules produces deviation from Fick’s law leading to a so-called restrained, or retarded diffusion. Using the Molecular Dynamics (MD) methodology, we have calculated the effects of the molecule-solid interaction by introducing the scale functions and used them in our finite element code PAK to find the equivalent bulk diffusivity and equivalent distance from the solid surface. This computational procedure is called numerical homogenization, where the equivalent diffusion parameters are evaluated from the mass release curves. Verification of our multiscale modeling of the retarded diffusion is demonstrated on glucose transport through a device with nanochannels (NDS). Our methodology is illustrated by two examples – one to show the effects of the molecular interaction within a space with solid spheres and several porosities, and another where the equivalent diffusion coefficients are computed using images of pancreatic tumor tissue.

References

1.
J. G, M D. Diffusion in supersaturated solutions. II. Glucose solutions *J Am Chem Soc*. 1953;75(16):3900–4.
2.
A. G, H. S, D. F, A. Z, S. GJ, L. H, et al. Nanochannel technology for constant delivery of chemotherapeutics: beyond metronomic administration. *Pharm Res*. 2011;28(2):292–300.
3.
J. KE, J. TM, V. C, M. TR, R. C, D. C, et al. Transport properties of pancreatic cancer describe gemcitabine delivery and response. *J Clin Invest*. 2014;124(4).
4.
M. K, M. M, A Z. *Computational Models in Biomedical Engineering – Finite Element Models Based on Smeared Physical Fields: Theory, Solutions, and Software*. 2022.
5.
M. K, M. M, N. K, K. K, M. F, A Z. A multiscale MD–FE model of diffusion in composite media with internal surface interaction based on numerical homogenization procedure. *Comput Methods Appl Mech Engrg*. 2014;269:123–38.
6.
M. K, M. M, N. K, J. KE, B. FJ, M. F, et al. Mass release curves as the constitutive curves for modeling diffusive transport within biological tissue. *Comput Biol Med*. 2018;92:156–67.
7.
M. K, M. M, S. W, E. B, M. F, A Z. Mass partitioning effects in diffusion transport. *Phys Chem Chem Phys*. 2015;17:20630–5.
8.
M. K, M. M, N. K, M. F, A Z. On diffusion in nanospace. *J Serbian Soc Comp Mechanics*. 2011;5(1):84–109.
9.
M. K, A. Z, M. M, V. I, N. K, M. R, et al. Transport in biological systems. *J Serbian Soc Comp Mechanics*. 2011;5(2):101–28.
10.
N. K, A. K, M K. Numerical determination of the solvent diffusion coefficient in a concentrated polymer solution. *Comm Num Meth Eng*. 2006;22:1003–13.
11.
N. K, J. PM, G. LG, K. RW, M. K, L KD. Ion electrodiffusion governs silk electrogelation. *Soft Matter*. 2012;8:6897–905.
12.
A. Z, M. K, M. M, N. K, F. H, M. F, et al. Hierarchical modeling of diffusive transport through nanochannels by coupling molecular dynamics with finite element method. *J Comput Phys*. 2011;230:5722–31.

Citation

Article metrics

Google scholar: See link

The statements, opinions and data contained in the journal are solely those of the individual authors and contributors and not of the publisher and the editor(s). We stay neutral with regard to jurisdictional claims in published maps and institutional affiliations.