The optimized geometries at the RHF/6-311++G** level, the relatively stable energy at the MPW1PW91/6-311++G** level and the structural characters of anions have been acquired, indicating the stability is related to ...The optimized geometries at the RHF/6-311++G** level, the relatively stable energy at the MPW1PW91/6-311++G** level and the structural characters of anions have been acquired, indicating the stability is related to the chemical bonding of μ2?P atoms and the distri- bution of negative charges. The configurations of cage units P8 4- and P9 5- are stable due to the less torsion, but their ES values are relatively higher than that of P7 3- with more μ2?P atoms and the isolated stability is lower than that of P7 . They potentially play an important role as intermediate 3- in chemical reaction of producing complicated polyphosphides. Based on the related electronic properties, a stable polyanion must have low valence electron concentration, no (μ2?P)?(μ2?P) bond and a little dispersive charge. The earmark IR frequencies of cage units have been assigned to the vibration models in the end.展开更多
In this paper, we present an efficient energy stable scheme to solve a phase field model incorporating contact line condition. Instead of the usually used Cahn-Hilliard type phase equation, we adopt the Allen-Cahn typ...In this paper, we present an efficient energy stable scheme to solve a phase field model incorporating contact line condition. Instead of the usually used Cahn-Hilliard type phase equation, we adopt the Allen-Cahn type phase field model with the static contact line boundary condition that coupled with incompressible Navier-Stokes equations with Navier boundary condition. The projection method is used to deal with the Navier-Stokes equa- tions and an auxiliary function is introduced for the non-convex Ginzburg-Landau bulk potential. We show that the scheme is linear, decoupled and energy stable. Moreover, we prove that fully discrete scheme is also energy stable. An efficient finite element spatial discretization method is implemented to verify the accuracy and efficiency of proposed schemes. Numerical results show that the proposed scheme is very efficient and accurate.展开更多
Here we focus on the numerical simulation of the phase separation about macromolecule microsphere composite (MMC) hydrogel. The model is based on time-dependent Ginzburg- Landau (TDGL) equation with the reticular ...Here we focus on the numerical simulation of the phase separation about macromolecule microsphere composite (MMC) hydrogel. The model is based on time-dependent Ginzburg- Landau (TDGL) equation with the reticular free energy. An unconditionally energy stable difference scheme is proposed based on the convex splitting of the corresponding energy functional. In the numerical experiments, we observe that simulating the whole process of the phase separation requires a considerably long time. We also notice that the total free energy changes significantly in initial stage and varies slightly in the following time. Based on these properties, we apply the adaptive time stepping strategy to improve the computational efficiency. It is found that the application of time step adaptivity can not only resolve the dynamical changes of the solution accurately but also significantly save CPU time for the long time simulation.展开更多
This paper studies the numerical simulations for the Cahn-Hilliard equation which describes a phase separation phenomenon.The numerical simulation of the Cahn-Hilliardmodel needs very long time to reach the steady sta...This paper studies the numerical simulations for the Cahn-Hilliard equation which describes a phase separation phenomenon.The numerical simulation of the Cahn-Hilliardmodel needs very long time to reach the steady state,and therefore large time-stepping methods become useful.The main objective of this work is to construct the unconditionally energy stable finite difference scheme so that the large time steps can be used in the numerical simulations.The equation is discretized by the central difference scheme in space and fully implicit second-order scheme in time.The proposed scheme is proved to be unconditionally energy stable and mass-conservative.An error estimate for the numerical solution is also obtained with second order in both space and time.By using this energy stable scheme,an adaptive time-stepping strategy is proposed,which selects time steps adaptively based on the variation of the free energy against time.The numerical experiments are presented to demonstrate the effectiveness of the adaptive time-stepping approach.展开更多
In this paper,we will develop a first order and a second order convex splitting,and a first order linear energy stable fully discrete local discontinuous Galerkin(LDG)methods for the modified phase field crystal(MPFC)...In this paper,we will develop a first order and a second order convex splitting,and a first order linear energy stable fully discrete local discontinuous Galerkin(LDG)methods for the modified phase field crystal(MPFC)equation.In which,the first order linear scheme is based on the invariant energy quadratization approach.The MPFC equation is a damped wave equation,and to preserve an energy stability,it is necessary to introduce a pseudo energy,which all increase the difficulty of constructing numerical methods comparing with the phase field crystal(PFC)equation.Due to the severe time step restriction of explicit timemarchingmethods,we introduce the first order and second order semi-implicit schemes,which are proved to be unconditionally energy stable.In order to improve the temporal accuracy,the semi-implicit spectral deferred correction(SDC)method combining with the first order convex splitting scheme is employed.Numerical simulations of the MPFC equation always need long time to reach steady state,and then adaptive time-stepping method is necessary and of paramount importance.The schemes at the implicit time level are linear or nonlinear and we solve them by multigrid solver.Numerical experiments of the accuracy and long time simulations are presented demonstrating the capability and efficiency of the proposed methods,and the effectiveness of the adaptive time-stepping strategy.展开更多
The main purpose of this work is to contrast and analyze a large timestepping numerical method for the Swift-Hohenberg(SH)equation.This model requires very large time simulation to reach steady state,so developing a l...The main purpose of this work is to contrast and analyze a large timestepping numerical method for the Swift-Hohenberg(SH)equation.This model requires very large time simulation to reach steady state,so developing a large time step algorithm becomes necessary to improve the computational efficiency.In this paper,a semi-implicit Euler schemes in time is adopted.An extra artificial term is added to the discretized system in order to preserve the energy stability unconditionally.The stability property is proved rigorously based on an energy approach.Numerical experiments are used to demonstrate the effectiveness of the large time-stepping approaches by comparing with the classical scheme.展开更多
In this paper,based on the Lagrange Multiplier approach in time and the Fourierspectral scheme for space,we propose efficient numerical algorithms to solve the phase field crystal equation.The numerical schemes are u...In this paper,based on the Lagrange Multiplier approach in time and the Fourierspectral scheme for space,we propose efficient numerical algorithms to solve the phase field crystal equation.The numerical schemes are unconditionally energy stable based on the original energy and do not need the lower bound hypothesis of the nonlinear free energy potential.The unconditional energy stability of the three semi-discrete schemes is proven.Several numerical simulations in 2D and 3D are demonstrated to verify the accuracy and efficiency of our proposed schemes.展开更多
文摘The optimized geometries at the RHF/6-311++G** level, the relatively stable energy at the MPW1PW91/6-311++G** level and the structural characters of anions have been acquired, indicating the stability is related to the chemical bonding of μ2?P atoms and the distri- bution of negative charges. The configurations of cage units P8 4- and P9 5- are stable due to the less torsion, but their ES values are relatively higher than that of P7 3- with more μ2?P atoms and the isolated stability is lower than that of P7 . They potentially play an important role as intermediate 3- in chemical reaction of producing complicated polyphosphides. Based on the related electronic properties, a stable polyanion must have low valence electron concentration, no (μ2?P)?(μ2?P) bond and a little dispersive charge. The earmark IR frequencies of cage units have been assigned to the vibration models in the end.
基金R. Chen is partially supported by the Fundamental Research Funds for Central Universities 24820182018RC25-500418780 and by the China Postdoctoral Science Foundation grant No. 2016M591122. X. Yang is partially supported by NSF DMS-1200487, NSF DMS-1418898, AFOSR FA9550-12-1-0178. H. Zhang is partially supported by NSFC/RGC Joint Research Scheme No. 11261160486, NSFC grant No. 11471046, 11571045.
文摘In this paper, we present an efficient energy stable scheme to solve a phase field model incorporating contact line condition. Instead of the usually used Cahn-Hilliard type phase equation, we adopt the Allen-Cahn type phase field model with the static contact line boundary condition that coupled with incompressible Navier-Stokes equations with Navier boundary condition. The projection method is used to deal with the Navier-Stokes equa- tions and an auxiliary function is introduced for the non-convex Ginzburg-Landau bulk potential. We show that the scheme is linear, decoupled and energy stable. Moreover, we prove that fully discrete scheme is also energy stable. An efficient finite element spatial discretization method is implemented to verify the accuracy and efficiency of proposed schemes. Numerical results show that the proposed scheme is very efficient and accurate.
文摘Here we focus on the numerical simulation of the phase separation about macromolecule microsphere composite (MMC) hydrogel. The model is based on time-dependent Ginzburg- Landau (TDGL) equation with the reticular free energy. An unconditionally energy stable difference scheme is proposed based on the convex splitting of the corresponding energy functional. In the numerical experiments, we observe that simulating the whole process of the phase separation requires a considerably long time. We also notice that the total free energy changes significantly in initial stage and varies slightly in the following time. Based on these properties, we apply the adaptive time stepping strategy to improve the computational efficiency. It is found that the application of time step adaptivity can not only resolve the dynamical changes of the solution accurately but also significantly save CPU time for the long time simulation.
基金We would like to thank Prof.Houde Han of Tsinghua University and Prof.Qiang Du of Penn State University for their helpful discussions.Z.R.Zhang was supported by National NSF of China under Grant 10601007Z.H.Qiao was supported by the FRG grants of the Hong Kong Baptist University under Grant No.FRG2/09-10/034.
文摘This paper studies the numerical simulations for the Cahn-Hilliard equation which describes a phase separation phenomenon.The numerical simulation of the Cahn-Hilliardmodel needs very long time to reach the steady state,and therefore large time-stepping methods become useful.The main objective of this work is to construct the unconditionally energy stable finite difference scheme so that the large time steps can be used in the numerical simulations.The equation is discretized by the central difference scheme in space and fully implicit second-order scheme in time.The proposed scheme is proved to be unconditionally energy stable and mass-conservative.An error estimate for the numerical solution is also obtained with second order in both space and time.By using this energy stable scheme,an adaptive time-stepping strategy is proposed,which selects time steps adaptively based on the variation of the free energy against time.The numerical experiments are presented to demonstrate the effectiveness of the adaptive time-stepping approach.
基金Research of R.Guo is supported by NSFC grant No.11601490Research of Y.Xu is supported by NSFC grant No.11371342,11626253,91630207.
文摘In this paper,we will develop a first order and a second order convex splitting,and a first order linear energy stable fully discrete local discontinuous Galerkin(LDG)methods for the modified phase field crystal(MPFC)equation.In which,the first order linear scheme is based on the invariant energy quadratization approach.The MPFC equation is a damped wave equation,and to preserve an energy stability,it is necessary to introduce a pseudo energy,which all increase the difficulty of constructing numerical methods comparing with the phase field crystal(PFC)equation.Due to the severe time step restriction of explicit timemarchingmethods,we introduce the first order and second order semi-implicit schemes,which are proved to be unconditionally energy stable.In order to improve the temporal accuracy,the semi-implicit spectral deferred correction(SDC)method combining with the first order convex splitting scheme is employed.Numerical simulations of the MPFC equation always need long time to reach steady state,and then adaptive time-stepping method is necessary and of paramount importance.The schemes at the implicit time level are linear or nonlinear and we solve them by multigrid solver.Numerical experiments of the accuracy and long time simulations are presented demonstrating the capability and efficiency of the proposed methods,and the effectiveness of the adaptive time-stepping strategy.
基金supported by the Fundamental Research Funds for the CentralUniversities andNationalNSF of China under grantNos.11271048,1130021 and 11571054.
文摘The main purpose of this work is to contrast and analyze a large timestepping numerical method for the Swift-Hohenberg(SH)equation.This model requires very large time simulation to reach steady state,so developing a large time step algorithm becomes necessary to improve the computational efficiency.In this paper,a semi-implicit Euler schemes in time is adopted.An extra artificial term is added to the discretized system in order to preserve the energy stability unconditionally.The stability property is proved rigorously based on an energy approach.Numerical experiments are used to demonstrate the effectiveness of the large time-stepping approaches by comparing with the classical scheme.
基金The work of Q.Zhuang is supported by the National Natural Science Foundation of China(No.11771083)The research of S.Zhai is supported in part by the Natural Science Foundation of China(No.11701196)+3 种基金the Natural Science Foundation of Fujian Province(No.2020J01074)The work of Z.Weng is supported in part by the Natural Science Foundation of China(No.11701197)Supported by the Fundamental Research Funds for the Central Universities(No.ZQN-702)the Key Laboratory of Intelligent Computing and Information Processing of Ministry of Education(Xiangtan University)(No.2020ICIP03).
文摘In this paper,based on the Lagrange Multiplier approach in time and the Fourierspectral scheme for space,we propose efficient numerical algorithms to solve the phase field crystal equation.The numerical schemes are unconditionally energy stable based on the original energy and do not need the lower bound hypothesis of the nonlinear free energy potential.The unconditional energy stability of the three semi-discrete schemes is proven.Several numerical simulations in 2D and 3D are demonstrated to verify the accuracy and efficiency of our proposed schemes.