第八章 边界条件

合集下载

boundary condition 条件

boundary condition 条件

boundary condition 条件引言:在数学和物理学中,边界条件是问题解决过程中非常重要的一部分。

它们定义了问题所在的领域,并确定了问题的解决方案。

本文将详细介绍边界条件的概念以及在不同领域中的应用。

1. 什么是边界条件?边界条件是指问题定义的区域边界上的条件或限制。

在数学中,边界条件是用来约束问题解的特定点或区域的条件。

在物理学中,边界条件是用来限定问题中参与运算的物理系统与其周围环境之间的相互作用。

2. 数学中的边界条件在微积分和偏微分方程中,边界条件用来限定定义域。

例如,在求解一维热传导问题时,可以通过指定热量的输入和输出来定义系统的边界条件。

常见的边界条件包括固定边界条件(温度或导数固定)和自由边界条件(热量或能量流量固定)。

3. 物理学中的边界条件在物理学中,边界条件将物理系统与其周围环境进行联系。

例如,在流体力学中,可以通过指定壁面上的速度或压力分布来定义流体的运动。

对于静电场问题,可以通过指定电势值或电场强度来定义电荷的分布。

4. 工程中的边界条件在工程领域中,边界条件用于模拟和优化各种系统的功能和行为。

例如,在结构工程中,边界条件可以用来模拟外部加载(例如风载荷、地震力等)对建筑物的影响。

在电气工程中,边界条件可以用来模拟电流和电压在电路中的传输和分布。

5. 边界条件对解的影响边界条件的选择和应用会对问题的解产生重要影响。

不同的边界条件可以导致不同的解,从而得到不同的结果。

因此,正确选择和应用边界条件是问题求解过程中的一个关键步骤。

6. 边界条件的设置方法在实际问题中,确定边界条件可能并不总是直观或容易的。

一种常用的方法是根据问题的物理意义和要求来选择合适的边界条件。

此外,使用数值方法也可以帮助确定边界条件。

通过将问题离散化为有限元或有限差分网格,并将边界条件应用于离散化的边界上,可以有效地求解复杂的问题。

结论:边界条件是定义问题域和约束解的一种方法。

无论是在数学、物理还是工程领域,正确选择和应用边界条件对于解决问题和获得准确结果至关重要。

《边界条件教程》课件

《边界条件教程》课件
总结词
Dirichlet边界条件是一种常见的边界条件,它指定了函数在边界上的值。
详细描述
在解决偏微分方程时,常常会遇到各种边界条件。其中,Dirichlet边界条件规定 了函数在边界上的取值,即要求函数在边界上达到特定的值。这种边界条件通常 用于控制流动、热传导等问题,以确保物理现象的合理性和实际意义。
Neumann边界条件
总结词
Neumann边界条件规定了函数在边界上的导数值。
详细描述
与Dirichlet边界条件不同,Neumann边界条件关注的是函数在边界上的导数。这种边界条件通常用于描述物理 现象的流出或流入,例如流体流动、热传导等。在解决偏微分方程时,Neumann边界条件可以确保物理量的连 续性和自然边界条件。
在有限差分法中实现边界条件
1 2 3
反射边界条件
在有限差分法中,对于反射边界,可以通过设置 边界上的网格点与相邻网格点的物理量相等来实 现。
吸收边界条件
对于吸收边界,可以通过设置边界上的网格点物 理量与相邻网格点物理量相同,但方向相反来实 现。
周期性边界条件
对于周期性边界条件,可以通过设置边界上的网 格点物理量与相邻网格点物理量相同来实现。
解的误差分析
评估边界条件对解的误差的影响,了解误差来源和误差传播机制。
解的敏感性和鲁棒性
分析边界条件对解的敏感性和鲁棒性的影响,了解解的稳定性和可 靠性。
05 边界条件的实际应用
在流体动力学中的应用
总结词
描述边界条件在流体动力学中的重要性及应用。
详细描述
在流体动力学中,边界条件是描述流体与固体边界相互作用的关键因素。它们 决定了流体在边界上的行为,如流动速度、压力和温度等。边界条件的应用范 围广泛,包括航空航天、船舶、汽车和能源等领域。

边界条件解析课件

