Numerical Analysis for Nematic Electrolytes

Ľubomír Baňas,Robert Lasarzik,Andreas Prohl

DOI: https://doi.org/10.48550/arXiv.2004.10836

2020-04-23

Abstract:We consider a system of nonlinear PDEs modeling nematic electrolytes, and construct a dissipative solution with the help of its implementable, structure-inheriting space-time discretization. Computational studies are performed to study the mutual effects of electric, elastic, and viscous effects onto the molecules in a nematic electrolyte.

Numerical Analysis,Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the electro - dynamic modeling and numerical analysis of nematic electrolytes. Specifically, the researchers constructed a nonlinear partial differential equation (PDE) system to describe the physical phenomena in nematic electrolytes and constructed a dissipative solution by implementing a structure - preserving space - time discretization method. The main focuses are:

1. **Electrical, elastic and viscous effects**: Study the influence of these effects on liquid crystal molecules.
2. **Solvability of the model**: Prove the solvability of the proposed model in a simplified form and approximate it by the finite element method.
3. **Numerical experiments**: Verify the theoretical results through computational research.

### Specific problem description

The nonlinear PDE system discussed in the paper includes:
- **Ericksen - Leslie equations**: Used to describe the average orientation distribution of molecules in the liquid crystal phase.
- **Nernst - Planck - Poisson system**: Used to simulate the phenomena caused by electrolytes.

The specific PDE system is as follows:
\[
\begin{aligned}
& \partial_t \mathbf{v}+(\mathbf{v} \cdot \nabla) \mathbf{v}-\nu \Delta \mathbf{v}+\nabla \cdot(\mathbf{d} \otimes \Delta \mathbf{d}+\epsilon_a (\mathbf{d} \cdot \nabla \Phi) \nabla \Phi)+(n_+ - n_-) \nabla \Phi+\nabla p = 0, \\
& \nabla \cdot \mathbf{v} = 0, \\
& \partial_t \mathbf{d}+(\mathbf{v} \cdot \nabla) \mathbf{d}-(\mathbf{I}-\mathbf{d} \otimes \mathbf{d})(\Delta \mathbf{d}+\epsilon_a (\mathbf{d} \cdot \nabla \Phi) \nabla \Phi) = 0, \\
& |\mathbf{d}| = 1, \\
& -\nabla \cdot(\epsilon(\mathbf{d}) \nabla \Phi) = n_+ - n_-, \\
& \partial_t n_\pm+(\mathbf{v} \cdot \nabla) n_\pm-\nabla \cdot(\epsilon(\mathbf{d})(\nabla n_\pm + n_\pm \nabla \Phi)) = 0,
\end{aligned}
\]
where:
- \(\mathbf{v}\) is the macroscopic velocity of the liquid crystal fluid;
- \(\mathbf{d}\) is the local orientation of liquid crystal molecules;
- \(\Phi\) is the electric potential;
- \(n_\pm\) are the concentrations of positive and negative ions;
- \(p\) is the pressure of the nematic electrolyte;
- \(\epsilon(\mathbf{d})=\mathbf{I}+\epsilon_a \mathbf{d} \otimes \mathbf{d}\) is the dielectric constant matrix.

### Main challenges

1. **Structure - preserving discretization scheme**: It is necessary to construct a discretization scheme that can inherit the properties of the original system (such as incompressibility and unit - length constraint) and ensure stability.
2. **Convergence and the concept of solution**: In the three - dimensional case, it is very difficult to directly construct a weak solution, so the concept of dissipative solution is introduced.

### Main contributions of the paper

- Proposed a structure - preserving space - time discretization method that can generate a series of approximations.