[Paper Review] Global existence of weak solutions of a model for electrolyte solutions - Part 1: Two-component case

Matthias Herz, Peter Knabner|arXiv (Cornell University)|May 24, 2016

Navier-Stokes equation solutions30 references3 citations

TL;DR

This paper establishes the global existence and uniqueness of weak solutions for the Darcy–Poisson–Nernst–Planck system modeling two-component electrolyte solutions with arbitrary valencies $ z_1 > 0 > z_2 $, using weighted test functions $ \varphi = |z_l|c_l $ to derive uniform a priori estimates, thereby overcoming limitations of prior results that required symmetry, electroneutrality, or volume-additivity constraints in two and three spatial dimensions.

ABSTRACT

This paper analytically investigates the Darcy-Poisson-Nernst-Planck system. This system is a mathematical model for electrolyte solutions. In this paper, we consider electrolyte solutions, which consist of a neutral fluid and two suspended oppositely charged chemical species with arbitrary valencies z_1 > 0 > z_2. We prove global existence and uniqueness of weak solutions in two space dimensions and three space dimensions. So far, most of the existence results have been proven for symmetric electrolyte solutions. These solutions consist of a neutral fluid and two suspended charged chemical species with symmetric valencies +/-z. As many electrolyte solutions in biological applications and hydrodynamical applications are not symmetric, the presented extension of the previous existence results is an important step.

Motivation & Objective

• To extend existing existence results for electrolyte models beyond symmetric valencies $ \pm z $ to arbitrary opposite valencies $ z_1 > 0 > z_2 $.
• To remove restrictive assumptions such as electroneutrality or volume-additivity constraints in the analysis of electrolyte systems.
• To establish global existence and uniqueness of weak solutions in both two and three spatial dimensions for the Darcy–Poisson–Nernst–Planck system.
• To develop uniform a priori estimates for the solute concentrations independent of the electric field, which is nontrivial for asymmetric valencies.
• To provide a rigorous analytical foundation for modeling complex electrolyte systems in biological and hydrodynamic applications where symmetric valency assumptions do not hold.

Proposed method

• Employed weighted test functions $ \varphi = |z_l|c_l $ instead of standard $ \varphi = c_l $ to derive a priori estimates independent of the electric field.
• Applied Moser iteration to obtain uniform $ L^\infty $ bounds on the solute concentrations $ c_l $, ensuring regularity and compactness.
• Used the generalized Schauder fixed point theorem in a Galerkin approximation framework to construct approximate solutions.
• Combined Aubin-Lions lemma with weak* convergence to pass to the limit in nonlinear terms, including reaction and drift integrals.
• Established local existence via a time-discretization and continuation argument, with subinterval time steps chosen small enough to preserve a priori bounds.
• Used exponential dependence of estimates on time to prevent blow-up and enable global continuation beyond any finite time horizon.

Experimental results

Research questions

• RQ1 Can global existence and uniqueness of weak solutions be established for the Darcy–Poisson–Nernst–Planck system with two oppositely charged species of arbitrary valencies?
• RQ2 How can uniform a priori estimates for solute concentrations be derived when the valencies are asymmetric and not symmetrically paired?
• RQ3 Is it possible to avoid the electroneutrality or volume-additivity constraints that are commonly imposed in prior analytical treatments?
• RQ4 What mathematical techniques are required to control the coupling between electrostatics, fluid flow, and ion transport in the absence of symmetry?
• RQ5 Can the solution be extended globally in time despite the nonlinear coupling and lack of symmetry in the system?

Key findings

• Global existence and uniqueness of weak solutions are proven for the two-component Darcy–Poisson–Nernst–Planck system in both two and three spatial dimensions.
• The $ L^\infty $ bound on the solute concentrations $ c_l $ is uniformly controlled via Moser iteration, ensuring regularity and compactness.
• The use of weighted test functions $ \varphi = |z_l|c_l $ enables derivation of a priori estimates independent of the electric field, overcoming a key difficulty in asymmetric cases.
• The solution construction avoids the need for electroneutrality or volume-additivity constraints, broadening applicability to real-world electrolyte systems.
• The continuation argument relies on exponential dependence of a priori estimates on time, which prevents blow-up and allows extension to arbitrary large time intervals.
• This is the first global existence and uniqueness result for such systems with arbitrary valencies in three dimensions without restrictive physical assumptions.