Rekha Khot,Andres E. Rubiano,Ricardo Ruiz-Baier
Abstract:This paper aims first to perform robust continuous analysis of a mixed nonlinear formulation for stress-assisted diffusion of a solute that interacts with an elastic material, and second to propose and analyse a virtual element formulation of the model problem. The two-way coupling mechanisms between the Herrmann formulation for linear elasticity and the reaction-diffusion equation (written in mixed form) consist of diffusion-induced active stress and stress-dependent diffusion. The two sub-problems are analysed using the extended Babuška--Brezzi--Braess theory for perturbed saddle-point problems. The well-posedness of the nonlinearly coupled system is established using a Banach fixed-point strategy under the smallness assumption on data. The virtual element formulations for the uncoupled sub-problems are proven uniquely solvable by a fixed-point argument in conjunction with appropriate projection operators. We derive the a priori error estimates, and test the accuracy and performance of the proposed method through computational simulations.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the robust analysis and numerical simulation of the coupled - stress - assisted diffusion problem. Specifically, the authors aim to:
1. **Conduct continuous analysis**: Conduct a robust continuous analysis of a mixed nonlinear formula that describes the stress - assisted diffusion process of solutes interacting with elastic materials.
2. **Propose and analyze the Virtual Element Method (VEM)**: Propose and analyze a discretization scheme based on the Virtual Element Method for the above - mentioned model problem.
### Problem Background
In solid materials, the diffusion process can lead to the generation of mechanical stress, a phenomenon known as chemical or diffusion - induced stress. This stress - diffusion interaction is very important in many practical applications, such as lithium - ion batteries, hydrogen diffusion in silicone rubber, polymer coatings, semiconductor manufacturing, oxidation of silicon nanostructures, and enhancement of the electrical conductivity of soft tissues. In these processes, the local heterogeneity and anisotropy of stress will affect the diffusion pattern.
### Research Content
The authors focus on a simple stress - dependent diffusion coefficient model and assume that it reverts to Fick diffusion in the absence of stress. The study considers linear elastostatics and assumes that the total Cauchy stress is also caused by (possibly nonlinear) diffusion - induced isotropic stress. On the other hand, the diffusion coefficient controlling the solute distribution in the reaction - diffusion equation is regarded as a nonlinear function of the Cauchy stress.
### Method
To handle the case of nearly incompressible materials, the author adopts a mixed formula of displacement and Herrmann pressure, which makes the stability still robust when the Lame parameter \(\lambda\) tends to infinity. The robustness of the Lame parameter \(\mu\) is achieved by appropriately scaling the inf - sup condition of the divergence operator. In addition, the reaction - diffusion sub - problem also requires robustness with respect to the Lame parameter because the diffusion coefficient depends on the Cauchy stress, which in turn depends on the Lame coefficients.
### Numerical Method
The authors use the Virtual Element Method (VEM) to discretize the problem and prove the unique solvability of the discrete problem. They also derive a priori error estimates and test the accuracy and performance of the proposed method through computational simulation.
### Main Contributions
1. **Theoretical Analysis**: Establish the well - posedness of the nonlinear coupled system and prove the well - posedness of the system when the data is small enough using the Banach fixed - point strategy.
2. **Numerical Method**: Propose the Virtual Element Method and prove its robustness and convergence in different parameter ranges.
In summary, this paper is dedicated to solving the robust analysis and efficient numerical solution methods of the stress - assisted diffusion problem, especially providing new ideas and tools when dealing with complex material behaviors.