HyperMesh和自定义导出模板简化CFD分析

2016-10-18 by:CAE仿真在线来源:互联网

项目介绍

亚利桑那州立大学(ASU)的物质运输和能源工程学院(SEMTE),希望通过使用商业化的通用前处理软件HyperMesh,来完成复杂几何外形的网格划分,以便在谱元法计算流体程序Nek5000中使用。由于CFD计算程序需要使用3D六面体单元,在使用Nek5000程序时就会遇到划六面体网格的挑战。

SEMTE的科研人员创建了一项研究项目,来开发HyperMesh网格和Nek5000代码间转化的工具。这个转化工具可以将HyperMesh网格转化为Nek5000程序可识别的格式。在使用这个导出模板后,整个流程具有了更好的用户体验并且更加不容易出错。首先在HyperMesh中生成网格,然后在导出模板中整理好所有的数据并将其分类,以便网格能够以Nek5000的正确格式导入到用户模块中。

挑战

综合模拟计算流体实验研究团队重点开发并利用了许多工具应用于并行计算机,研究分析复杂的工程问题和物理系统。他们用开源的计算流体求解器Nek5000来进行他们的研究。Nek5000是基于谱元法的一个求解器,由阿贡国家实验室提供技术支持。Nek5000被广泛的应用于许多问题中,比如反应堆堆芯热工水力学、血管流动、大气和海洋建模以及燃烧等领域。SEMTE的科研人员使用这个程序来研究一些领域的基本的物理问题,比如热传递问题和湍流问题,它们目前的商业应用主要是在航空航天和风能行业。该团队(综合模拟计算流体实验研究团队)正在开发各种应用来提高Nek5000的功能。

Nek5000采用了谱元法,将有限元求解器的几何灵活性与全局谱元法的高精度特性相结合。这个程序的优点在于它能够以很高的精度来研究十分复杂的流动形式,包括几何上的复杂性和多组分/多物理流动。

对于高保真方法预测的湍流流动,比如直接数值模拟法,在这种方法下,N-S方程在足够精细的网格上被离散,所有尺度的运动都得到直接的求解,不需要模型。这是一种优秀的方法,它具有最小的离散度和耗散误差。使用Nek5000最大的困难在于它要求采用3D六面体单元。通常来说,大多数CFD软件使用四面体网格,因为对于一般的几何外形,四面体网格更加容易生成。这给亚利桑纳州立大学的研究人员们带来了挑战——如何利用多功能、先进的网格划分软件(如HyperMesh)。在使用Nek5000的时候,生成高质量的六面体网格。

“我很高兴能使用HyperMesh,我认为它是我用过的最好的网格生成工具。从很多方面来讲,它都是一个非常有价值的软件。比如它出色的参考资料,这些资料都十分容易获得,对于其他的商用软件而言就不一定了。使用HyperMesh后,我们在许多问题上的工作效率都变快了,划分网格的流程也十分简单。事实上,在过去的一年里,我们团队里面的硕士生有一半使用过HyperMesh自定义导出模块,为研究流经几何结构的流动生成网格,否则这些网格将不能被识别。”

Philip Sakievich

博士生、助理研究员

亚利桑那州立大学航空航天与机械工程系

解决方案

为了应对这项挑战,综合模拟计算流体实验室的研究人员立了一个课题来研究开发转换工具,使从商用软件HyperMesh中导出的网格转换为Nek5000程序能够使用的格式。在这个课题开始前,研究人员还没有一个划分网格的流程。大多数的网格划分都是由15-20年前开发的定制软件来处理的,而且这些软件从那以后只有较少的更新。另外一些用户针对一些特定的问题,在一些软件比如Matlab中,创建他们自己的网格划分工具。当处理不规则几何外形的时候就变得更加复杂了。在建立一个标准的流程的时候,一个成熟的工具比如HyperMesh是很有必要的,因为它允许用户通过几何划分网格。 SEMTE的研究人员选择HyperMesh,是因为它中立于求解器,出众的参考文档,以及开放的架构。

“使用HyperMesh后,我们团队能生成高质量的网格。此外,由于它能读取所有的常用网格格式,其它团队用另外一些软件划分的网格,我们可以用它进行转换。这极大的增加了我们与其他研究人员的合作能力,并且这使得我们能集中更多精力在我们关心的实际研究问题上。这都是我们团队获得的好处。”Philip Sakievich说道。

有了工具后,工程师们开始研究一些很小的问题,这些问题由5到10个单元组成。在项目的过程中,他们一直在开发和调试程序,随后在更大的范围去测试。经过几周的深入研究和编写程序,这些研究人员最终确定了导出模板,编写了一个总体上对用户更友好、更不容易出错的流程。网格在HyperMesh中生成以后,导出模板将所有的数据组织分类,以便它能在Nek5000中以需要的格式在用户模板中导入。转换工具先将每个单元的节点坐标从原始的HyperMesh格式转换成Nek5000的数据格式。用户也通过HyperMesh施加了边界条件,这些数据将在几何转换完成后,转换成Nek5000的格式。然后所有的这些数据都被写入到一个文件中,这个文件就可以提交到Nek5000中进行求解。工程师们也为网格加入了边界点,以便增加几何的可变形特性,并且转换器会自动识别这个基于HyperMesh的单元类型的信息。