学习自 https://confluence.cornell.edu/pages/viewpage.action?pageId=333371302
FLUENT算例——Vertical Axis Wind Turbine (Part 1) 垂直轴风力机(第1部分)
以ANSYS 17.0为例。
该算例分为两个部分,第一部分将采用运动参考系(Moving Frame of Reference(MRF))来分析,第二部分将采用滑移网格(Sliding Mesh(SM))来作分析。
学习目标
- 获取通过垂直轴风力机的流动特性。
- 应用FLUENT中的运动参考系(MRF)和滑移网格(SM)的理念将非定常问题转化为定常问题求解。
问题描述
如上图所示,流过某垂直轴风力机(Vertical Axis Wind Turbine(VAWT))的均匀来流速度为V=10m/s,VAWT直径为12cm且有3个均布叶片,每个叶片弦长为2cm。为简便计,假设其以40RPM的恒定角速度做旋转运动,每个叶片的中心到轮毂中心的距离为0.04m。
实际上,这是对上图所示的三维实际风机的简化二维计算模型。
注意,这是达里厄型垂直轴风力发电机(Darrieus VAWT),其属于升力型风力机。与其对应的是萨沃纽斯型垂直轴风力发电机(Savonius VAWT),这是阻力型风力机。感兴趣的读者可以参考风力机的相关资料。
注意
这里咱们假设风力机的转动是和流体无关的,这和实际情况并不相符,因为流动实际上推动了风力机的旋转运动,即,对于确定的几何形状和质量(转动惯量)的风力机而言,对应外部来流速度的只有唯一的稳定转速。然而,数值模拟由流动所导致的几何体的运动是非常困难的事情,而且需要用到“6DOF solver”(6自由度求解器),首先未必能算,其次算得很慢,再次算得未必靠谱。所以呢,咱们这儿就简单点,同时给定风力机旋转速度和来流的速度,先算个稳态的问题(part1-MRF),再算个瞬态的问题(part2-SM)。
1 预分析和准备工作
预分析
- 数学模型:关于控制方程、边界条件、数学模型中的相关假设;
- FLUENT中的分析过程:对FLUENT中的求解过程做简要概述;
- 手动计算预期结果:对计算结果的展示与后处理,得出风力机性能的相关数据。
数学模型
控制方程
控制方程为运动参考系下的质量守恒和动量守恒(Navier-Stokes方程)。
质量守恒
在以恒定角速度旋转的运动参考系下的Navier-Stokes方程
其中
FLUENT将求解在运动参考系下的上述方程,这个近似所带来的好处就是咱们不再需要处理运动网格了。
注意咱们要求解湍流问题,那么就用 k − ϵ k-\epsilon k−ϵ模型好了。
边界条件
创建的流体计算域应该比风力机的几何结构大数倍才好,该计算域正是风力机所扰动(波及,影响)到的流场区域。这里用一个外部大圆来作为整个计算域,实际上随便用个几何形状来做为“远场”区域都可以,但是用圆的话相对比较简单些。
边界条件为:
- Inlet(far-field):速度进口,x方向速度恒定10m/s,湍流强度5%,湍流粘性比为1。
- Outlet(far-field):压力出口,绝对压力101325Pa,或1个大气压。
- Blades:壁面,速度为0(无滑移边界条件)。
FLUENT中的分析过程
FLUENT将采用有限体积方法求解,先把控制域划分成多个控制体(或单元),然后在每个单元上去积分控制方程,保证每个单元内质量、动量、能量的生成量和通过边界流入的通量是相互平衡的,从而将原本的偏微分控制方程组转化为控制单元的非线性代数方程组,再将其转化为线性化方程组,接下来用迭代方法去求解这些方程组,直到残差满足收敛指标为止。
速度、压力、角速度和湍动能k将参与计算,而其他的值可以通过后处理由这些值推出来,比如剪切力、涡量场等。
CFD的理论相当复杂,具体算法可参考计算流体动力学的相关书籍。
手动计算预期结果
叶尖速比Tip Speed Ratio (TSR)
TSR由下式定义:
λ = 叶 尖 处 的 速 度 来 流 风 速 = r × ω V \lambda=\frac{叶尖处的速度}{来流风速}=\frac{r\times\omega}{V} λ=来流风速叶尖处的速度=V