首页/文章/ 详情

一个关于线圈加热的CFD仿真与实验对比的故事

1月前浏览2051

作为工程师,通常希望将模拟结果与实验结果进行比较,特别是在有限元分析(FEA)和计算流体力学(CFD)的学习过程中。追求的不是那些色彩斑斓、吸引眼球的图像和动画,而是真实可靠的模拟结果。有时,特定CFD问题的模拟结果的有效性难以自我验证。

测量自然对流引起的空气流速是一个有趣的研究课题。实验室中不打算通过点火或烧水来产生热源。这些热源产生的流动可能较大,便于使用热线风速计进行测量。然而,需要的是一个小型且可重复的实验设置。

焦耳热源是一个不错的选择。可以使用扬声器驱动器的音圈线圈,通过直流电作为热源。在没有风扇或空调引起的流动的“静态”房间空气中进行了测试。风速计探头距离线圈约1厘米,该仪表能够测量从0.01到25米/秒的速度,分辨率为0.01米/秒。但测试结果并不理想,热空气可能会干扰探头的热特性,导致无法进行CFD和比较。

自然对流的流动可能太小(<0.01米/秒)或者需要集中才能被仪表检测到。因此,探头需要远离热源。使用了一个纸板圆筒管,像烟囱一样,来约束流动。线圈位于管子的底部附近,风速计探头在管子的顶部(出口)测量流速。线圈表面和出口的温度由热电偶数据记录器测量。

从图表中可以看出,线圈表面温度达到了88.4摄氏度,出口空气温度约为37.6摄氏度。其他传感器用于测量环境温度。出口处的气流速度是0.33米/秒,这个读数是可靠的,因为它远高于仪表的分辨率和测量范围的最小值。现在有了CFD验证所需的实验数据。

为了简化CFD分析,目前不考虑焦耳加热,因为没有准确的线圈参数。在这种情况下,热生成的准确性将受到怀疑。线圈被假设为一个等温物体,其温度由测量的线圈表面温度给出。

模型的网格在下图中进行了说明。纸板管被简单地表示为无滑移壁面,并没有进行“物理”建模。线圈绕组作为一个铜实体CFD组件。“背景”流体是空气,CFD中假设了理想气体、可压缩流求解器和层流。

下面图表显示了CFD和实验数据。由于CFD中的线圈温度实际上是测量数据,两条曲线是相同的。出口温度测量和CFD数据非常吻合。放大的曲线图显示了测量中的缓慢上升,下降速率几乎相同。测量数据的移动平均值低于CFD数据。

在流速上观察到较大的偏差。测量和CFD数据的初始上升速率几乎相同,但随后测量数据的增加速率要低得多。CFD数据达到了稳态,而测量数据仍在上升。尽管如此,模拟和测量结果相差不大,它们处于同一数量级,这是一个令人鼓舞的结果。

实验和CFD之间的偏差可能由实验设置和测量误差造成。出口处的探头可能受到了较温暖的温度影响(在仪表的操作温度范围内),这可能对准确性产生了负面影响。另一个可能性是模拟中没有考虑纸板管的影响。使用的不正确的边界条件和参数也可能是造成差异的因素。这些与CFD相关的问题可以相对容易地解决,并应首先解决。

CFD分析中流体的流动、温度和密度下图中进行了展示。



来源:CFD饭圈
FluentCFX燃烧Polyflow理论科普ParaViewParticleWorks
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-15
最近编辑:1月前
CFD饭圈
硕士 分享CFD文章,感谢关注
获赞 22粉丝 21文章 376课程 0
点赞
收藏
作者推荐

【经典教材翻译】21-螺旋桨的叶片元素理论

