1
|
Ronsin OJJ, Jang D, Egelhaaf HJ, Brabec CJ, Harting J. Phase-Field Simulation of Liquid-Vapor Equilibrium and Evaporation of Fluid Mixtures. ACS APPLIED MATERIALS & INTERFACES 2021; 13:55988-56003. [PMID: 34792348 DOI: 10.1021/acsami.1c12079] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
In solution processing of thin films, the material layer is deposited from a solution composed of several solutes and solvents. The final morphology and hence the properties of the film often depend on the time needed for the evaporation of the solvents. This is typically the case for organic photoactive or electronic layers. Therefore, it is important to be able to predict the evaporation kinetics of such mixtures. We propose here a new phase-field model for the simulation of evaporating fluid mixtures and simulate their evaporation kinetics. Similar to the Hertz-Knudsen theory, the local liquid-vapor (LV) equilibrium is assumed to be reached at the film surface and evaporation is driven by diffusion away from this gas layer. In the situation where the evaporation is purely driven by the LV equilibrium, the simulations match the behavior expected theoretically from the free energy: for evaporation of pure solvents, the evaporation rate is constant and proportional to the vapor pressure. For mixtures, the evaporation rate is in general strongly time-dependent because of the changing composition of the film. Nevertheless, for highly nonideal mixtures, such as poorly compatible fluids or polymer solutions, the evaporation rate becomes almost constant in the limit of low Biot numbers. The results of the simulation have been successfully compared to experiments on a polystyrene-toluene mixture. The model allows to take into account deformations of the liquid-vapor interface and, therefore, to simulate film roughness or dewetting.
Collapse
Affiliation(s)
- Olivier J J Ronsin
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Fürther Straße 248, 90429 Nürnberg, Germany
- Department of Chemical and Biological Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429 Nürnberg, Germany
| | - DongJu Jang
- ZAE Bayern─Solar Factory of the Future, Energy Campus Nürnberg, Fürther Straße 250, 90429 Nürnberg, Germany
| | - Hans-Joachim Egelhaaf
- ZAE Bayern─Solar Factory of the Future, Energy Campus Nürnberg, Fürther Straße 250, 90429 Nürnberg, Germany
- Institute of Materials for Electronics and Energy Technology (i-MEET), Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Immerwahrstrasse 2, 91058 Erlangen, Germany
| | - Christoph J Brabec
- Institute of Materials for Electronics and Energy Technology (i-MEET), Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Immerwahrstrasse 2, 91058 Erlangen, Germany
| | - Jens Harting
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Fürther Straße 248, 90429 Nürnberg, Germany
- Department of Chemical and Biological Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429 Nürnberg, Germany
- Department of Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429 Nürnberg, Germany
| |
Collapse
|
2
|
Tönsmann M, Scharfer P, Schabel W. Critical Solutal Marangoni Number Correlation for Short-Scale Convective Instabilities in Drying Poly(vinyl acetate)-Methanol Thin Films. Polymers (Basel) 2021; 13:polym13172955. [PMID: 34502995 PMCID: PMC8433935 DOI: 10.3390/polym13172955] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Revised: 08/27/2021] [Accepted: 08/28/2021] [Indexed: 12/22/2022] Open
Abstract
A new empiric correlation for the critical solutal Marangoni number as function of the Péclet and Schmidt numbers is proposed. It is based on previously published experimental flow field data in drying poly(vinyl acetate)-methanol films with an initial thickness in the range of 20–100 μm and an initial solvent load of 1 to 2 gMeOH/gPVAc, as well as newly derived concentration profile measurements and 1D drying simulations. The analysis accounts for realistic transient material properties and describes the occurrence of short-scale convective Marangoni (in)stabilities during the entire drying process with an accuracy of 9%. In addition, the proposed correlation qualitatively follows trends known from theory. As convective Marangoni instabilities in drying polymer films may induce surface deformations, which persist in the dry film, the correlation may facilitate future process design for either thin films with uniform thickness or deliberate self-assembly.
Collapse
|
3
|
Gracia‐Medrano‐Bravo V, Scharfer P, Schabel W. On the drying kinetics of non‐spherical particle‐filled polymer films: A numerical study. AIChE J 2021. [DOI: 10.1002/aic.17398] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Philip Scharfer
- Thin Film Technology (TFT) Karlsruhe Institute of Technology (KIT) Karlsruhe Germany
| | - Wilhelm Schabel
- Thin Film Technology (TFT) Karlsruhe Institute of Technology (KIT) Karlsruhe Germany
| |
Collapse
|
4
|
Börnhorst T, Scharfer P, Schabel W. Drying Kinetics from Micrometer- to Nanometer-Scale Polymer Films: A Study on Solvent Diffusion, Polymer Relaxation, and Substrate Interaction Effects. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:6022-6031. [PMID: 33947182 DOI: 10.1021/acs.langmuir.1c00641] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The drying behavior of two different polymers [polyvinyl pyrrolidone (PVP) and polyisobutylene (PIB)] with different glass transition temperatures are investigated and compared as a function of film thickness from micrometer (∼3 μm) to nanometer scale (∼10 nm). The focus of this study is to distinguish between solvent diffusion, polymer relaxation, and substrate confinement of polymer chain mobility toward the interface as the dominating mechanism of drying kinetics. Relaxation kinetics becomes more dominant when the film thickness is reduced, which is shown experimentally for the first time for nanometer-scale film thicknesses. Identical drying curves regardless of the film thickness of PVP/methanol indicate the limitation of solvent transport by relaxation kinetics. The viscoelastic relaxation behavior of the polymer/solvent film is modeled by a Maxwell element. The results are in accordance with the experimental drying curves and allow for the determination of the characteristic relaxation time. Relaxation limitation becomes relevant at high diffusion Deborah numbers when the relaxation time-which is a function of the deployed material and the polymer/solvent composition-is higher than the characteristic diffusion time in the film. The latter is a function of the polymer/solvent composition and the thickness of the film. Drying curves of PIB/toluene films show additional effect in a substrate-near region of about 5 nm in which polymer chain mobility is confined, resulting in decelerated solvent diffusion. Although this effect near the substrate interface is expected to be present regardless of the film thickness, it becomes more dominant when the substrate-near region represents a significant fraction of the total film thickness. The key to the derived methodology for characterization of the polymer/solvent drying process is to vary dry film thickness from micrometers to a few nanometers which allows us to determine the dominating mechanism of drying kinetics.
Collapse
Affiliation(s)
- Tobias Börnhorst
- Thin Film Technology (TFT), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Philip Scharfer
- Thin Film Technology (TFT), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Wilhelm Schabel
- Thin Film Technology (TFT), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| |
Collapse
|
5
|
Transient Three-Dimensional Flow Field Measurements by Means of 3D µPTV in Drying Poly(Vinyl Acetate)-Methanol Thin Films Subject to Short-Scale Marangoni Instabilities. Polymers (Basel) 2021; 13:polym13081223. [PMID: 33920103 PMCID: PMC8068913 DOI: 10.3390/polym13081223] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2021] [Revised: 04/06/2021] [Accepted: 04/08/2021] [Indexed: 12/17/2022] Open
Abstract
Convective Marangoni instabilities in drying polymer films may induce surface deformations, which persist in the dry film, deteriorating product performance. While theoretic stability analyses are abundantly available, experimental data are scarce. We report transient three-dimensional flow field measurements in thin poly(vinyl acetate)-methanol films, drying under ambient conditions with several films exhibiting short-scale Marangoni convection cells. An initial assessment of the upper limit of thermal and solutal Marangoni numbers reveals that the solutal effect is likely to be the dominant cause for the observed instabilities.
Collapse
|