3多维非稳态对流扩散问题
热流问题的数值计算
整理,得
a PTP a E TE aW TW a N TN a S TS b
其中
rP r aE (x) e / e
rP r aW (x) w / w
0 P
rn x aN (y ) n / n
rs x aS (y ) s / s
( c) P V a t
0 0 0 a P a E aW a N a S a P S P V b S C V a P TP
aE , aW , aS , aN 都是相邻两节点间导热热阻的倒
0 a 数, 具有热惯性意义, P热惯性越大,上一
时层的温度对下一时层的影响越大。
4.3 源项及边界条件的处理
离散方程为:
a PTP a E TE aW TW a N TN a S TS b
其中
r aE re ( ) e / e
r aW rw ( ) w / w
0 P
rn aN (r ) n / n
rs aS (r ) s / s
( c) P V a V 0.5(rn rs )r t
经整理,得
dT dT e Ae ( ) e w Aw ( ) w ( S C S p TP ) AP(3) x 0 dx dx
令
Ae e aE (x) e
Aw w aW (x) w
b S C AP x 则式(3)变为:
a P a E aW S P AP x
因为离散方程都可表示为
a PTP anbTnb b
aP anb S P V
线性代数方程迭代求解收敛的一个充分必 要条件是对角占优,即:
扩散方程 稳态扩散与非稳态扩散.
一、扩散方程稳态扩散与非稳态扩散1.稳态扩散下的菲克第一定律(一定时间内,浓度不随时间变化dc/dt=0)单位时间内通过垂直于扩散方向的单位截面积的扩散物质流量(扩散通量)与该面积处的浓度梯度成正比即J=-D(dc/dx)其中D:扩散系数,cm2/s,J:扩散通量,g/cm2·s ,式中负号表明扩散通量的方向与浓度梯度方向相反。
可见,只要存在浓度梯度,就会引起原子的扩散。
x轴上两单位面积1和2,间距dx,面上原子浓度为C1、C2则平面1到平面2上原子数n1=C1dx ,平面2到平面1上原子数n2=C2dx若原子平均跳动频率f, dt时间内跳离平面1的原子数为n1f·dt跳离平面2的原子数为n2fdt,但沿一个方向只有1/2的几率,则单位时间内两者的差值即扩散原子净流量。
令,则上式2.扩散系数的测定:其中一种方法可通过碳在γ-Fe中的扩散来测定纯Fe的空心园筒,心部通渗碳气氛,外部为脱碳气氛,在一定温度下经过一定时间后,碳原子从内壁渗入,外壁渗出达到平衡,则为稳态扩散单位时单位面积中碳流量:A:圆筒总面积,r及L:园筒半径及长度,q:通过圆筒的碳量则:即:则:q可通过炉内脱碳气体的增碳求得,再通过剥层法测出不同r处的碳含量,作出C-lnr曲线可求得D。
第一定律可用来处理扩散中浓度不因时间变化的问3.菲克第二定律:解决溶质浓度随时间变化的情况,即dc/dt≠0两个相距dx垂直x轴的平面组成的微体积,J1、J2为进入、流出两平面间的扩散通量,扩散中浓度变化为,则单元体积中溶质积累速率为(Fick第一定律)(Fick第一定律)(即第二个面的扩散通量为第一个面注入的溶质与在这一段距离内溶质浓度变化引起的扩散通量之和)若D不随浓度变化,则故:4.Fick第二定律的解:很复杂,只给出两个较简单但常见问题的解a. 无限大物体中的扩散设:1)两根无限长A、B合?金棒,各截面浓度均匀,浓度C2>C12)两合金棒对焊,扩散方向为x方向3)合金棒无限长,棒的两端浓度不受扩散影响4)扩散系数D是与浓度无关的常数根据上述条件可写出初始条件及边界条件初始条件:t=0时, x>0则C=C1,x<0, C=C2边界条件:t≥0时, x=∞,C=C1, x=-∞, C=C2令,代入则,则菲克第二定律为即(1)令代入式(1)则有(2)若代入(2)左边化简有而积分有(3)令,式(3)为由高斯误差积分:应用初始条件t=0时x>0, c=c1,x<0, c=c2,从式(4)求得(5)则可求得(6)将(5)和(6)代入(4)有:上式即为扩散偶经过时间t扩散之后,溶质浓度沿x方向的分布公式,其中为高斯误差函数,可用表查出:根据不同条件,无限大物体中扩散有不同情况(1)B金属棒初始浓度,则(2)扩散偶焊接面处溶质浓度c0,根据x=0时,,则,若B棒初始浓度,则。
非稳态扩散名词解释
非稳态扩散名词解释非稳态扩散:非稳态扩散又称自发过程、非平衡扩散,是不需要外界能量的作用就能自发进行的一类扩散过程。
自然界中的扩散多数是非稳态扩散。
本文主要介绍非稳态扩散。
⑴非稳态扩散定义:在无限时间内,随着流体通过的断面面积变化而产生的浓度变化,称为非稳态扩散。
它是依靠单位时间通过的总面积或所有截面的总通量来描述的。
⑵非稳态扩散过程的分类(1)化学平衡理论解释的平衡分布过程(2)熵值理论解释的平衡分布过程(量子力学认为所有物质的微观运动形式都服从统计规律)。
⑶非稳态扩散机理研究意义:⑴使对象更加复杂;⑵导致各个部分的性质发生改变;⑶导致结构变得更加复杂;⑷使人们可以获取到更加丰富的信息。
非稳态扩散:非稳态扩散又称自发过程、非平衡扩散,是不需要外界能量的作用就能自发进行的一类扩散过程。
自然界中的扩散多数是非稳态扩散。
非稳态扩散分为两大类:吸附扩散和分子扩散。
吸附扩散是分子或颗粒物质因受其他物质吸引,相互接近而引起的一种扩散过程。
其特点是扩散的浓度比在空气中低,扩散的传质系数比在空气中大。
分子扩散是指由于温度差或化学反应等引起的扩散。
其特点是扩散的浓度比在空气中高,扩散的传质系数比在空气中小。
⑷扩散过程与状态变化特点:①属于等温、等压、等体积过程;②扩散速率不受浓度差的影响;③扩散的方向性;④存在固定的扩散系数;⑤有特定的传质系数。
⑸浓度梯度与传质系数:在某一瞬时,物料中每一点上的浓度梯度是该点处各个浓度单位的相应值的连乘积。
传质系数k是单位时间内从扩散体系一侧通过单位截面积物料的量,也就是单位时间内每单位面积上的物料浓度梯度除以物料的体积V。
⑹不稳态扩散与稳态扩散:稳态扩散与不稳态扩散之间的区别是前者有外力推动。
而后者则没有,这是非稳态扩散与稳态扩散最根本的区别。
⑺稳态扩散与不稳态扩散的比较:在稳态扩散中,由于断面面积变化而产生的浓度变化可忽略,所以其分析的方法与其他类型的扩散相同。
而在不稳态扩散中,外力推动必须满足浓度变化的情况下才有意义,故其分析方法不同。
对流扩散方程解析解
对流扩散方程解析解
《流体力学中的湍流扩散方程解析解》
一、什么是湍流扩散方程?
湍流扩散方程是描述物理流体扩散过程数学模型,是由流体力学中的湍流动力学概念推导出来的一种方程,是一种常用的偏微分方程。
它是一种描述在空间中湍流的扩散过程的数学方程,其目的是描述物质和能量在湍流中的传播。
二、湍流扩散方程的公式:
湍流扩散方程的公式为:
∂C/∂t = D∇2C
左侧的第一项是物质的局部变化率,t 代表时间;右侧的第一项用来描述物质在空间中的传播,D 为扩散系数,∇2C 为Laplace 算子。
三、湍流扩散方程的解析解:
1.快速波动方法:即快速 Fourier 过程,是一种快速处理湍流扩散方程的方法,其大致操作是用离散傅立叶变换把扩散方程转化为一个秩为 0
的傅立叶方程,然后使用傅立叶级数解决得出结果;
2.有限差分方法:给定的湍流扩散方程先采用有限的体积分解,即在时间及空间的二维平面上将扩散方程的计算区域划分成均匀的小单元,然后在每个区间内建立一个线性的有限差分矩阵,把扩散方程就变为简单的线性方程组;
3.格式方法:即 Finite Element 方法,用此方法可以把湍流扩散方程从不同的坐标方程中任意变换到球形坐标系,然后用有限元计算机程序解决;
4.积分方法:则是用数值积分的方法解决湍流扩散方程,包括 Runge-Kutta 方法、Adams 方法及其它积分的方法。
四、总结
湍流扩散方程是描述物理流体扩散过程的数学模型,是由流体力学中的湍流动力学概念推导出来的一种方程。
解决该方程有几种方法,即快速波动方法、有限差分方法、格式方法及积分方法。
以上是关于湍流扩散方程解析解的相关介绍,希望能够帮助到大家。
第四章扩散习题答案
第四章扩散习题答案4-1解:本题为非稳态扩散问题,采用菲克第二定律。
由题意知,初始条件和边界条件为:t=0时,x>0 C = C 0; t ≥0时,x=0 C = C s ;t ≥0时,x=∞ C = C 0所以,可应用误差函数解来求解。
c(x,t)= C s -(C s - C 0)erf (Dtx2),其中,Cs=0, C 0=0.85%, 取C=0.8% 则:(C s –C)/( C s –C 0)= erf(Dt x2)∴(0-0.8%)/( 0-0.85%)= erf(3600*10*1.1211-x)∴erf(3600*10*1.1211-x)=0.94 查表知,3600*10*1.1211-x=1.33∴x=0.00053m=0.53mm∴ 应车去 0.53mm 的深度。
4-2解:本题为非稳态扩散问题,采用菲克第二定律。
由题意知,初始条件和边界条件为:t=0时,x>0 C = C 0; t ≥0时,x=0 C = C s ;t ≥0时,x=∞ C = C 0所以,可应用误差函数解来求解。
c(x,t)= C s -(C s - C 0)erf (Dt x 2),其中,Cs=1.1%, C 0=0 ①β=Dt x2=10*10*5.825.02-=0.33∴ erf(β)=0.3593C= C s -( C s – C 0) erf(β)=1.1%-1.1%×0.5393=0.705%②β= Dt x2=10*10*5.820.12-=0.66∴ erf(β)=0.6494C= C s -( C s -C 0) erf(β)=1.1%-1.1%×0.6494=0.386%③β= Dt x2=10*10*5.822.12-=0.79∴ erf(β)=0.7361C= C s -( C s -C 0) erf(β)=1.1%-(1.1%-0)×0.7361=0.29%④β= Dt x2=10*10*5.825.12-=0.98∴ erf(β)=0.8342C= C s -( C s -C 0) erf(β)=1.1%-(1.1%-0)×0.8342=0.18%⑤β= Dt x2=10*10*5.8222-=1.31∴ erf(β)=0.9361C= C s -( C s -C 0) erf(β)=1.1%-(1.1%-0)×0.9361=0.07%渗层内碳的浓度分布曲线如下:4-3解:本题为非稳态扩散问题,采用菲克第二定律。
非稳态扩散名词解释
非稳态扩散名词解释非稳态扩散是指物理、材料、化学及生物系统中的扩散过程,它对应于热力学的熵的非稳态增加。
它的本质是物质或能量从一个地方向另一个地方转移,以满足不均衡态的需要和条件。
而且它还可以在某些情况下产生有益的结果,从而改善系统的性能和功能。
非稳态扩散涉及一系列不同的现象,如传热、传热和传质等。
扩散在热力学中是指物质穿过温度或其他物理变量时,将温度或其他物理变量不均匀分布于各处的过程。
热扩散可以在物体表面和内部之间传递热量,而扩散系数描述的是热传递的强度。
扩散过程的速率取决于物质的种类、浓度和热能等。
热扩散是通过物体温度从高温区至低温区的传热而实现的,它对于定义散热器的有效性至关重要。
传质扩散是指在不同区域之间物质的扩散过程,即使温度在整个系统中是均衡的,也可以发生物质的扩散。
传质扩散的驱动力是物质的溶质浓度差异。
当物质的浓度在某一位置大于其他位置时,物质就会从高浓度区域流向低浓度区域。
这种可能会导致不良的影响,如污染和气味的扩散,因此应该尽可能多地阻止传质扩散。
生物系统中的非稳态扩散是指物质或信息在生物体之间以及内部生物体内部的扩散。
在生物体之间,细胞间信号传导是一种重要的非稳态扩散过程,其涉及受体细胞、信号转导蛋白与接受信号的分子等,其中都包含扩散的本质。
生物体内部的非稳态扩散是指细胞的活性蛋白的扩散,涉及到细胞间质和细胞膜的传递,从而影响各种细胞活性及细胞内功能的分布。
从上述内容可以看出,非稳态扩散是一种普遍存在的物理现象,它可以发生在物理、材料、化学和生物系统中,对于系统的功能性和性能有着重要的影响。
非稳态扩散涉及不同的现象,包括热传导、传热和传质扩散,以及生物体内部的细胞间信号传导和细胞内蛋白的扩散。
因此,了解非稳态扩散的机制和其对系统性能影响的理解,对于提高材料或生物系统的效率以及可靠性具有重要的意义。
扩散方程稳态扩散与非稳态扩散
一、扩散方程稳态扩散与非稳态扩散1.稳态扩散下的菲克第一定律(一定时间,浓度不随时间变化dc/dt=0)单位时间通过垂直于扩散方向的单位截面积的扩散物质流量(扩散通量)与该面积处的浓度梯度成正比即J=-D(dc/dx)其中D:扩散系数,cm2/s,J:扩散通量,g/cm2·s ,式中负号表明扩散通量的方向与浓度梯度方向相反。
可见,只要存在浓度梯度,就会引起原子的扩散。
x轴上两单位面积1和2,间距dx,面上原子浓度为C1、C2则平面1到平面2上原子数n1=C1dx ,平面2到平面1上原子数n2=C2dx若原子平均跳动频率f, dt时间跳离平面1的原子数为n1f·dt跳离平面2的原子数为n2fdt,但沿一个方向只有1/2的几率,则单位时间两者的差值即扩散原子净流量。
令,则上式2.扩散系数的测定:其中一种方法可通过碳在γ-Fe中的扩散来测定纯Fe的空心园筒,心部通渗碳气氛,外部为脱碳气氛,在一定温度下经过一定时间后,碳原子从壁渗入,外壁渗出达到平衡,则为稳态扩散单位时单位面积中碳流量:A:圆筒总面积,r及L:园筒半径及长度,q:通过圆筒的碳量则:即:则:q可通过炉脱碳气体的增碳求得,再通过剥层法测出不同r处的碳含量,作出C-lnr曲线可求得D。
第一定律可用来处理扩散中浓度不因时间变化的问3.菲克第二定律:解决溶质浓度随时间变化的情况,即dc/dt≠0两个相距dx垂直x轴的平面组成的微体积,J1、J2为进入、流出两平面间的扩散通量,扩散中浓度变化为,则单元体积中溶质积累速率为(Fick第一定律)(Fick第一定律)(即第二个面的扩散通量为第一个面注入的溶质与在这一段距离溶质浓度变化引起的扩散通量之和)若D不随浓度变化,则故:4.Fick第二定律的解:很复杂,只给出两个较简单但常见问题的解a. 无限大物体中的扩散设:1)两根无限长A、B合?金棒,各截面浓度均匀,浓度C2>C12)两合金棒对焊,扩散方向为x方向3)合金棒无限长,棒的两端浓度不受扩散影响4)扩散系数D是与浓度无关的常数根据上述条件可写出初始条件及边界条件初始条件:t=0时, x>0则C=C1,x<0, C=C2边界条件:t≥0时, x=∞,C=C1, x=-∞, C=C2令,代入则,则菲克第二定律为即(1)令代入式(1)则有(2)若代入(2)左边化简有而积分有(3)令,式(3)为由高斯误差积分:应用初始条件t=0时x>0, c=c1,x<0, c=c2,从式(4)求得(5)则可求得(6)将(5)和(6)代入(4)有:上式即为扩散偶经过时间t扩散之后,溶质浓度沿x方向的分布公式,其中为高斯误差函数,可用表查出:根据不同条件,无限大物体中扩散有不同情况(1)B金属棒初始浓度,则(2)扩散偶焊接面处溶质浓度c0,根据x=0时,,则,若B棒初始浓度,则。
第五章对流扩散问题(假扩散)
该问题的数值解如下:
1 n i
u t n u t n (1 ) i i 1 x x
MUD : du d ux d (( ) ) dx dx 2 dx
第五章 对流扩散问题———假扩散
由图可以看出,在区间 P 2 ,中心差分格式预报的 P 值优于迎风格式的预报值。对比这两种格式,其扩散项 的处理是完全相同的,所不同的仅仅是对流项的处理上 ,在中心差分格式中对流项的差分格式具有二阶精度, 而在迎风格式中对流项的差分格式只具有一阶精度。在 区间 P 2 ,两种格式预报 P值所表现出的差异性恰恰是 这两种格式精度不同的体现。观察上图,迎风格式所预 报的 P值具有该高不高和该低不低的特点,这一特点正 是由一阶精度迎风格式所引起的扩散系数为 ux / 2 的 假扩散项造成的。也反映了假扩散项的影响。
n n ux ut 2 n ( ) i u( ) i (1 )( 2 ) i O( x 2 , t 2 ) t x 2 x x
由此可以看出,我们前边得到的差分方程所逼近的是 一个非稳态对流扩散问题,而非原型问题所要求的非 稳态对流问题。
第五章 对流扩散问题———假扩散
1 n (1 P
ut n ut n ) P W x x
用编号法表示
1 n i
1 n , n 在点 (i, n) i i 1
u t n u t n (1 ) i i 1 x x
做Taylor展开
n n u 2 n 1 2 n ( ) i u( ) i ( 2 ) i x ( 2 ) i t O( x 2 , t 2 ) t x 2 x 2! t
第五章 对流扩散问题———假扩散
将 i 1 和 i 1 台劳 级数展开代入
对流扩散与相间传质讲解
2019/6/8
对流扩散与相间传质
2/27
JA
( D
DE
)
dcA dz
对流传质通量
在湍流主体中
DE D 在层流内层中
NA
DE RT
dpA dz
NA
D RT
dpA dz
D 较小
dpA dz
较大,DE
0
在缓冲层内
NA
D DE RT
dpA dz
溶
质p
A
在 气
pi
相
中
的
分
压
2019/6/8
气相主体
气膜 液膜
液
膜
边
液相主体
界
气
膜
边
界
相界面 传质方向
图32 2 双膜模型
对流扩散与相间传质
溶
质
A 在
液
相
ci
中 的
c
摩 尔
浓
度
6/27
双膜模型的理论要点是:
①在气-液两相接触面附近,分别存在着呈滞流流动的稳定气膜和液膜。 溶质连续稳定地通过两膜,膜的厚度随流体流动状态而变化; ②气-液两相在相界面上呈平衡状态,即相界面上不存在传质阻力。如以低浓度 气体溶解为例,则平衡关系服从Henry定律,即有 ci Hpi 或 c Hp ,其中H 为溶解度系数,单位随 c 和 p 的单位而定; ③膜层以外的气、液相主体,由于流体的充分湍动,分压或浓度均匀化,无分压 或浓度梯度。
溶质由气相主体 湍流扩散气膜边界 p 分子扩散相界面气侧 pi 无阻力溶解 相界面液侧 ci 分子扩散液膜边界 c 湍流扩散液相主体
对流扩散方程解析解
对流扩散方程解析解对流扩散方程(Convection-DiffusionEquation)是在求解流体,如气体或液体的输运问题时需要使用的普通微分方程。
它表示物质被三种因素作用所引起的质量流动:对流、扩散和反应。
在本文中,我们将讨论对流扩散方程的解析解,以及它在工程中的重要作用。
首先,要理解对流扩散方程,我们必须从它的数学形式开始。
它可以用以下形式表示:$$frac{partial c}{partial t}+ vec{u} cdotabla c-Dabla^2 c=R$$在这里,$c$表示物质的浓度,$vec{u}$表示流体的速度,$D$表示物质的扩散系数,$R$表示反应的密度。
对流扩散方程的解析解是一种运用数学方法来求解这个方程的方法。
它主要是利用积分变换法(Integral Transform Method),将复杂的运动学问题转化为一组常微分方程求解。
解析解方法在解决一定类型的常微分方程时尤其有用,特别是当一个系统的边界条件是确定的时。
解析解的优势在于它可以提供直观的解,方便比较和评估结果,便于理解物理机理。
它也可以提供准确的结果,并可以用于组合的求解方法中。
在工程领域,对流扩散方程解析解的应用非常重要。
它可以被应用于温度或物质浓度输运,以及其他类似现象的计算。
例如,对流扩散方程可以用来模拟一定范围内扩散方式的热量传输,从而推测温度场分布;也可以用来模拟入口流场和出口的物质浓度的变化;它还可以用来描述各种物质在工程系统内的扩散问题。
再者,解析解方法也被广泛应用于制药行业。
对流扩散方程可以用来模拟药物在体内的运动,从而计算出最佳控制方案,以达到药物最佳疗效。
这不仅可以为药物分布模型提供依据,还可以用来估算药物组分以及药物与体细胞的相互作用等工程相关问题,从而帮助制药公司最大程度地提高药品安全性和疗效。
最后,对流扩散方程的解析解是一种非常有效的数学方法,它可以帮助我们更加清晰地理解流体输运问题,并可以提供准确可靠的结果。
第六章对流与扩散
该格式计算量比指数小,且与指数格式的解差别很小。
§ 6-3
通用表达式
为了在讨论中引入 PE J* J u d * 记 J x x ( ) d( ) d i i+1 P d x i+1/2 d( ) 1 界面i+ 上的值可以有界面两侧节点值表示
第六章
对流扩散方程的差分格式
导热型方程:(原始或经过变换的)
二阶导数项(扩散),源项
对流扩散方程:(动量或能量)
二阶导数项(扩散),源项 一阶导数项(对流),压力梯度。 一维稳态无内热源的对流扩散方程:
d d d ( u ) ( ) 密度, 扩散系数。 dx dx dx
对流热能量方程
aE Pe De
aE 1 1 Pe De 2
指数
aE 0 De
二.混合格式
虽然指数格式是精确解,但计算过繁,通过对 随 Pe 变化及其三条切线 aE Pe 0 De aE Pe Pe De aE 1 Pe 0 1 Pe De 2 斯帕尔丁提出 aPP aEE aww
F u J * 而 P ,J D ( ) D x * 根据通量守恒 Je J De Je D J* 0
P{De B(P e ) D A(P )} De A(Pe )E e D BP W
aE De A(Pe ) De{A(| Pe |) [| Pe ,0 |]}
Pe 10
aE Pe DE aE (1 0.1Pe )5 Pe DE aE (1 0.1Pe )5 DE aE 0 DE
10 Pe 0
0 Pe 10
Pe 10
(f)
三维对流扩散问题的有限元计算论文
连续性方程和运动方程)求出。另外一种是 u 和流场中的密度 或速度 v 有直接的
0 kij
u x j
ni
1u
2
0
式中0 、1 、2 是已知函数。 初始条件为:在 t = 0 时,给出 u 的分布
u(xi , 0) u (xi )
方程(1.1)式中的速度 v 一般有两种情况:一种是 u 对于流场中的密度 或速度 v 没
有直接影响,或者影响很小可以忽略,此时 v 可以独立通过流体力学的主流方程(即
V (t t) vn
S V (t)
u u
图 1.1 对流过程 Figure 1.1 Process of convection
d dt
udV
V
V
u t
dV
Ñ uvndS
s
其中 vn v n 是流体速度在 S 面上的法向分量。利用 Green-Gauss 公式
Ñ uvndSቤተ መጻሕፍቲ ባይዱ div(uv)dV
Keywords: Convection-diffusion equation, Finite element method, Object-oriented method, Programming
II
重庆大学硕士学位论文
1 绪论
1绪 论
对流扩散方程是一类基本的运动方程,它可以用来描述河流污染、大气污染、 核废物污染中污染物质的分布,流体的流动和流体中热的传导等众多物理现象。 对于同时伴有物质输运和分子扩散的物理过程以及粘性流体的流动,其数学模型 通常为对流扩散方程或含有此类方程的偏微分方程组的定解问题。对流扩散问题 数值计算方法的研究具有重要的理论和实际意义,可用于环境科学、能源开发、 流体力学和电子科学等许多领域。对流扩散问题的有效数值解法一直是计算数学 中重要的研究内容,求解对流扩散方程的数值方法主要是有限差分法(FDM)、有限 元法(FEM)、有限体积法(FVM)、有限解析法(FAM)、边界元法(BEM)、谱方法(SM) 等多种方法。但是对于对流占优问题,用通常的差分法或有限元法进行求解将出 现数值震荡。为了克服数值震荡,80 年代,J.Douglas,Jr.和 T.F.Russell 等提出特征 修正技术求解对流扩散占优的对流扩散问题,与其它方法相结合,提出了特征有 限元方法、特征有限差分方法、特征混合元方法;T.J.Hughes 和 A.Brooks 提出过 一种沿流线方向附加人工黏性的间断有限元法,称为流线扩散方法(SDM)。有限差 分法、有限元法、有限体积法是工程应用中的主要方法。本文所研究是低 Peclet 数下,粘性不可压缩流体的流动中, 三维定常的对流扩散问题的有限元解法。验 证 Peclet 数较小时,三维定常的对流扩散方程的有限元解法的有效性,用面向对 象的 C++程序设计语言开发出有应用价值的通用程序,对相关的对流扩散问题进 行数值模拟。对于对流占优对流扩散问题的求解,采用本文方法,要得到稳定解, 则要通过加密有限元网格来实现。
第五章对流扩散问题(多维对流扩散问题)
n 1 P
U P x VP y ) S P
xx yy 2 A x 2 B y g
A 0.5 RU P , B 0.5 RVP , R n 1 g (P P ) SP t
第三章 对流扩散问题———多维对流扩散问题
给线化的方程配置边界条件
R n 1 ( C nb nb C P P C P S P ) t n 1 P ( 1 C P R / t )
8
第三章 对流扩散问题———多维对流扩散问题 如果把上式中的h和k改为x和y,则就获得了以前边 所配一类边界条件下对流扩散问题的精确解
R n 1 ( Cnb ( x, y ) nb CP ( x, y ) P C P ( x, y )S P ) t n 1 ( x , y ) (1 C P ( x, y )R / t )
De Dw 1, W 50, Fe Fw 4 Ds Dn 1, Fs Fn 4 E 200, S 50, N 200,
则,
a E 1 4 / 2 1 a W 1 4 / 2 3 a N 1 4 / 2 1 a S 1 4 / 2 3 a P 1 3 1 3 4
分离变量法
R n 1 ( C nb nb C P P C P S P ) t n 1 P ( 1 C P R / t )
8
第三章 对流扩散问题———多维对流扩散问题 其中,
CEC exp( Ah ) exp(Ah ) PB , C WC PB 2 cosh( Ah ) 2 cosh( Ah ) exp(Ah Bk ) (1 PA PB ) 4 cosh( Ah ) cosh( Bk ) exp(Ah Bk ) (1 PA PB ) 4 cosh( Ah ) cosh( Bk )
流体的稳态与非稳态
流体的稳态与非稳态流体是一种特殊的物质形态,它具有流动性和变形性。
在许多物理和工程学领域,研究流体的稳态与非稳态是非常重要的。
稳态指的是流体在一定条件下保持恒定的流动状态,不随时间变化;非稳态则指流体的流动状态会随时间发生变化。
一、稳态流动稳态流动是指流体在空间和时间上保持恒定的流动状态。
稳态条件下,流体的速度场、压力场和密度场均不随时间变化,且在空间各点上保持恒定。
稳态流动通常发生在稳定的管道、河流以及一些工程设备中。
稳态流动具有以下特点:1. 守恒性质:在稳态条件下,流体在整个流动过程中,质量、动量和能量都保持不变。
这意味着在稳态流动中,流体在任何两点之间的质量流量、动量和能量传递都是平衡的。
2. 流速分布均匀:稳态流动时,流速在整个流场内分布均匀,没有明显的变化。
这是因为在稳态流动中,速度场保持不变,流体沿着固定路径流动,流速沿流动方向不发生变化。
3. 马克氏数恒定:在稳态流动中,流体的马克氏数为常数。
马克氏数是流体力学中一个重要的无量纲参数,它描述了流体的运动状态和热力学性质之间的关系。
二、非稳态流动非稳态流动是指流体在空间和时间上不断变化的流动状态。
非稳态流动常常发生在开关阀门、开始或停止流动等瞬态过程中,以及在激波、涡旋等复杂流动现象中。
非稳态流动具有以下特点:1. 时间变化:非稳态流动涉及流体的时间变化,即流速、压力和密度等在一段时间内发生变化。
在非稳态流动中,流场的速度和压力场随时间的推移而改变。
2. 波动与脉动:非稳态流动时,流体中的流速存在波动和脉动现象。
波动是指流速的周期性变化,脉动则是指流速的非周期性变化。
这些波动和脉动可能会导致流动的非均匀性和紊乱性。
3. 动量传播:在非稳态流动中,流体中的质量、动量和能量会随着流动的变化而传播,并在流场中发生扩散、对流和传导过程。
三、稳态与非稳态的区别稳态流动与非稳态流动之间存在明显的区别:1. 流动状态:稳态流动是持续且恒定的,而非稳态流动是变化的和不稳定的。
稳态扩散和非稳态扩散名词解释
稳态扩散和非稳态扩散是物理学和化学领域中常用的术语,用来描述物质在空间中的扩散行为。
稳态扩散是指物质在均匀介质中的扩散过程,而非稳态扩散则是指物质在非均匀介质中的扩散过程。
本文将通过对稳态扩散和非稳态扩散的名词解释,解析其物理意义、数学表达和实际应用,帮助读者更好地理解这两个概念。
一、稳态扩散的名词解释1.1 稳态扩散的物理意义稳态扩散是指当物质在均匀介质中的浓度分布达到稳定状态时的扩散过程。
在稳态扩散中,物质的浓度分布不再发生变化,达到了动态平衡状态。
这种扩散过程通常以弥散系数来描述,可以用弥散方程进行数学建模,是一种重要的物质传输过程。
1.2 弥散系数的数学表达在稳态扩散中,物质的弥散系数是一个重要的物理参数,用来描述物质在均匀介质中扩散的能力。
它通常用符号D表示,是一个与时间和空间无关的常数。
弥散系数与介质的性质、温度和压力等因素有关,不同的物质在不同的介质中有不同的弥散系数。
1.3 稳态扩散的实际应用稳态扩散在化学工程、环境科学和生物医学等领域有着广泛的应用。
在化学工程中,稳态扩散常用来描述气体或液体在反应器中的传质过程;在环境科学中,稳态扩散被用来研究大气、水体和土壤中的污染物传播行为;在生物医学中,稳态扩散可用于分子在细胞内的扩散和运输研究。
二、非稳态扩散的名词解释2.1 非稳态扩散的物理意义非稳态扩散是指当物质在非均匀介质中的浓度分布随时间和空间发生变化的扩散过程。
在非稳态扩散中,物质的浓度分布不断改变,未达到动态平衡状态,通常需要考虑时间和空间的变化。
2.2 非稳态扩散的数学表达在非稳态扩散中,物质的扩散过程通常需要考虑时间和空间的变化,所以需要使用偏微分方程进行描述。
这类方程常常包括时间导数和空间导数,需要通过适当的数值或解析方法进行求解,得到物质浓度随时间和空间的变化规律。
2.3 非稳态扩散的实际应用非稳态扩散在材料科学、地球科学和生物学等领域有着重要的应用价值。
在材料科学中,非稳态扩散常用来研究材料中的晶体生长和变形过程;在地球科学中,非稳态扩散被用来描述地表和地下水体中的渗透和溶质迁移过程;在生物学中,非稳态扩散可用于描述细胞内物质的输运和信号传导过程。
5.3多维非稳态对流扩散问题
考虑非常数源项时的数值结果
评述
• 这种虚假扩散现象是一种特殊情况,但可 能在不同的离散格式、不同源项分布情况 下出现。
• 如何减少这种虚假扩散,还有待深入研究 ,但对流项采用高阶精度离散格式,对减 轻相应的影响显然是有益的
虚假扩散无处不在, 无孔不入!!
采用高阶的对流项差分格式 总是有一定好处的
c(unj
un j1
)
(c at x)
修正的偏微分方程(MPDE)
• 迎风格式的泰勒展开
ut
aux
t 2
utt
ax 2
uxx
t 2 6
uttt
ax2 6
uxxx
• MPDE
自循环消元过程
ut
aux
ax 2
(1
c)uxx
ax2 6
(2c
1)(c
1)uxxx
虚假的流向扩散
• MPDE中的二阶空间导数代表扩散作用(粘 性效应),相当在原始方程中增加了扩散作 用(人工粘性作用),这引入了原始方程中没 有的一种虚假扩散。
5.4 对流扩散方程离散格式的 虚假扩散问题
1. 人工粘性引起流向扩散 2. 网格取向引起交叉扩散 3. 非常数源项带来的虚假扩散
5.4.1人工粘性所引起的 流向扩散
修正的偏微分方程(MPDE)
• 一维对流方程(波动方程)
u a u 0 (a 0) t x
• 一阶迎风格式
un1 j
unj
• 根据采用的三点离散格式不同,选定A(|P|) 函数形式不同,参见前一节的表格
5.3.2 三维非稳态对流扩散方程
• 离散结果 aPP aEE aWW aNN aSS aTT aBB b
非稳态扩散名词解释
非稳态扩散名词解释
非稳态扩散:一种“活跃”的化学扩散
在微观层面上,物质扩散体现为分子、原子或颗粒在不断变化的温度和压力的
作用下,穿过物质中的液体、气体,或者在它们的表面上的不断运动。
扩散,也可以在物理或化学上被描述为一种“运动”,而其中物质在运动过程中发生变化。
综上所述,扩散是宏观层面上和微观层面上不可缺少的一种运动过程。
在物理学中,扩散分为稳态扩散和非稳态扩散两类。
非稳态扩散是一种“活跃”的化学扩散,它会对它所经过部位的物态和选择性
形成一定的影响。
非稳态扩散过程中所演变出的物质的特性和空间散布情况反映出来,它的动力学特性一般由粒子的时空变化表征。
非稳态扩散一般不受外力的影响,它比一般的扩散更加活跃,所以在物质的多种化学活动中,它都起着重要的作用。
此外,非稳态扩散特别常见于微小体系,其会受到表面能和粒子相互作用的影响。
非稳态扩散比起稳态来说有较大的自主性,它能够更有效地传播激发信号,在传输消息中起到重要的作用。
它可以被用来驱动多时间尺度短、热休克信号跳跃或复杂化学反应等现象,发生在不稳定的体系中。
总而言之,非稳态扩散是一种多变现象,先于某一物质的变化而发生,而它的
演变轨迹由粒子的时空变化所誊描,它的影响力充满多样性,因为它可以驱动过多个时间尺度复杂的化学变化现象。
未稳态扩散可以说是化学活动中要有的不可缺
少的部分,它所提供的经历可以帮助研究者更深入地分析和预测化学系统的运行及其复杂程度。
扩散方程稳态扩散与非稳态扩散
扩散⽅程稳态扩散与⾮稳态扩散⼀、扩散⽅程稳态扩散与⾮稳态扩散1.稳态扩散下的菲克第⼀定律(⼀定时间内,浓度不随时间变化dc/dt=0)单位时间内通过垂直于扩散⽅向的单位截⾯积的扩散物质流量(扩散通量)与该⾯积处的浓度梯度成正⽐即J=-D(dc/dx)其中D:扩散系数,cm2/s,J:扩散通量,g/cm2·s ,式中负号表明扩散通量的⽅向与浓度梯度⽅向相反。
可见,只要存在浓度梯度,就会引起原⼦的扩散。
x轴上两单位⾯积1和2,间距dx,⾯上原⼦浓度为C1、C2则平⾯1到平⾯2上原⼦数n1=C1dx ,平⾯2到平⾯1上原⼦数n2=C2dx若原⼦平均跳动频率f, dt时间内跳离平⾯1的原⼦数为n1f·dt跳离平⾯2的原⼦数为n2fdt,但沿⼀个⽅向只有1/2的⼏率,则单位时间内两者的差值即扩散原⼦净流量。
令,则上式2.扩散系数的测定:其中⼀种⽅法可通过碳在γ-Fe中的扩散来测定纯Fe的空⼼园筒,⼼部通渗碳⽓氛,外部为脱碳⽓氛,在⼀定温度下经过⼀定时间后,碳原⼦从内壁渗⼊,外壁渗出达到平衡,则为稳态扩散单位时单位⾯积中碳流量:A:圆筒总⾯积,r及L:园筒半径及长度,q:通过圆筒的碳量则:即:则:q可通过炉内脱碳⽓体的增碳求得,再通过剥层法测出不同r处的碳含量,作出C-lnr曲线可求得D。
第⼀定律可⽤来处理扩散中浓度不因时间变化的问3.菲克第⼆定律:解决溶质浓度随时间变化的情况,即dc/dt≠0两个相距dx垂直x轴的平⾯组成的微体积,J1、J2为进⼊、流出两平⾯间的扩散通量,扩散中浓度变化为,则单元体积中溶质积累速率为(Fick第⼀定律)(Fick第⼀定律)(即第⼆个⾯的扩散通量为第⼀个⾯注⼊的溶质与在这⼀段距离内溶质浓度变化引起的扩散通量之和)若D不随浓度变化,则故:4.Fick第⼆定律的解:很复杂,只给出两个较简单但常见问题的解a. ⽆限⼤物体中的扩散设:1)两根⽆限长A、B合?⾦棒,各截⾯浓度均匀,浓度C2>C12)两合⾦棒对焊,扩散⽅向为x⽅向3)合⾦棒⽆限长,棒的两端浓度不受扩散影响4)扩散系数D是与浓度⽆关的常数根据上述条件可写出初始条件及边界条件初始条件:t=0时, x>0则C=C1,x<0, C=C2边界条件:t≥0时, x=∞,C=C1, x=-∞, C=C2令,代⼊则,则菲克第⼆定律为即(1)令代⼊式(1)则有(2)若代⼊(2)左边化简有⽽积分有(3)令,式(3)为由⾼斯误差积分:应⽤初始条件t=0时x>0, c=c1,x<0, c=c2,从式(4)求得(5)则可求得(6)将(5)和(6)代⼊(4)有:上式即为扩散偶经过时间t扩散之后,溶质浓度沿x⽅向的分布公式,其中为⾼斯误差函数,可⽤表查出:根据不同条件,⽆限⼤物体中扩散有不同情况(1)B⾦属棒初始浓度,则(2)扩散偶焊接⾯处溶质浓度c0,根据x=0时,,则,若B棒初始浓度,则。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
• 对一阶精度迎风格式加入适量的逆耗散,以减小 扩散系数 • 离散格式中包含更多邻节点个数
5.4.3 非常数源项引起的 虚假扩散
对称边界
• 由对称性,有
u v 0, 0, 0 y y
固壁边界
• 对粘性流体,壁面无渗透,其壁面速度为 零,即 u v 0 • 对于 ,可提1、2、3类边界条件。
出口边界
• 难点:除非实测,不可能获得出口截面信息
• 出口截面局部坐标单向化:假定出口截面节 点对它近邻的内节点无影响,从而令边界节 点对内节点的影响系数为零
aN Dn A( Pn ) Dn A Pn F , 0 n
aS Ds B( Ps ) Ds A Ps F , 0 s
• 根据采用的三点离散格式不同,选定A(|Βιβλιοθήκη |) 函数形式不同,参见前一节的表格
5.4.2 网格取向效应引起的 交叉扩散
由于网格线和流线之间并非平行或垂直 ,而是有一定角度的交叉而导致的扩散
虚假扩散逐渐抹平阶梯分布
0
0
来流与网格线平行和交叉时 迎风格式计算结果
P W
P W S 2
对 策
• 尽量减小流线与网格线间的倾斜和交叉。采用自 适应网格,如“旋转坐标”技术。 • 改进对流项格式设计方案,采用高阶精度迎风格式
• 连续方程积分结果
0 P P
t
xy ( Fe Fw ) ( Fn Fs ) 0
两式相减合并
0 0 P (P P )
t ( SC SPP )xy
xy ( J e FeP ) ( J w FwP ) ( J n FnP ) ( J s FsP )
• 需要注意:一定要跟连续方程联立,才能 得到正确的结果,才能适用于可压和不可 压的情况
最终表达式
aPP aEE aWW aNN aSS b
aE De A( Pe ) De A Pe F , 0 e aW Dw B( Pw ) Dw A Pw F , 0 w
J x t s w x dxdydt J x e J x w yt ( J e J w )t t t n e J y J y n J y s xt ( J n J s )t t s w y dxdydt
5.3.2 三维非稳态对流扩散方程
• 离散结果
aPP aEE aWW aNN aSS aT T aBB b
• 系数表达式见课本
5.3.3 多维对流扩散问题的 边界条件处理
几种可能的边界条件
• 以有回流的突扩通道为例
入口边界
• 一般规定入口边界上的函数值 和 流速 u 和 v 的分布
5.4 对流扩散方程离散格式的 虚假扩散问题
1. 人工粘性引起流向扩散 2. 网格取向引起交叉扩散 3. 非常数源项带来的虚假扩散
5.4.1人工粘性所引起的 流向扩散
修正的偏微分方程(MPDE)
• 一维对流方程(波动方程)
u u a 0 (a 0) t x
• 一阶迎风格式
1 n n n un u c ( u u j j j j 1 )
(c at x)
修正的偏微分方程(MPDE)
• 迎风格式的泰勒展开
t ax t 2 ax 2 ut aux utt uxx uttt uxxx 2 2 6 6
自循环消元过程
• MPDE
ax ax 2 ut aux (1 c)uxx (2c 1)(c 1)uxxx 2 6
• 连续方程
( u) ( v ) 0 t x y
引入通量密度
• 对流扩散总通量密度:
J x u , x J y v y
• 质量通量密度:
Fx u, Fy v
用通量表示的控制方程
• 控制方程:
( ) J x J y S t x y
• 连续方程
( ) Fx Fy 0 t x y
2.控制容积积分离散
• 非稳态项:假设沿空间为均匀分布
( ) 0 dtdxdy ( ) ( ) xy P P t
s
n
e
t t
w t
• 对流、扩散通量项:时间积分取隐式,空 间取均匀分布
5.3 多维非稳态对流扩散问题
5.3.1 二维非稳态对流扩散 方程的离散
1.直角坐标系下的对流扩散 方程和连续方程
• 控制方程
( ) ( u ) ( v ) S t x y x x y y
虚假的流向扩散
• MPDE中的二阶空间导数代表扩散作用(粘 性效应),相当在原始方程中增加了扩散作 用(人工粘性作用),这引入了原始方程中没 有的一种虚假扩散。 • 流向扩散(streamwise diffusion): 只要求解 函数顺流向存在不为零的一阶导数时,它 使方程的真解被光滑,导致数值计算误差
t t n e
源项
• 线化为
S SC SP ( SP 0)
• 时间、空间均取均匀分布
s
n
e
t t
w t
Sdtdxdy ( SC SP )xyt
积分结果
( ) P ( )0 P xy ( J e J w ) ( J n J s ) ( SC SPP )xy t