首页 / 专利库 / 材料和表面特性 / 剪切粘度 / 一种风力机结冰数值模拟方法

一种机结数值模拟方法

阅读:3发布:2021-11-26

专利汇可以提供一种机结数值模拟方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及 风 力 机测试技术领域,公开了一种 风力 机结 冰 数值模拟方法,在动量-叶素理论模型 框架 下提出了基于粘性无粘耦合及浸入边界方法的覆冰风力机准三维计算模型;粘性无粘耦合方法对于附着流条件下的 翼型 扰流问题,具有计算效率高、计算结果准确的特点;基于 正交 笛卡尔网格的浸入边界方法,可以在避免网格重构的同时兼顾流动分离计算的准确性。本发明的覆冰风力机三维计算模型依托动量-叶素理论模型分析,将三维风力机结冰问题化解为沿展向不同 位置 的二维翼型结冰问题,在此 基础 上充分结合粘性无粘耦合方法及浸入边界方法的优势,实现风力机结冰的准确、高效模拟,为覆冰风力机的工程分析奠定重要的基础。,下面是一种机结数值模拟方法专利的具体信息内容。

1.一种机结数值模拟方法,其特征在于,包括如下步骤:
S1:获取风力机外形参数,及流场工况(入流风速、风向、风剪切剖面、湍流度等),通过动量叶素理论模型计算,获得光滑叶片沿展向不同二维翼型截面位置处的α及相对入流速度Vrel分布;
其中,V为平自由来流速度,Ω为风力机旋转速度,r为当前翼型截面的展向位置, 为当前翼型截面的扭角。a和b分别为轴向和周向诱导因子,可通过经典的动量-叶素理论模型计算得出;
S2:根据S1中计算出的当地二维翼型截面的入流风速和迎角,利用粘性无粘耦合模型计算当地二维翼型的流场解u、v:
其中μ为偶极子强度、σ为点源强度、r为速度势、n为几何形面的法向量,Sb和Sw固壁和尾涡几何形面,s为面积;
其中φ为流场中的速度势,
S3:判断S2中计算是否收敛,若收敛则利用二维翼型的流场解u、v开展后续计算,否则采用浸入边界数值计算模型,求解二维翼型的流场解u、v,计算公式为:
其中 为速度向量、t为时间、ρ为流体密度、p为压力、μ为流体动力粘度,f为固壁边界的作用力函数;
S4:基于S2或S3的计算结果,根据雾参数及翼型的几何外形,计算液滴在当前二维翼型截面内的运动轨迹,计算公式如下:
其中Cd为球形液滴假设条件下的阻力系数; 为液滴惯性参数;
为液滴相对速度雷诺数; 为液滴相对于空气的无量纲运动速度;
分别为无量纲形式的液滴x、y向的速度分量及无量纲时间;ud、vd为
液滴x、y向的速度分量;V∞为无穷远来流速度;t为时间变量;c为翼型弦长;ρa、ρd为空气和液滴的密度;g为重力加速度,α为当前翼型的迎角。
根据液滴轨迹计算结果可以求得当前翼型的收集系数β为
其中ds∞为无穷远来流中的相邻液滴在当前二维翼型截面内的几何间距,dsi为该相邻液滴撞击在翼型表面撞击点的几何曲线间距。
S5:根据S4的计算结果,利用经典Messinger模型进行当前二维翼型的结冰计算,获取翼型曲面上不同位置的结冰冻结因子fr。当前曲面位置的结冰厚度hi的计算公式为:
其中dli为当前计算单元的曲线长度;LWC为空气中的液态水含量;hc为对流换热系数;R为气体常数;Cp为空气的定压比热容;ρa为空气密度;Le为刘易斯数(Lewis Number);psat为饱和蒸气压;Ts和T∞为当前控制单元及无穷远来流的温度
S6:根据S5的计算结果,在获取了冰形厚度分布后,更新覆冰翼型的几何外形,并基于S2和S3中的方法进行覆冰翼型的流场计算,获得覆冰翼型的气动力系数,实现覆冰翼型的气动性能评估。
2.根据权利要求1所述的风力机结冰数值模拟方法,其特征在于:步骤S2中的浸入边界方法的控制方程计算时,设在计算空间Ω内,存在固体边界Γ,为满足Γ边界的无滑移条件,浸入边界模型的通过在控制方程为添加边界作用力源项,来模拟固壁边界与流场之间的相互作用;添加边界作用力源项的积分形式的N-S方程如下所示:
其中 为流场中速度向量,p为压力,ρ为流体密度,μ为动力粘度,为边界法向量; 为作用在固壁边界的虚拟作用力,用以表征固壁对流场的作用。

说明书全文

一种机结数值模拟方法

技术领域

[0001] 本发明涉及风力机测试技术领域,具体涉及一种风力机结冰数值模拟方法。

背景技术

