论文部分内容阅读
自然界中广泛存在着非平衡态现象,如何对这些现象建立满足热力学定律的模型一直是一个重要的问题。除了满足热力学定律外,我们还希望建立的模型要尽可能地满足其它物理定律,比如体积守恒、质量守恒、原子守恒等。近年来,相场成为描述系统从非平衡态向平衡态动力学演化的一个重要方法,我们将满足热力学定律及其他物理定律的相场模型称为热力学一致相场模型。相变量可能代表系统内某组分的体积分数、质量分数、原子密度等物理量,如何在模型的构建中保证这些物理量在整个区域上的守恒性成为了建模的重点。Allen-Cahn模型及Cahn-Hillard模型是建模过程中比较常用的两种模型,Allen-Cahn模型不能保证相变量所描述的物理量在整个区域上积分后守恒(如体积守恒、质量守恒及原子数守恒等),而Cahn-Hillard方程却可以严格地保证这一守恒性。但是Cahn-Hillard方程中空间算子的阶数要高于Allen-Cahn方程中空间算子的阶数,这导致数值求解Cahn-Hilard方程需要的时间成本要高于求解Allen-Cahn方程。在本文中,通过在经典的Allen-Cahn模型后添加非局部约束以使得修正后的Allen-Cahn模型满足某种守恒性,我们称之为带非局部约束的Allen-Cahn模型。我们设计了两种带非局部约束的Allen-Cahn模型,可以证明这两种带非局部约束的Allen-Cahn模型都是能量耗散的,并且能够保证守恒性。我们对这一模型在近平衡态的动力学进行了线性稳定性分析,并比较了它与Allen-Cahn模型及Cahn-Hillard模型失稳条件的异同。在算法设计部分,通过使用能量二次化方法(EQ)及标量辅助变量(SAV)方法对带非局部约束的Allen-Cahn模型设计了一系列线性的、二阶精度的、能量稳定的时间半离散的数值格式;之后,我们用差分方法进行空间离散以得到一个全离散的数值格式。我们证明了所有的时间半离散及全离散数值格式都满足无条件能量稳定性及解的存在唯一性。通过对这些数值格式的实现,我们对不互溶两相流问题,晶体生长问题以及聚合物溶液的相分离问题进行了数值模拟,并且设计了多个算例对数值格式进行了评估。针对大步长模拟,我们提供了一些提高精度的技巧。本文的另一个重点是探究溶液中化学反应对输运过程的影响。我们从热力学基本定律出发对多组分化学反应系统中的输运过程及化学反应过程进行了刻画,可以证明,在忽略界面效应、处于均匀状态的系统中,系统的Helmholtz自由能是耗散的;而在非均匀的化学反应系统中,在适当的边界条件及假设下,也可以保证相应的自由能是耗散的。应用广义Onsager原理,反应-输运方程可以写作为一个Allen-Cahn形式的方程和一个Cahn-Hillard形式方程的加和。最后我们以多组分的聚合物溶液为例通过数值模拟探究了化学反应对输运过程的影响。