预测螺旋桨(以及风扇或风力涡轮机)性能的一种相对简单的方法是使用叶片元素理论BladeElementTheory。在这种方法中,螺旋桨被分成沿长度的若干独立部分。在每个部分应用力平衡,涉及2D截面升力和阻力与该部分产生的推力和扭矩。同时,应用轴向和角动量的平衡。这产生一组非线性方程,可以通过迭代求解每个叶片部分。得到的各部分推力和扭矩的值可以累加以预测螺旋桨的整体性能。该理论不包括由螺旋桨脱落的尖端涡流引起的3-D流速或由螺旋桨旋转引起的角加速度引起的径向流分量等二次效应。与实际螺旋桨结果相比,该理论会过高预测推力,过低预测扭矩,导致理论效率比实际性能高出5%到10%。当叶片上的流动变得停滞或有相当一部分螺旋桨叶片处于风车配置状态而其他部分仍在产生推力时,一些流的假设也会在极端条件下失效。尽管存在上述限制,该理论在比较研究中非常有用,例如优化给定巡航速度的叶片俯仰设置或确定螺旋桨的最佳叶片实度。它仍然是在广泛的操作条件下获得推力、扭矩和效率的良好一级预测的最佳工具。1-叶片元素细分螺旋桨叶片可以如上图所示细分为离散的部分。如果假设每个部分只有轴向和角速度分量,并且来自其他部分的诱导流输入可以忽略不计,则可以独立分析每个部分的流动。因此,在上图所示的AA部分(半径=r),叶片上的流动将包括以下组成部分。其中V0是螺旋桨盘的轴向流,V2是角速度矢量,V1是部分局部流速矢量,向量V0和V2的总和。由于螺旋桨叶片将设置在给定的几何俯仰角(θ),局部速度矢量将在部分上产生攻角。可以使用标准的2-D翼型属性计算部分的升力和阻力。(注意:螺旋桨使用改变的参考线:零升力线不是部分弦线)。可以计算部分对推力和扭矩的贡献的升力和阻力分量。推力和升力方向之间的角度差定义为φ=θ−α因此,这个叶片元素的元素推力和扭矩可以写成用给定α的部分数据(CL和CD)替换,得到以下方程。每个叶片,其中ρ是空气密度,c是叶片弦长,这样叶片元素的升力产生面积就是c.dr。如果螺旋桨叶片的数量是B,那么,2-流入因素应用这一理论时,主要的复杂性出现在试图确定两个流分量V0&gt;和V2的大小时。V0大致等于飞机的前进速度(V∞),但由于螺旋桨自身的诱导轴向流进入尾流而增加。V2大致等于叶片部分的角速度(Ωr),但由于螺旋桨诱导的流动的旋转性质而略有减少。要准确计算V0和V2,必须应用轴向和角动量平衡来预测给定叶片元素上的诱导流效应。如下面的图表所示,诱导分量可以定义为增加或减少主要流分量的因子。通过AA部分的典型流管将具有速度因此,对于前一节流程图中显示的V0和V2的速度,其中a是轴向流入因子,其中b是角流入因子theangularinflowfactor(涡流因子swirlfactor)叶片部分的局部流速和攻角因此是3-轴向和角动量动量守恒可以对轴向和圆周方向应用动量流守恒的控制原理。对于轴向,从上游开始,通过螺旋桨在AA部分,然后进入尾流的流管,沿流管的动量变化必须等于这个叶片元素产生的推力。为了消除螺旋桨旋转引起的非稳态效应,使用的流管是覆盖螺旋桨叶片元素扫过的整个螺旋桨盘面积的流管,所有变量都假定为时间平均值。通过应用伯努利方程和动量守恒,对于管的三个单独部分,从自由流到盘的前面,从盘的后面到远下游的尾流以及平衡压力和面积与推力,可以证明盘处的轴向速度将是自由流和尾流速度的平均值。这意味着因此对于角动量通过结合角动量守恒和轴向速度变化,可以证明尾流中的角速度将是螺旋桨盘处的两倍。因此由于这些动量方程的最终形式仍然包含元素推力和扭矩的变量,它们不能直接用来求解流入因子。然而,现在存在一个非线性方程组(1),(2),(3),(4),(5)和(6),包含四个主要未知变量ΔT,ΔQ,a,b,因此可以对这个系统进行迭代求解。4-叶片元素理论的迭代求解过程叶片元素流动的求解方法是首先对流入因子a和b进行一些初始猜测。使用这些来找到叶片上的流动角(方程(3),(4)),然后使用叶片部分属性来估计元素推力和扭矩(方程(1),(2))。有了这些近似的推力和扭矩值,可以使用方程(5)和(6)来给出改进的流入因子a和b的估计。这个过程可以重复,直到a和b的值在指定的容差内收敛。应该注意到,这个非线性方程组的收敛并不是保证的。通常,应用一些收敛增强技术(例如Crank-Nicholson欠松弛)就可以得到结果,当使用线性翼型部分属性时。当使用非线性属性时,即包括失速效应时,获得收敛将更加困难。对于最终的流入因子a和b的值,将从方程(1)和(2)获得元素推力和扭矩的准确预测。螺旋桨推力和扭矩系数及效率通过累加所有径向叶片元素值,将获得整体螺旋桨推力和扭矩。然后可以计算非维数推力和扭矩系数以及它们被计算时的前进比。和其中n是螺旋桨的旋转速度,以每秒转数表示,D是螺旋桨直径。然后在这些飞行条件下螺旋桨的效率将是其中前进比J=V∞∕nD。邀您关注▽纯粹CFD:软件教程、行业应用、专业理论、基础科普、研究前沿、严选培训广告▽只聊CFD相关的大小事,信手天成,娓娓道来来源:CFD饭圈

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