We propose an adaptive threshold dynamics method for wetting problems in three space dimensions.The method is based on solving a linear heat equation and a thresholding step in each iteration.The heat equation is disc...We propose an adaptive threshold dynamics method for wetting problems in three space dimensions.The method is based on solving a linear heat equation and a thresholding step in each iteration.The heat equation is discretized by a cell-centered finite volume method on an adaptively refined mesh.An efficient technique for volume conservation is developed on the nonuniform meshes based on a quick-sorting operation.By this method,we compute some interesting wetting problems on complicated surfaces.Numerical results verify some recent theory for the apparent contact angle on rough and chemically patterned surfaces.展开更多
The Merriman-Bence-Osher threshold dynamics method is an efficient algorithm to simulate the motion by mean curvature.It has the advantages of being easy to implement and with high efficiency.In this paper,we propose ...The Merriman-Bence-Osher threshold dynamics method is an efficient algorithm to simulate the motion by mean curvature.It has the advantages of being easy to implement and with high efficiency.In this paper,we propose a threshold dynamics method for dislocation dynamics in a slip plane,in which the spatial operator is essentially an anisotropic fractional Laplacian.We show that this threshold dislocation dynamics method is able to give two correct leading orders in dislocation velocity,including both the O(logε)local curvature force and the O(1)nonlocal force due to the long-range stress field generated by the dislocations as well as the force due to the applied stress,whereεis the dislocation core size,if the time step is set to beΔt=ε.This generalizes the available result of threshold dynamics with the corresponding fractional Laplacian,which is on the leading order O(logΔt)local curvature velocity under the isotropic kernel.We also propose a numerical method based on spatial variable stretching to correct the mobility and to rescale the velocity for efficient and accurate simulations,which can be applied generally to any threshold dynamics method.We validate the proposed threshold dislocation dynamics method by numerical simulations of various motions and interaction of dislocations.展开更多
基金supported in part by NSFC grants DMS-11971469,DMS-11771290the National Key R&D Program of China under Grant 2018YFB0704304,Grant 2018YFB0704300.
摘要We propose an adaptive threshold dynamics method for wetting problems in three space dimensions.The method is based on solving a linear heat equation and a thresholding step in each iteration.The heat equation is discretized by a cell-centered finite volume method on an adaptively refined mesh.An efficient technique for volume conservation is developed on the nonuniform meshes based on a quick-sorting operation.By this method,we compute some interesting wetting problems on complicated surfaces.Numerical results verify some recent theory for the apparent contact angle on rough and chemically patterned surfaces.
基金supported by the Hong Kong Research Grants Council Collaborative Research Fund C1005-19Gthe Project of Hetao Shenzhen-HKUST Innovation Cooperation Zone HZQB-KCZYB-2020083supported by Shenzhen Fund 2021 Basic Research General Programme(project code:JCYJ20210324115400002).
摘要The Merriman-Bence-Osher threshold dynamics method is an efficient algorithm to simulate the motion by mean curvature.It has the advantages of being easy to implement and with high efficiency.In this paper,we propose a threshold dynamics method for dislocation dynamics in a slip plane,in which the spatial operator is essentially an anisotropic fractional Laplacian.We show that this threshold dislocation dynamics method is able to give two correct leading orders in dislocation velocity,including both the O(logε)local curvature force and the O(1)nonlocal force due to the long-range stress field generated by the dislocations as well as the force due to the applied stress,whereεis the dislocation core size,if the time step is set to beΔt=ε.This generalizes the available result of threshold dynamics with the corresponding fractional Laplacian,which is on the leading order O(logΔt)local curvature velocity under the isotropic kernel.We also propose a numerical method based on spatial variable stretching to correct the mobility and to rescale the velocity for efficient and accurate simulations,which can be applied generally to any threshold dynamics method.We validate the proposed threshold dislocation dynamics method by numerical simulations of various motions and interaction of dislocations.