多尺度模拟方法概述-计算传热学作业.pdfVIP

《计算传热学》学期作业 多尺度模拟方法概述 摘要：本文简单介绍多尺度模拟的思想，应用及存在的问题。 关键词：数值模拟；多尺度模拟 世界的本质是多尺度的，在不同的尺度下物质表现出不同的特征。如流体在 分子尺度下表现为离散的不确定的粒子，而在宏观尺度下表现为连续的确定性的 介质。在不同的时间和空间尺度下由于其尺度特性的不同，往往所采用的方法也 不同,如图1[1]所示。 图1 各种空间时间尺度下适用的模拟方法 文献[2]利用Kn数来鉴定何种特征尺度下流体流动适合用何种方法。Kn数的 物理意义是分子平均自由程与特征长度的比值。 -3 Kn10 ，流动符合连续介质假设，可用N-S方程； -3 -1 10 Kn10 ，边界是滑移边界，速度和温度有跳跃，控制方程为N-S方程； -1 10 Kn10 ，过渡流动，N-S方程不再适用，可用格子Boltzmann方法； Kn10 ，分子流动，可用分子动力学模拟方法。 模拟方法大致可分为宏观方法，介观方法，微观方法。宏观方法即流动符合 连续介质假设，传热的空间尺度和时间尺度符合傅立叶导热定律；微观方法是从 分子运动碰撞理论来建立方程；介观方法是介于微观方法和宏观方法之间。这三 种方法各有优缺点。宏观方法不能揭示微观的物理现象，但是方法成熟，应用方 便。微观或介观方法更适合描述极端尺度的物理现象，但是计算量巨大，方法不 成熟，工程应用极少。如果在采用宏观方法的过程中，可将微观尺度的信息带入， 建立一种微观——宏观耦合的多尺度模拟方法可以结合两者的优点，又可以削弱 两者的缺点。 [3] 多尺度问题表现 为: 已知一个模型的宏观描述, 但这种宏观描述在某些局 部区域失效, 必须要用低尺度微观非线性描述代替。模型的微观特性既受制于宏 观上的作用因素, 又可能显著影响宏观性能。但微观结构, 性能与状态何时、以 怎样的途径去影响宏观性能并不清楚。 假定一个给定系统的微观行为可以使用微观模型变量u表示, 系统的宏观行 为用宏观模型变量U表示, 那么宏观模型变量U与微观模型变量u可以通过压缩乘 子Q或者重构算子R联系起来： U=Qu RU=u 多尺度模拟的难度在于两种尺度的耦合，即如何建模。建模的策略有两种[4- 6] ：一种策略是先在较低的尺度上建模, 然后将结果放入高尺度模型中, 这是一 个从小尺度到大尺度的递阶过程。但低尺度建模的理论是一个重要问题。采用这 种策略的方法一般称作信息传递的多尺度方法或递阶的多尺度方法另一种策略 是在不同尺度上同时建模, 将区域分成不同尺度定律控制的区域, 这些区域可以 重叠也可以不重叠,在交界处实现连接。在这种策略中, 区域之间的连接也是一 个重要问题采用这种策略的方法一般称作并发(一致) 的多尺度方法。 国内外许多学着都致力于开发多尺度模拟方法，主要是介观宏观耦合和微观 宏观耦合。多尺度模拟可用于分析材料、化学、能源工程等领域的问题，特别是 微小装置的结构、流动和传热问题。随着微纳米科学技术的发展诞生出一个新的 技术领域，微/纳机电系统(Micro/Nano ElectroMechanical System，M/NEMS) 。微 机电系统在工业、通信、环境、生物、医疗和航空航天等领域有着十分广阔的应 用前景。 对于 M/NEMS 尺度来说，分子动力学模拟虽可提供原子尺度信息，但只能 考虑几百万个原子，处理的规模太小；而连续介质力学模拟不能提供接触区域（通 常只有几层原子）微观结构的变化；因而不利于人们全面地揭示微/纳尺度下各 种现象的相关性。多尺度模拟在一个系统的不同区域内采用不同的模型。例如， 在发生较大变形的区域采用量子力学或分子动力学模型，在Kn数较大的区域采 用分子动力学模拟或格子Boltzmann方法，以获得该区域的原子尺度信息；在变 形较小、Kn数较小的区域采用连续介质模型，以减