TY - JOUR
T1 - Computing dynamics of thin films via large scale GPU-based simulations
AU - Lam, Michael Angelo Y.H.
AU - Cummings, Linda J.
AU - Kondic, Lou
N1 - Funding Information:
Research supported by the National Science Foundation under grant DMS-1211713.Research supported by the National Science Foundation under grant DMS-1815613.Research supported by the National Science Foundation under grant CBET-1604351.Research supported by the National Aeronautics and Space Administration under grant No. NNX16AQ79G.
Publisher Copyright:
© 2018 The Author(s)
PY - 2019/3
Y1 - 2019/3
N2 - We present the results of large scale simulations of 4th order nonlinear partial differential equations of diffusion 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 computational 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 [28]. 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.
AB - We present the results of large scale simulations of 4th order nonlinear partial differential equations of diffusion 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 computational 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 [28]. 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.
KW - Film instabilities
KW - Finite difference simulations
KW - GPU computing
KW - Long-wave approximation
KW - Thin films
UR - http://www.scopus.com/inward/record.url?scp=85064808338&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85064808338&partnerID=8YFLogxK
U2 - 10.1016/j.jcpx.2018.100001
DO - 10.1016/j.jcpx.2018.100001
M3 - Article
AN - SCOPUS:85064808338
SN - 2590-0552
VL - 2
JO - Journal of Computational Physics: X
JF - Journal of Computational Physics: X
M1 - 100001
ER -