Finite element methods for multicomponent convection-diffusion

Francis R A Aznaran,Patrick E Farrell,Charles W Monroe,Alexander J Van-Brunt

DOI: https://doi.org/10.1093/imanum/drae001

2024-04-27

IMA Journal of Numerical Analysis

Abstract:Abstract We develop finite element methods for coupling the steady-state Onsager–Stefan–Maxwell (OSM) equations to compressible Stokes flow. These equations describe multicomponent flow at low Reynolds number, where a mixture of different chemical species within a common thermodynamic phase is transported by convection and molecular diffusion. Developing a variational formulation for discretizing these equations is challenging: the formulation must balance physical relevance of the variables and boundary data, regularity assumptions, tractability of the analysis, enforcement of thermodynamic constraints, ease of discretization and extensibility to the transient, anisothermal and nonideal settings. To resolve these competing goals, we employ two augmentations: the first enforces the definition of mass-average velocity in the OSM equations, while its dual modifies the Stokes momentum equation to enforce symmetry. Remarkably, with these augmentations we achieve a Picard linearization of symmetric saddle point type, despite the equations not possessing a Lagrangian structure. Exploiting structure mandated by linear irreversible thermodynamics, we prove the inf-sup condition for this linearization, and identify finite element function spaces that automatically inherit well-posedness. We verify our error estimates with a numerical example, and illustrate the application of the method to nonideal fluids with a simulation of the microfluidic mixing of hydrocarbons.

mathematics, applied

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to effectively combine the steady - state Onsager - Stefan - Maxwell (OSM) equations with the compressible Stokes flow equations in multi - component fluids to describe the convection and molecular diffusion processes of different chemical substance mixtures in the same thermodynamic phase under low - Reynolds - number conditions. Specifically, the authors developed a finite - element method to discretize these equations, which involves several challenges:

1. **Physical relevance of variables and boundary data**: It is necessary to ensure that the selected variables and boundary conditions can accurately reflect physical phenomena.
2. **Regularity assumptions**: Mathematical assumptions required to ensure the existence and uniqueness of solutions.
3. **Tractability of analysis**: The mathematical form of the equations should be convenient for theoretical analysis.
4. **Enforcement of thermodynamic constraints**: Ensure that the system satisfies the basic principles of thermodynamics.
5. **Simplicity and extensibility of discretization**: The method should be easy to implement and be able to be extended to transient, non - isothermal and non - ideal cases.

To address these challenges, the authors adopted two enhancement techniques:
- The first enhancement technique is used to enforce the definition of mass - average velocity in the OSM equations.
- The second enhancement technique modifies the Stokes momentum equation to enforce symmetry.

Through these enhancement techniques, the authors successfully achieved Picard linearization with a symmetric saddle - point structure, although the original equations do not have a Lagrangian structure. Using the linear irreversible thermodynamic structure, the authors proved the inf - sup condition of this linearization and determined the finite - element function spaces that automatically inherit continuous well - posedness.

In addition, the authors also verified the error estimates through numerical examples and demonstrated the application of this method in non - ideal fluids, such as the simulation of hydrocarbon mixing in microfluidic devices.

In summary, the main contribution of this paper lies in proposing a new variational formula for describing the system equations of non - ideal miscible isothermal multi - component fluids and developing a structure - preserving finite - element discretization method that can automatically inherit continuous well - posedness.