边界条件解析课件
边界条件解析课件
目录
CONTENTS
• 边界条件概述 • 常见边界条件类型解析 • 边界条件的设置与调整 • 边界条件在各领域的应用 • 边界条件的未来发展与展望
01 边界条件概述
定义与分类
定义
边界条件是指在求解数学问题时,对所求问题的定义域或解的取值范围进行限 制的一些条件。
分类
根据不同的分类标准,边界条件可以分为多种类型,如根据约束类型可分为显 式和隐式边界条件,根据对解的影响可分为本质边界条件和非本质边界条件等。
详细描述
在物理、工程和自然界中,很多现象具有周期性变化的特性。 周期性边界条件是指在系统的边界上,变量的值按照一定的 时间或空间周期重复。这种边界条件通常用于模拟具有周期 性变化的现象,例如振动、波动等。
反射性边界条件
总结词
反射性边界条件是指系统在边界上的变量值被反射回来,不穿过边界。
详细描述
反射性边界条件通常用于模拟具有反射特性的现象,例如声音在硬表面上的反射、 光线在镜面上的反射等。在边界上,变量的值被反射回来,不穿过边界,保持了 系统的完整性。
总结词
在数据分析中,边界条件可以用来限制数据范围和分 析结果。
详细描述
在数据分析中,边界条件可以用来限制数据范围和分析 结果。例如,在统计分析中,边界条件可以用来排除异 常值和离群点,以确保分析结果的准确性和可靠性。此 外,在数据挖掘中,边界条件也可以用来限制搜索范围 和优化算法性能。
05 边界条件的未来发展与展 望
详细描述
常见的问题包括边界条件的设定不正确、不 合理或不符合实际情况等。针对这些问题, 可以采取相应的解决方案,如重新审查物理 定律和数学原理、增加额外的数据或信息来 源、寻求专业人士的帮助等。同时,应注意 及时总结和归纳经验教训,不断完善边界条 件的设置和调整过程。

流体力学第八章答案

流体力学第八章答案

流体力学第八章答案【篇一:流体力学第8、10、11章课后习题】>一、主要内容(一)边界层的基本概念与特征1、基本概念:绕物体流动时物体壁面附近存在一个薄层,其内部存在着很大的速度梯度和漩涡,粘性影响不能忽略,我们把这一薄层称为边界层。

2、基本特征:(1)与物体的长度相比,边界层的厚度很小;(2)边界层内沿边界层厚度方向的速度变化非常急剧,即速度梯度很大;(3)边界层沿着流体流动的方向逐渐增厚;(4)由于边界层很薄,因而可以近似地认为边界层中各截面上压强等于同一截面上边界层外边界上的压强;(5)在边界层内粘性力和惯性力是同一数量级;(6)边界层内流体的流动与管内流动一样,也可以有层流和紊流2种状态。

(二)层流边界层的微分方程(普朗特边界层方程)??v?vy?2v1?p?vy?????vx?x?y??x?y2????p??0?y???v?vy???0?x?y??其边界条件为:在y?0处,vx?vy?0 在y??处,vx?v(x)(三)边界层的厚度从平板表面沿外法线到流速为主流99%的距离,称为边界层的厚度,以?表示。

边界层的厚度?顺流逐渐加厚,因为边界的影响是随着边界的长度逐渐向流区内延伸的。

图8-1 平板边界层的厚度1、位移厚度或排挤厚度?1?1?2、动量损失厚度?2?vx1?(v?v)dy?(1?)dy x??00vv?2?1?v2???vx(v?vx)dy???vxv(1?x)dy vv(四)边界层的动量积分关系式??2???p?vdy?v?vdy?????wdx xx??00?x?x?x对于平板上的层流边界层,在整个边界层内每一点的压强都是相同的,即p?常数。

这样,边界层的动量积分关系式变为?wd?2d?vdy?vvdy?? x?x??00dxdx?二、本章难点(一)平板层流边界层的近似计算根据三个关系式:(1)平板层流边界层的动量积分关系式;(2)层流边界层内的速度分布关系式;(3)切向应力关系式。

第八章 边界条件

第八章  边界条件

第八章边界条件任何数值模拟都可以认为仅仅是在物理区域或系统的一部分中进行的。

区域的断层产生了人工边界,在这个断层中有我们处理的物理量。

此外,还有暴露在流体中的自然边界。

边界条件的数值处理需要特别注意。

在实际的系统中处理不当模拟就会出现偏差。

与此同时,稳定性和求解方案中的合成速度同样对数值模拟有消极的影响。

下边边界条件的类型是我们在欧拉方程和N-S方程中数值计算最常见的几种:·固体壁面·外表面的远场和流体内部流出或流出的表面·对称面·平整切割和周期性边界。

·平板间的边界这些边界条件的处理问题在以后几节中会进行详细的介绍。

