【纤维粘弹性流体悬浮液动力学若干问题】.pdf

浙江大学硕士学位论文纤维粘弹性流体悬浮液动力学若干问题姓名:王立刚导师:范西俊 TavTa 万物皆流一九九四年十二月

浙江大学申请硬士学位论文纤维粘弹性流体悬浮液动力学若干问题 S1粘弹性流动的数值模拟对徽分型本构方程菜用Galerkin有限元方法计算得到的速度场产生失真振荡 [13,14]。这是由数值误差造成的。W数增大时,本构方程的对流项占主导地位,在流场中控制方程的类型会发生转变[15]:定常时由椭圆型转变为双曲型,非定常时由抛物型转变为双曲型[16]。在双曲型问题中构造具有高阶精度和较好稳定性的有限元算法是一项十分困难的工作,开始人们想到将应力沿实特征线(即流线)积分的流线有限元或边界元方法,将应力场作为拟体力代人守恒方程的下一次送代、这种方法对计算机内存要求高,难以

采用时间相关法和纤维Oldroyd-B流体半浓悬浮液的本构方程,对绕圆管内一球布、速度场、流函数和无量纲阻力系数。时间相关法的优点是不要求给定进口截面的纤维取向分布和应力张量分布,能方便地用于涡流区域的计算,积分方程所需的初值可任意给定。缺点是精度不够高,计算的W数只达到0。守恒方程是椭圆型的,但纤维结构张量的变化方程和高分子对偏应力张量贡献的本构方程是双曲型的。与其它粘弹性流动问题相比,粘弹性流体纤维悬浮液流动问题的控制方程是强双曲型的。而且流动中应力奇点的存在,它的作用尚不知道、纤维的存在使奇异问题更为严重,粘弹性流体纤维悬浮液的数值模拟比粘弹性流体更为因难 S1.