-
公开(公告)号:CN113158339B
公开(公告)日:2022-10-18
申请号:CN202110408433.0
申请日:2021-04-16
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/15 , G06F30/28 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明公开了一种针对SST湍流模型的湍流长度尺度修正方法,本修正方法以无量纲速度散度λl的值为基本自变量来确定修正源项的大小,通过控制函数tanh(h2(η‑h3))‑1实现了对修正源项作用区域的控制。本发明方法不依赖于壁面距离这一参数,而是根据流场中速度散度的强度大小来确定修正源项的大小,可以有效避免现有代数方法的不足。
-
公开(公告)号:CN109506807A
公开(公告)日:2019-03-22
申请号:CN201811421190.9
申请日:2018-11-26
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
Abstract: 本发明公开了一种稳态条件下的高温结构内部温度及壁厚同时测量方法,解决了稳态条件下的高温结构内部温度及壁厚无法同时测量的问题,该方法基于超声检测信号,将结构厚度和内部温度的同时测量转化为热传导问题热边界条件和结构厚度的多参数识别问题。其中,采用增大初始时间和时间步,使得瞬态求解得到的结构内部非均匀温度场与稳态结果保持一致,将稳态传热条件下的一个有效测量数据转化为瞬态条件下的多个有效测量数据,以有效解决多参数识别中输入信息不足带来的矩阵奇异问题。该方法通过求解热传导反问题,可快速、无损、非接触地测量稳态传热条件下相关的结构内部温度和厚度。
-
公开(公告)号:CN117236144A
公开(公告)日:2023-12-15
申请号:CN202311515846.4
申请日:2023-11-15
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/23 , G06F30/15 , G06F119/08
Abstract: 本发明公开了基于安装型面的正交各向异性防热材料导热主轴确定方法,涉及热传导数值计算领域,包括:获得安装型面离散单元的表示集合和离散单元对应的外法向量集合;建立位置与离散单元的外法向量之间的射线,获得射线上任意一点的第一表达式和离散单元上任意一点的第二表达式;基于第一表达式和第二表达式建立联立方程式,判断联立方程式是否存在物理解,若存在则表示该射线与离散单元有交点;更换离散单元返回获得第一表达式和第二表达式,直至所有离散单元均被遍历完成累计获得多个交点,基于交点与位置之间的距离大小以及位置获得防热结构导热主轴中的第一主轴方向,本发明解决了复杂防热结构正交各向异性防热材料导热主轴的确定问题。
-
公开(公告)号:CN116611173B
公开(公告)日:2023-10-03
申请号:CN202310869193.3
申请日:2023-07-17
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/15 , G06F30/28 , G06F113/08 , G06F119/08 , G06F119/14 , G06F119/12
Abstract: 本发明公开了一种多层级自适应耦合时间步长的飞行器累积热变形计算方法,包括以下步骤:S1、针对固定飞行器结构及材料,给定其长航时巡航计算状态,划分好飞行器流场计算网格G1和结构场计算网格R1;S2、进行t=0时刻气动力/热环境数据的计算求解;S3、在气动力/热环境数据基础上结合热壁热流修正方法开展第一层级的累积热变形计算,获得累积热变形的宏观变化特征;S4、根据该宏观变化特征,在温升变化剧烈的区域选择小的时间步,在温升变化缓慢的区域选择大的时间步开展累积热变形计算,获得新的热变形特征;S5、根据新的温升特征重新进行耦合时间步的选取,重复迭代开展高精度的累积热变形计算,直至热变形计算收敛。
-
公开(公告)号:CN116610905B
公开(公告)日:2023-09-22
申请号:CN202310892878.X
申请日:2023-07-20
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
Abstract: 本发明公开了一种基于各向异性尺度修正的反距离权重数据插值方法,包括以下步骤:针对固定的飞行器翼面或者舵面构型S,将已知物理量的数据点、待插值的数据点表示为集合;由厚度方向、展向和弦向三个方向构成插值坐标系;统计插值坐标系下已知物理量的数据点集合和待插值的数据点集合在厚度方向、展向和弦向三个方向上的最小值和最大值;在插值坐标系下根据最小值和最大值对已知物理量和待插值数据点进行各向异性修正得到数据点;在插值坐标系下进行插值;计算得到待插值点的物理量,用于飞行器薄的翼面、舵面结构进行多场耦合。本发明可以各向异性修正翼面和舵面各方向的坐标,增大薄层方向的距离权重系数,提高数据插值精度。
-
公开(公告)号:CN113158340B
公开(公告)日:2022-10-18
申请号:CN202110408435.X
申请日:2021-04-16
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/15 , G06F30/28 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明公开了一种针对k‑epsion湍流模型的湍流长度尺度修正方法,本修正方法以无量纲速度散度λl的值为基本自变量来确定修正源项的大小,通过控制函数tanh(h2(η‑h3))‑1实现了对修正源项作用区域的控制。本发明方法不依赖于壁面距离这一参数,而是根据流场中速度散度的强度大小来确定修正源项的大小,可以有效避免现有代数方法的不足。
-
公开(公告)号:CN112989485B
公开(公告)日:2022-04-12
申请号:CN202110065551.6
申请日:2021-01-18
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/15 , G06F119/08
Abstract: 本发明提供的一种基于热壁修正的沿弹道热流插值方法中,利用热壁修正公式,实现前后锚点壁面温度的统一,获得前后锚点统一壁温下的壁面热流,之后利用Fay‑Riddell求得的驻点热流值将前后锚点壁面热流无量纲化,再利用反时间距离权重法插值获得全部结构传热计算时刻点上的统一壁温下的无量纲壁面热流,最后利用对应结构传热计算时刻点上的驻点热流值进行有量纲化和利用热壁修正公式将统一壁面温度修正回各自的非均匀壁温,完成锚点壁面热流向结构传热计算时刻点壁面热流的插值过程。本发明考虑了弹道以及壁面温度变化对于热流插值结果的影响,能够减小热流插值误差,提高热流插值计算精度。
-
公开(公告)号:CN111931295A
公开(公告)日:2020-11-13
申请号:CN202010963871.9
申请日:2020-09-13
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/15 , G06F111/10 , G06F113/28 , G06F30/28
Abstract: 本发明公开了一种全弹道整体迭代的气动热/传热耦合计算方法。该方法迭代过程为重复热环境-传热的单向计算,操作简单易于实现;该方法按照一定规则在全弹道上进行锚点选取,锚点之间的热环境可并行计算,流场计算热环境和结构场计算的温度分布通过特定插值相互耦合求解,沿弹道整体迭代若干轮后计算精度可满足特定需求。相较于沿时间方向依次进行耦合计算的方法而言,全弹道整体迭代的气动热/传热耦合计算方法计算效率提升,提升倍数为锚点数/整体迭代轮数;由于沿全弹道整体迭代为误差减小过程,相较于误差增大过程的沿时间方向依次耦合计算,收敛误差至最小范围,计算精度提升。
-
公开(公告)号:CN111931295B
公开(公告)日:2024-08-20
申请号:CN202010963871.9
申请日:2020-09-13
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/15 , G06F111/10 , G06F113/28 , G06F30/28
Abstract: 本发明公开了一种全弹道整体迭代的气动热/传热耦合计算方法。该方法迭代过程为重复热环境‑传热的单向计算,操作简单易于实现;该方法按照一定规则在全弹道上进行锚点选取,锚点之间的热环境可并行计算,流场计算热环境和结构场计算的温度分布通过特定插值相互耦合求解,沿弹道整体迭代若干轮后计算精度可满足特定需求。相较于沿时间方向依次进行耦合计算的方法而言,全弹道整体迭代的气动热/传热耦合计算方法计算效率提升,提升倍数为锚点数/整体迭代轮数;由于沿全弹道整体迭代为误差减小过程,相较于误差增大过程的沿时间方向依次耦合计算,收敛误差至最小范围,计算精度提升。
-
公开(公告)号:CN117407635A
公开(公告)日:2024-01-16
申请号:CN202311358576.0
申请日:2023-10-18
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F17/11 , G06F30/20 , G06F119/08 , G06F113/08
Abstract: 本发明公开了一种基于结霜相似律的平板结霜厚度预测方法,涉及结霜行为预测领域,包括:基于霜层厚度和霜层平衡厚度建立无量纲霜层厚度的第一计算公式;基于结霜时间和结霜特征时间建立无量纲结霜时间的第二计算公式;将第一计算公式转化为第一关系式;将第二计算公式转化为第二关系式;获得第一低温平板在第一来流速度和第一低温平板长度下的霜层生长情况;利用第一关系式和第二关系式,基于第一低温平板在第一来流速度和第一低温平板长度下的霜层生长情况,预测第二低温平板在第二来流速度和第二低温平板长度条件下第二时刻的霜层厚度,本发明目的为减少低温平板上霜层厚度预测的计算量,提高预测效率。
-
-
-
-
-
-
-
-
-