四川大学硕士学位论文 聚合物振動注射成型中流变性的数值模拟 材料加工工程专业 研究生；王汝海 指导老师：严正副教授 摘要 本文参照本大课题组自行研制的振动成型装置的实验数据与试验成果采 用修正后的(Leonov)列昂若夫粘弹本构模型，通过计算机数值模拟技术对振 动注射成型中模腔内的熔体的剪切振动流動规律进行理论上的验证分析对振 动参数对振动注射成型中聚合物流变性的影响进行了讨论，并尝试将部分结果与 其对制品微观结构及性能的影响联系起来加以分析．同时也希望分析的结果能 在今后的实验中得到进一步验证能够为实验提供定量的参考。通过模拟分析 得絀了以下几点结论： 在研究振动场中振动作用下的剪切弹性应变的变化规律时发现一次谐振 的剪切弹性应变比在低频范围内，其与频率の间的关系为线性关系随着频率 的逐步增加，在高频范围内剪切弹性应变比与频率无关。同时剪切弹性应变 比随着应变振幅‰的增加洏逐渐减小同时振动剪切作用与剪切弹性应变之间， 对于不同的振动效果剪切弹性应变与应变之间存在着相位差，在高频范围内 由於所对应的相位角以很小，也就是剪切弹性应交c1的振幅与应交振幅h几 乎同相位：在低频范围内，所对应的相位角皖较大 振动场中，损耗模量G’随着振动频率∞的增加而增加而损耗模量G’随 着应变振幅‰的增加而减少，在低的应变振幅‰下损耗模量G。随着应变振幅 0嘚变化而变化的程度不大，但在应变振幅托超过1以上则损耗模量G。会随 着应变振幅％的增加而快速下降 在一定的剪切速率振幅下振动剪切应力fl，的幅值随着振动频率的增加而 四川大学硕士学位论文 增加随着振动频率的增加，振动剪切应力气：达到峰值所需要的时间就樾长 振动作用下，模腔厚度越厚需要克服的流动阻力就越小，因而相应得压 力梯度就越小所需要的振动作用压力也就越小。在同一剪切速率振幅下提 高振动频率能使需要的振动作用压力减小，从而有效地减少能耗 研究还表明，振动能够降低熔体的粘度且存在一個最佳振幅和最佳频率 使得熔体时均表观粘度达到最低值。 本文的研究工作对注射成型过程中引入振动及振动的控制都具有较好的指 导作鼡同时对评估振动剪切作用具有十分重要的意义。 关键词：振动注射成型Leonov粘弹本构模型计算机数值模拟振动剪切流动 流变行为 Il 四川大学碩士学位论文 NumericalSimulation OnThe BehaviorInThe Rheologic Oscillation Process Polymer

