比拟理论与无量纲化导出特征数

$\begin{cases} \frac{\partial U}{\partial X}+\frac{\partial V}{\partial Y}=0\\ U\frac{\partial U}{\partial X}+V\frac{\partial U}{\partial Y}=\frac{1}{Re}\frac{\partial^2 U}{\partial Y^2}\\ U\frac{\partial \Theta}{\partial X}+V\frac{\partial \Theta}{\partial Y}=\frac{1}{Re}\frac{\partial ^2 \Theta}{\partial Y^2}\\ \end{cases}\tag{Pr=1}$
$\frac{1}{X}Nu_X=\frac{\partial \Theta}{\partial Y}|_{y=0}$
$\frac{\partial U}{\partial Y}|_{y=0}=\frac{C_f}{2}Re_x$

忽略体积力，压力的微分方程组如下：

$\begin{cases} \frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0\\ u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\nu(\frac{\partial^2 u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2})\\ u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}=\nu(\frac{\partial^2 v}{\partial x^2}+\frac{\partial^2 v}{\partial y^2})\\ u\frac{\partial t}{\partial x}+v\frac{\partial t}{\partial y}=a(\frac{\partial ^2t}{\partial x^2}+\frac{\partial ^2t}{\partial y^2})\\ h_x=-\frac{\lambda}{t_w-t_f}(\frac{\partial t}{\partial y})_{y=0} \end{cases}\tag{微分方程组}$

1.无量纲化

$X=\frac{x}{l};Y=\frac{y}{l};P=\frac{p}{\rho u^2_\infty};U=\frac{u}{u_\infty};V=\frac{v}{u_\infty};\Theta=\frac{t-t_w}{t_f-t_w}\tag{无量纲化}$

①连续性方程无量纲化