对于文献中进一步涉及的边界条件,比如壁面上的热辐射或者是自由表面上的(热辐射),读者可以在3.4节中了解。

8.1 虚拟单元的概念在我们讨论边界条件时,我们需要提到虚拟单元(也可以被称作虚拟点)这个概念。

在规则的网格中这种方法非常的流行。

然而,在不规则网格中,虚拟单元仍然有很多的优点。

虚拟单元是在物理区域外部附加层上的一些网格点。

这个可以由图8.1中的二维规则网格中看到。

正如我们看到的,整个计算区域被两层虚拟单元包围着(由虚线标出),虚拟单元(点)通常不会像区域内的网格一样产生(除过多平板的网格)。

尽管它仍然有几何形状,比如体积或者表面的矢量,但是它仅仅是虚拟的。

利用虚拟单元可以简化计算沿边界的通量,梯度,散度等等。

这是由于在边界上可以将空间离散的模型进行扩展。

正如图8.1中我们看到的,在物理区域内同样可以进行离散。

因此,我们可以在所有的“物理”网格点中求解控制方程。

这种方法可以使离散工作非常简单。

此外,所有规则的网格点可以存在在一个单独的区域内,这在矢量计算中非常很有用。

虚拟单元不但包含有守恒变量,同时也有几何量。

很明显的是,虚拟单元层必须完全覆盖物理区域外。

几何量通常由边界的控制体积来求得。

在多网格平板中(3.1节),所有的流体变量和几何变量可以从相邻的平板求得。

弹性力学-边界条件

弹性力学-边界条件
y
yx

x
P y
fx
n
l cosn, x cos m cosn, y sin
xy

由 x s m xy s f x xy s m y s f y
fy
x s cos yx s sin 0

h 2 h 2
h 2 h 2
f x ydy M
则边界条件可以写成(P.23 (b)):

x x l
dy Fx ,

xy x l
dy Fy ,

x x l
ydy M
悬臂梁的例子:
y
h 2 h 2
y y x
h 2 h 2
x
P
L
L
对边界条件的积分为: (P.23 (b)):
x yx
xy l fx y s m fy
x 上面:l=0,m=-1 左面: 右面: l=-1 l=1 m=0 m=0 下面:l=0,m=1 y
(2).上下两面 ( ) f l 0 m 1 ( ) f
二、应力边界条件 在边界上的楔形体(单位厚度)如图所示: 弹性体内单元体斜面上的 y 应力分量与坐标面应力的 yx 关系有(静力平衡) f xn X x p x x xy l p y m y yx
• 所得到的应力分量必须在所有边界上各点处严 格满足应力边界条件,才是所论问题的解答。 • 在小边界上,如果不能严格满足边界条件,可 以用圣维南原理在静力等效意义上满足(积分 意义上的)边界条件。 • 根据这个原理:两组面力其分布尽管不同,但 如果两者的合力与合力矩相同(静力等效),此 时它们所产生的作用结果仅仅在局部有比较大 的差异,远离这个局部,结果基本相同。

边界条件定义

边界条件定义

边界条件定义边界条件是指在一个问题或系统中,所设定的特定条件或限制,用于测试或确定系统的行为。

边界条件在各个领域都有应用,例如软件开发、数学、物理学等等。

在软件开发中,边界条件是指在测试或运行程序时,需要考虑的各种极端情况。

下面将从不同领域的角度,讨论边界条件的定义和应用。

在数学中,边界条件是指在函数或方程中所设定的特定约束条件。

例如,在求解微分方程时,需要给定初始条件或边界条件,以确定唯一的解。

边界条件可以是函数在某一点的值,或者函数在某一区间的行为。

这些边界条件的设定,对于解的存在性和唯一性具有重要影响。

在物理学中,边界条件指的是在物理系统中所设定的限制条件。

例如,在求解波动方程时,需要考虑波函数在边界处的行为。

边界条件可以是波函数在某一点的值,或者波函数在某一区域的导数。

这些边界条件的设定,对于确定系统的行为和性质具有重要影响。

在计算机科学中,边界条件是指在程序设计或算法实现中所设定的特定限制条件。

例如,在编写排序算法时,需要考虑数组的边界条件,即数组的起始位置和结束位置。

边界条件的设定,可以避免数组越界和程序崩溃的情况发生。

边界条件的考虑也可以提高程序的效率和性能。

除了数学、物理学和计算机科学,边界条件在其他领域也有广泛的应用。