[0002] 寒冷气候条件下的风能资源具有风能密度大、分布区域广等优势。截止到2011年,全球范围内寒冷气候条件下已探明的风能资源储量已高达到60GW。在我国,随着“三北地区”及七大风电基地的相继开发,风电市场已逐渐向寒冷气候特征严重的地区扩展。然而,在寒冷气候条件下风力机组覆冰问题成为影响发电量及机组安全的重要原因。
[0003] 覆冰计算通常包括四个模,即流场计算、液滴运动轨迹及收集率计算、结冰模型计算及固壁边界重构,其中流场计算是求解液滴输运特征的重要基础。目前,在绝大多数基于CFD流场解的结冰分析模型中,流场计算通常需要消耗较大的计算资源,影响结冰计算的效率。TURBICE和LEWINT是目前广泛适用于风力机覆冰计算的二维求解模型,其中LEWINT的核心模块为NASA开发的LEWICE模型。为了加速计算,TURBICE和LEWINT模型在求解流场时均采用了基于势流理论的二维面元方法(PanelMethod),这意味着当流场中存在较严重流动分离情况时,该数值模型会带来较大的计算误差。
[0004] 三维计算流体力学(CFD)的快速发展为分析复杂流场环境下的结冰问题提供了有效手段。为了提高风力机三维覆冰分析的准确性,不同的学者利用CFD方法直接求解风力机三维流场,获得了较精确的覆冰结果。然而,复杂固壁外形条件下网格重构困难,同时CFD较长的计算时间也是限制CFD技术在风力机覆冰工程应用中的主要瓶颈

发明内容

[0005] 基于以上问题,本发明提供一种风力机结冰数值模拟方法,实现风力机结冰的准确、高效模拟,为覆冰风力机的工程分析奠定重要的基础。
[0006] 为解决以上技术问题,本发明提供了一种风力机结冰数值模拟方法,包括如下步骤:
[0007] S1:获取风力机外形参数,及流场工况(入流风速、风向、风剪切剖面、湍流度等),通过动量叶素理论模型计算,获得光滑叶片沿展向不同二维翼型截面位置处的α及相对入流速度Vrel分布;
[0008]
[0009]
[0010] 其中,V为平自由来流速度,Ω为风力机旋转速度,r为当前翼型截面的展向位置,为当前翼型截面的扭角。a和b分别为轴向和周向诱导因子,可通过经典的动量-叶素理论模型计算得出;
[0011] S2:根据S1中计算出的当地二维翼型截面的入流风速和迎角,利用粘性无粘耦合模型计算当地二维翼型的流场解u、v:
[0012]
[0013] 其中μ为偶极子强度、σ为点源强度、r为速度势、n为几何形面的法向量,Sb和Sw固壁和尾涡几何形面,s为面积;
[0014]
[0015]
[0016]
[0017] 其中φ为流场中的速度势,
[0018] S3:判断S2中计算是否收敛,若收敛则利用二维翼型的流场解u、v开展后续计算,否则采用浸入边界数值计算模型,求解二维翼型的流场解u、v,计算公式为:
[0019]
[0020] 其中 为速度向量、t为时间、ρ为流体密度、p为压力、μ为流体动力粘度,f为固壁边界的作用力函数;
[0021] S4:基于S2或S3的计算结果,根据雾参数及翼型的几何外形,计算液滴在当前二维翼型截面内的运动轨迹,计算公式如下:
[0022]
[0023]
[0024] 其中Cd为球形液滴假设条件下的阻力系数; 为液滴惯性参数;为液滴相对速度雷诺数; 为液滴相对于空气的无量纲运动速度;
分别为无量纲形式的液滴x、y向的速度分量及无量纲时间;ud、vd为
液滴x、y向的速度分量;V∞为无穷远来流速度;t为时间变量;c为翼型弦长;ρa、ρd为空气和液滴的密度;g为重力加速度,α为当前翼型的迎角。
[0025] 根据液滴轨迹计算结果可以求得当前翼型的收集系数β为
[0026]
[0027] 其中ds∞为无穷远来流中的相邻液滴在当前二维翼型截面内的几何间距,dsi为该相邻液滴撞击在翼型表面撞击点的几何曲线间距。
[0028] S5:根据S4的计算结果,利用经典Messinger模型进行当前二维翼型的结冰计算,获取翼型曲面上不同位置的结冰冻结因子fr。当前曲面位置的结冰厚度hi的计算公式为:
[0029]
[0030]
[0031]
[0032] 其中dli为当前计算单元的曲线长度;LWC为空气中的液态水含量;hc为对流换热系数;R为气体常数;Cp为空气的定压比热容;ρa为空气密度;Le为刘易斯数(LewisNumber);psat为饱和蒸气压;Ts和T∞为当前控制单元及无穷远来流的温度
[0033] S6:根据S5的计算结果,在获取了冰形厚度分布后,更新覆冰翼型的几何外形,并基于S2和S3中的方法进行覆冰翼型的流场计算,获得覆冰翼型的气动力系数,实现覆冰翼型的气动性能评估。
[0034] 进一步地,步骤S2中的浸入边界方法的控制方程计算时,设在计算空间Ω内,存在固体边界Γ,为满足Γ边界的无滑移条件,浸入边界模型的通过在控制方程为添加边界作用力源项,来模拟固壁边界与流场之间的相互作用;添加边界作用力源项的积分形式的N-S方程如下所示:
[0035]
[0036]
[0037] 其中 为流场中速度向量,p为压力,ρ为流体密度,μ为动力粘度,为边界法向量;为作用在固壁边界的虚拟作用力,用以表征固壁对流场的作用。
[0038] 与现有技术相比,本发明的有益效果为:覆冰风力机三维计算模型依托动量-叶素理论模型分析,将三维风力机结冰问题化解为沿展向不同位置的二维翼型结冰问题,在此基础上充分结合粘性无粘耦合方法及浸入边界方法的优势,实现风力机结冰的准确、高效模拟,为覆冰风力机的工程分析奠定重要的基础。附图说明
[0039] 图1为实施例中风力机结冰数值模拟方法的流程图

