Computing dynamics of thin films via large scale GPU-based simulations
Computing dynamics of thin films via large scale GPU-based simulations
复制标题
通过基于 GPU 的大规模模拟计算薄膜动力学
DOI：
10.1016/j.jcpx.2018.100001
复制
发表时间：
2019
影响因子：
--
通讯作者：
Kondic, Lou
中科院分区：
文献类型：
--
作者：
Lam, Michael-Angelo Y.-H.;Cummings, Linda J.;Kondic, Lou
我们给出了四阶非线性扩散型偏微分方程的大规模模拟结果，这些方程通常在基材上的薄流体膜动力学建模时遇到。仿真基于交替方向隐式（ADI）方法，主要部分计算工作在GPU计算环境下进行。高效和准确的计算允许在三维空间（3D）和长计算时间的大型计算域上进行模拟。我们将所开发的方法应用于纳米级厚度的流体薄膜的不稳定性问题。模拟的大尺度最大限度地减少了边界的影响，也允许模拟在已发表的实验中遇到的大小域。因此，我们可以以前所未有的详细程度分析不稳定性的发展。一个特别的重点是分析不稳定性发展的方式，特别是关于线性不稳定膜的spinodal和成核类型的脱湿的差异，以及亚稳膜的不稳定性。3D模拟允许考虑一些最近的结果，这些结果是以前在二维几何中获得的。一些新的结果包括使用傅里叶变换和拓扑不变量（贝蒂数）来区分spinodal和成核类型的不稳定性的结果，用精确的术语描述导致卫星滴形成的复杂过程，以及区分线性不稳定和亚稳定状态下不断发展的膜锋的形状。我们还讨论了向列液晶和聚合物薄膜的模拟结果与现有实验结果的直接比较。
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.
登录
查看更多内容
DOI：
10.1016/j.jcp.2015.03.051
发表时间：
2014-09
期刊：
J. Comput. Phys.
影响因子：
--
作者：
K. Mahady;S. Afkhami;L. Kondic
通讯作者：
K. Mahady;S. Afkhami;L. Kondic
直接计算流动（包括流固相互作用）的数值方法：模拟接触角、薄膜破裂和去湿
DOI：
--
发表时间：
2015
期刊：
影响因子：
--
作者：
K. Mahady;S. Afkhami;L. Kondic
通讯作者：
L. Kondic
影响因子：
4.1
作者：
Seric, Ivana;Afkhami, Shahriar;Kondic, Lou
通讯作者：
Kondic, Lou
向列液晶在倾斜基板上的三维涂覆流动†
影响因子：
1.9
作者：
M. Lam;L. Cummings;T. Lin;L. Kondic
通讯作者：
L. Kondic
nCB 液晶薄膜中成核引起的波动不稳定性
DOI：
--
发表时间：
2001
期刊：
影响因子：
--
作者：
S. Schlagowski;K. Jacobs;S. Herminghaus
通讯作者：
S. Herminghaus