例如,在经济学中,边界条件是指经济模型中所设定的特定限制条件,用于分析和预测经济现象。

在生物学中,边界条件是指生物系统中所设定的特定约束条件,用于研究生物过程和现象。

边界条件的设定需要考虑到问题的特性和目标,以及系统的实际情况。

边界条件的选择应该合理、准确,能够准确反映问题的本质和复杂性。

同时,边界条件的设定也需要符合问题的要求和约束,以保证系统的稳定性和可靠性。

边界条件是问题或系统中所设定的特定条件或限制,用于测试或确定系统的行为。

边界条件的设定在各个领域都有重要的应用,对于解决问题和研究系统行为具有关键作用。

边界条件的设定需要考虑问题的特性和目标,以及系统的实际情况,从而保证系统的稳定性和可靠性。

第八章 非齐次边界条件处理(3节).

第八章 非齐次边界条件处理(3节).
§8.3非齐次边界条件的处理
本节中心内容
非齐次边界转化为齐次边界的问题;
本节基本要求
掌握非齐次边界齐次化的方法 着重掌握求解四种非齐次边界问题的解题思 路、解题步骤。 掌握求解非齐次边界问题的特殊方法
以前处理方程都是对齐次边界条件,而生活实践中大多数 是对于非齐次边界条件如何处理? 一般处理方法是要经过 代换转换为齐次的。
2 0 1 ( 0 0 ) l l
l
2 Cn { 0 ( x) ( x 0) 0 ( x) [ x (l o)]} l 0 n x con dx l 2 2 n [ 0 0 cos n ] 0 [1 () ] l l 0 (n为奇数) 4 0 (n为偶数, n 2k , k 1, 2,3,...) l
四、
X ( x) X ( x) 0 X (0) 0 X (l ) 0
1 n 2 x (n 0,1, 2,) X n ( x) Cnco s l
其本征值和本征函数分别为
1 2 2 (n ) 2 l2
五、
( ) ( ) 0 ( 2 ) ( )
一、齐次化的一般方法:
1、第一种边界条件:
x 0 , x l两端都是第二类
非齐次边界条件
utt a 2u xx 0
u x0 (t ), u xl (t )
u
y 0
( x), u t
t 0
( x)
此时,边界条件为非齐次的。选取一个函数
v( x, t )
( x ) t
t 0
尽管
的方程一般是非齐次的, w( x, t )
u(t ) u xl (t )
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第八章边界条件任何数值模拟都可以认为仅仅是在物理区域或系统的一部分中进行的。

区域的断层产生了人工边界,在这个断层中有我们处理的物理量。

此外,还有暴露在流体中的自然边界。

边界条件的数值处理需要特别注意。

在实际的系统中处理不当模拟就会出现偏差。

与此同时,稳定性和求解方案中的合成速度同样对数值模拟有消极的影响。

下边边界条件的类型是我们在欧拉方程和N-S方程中数值计算最常见的几种:·固体壁面·外表面的远场和流体内部流出或流出的表面·对称面·平整切割和周期性边界。

·平板间的边界这些边界条件的处理问题在以后几节中会进行详细的介绍。

对于文献中进一步涉及的边界条件,比如壁面上的热辐射或者是自由表面上的(热辐射),读者可以在3.4节中了解。

8.1 虚拟单元的概念在我们讨论边界条件时,我们需要提到虚拟单元(也可以被称作虚拟点)这个概念。

在规则的网格中这种方法非常的流行。

然而,在不规则网格中,虚拟单元仍然有很多的优点。

虚拟单元是在物理区域外部附加层上的一些网格点。

这个可以由图8.1中的二维规则网格中看到。

正如我们看到的,整个计算区域被两层虚拟单元包围着(由虚线标出),虚拟单元(点)通常不会像区域内的网格一样产生(除过多平板的网格)。

尽管它仍然有几何形状,比如体积或者表面的矢量,但是它仅仅是虚拟的。

利用虚拟单元可以简化计算沿边界的通量,梯度,散度等等。

这是由于在边界上可以将空间离散的模型进行扩展。

正如图8.1中我们看到的,在物理区域内同样可以进行离散。

因此,我们可以在所有的“物理”网格点中求解控制方程。

这种方法可以使离散工作非常简单。

此外,所有规则的网格点可以存在在一个单独的区域内,这在矢量计算中非常很有用。

虚拟单元不但包含有守恒变量,同时也有几何量。

很明显的是,虚拟单元层必须完全覆盖物理区域外。

