原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
流体界面不稳定性数值模拟中界面的一种处理方法
尤迎玖·，李晓林一，王继海，，蔚喜军+
·北京应用物理与计算数学研究所，计算物理实验室
北京8009信箱，100088
ofMathematicalSciencesIndianaUniv·PurdueUnivat
”Department Indianapolis
Indianapolis，IN46202，USA
摘要
在界面不稳定性的数值模拟中，Level
set方法局限于密度不同的单种介质或绝热指数
相同的两种介质的情形。针对于此，本文在Levelset方法的基础上引入了Front
Tracking
方法的思想，即在界面上强制加入接触间断边界条件，从而避免了界面上参数的非物理振
算，得到了扰动随时间发展的图象。
关键词：界面不稳定性，LevelSet方法，TVD格式
1．引言
在界面不稳定性的数值模拟中，所用的数值格式对随时间发展的界面上的接触间断应
铲
’ 有较高的分辨能力。传统的二阶TVD格式较好地解决了激波传播问胚．但由于其格式耗散
口
^^ 较大，对接触间断的求解往往难以得到较高的精度。因此用这种格式计算界面不稳定性问
矿 题，几个时间步长后初始极薄的扰动界面就会弥散到几个网格上。加了人工压缩项的TVD
格式，即TVD／AC格式[【1能较好地解决这个问题，计算中数值耗散并不随时间步长的增加
冉- 而增大。
在过去几年中，对运动界面的研究使FrontTracking方法[2】[3]和LevelSet[4][5]方法得
到了充分的发展。Front
set方法是一种Euler方法，
Tracking方法在数值计算中非常繁琐而且计算量极大。而LeveI
都可以重新构造界面。但流体质点上的Level函数并不是实际的物理量，它仅仅需要满足
一个额外的方程即可，因此Level函数的加入不对流体力学方程组产生任何影响。这种方
法思路清晰且容易实现，因而得到了广泛的应用。
Set方法计
遗憾的是，当两种理想流体的绝热指数不同时，用TVD／AC格式附以Level
算所得界面上物理量出现了非物理振荡[1】。对于这个问题，目前的看法和认识有分歧，有
关的研究也较少，但这个问题的解决对界面不稳定性的研究有着十分重要的意义。本文认
Set方法并不
为，由于Level函数的引入并不影响实际流场中的物理量的求解．因而Level
是界面上物理量产生非物理振荡的原因，其真正的原因在于根据差分格式计算界面一侧的
点时用到界面另一侧满足不同：次态方程的点上的物理量。
175
鉴于此，本文尝试在处理界面时，引入Front
Tracking的思想，在界面上强制给出接触
文档评论（0）