No Arabic 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 equilibrium is assumed to be reached at the film surface and evaporation is driven by diffusion away from this gas layer. The simulations match the behavior expected theoretically from the free energy, as well as experimental results. 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 non-ideal mixtures, such as poorly compatible fluids or polymer solutions, the evaporation rate becomes almost constant in the limit of low Biot numbers.
The performance of solution-processed solar cells strongly depends on the geometrical structure and roughness of the photovoltaic layers formed during film drying. During the drying process, the interplay of crystallization and liquid-liquid demixing leads to the structure formation on the nano- and microscale and to the final rough film. In order to better understand how the film structure can be improved by process engineering, we aim at theoretically investigating these systems by means of phase-field simulations. We introduce an evaporation model based on the Cahn-Hilliard equation for the evolution of the fluid concentrations coupled to the Allen-Cahn equation for the liquid-vapour phase transformation. We demonstrate its ability to match the experimentally measured drying kinetics and study the impact of the parameters of our model. Furthermore, the evaporation of solvent blends and solvent-vapour annealing are investigated. The dry film roughness emerges naturally from our set of equations, as illustrated through preliminary simulations of spinodal decomposition and film drying on structured substrates.
A phase-field crystal model based on the density-field approach incorporating high-order interparticle direct correlations is developed to study vapor-liquid-solid coexistence and transitions within a single continuum description. Conditions for the realization of the phase coexistence and transition sequence are systematically analyzed and shown to be satisfied by a broad range of model parameters, demonstrating the high flexibility and applicability of the model. Both temperature-density and temperature-pressure phase diagrams are identified, while structural evolution and coexistence among the three phases are examined through dynamical simulations. The model is also able to produce some temperature and pressure related material properties, including effects of thermal expansion and pressure on equilibrium lattice spacing, and temperature dependence of saturation vapor pressure. This model can be used as an effective approach for investigating a variety of material growth and deposition processes based on vapor-solid, liquid-solid, and vapor-liquid-solid growth.
We present general calculations allowing to express the thermodynamical coefficients and thermophysical properties (compressibility, thermal coefficients and heat capacities) of a material composed of a mixture of two constituents or phases, regardless of the equations of state considered for each of the constituant or for the mixture. We consider mixtures under complete thermodynamical equilibrium, either with mass exchange between the two constituents (phase change), such as an oil and gas mixture (black-oil) and water-vapor system or with two immiscible phases, such as air and water mixtures (foam). Keywords Thermodynamic coefficient, Thermophysical properties, thermal coefficients, speed of sound, two phase systems, liquid vapor equilibrium
We investigate numerically the behaviour of a phase-separating mixture of a blue phase I liquid crystal with an isotropic fluid. The resulting morphology is primarily controlled by an inverse capillary number, $chi$, setting the balance between interfacial and elastic forces. When $chi$ and the concentration of the isotropic component are both low, the blue phase disclination lattice templates a cubic array of fluid cylinders. For larger $chi$, the isotropic phase arranges primarily into liquid emulsion droplets which coarsen very slowly, rewiring the blue phase disclination lines into an amorphous elastic network. Our blue phase/simple fluid composites can be externally manipulated: an electric field can trigger a morphological transition between cubic fluid cylinder phases with different topologies.
Multicomponent systems are ubiquitous in nature and industry. While the physics of few-component liquid mixtures (i.e., binary and ternary ones) is well-understood and routinely taught in undergraduate courses, the thermodynamic and kinetic properties of $N$-component mixtures with $N>3$ have remained relatively unexplored. An example of such a mixture is provided by the intracellular fluid, in which protein-rich droplets phase separate into distinct membraneless organelles. In this work, we investigate equilibrium phase behavior and morphology of $N$-component liquid mixtures within the Flory-Huggins theory of regular solutions. In order to determine the number of coexisting phases and their compositions, we developed a new algorithm for constructing complete phase diagrams, based on numerical convexification of the discretized free energy landscape. Together with a Cahn-Hilliard approach for kinetics, we employ this method to study mixtures with $N=4$ and $5$ components. We report on both the coarsening behavior of such systems, as well as the resulting morphologies in three spatial dimensions. We discuss how the number of coexisting phases and their compositions can be extracted with Principal Component Analysis (PCA) and K-Means clustering algorithms. Finally, we discuss how one can reverse engineer the interaction parameters and volume fractions of components in order to achieve a range of desired packing structures, such as nested `Russian dolls and encapsulated Janus droplets.