CN107562984A - 极限流场的保向共轭映射方法和装置 - Google Patents
极限流场的保向共轭映射方法和装置 Download PDFInfo
- Publication number
- CN107562984A CN107562984A CN201710584904.7A CN201710584904A CN107562984A CN 107562984 A CN107562984 A CN 107562984A CN 201710584904 A CN201710584904 A CN 201710584904A CN 107562984 A CN107562984 A CN 107562984A
- Authority
- CN
- China
- Prior art keywords
- curved surface
- point
- field
- limit flow
- pneumatic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Processing Or Creating Images (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种极限流场的保向共轭映射方法和装置。该方法包括:根据气动曲面的吸力面数据,获取所述气动曲面的极限流场;其中,所述气动曲面为三维曲面;将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场;将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场;其中,所述参数平面为将所述气动曲面进行参数化后得到的二维平面。本发明的方法,使得基于极限流场的气动曲面设计与加工从三维空间转化为平面,使得设计与加工更加精确和便捷,提高了气动曲面设计的精确度。
Description
技术领域
本发明涉及流体力学技术,尤其涉及一种极限流场的保向共轭映射方法和装置。
背景技术
极限流场是气动力学领域的概念,工程上选取靠近物面特定距离处的流场分布作为极限流场,以得到极限流线。极限流场是最为贴近物面(例如叶片表面)的流场分布,能从极大程度上反映叶片表面的物理环境以及曲面的几何特征。通过极限流线的形态可以判断流动的分离、再附和旋涡等现象,这与叶片气动性能密切相关。研究黏性气体的分离、流动情况以及极限流场的性态,对飞机、导弹和高超声速机动飞行器等航空航天高精度器件的设计,具有重要意义。
与气动性能直接相关的器件,外观多为较为复杂的三维曲面,其主要表现为:曲率分布不均且变化率较大。对于较为复杂的自由曲面而言,表面的极限流场也会相对紊乱。除此之外,由于三维空间中对曲面数据散点与离散矢量场的拟合过程中存在误差,而对于大部分气动曲面而言,设计或加工中精度的控制至关重要。例如,航空航天领域的处于复杂的空气流场中压气机叶片和机翼等部件,这些流场对压气机叶片和机翼的设计起着至关重要的作用。
但是,目前工程上难以利用极限流场的性质对三维曲面形态的叶片或者机翼进行精确的设计与加工。
发明内容
本发明提供一种极限流场的保向共轭映射方法和装置,用以解决现有技术中难以利用极限流场的性质对三维曲面形态的叶片或者机翼进行精确的设计与加工的技术问题。
第一方面,本发明提供一种极限流场的保向共轭映射方法,包括:
根据气动曲面的吸力面数据,获取所述气动曲面的极限流场;其中,所述气动曲面为三维曲面;
将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场;
将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场;其中,所述参数平面为将所述气动曲面进行参数化后得到的二维平面。
可选的,所述气动曲面为S(u1,u2),所述将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场,具体包括:
对||pk-qk||进行优化,得到并通过得到所述气动曲面上与所述极限流场的数据点pk的欧式距离最小的位置点qk;所述k为所述极限流场中的数据点的个数;
根据所述数据点pk和所述位置点qk,得到
根据公式得到dui/dt;其中,为所述气动曲面S(u1,u2)在qk点处关于参数u1的偏导数,为所述气动曲面S(u1,u2)在qk点处关于参数u2的偏导数;所述dpk/dt为所述数据点pk的极限流场方向向量;
根据公式dqk/dt=∑i=1,2Si·dui/dt,得到所述气动曲面上qk点处的流线场方向向量,并根据所述气动曲面上qk点处的流线场向量得到所述气动曲面的流线场。
可选的,所述将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场,具体包括:
根据公式计算得到所述参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b);其中,
根据所述参数平面上的点(u1,u2)处的方向向量(ξ,η),得到所述参数平面上的向量场。
进一步地,所述根据所述参数平面上的点(u1,u2)处的方向向量(a,b),得到所述参数平面上的向量场,具体包括:
对(ξ,η)=(a,b)进行归一化,得到归一化后的
根据所述参数平面上的点(u1,u2)处归一化后的(ξ,η),得到所述参数平面上的向量场。
可选的,所述极限流场为n*6的矩阵,或者6*n的矩阵,所述n为所述极限流场中的数据点的个数;
当所述极限流场为n*6的矩阵时,所述n*6的矩阵的任一行的前三列为所述行对应的数据点的位置坐标信息,所述任一行的后三列为所述数据点在在x、y、z轴上的矢量;
当所述极限流场为6*n的矩阵时,所述6*n的矩阵的任一列的前三行为所述列对应的数据点的位置坐标信息,所述任一列的后三行为所述数据点在在x、y、z轴上的矢量。
第二方面,本发明提供一种极限流场的保向共轭映射装置,包括:
获取模块,用于根据气动曲面的吸力面数据,获取所述气动曲面的极限流场;其中,所述气动曲面为三维曲面;
第一映射模块,用于将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场;
第二映射模块,用于将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场;其中,所述参数平面为将所述气动曲面进行参数化后得到的二维平面。
可选的,所述气动曲面为S(u1,u2),所述第一映射模块,具体用于:
对||pk-qk||进行优化,得到并通过得到所述气动曲面上与所述极限流场的数据点pk的欧式距离最小的位置点qk;所述k为所述极限流场中的数据点的个数;
根据所述数据点pk和所述位置点qk,得到
根据公式得到dui/dt;其中,为所述气动曲面S(u1,u2)在qk点处关于参数u1的偏导数,为所述气动曲面S(u1,u2)在qk点处关于参数u2的偏导数;所述dpk/dt为所述数据点pk的极限流场方向向量;
根据公式dqk/dt=∑i=1,2Si·dui/dt,得到所述气动曲面上qk点处的流线场方向向量,并根据所述气动曲面上qk点处的流线场向量得到所述气动曲面的流线场。
可选的,所述第二映射模块,具体用于:
根据公式计算得到所述参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b);其中,
根据所述参数平面上的点(u1,u2)处的方向向量(ξ,η),得到所述参数平面上的向量场。
进一步地,所述第二映射模块,具体用于:
对(ξ,η)=(a,b)进行归一化,得到归一化后的
根据所述参数平面上的点(u1,u2)处归一化后的(ξ,η),得到所述参数平面上的向量场。
可选的,所述极限流场为n*6的矩阵,或者6*n的矩阵,所述n为所述极限流场中的数据点的个数;
当所述极限流场为n*6的矩阵时,所述n*6的矩阵的任一行的前三列为所述行对应的数据点的位置坐标信息,所述任一行的后三列为所述数据点在在x、y、z轴上的矢量;
当所述极限流场为6*n的矩阵时,所述6*n的矩阵的任一列的前三行为所述列对应的数据点的位置坐标信息,所述任一列的后三行为所述数据点在在x、y、z轴上的矢量。
本发明提供的极限流场的保向共轭映射方法和装置,首先根据气动曲面的吸力面数据,获取气动曲面的极限流场,然后将该极限流场映射至气动曲面的曲面切平面上,得到气动曲面的流线场,进而将气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到参数平面上的向量场,从而使得三维气动曲面上的流线场与参数平面上向量场存在一一对应关系,进而使得基于极限流场的气动曲面设计与加工从三维空间转化为平面,使得设计与加工更加精确和便捷,提高了气动曲面设计的精确度。
附图说明
图1为本发明提供的极限流场的处理系统的结构示意图;
图2为本发明提供的极限流场的保向共轭映射方法实施例一的流程示意图;
图3a为本发明提供的极限流场的示意图;
图3b为本发明提供的极限流场到气动曲面的映射示意图;
图3c为本发明提供的流线场到参数平面的映射示意图;
图4为本发明提供的极限流场的保向共轭映射方法实施例二的流程示意图;
图5为本发明提供的极限流场的保向共轭映射方法实施例三的流程示意图;
图6为本发明提供的极限流场的保向共轭映射装置实施例的结构示意图。
具体实施方式
本发明提供的极限流场的保向共轭映射方法,可以适用于图1所示的极限流场的处理系统。如图1所示,该处理系统包括:输入装置10、处理设备11和输出设备12。其中,输入装置可以向用户提供的数据的输入接口,使得用户可以通过输入装置向处理设备输入机翼三维曲面的数据或者压气机叶片的三维曲面数据等吸力面数据;处理设备可以为计算机,还可以为服务器,还可以是其他的具有处理和计算能力的设备;输出设备可以为能够将处理设备输出的结果通知给用户的设备,可选的,该输出设备可以为一显示设备,从而可以将处理设备计算的结果可视化的显示给用户。
由于与气动性能直接相关的器件,外观多为较为复杂的三维曲面,其表面的极限流场相对紊乱。除此之外,由于三维空间中对曲面数据散点与离散矢量场的拟合过程中存在误差,而对于大部分气动曲面而言,设计或加工中精度的控制至关重要。例如,航空航天领域的处于复杂的空气流场中压气机叶片和机翼等部件,这些流场对压气机叶片和机翼的设计起着至关重要的作用。但是,目前工程上难以利用极限流场的性质对三维曲面形态的叶片或者机翼进行精确的设计与加工。
因此,本发明提供的极限流场的保向共轭映射方法,通过将三维曲面上的极限流场通过保向共轭映射,得到曲面参数域的向量场(二维),从而使得三维曲面上的流线场与参数平面上向量场存在一一对应关系,进而使得基于极限流场的气动曲面设计与加工从三维空间转化为平面,直接在参数平面上进行设计与加工更加便捷,同时保证了曲面上任意点的插值性。
下面以具体地实施例对本发明的技术方案进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。
图2为本发明提供的极限流场的保向共轭映射方法实施例一的流程示意图。本实施例的执行主体可以是任一的处理设备,该处理设备可以是计算机等具体处理计算功能的设备,下述方法实施例以执行主体是计算机为例来进行说明。本实施例涉及的是计算机通过将三维曲面上的极限流场通过保向共轭映射,得到二维的参数平面上的向量场,使得三维曲面上的流线场与参数平面上向量场存在一一对应关系,进而使得基于极限流场的气动曲面设计与加工从三维空间转化为平面,使得设计与加工更加精确和便捷的具体过程。
如图1所示,该方法包括如下步骤:
S101:根据气动曲面的吸力面数据,获取所述气动曲面的极限流场;其中,所述气动曲面为三维曲面。
具体的,在进行压气机叶片或者机翼的设计时,首先需要获取叶片或者机翼等三维气动曲面的吸力面数据,该吸力面数据可以为能够反映叶片或者机翼等物体表面的物理环境以及曲面的几何特征的数据。当计算机获得气动曲面的吸力面数据之后,可以对该吸力面数据进行预处理。以压气机叶片的吸力面数据为例,预处理过程可以包括:将得到的吸力面数据统一尺度,然后通过适当的坐标变换,将吸力面数据置于合适的坐标系下,并且去掉冗余或噪声数据,最后以样条函数形式表示该叶片的气动曲面。可选的,该气动曲面的样条函数形式可以为S(u1,u2)。
当计算机对所获取的气动曲面的吸力面数据进行预处理后,将处理后的吸力面数据导入计算流体动力学(Computational Fluid Dynamics,简称CFD)软件中,获得气动曲面表面的极限流场。可选的,极限流场可以参见图3a所示。需要说明的是,本实施例涉及的气动曲面的极限流场,均是以靠近气动曲面特定距离处的流场分布作为极限流场。可选的,上述CFD软件可以是ANSYS软件。
可选的,上述所得到的气动曲面的极限流场可以为n*6的矩阵,还可以为6*n的矩阵,n为极限流场中的数据点的个数。
当极限流场为n*6的矩阵时,该矩阵的一行为极限流场中的一个数据点对应的信息。对于该n*6的矩阵的某一行而言,该行的前三列为该行对应的数据点的位置坐标信息,后三列为该数据点在在x、y、z轴上的矢量。
当极限流场为6*n的矩阵时,该矩阵的一列为极限流场中的一个数据点对应的信息。对于该6*n的矩阵的某一列而言,该列的前三行为该列对应的数据点的位置坐标信息,后三行为该数据点在在x、y、z轴上的矢量。
S102:将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场。
具体的,当获取气动曲面的极限流场之后,可以将该极限流场映射至气动曲面的曲面切片面上,得到气动曲面的流线场。需要说明的是,上述极限流场中包括n个数据点pk,每个数据点pk,均会对应该数据点pk处的极限流场方向向量,可以记作dpk/dt,k为极限流场中的数据点的个数,k=1,2,…,n。。另外,对于极限流场中的每个数据点,在气动曲面上均会有一个点qk与pk所对应,该qk为气动曲面上与极限流场的数据点pk的欧式距离最小的位置点,每一个qk均会对应一个曲面切平面,每个qk在其所对应的曲面切平面上均会具有该点的流线场方向向量,所有qk点处的流线场方向向量共同构成气动曲面的流线场。也就是说,在获取气动曲面的流线场时,只要能够获得每个qk点处的流线场方向向量,就可以得到三维气动曲面的流线场。
可选的,极限流场到气动曲面的映射示意图,可以参见图3b所示。
S103:将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场;其中,所述参数平面为将所述气动曲面进行参数化后得到的二维平面。
具体的,上述气动曲面进行参数化之后可以得到一个二维的参数平面,具体的参数化过程可以参见现有技术,在此不再赘述。当计算机得到三维气动曲面的流线场之后,就可以通过保向共轭映射的方式,将该流线场映射至该参数平面,从而得到参数平面上的向量场。需要说明的是,保向共轭的映射方式,可以使得两个连续流之间不仅保持轨道的拓扑结构,而且保持参数之间具有一定的关联关系或者等价关系。本实施例中,气动曲面上的某一个点qk点均会对应至参数平面上的某一个点,也就是说,气动曲面上的每个点的流线场方向向量,与参数平面上的每个点的方向向量之前均具有一定的关联关系。当计算机得到参数平面上的每个点的方向向量之后,就可以得到参数平面上的向量场,从而使得三维气动曲面上的流线场与参数平面上向量场存在一一对应关系,进而使得基于极限流场的气动曲面设计与加工从三维空间转化为平面,使得设计与加工更加精确和便捷,提高了气动曲面设计的精确度。可选的,气动曲面的流线场到参数平面的映射,可以参加图3c所示。
本发明提供的极限流场的保向共轭映射方法,首先根据气动曲面的吸力面数据,获取气动曲面的极限流场,然后将该极限流场映射至气动曲面的曲面切平面上,得到气动曲面的流线场,进而将气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到参数平面上的向量场,从而使得三维气动曲面上的流线场与参数平面上向量场存在一一对应关系,进而使得基于极限流场的气动曲面设计与加工从三维空间转化为平面,使得设计与加工更加精确和便捷,提高了气动曲面设计的精确度。
图4为本发明提供的极限流场的保向共轭映射方法实施例二的流程示意图。本实施例涉及的是获取气动曲面流线场的具体过程,在上述实施例的基础上,进一步地,上述S102具体可以包括:
S201:对||pk-qk||进行优化,得到并通过得到所述气动曲面上与所述极限流场的数据点pk的欧式距离最小的位置点qk;所述k为所述极限流场中的数据点的个数。
S202:根据所述数据点pk和所述位置点qk,得到
具体的,当计算机获取气动曲面的极限流场之后,根据极限流场中的数据点pk以及对||pk-qk||进行优化,得到进而通过得到上述气动曲面上与极限流场的数据点pk的欧式距离最小的位置点qk,该步骤实际上是求向量的最小模长,从而得到qk。可选的,可以通过matlab自带的优化工具箱求解得到。
S203:根据公式得到dui/dt;其中,为所述气动曲面S(u1,u2)在qk点处关于参数u1的偏导数,Su2为所述气动曲面S(u1,u2)在qk点处关于参数u2的偏导数。
S204:根据公式dqk/dt=∑i=1,2Si·dui/dt,得到所述气动曲面上qk点处的流线场方向向量,并根据所述气动曲面上qk点处的流线场向量得到所述气动曲面的流线场。
具体的,由于极限流场的表现形式为(pk,dpk/dt),pk为极限流场中的第k个数据点的位置坐标向量(三维),dpk/dt为该数据点对应的极限流场方向向量,因此,计算机在通过S201的步骤得到数据点pk对应的qk之后,将pk的极限流场方向向量dpk/dt平移到三维气动曲面上的点qk,即计算pk的极限流场方向向量dpk/dt对应至qk点处曲面切平面上的流线场方向向量,即向qk点处的曲面切平面做正交投影,通过公式得到dui/dt(i=1,2)或者duj/dt(j=1,2),然后根据公式dqk/dt=∑i=1,2Si·dui/dt,得到所述气动曲面上qk点处的流线场方向向量dqk/dt。需要说明的是,上述i=1,2实际上是i=1和2时两个方程式的简写,其中包含两个未知数,分别是du1/dt、du2/dt,i取1和2时可以得到两个方程,联立求解从而得到du1/dt、du2/dt的一组解。
当得到气动曲面上的每个qk点处的流线场方向向量dqk/dt,这些qk点处的流线场方向向量就可以构成上述气动曲面的流线场。其中,上述气动曲面的极限流场的数据点信息可以记作(pk,dpk/dt),对应三维曲面上流线场数据可以记作(qk,dqk/dt),三维曲面表达式可以记作S(u1,u2)。
本实施例提供的极限流场的保向共轭映射方法,根据极限流场中的数据点找到三维气动曲面上距离该数据点的欧式距离最短的位置点qk,然后根据该数据点处的极限流场方向向量dpk/dt以及公式得到dui/dt,进而根据公式dqk/dt=∑i=1,2Si·dui/dt得到气动曲面上qk点处的流线场方向向量dqk/dt,并根据气动曲面上qk点处的流线场向量得到气动曲面的流线场。本实施例可以得到气动曲面的极限流场对应的流线场,为计算机将三维的极限流场转换为二维参数平面上的向量场提供了良好的铺垫和基础,提高了计算机将三维的极限流场转换为二维参数平面上的向量场的转换效率。
图5为本发明提供的极限流场的保向共轭映射方法实施例三的流程示意图。本实施例涉及的是将气动曲面的流线场通过保向共轭映射的方式映射至二维的参数平面,得到参数平面上的向量场的具体过程。在上述实施例的基础上,进一步地,上述S103可以包括如下步骤:
S301:根据公式计算得到所述参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b);其中,
具体的,本实施例中(u1,u2)表示二维参数平面上的一个点,即三维气动曲面上的三维坐标的点(例如qk点)对应到二维参数平面上的点(例如(u1,u2)表示的点)时,可以通过二维坐标表示。本实施例中的t(u1,u2)实际上为qk点处的流线场方向向量,t=t(u1,u2)。
当计算机得到气动曲面的流线场之后,即计算机得到了气动曲面上的每一个qk点处的流线场方向向量dqk/dt,然后计算机根据公式从而可以得到参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b),即可以通过以下方程组得到参数平面上的点(u1,u2)处的方向向量(ξ,η),具体的方程组如下:
可选的,可以将气动曲面的第一基本形式以惯用记法E、F、G表示,其中E=Su1·Su1,F=Su1·Su2,G=Su2·Su2。通过matlab求解,获得下列解表达式的数值形式:故而,就得到了参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b)。
S302:根据所述参数平面上的点(u1,u2)处的方向向量(ξ,η),得到所述参数平面上的向量场。
基于上述过程,计算机可以得到参数平面上每个点在参数平面上的方向向量,因此,结合每个点的处的方向向量(ξ,η),就可以得到参数平面上的向量场,进而实现将三维的极限流场转换为二维的参数平面上的向量场,从而使得三维气动曲面上的流线场与参数平面上向量场存在一一对应关系,进而使得基于极限流场的气动曲面设计与加工从三维空间转化为平面,使得设计与加工更加精确和便捷,提高了气动曲面设计的精确度。
可选的,为了提高三维空间中对曲面数据散点与离散矢量场的拟合效果,本实施例还可以对上述参数平面上的点(u1,u2)处的方向向量(ξ,η)进行归一化,从而得到归一化后的基于该归一化的过程,计算机可以得到参数平面上每个点在参数平面上的归一化后的方向向量,因此,结合每个点处的归一化后的(ξ,η),就可以得到参数平面上的向量场,提高了三维空间中对曲面数据散点与离散矢量场的拟合度。
本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
图6为本发明提供的极限流场的保向共轭映射装置实施例的结构示意图。该装置可以集成在上述处理设备中,还可以为独立的处理设备。该装置可以通过软件、硬件或者软硬件结合的方式实现。如图所示,该装置包括获取模块21、第一映射模块22和第二映射模块23。
具体的,获取模块21,用于根据气动曲面的吸力面数据,获取所述气动曲面的极限流场;其中,所述气动曲面为三维曲面;
第一映射模块22,用于将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场;
第二映射模块23,用于将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场;其中,所述参数平面为将所述气动曲面进行参数化后得到的二维平面。
可选的,所述气动曲面为S(u1,u2),所述第一映射模块22,具体用于:
对||pk-qk||进行优化,得到并通过得到所述气动曲面上与所述极限流场的数据点pk的欧式距离最小的位置点qk;所述k为所述极限流场中的数据点的个数;
根据所述数据点pk和所述位置点qk,得到
根据公式i=1,2,得到dui/dt;其中,为所述气动曲面S(u1,u2)在qk点处关于参数u1的偏导数,为所述气动曲面S(u1,u2)在qk点处关于参数u2的偏导数;所述dpk/dt为所述数据点pk的极限流场方向向量;
根据公式dqk/dt=∑i=1,2Si·dui/dt,得到所述气动曲面上qk点处的流线场方向向量,并根据所述气动曲面上qk点处的流线场向量得到所述气动曲面的流线场。
可选的,所述第二映射模块23,具体用于:
根据公式计算得到所述参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b);其中,
根据所述参数平面上的点(u1,u2)处的方向向量(ξ,η),得到所述参数平面上的向量场。
进一步地,所述第二映射模块23,具体用于:
对(ξ,η)=(a,b)进行归一化,得到归一化后的
根据所述参数平面上的点(u1,u2)处归一化后的(ξ,η),得到所述参数平面上的向量场。
可选的,所述极限流场为n*6的矩阵,或者6*n的矩阵,所述n为所述极限流场中的数据点的个数;
当所述极限流场为n*6的矩阵时,所述n*6的矩阵的任一行的前三列为所述行对应的数据点的位置坐标信息,所述任一行的后三列为所述数据点在在x、y、z轴上的矢量;
当所述极限流场为6*n的矩阵时,所述6*n的矩阵的任一列的前三行为所述列对应的数据点的位置坐标信息,所述任一列的后三行为所述数据点在在x、y、z轴上的矢量。
本发明实施例提供的极限流场的保向共轭映射装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种极限流场的保向共轭映射方法,其特征在于,包括:
根据气动曲面的吸力面数据,获取所述气动曲面的极限流场;其中,所述气动曲面为三维曲面;
将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场;
将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场;其中,所述参数平面为将所述气动曲面进行参数化后得到的二维平面。
2.根据权利要求1所述的方法,其特征在于,所述气动曲面为S(u1,u2),所述将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场,具体包括:
对||pk-qk||进行优化,得到并通过得到所述气动曲面上与所述极限流场的数据点pk的欧式距离最小的位置点qk;所述k为所述极限流场中的数据点的个数;
根据所述数据点pk和所述位置点qk,得到
根据公式得到dui/dt;其中,为所述气动曲面S(u1,u2)在qk点处关于参数u1的偏导数,为所述气动曲面S(u1,u2)在qk点处关于参数u2的偏导数;所述dpk/dt为所述数据点pk的极限流场方向向量;
根据公式dqk/dt=∑i=1,2Si·dui/dt,得到所述气动曲面上qk点处的流线场方向向量,并根据所述气动曲面上qk点处的流线场向量得到所述气动曲面的流线场。
3.根据权利要求2所述的方法,其特征在于,所述将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场,具体包括:
根据公式计算得到所述参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b);其中,
根据所述参数平面上的点(u1,u2)处的方向向量(ξ,η),得到所述参数平面上的向量场。
4.根据权利要求3所述的方法,其特征在于,所述根据所述参数平面上的点(u1,u2)处的方向向量(a,b),得到所述参数平面上的向量场,具体包括:
对(ξ,η)=(a,b)进行归一化,得到归一化后的
根据所述参数平面上的点(u1,u2)处归一化后的(ξ,η),得到所述参数平面上的向量场。
5.根据权利要求1-4任一项所述的方法,其特征在于,所述极限流场为n*6的矩阵,或者6*n的矩阵,所述n为所述极限流场中的数据点的个数;
当所述极限流场为n*6的矩阵时,所述n*6的矩阵的任一行的前三列为所述行对应的数据点的位置坐标信息,所述任一行的后三列为所述数据点在在x、y、z轴上的矢量;
当所述极限流场为6*n的矩阵时,所述6*n的矩阵的任一列的前三行为所述列对应的数据点的位置坐标信息,所述任一列的后三行为所述数据点在在x、y、z轴上的矢量。
6.一种极限流场的保向共轭映射装置,其特征在于,包括:
获取模块,用于根据气动曲面的吸力面数据,获取所述气动曲面的极限流场;其中,所述气动曲面为三维曲面;
第一映射模块,用于将所述极限流场映射至所述气动曲面的曲面切平面上,得到所述气动曲面的流线场;
第二映射模块,用于将所述气动曲面的流线场通过保向共轭映射的方式,映射至参数平面,得到所述参数平面上的向量场;其中,所述参数平面为将所述气动曲面进行参数化后得到的二维平面。
7.根据权利要求6所述的装置,其特征在于,所述气动曲面为S(u1,u2),所述第一映射模块,具体用于:
对||pk-qk||进行优化,得到并通过得到所述气动曲面上与所述极限流场的数据点pk的欧式距离最小的位置点qk;所述k为所述极限流场中的数据点的个数;
根据所述数据点pk和所述位置点qk,得到
根据公式得到dui/dt;其中,为所述气动曲面S(u1,u2)在qk点处关于参数u1的偏导数,为所述气动曲面S(u1,u2)在qk点处关于参数u2的偏导数;所述dpk/dt为所述数据点pk的极限流场方向向量;
根据公式dqk/dt=∑i=1,2Si·dui/dt,得到所述气动曲面上qk点处的流线场方向向量,并根据所述气动曲面上qk点处的流线场向量得到所述气动曲面的流线场。
8.根据权利要求7所述的装置,其特征在于,所述第二映射模块,具体用于:
根据公式计算得到所述参数平面上的点(u1,u2)处的方向向量(ξ,η)=(a,b);其中,
根据所述参数平面上的点(u1,u2)处的方向向量(ξ,η),得到所述参数平面上的向量场。
9.根据权利要求8所述的装置,其特征在于,所述第二映射模块,具体用于:
对(ξ,η)=(a,b)进行归一化,得到归一化后的
根据所述参数平面上的点(u1,u2)处归一化后的(ξ,η),得到所述参数平面上的向量场。
10.根据权利要求6-9任一项所述的装置,其特征在于,所述极限流场为n*6的矩阵,或者6*n的矩阵,所述n为所述极限流场中的数据点的个数;
当所述极限流场为n*6的矩阵时,所述n*6的矩阵的任一行的前三列为所述行对应的数据点的位置坐标信息,所述任一行的后三列为所述数据点在在x、y、z轴上的矢量;
当所述极限流场为6*n的矩阵时,所述6*n的矩阵的任一列的前三行为所述列对应的数据点的位置坐标信息,所述任一列的后三行为所述数据点在在x、y、z轴上的矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710584904.7A CN107562984B (zh) | 2017-07-18 | 2017-07-18 | 极限流场的保向共轭映射方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710584904.7A CN107562984B (zh) | 2017-07-18 | 2017-07-18 | 极限流场的保向共轭映射方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107562984A true CN107562984A (zh) | 2018-01-09 |
CN107562984B CN107562984B (zh) | 2019-12-31 |
Family
ID=60973634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710584904.7A Active CN107562984B (zh) | 2017-07-18 | 2017-07-18 | 极限流场的保向共轭映射方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107562984B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003091739A (ja) * | 2001-09-17 | 2003-03-28 | Namco Ltd | 画像生成システム、プログラム及び情報記憶媒体 |
CN105298925A (zh) * | 2015-11-04 | 2016-02-03 | 北京航空航天大学 | 槽道展向高度渐缩变化的叶根开槽压气机静子 |
CN105631158A (zh) * | 2016-01-14 | 2016-06-01 | 北京航空航天大学 | 一种基于流线场理论的压气机叶片吸力面建模方法 |
-
2017
- 2017-07-18 CN CN201710584904.7A patent/CN107562984B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003091739A (ja) * | 2001-09-17 | 2003-03-28 | Namco Ltd | 画像生成システム、プログラム及び情報記憶媒体 |
CN105298925A (zh) * | 2015-11-04 | 2016-02-03 | 北京航空航天大学 | 槽道展向高度渐缩变化的叶根开槽压气机静子 |
CN105631158A (zh) * | 2016-01-14 | 2016-06-01 | 北京航空航天大学 | 一种基于流线场理论的压气机叶片吸力面建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107562984B (zh) | 2019-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102203782B (zh) | 求解拉格朗日形式的欧拉方程的数值方法 | |
Ruiz-Gironés et al. | Generation of curved high-order meshes with optimal quality and geometric accuracy | |
Chen | Aerothermodynamics of turbomachinery: analysis and design | |
CN103077259A (zh) | 高超声速导弹多场耦合动力学一体化仿真分析方法 | |
Ramírez et al. | New high-resolution-preserving sliding mesh techniques for higher-order finite volume schemes | |
CN103729564B (zh) | 一种基于粒子图像测速技术的压力场计算方法和装置 | |
Müller et al. | Low Reynolds number drag and particle collision efficiency of a cylindrical fiber within a parallel array | |
CN109408836A (zh) | 利用Boltzmann方程进行流体仿真的方法 | |
CN109325303A (zh) | 一种船舶回转运动模拟方法 | |
CN109740182A (zh) | 一种基于再生核粒子的无网格物理变形仿真方法 | |
Chakrabarti et al. | Representing rectangular jet dynamics through azimuthal fourier modes | |
CN114996866A (zh) | 一种用于射流式脉冲发动机内部的数学模型和离散方法 | |
Watanabe et al. | Moving computational domain method and its application to flow around a high-speed car passing through a hairpin curve | |
Hazra et al. | Simultaneous pseudo-timestepping for aerodynamic shape optimization problems with state constraints | |
Mölder | Reflection of curved shock waves | |
Takashima et al. | Wedge-cone waverider configuration for engine-airframe interaction | |
CN107562984A (zh) | 极限流场的保向共轭映射方法和装置 | |
Balan et al. | Hp-adaptivity on anisotropic meshes for hybridized discontinuous Galerkin scheme | |
Jin | Numerical simulation of wind turbine wakes based on actuator line method in NEK5000 | |
CN116467957A (zh) | 一种基于piv技术的超声速复杂流场的压力场重构方法 | |
Driant et al. | Numerical simulation and wind tunnel measurements on a tricycle wheel sub-system | |
Hao et al. | Numerical studies of acoustic diffraction by rigid bodies | |
Yang et al. | An immersed boundary method based on parallel adaptive Cartesian grids for high Reynolds number turbulent flow | |
Ma et al. | Aerodynamic configuration of the HCW based on the lifting body | |
Shirayama et al. | Flow visualization in computational fluid dynamics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |