高等流体力学-第五讲

1北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论本讲主要内容本讲主要内容本讲主要内容本讲主要内容一、基本概念二、描述扩散运动的基本运动方程三、扩散系数及其分析确定方法四、扩散运动的解析解五、岸边排放与中心排放污染带的计算2北京工业大学市政学科部——马长明高等流体(水）力学讲稿一、基本概念1111、扩散现象烟囱排烟；河流排污；水面蒸发；食糖与食盐的溶解等。2222、传输过程流体中所含有物质（如各种污染物，也包括动量、能量和热量）在流场中某一处到另一处转移的过程。3333、扩散（DiffusionDiffusionDiffusionDiffusion）是一类传输过程，指物质由含量高处向含量低处的传输过程。（1111）按扩散的机制可分为：分子扩散、对流传输扩散及紊动扩散。1111）分子扩散（MolecularDiffusionMolecularDiffusionMolecularDiffusionMolecularDiffusion）：由分子运动产生。2222）紊动扩散（TurbulentDiffusionTurbulentDiffusionTurbulentDiffusionTurbulentDiffusion）：由流体质点的紊动产生的扩散。3333）移流传输（AdvectionAdvectionAdvectionAdvection）：扩散物随同流体质点的时均运动而转移。第五讲扩散理论3北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论（2222）按研究问题的类型分：1111）剪切流中的离散（DispersionDispersionDispersionDispersion）：由于剪切流中速度分布不均匀产生含有随流散开的作用，也称弥散。离散中包含有移流扩散和紊动扩散。2222）射流扩散：指从各种排泄口喷出流入周围另一流体域内运动的一股流体。包括移流扩散和紊动掺混扩散。其中射入同种性质的流体内称淹没射流，射入不同性质的流体内为非淹没射流。按射流的原动力还可分为：动量射流、浮力羽流、浮射流。动量射流（JetsJetsJetsJets）：射流以出流的动量（MomentumMomentumMomentumMomentum）为原动力，该动量对射流运动起主要作用。浮力羽流（PlumesPlumesPlumesPlumes）：浮力（BouyancyForcesBouyancyForcesBouyancyForcesBouyancyForces）是原动力，产生的运动形态呈羽毛状。如烟气，水体中泄入污染液体后的运动等。浮射流（BuoyantJetsBuoyantJetsBuoyantJetsBuoyantJets）：原动力既有动量又有浮力。3333）分层流（StratifiedFlowingStratifiedFlowingStratifiedFlowingStratifiedFlowing）：在重力场中密度不均匀的流体形成有层次的流动。4北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论4444、污水泄入河道中问题研究的阶段划分第一阶段：污水离开排污口与周围水体掺混，以射流和浮力羽流形式扩散，一般按三维运动问题处理；第二阶段：污水还没有扩展到河流全断面，污水随河水一起运动，按二维紊动离散问题研究；第三阶段：污水已扩展到河流全断面，并且全断面完全混合，污水沿纵向继续离散，可按一维纵向离散问题分析。5555、扩散理论的基本假设（1111）扩散质的存在不改变流体质点的流动特性。即将扩散质视为标志物质或称示踪剂（TracerTracerTracerTracer）。（2222）流体质点上所带的扩散质在运动过程中保持不变，即流体质点间不发生扩散质的转移（不计分子扩散），扩散质的扩散完全是由于带有扩散质的流体质点发生掺混的结果。（3333）对不可压流体，携带扩散质流体质点的总体积在扩散过程中保持不变，扩散结果反映在携带扩散质的流体质点所占空间位置和轮廓随时间而变化。5北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论二、描述扩散运动的基本方程1111、分子扩散的费克定律（第一定律）费克（AdolphFick，1855）提出假设：盐分在其溶液中扩散的物理定律应等同于傅立叶（Fourier，1822）提出的热传导定律相同，即：其中：费克定律说明：在扩散溶液浓度场中的时空点上，单位时间内通过单位面积的扩散质的质量与该点处扩散溶液浓度的梯度成正比，比例系数为该种扩散溶液的分子扩散系数；方向与浓度梯度方向相反。分子扩散系数由扩散质及溶解质的种类、温度、压强决定，与溶液的运动形态无关，是物性参数。iiiimmmmxxxxCCCCDDDDqqqq∂∂−=—CCCC扩散液的浓度；量纲：3−MLMLMLML—qqqq在xxxxiiii方向的单位面积的扩散质量通量；量纲：12−−TTTTMLMLMLML—mmmmDDDD分子扩散系数；量纲：12−TTTTLLLL常见扩散质在水中的扩散系数表对三维情况，以矢量表示：CCCCDDDDqqqqmmmm∇−=r6北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论2222、分子扩散方程（费克第二定律）分析分子扩散应满足的控制方程，以一维情况为例，流体静止。取控制体如图，由质量守恒定律，可得：dtdtdtdtxxxxxxxxqqqqqqqqdtdtdtdtqqqqdtdtdtdtttttxxxxCCCC⋅⋅∆∂∂+−⋅⋅=⋅∂⋅∆∂1)(1)1(0=∂∂+∂∂xxxxqqqqttttCCCC整理得：将费克第一定律代入，可得：22xxxxCCCCDDDDttttCCCCmmmm∂∂=∂∂22xxxxqqqqDDDDttttqqqqmmmm∂∂=∂∂或CCCCDDDDzzzzCCCCyyyyCCCCxxxxCCCCDDDDttttCCCCmmmmmmmm2222222∇=⎟⎠⎞⎜⎝⎛∂∂+∂∂+∂∂=∂∂对三维情况：上述方程称之为分子扩散方程（DiffusionEquationsDiffusionEquationsDiffusionEquationsDiffusionEquations）。7北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论3333、移流扩散方程（AdvectiveDiffusionEquationsAdvectiveDiffusionEquationsAdvectiveDiffusionEquationsAdvectiveDiffusionEquations）取控制体如图，以xxxx1111方向为例。假设：层流运动时溶液的扩散与流体静止时的分子扩散相同。由质量守恒定律，可得：dtdtdtdtdxdxdxdxdxdxdxdxqqqqqqqqdtdtdtdtdxdxdxdxqdxqdxqdxqdxdtdtdtdtttttdxdxdxdxdxdxdxdxCdxCdxCdxCdx3232321)()(∆+−=⋅∂∂1xxxxxxxxqqqqqqqq∆∂∂=∆⎟⎟⎠⎞⎜⎜⎝⎛∂∂−+=11xxxxCCCCDDDDCuCuCuCuqqqqmmmm及21211)(xxxxCCCCDDDDxxxxCuCuCuCuttttCCCCmmmm∂∂+∂∂−=∂∂整理可得：对三维流动：⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂=∂∂+∂∂+∂∂+∂∂232222212332211)()()(xxxxCCCCxxxxCCCCxxxxCCCCDDDDxxxxCuCuCuCuxxxxCuCuCuCuxxxxCuCuCuCuttttCCCCmmmm8北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论4444、紊动扩散方程对紊动扩散，瞬时浓度和流速都可分解时均值与脉动值之和，即：代入移流扩散方程后取时均值，有：仿照分子扩散系数的表示形式，引入紊动扩散系数DDDDijijijij，令：紊动扩散方程可表示为（考虑源、汇项后）：⎟⎟⎠⎞⎜⎜⎝⎛′′−∂∂∂∂=∂∂+∂∂iiiiiiiimmmmiiiiiiiiiiiiuuuuCCCCxxxxCCCCDDDDxxxxxxxxCCCCuuuuttttCCCCuuuuuuuuuuuu′+=CCCCCCCCCCCC′+=和iiiiuuuuCCCC′′−其中：是由紊动脉动量产生的沿iiii方向的质量扩散通量。jjjjijijijijiiiixxxxCCCCDDDDuuuuCCCC∂∂=′′−ccccjjjjijijijijiiiimmmmiiiiiiiiiiiiFFFFxxxxCCCCDDDDxxxxCCCCDDDDxxxxxxxxCCCCuuuuttttCCCC+⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂∂∂=∂∂+∂∂9北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论紊动扩散质量通量与紊动扩散系数DDDDijijijij（TurbulentDiffusionTurbulentDiffusionTurbulentDiffusionTurbulentDiffusionCoefficientsCoefficientsCoefficientsCoefficients）可用矩阵表示为：DDDDijijijij应是空间坐标的函数，当选择坐标使其与二阶张量的主轴方向一致时，九个量中仅有三个主值，即：DDDD11111111，DDDD22222222，DDDD33333333不为零。当满足各向同性条件下，有：⎟⎟⎟⎟⎟⎟⎟⎠⎞⎜⎜⎜⎜⎜⎜⎜⎝⎛∂∂∂∂∂∂⎟⎟⎟⎠⎞⎜⎜⎜⎝⎛−=⎟⎟⎟⎠⎞⎜⎜⎜⎝⎛′′′′′′321333231232221131211321xxxxCCCCxxxxCCCCxxxxCCCCDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDuuuuCCCCuuuuCCCCuuuuCCCCTTTTDDDDDDDDDDDDDDDD===332211iiiiiiiiTTTTiiiiiiiixxxxxxxxCCCCDDDDxxxxCCCCuuuuttttCCCC∂∂∂=∂∂+∂∂2mmmmTTTTDDDDDDDD因为，分子扩散可忽略不计，紊动扩散方程为：10北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论三、扩散系数及其分析确定1111、分子扩散系数与概率统计量间的关系（1111）分子扩散方程的基本解问题：（考虑一维问题）在t=0t=0t=0t=0时刻，坐标原点处（x=0x=0x=0x=0）放置质量为MMMM的扩散质，确定浓度沿xxxx轴的扩散过程。22xxxxCCCCDDDDttttCCCCmmmm∂∂=∂∂基本方程：求解方法：1111）量纲分析相似解法；2222）数理方程解法定解条件：由质量守恒，在任何时刻，有：∫+∞∞−=MMMMCdxCdxCdxCdx或))))(((()))),,,,((((xMxCδ=011北京工业大学市政学科部——马长明高等流体(水）力学讲稿第五讲扩散理论量纲分析相似解法量纲分析：CCCC（xxxx，tttt）必然与MMMM、DDDDmmmm和tttt有关，有量纲关系代入扩散方程，得：得出：分离变量法结果：)4(4)(4),(ttttDDDDxxxxffffttttDDDDMMMMffffttttDDDDMMMMttttxxxxCCCCmmmmmmmmmmmmππππηηηηππππ==02=+ffffdddddfdfdfdfηηηηηηηη20)(ηηηηηηηη−=eeeeccccffff再由定解条件确定系数cccc0000。∫+∞∞−=MMMMCdxCdxCdxCdx)4exp(4),(2ttttDDDDxxxxttttDDDDMMMMttttxxxxCCCCmmmmmmmm−=ππππ∫∞+∞−−−=ξξξξξξξξππππξξξξddddttttDDDDxxxxttttDDDDCCCCttttxxxxCCCCmmmmmmmm]4)(exp[4)0,(),(212北京工业大学市政学科部——马长明高等流体(水）力学讲稿