第五周：CFD偏微分方程部分+Ansys fluent 求解器设置学习周报

摘要

本周内主要学习了CFD中偏微分方程的基础概念和在Ansys fluent 中求解器设置的实例操作，对于偏微分方程在CFD中的应用和基本理论概念有了一定的了解，包括：何为特征线，线性单波方程的边界条件如何设定，如何使用特征线计算物理量，如何判断偏微分方程的类型等内容，并且对Ansys fluent的求解器进行了实际操作，明白了Ansys功能的强大。

Abstract

This week，I totally learned about the basic concept of partial differential equations in CFD，and the try to run the operation of solver settings in Ansys fluent.and gained a certain understanding of the application and basic theoretical concepts of partial differential equations in CFD, including: what is a eigenline, how to set the boundary conditions of a linear single-wave equation, how to use eigenlines to calculate physical quantities, how to judge the type of partial differential equation, etc. ,and the Ansys fluent solver was actually operated. Understood the power of Ansys.

1.常系数线性单波方程

  定义为某波以常速度c向右传播。波形与振幅均保持不变。其中x-ct=const为方程特征线。为自变量空间中一条曲线，在该曲线上物理量方程可以简化。

关于该偏微分方程的边界条件设置，见下图：

其核心还是在于关注流体在流场中的流动发展的影响，找出主要的影响地点。

2.特征线

特征线概念为在一个规定的特征方向上，在该方向内的数学性质，为在子空间中最基本的性质，即形状不发生改变，改变的只有模。（方向固定）可以认为在特征线上把函数变为矢量的过程。在特征线（面）上被算子作用后，只有长度发生变化，即只有拉伸，其基本的方向不会发生改变。利用“特征相容关系”可以在特征线上实现方程的简化，将偏微分方程转换为常微分方程。

2.1利用特征线计算物理量

（1）设定积分部长Δs

（2）在边界上选取初始点，由边界条件确定该点的物理量值u0

（3）根据特征线及特征相容关系数值积分，求出特征线下一个点的坐标和函数值，在进行递推计算，直到计算出整条特征线的（离散）坐标及物理量的（离散）值。

3.线性偏微分方程组

二阶偏微分方程确定类型时，尝试将其转化为一阶方程组进行判断：

关于欧拉方程组的学习可见下图：

4.使用“雷诺运输定律”进行守恒方程的推导（以质量守恒方程为例）

5.CFD常用方程

（1）不可压缩方程

（2）无粘方程，即欧拉方程

（3）Thinlayer-N-S方程

将交叉导数项去掉，由边界层性质为流动的法向剪切力为主导，粘性项中的流向和窄向剪切力影响较小。

（4）RANS（雷诺平均方程）

为湍流模式方程。

（5）PNS方程

将NS方程抛物化，进行简化计算。（边界层方程）

其中关于小扰动时的流动状态，还存在三个方程：扰动方程、小扰动方程、线化NS方程。

6.Ansys fluent 求解器设置操作

因本次只为在软件中熟悉求解器的界面操作和相关属性，故其设计参数均选择默认。求解器可分散三步进行：选择模型→指定介质参数→指定边界条件，完成这三步后即可开始进行计算

求解器分为三个部分：模型展示、选项界面、参数设置界面

将模型设置为湍流模型，流体为空气，边界条件及其余选项均为默认，

上图为残差值设置，每个物理量在计算进行迭代的过程中两次计算中的差值到达0.001时，认为已经收敛，停止计算。进行初始化后即可开始计算。在进行初始化后，设定迭代数为100，开始计算。

计算结果如下：

总结

在本周的学习中，对CFD的基础理论学习进行了丰富，了解到其一个基本的理论支持，学习了相关数学概念，让我在看资料视频时的阻力减少，CFD为交叉理论学科，需要补充的其他学科知识点还有很多，在以后的学习中我会针对性对数学理念进行补充学习。熟悉了Ansys fluent 求解器的设置及运行操作让我对Ansys的功能理解更进一步。