Abstract
In this paper we consider the system of partial differential equations describing the stress-assisted diffusion of a solute into an elastic material, and introduce and analyze a Banach spaces-based variational approach yielding a new mixed-primal finite element method for its numerical solution. The elasticity model involved, which is initially defined according to the constitutive relation given by Hooke’s law, and whose momentum equation holds with a concentration-depending source term, is reformulated by using the non-symmetric pseudostress tensor and the displacement as the only unknowns of the associated mixed scheme, in addition to assuming a Dirichlet boundary condition for the latter. In turn, the diffusion equation, whose diffusivity function and source term depend on the stress and the displacement of the solid, respectively, is set in primal form in terms of the concentration unknown and a Dirichlet boundary condition for it as well. The resulting coupled formulation is rewritten as an equivalent fixed point operator equation, so that its unique solvability is established by employing the classical Banach theorem along with the corresponding Babuška-Brezzi theory and the Lax-Milgram theorem. The aforementioned dependence of the diffusion coefficient and the subsequent treatment of this term in the continuous analysis, suggest to better look for the solid unknowns in suitable Lebesgue spaces. The discrete analysis is performed similarly, and the Brouwer theorem yields existence of a Galerkin solution. A priori error estimates are derived, and rates of convergence for specific finite element subspaces satisfying the required discrete inf-sup conditions, are established in 2D. Finally, several numerical examples illustrating the performance of the method and confirming the theoretical convergence, are reported.




