专利名称：固体推进剂配方设计用羽流场数据结构网格化处理方法
技术领域：
本发明属于固体推进剂羽流场数据处理技术领域，尤其是涉及一种固体推进剂配方设计用羽流场数据结构网格化处理方法。
背景技术：
随着机械化、信息化技术的发展，现代武器系统越来越多的采用雷达、光电等自动跟踪、控制和制导手段，武器系统的自动化程度、打击精度等方面得到显著提高，极大的提升了弹箭武器系统的作战效能。但与此同时，弹箭武器发射和运载过程中的烟焰问题危害性日益显现，对武器系统信息化、智能化、隐身化性能影响很大，对武器系统快速、精确打击、隐身等功能造成了严重的不利影响，甚至已成为制约我军弹箭武器发挥正常效能、新装备订货以及新型武器装备发展的主要技术瓶颈。尤其是近年来，重型反坦克导弹研制过程中出现推进剂装药烟雾遮挡制导信号的问题，解决问题耗费巨大。固体推进剂是弹箭武器系统的发射能源，其低特征信号性能是保证武器系统战场使命完成的重要因素。羽流是从火箭发动机喷管喷射出来的羽毛状的高速高温燃气流，火箭羽流是一种气体分子浓度大、电子密度和电子碰撞频率都很高的弱等离子体，其与雷达微波之间会发生相互作用，使微波信号功率大为降低，影响导弹的制导。而羽流场是火箭发动机喷管喷射出来的羽毛状的高速高温燃气流的流体运动所占据的空间。为评价固体推进剂羽流特性，世界上一些发达国家都建立了相关测试评估方法，建设了各种实验测试设施来对固体推进剂特征信号性能进行检测及表征。国内在对国外测试设施研究的基础上，已建立起一套具有独立知识产权的低成本、易控制的固体推进剂羽流特性检测系统。但在低特征信号推进剂研制及评价过程中，解决问题需要耗费大量的人力物力，研制周期大幅度延长，因而带来巨大的经济损失和政治影响，现如今迫切需要虚拟试验对固体推进剂相关性能进行快速评价。而虚拟试验是借助于计算机的高速解算，按实际试验要求对基于推进剂燃烧和燃烧产物流动的数学模型进行模拟试验，不仅可以作为真实试验的前期准备工作或在一定程度上替代传统的试验(如一些极限工况)；并且能大幅减少真实试验次数，降低试验费用，缩短试验周期； 同时，具有较好的交互性，使得各种试验信息及时反馈；另外，虚拟试验不受气象条件、场地、时间和次数的限制，试验过程可方便实现回放、再现和重复。目前，对推进剂燃烧和燃烧产物流动的数学模型进行模拟试验时，通常都采用FLUNT软件对固体推进剂燃烧后的羽流场进行模拟，且采用FLUNT软件完成所设计固体推进剂的羽流场计算后，便能自动输出*.0ut文件供微波衰减计算、红外福射计算以及光电特性计算使用。但由于FULENT软件存盘文件dat中保存的是基于非结构网格的流场数据，而且包含了喷管区域和喷流区域。由于微波衰减计算、红外辐射计算以及光电特性计算等后续计算中只针对喷流区域，因而需要将FULENT软件计算后所获得喷流区域的流场数据提取出来进行结构网格转换，并转换成基于结构网格的对应数据以便后续计算使用。此外，由于FULENT只能保存粒子轨迹历史数据，而红外辐射计算和光电特性计算需要的是不同直径粒子的数密度和温度，因而还需通过对FULENT软件计算后所获得喷流区域的流场数据进行处理，并相应得出喷流区域中每一条粒子轨迹上各结构网格内不同直径粒子的数密度和温度。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种固体推进剂配方设计用羽流场数据结构网格化处理方法，其方法步骤简单、设计合理且实现方便、使用效果好，能简便将基于非结构网格的流场数据转换成基于结构网格的对应数据以便后续计算使用。为解决上述技术问题，本发明采用的技术方案是:一种固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于该方法包括以下步骤:步骤一、羽流场计算及数据存储:采用数据处理器调用FULENT软件且对所设计固体推进剂进行羽流场计算后，将羽流场计算结果自动输出并存储至与所述数据处理器相接的数据存储单元内；其中，所输出的羽流场计算结果中包括气相产物羽流场计算结果和凝聚相产物羽流场计算结果；步骤二、燃烧产物信息输入:通过参数输入单元输入所设计固体推进剂燃烧后所产生燃烧产物的数量m以及各燃烧产物的属性信息和平衡组分，各燃烧产物的属性信息均包括化学式和相态，其中相态为气相或凝聚相个燃烧产物的平衡组分之和=1，且m个燃烧产物中所有凝聚相产物的平衡组分之和记为nx ;同时，通过所述参数输入单元，对所设计固体推进剂燃烧后所产生燃烧产物中Al2O3颗粒的平均粒径Dn进行设定；所述参数输入单元与所述数据处理器相接；
步骤三、气相产物羽流场数据结构网格化处理，其处理过程如下:步骤301、喷流区域气相产物流场数据读取:采用所述数据处理器从步骤一中所输出的气相产物羽流场计算结果中读取所设计固体推进剂的喷流区域中所有非结构网格节点的流场数据；所述喷流区域为所述发动机喷管出口后方的矩形区域；步骤302、轴向坐标轴上非结构网格点提取:采用所述数据处理器从步骤301中所述喷流区域中的所有非结构网格节点中提取位于轴向坐标轴上的所有非结构网格点，本步骤中所提取的位于轴向坐标轴上的非结构网格点总数量为Nx;其中，轴向坐标轴为所述发动机喷管的中心轴线所在的横坐标轴，位于轴向坐标轴上的非结构网格点的径向坐标yh=0且其轴向坐标xh彡0，其中h为正整数，且h=l，2，-,Nx；步骤303、径向坐标轴上非结构网格点提取:采用所述数据处理器从步骤301中所述喷流区域中的所有非结构网格节点中提取位于径向坐标轴上的所有非结构网格点，所提取的位于径向坐标轴上的非结构网格点总数量为NYa ;其中，径向坐标轴为所述发动机喷管出口处所在的纵坐标轴且发动机喷管出口处的轴向坐标值为0，位于径向坐标轴上的非结构网格点的轴向坐标Xkl=O且其径向坐标ykl彡0，其中kl为正整数，且kl=l，2，…，Ny出；步骤304、构建结构网格图:将Nx条直线X=Xh和Ny出条直线y=ykl正交后，构建出一个包含(Nx — I) X (ΝΥ -1)个矩形网格的结构网格图；步骤305、气相产物羽流场数据结构网格化处理:采用所述数据处理器对步骤304中所构建结构网格图中各矩形网格的四个顶点的气相产物流场数据分别进行重新赋值；所有矩形网格中各顶点的气相产物流场数据重新赋值方法均相同，其中对于所构建结构网格图中任一个矩形网格的任一顶点的气相产物流场数据进行重新赋值时，所述数据处理器先从步骤301中所述喷流区域中的所有非结构网格节点中找出与当前被赋值顶点距离最近的非结构网格节点，并将所找出的非结构网格节点的气相产物流场数据赋给当前被赋值的顶点；步骤四、凝聚相产物羽流场数据结构网格化处理:所述数据处理器利用步骤304中所构建的结构网格图对M个不同粒径粒子的粒子轨迹数据分别进行结构网格化处理，过程如下:步骤401、结构网格化处理初始参数设定:采用所述参数输入单元，对M的取值和M个不同粒径粒子的粒径Dm分别进行设定；其中，r为正整数，且r=l，2，…，M ;步骤402、粒子轨迹数据读取:采用所述数据处理器从步骤一中所输出的凝聚相产物羽流场计算结果中读取所设计固体推进剂的所有粒子轨迹数据；其中，所读取的凝聚相产物羽流场计算结果包括粒子质量文件、粒子温度随轨迹变化的文件、粒子直径随轨迹变化的文件和粒子轨迹的时间步长文件；步骤403、发动机喷管入口矩形网格数量获取及各喷管入口矩形网格的上下边界确定:首先，采用所述数据处理器从步骤一中所输出的气相产物羽流场计算结果中读取所设计固体推进剂的发动机喷管区域中所有非结构网格节点的流场数据；之后，采用所述数据处理器从所读取的发动机喷管区域中的所有非结构网格节点中提取位于直线X=-Ad上的所有非结构网格节点，所提取的位于直线X=- Λ d上的非结构网格点总数量为Ny入；其中，位于直线χ=_Δ d上的非结构网格节点的轴向坐标Xk2=-A d且其径向坐标yk2彡O,其中k2为正整数，且k2=l，2，…，ΝΥλ; Ad为所述发动机喷管入口至喷管出口之间的距离；所获取的喷管入口矩形网格数量为(ΝΥλ_1)个，各喷管入口矩形网格的上下边界分别为上下相邻的两条直线η ；步骤404、凝聚相产物羽流场数据结构网格化处理:采用数据处理器并以步骤403中所述的(ΝΥλ-1)个喷管入口矩形网格分别作为起始网格进行粒子轨迹结构网格化处理，每个喷管入口矩形网格作为起始网格进行粒子轨迹结构网格化处理的过程均相同；并且，以任一个喷管入口矩形网格作为起始网格进行粒子轨迹结构网格化处理时，均对以当前所处理喷管入口矩形网格作为起始网格的M个不同粒径粒子轨迹分别进行结构网格化处理，且M个不同粒径粒子轨迹的结构网格化处理方法均相同；其中，对M个不同粒径粒子轨迹中的任一条粒子轨迹进行结构网格化处理时，均分别计算出当前所处理粒子轨迹在所述结构网格图中途经的所有矩形网格内的粒子轨迹网格化数据，且每个矩形网格内的粒子轨迹网格化数据均包括粒子质量、粒子数密度、粒子平均直径和粒子平均温度。上述固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征是:步骤404中对当前所处理粒子轨迹在所述结构网格图中途经的所有矩形网格内的粒子轨迹网格化数据进行计算时，按照所途经矩形网格的布设位置前后顺序由先至后进行计算。上述固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征是:步骤401中Μ=8，且8个不同粒径粒子的粒径Dm分别为Dnl、Dn2、Dn3> Dn4, Dn5, Dn6, Dn7和Dn8，其中，Dnl < Dn2 < Dn3 < Dn4 < Dn < Dn5 < Dn6 < Dn7 < Dn8，其中，Dn 为步骤二中所输入的 Al2O3 颗粒的平均粒径。上述固体推进剂配方设 计用羽流场数据结构网格化处理方法，其特征是:步骤一中对所设计固体推进剂进行羽流场计算时，所采用凝聚相产物的控制方程为拉格朗日颗粒轨道模型。上述固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征是:步骤404中对M个不同粒径粒子轨迹中的任一条粒子轨迹进行结构网格化处理时，其处理过程如下:步骤4041、发动机喷管入口起始网格确定及粒子质量流 η计算:所述数据处
理器从步骤402中所读取的所有粒子轨迹数据中找出当前所处理粒子轨迹的粒子轨迹数据，并根据所找出的粒子轨迹数据中的发动机喷管入口处轨迹点的径向坐标和步骤403中所确定的各喷管入口矩形网格的上下边界，对当前所处理粒子轨迹的喷管入口起始网格进行确定；所找出粒子轨迹数据中包括粒子质量文件、粒子温度随轨迹变化的文件、粒子直径随轨迹变化的文件和粒子轨迹的时间步长文件内所存储的数据；待当前所处理粒子轨迹的喷管入口起始网格确定后，所述数据处理器找出所确定喷管入口起始网格的顶点一和顶点二的气相产物流场数据，并自所找出的顶点一的气相产物流场数据中找出顶点一处的气相密度Pgl和气相轴向速度Ugl，且自所找出的顶点二的气相产物流场数据中找出顶点二处的气相密度P g2和气相轴向速度u g2，其中顶点一为当前所处理粒子轨迹的喷管入口起始网格的左上方顶点，且顶点二为当前所处理粒子轨迹的
喷管入口起始网格的左下方顶点；之后，所述数据处理器根据公式
权利要求
1.一种固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于该方法包括以下步骤: 步骤一、羽流场计算及数据存储:采用数据处理器调用FULENT软件且对所设计固体推进剂进行羽流场计算后，将羽流场计算结果自动输出并存储至与所述数据处理器相接的数据存储单元内；其中，所输出的羽流场计算结果中包括气相产物羽流场计算结果和凝聚相产物羽流场计算结果； 步骤二、燃烧产物信息输入:通过参数输入单元输入所设计固体推进剂燃烧后所产生燃烧产物的数量m以及各燃烧产物的属性信息和平衡组分，各燃烧产物的属性信息均包括化学式和相态，其中相态为气相或凝聚相，m个燃烧产物的平衡组分之和=1，m个燃烧产物中所有凝聚相产物的平衡组分之和记为nx ;同时，通过所述参数输入单元，对所设计固体推进剂燃烧后所产生燃烧产物中Al2O3颗粒的平均粒SDn进行设定；所述参数输入单元与所述数据处理器相接； 步骤三、气相产物羽流场数据结构网格化处理，其处理过程如下: 步骤301、喷流 区域气相产物流场数据读取:采用所述数据处理器从步骤一中的气相产物羽流场计算结果中读取所设计固体推进剂的喷流区域中所有非结构网格节点的流场数据；所述喷流区域为所述发动机喷管出口后方的矩形区域； 步骤302、轴向坐标轴上非结构网格点提取:采用所述数据处理器从步骤301中所述喷流区域中的所有非结构网格节点中提取位于轴向坐标轴上的所有非结构网格点，本步骤中所提取的位于轴向坐标轴上的非结构网格点总数量为Nx ;其中，轴向坐标轴为所述发动机喷管的中心轴线所在的横坐标轴，位于轴向坐标轴上的非结构网格点的径向坐标yh=0且其轴向坐标xh彡0，其中h为正整数，且h=l，2，-,Nx； 步骤303、径向坐标轴上非结构网格点提取:采用所述数据处理器从步骤301中所述喷流区域中的所有非结构网格节点中提取位于径向坐标轴上的所有非结构网格点，所提取的位于径向坐标轴上的非结构网格点总数量为NYa ;其中，径向坐标轴为所述发动机喷管出口处所在的纵坐标轴且发动机喷管出口处的轴向坐标值为0，位于径向坐标轴上的非结构网格点的轴向坐标Xkl=O且其径向坐标ykl彡0，其中kl为正整数，且kl=l，2，…，Ny出； 步骤304、构建结构网格图:将Nx条直线X=Xh和Ny a条直线y=ykl正交后，构建出一个包含(Nx—l)X (队出一 I)个矩形网格的结构网格图； 步骤305、气相产物羽流场数据结构网格化处理:采用所述数据处理器对步骤304中所构建结构网格图中各矩形网格的四个顶点的气相产物流场数据分别进行重新赋值；所有矩形网格中各顶点的气相产物流场数据重新赋值方法均相同，其中对于所构建结构网格图中任一个矩形网格的任一顶点的气相产物流场数据进行重新赋值时，所述数据处理器先从步骤301中所述喷流区域中的所有非结构网格节点中找出与当前被赋值顶点距离最近的非结构网格节点，并将所找出的非结构网格节点的气相产物流场数据赋给当前被赋值的顶占.步骤四、凝聚相产物羽流场数据结构网格化处理:所述数据处理器利用步骤304中所构建的结构网格图对M个不同粒径粒子的粒子轨迹数据分别进行结构网格化处理，过程如下: 步骤401、结构网格化处理初始参数设定:采用所述参数输入单元对M的取值和M个不同粒径粒子的粒径Dm分别进行设定；其中，r为正整数，且r=l，2，…，M ; 步骤402、粒子轨迹数据读取采用所述数据处理器从步骤一中所输出的凝聚相产物羽流场计算结果中读取所设计固体推进剂的所有粒子轨迹数据；其中，所读取的凝聚相产物羽流场计算结果包括粒子质量文件、粒子温度随轨迹变化的文件、粒子直径随轨迹变化的文件和粒子轨迹的时间步长文件； 步骤403、发动机喷管入口矩形网格数量获取及各喷管入口矩形网格的上下边界确定首先，采用所述数据处理器从步骤一中所输出的气相产物羽流场计算结果中读取所设计固体推进剂的发动机喷管区域中所有非结构网格节点的流场数据；之后，采用所述数据处理器从所读取的发动机喷管区域中的所有非结构网格节点中提取位于直线X=-Ad上的所有非结构网格节点，所提取的位于直线X=- Λ d上的非结构网格点总数量为Ny入；其中，位于直线X=- Λ d上的非结构网格节点的轴向坐标Xk2=- Δ d且其径向坐标yk2彡O，其中k2为正整数，且k2=l，2，…，ΝΥλ; Ad为所述发动机喷管入口至喷管出口之间的距离；所获取的喷管入口矩形网格数量为(ΝΥλ-1)个，各喷管入口矩形网格的上下边界分别为上下相邻的两条直线η ； 步骤404、凝聚相产物羽流场数据结构网格化处理采用数据处理器并以步骤403中所述的(ΝΥλ-1)个喷管入口矩形网格分别作为起始网格进行粒子轨迹结构网格化处理，每个喷管入口矩形网格作为起始网格进行粒子轨迹结构网格化处理的过程均相同；并且，以任一个喷管入口矩形网格作为起始网格进行粒子轨迹结构网格化处理时，均对以当前所处理喷管入口矩形网格作为起始网格的M个不同粒径粒子轨迹分别进行结构网格化处理，且M个不同粒径粒子轨迹的结构网格化处理方法均相同；其中，对M个不同粒径粒子轨迹中的任一条粒子轨迹进行结构网格化处理时，均分别计算出当前所处理粒子轨迹在所述结构网格图中途经的所有矩形网格内的粒子轨迹网格化数据，且每个矩形网格内的粒子轨迹网格化数据均包括粒子质量、粒子数密度、粒子平均直径和粒子平均温度。
2.按照权利要求I所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于步骤404中对当前所处理粒子轨迹在所述结构网格图中途经的所有矩形网格内的粒子轨迹网格化数据进行计算时，按照所途经矩形网格的布设位置前后顺序由先至后进行计算。
3.按照权利要求I或2所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于步骤401中Μ=8，且8个不同粒径粒子的粒径Dm分别为Dnl、D 2、Dn3、Dn4、Dn5、Dn6>Dn7 和 Dn8，其中，Dnl < Dn2 < Dn3 < Dn4 < Dn < Dn5 < Dn6 < Dn7 < Dn8，其中，Dn 为步骤二中所输入的Al2O3颗粒的平均粒径。
4.按照权利要求I或2所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于步骤一中对所设计固体推进剂进行羽流场计算时，所采用凝聚相产物的控制方程为拉格朗日颗粒轨道模型。
5.按照权利要求4所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于步骤404中对M个不同粒径粒子轨迹中的任一条粒子轨迹进行结构网格化处理时，其处理过程如下 步骤4041、发动机喷管入口起始网格确定及粒子质量流率计算所述数据处理器从步骤402中所读取的所有粒子轨迹数据中找出当前所处理粒子轨迹的粒子轨迹数据，并根据所找出的粒子轨迹数据中的发动机喷管入口处轨迹点的径向坐标和步骤403中所确定的各喷管入口矩形网格的上下边界，对当前所处理粒子轨迹的喷管入口起始网格进行确定；所找出粒子轨迹数据中包括粒子质量文件、粒子温度随轨迹变化的文件、粒子直径随轨迹变化的文件和粒子轨迹的时间步长文件内所存储的数据； 待当前所处理粒子轨迹的喷管入口起始网格确定后，所述数据处理器找出所确定喷管入口起始网格的顶点一和顶点二的气相产物流场数据，并自所找出的顶点一的气相产物流场数据中找出顶点一处的气相密度P gl和气相轴向速度Ugl，且自所找出的顶点二的气相产物流场数据中找出顶点二处的气相密度P g2和气相轴向速度Ug2，其中顶点一为当前所处理粒子轨迹的喷管入口起始网格的左上方顶点，且顶点二为当前所处理粒子轨迹的喷管入口起始网格的左下方顶点；之后，所述数据处理器根据公式
6.按照权利要求5所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于:步骤III中对当前所计算矩形网格内的粒子平均质量％和粒子平均温度Tp进行计算之前，先根据当前所计算矩形网格的上下边界径向坐标和左右边界轴向坐标，并结合步骤4041中所找出粒子轨迹数据中各轨迹点的轴向坐标和径向坐标，对当前所处理粒子轨迹在当前所计算矩形网格内的轨迹点总数量K和K个轨迹点的轴向坐标与径向坐标分别进行确定。
7.按照权利要求4所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于:步骤一中进行羽流场计算时，所述数据处理器先根据预先设定的喷管几何参数与喷流计算域范围，建立对所述发动机喷管的内外羽流场进行数值计算的二维轴对称模型；之后，所述数据处理器调用CFD前处理器，生成所设计固体推进剂的羽流场计算域网格图，且所述CFD前处理器为GAMBIT软件；然后，所述数据处理器调用FULENT软件，并结合预先设计的喷管几何参数、喷流计算域范围和燃烧室工作参数，对所设计固体推进剂进行羽流场计算，并输出羽流场计算结果； 其中，对燃烧室工作参数进行设定时，其设定过程如下:首先，通过所述参数输入单元对发动机的燃烧室内压强P。、燃烧室绝热温度、环境压力和环境温度Tif分别进行设定；之后，再对燃烧产物平衡组成进行设定；然后，通过所述参数输入单元对所设计固体推进剂燃烧后所产生Al2O3颗粒的平均粒径Dn进行设定； 对燃烧产物平衡组成进行设定时，所设定的燃烧产物平衡组成中包括所设计固体推进剂燃烧后所产生的(m-Q)个气相产物的平衡组分和Al2O3颗粒的平衡组分；所设定的(m-Q)个气相产物的平衡组分分别为步骤二中所输入的所设计固体推进剂燃烧后所产生(m-Q)个气相产物的平衡组分；当Q=I时，所设定的Al2O3颗粒的平衡组分为步骤二中所输入的所设计固体推进剂燃烧后所产生Al2O3颗粒的平衡组分；当Q > I时，所设定的Al2O3颗粒的平衡组分为步骤二中所述的nx。
8.按照权利要求1或2所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于:步骤302中提取出位于轴向坐标轴上的所有非结构网格点后，所述数据处理器将所提取的位于轴向坐标轴上的所有非结构网格点按照轴向坐标由小至大的顺序进行排列，并由左至右标注于所述轴向坐标轴上；步骤303中提取出位于径向坐标轴上的所有非结构网格点后，所述数据处理器将所提取的位于径向坐标轴上的所有非结构网格点按照径向坐标由小至大的顺序进行排列，并由下至上标注于所述径向坐标轴上。
9.按照权利要求1或2所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于:步骤一中所述气相产物羽流场计算结果为存储在*_tec.dat文件内的计算结果；步骤402中所述粒子质量文件、粒子温度随轨迹变化的文件、粒子直径随轨迹变化的文件和粒子轨迹的时间步长文件，分别为FULENT软件计算完成后自动保存的*_mass.fvp> *_temp.fvp> *_diam.fvp 和 *_time.fvp 文件。
10.按照权利要求1或2所述的固体推进剂配方设计用羽流场数据结构网格化处理方法，其特征在于:步骤301中对所设计固体推进剂的喷流区域中所有非结构网格节点的流场数据进行读取时，首先建立F个动态一维数组，其中F=F1+F2，其中Fl=7，F2为所设计固体推进剂燃烧后所产生气相产物的数量^个动态一维数组分别用于存储N个非结构网格点的F个流场变量数据，并且F 个流场变量数据分别为轴向坐标、径向坐标、温度、压强、F2个气相产物的平衡组分、凝聚相浓度、气相轴向速度和气相密度。
全文摘要
本发明公开了一种固体推进剂配方设计用羽流场数据结构网格化处理方法，包括步骤一、羽流场计算及数据存储；二、燃烧产物信息输入；三、气相产物羽流场数据结构网格化处理，过程如下喷流区域气相产物流场数据读取、构建结构网格图和气相产物羽流场数据结构网格化处理；四、凝聚相产物羽流场数据结构网格化处理，过程如下初始参数设定、粒子轨迹数据读取、发动机喷管入口处矩形网格数量获取及各喷管入口矩形网格的上下边界确定和凝聚相产物羽流场数据结构网格化处理。本发明方法步骤简单、设计合理且实现方便、使用效果好，能简便、快速将基于非结构网格的流场数据转换成基于结构网格的对应数据以便后续计算使用。
文档编号G06F17/50GK103218504SQ20131017879
公开日2013年7月24日 申请日期2013年5月15日 优先权日2013年5月15日
发明者赵凤起, 肖川, 向红军, 李猛, 徐司雨, 李丽, 梁勇, 孙美 申请人:西安近代化学研究所