圆管湍流流速分布和带表面活性剂不可压缩两相流的数值模拟

来源 :湘潭大学 | 被引量 : 0次 | 上传用户:buffisher
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在过去的数十年里,尽管计算机的速度和内存都有很大提高,但很多实际问题的数值模拟仍需在很多简化下才可完成。也就是说,现有的计算机能力对付过多的自由度仍然有很多困难。在这种情形下,自适应算法应运而生,并在理论和实际应用上得到了广泛的重视。自适应网格方法是一种用来解决微分方程近似解的重要计算方法。把它应用在边界层和内层问题的近似解上是非常有效的。针对奇异摄动问题的自适应方法,我们做了一些基础性的研究。利用插值的思想,得到了基于导数误差的自适应移动有限元方法。针对椭圆型奇异摄动问题差分格式,在Bakhvalov-Shishkin网格下,我们通过构造网格函数并利用比较定理等证明了数值解一致一阶收敛于真解。  我们将自适应移动网格方法成功地应用于解决化工过程模拟与优化中涉及到一类流体计算问题。Luo等人提出了一个基于修正的Prandtl混合长概念的圆管湍流速度分布模型。这是一个控制问题,其数学模型是带边界层的非线性奇异摄动微分方程。我们通过迭代方法求解控制问题,利用自适应移动网格方法抓住奇异摄动微分方程的边界层,用Newton-Raphson方法求解非线性方程组,有效地求解了该问题。通过在流体速度,涡流粘性,摩擦因子等方面与实验结果的比较,表明这个新的模型比之前的一些模型更适合于描述圆管湍流的情况。  带表面活性剂的两相流是实验流体力学、理论流体力学和计算流体力学中的十分活跃且极具挑战性的重要课题。我们提出一种基于水平集方法求解带表面活性剂的两相流问题,在该方法中,我们用水平集函数来跟踪流体的交界面,把流体压力张量在界面处的Laplace-Young间断条件作为奇异界面力通过水平集函数写到流体的动量方程中。在计算方法中,我们首先用一种包含界面的小支集的光滑函数来近似表示奇异界面力的Delta函数,然后用一种二阶精度的投影方法求解正则化了的Navier-Stokes方程组,由此得到的速度场用来计算下一个时间步的水平集函数,以及活性剂的浓度。我们用高阶TVDRunge-Kutta方法和高阶WENO方法来求解水平集函数的运输方程,用一种二阶精度的Euler方法来求解活性剂在流体交界面上的运输扩散方程。我们的这一套方法简单易行,可直接推广到三维问题。我们针对Navier-Stokes剪切流中有一个带不可溶活性剂的小液滴的情况,进行了大量的数值模拟实验,我们的计算结果符合物理规律,特别是我们发现Marangoni力会阻碍小滴的拉伸。
其他文献
自然邻接插值是一种广泛应用于多元数据拟合的插值方法.自然元方法就是以自然邻接插值(形函数)作为试探函数和检验函数应用在Galerkin过程的数值计算方法,其依赖于Voronoi图
本文研究了非线性系统的自适应控制问题、状态观测器设计问题以及基于观测器的输出反馈镇定问题.主要研究成果如下:   首先,本文研究一类含有时变周期参数的非线性离散时间
空间逻辑是当前计算机理论、人工智能理论领域中正在兴起,而又具强大生命力的学科。在很短的时间内,国内外形成很多科研团体,各个国家投入大量的科研经费。新的文献如此之多,
作为非线性科学研究领域的一个重要课题,无穷维动力系统已成功地应用到许多重要方面,有着广泛的应用前景.尤其是,无穷维动力系统中的格点系统,受到了科学家的广泛关注,已被成
2011~2012年在12个省份采集129份烟叶样品,检测其物理指标和常规化学成分含量。数据统计结果显示,烤烟产地对烟叶的单叶重、含梗率、平衡含水率、抗张强度、叶片厚度、填充值
本文主要研究了n元集合上的拓扑个数T(n)和To拓扑个数To(n)的一些性质。我们首先证明了T(n)模2的周期性,即T(n)的奇偶性。利用б-代数和集合的划分之间的对应关系我们得到了T(
时滞微分方程理论研究的一个重要问题是研究解的定性性质,如解的稳定性、周期解、概周期解等.概周期解问题一直是微分方程理论的一个重要分支,在理论和实际应用中有着非常重
学位
奇异摄动问题出现在应用数学的许多分支里面。近年来,奇异摄动方程的解析解以及数值解的研究引起了许多学者的注意,并涌现了大量的这方面的文章。但是,这些文章大多都是讨论
近年来,对于分数积分方程和发展方程的研究获得了许多新的结果.但是,相对于整数阶微分方程而言,分数阶方程在理论研究方面还很不完善,有许多领域尚未涉及,需要进一步研究.  
顾客满意度是指顾客对其要求已被满足的程度的感觉,它是企业服务工作的目标,也是企业成功的动力,几乎可以说是左右企业生存发展重要因素之一。顾客满意度作为顾客满意的量化