Introduction

The Chemical FEM physics is the responsible for providing an implementation for the multispecies reactive solute transport equation that can be used together with the Fem process to solve the steady state and transient reactive solute transport through a porous media such as rock, in 2D and 3D. The main equation solved by this physics is:

\[ k_i \frac{\partial C_i}{\partial t} = \nabla\cdot(D \nabla C_i)-v\nabla C_i + R_i \]

where:

\( C_i \) is the concentration of species in \(mol/m^3\)
\( D \) is the diffusion in \(m^2/s\)
\( v \) is the convective flux in \(m/s\)
\( R_i\) is the reaction of species in \(mol/m^3*s\)
\( k-i \) is the retardation fator of species

In order to solve the chemical problem we need to state appropiate initial and boundary conditions. Two types of boundary conditions be emplemented in GeMA, Dirichlet (Prescribed concentration of species in \(mol/m3\)) and Neumann (prescibed concentration flow in \(mol/m^3*s\)).

Figure 1: multispecies reactive solute transport process

Formulation Reference

The Chemical FEM physics plugin supports the following formulations:

Continuum Elements
Pipe (1D) Element

The FEM discretization of multispecies reactive solute transport equation used by the Chemical Fem physics plugin is based on the book Fundamentals of the Fundamentals of the Finite Element Method for Heat and Fluid Flow by Roland W. Lewis, Perumal Nithiarasu and Kankanhalli N. Seetharamu, 2004.

Plugin options

A reference manual documenting the set of state variables and material properties expected by the plugin, along with all of its supported configuration and result options can be found here.

Examples

Several example simulations involving steady state and transient flow calculations are available at the examples directory provided with the GeMA instalation. A summary of those examples can be found here.

GeMA project main page