Rainey Lyons,Emilio N. M. Cirillo,Adrian Muntean
Abstract:We study a nonlinear coupled parabolic system with non-local drift terms modeling at the continuum level the inter-species interaction within a ternary mixture that allows the evaporation of one of the species. In the absence of evaporation, the proposed system coincides with the hydrodynamic limit of a stochastic interacting particle system of Blume-Capel-type driven by the Kawasaki dynamics. Similar governing dynamics are found in models used to study morphology formation in the design of organic solar cells, thin adhesive bands, and other applications. We investigate the well-posedness of the target system and present preliminary numerical simulations which incorporate evaporation into the model. We employ a finite volumes scheme to construct approximations of the weak solution and illustrate how the evaporation process can affect the shape and connectivity of the evolving-in-time morphologies.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the phase separation and morphology formation in ternary mixtures under evaporation conditions. Specifically, the author studied a nonlinear coupled parabolic system that contains non - local drift terms for modeling the interactions between species in a ternary mixture at the continuous level, where one species is allowed to evaporate. The goal of the paper is to explore the well - posedness of this system and, through numerical simulations, show how evaporation affects the shape and connectivity of the evolving morphology over time.
### Background and Motivation
In materials science, especially when designing thin films for organic solar cells, the phase separation process is very important. These applications usually involve two different solutes mixed in one solvent, allowing phase separation. The evaporation process plays a crucial role in this case because the evaporation of the solvent will affect the finally formed morphological structure. By controlling the evaporation mechanism, a specific morphology can be actively selected as the final configuration.
### Model Description
The paper first discusses the model without evaporation conditions, that is, two coupled non - local nonlinear parabolic equations:
\[
\begin{cases}
\partial_t m=\nabla\cdot(\nabla m - 2\beta(\phi - m^2)(\nabla J * m)),&(t, x)\in(0, T)\times\Omega\\
\partial_t\phi=\nabla\cdot(\nabla\phi - 2\beta m(1 - \phi)(\nabla J * m)),&(t, x)\in(0, T)\times\Omega
\end{cases}
\]
where \(t\) and \(x\) represent the time and space variables respectively, \(m\) represents the average spin density (also known as magnetization), and \(\phi\) represents the solute volume concentration. In addition, \(\Omega\subset\mathbb{R}^d\) is a cube with periodic boundary conditions, \(T\) is some positive time, \(\beta> 0\) is a constant representing the inverse temperature, and \(J\in C^2_+(\bar{\Omega})\) is a symmetric compact - supported potential function and satisfies \(\int_{\Omega}J(r)dr = 1\).
### Introduction of the Evaporation Process
To simulate the evaporation process, the author introduced a linear evaporation term into the above system, obtaining a new system:
\[
\begin{cases}
\partial_t m=\nabla\cdot(\nabla m - 2\beta(\phi - m^2)(\nabla J * m)),&(t, x)\in(0, T)\times\Omega\\
\partial_t\phi=\nabla\cdot(\nabla\phi - 2\beta m(1 - \phi)(\nabla J * m))+F(\phi),&(t, x)\in(0, T)\times\Omega
\end{cases}
\]
where \(F: [0, 1]\to\mathbb{R}\) is a bounded, non - increasing, Lipschitz - continuous function and \(F(1) = 0\).
### Main Results
The main result of the paper is to prove the well - posedness of this system and maintain the inequality \(0\leq|m|\leq\phi\leq1\) on \(Q_T\). In particular, the author focuses on the existence and uniqueness of weak solutions and gives specific definitions and theorems.
### Numerical Simulation
The author also carried out preliminary numerical simulations to show how the evaporation process affects the shape and connectivity of the evolving morphology over time. These simulations use the finite - volume method to construct approximations of weak solutions.