Computing dynamics of thin films via large scale GPU-based simulations
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...
Main Authors: | , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Elsevier
2019-03-01
|
Series: | Journal of Computational Physics: X |
Online Access: | http://www.sciencedirect.com/science/article/pii/S2590055218300015 |
id |
doaj-1f926efa789b4082b2f2378c58f8aa4e |
---|---|
record_format |
Article |
spelling |
doaj-1f926efa789b4082b2f2378c58f8aa4e2020-11-25T00:46:47ZengElsevierJournal of Computational Physics: X2590-05522019-03-012Computing dynamics of thin films via large scale GPU-based simulationsMichael-Angelo Y.-H. Lam0Linda J. Cummings1Lou Kondic2Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, NJ, 07102, USADepartment of Mathematical Sciences, New Jersey Institute of Technology, Newark, NJ, 07102, USACorresponding author.; Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, NJ, 07102, USAWe 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. Keywords: Thin films, Long-wave approximation, Finite difference simulations, GPU computing, Film instabilitieshttp://www.sciencedirect.com/science/article/pii/S2590055218300015 |
collection |
DOAJ |
language |
English |
format |
Article |
sources |
DOAJ |
author |
Michael-Angelo Y.-H. Lam Linda J. Cummings Lou Kondic |
spellingShingle |
Michael-Angelo Y.-H. Lam Linda J. Cummings Lou Kondic Computing dynamics of thin films via large scale GPU-based simulations Journal of Computational Physics: X |
author_facet |
Michael-Angelo Y.-H. Lam Linda J. Cummings Lou Kondic |
author_sort |
Michael-Angelo Y.-H. Lam |
title |
Computing dynamics of thin films via large scale GPU-based simulations |
title_short |
Computing dynamics of thin films via large scale GPU-based simulations |
title_full |
Computing dynamics of thin films via large scale GPU-based simulations |
title_fullStr |
Computing dynamics of thin films via large scale GPU-based simulations |
title_full_unstemmed |
Computing dynamics of thin films via large scale GPU-based simulations |
title_sort |
computing dynamics of thin films via large scale gpu-based simulations |
publisher |
Elsevier |
series |
Journal of Computational Physics: X |
issn |
2590-0552 |
publishDate |
2019-03-01 |
description |
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. Keywords: Thin films, Long-wave approximation, Finite difference simulations, GPU computing, Film instabilities |
url |
http://www.sciencedirect.com/science/article/pii/S2590055218300015 |
work_keys_str_mv |
AT michaelangeloyhlam computingdynamicsofthinfilmsvialargescalegpubasedsimulations AT lindajcummings computingdynamicsofthinfilmsvialargescalegpubasedsimulations AT loukondic computingdynamicsofthinfilmsvialargescalegpubasedsimulations |
_version_ |
1725263196097347584 |