$\frac{\partial U}{\partial x}= \frac{\partial \frac{u}{u_\infty}}{\partial x}=\frac{1}{u_\infty}\frac{\partial u}{\partial x}\\ 连续性方程两边除以u_\infty\\ \frac{\partial U}{\partial x}+\frac{\partial V}{\partial x}=0\\ \because X=x/l; l\cdot dX=dx;\\ \therefore \frac{\partial U}{\partial x}=\frac{\partial U}{\partial X\cdot l};\\ l\neq 0等式两边除以\\ \frac{1}{l}\frac{\partial U}{\partial x}+\frac{1}{l}\frac{\partial V}{\partial x}=0\\$
$\color{#F00}{\frac{\partial U}{\partial X}+\frac{\partial V}{\partial Y}=0\tag{1}}$

②x方向上动量微分方程无量纲化

$动量微分方程两边除以u_\infty^2\\ \frac{u}{u_\infty}\frac{\frac{\partial u}{u_\infty}}{\partial x}+\frac{v}{u_\infty}\frac{\frac{\partial u}{u_\infty}}{\partial y}=\frac{\nu}{u_\infty}(\frac{\partial ^2(\frac{u}{u_\infty})}{\partial x^2}+\frac{\partial ^2(\frac{u}{u_\infty})}{\partial x^2})\\ 变化为：U\frac{\partial U}{\partial x}+V\frac{\partial U}{\partial y}=\frac{\nu}{u_\infty}(\frac{\partial^2 U}{\partial x^2}+\frac{\partial^2 U}{\partial y^2})\\ 注意一下，偏微分中\frac{\partial^2 u}{\partial x^2}=\frac{\partial }{\partial x}(\frac{\partial u}{\partial x})，因此对u的偏导项除以一个u_\infty就可以\\ 而对于\partial x^2，相当于平方项需要除以l^2\\ \because X=x/l; x=X\cdot l;x^2=X^2\cdot l^2\\ 微分方程式变为：\\ U\frac{\partial U}{\partial X}\cdot \frac{1}{l}+V\frac{\partial U}{\partial Y}\cdot \frac{1}{l}=\frac{\nu}{u_\infty\cdot l^2}(\frac{\partial^2 U}{\partial X^2}+\frac{\partial^2 U}{\partial Y^2})\\ 上式两边乘以l\\$
$\color{#F00}{U\frac{\partial U}{\partial X}+V\frac{\partial U}{\partial Y}=\frac{\nu}{u_\infty\cdot l}(\frac{\partial^2 U}{\partial X^2}+\frac{\partial^2 U}{\partial Y^2})=\frac{1}{Re}(\frac{\partial^2 U}{\partial X^2}+\frac{\partial^2 U}{\partial Y^2})\tag{2}}$

在(2)式中， $\frac{\nu}{u_\infty\cdot l}=\frac{1}{Re}$

雷诺数通过对x方向上的动量方程无量纲化导出

（该方法是对已经学过的微分方程处理得到一个已知的雷诺数的范例，其他微分方程也是同样的套路）

③y方向上的动量方程同样可以得到

$\color{#F00}{U\frac{\partial V}{\partial X}+V\frac{\partial V}{\partial Y}=\frac{\nu}{u_\infty\cdot l}(\frac{\partial^2 V}{\partial X^2}+\frac{\partial^2 V}{\partial Y^2})=\frac{1}{Re}(\frac{\partial^2 V}{\partial X^2}+\frac{\partial^2 V}{\partial Y^2})\tag{3}}$

④能量微分方程的无量纲化

同样的手段，经过三步分别对u,x,t无量纲化
$U\frac{\partial \Theta}{\partial X}+V\frac{\partial \Theta}{\partial Y}=\frac{a}{u_\infty l}(\frac{\partial \Theta^2}{\partial X^2}+\frac{\partial \Theta^2}{\partial Y^2})$

$\frac{a}{u_\infty l}是什么无量纲数？$

$已知\frac{\nu}{u_\infty l}=\frac{1}{Re}\\ \frac{a}{u_\infty l}=\frac{a \cdot \nu}{u_\infty l\cdot \nu}\\ \because Pr=\frac{\nu}{a}\\ \therefore \frac{a}{u_\infty l}=\frac{\nu \cdot a}{u_\infty l\cdot \nu}=\frac{1}{Re\cdot Pr}$
$\color{#F00}{U\frac{\partial \Theta}{\partial X}+V\frac{\partial \Theta}{\partial Y}=\frac{1}{Re\cdot Pr}(\frac{\partial \Theta^2}{\partial X^2}+\frac{\partial \Theta^2}{\partial Y^2})\tag{4}}$

二、数据处理

①由(1)式 $\frac{\partial U}{\partial X}+\frac{\partial V}{\partial Y}=0$ 的实验值进行数据处理：

$f_1(X,Y,U,V)=0\tag{5}$

②由(2)式 $U\frac{\partial U}{\partial X}+V\frac{\partial U}{\partial Y}=\frac{\nu}{u_\infty\cdot l}(\frac{\partial^2 U}{\partial X^2}+\frac{\partial^2 U}{\partial Y^2})=\frac{1}{Re}(\frac{\partial^2 U}{\partial X^2}+\frac{\partial^2 U}{\partial Y^2})$

U是这四个实验值的函数：

$U=f_2(X,Y,V,Re)\tag{6}$

③由(3)式 $U\frac{\partial V}{\partial X}+V\frac{\partial V}{\partial Y}=\frac{\nu}{u_\infty\cdot l}(\frac{\partial^2 V}{\partial X^2}+\frac{\partial^2 V}{\partial Y^2})=\frac{1}{Re}(\frac{\partial^2 V}{\partial X^2}+\frac{\partial^2 V}{\partial Y^2})$

V是这四个实验值的函数：

$V=f_3(X,Y,U,Re)\tag{7}$

④由(4)式 $U\frac{\partial \Theta}{\partial X}+V\frac{\partial \Theta}{\partial Y}=\frac{1}{Re\cdot Pr}(\frac{\partial \Theta^2}{\partial X^2}+\frac{\partial \Theta^2}{\partial Y^2})$

$\Theta$ 是六个实验值的函数
$\Theta=f_4(X,Y,U,V,Re,Pr)\tag{8}$

⑤一般V比U小，把V的表达式(6)带入(7)然后得到

$V=f_5(X,Y,Re)\tag{9}$

⑥(9)式带回(6)，消去一个因变量

$U=f_6(X,Y,Re)\tag{10}$

可以看到，无量纲的速度U和V就是坐标X，Y与雷诺数的函数

⑦(9)+(10)带回(8)

$\Theta=f_7(X,Y,Re,Pr)\tag{11}$

无量纲的温度分布最终需要得到对流换热系数，其是坐标X，Y与Re、Pr数的函数

⑧带回对流换热微分方程: $h_x=-\frac{\lambda}{t_w-t_f}(\frac{\partial t}{\partial y})_{y=0}$

$\Theta=\frac{t-t_w}{t_f-t_w}\\ \therefore h_x=-\frac{\lambda}{-(t_f-t_w)}(\frac{\partial(t-t_w)}{\partial(Y\cdot l)})_{y=0}=\frac{\lambda}{l}(\frac{\partial{(\frac{t-t_w}{t_f-t_w}})}{\partial Y})_{Y=0}=\frac{\lambda}{l}(\frac{\partial \Theta}{\partial Y})_{Y=0}$
$又Nu_x=h_x\cdot x/\lambda$

$Nu_x=\frac{x}{l}(\frac{\partial \Theta}{\partial Y})_{Y=0}=X(\frac{\partial \Theta}{\partial Y})_{Y=0}\tag{12}$

⑨局部努塞尔数

$Nu_x=f_8(X,Y,\Theta)|_{Y=0}$
因为Y=0确定，所以局部努塞尔数是X,Re,Pr的函数
$Nu_x=f_9(X,Re,Pr)$

⑩积分平均掉X，努塞尔数

$Nu=f_{10}(Re,Pr)$

最后编辑于：2022.11.30 10:16:24

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 221,576评论 6赞 515
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 94,515评论 3赞 399
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 168,017评论 0赞 360
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 59,626评论 1赞 296
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 68,625评论 6赞 397
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 52,255评论 1赞 308
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 40,825评论 3赞 421
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 39,729评论 0赞 276
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 46,271评论 1赞 320
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 38,363评论 3赞 340
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 40,498评论 1赞 352
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 36,183评论 5赞 350
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 41,867评论 3赞 333
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 32,338评论 0赞 24
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 33,458评论 1赞 272
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 48,906评论 3赞 376
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 45,507评论 2赞 359