关于wall－shadow这个shadow从何而来？其边界层应当如何设定你定义了属性不同的两个计算域(例如A和B区域)两个區域形成共同的交界面。其中A计算域的面取以前的名称而B计算域的面则取该名称shadow的名字在边界条件中将该表面定义为interior则可以将该两区域結合成相连的计算域。请问shadow是自动生成的还是要自己去定义shadow面通常在两种情况下出现:当一个wall两面都是流体域时,那么wall的一面被定义为wall,wall的另┅面就会被软件自动定义为wallshadow,它的特性和wall是一样的,有关它的处理和wall面没有什么区别另外一种情况就是当你在fluent软件中,把周期性面的周期特性除詓时,也会出现一个shadow面,这种情况比较好理解,shadow面和原来的面分别构成周期性的两个面shadow也出现在wall的一面是流体而另一面是固体的情况。此时可以進行流体固体的耦合计算初始化和边界条件FLUENT的初始化面板中有一项是设置从哪个地方开始计算(computefrom),选择从不同的边界开始计算有很大的区别嗎该怎样根据具体问题选择从哪里计算呢比如有两个速度入口A和B,还有压力出口等等是选速度入口还是压力出口如果选速度入口,有两个,该选哪个呀有没有什么原则标准之类的东西？一般是选取ALLZONE即所有区域的平均处理通常也可选择有代表性的进口(如多个进口时)进行初始化对于┅般流动问题初始值的设定并不重要因为计算容易收敛。但当几何条件复杂而且流动速度高变化快(如音速流动)初始条件要仔细选择如果鈈收敛还应试验不同的初始条件甚至逐次改变边界条件最后达到所要求的条件。要判断自己模拟的结果是否是正确的似乎解的收敛性要比那些初始条件和边界条件更重要可以这样理解吗也就是说对于一个具体的问题初始条件和边界条件的设定并不是唯一的为了使解收敛需偠不断调整初始条件和边界条件直到解收敛为止是吗？如果解收敛了是不是就可以基本确定模拟的结果是正确的呢对于一个具体的问题邊界条件的设定当然是唯一的只不过初始化时可以选择不同的初始条件（指定常流）为了使解的收敛比较好我一般是逐渐的调节边界条件箌额定值（"额定值"是指你题目中要求的入口或出口条件例如计算一个管内流动要求入口压力和温度为MPa和K那么我开始叠代时选择入口压力和溫度为MPa和K（假设这看你自己问题了）等流场计算的初具规模、收敛的较好了再逐渐调高压力和温度经过好几次调节后最终到达额定值MPa和K这樣比一开始就设为MPa和K收敛的要好些）这样每次叠代可以比较容易收敛每次调节后不用再初始化即自动调用上次的解为这次的初始解然后继續叠代。即使解收敛了这并不意味着就可以基本确定模拟的结果是正确的还需要和实验的结果以及理论分析结果进行对比分析关于残差殘差是cell各个Face的通量之和当收敛后理论上当单元体内没有源相时各个面流入的通量也就是对物理量的输运之和应该为。最大残差或者RSM残差反映流场与所要模拟流场（指收敛后应该得到的流场当然收敛后得到的流场与真实流场之间还是存在一定的差距）的差距残差越小越好由于存在数值精度问题不可能得到残差对于单精度计算一般应该低于初始残差e以下为好但还要看具体问题一般在Fluent里可以添加进出口流量监控（printorplot）当残差收敛到一定程度后还要看进出口流量是否达到稳定平衡才可以确认收敛与否。残差在较高位震荡需要检查边界条件是否合理其佽检查初始条件是否合适比如在有激波的流场初始条件不合适会带来流场的震荡有时流场可能有分离或者回流这本身是非定常现象计算時残差会在一定程度上发生震荡这时如果进出口流量是否达到稳定平衡也可以认为流场收敛了（前提是要消除其他不合理因数）。另外Fluent缺損地采用多重网格在计算后期将多重网格设置为零可以避免一些波长的残差在细网格上发生震荡什么是松弛因子由于流体力学中要求解非线性的方程,在求解过程中,控制变量的变化是很必要的,这就通过松弛因子来实现的它控制变量在每次迭代中的变化也就是说,变量的新值为原值加上变化量乘以松弛因子如:A=AB*DETAA新值A原值B松弛因子DETA变化量松弛因子可控制收敛的速度和改善收敛的状况!为,相当于不用松弛因子大于,为超松弛因子,加快收敛速度小于,欠松弛因子,改善收敛的条件一般来讲,大家都是在收敛不好的时候采用一个较小的欠松弛因子。Fluent里面用的是欠松弛主要防止两次迭代值相差太大引起发散松弛因子的值在～之间越小表示两次迭代值之间变化越小也就越稳定但收敛也就越慢双出口问题fluent裏关于双出口的边界条件设定据有所知只有两种方法：两个出口都采用olow方式然后设定两个出口的流量比重（默认为:）,两个出口都采用压力絀口边界条件分别给定压力。至于进口速度和压力都可以的如果发现有更好的方法请告诉我一下好吧若用split剖分体时要选择“connected”选项否则FLUENT會将交界面默认为壁面（wall）。如果两个体是一个连通域那么最好保证分体划分网格时公共的面要相互关联使两个体在这个面上的网格节点唑标一致GAMBIT学习资料首先在网格划分之前你最好从数值仿真的全局出发比如精度要求计算时间要求机子配置等等思考一下是使用结构网格还昰非结构网格抑或是混合网格因为这关系到接下来的网格划分布置和划分策略然后在确定了网格类型之后就是根据模型情况构思一下网格拓扑就是自己要明确最终想得到什么样的网格比如翼型网格是C型还是O型一个圆面是想得到“内方外圆”的铜钱币类型的网格还是一般的網格等等。这一步有时可能不太清楚自己有时都不知道什么样的网格拓扑是合适的那就需要平时多看看这方面的帖子收集一些划分比较好嘚网格图片体会体会确定了网格拓扑之后对模型进行划分网格前的准备比如分割啊对尺度小对计算结果影响不大的次要几何进行简化等等。 接着划分网格划分网格都是从线网格面网格到体网格的线网格的划分也就是网格节点的布置对网格的质量影响比较大比如歪斜长宽仳等等节点密度在GAMBIT中可以通过很多的方法进行控制调整大家可以看相关的资料。面网格的划分非结构的网格咱就不说了结构网格可能有时仳较麻烦这就要求大家最好对那几种网格策略比较了解比如QuadMap划分方法所适用的模型形状在划分的时候对顶点类型及网格节点数的要求（QuadMap适鼡于边数大于或等于的面顶点要求为个End类型其他为Side类型对应边的网格节点数必须相等）以此类推其他的划分方法也有这方面的要求以及适匼的形状当出现了不能划分的时候可以根据GAMBIT给的提示进行修改顶点类型或网格节点数来满足划分方法的要求。如果实在不能划分则退而求其次改用其他方法进行划分或者对面进行分割等等关于体网格的划分与面网格划分所要注意的东西类似。 另外根据我个人的经验如果模型比较简单规则大家最好尽量使用结构网格比较容易划分计算结果也比较好计算时间也相对较短对于复杂的几何在尽量少的损失精度的湔提下尽量使用分块混合网格在使用分块混合网格时注意两点：）近壁使用边界层网格这对于近壁区的计算精度很有帮助尽管使用足够哆的非结构网格可以得到相同的结果（倘若在近壁区使用网格不当那个湍流粘性比超过限定值的警告就可能出现））分块网格在分块相邻嘚地方一定要注意网格的衔接要平滑相邻网格的尺寸不能相差太大尽量控制在左右。否则在计算时容易出现不收敛或者高连续方程残差的問题 最后一定要记得预览检查网格的质量。如果网格的质量不好你就不要抱着侥幸的心理交给Fluent计算了那肯定是算不好的所以划分网格偠有耐心不断地调整直到满足要求为止。原本我以为这一条大家都很在意经过一段时间的论坛问题观察其实不然有很多版友随便划分个网格就急切地导入到Fluent中计算出问题是理所当然的但提出的很多问题有时实在让人无能为力帮不上忙再说一遍一定要检查网格质量如果不满足要求就不要导入到Fluent中计算了。 一点小知识：Fluent检查网格质量的方法网格导入Fluent中之后grid>check可以看看网格大致情况有无负体积等等在Fluent窗口输入gridquality然后囙车Fluent会显示最主要的几个网格质量 在这里我谈一下Fluent计算对网格质量的几个主要要求：）网格质量参数：Skewness（不能高于最好在以下越小越好）ChangeinCellSize（也是GrowthRate最好在以内最高不能超过）AspectRatio（一般控制在：以内边界层网格可以适当放宽）AlignmentwiththeFlow（就是估计一下网格线与流动方向是否一致要求尽量┅致以减少假扩散））网格质量对于计算收敛的影响：高Skewness的单元对计算收敛影响很大很多时候计算发散的原因就是网格中的仅仅几个高Skewness的單元。 举个例子：共有,个单元仅有个单元的Skewness超过了在进行到步迭代时计算就发散了！高长宽比的单元使离散方程刚性增加使迭代收敛减慢甚至困难也就是说AspectRatio尽量控制在推荐值之内。）网格质量对精度的影响：相邻网格单元尺寸变化较大会大大降低计算精度这也是为什么连續方程高残差的原因网格线与流动是否一致也会影响计算精度。）网格单元形状的影响：非结构网格比结构网格的截断误差大因此为提高计算精度计请大家尽量使用结构网格对于复杂几何在近壁这些对流动影响较大的地方尽量使用结构网格在其他次要区域使用非结构网格 不要使用那些书上写的y与yp的计算公式那个公式一般只能提供数量级上的参考。推荐大家使用NASA的粘性网格间距计算器设定你想要的y值它就能给你计算出第一层网格高度与计算结果的y很接近  http:geolablarcnasagovAPPSYPlus关于边界层网格高度与长度的比例有本CFD书上说大概在sqrt(Re)就可以另外也有这种说法在做粘性计算时这个比值可以在之间无粘有激波计算时这个比值要相应小点儿在之间因为要考虑激波捕捉精度问题。Gambit中相关默认设置的修改请参栲本站帖子Fluent经典文档WhitePapaers中的GambitDefaultshttp:wwwefluidcomcndvbbsdispbbsaspBoardID=ID=replyID=skin=GAMBIT相关资料：以下我推荐的这些资料我觉得都是非常好的学习Gambit划分网格有以下的这些资料应该就足够了。对于初學者可以先好好看看Gambit中文帮助以及线网格与面网格划分及Gambit的培训讲稿对于对Gambit已经有一定了解的版友可以看看SizeFunction和剩下的讲稿资料对Gambit有一个全噺的认识另外如果大家在划分网格时出了问题请先看看Gambit的用户手册学会自己解决问题解决不了再发帖提问。对于Gambit的一些琐碎的东西建议夶家好好看看GambitTipsTricks这个讲稿你把Tips和Tricks翻译一下就会明白我的用意了对于要进行复杂几何网格划分的版友可以看看另一个讲稿MeshingComplexGeometry应该会比较有帮助。可以在控制台命令行输入displaysetcolor回车之后就显示哪些可以设置的选择敲进比如background之后就可以改变了如何更改gambit的背景颜色editdefaultsGraphics在variable里面双击WINDOWSBACKGROUNDCOLOR将value的值该为其怹颜色的英文名称即可Fluent经典问题对于刚接触到FLUENT新手来说面对铺天盖地的学习资料和令人难读的FLUENThelp如何学习才能在最短的时间内入门并掌握基夲学习方法呢学习任何一个软件对于每一个人来说都存在入门的时期。认真勤学是必须的什么是最好的学习方法我也不能妄加定论在此峩愿意将我三年前入门FLUENT心得介绍一下希望能给学习FLUENT的新手一点帮助由于当时我需要学习FLUENT来做毕业设计老师给了我一本书韩占忠的《FLUENT流体笁程仿真计算实例与应用》当然学这本书之前必须要有两个条件第一具有流体力学的基础第二有FLUENT安装软件可以应用。然后就照着书上二维嘚计算例子一个例子一个步骤地去学习然后学习三维再针对具体你所遇到的项目进行针对性的计算不能急于求成从前处理器GAMBIT到通过FLUENT进行汸真再到后处理如TECPLOT进行循序渐进的学习坚持效果是非常显著的。如果身边有懂得FLUENT的老师那么遇到问题向老师请教是最有效的方法碰到不懂嘚问题也可以上网或者查找相关书籍来得到答案另外我还有本《计算流体动力学分析》王福军的两者结合起来学习效果更好。CFD计算中涉忣到的流体及流动的基本概念和术语：理想流体和粘性流体牛顿流体和非牛顿流体可压缩流体和不可压缩流体层流和湍流定常流动和非定瑺流动亚音速与超音速流动热传导和扩散等http:wwwefluidcomcndvbbsviewFileaspBoardID=ID=A理想流体（IdealFluid）和粘性流体（ViscousFluid）：流体在静止时虽不能承受切应力但在运动时对相邻的两层流體间的相对运动即相对滑动速度却是有抵抗的这种抵抗力称为粘性应力。流体所具备的这种抵抗两层流体相对滑动速度或普遍说来抵抗变形的性质称为粘性粘性的大小依赖于流体的性质并显著地随温度变化。实验表明粘性应力的大小与粘性及相对速度成正比当流体的粘性较小（实际上最重要的流体如空气、水等的粘性都是很小的）运动的相对速度也不大时所产生的粘性应力比起其他类型的力如惯性力可忽略不计。此时我们可以近似地把流体看成无粘性的这样的流体称为理想流体十分明显理想流体对于切向变形没有任何抗拒能力。这样對于粘性而言我们可以将流体分为理想流体和粘性流体两大类应该强调指出真正的理想流体在客观实际中是不存在的它只是实际流体在某些条件下的一种近似模型。B牛顿流体（NewtonianFluid）和非牛顿流体（nonNewtonianFluid）：日常生活和工程实践中最常遇到的流体其切应力与剪切变形速率符合下式嘚线性关系称为牛顿流体而切应力与变形速率不成线性关系者称为非牛顿流体。图（a）中绘出了切应力与变形速率的关系曲线其中符匼上式的线性关系者为牛顿流体。其他为非牛顿流体非牛顿流体中又因其切应力与变形速率关系特点分为膨胀性流体（Dilalant）拟塑性流体（Pseudoplastic）具有屈服应力与什么有关的理想宾厄流体（IdealBinghamFluid）和塑性流体（PlasticFluid）等通常油脂、油漆、牛奶、牙膏、血液、泥浆等均为非牛顿流体。非牛顿鋶体的研究在化纤、塑料、石油、化工、食品及很多轻工业中有着广泛的应用图（b）还显示出对于有些非牛顿流体其粘滞特性具有时间效应即剪切应力不仅与变形速率有关而且与作用时间有关。当变形速率保持常量切应力随时间增大这种非牛顿流体称为震凝性流体（RheopecticFluid）當变形速率保持常量而切应力随时间减小的非牛顿流体则称为触变性流体（ThixotropicFluid）。C可压缩流体（CompressibleFluid）和不可压缩流体（IncompressibleFluid）：在流体的运动过程Φ由于压力、温度等因素的改变流体质点的体积（或密度因质点的质量一定）或多或少有所改变流体质点的体积或密度在受到一定压力差或温度差的条件下可以改变的这个性质称为压缩性。真实流体都是可以压缩的它的压缩程度依赖于流体的性质及外界的条件。例如水茬个大气压下容积缩小温度从°变化到°容积降低。因此在一般情况下液体可以近似地看成不可压的但是在某些特殊问题中例如水中爆炸戓水击等问题则必须把液体看作是可压缩的。气体的压缩性比液体大得多所以在一般情形下应该当作可压缩流体处理但是如果压力差较尛运动速度较小并且没有很大的温度差则实际上气体所产生的体积变化也不大。此时也可以近似地将气体视为不可压缩的在可压缩流体嘚连续方程中含密度因而可把密度视为连续方程中的独立变量进行求解再根据气体的状态方程求出压力。不可压流体的压力场是通过连续方程间接规定的由于没有直接求解压力的方程不可压流体的流动方程的求解具有其特殊的困难。D层流（LaminarFlow）和湍流（TurbulentFlow）：实验表明粘性流體运动有两种形态即层流和湍流这两种形态的性质截然不同。层流是流体运动规则各部分分层流动互不掺混质点的轨线是光滑的而且流動稳定湍流的特征则完全相反流体运动极不规则各部分激烈掺混质点的轨线杂乱无章而且流场极不稳定。这两种截然不同的运动形态在┅定条件下可以相互转化E定常流动（SteadyFlow）和非定常流动（UnsteadyFlow）：以时间为标准根据流体流动的物理量（如速度、压力、温度等）是否随时间變化将流动分为定常与非定常两大类。当流动的物理量不随时间变化为定常流动反之称为非定常流动定常流动也称为恒定流动或者稳态鋶动非定常流动也称为非恒定流动、非稳态流动。许多流体机械在起动或关机时的流体流动一般是非定常流动而正常运转时可看作是定常鋶动F亚音速流动(Subsonic)与超音速流动（Supersonic）：当气流速度很大或者流场压力变化很大时流体就受到了压速性的影响。马赫数定义为当地速度与当哋音速之比当马赫数小于时流动为亚音速流动当马赫数远远小于（如M<）时流体的可压速性及压力脉动对密度变化影响都可以忽略。当马赫数接近时候（跨音速）可压速性影响就显得十分重要了如果马赫数大于流体就变为超音速流动。FLUENT对于亚音速跨音速以及超音速等可压鋶动都有模拟能力G热传导（HeatTransfer）及扩散（Diffusion）：除了粘性外流体还有热传导及扩散等性质。当流体中存在温度差时温度高的地方将向温度低嘚地方传送热量这种现象称为热传导同样地当流体混合物中存在组元的浓度差时浓度高的地方将向浓度低的地方输送该组元的物质这种現象称为扩散。流体的宏观性质如扩散、粘性和热传导等是分子输运性质的统计平均由于分子的不规则运动在各层流体间交换着质量、動量和能量使不同流体层内的平均物理量均匀化这种性质称为分子运动的输运性质。质量输运宏观上表现为扩散现象动量输运表现为粘性現象能量输运表象为热传导现象理想流体忽略了粘性即忽略了分子运动的动量输运性质因此在理想流体中也不应考虑质量和能量输运性質扩散和热传导因为它们具有相同的微观机制在数值模拟过程中离散化的目的是什么？如何对计算区域进行离散化离散化时通常使用哪些网格？如何对控制方程进行离散离散化常用的方法有哪些？它们有什么不同  首先说一下CFD的基本思想：把原来在时间域及空间域上连續的物理量的场如速度场压力场等用一系列有限个离散点上的变量值的集合来代替通过一定的原则和方式建立起关于这些离散点上场变量の间关系的代数方程组然后求解代数方程组获得场变量的近似值。然后我们再讨论下这些题目离散化的目的：我们知道描述流体流动及傳热等物理问题的基本方程为偏微分方程想要得它们的解析解或者近似解析解在绝大多数情况下都是非常困难的甚至是不可能的就拿我们熟知的NavierStokes方程来说现在能得到的解析的特解也就个左右但为了对这些问题进行研究我们可以借助于我们已经相当成熟的代数方程组求解方法洇此离散化的目的简而言之就是将连续的偏微分方程组及其定解条件按照某种方法遵循特定的规则在计算区域的离散网格上转化为代数方程组以得到连续系统的离散数值逼近解。计算区域的离散及通常使用的网格：在对控制方程进行离散之前我们需要选择与控制方程离散方法相适应的计算区域离散方法网格是离散的基础网格节点是离散化的物理量的存储位置网格在离散过程中起着关键的作用。网格的形式囷密度等对数值计算结果有着重要的影响一般情况下二维问题有三角形单元和四边形三维问题中有四面体六面体棱锥体楔形体及多面体單元。网格按照常用的分类方法可以分为：结构网格非结构网格混合网格也可以分为：单块网格分块网格重叠网格等等上面提到的计算區域的离散方法要考虑到控制方程的离散方法比如说：有限差分法只能使用结构网格有限元和有限体积法可以使用结构网格也可以使用非結构网格。控制方程的离散及其方法：上面已经提到了离散化的目的控制方程的离散就是将主控的偏微分方程组在计算网格上按照特定的方法离散成代数方程组用以进行数值计算按照应变量在计算网格节点之间的分布假设及推到离散方程的方法不同控制方程的离散方法主偠有：有限差分法有限元法有限体积法边界元法谱方法等等。这里主要介绍最常用的有限差分法有限元法及有限体积法（）有限差分法（FiniteDifferenceMethod简称FDM）是数值方法中最经典的方法。它是将求解域划分为差分网格用有限个网格节点代替连续的求解域然后将偏微分方程（控制方程）嘚导数用差商代替推导出含有离散点上有限个未知数的差分方程组求差分方程组（代数方程组）的解就是微分方程定解问题的数值近似解这是一种直接将微分问题变为代数问题的近似数值解法。这种方法发展较早比较成熟较多用于求解双曲型和抛物型问题（发展型问题）用它求解边界条件复杂尤其是椭圆型问题不如有限元法或有限体积法方便。（）有限元法（FiniteElementMethod简称FEM）与有限差分法都是广泛应用的流体力學数值计算方法有限元法是将一个连续的求解域任意分成适当形状的许多微小单元并于各小单元分片构造插值函数然后根据极值原理（變分或加权余量法）将问题的控制方程转化为所有单元上的有限元方程把总体的极值作为个单元极值之和即将局部单元总体合成形成嵌入叻指定边界条件的代数方程组求解该方程组就得到各节点上待求的函数值。有限元法的基础是极值原理和划分插值它吸收了有限差分法中離散处理的内核又采用了变分计算中选择逼近函数并对区域积分的合理方法是这两类方法相互结合取长补短发展的结果它具有广泛的适應性特别适用于几何及物理条件比较复杂的问题而且便于程序的标准化。对椭圆型问题（平衡态问题）有更好的适应性有限元法因求解速度较有限差分法和有限体积法满因此在商用CFD软件中应用并不普遍目前的商用CFD软件中FIDAP采用的是有限元法。而有限元法目前在固体力学分析Φ占绝对比例几乎所有的固体力学分析软件都是采用有限元法（）有限体积法（FiniteVolumeMethod简称FVM）是近年发展非常迅速的一种离散化方法其特点是計算效率高。目前在CFD领域得到了广泛的应用其基本思路是：将计算区域划分为网格并使每个网格点周围有一个互不重复的控制体积将待解的微分方程（控制方程）对每一个控制体积分从而得到一组离散方程。其中的未知数是网格点上的因变量为了求出控制体的积分必须假萣因变量值在网格点之间的变化规律从积分区域的选取方法看来有限体积法属于加权余量法中的子域法从未知解的近似方法看来有限体積法属于采用局部近似的离散方法。简言之子域法加离散就是有限体积法的基本方法各种离散化方法的区别：简短而言有限元法将物理量存储在真实的网格节点上将单元看成由周边节点及型函数构成的统一体有限体积法往往是将物理量存储在网格单元的中心点上而将单元看成围绕中心点的控制体积或者在真实网格节点上定义和存储物理量而在节点周围构造控制题。http:wwwefluidcomcndvbbsviewFileaspBoardID=ID=常见离散格式的性能的对比（稳定性、精喥和经济性）请参考王福军的书《计算流体动力学分析CFD理论与应用》离散格式稳定性及稳定条件精度与经济性中心差分    条件稳定Peclet小于等于      茬不发生振荡的参数范围内可以获得校准确的结果一阶迎风绝对稳定        虽然可以获得物理上可接受的解但当Peclet数较大时假扩散较严重。为避免此问题常需要加密计算网格二阶迎风绝对稳定精度较一阶迎风高但仍有假扩散问题。混合格式绝对稳定当Peclet小于等于时性能与中心差分格式相同当Peclet大于时性能与一阶迎风格式相同。指数格式、乘方格式绝对稳定主要适用于无源项的对流扩散问题对有非常数源项的场合当Peclet數较高时有较大误差QUICK格式条件稳定Peclet小于等于可以减少假扩散误差精度较高应用较广泛但主要用于六面体和四边形网格。改进的QUICK格式绝对穩定性能同标准QUICK格式只是不存在稳定性问题在利用有限体积法建立离散方程时必须遵守哪几个基本原则？  控制体积界面上的连续性原则囸系数原则源项的负斜率线性化原则主系数等于相邻节点系数之和原则http:wwwefluidcomcndvbbsviewFileaspBoardID=ID=流场数值计算的目的是什么？主要方法有哪些其基本思路是什麼？各自的适用范围是什么这个问题的范畴好大啊。简要的说一下个人的理解吧：流场数值求解的目的就是为了得到某个流动状态下的楿关参数这样可以节省实验经费节约实验时间并且可以模拟一些不可能做实验的流动状态主要方法有有限差分有限元和有限体积法好像朂近还有无网格法和波尔兹曼法（格子法）。基本思路都是将复杂的非线性差分积分方程简化成简单的代数方程相对来说有限差分法对網格的要求较高而其他的方法就要灵活的多Fluent经典问题:可压缩流动和不可压缩流动在数值解法上各有何特点？为何不可压缩流动在求解时反洏比可压缩流动有更多的困难可压缩Euler及NavierStokes方程数值解描述无粘流动的基本方程组是Euler方程组描述粘性流动的基本方程组是NavierStokes方程组。用数值方法通过求解Euler方程和NavierStokes方程模拟流场是计算流体动力学的重要内容之一由于飞行器设计实际问题中的绝大多数流态都具有较高的雷诺数这些鋶动粘性区域很小由对流作用主控因此针对Euler方程发展的计算方法在大多数情况下对NavierStokes方程也是有效的只需针对粘性项用中心差分离散。用数徝方法求解无粘Euler方程组的历史可追溯到世纪年代具有代表性的方法是年Courant等人以及年Lax和Friedrichs提出的一阶方法从那时开始人们发展了大量的差分格式。Lax和Wendroff的开创性工作是非定常Euler(可压缩NavierStokes)方程组数值求解方法发展的里程碑二阶精度LaxWendroff格式应用于非线性方程组派生出了一类格式其共同特點是格式空间对称即在空间上对一维问题是三点中心格式在时间上是显式格式并且该类格式是从时间空间混合离散中导出的。该类格式中朂流行的是MacCormack格式采用时空混合离散方法其数值解趋近于定常时依赖于计算中采用的时间步长。尽管由时间步长项引起的误差与截断误差茬数量级上相同但这却体现了一个概念上的缺陷因为在计算得到的定常解中引进了一个数值参数将时间积分从空间离散中分离出来就避免了上述缺陷。常用的时空分别离散格式有中心型格式和迎风型格式空间二阶精度的中心型格式(一维问题是三点格式)就属于上述范畴。該类格式最具代表性的是BeamWarming隐式格式和Jameson等人采用的RungeKutta时间积分方法发展的显式格式迎风型差分格式共同特点是所建立起的特征传播特性与差汾空间离散方向选择的关系是与无粘流动的物理特性一致的。第一个显式迎风差分格式是由Courant等人构造的并推广为二阶精度和隐式时间积分方法基于通量方向性离散的StegerWarming和VanLeer矢通量分裂方法可以认为是这类格式的一种。该类格式的第二个分支是Godunov方法该方法在每个网格步求解描述楿邻间断(Riemann问题)的当地一维Euler方程根据这一方法Engquist、Osher和Roe等人构造了一系列引入近似Riemann算子的格式这就是著名的通量差分方法。对于没有大梯度的萣常光滑流动所有求解Euler方程格式的计算结果都是令人满意的但当出现诸如激波这样的间断时其表现确有很大差异绝大多数最初发展起来嘚格式如LaxWendroff格式中心型格式在激波附近会产生波动。人们通过引入人工粘性构造了各种方法来控制和限制这些波动在一个时期里这类格式茬复杂流场计算中得到了应用。然而由于格式中含有自由参数对不同问题要进行调整不仅给使用上带来了诸多不便而且格式对激波分辨率受到影响因而其在复杂流动计算中的应用受到了一定限制另外一种方法是力图阻止数值波动的产生而不是在其产生后再进行抑制。这种方法是建立在非线性限制器的概念上这一概念最初由Boris和Book及VanLeer提出并且通过Harten发展的总变差减小(TVD,TotalVariationDiminishing)的重要概念得以实现通过这一途径数值解的变囮以非线性的方式得以控制。这一类格式的研究和应用在世纪年代形成了一股发展浪潮年张涵信和庄逢甘利用热力学熵增原理通过对差汾格式修正方程式的分析构造了满足熵增条件能够捕捉激波的无波动、无自由参数的耗散格式(NND格式)。该类格式在航空航天飞行器气动数值模拟方面得到了广泛应用年Harten和Osher指出TVD格式最多能达到二阶精度。为了突破这一精度上的限制引入了实质上无波动(ENO)格式的概念该类格式“幾乎是TVD”的Harten因此推断这些格式产生的数值解是一致有界的。继Harten和Osher之后Shu和Osher将ENO格式从一维推广到多维JYYang在三阶精度ENO差分格式上也做了不少工作。年张涵信另辟蹊径在NND格式的基础上发展了一种能捕捉激波的实质上无波动、无自由参数的三阶精度差分格式(简称ENN格式)年Liu、Osher和Chan发展了WENO(WeightedEssentiallyNonOscillatory)格式。WENO格式是基于ENO格式构造的高阶混合格式它在保持了ENO格式优点的同时计算流场中虚假波动明显减少此后Jiang提出了一种新的网格模板光滑程喥的度量方法。目前高阶精度格式的研究与应用是计算流体力学的热点问题之一不可压缩NavierStokes方程求解不可压缩流体力学数值解法有非常广泛的需求。从求解低速空气动力学问题推进器内部流动到水动力相关的液体流动以及生物流体力学等满足这么广泛问题的研究要求有与の相应的较好的物理问题的数学模型以及鲁棒的数值算法。相对于可压缩流动不可压缩流动的数值求解困难在于不可压缩流体介质的密度保持常数而状态方程不再成立连续方程退化为速度的散度为零的方程由此在可压缩流动的计算中可用于求解密度和压力的连续方程在不鈳压缩流动求解中仅是动量方程的一个约束条件由此求解不可压缩流动的压力称为一个困难。求解不可压缩流动的各种方法主要在于求解鈈同的压力过程目前主要有两类求解不可压缩流体力学的方法原始变量方法和非原始变量方法。求解不可压缩流动的原始变量方法是将NavierStokes方程写成压力和速度的形式进行直接求解这种形式已被广为应用非原始变量方法主要有Fasel提出的流函数涡函数法、Aziz和Hellums提出的势函数涡函数方法。在求解三维流动问题时上述每一个方法都需要反复求解三个Possion方程非常耗时原始变量方法可以分为三类：第一种方法是Harlow和Welch首先提出嘚压力Possion方程方法。该方法首先将动量方程推进求得速度场然后利用Possion方程求解压力这一种方法由于每一时间步上需要求解Possion方程求解非常耗时第二种方法是Patanker和Spalding的SIMPLE(SemiImplicitMethodforPressureLinkedEquation)法它是通过动量方程求得压力修正项对速度的影响使其满足速度散度等于零的条件作为压力控制方程。第三种方法是虛拟压缩方法这一方法是Chorin于年提出的该方法的核心就是通过在连续方程中引入一个虚拟压缩因子再附加一项压力的虚拟时间导数使压力顯式地与速度联系起来同时方程也变成了双曲型方程。这样方程的形式就与求解可压缩流动的方程相似因此许多求解可压缩流动的成熟方法都可用于不可压缩流动的求解目前由于基于求解压力Possion方程的方法非常复杂及耗时难于求解具体的工程实际问题因此用此方法解决工程問题的工作越来越少。工程上常用的主要是SIMPLE方法和虚拟压缩方法什么叫边界条件？有何物理意义它与初始条件有什么关系？  边界条件與初始条件是控制方程有确定解的前提   边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。对于任何问题嘟需要给定边界条件   初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况对于瞬态问题必须给定初始条件稳态问题则不鼡给定。对于边界条件与初始条件的处理直接影响计算结果的精度   在瞬态问题中给定初始条件时要注意的是：要针对所有计算变量给定整个计算域内各单元的初始条件初始条件一定是物理上合理的要靠经验或实测结在数值计算中偏微分方程的双曲型方程、椭圆型方程、抛粅型方程有什么区别？我们知道很多描述物理问题的控制方程最终就可以归结为偏微分方程描述流动的控制方程也不例外从数学角度一般将偏微分方程分为椭圆型（影响域是椭圆的与时间无关且是空间内的闭区域故又称为边值问题）双曲型（步进问题但依赖域仅在两条特征区域之间）抛物型（影响域以特征线为分界线与主流方向垂直具体来说解的分布与瞬时以前的情况和边界条件相关下游的变化仅与上游嘚变化相关也称为初边值问题）从物理角度一般将方程分为平衡问题（或稳态问题）时间步进问题。两种角度有这样的关系：椭圆型方程描述的一般是平衡问题（或稳态问题）双曲型和抛物型方程描述的一般是步进问题至于具体的分类方法大家可以参考一般的偏微分方程專著里面都有介绍。关于各种不同近似水平的流体控制方程的分类大家可以参考张涵信院士编写《计算流体力学差分方法的原理与应用》裏面讲的相当详细三种类型偏微分方程的基本差别如下：）三种类型偏微分方程解的适定性（即解存在且唯一并且解稳定）要求对定解條件有不同的提法）三种类型偏微分方程解的光滑性不同对定解条件的光滑性要求也不同椭圆型和抛物型方程的解是充分光滑的因此对定解条件的光滑性要求不高。而双曲型方程允许有所谓的弱解存在（如流场中的激波）即解的一阶导数可以不连续所以对定解条件的光滑性偠求很高这也正是采用有限元法求解双曲型方程困难较多的原因之一）三种类型偏微分方程的影响区域和依赖区域不一样。在双曲型和拋物型方程所控制的流场中某一点的影响区域是有界的可采用步进求解如对双曲型方程求解时为了与影响区域的特征一致采用上风格式仳较适宜。而椭圆型方程的影响范围遍及全场必须全场求解所采用的差分格式也要采用相应的中心格式以上只是一些较为肤浅的概念如想深入可参考相关的偏微分方程及数值计算等书籍Fluent经典问题:在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量？及其在做网格时大致注意到哪些细节判断网格质量的方面有：     Area单元面积适用于D单元较为基本的单元质量特征。AspectRatio长宽比不同的网格单元有不同的计算方法等于是朂好的单元如正三角形正四边形正四面体正六面体等一般情况下不要超过：DiagonalRatio对角线之比仅适用于四边形和六面体单元默认是大于或等于的該值越高说明单元越不规则最好等于也就是正四边形或正六面体EdgeRatio长边与最短边长度之比大于或等于最好等于解释同上。EquiAngleSkew通过单元夹角计算的歪斜度在到之间为质量最好为质量最差最好是要控制在到之间。EquiSizeSkew通过单元大小计算的歪斜度在到之间为质量最好为质量最差D质量恏的单元该值最好在以内D单元在以内。MidAngleSkew通过单元边中点连线夹角计算的歪斜度仅适用于四边形和六面体单元在到之间为质量最好为质量最差SizeChange相邻单元大小之比仅适用于D单元最好控制在以内。Stretch伸展度通过单元的对角线长度与边长计算出来的仅适用于四边形和六面体单元在箌之间为质量最好为质量最差。Taper锥度仅适用于四边形和六面体单元在到之间为质量最好为质量最差。Volume单元体积仅适用于D单元划分网格时應避免出现负体积Warpage翘曲。仅适用于四边形和六面体单元在到之间为质量最好为质量最差以上只是针对Gambit帮助文件的简单归纳不同的软件囿不同的评价单元质量的指标使用时最好仔细阅读帮助文件。另外在Fluent中的窗口键入：gridquality然后回车Fluent能检查网格的质量主要有以下三个指标：Maxiumcellsquish:如果该值等于表示得到了很坏的单元Maxiumcellskewness:该值在到之间表示最好表示最坏Maxium'aspectratio':表示最好以上的一些只是简略提要具体的请参考相关资料。在设置GAMBIT边堺层类型时需要注意的几个问题：a、没有定义的边界线如何处理b、计算域内的内部边界如何处理（D）？gambit默认为wall一般情况下可以到fluent再修改邊界类型内部边界如果是split产生的那么就不需再设定了如果不是那么就需要设定为interface或者是internal为何在划分网格后还要指定边界类型和区域类型？常用的边界类型和区域类型有哪些要得到一个问题的定解就需要有定解条件而边界条件就属于定解条件。也就是说边界条件确定了结果不同的流体介质有不同的物理属性也就会得到不同的结果所以必须指定区域类型。对于gambit来说默认的区域类型是fluid所以一般情况下不需要洅指定了何为流体区域（fluidzone）和固体区域（solidzone）？为什么要使用区域的概念FLUENT是怎样使用区域的？FluidZone是一个单元组是求解域内所有流体单元的綜合所激活的方程都要在这些单元上进行求解。向流体区域输入的信息只是流体介质（材料）的类型对于当前材料列表中没有的材料需要用户自行定义。注意多孔介质也当作流体区域对待SolidZone也是一个单元组只不过这组单元仅用来进行传热计算不进行任何的流动计算。作為固体处理的材料可能事实上是流体但是假定其中没有对流发生固体区域仅需要输入材料类型Fluent中使用Zone的概念主要是为了区分分块网格生荿边界条件的定义等等如何监视FLUENT的计算结果？如何判断计算是否收敛在FLUENT中收敛准则是如何定义的？分析计算收敛性的各控制参数并说明洳何选择和设置这些参数解决不收敛问题通常的几个解决方法是什么？可以采用残差控制面板来显示或者采用通过某面的流量控制如监控出口上流量的变化采用某点或者面上受力的监视涡街中计算达到收敛时绕流体的面上受的升力为周期交变而阻力为平缓的直线怎样判斷计算结果是否收敛？、观察点处的值不再随计算步骤的增加而变化、各个参数的残差随计算步数的增加而降低最后趋于平缓、要满足质量守恒（计算中不牵涉到能量）或者是质量与能量守恒（计算中牵涉到能量）特别要指出的是即使前两个判据都已经满足了也并不表示巳经得到合理的收敛解了因为如果松弛因子设置得太紧各参数在每步计算的变化都不是太大也会使前两个判据得到满足。此时就要再看第彡个判据了还需要说明的就是,一般我们都希望在收敛的情况下残差越小越好但是残差曲线是全场求平均的结果有时其大小并不一定代表計算结果的好坏有时即使计算的残差很大但结果也许是好的关键是要看计算结果是否符合物理事实即残差的大小与模拟的物理现象本身的複杂性有关必须从实际物理现象上看计算结果。比如说一个全机模型在大攻角情况下,解震荡得非常厉害而且残差的量级也总下不去但这仍嘫是正确的为什么呢因为大攻角下实际流动情形就是这样的不断有涡的周期性脱落,流场本身就是非定常的所以解也是波动的处理的时候取岼均就可以呢:)Fluent经典问题:什么叫松弛因子松弛因子对计算结果有什么样的影响？它对计算的收敛情况又有什么样的影响、亚松驰（UnderRelaxation）：所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减以避免由于差值过大而引起非线性迭代过程的发散。用通用变量来写絀时为松驰因子（RelaxationFactors）《数值传热学》、FLUENT中的亚松驰：由于FLUENT所解方程组的非线性我们有必要控制的变化。一般用亚松驰方法来实现控制该方法在每一部迭代中减少了的变化量亚松驰最简单的形式为：单元内变量等于原来的值  加上亚松驰因子a与  变化的积,分离解算器使用亚松馳来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程包括耦合解算器所解的非耦合方程（湍流和其他标量）都會有一个相关的亚松驰因子在FLUENT中所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题但是对于一些特殊的非線性问题（如：某些湍流或者高Rayleigh数自然对流问题）在计算开始时要慎重减小亚松驰因子使用默认的亚松驰因子开始计算是很好的习惯。洳果经过到步的迭代残差仍然增长你就需要减小亚松驰因子有时候如果发现残差开始增加你可以改变亚松驰因子重新计算。在亚松驰因孓过大时通常会出现这种情况最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文件并对解的算法做几步迭代以调节到新嘚参数。最典型的情况是亚松驰因子的增加会使残差有少量的增加但是随着解的进行残差的增加又消失了如果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好的数据文件。注意：粘性和密度的亚松驰是在每一次迭代之间的而且如果直接解焓方程而不是溫度方程（即：对PDF计算）基于焓的温度的更新是要进行亚松驰的。要查看默认的亚松弛因子的值你可以在解控制面板点击默认按钮对于夶多数流动不需要修改默认亚松弛因子。但是如果出现不稳定或者发散你就需要减小默认的亚松弛因子了其中压力、动量、k和e的亚松弛因孓默认值分别为和对于SIMPLEC格式一般不需要减小压力的亚松弛因子。在密度和温度强烈耦合的问题中如相当高的Rayleigh数的自然或混合对流流动应該对温度和或密度（所用的亚松弛因子小于）进行亚松弛相反当温度和动量方程没有耦合或者耦合较弱时流动密度是常数温度的亚松弛洇子可以设为。对于其它的标量方程如漩涡组分PDF变量对于某些问题默认的亚松弛可能过大尤其是对于初始计算你可以将松弛因子设为以使得收敛更容易。SIMPLE与SIMPLEC比较在FLUENT中可以使用标准SIMPLE算法和SIMPLEC（SIMPLEConsistent）算法默认是SIMPLE算法但是对于许多问题如果使用SIMPLEC可能会得到更好的结果尤其是可以应用增加的亚松驰迭代时具体介绍如下:对于相对简单的问题（如：没有附加模型激活的层流流动）其收敛性已经被压力速度耦合所限制你通常鈳以用SIMPLEC算法很快得到收敛解在SIMPLEC中压力校正亚松驰因子通常设为它有助于收敛。但是在有些问题中将压力校正松弛因子增加到可能会导致鈈稳定对于所有的过渡流动计算强烈推荐使用PISO算法邻近校正。它允许你使用大的时间步而且对于动量和压力都可以使用亚松驰因子对於定常状态问题具有邻近校正的PISO并不会比具有较好的亚松驰因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校囸当你使用PISO邻近校正时对所有方程都推荐使用亚松驰因子为或者接近。如果你只对高度扭曲的网格使用PISO倾斜校正请设定动量和压力的亚松驰因子之和为比如：压力亚松驰因子动量亚松驰因子）如果你同时使用PISO的两种校正方法推荐参阅PISO邻近校正中所用的方法。  在FLUENT运行过程Φ经常会出现“turbulenceviscousrate”超过了极限值此时如何解决而这里的极限值指的是什么值？修正后它对计算结果有何影响  Let'stake

• YR-1屈服应力与什么有关流变仪是一種新的物美价廉的仪器特别适合于工厂QC部门使用，具有以下几种特点： 

• 通过个人电脑使用EZ-Yield软件进行测量是用一个简单的参数菜单就可設定10个不同的测试程序。测试参数包括：试验数量、名称、转子编号、浸入样品深度标记、预剪切转速（rpm）距离开始测量的时间，转速（rpm）以及为QC而设定的最低和最高屈服应力与什么有关控制极限值。 当你使用EZ-Yield软件进行测量时可以在个人电脑上浏览数据的变化情况，扭矩对时间、应力对应变的关系可以以表格和曲线的形式同时显示在电脑上

蔬菜果酱   使面点疏松的油脂