A computational model for diffusion-reaction-deformation coupled problems and its finite element implementation

Abstract

Purpose

A diffusion-reaction-deformation coupled model is employed and implemented as a user-defined element (UEL) subroutine in the commercial finite element software package ABAQUS.

Design/methodology/approach

Chemical reaction and diffusion are treated as two distinct processes by introducing the extent of reaction and the diffusion concentration as two kinds of independent variables, for which the independent governing equations for chemical reaction and diffusion processes are obtained. Furthermore, an exponential form of chemical kinetics, instead of the linearly phenomenological relation, between the reaction rate and the chemical affinity is used to describe reaction process. As a result, complex chemical reaction can be simulated, no matter it is around or away from equilibrium.

Findings

Two numerical examples are presented, one for validation of the model and another for the modeling of the deflection of a plane caused by a chemical reaction.

Originality/value

1. Independent governing equations for diffusion and reaction processes are given. 2. An exponential relation between the reaction rate and its driving force is employed. 3. The UEL subroutine is used to implement the finite element procedure.

Citation

Qin, B., Zhou, Y. and Zhong, Z. (2022), "A computational model for diffusion-reaction-deformation coupled problems and its finite element implementation", Engineering Computations, Vol. 39 No. 3, pp. 837-857. https://doi.org/10.1108/EC-02-2021-0084