论文部分内容阅读
随着塑料和复合材料加工工业的发展以及消费者对成品质量要求的提高,需要对加工中聚合物流体等的流动规律和微观结构进行深入的研究。而采用分子动力学模型及模拟方法可以对分子构形进行粗略的描述,并可以模拟出流动诱导过程中分子构形的一些重要特性。流变学的研究方法有两种途径:宏观唯象法和微观分子理论。分子论方法是基于分子水平的非牛顿流体流动的模拟手段,由统计学的精神得知,要得到有意义的统计量(如流体的各种应力)必须模拟足够数量的分子与足够长的时间,这对电脑的运算能力产生很大的挑战。因此比较实际可行的方法是选择适当的分子分子模拟方法,不但能掌握主要的分子动态特性,又能大幅度的节省电脑的运算时间。CONNFFESSIT方法(布朗动力学)是基于分子论方法,直接利用微结构模型对聚合物流动过程进行动力学模拟。它有着不需要封闭的本构方程,容易在不同聚合物模型中转换等优点,在聚合物模拟领域有着一定的发展前景。本文把方差缩减方法运用到CONNFFESSIT方法中,模拟了库埃特流等算例,并对结果进行比较分析。本文的主要工作如下:(1)系统的分析了这种研究方法,从布朗动力学(CONNFFESSIT算法)的基础出发,建立库埃特流的基本控制方程,采用哑铃分子模型,通过布朗动力学方法,推导出库埃特流动下聚合物分子构形的控制方程。(2)在聚合物分子构形的控制方程的基础上,用系综平均的方法得到其应力方程,并结合宏观介质力学的动量方程,分析方程的求解,编写计算机模拟程序,得到库埃特流的应力参数和速度参数。(3)把CONNFFESSIT算法得到的结果与采用传统的Oldroyd-B模型的结果进行比较后,进一步对这种算法采用构形场法和控制变量法两种进行方差缩减,推导出方差缩减后CONNFFESSIT算法的控制方程,分析方程的求解,得到模拟结果。发现这种算法计算结果更加稳定,标准差小,精度高,并且所得到的模拟结果与采用传统的Oldroyd-B模型的结果非常吻合。