2010, Münch, Andreas, Please, Colin P., Wagner, Barbara

We consider a mathematical model of spin coating of a single polymer blended in a solvent. The model describes the one-dimensional development of the thin layer of the mixture as the layer thins due to flow created by a balance of viscous forces and centrifugal forces and due to evaporation of the solvent. In the model both the diffusivity of the solvent in the polymer and the viscosity of the mixture are very rapidly varying functions of the solvent volume fraction. Guided by numerical solutions an asymptotic analysis reveals a number of different possible behaviours of the thinning layer dependent on the nondimensional parameters describing the system. The main practical interest is in controlling the appearance and development of a ``skin'' on the polymer where the solvent concentration reduces rapidly on the outer surface leaving the bulk of the layer still with high concentrations of solvent. The critical parameters controlling this behaviour are found to be eps the ratio of the diffusion to advection time scales, delta the ratio of the evaporation to advection time scales and exp(-gamma), the ratio of the diffusivity of the initial mixture and the pure polymer. In particular, our analysis shows that for very small evaporation with delta ll exp(-3/(4gamma)) eps^3/4 skin formation can be prevented

2021, Peschka, Dirk, Heltai, Luca

We present a mathematical and numerical framework for the physical problem of thin-film fluid flows over planar surfaces including dynamic contact angles. In particular, we provide algorithmic details and an implementation of higher-order spatial and temporal discretisation of the underlying free boundary problem using the finite element method. The corresponding partial differential equation is based on a thermodynamic consistent energetic variational formulation of the problem using the free energy and viscous dissipation in the bulk, on the surface, and at the moving contact line. Model hierarchies for limits of strong and weak contact line dissipation are established, implemented and studied. We analyze the performance of the numerical algorithm and investigate the impact of the dynamic contact angle on the evolution of two benchmark problems: gravity-driven sliding droplets and the instability of a ridge.