Similar content being viewed by others
Data Availability
Enquiries about data availability should be directed to the authors.
References
Álvarez, M., Gatica, G.N., Ruiz-Baier, R.: An augmented mixed-primal finite element method for a coupled flow-transport problem. ESAIM Math. Model. Numer. Anal. 49(5), 1399–1427 (2015)
Álvarez, M., Gatica, G.N., Ruiz-Baier, R.: A mixed-primal finite element approximation of a steady sedimentation-consolidation system. Math. Models Methods Appl. Sci. 26(5), 867–900 (2016)
An, Y., Jiang, H.: A finite element simulation on transient large deformation and mass diffusion in electrodes for lithium ion batteries. Model. Simul. Materials Sci. Engrg. 21(7), 074007 (2013)
Barrientos, M.A., Gatica, G.N., Stephan, E.P.: A mixed finite element method for nonlinear elasticity: two-fold saddle point approach and a-posteriori error estimate. Numer. Math. 91(2), 197–222 (2002)
Bernardi, C., Canuto, C., Maday, Y.: Generalized inf-sup conditions for Chebyshev spectral approximation of the Stokes problem. SIAM J. Numer. Anal. 25(6), 1237–1271 (1988)
Boffi, D., Brezzi, F., Fortin, M.: Mixed Finite Element Methods and Applications. Springer Series in Computational Mathematics, 44. Springer, Heidelberg, (2013)
Brezzi, F., Fortin, M.: Mixed and Hybrid Finite Element Methods. Springer-Verlag, New York (1991)
Cherubini, C., Filippi, S., Gizzi, A., Ruiz-Baier, R.: A note on stress-driven anisotropic diffusion and its role in active deformable media. J. Theoret. Biol. 430(7), 221–228 (2017)
Ciarlet, P.: Linear and Nonlinear Functional Analysis with Applications. Society for Industrial and Applied Mathematics, Philadelphia, PA (2013)
Colmenares, E., Gatica, G.N., Oyarzúa, R.: Analysis of an augmented mixed-primal formulation for the stationary Boussinesq problem. Numer. Methods Partial Differential Equations 32(2), 445–478 (2016)
Durán, R.G.: Error analysis in \(\rm L ^p\), \(1 \le p \le \infty \), for mixed finite element methods for linear and quasi-linear elliptic problems. RAIRO Modél. Math. Anal. Numér. 22(3), 371–387 (1988)
Ern, A., Guermond, J.-L.: Theory and Practice of Finite Elements. Applied Mathematical Sciences, 159. Springer-Verlag, New York, (2004)
Foster, J.M., Chapman, S.J., Richardson, G., Protas, B.: A mathematical model for mechanically-induced deterioration of the binder in lithium-ion electrodes. SIAM J. Appl. Math. 77(6), 2172–2198 (2017)
Fromm, S.J.: Potential space estimates for Green potentials in convex domains. Proc. Amer. Math. Soc. 119(1), 225–233 (1993)
Gatica, G.N.: A Simple Introduction to the Mixed Finite Element Method. Theory and Applications. SpringerBriefs in Mathematics. Springer, Cham (2014)
Gatica, G.N., Gatica, L.F., Sequeira, F.A.: A priori and a posteriori error analyses of a pseudostress-based mixed formulation for linear elasticity. Comput. Math. Appl. 71(2), 585–614 (2016)
Gatica, G.N., Gómez-Vargas, B., Ruiz-Baier, R.: Analysis and mixed-primal finite element discretisations for stress-assisted diffusion problems. Comput. Methods Appl. Mech. Engrg. 337, 411–438 (2018)
Gatica, G.N., Gómez-Vargas, B., Ruiz-Baier, R.: Formulation and analysis of fully-mixed methods for stress-assisted diffusion problems. Comput. Math. Appl. 77(5), 1312–1330 (2019)
Gatica, G.N., Inzunza, C.: On the well-posedness of Banach spaces-based mixed formulations for the nearly incompressible Navier-Lamé and Stokes equations. Comput. Math. Appl. 102, 87–94 (2021)
Gatica, G.N., Márquez, A., Oyarzúa, R., Rebolledo, R.: Analysis of an augmented fully-mixed approach for the coupling of quasi-Newtonian fluids and porous media. Comput. Methods Appl. Mech. Engrg. 270, 76–112 (2014)
Gatica, G.N., Meddahi, S., Ruiz-Baier, R.: An \({{\rm L}}^{p}\) spaces-based formulation yielding a new fully mixed finite element method for the coupled Darcy and heat equations. IMA J. Numer. Anal., https://doi.org/10.1093/imanum/drab063
Lewicka, M., Mucha, P.B.: A local and global well-posedness results for the general stress-assisted diffusion systems. J. Elasticity 123(1), 19–41 (2016)
Manda, M.L., Shepard, R., Fair, B., Massoud, H.Z.: Stress-assisted diffusion of boron and arsenic in silicon. Mat. Res. Soc. Symp. Proc. 36, 71–76 (1985)
Roy, S., Vengadassalam, K., Wang, Y., Park, S., Liechti, K.M.: Characterization and modeling of strain assisted diffusion in an epoxy adhesive layer. Int. J. Solids Struct. 43, 27–52 (2006)
Si, H.: TetGen: A Quality Tetrahedral Mesh Generator and 3D Delaunay Triangulator v.1.5 User’s manual, Tech. Report 13, Weierstrass Institute for Applied Analysis and Stochastics, Berlin, (2013)
Taralova, V., Iliev, O., Efendiev, Y.: Derivation and numerical validation of a homogenized isothermal Li-ion battery model. J. Engr. Math. 101, 1–27 (2016)
Yost, F.G., Amos, D.E., Roming Jr., A.D.: Stress-driven diffusive voiding of aluminum conductor lines. Proc. Int. Rel. Phys. Symp., 193–201 (1989)
Funding
This research was partially supported by ANID-Chile through the projects Centro de Modelamiento Matemático (FB210005), Anillo of Computational Mathematics for Desalination Processes (ACT 210087), and the Becas Chile Programme for national students; by Centro de Investigación en Ingeniería Matemática (CI2MA), Universidad de Concepción; and by Universidad Nacional (Costa Rica), through the project 0140-20.
Author information
Authors and Affiliations
Contributions
All authors contributed equally to the study conception and development of this work. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors have no relevant financial or non-financial interests to disclose
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Gatica, G.N., Inzunza, C. & Sequeira, F.A. A Pseudostress-Based Mixed-Primal Finite Element Method for Stress-Assisted Diffusion Problems in Banach Spaces. J Sci Comput 92, 103 (2022). https://doi.org/10.1007/s10915-022-01959-9
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-022-01959-9