具体实施方式

[0040] 为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
[0041] 实施例:
[0042] S1:获取风力机外形参数,及流场工况(入流风速、风向、风剪切剖面、湍流度等),通过动量叶素理论模型计算,获得光滑叶片沿展向不同二维翼型截面位置处的攻角α及相对入流速度Vrel分布;
[0043]
[0044]
[0045] 其中,V为水平自由来流速度,Ω为风力机旋转速度,r为当前翼型截面的展向位置,为当前翼型截面的扭角。a和b分别为轴向和周向诱导因子,可通过经典的动量-叶素理论模型计算得出;
[0046] S2:根据S1中计算出的当地二维翼型截面的入流风速和迎角,利用粘性无粘耦合模型计算当地二维翼型的流场解u、v:
[0047]
[0048] 其中μ为偶极子强度、σ为点源强度、r为速度势、n为几何形面的法向量,Sb和Sw固壁和尾涡几何形面,s为面积;
[0049]
[0050]
[0051]
[0052] 其中φ为流场中的速度势,
[0053] S3:判断S2中计算是否收敛,若收敛则利用二维翼型的流场解u、v开展后续计算,否则采用浸入边界数值计算模型,求解二维翼型的流场解u、v,计算公式为:
[0054]
[0055] 其中 为速度向量、t为时间、ρ为流体密度、p为压力、μ为流体动力粘度,f为固壁边界的作用力函数;
[0056] S4:基于S2或S3的计算结果,根据云雾参数及翼型的几何外形,计算液滴在当前二维翼型截面内的运动轨迹,计算公式如下:
[0057]
[0058]
[0059] 其中Cd为球形液滴假设条件下的阻力系数; 为液滴惯性参数;为液滴相对速度雷诺数; 为液滴相对于空气的无量纲运动速度;
分别为无量纲形式的液滴x、y向的速度分量及无量纲时间;ud、vd为
液滴x、y向的速度分量;V∞为无穷远来流速度;t为时间变量;c为翼型弦长;ρa、ρd为空气和液滴的密度;g为重力加速度,α为当前翼型的迎角。
[0060] 根据液滴轨迹计算结果可以求得当前翼型的收集系数β为
[0061]
[0062] 其中ds∞为无穷远来流中的相邻液滴在当前二维翼型截面内的几何间距,dsi为该相邻液滴撞击在翼型表面撞击点的几何曲线间距。
[0063] S5:根据S4的计算结果,利用经典Messinger模型进行当前二维翼型的结冰计算,获取翼型曲面上不同位置的结冰冻结因子fr。当前曲面位置的结冰厚度hi的计算公式为:
[0064]
[0065]
[0066]
[0067] 其中dli为当前计算单元的曲线长度;LWC为空气中的液态水含量;hc为对流换热系数;R为气体常数;Cp为空气的定压比热容;ρa为空气密度;Le为刘易斯数(LewisNumber);psat为饱和蒸气压;Ts和T∞为当前控制单元及无穷远来流的温度;
[0068] S6:根据S5的计算结果,在获取了冰形厚度分布后,更新覆冰翼型的几何外形,并基于S2和S3中的方法进行覆冰翼型的流场计算,获得覆冰翼型的气动力系数,实现覆冰翼型的气动性能评估。
[0069] 进一步地,步骤S2中的浸入边界方法的控制方程计算时,设在计算空间Ω内,存在固体边界Γ,为满足Γ边界的无滑移条件,浸入边界模型的通过在控制方程为添加边界作用力源项,来模拟固壁边界与流场之间的相互作用;添加边界作用力源项的积分形式的N-S方程如下所示:
[0070]
[0071]
[0072] 其中 为流场中速度向量,p为压力,ρ为流体密度,μ为动力粘度,为边界法向量;为作用在固壁边界的虚拟作用力,用以表征固壁对流场的作用。
[0073] 如上即为本发明的实施例。上述实施例以及实施例中的具体参数仅是为了清楚表述发明验证过程,并非用以限制本发明的专利保护范围,本发明的专利保护范围仍然以其权利要求书为准,凡是运用本发明的说明书及附图内容所作的等同结构变化,同理均应包含在本发明的保护范围内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