Computing dynamics of thin films via large scale GPU-based simulations

DOI10.1016/J.JCPX.2018.100001MaRDI QIDQ6186179zbMATH OpenOpenAlexWikidataFDO

Authors M. A. Lam, Linda Cummings, L. Kondic

Publication date 9 January 2024

Published in Journal of Computational Physics: X (Search for Journal in Brave)

Copyright license Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International

Full work available at URL https://arxiv.org/abs/1806.10550

zbMATH Keywords

thin filmsGPU computinglong-wave approximationfinite difference simulationsfilm instabilities

Mathematics Subject Classification ID

Incompressible viscous fluids (76Dxx) Basic methods in fluid mechanics (76Mxx) Foundations, constitutive equations, rheology, hydrodynamical models of non-fluid phenomena (76Axx)

Abstract: We present the results of large scale simulations of 4th order nonlinear partial differential equations of dif- fusion type that are typically encountered when modeling dynamics of thin fluid films on substrates. The simulations are based on the alternate direction implicit (ADI) method, with the main part of the compu- tational work carried out in the GPU computing environment. Efficient and accurate computations allow for simulations on large computational domains in three spatial dimensions (3D) and for long computational times. We apply the methods developed to the particular problem of instabilities of thin fluid films of nanoscale thickness. The large scale of the simulations minimizes the effects of boundaries, and also allows for simulating domains of the size encountered in published experiments. As an outcome, we can analyze the development of instabilities with an unprecedented level of detail. A particular focus is on analyzing the manner in which instability develops, in particular regarding differences between spinodal and nucleation types of dewetting for linearly unstable films, as well as instabilities of metastable films. Simulations in 3D allow for consideration of some recent results that were previously obtained in the 2D geometry (J. Fluid Mech. 841, 925 (2018)). Some of the new results include using Fourier transforms as well as topological invariants (Betti numbers) to distinguish the outcomes of spinodal and nucleation types of instabilities, describing in precise terms the complex processes that lead to the formation of satellite drops, as well as distinguishing the shape of the evolving film front in linearly unstable and metastable regimes. We also discuss direct comparison between simulations and available experimental results for nematic liquid crystal and polymer films.

Cites work

• scientific article; zbMATH DE number 852532 (Why is no real title available?)
• scientific article; zbMATH DE number 274399 (Why is no real title available?)
• 3D phase-field simulations of interfacial dynamics in Newtonian and viscoelastic fluids
• A parallel partition method for solving banded systems of linear equations
• A volume of fluid method for simulating fluid/fluid interfaces in contact with solid boundaries
• ADI schemes for higher-order nonlinear diffusion equations.
• Calculation of two-phase Navier-Stokes flows using phase-field modeling
• Computing three-dimensional thin film flows including contact lines.
• Contact-line dynamics of a diffuse fluid interface
• Direct numerical simulation of variable surface tension flows using a volume-of-fluid method
• Evolution of a thin film of nematic liquid crystal with anisotropic surface energy
• Front propagation into unstable states
• Long-wave theory of bounded two-layer films with a free liquid–liquid interface: Short- and long-time evolution
• MODELING SHAPES AND DYNAMICS OF CONFINED BUBBLES
• Modeling flow of nematic liquid crystal down an incline
• Modelling spreading dynamics of nematic liquid crystals in three spatial dimensions
• Numerical simulation of moving contact line problems using a volume-of-fluid method
• On evaporation of sessile drops with moving contact lines
• On the breakup of fluid films of finite and infinite extent
• Sharp-interface limit of the Cahn-Hilliard model for moving contact lines
• Some Aspects of the Cyclic Reduction Algorithm for Block Tridiagonal Linear Systems
• Some Computer Organizations and Their Effectiveness
• Stability of thin fluid films characterised by a complex form of effective disjoining pressure
• Thin films flowing down inverted substrates: Two dimensional flow
• Thin films flowing down inverted substrates: three-dimensional flow
• Three-dimensional coating flow of nematic liquid crystal on an inclined substrate

This page was built for publication: Computing dynamics of thin films via large scale GPU-based simulations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6186179)