计算流体力学内核.pptVIP

计算流体力学内核王建文 有限单元法——有限元法的来源 有限元法的基本思想来源于Courant在1943年发表的论文。他尝试应用在一系列三角形区域上定义的分片连续函数和最小位能原理相结合，来求解St. Venant扭转问题。 Turner和Clough在1956年将桁架分析中的位移法推广到弹性力学平面问题，用于飞机结构的分析。并首次给出了用三角形单元求解平面应力问题的正确解答。 Clough在1960年第一次提出了有限元法（Finite Element Method, FEM)这个名称，使人们更清楚地认识到有限单元法的特性和功效。 有限单元法——有限单元法的要点 （1）将一个表示结构或连续体的求解域离散为若干个子域（单元），并通过它们边界上的结点相互联结成为组合体。 （2）用每个单元内所假设的近似函数来分片地表示求解域内待求的未知场变量。而每个单元中的近似函数本身就是由未知场函数在单元结点上的数值和与其对应的插值函数来表示。使原来待求场函数的无穷多自由度问题转换为求解场函数结点值的有限自由度问题。 （3）通过和原问题数学模型等效的变分原理或加权余量法，建立求解基本未知量（场函数的结点值）的代数方程组或常微分方程组。 有限单元法——微分方程的等效积分形式 例：三维稳态热传导方程 有限单元法——加权余量法 在求解域Ω中，若场函数是精确解，则满足A（u）= 0和B（u）= 0，此时同样满足等效积分形式。 但是对于复杂的实际问题，精确解往往很难找到，因此需要找到具有一定精度的近似解。 对于微分方程和边界条件所表达的物理问题，假设未知场函数u可以采用近似函数来表示。该近似函数可以是分片有效的。 近似函数是一族带有待定参数的已知函数，一般形式为： 有限单元法——加权余量法 在通常n取有限项的情况下，近似解是不能精确满足微分方程和边界条件的，它们将产生残差。 有限单元法——加权余量法 配点法：强迫余量在域内n个点上等于零。 计算流体力学——通用求解理论 量纲 控制方程 浮力 多组分流 附加变量 旋转力 源项 边界条件 自动时间步计算 自适应网格 微孔介质流 计算流体力学——通用求解理论 量纲：长度（L）、质量（M）、时间（T）、温度（Θ）和化学当量（A）。 计算流体力学——通用求解理论 计算流体力学——通用求解理论 计算流体力学——通用求解理论 计算流体力学——通用求解理论 计算流体力学——通用求解理论 计算流体力学——通用求解理论 计算流体力学——通用求解理论 计算流体力学——通用求解理论 5 湍流 湍流是流体微元在时空上的涨落。 湍流是流动不确定性的体现。 湍流是惯性力大于粘性力的高雷诺数流动。 在理论上Ｎ－Ｓ方程既可以描述层流，也可以描述湍流。但是湍动的尺度和微元体是同样的数量级，但是直接模拟法（Direct numerical simulation）需要占用大量的计算空间，现有计算机无法完成。 现有算法只能简化，提出多种湍流模型，使湍流表现出统计特性。 5 湍流—雷诺数划分 5 湍流—速度表示法 湍流平均特性：时间上，用大于湍流涨落的时间观察湍流，湍流流动表现出平均特性。 湍流的涨落特性：瞬时变化与平均特性的差值。 5 湍流—湍流模型 湍流模型基本都是通过引入时均量和脉动量，来修正初始非稳态N-S方程，形成雷诺平均N-S方程。 湍流模型的实质是获得湍流的平均量，而不必求解脉动量。 所有尺度的湍流场都可以通过模型来模拟。 用模型来模拟湍流场可以减少计算量，这是模型法的优点。 5 湍流—雷诺平均模型 5 湍流—湍流涡粘模型 涡粘理论认为湍流是由不断产生和消散的小涡组成，在这些小涡中雷诺应力和速度梯度和涡粘性成正比。 5 湍流—湍流涡粘模型 把涡粘假设带入雷诺平均模型： 5 湍流—湍流涡粘模型 涡粘模型根据涡粘度和涡散系数的表述方式分为不同的类型。 零方程模型是最简单的湍流涡粘方程，直接由平均速度计算全局涡粘系数，并根据经验公式计算长度尺寸，而不求解雷诺输运方程。 5 湍流—湍流涡粘模型 两方程模型分别求解速度输运方程和长度输运方程。是计算效率和计算精度之间的良好折中。 κ-ε和κ-ω两方程模型采用梯度扩散假设建立雷诺应力与平均速度梯度及湍流粘度之间的关系，即湍流粘度看作是湍流速度和湍流长度尺度的生成项。 两方程模型中，湍流速度通过湍动能方程（输运方程）求解，湍流长度通过湍流场的两个特性（通常是湍动能和湍动能耗散率）求解。 5 湍流—湍流涡粘模型 CFX中的κ-ε方程： κ是湍动能，源于湍流场中的速度波动，单位为m2/s2. ε是湍动能耗散率，源于湍流场中的速度波动耗散，单位为m2/s3. 5 湍流—湍流涡粘模型 5 湍流—湍流涡粘模型 CFX中的RNGκ