几何量通常由边界的控制体积来求得。

在多网格平板中(3.1节),所有的流体变量和几何变量可以从相邻的平板求得。

图8.1中灰色的虚拟单元产生了一个问题,由于不清楚怎么设置它的值(如果没有相邻网格板),通过环绕形的离散模型,我们不需要知道它的值。

但是,在计算梯度(粘性通量—见4.4节)或者在多网格中计算转移量时就很有比必要知道它的值。

通常,如图8.1中用箭头表示的一样,相邻的“规则的”虚拟单元的平均值是很有用的。

8.2 固体壁面8.2.1非粘性流动在非粘性流动中,流体流过表面,由于没有摩擦力,速度矢量与壁面相切,与就是说在壁面的法线方向没有作用力。

即:0=⋅n v 在壁面上, (8.1)n 表示壁面上的单位法相矢量。

因此,相对应的速度V 等于0。

即,对流通量的矢量简化成只有压力项:⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=0 0 )(w z w y w x c p n p n p n F ω (8.2)其中Pw 为壁面压力。

规则的单元中心方案在以单元中心的方案中,在单元的质心可以求出压力。

但是,在等式(8.2)中,边界表面单元的Pw 需要求解出。

我们可以很容易的通过内部区域推导出壁面的压力。

考虑到图8.2,我们可以简单的设Pw=P2。

通过任一个两两相连点,我们可以取得较高的精确度。

)3(2132p p p w -=(8.3)或者一个三点相连的推导公式是)31015(81432p p p p w +-=(8.4)为了解释网格的延伸,到壁面的距离可以替代常数系数。

【1】上式的推导公式8.3和8.4没有涉及到网格和表面的几何形状。

一个可供选择的方法—Rizzi 发明的被称作为法向动量的关系。

它是基于在非粘性流动中,壁面是流线型的。

在式8.1中沿着流线型壁面的法线流动方向的导数为0。

并且动量方程可以替换成如下的公式:p n n v v ∇⋅=∇⋅⋅)(ρ (8.5)等式(8.5)包括密度,速度和壁面法线方向的压力。

法向动量在【1】中给出了精确的结果。

但是在复杂的几何表面,法向动量的数值求解存在着问题。

推导中详细的描述和精确的比较可以在文献【1】中找到。

虚拟单元中守恒变量的值在壁面内部能够被线性的推导出:320321232W W W W W W -=-= (8.6)公式(8.6)中的参数和图8.2中的相对应。

如果虚拟单元在空间上离散后可以被利用,那么对流通量的计算就可以和公式(8.2)一致。

规则的单元顶点的方案通过重叠的控制体积(见4.2.2),对单元顶点的离散可以直接由方程(8.1)在边界条件上应用。

壁面的对流通量可以通过公式(8.2)来求解。

壁面的压力Pw 可以通过计算平均节点的值来表示,二维公式节点的值可以通过(4.26)计算,三维公式节点的值可以通过(4.27)计算。

分类公式(二维公式中的(4.30))现在仅仅解释了两个单元(在三维中是4个)。

在重叠的控制体积内,单元顶点的值可以通过许多不同的方法求解。

一种方法是通过公式(8.2)将壁面上每个表面的控制体积分离。

因此,通过图8.3,我们可以写出:⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=+++------0 )()()()()()(0 0 )()()()()()(0 )(2,4/12,2,4/12,2,4/12,2,4/12,12,4/12,12,4/12,12,,i w i z i w i y i w i x i w i z i w i y i w i x i w c p n p n p n p n p n p n F (8.7)公式(8.7)中的压力2,4/1)(-i w p 和2,4/1)(+i w p 可以通过线性的插值表示:[]2,12,2,4/1)()(341)(+++=i w i wi w p pp (8.8)对应的三维公式可以表示未划分的网格。

另外一种可能的方法是在各自的壁面节点利用公式(8.2)来进行计算。

壁面的压力Pw 可以简单的设它等于P2。

单位法向矢量可以通过节点2上所有的面的法向矢量的平均值来表示。

这种方法要求对速度矢量进行修正。

通过时间步进的求解,壁面上的速度矢量被投影到切平面。

【3】【4】2,2,2,2,2,)(])([i av i av i i corri n n v v v ⋅⋅-=)( (8.9) av n 是平均单位法向矢量。

通过这种方法,流动将会变成与壁面相切。

为了对虚拟点进行赋值,在内部流场中,利用与公式(8.6)相关的公式推导守恒变量是非常有效的。

不规则的以单元为中心的方案公式(8.1)中壁面的边界条件也同样适用于未构建以单元为中心的方案。

如果边界单元是四面体、六面体或者是棱形(壁面上有三角形的表面),压力可以通过公式(8.3)来推导。

相邻的单元(图8.2中的3号)可以通过5.2.1中表面的数值来代替。

对于三角形或者四面体单元,在文献【5】【6】中通过一层虚拟单元来表示。

壁面上边界单元的速度矢量可以表示虚拟单元中速度分量。

比如:在图8.2中的虚拟单元1速度可以如下来表示:n V v v 2212-= (8.10)在这里z y x n w n v n u V 2222++=是推导的速度,Tz y x n n n n ],,[=表示的是壁面的单位法向矢量。

我们假设虚拟单元中的压力和密度与边界单元的值相等。

(即就是2p p w =)。

规则的双重中线的方案公式(8.1)中的边界条件需要更多的注意离散化的双重中线问题。

图8.4和图8.5分别表示了二维和三维的情况。

公式(8.2)中的对流通量可以在壁面上的每一个控制体积上的表面分开来计算。

这种方法与第一种构建单元顶点的方案相同。

对于一个四面体单元(就像图8.4中的1-3-4-5),压力可以通过公式8.8来代替。

对于六面体,棱形或者棱锥形,控制体积的表面是四边形(就像图8.5的面1-4-5-6),推导的公式可以如下来表示:)339(1615641int p p p p p +++=(8.11)如果边界元素是四面体(或者二维中的三角形),壁面的压力应该通过有限体积法来计算。

【7】在壁面的1-2部分,例如,*21-表面的压力可以通过如下的公式计算:)5(6121int p p p +=(8.12)对于四面体,例如图8.5中的壁面1-2-3,压力可以如下表示:)6(81321int p p p p ++=(8.13)8.2.2 粘性流动对于经过固体壁面的粘性流体,介于表面和流体的中垂直于表面的为零。

因此,我们我们可以称之为非滑动的边界条件。

在一个固定的壁面,速度分量在表面变为:0===w v u (8.14)对于非滑动的边界条件,有两个基本的结果。

第一点,我们不需要解壁面上的动量方程,这个在单元顶点方案中已经应用。

第二点,对流通量通过非滑动的壁面的公式已经在公式(8.2)中给出,并且在公式(2.24)中的项已经简化成T k ∇=Θ,因此,对流通量中的壁面压力同样可以通过非粘性流动来求得。

但是,虚拟的单元处理的方法不同。

以单元为中心的方案公式(8.14)中非滑动边界条件的处理可以通过利用虚拟单元来简化。

在一个绝热的壁面(没有热通量通过壁面),我们可以取(如图8.2)2121212121 w w v v u u E E -=-=-===,,,ρρ (8.15)同样的对于单元0和3。

方法同样适用于构建的以单元为中心的方案的和非构建的以单元为中心的方案(参考文献【6】)。

如果壁面的温度已知,速度分量仍然像方程(8.15)中一样是相反的。

利用壁面的温度可以将周围的温度线性的表示。

由于压力的梯度垂直于壁面并且为零,边界元素的压力同样可以用虚拟单元来表示(例如210 p p p ==)。

虚拟单元中的密度和总能量可以通过推导的值来计算。

单元顶点的方案由于动量方程不需要求解,壁面上的对流通量没有作用。

公式(2.23)中的粘性通量仅仅对能量方程中垂直于壁面的温度梯度有作用。

对于一个绝热的壁面,n T w ⋅∇为零。

因此,我们没有必要求解任何壁面上的对流或者粘性通量。

为了防止避免节点上非零速度分量的产生,动量方程中剩余的项都为零。

在已知壁面温度的情况下,我们可以直接设壁面(例如图8.3中的节点)2,(i )的总能量为(假设为理想气体)w i pi T c E 2,2,)(ργρ=(8.16)在这里,w T 表示给定的壁面温度。

剩余的动量和能量方程都设为零。

对于未建立一单元为中心的方案同样适用。

另外一种方法,对于一些应用来说比较简略,并且根本没有解决壁面的控制方程。

同样的,密度和能量可以直接简化为RT p w i i 3,2,=ρ 和γρ3,2,)(i i p E =(8.17)方程(8.17)中假设垂直于壁面没有压力梯度(因此32p p =)。

由于所有的守恒变量都是一定的,剩余的所有的方程都等于零。

这个方法同样适用于未构建的网格。

相关文档
最新文档