资源描述
DearSongandZiniu,IwouldliketorecruitagraduatestudentwhohasundergraduateeducationfromTsinghuaUniversityunderyoutwo.Isitpossiblethatyoucanhelpmetofindone?IhopethathecanspendcoupleyearsinPekingUniversityforMasterDegreeandcomeUS(orstayatPKU)twoorthreeyearslaterforPh.D.IassumethatheshouldbequalifiedforthewaiveroftheentryexaminationfromTsinghua.Inthisway,IthinkthatwecantransferthewaivertoPekingUniversity.TheideaofgettingagraduatestudentfromTsinghuaissoobvious.Iwouldliketoencouragethemixingofacademicculture.Ifyouthinkthatisagoodidea,weshoulddothisassoonaspossible.Thetimeisup.Ofcourse,ifwecannotfindsuchacandidatethisyear,wecantrynextyear.Theresearchareaswillinclude:TurbulencetheoryandComputations,CFDandapplications,inparticularincompbustionandmultiphaseflows,nanofluidicsandmicrofluidics,pollentransportinbiocomplxsystemsandengineeringdataminingrelatedinformationtechnology.HopetohearfromyousoonandBestregards,Shiyi意义许多实际问题属于一维非定常流动：如汽车进排气管道中的流动。某些三维流动（如点爆炸）也可以看成一维流动，如柱面波和球面波问题许多高维问题局部存在一维效应研究一维非定常问题可以揭示许多流动现象，因为问题存在精确解VI:非定常流动问题特征考虑大扰动波或有限振幅波的传播，属于非线性问题该类问题存在一些特殊解，如膨胀波、压缩波和激波；需要了解这些解的运动规律及它们之间或与其它物体之间的相互作用规律VI:非定常流动研究内容基本方程及特征线法一维非定常均熵流动间断流波的反射与相互作用其它问题VI:非定常流动VI-1:基本方程与特征线法基本方程特征线方程相容关系式黎曼不变量与简单波VI1:基本方程方程组推导：控制体VI-1:基本方程几何说明截面积,控制体体积对于周长为的非圆截面管道，可以定义水利学直径因此可以看成是直径为D的当量圆截面管道所有流动参数都用截面平均值扩张角满足VI-1:基本方程质量守恒方程左边界质量流量右边界质量流量添质作用控制体质量变化率质量守恒关系式VI-1:基本方程连续性方程由得动量守恒：摩擦力表达式壁面剪切应力这里f 为摩阻系数（无量纲，基本为常数）管壁摩擦力管壁摩擦力沿轴线投影为利用，得VI-1:基本方程动量方程：总受力分析左边界受力：，右边界受力：侧面受力：合力VI-1:基本方程动量方程：动量流量从左边流进的动量从右边流出的动量从侧面添质带进的动量：动量流量引起的净增加VI-1:基本方程动量守恒方程外力和净动量流量引起控制体内动量增加因此，动量方程为VI-1:基本方程动量方程由得VI-1:基本方程能量方程：守恒率V-1:基本方程能量方程形式一般形式对于量热完全气体，可以得出V-1:基本方程广义一维流动控制方程组VI-1:基本方程矩阵形式矩阵形式如何求解？方程为拟线性偏微分方程组，传统方法无法求解。比较有效的方法是采用特征线法。方程沿某些曲线（特征线）可以化成常微分方程，从而可以积分出来。因此考虑波沿这些曲线的传播。特征线理论按得特征值相应的左特征向量为相容关系式1对于特征线，相容关系式为相容关系式2对于特征线，相容关系式为相容关系式3对于特征线，相容关系式为简单流动的相容关系式对于的简单流动，有VI-2:一维非定常均熵流动黎曼不变量与简单波膨胀波与压缩波的定义中心稀疏波VI2:均熵流动黎曼不变量对于量热完全气体()和等熵流动可以将第1和第3个相容关系式沿特征线积分，得沿第2条特征线，熵为常数；在熵为常数的前提下，沿第1和第3条特征线，分别不变。称为黎曼不变量。也就是说，沿给定的（第1和第3条）特征线，黎曼不变量为常数，但对于不同特征线，它们的值可以不同。为了方便，记VI2:均熵流动黎曼不变量的作用根据黎曼不变量的定义，沿特征线方向的任意点的解可以由某一初始点的解获得。特征值的正负反应了波的传播方向。不变量可以用于定义边界条件。例如，在入口边界，如果，那么对应的不变量要给定。VI2:均熵流动边界条件个数VI2:均熵流动左（右）边边界条件类型设有个正(负)特征值，需要给定个不变量作为边界条件。对应正(负)特征值的波称为入射波，对应负(正)特征值的波称为出口波。无反射边界条件：无论出口波如何，只给定入射波的值作为边界条件反射边界条件：出口波被部分反射为入口波VI2:均熵流动例：亚音速入流边界条件反射边界条件(一般边界条件)也可以按扰动量写成考虑到因此上述一般边界条件可以写成VI2:均熵流动病态条件特例如果在亚音速入口给定速度和压力，即令则一般边界条件成为因此，不存在与无关的参数，使(3)等价于无反射边界条件的一般形式（2）。VI2:均熵流动良态条件特例如果在亚音速入口给定密度和速度，即令则一般边界条件成为因此，取，上述条件便变为VI2:均熵流动简单波定义1：如果流动所涉及的两个或多个参数互为单值函数，那么流动称为简单波流动。它比均匀流区域稍复杂，比非简单波区域更简单。定义2：在简单波区域，其中一个黎曼不变量为常数（从而是两个，因为已经假定熵不变)。可以证明(p.253)，以上两个定义是等价的。数学上，简单波有更一般的定义。这些定义正好等价。VI2:均熵流动状态平面上的简单波平面(V,a)称为状态平面。均匀流对应状态平面上一点，简单波对应一条直线，非简单波有两束特征线（在状态平面上都为直线）。VI2:均熵流动简单波的连接（特征平面）与均匀流相连接的区域只可能是简单波区域。简单波区既可与均匀流区相接，也可以与非简单波区相接。均匀流简单波区非简单波区均匀流VI2:均熵流动左行简单波（站在运动质点上看）对于第一个黎曼不变量为常数的简单波区域，有在第I束特征线上，第二个黎曼不变量也为常数因此，在第I束特征线上，即第I束特征线为直线，从而有VI2:均熵流动左行简单波一般情况下，c=c(V)orc(a).如果c0为常数,则波是从一个点出发的，称为以为中心的中心简单波。简单波可分为膨胀波（稀疏波）与压缩波两种VI2:均熵流动左行简单波图示VI2:均熵流动左行中心简单波关系式假设中心简单波区域由定义，由（1）简单波的定义得另外由简单波特征线的性质得即，从而VI2:均熵流动左行中心简单波关系式续假设中心简单波左右均匀流状态分别为和，它们满足条件于是简单波区的定义为在简单波区任意点(x,t)各参数定义为VI2:均熵流动左行中心简单波的性质由于所以左行中心简单波必为膨胀波，从而被称为中心稀疏波(expansionfan).问题：中心稀疏波与坐标选择有关吗？(Z.N.Wu,Expansionfaninforcedandpassivelymovingframes,ActaMechanica,acceptedandtoappearin2019)VI2:均熵流动习题请结合前面介绍和书中254-256页，将上述内容推广到右行简单波的情况。VI2:均熵流动活塞后撤问题设活塞后撤速度为。引起右行膨胀波。VI2:均熵流动活塞后撤问题对于右行简单波(3-简单波)，有因此从而由得VI2:均熵流动活塞后撤问题于是由上式解出V得VI2:均熵流动VI-3:有间断的流动激波与接触间断兰金于戈尼奥(Rankine-Hugoniot)关系式（激波跳跃关系式）。以一般间断作为初始条件的问题：黎曼问题波的相互作用VI3:有间断的流动激波的产生：活塞效应左边活塞向右加速至某一固定速度越靠左，扰动传播速度就越大，从而左边的波试图超越右边的波，即波形越来越陡峭。VI3:有间断的流动激波的产生续波形越陡峭，局部速度梯度（还有温度梯度）就越大，从而粘性耗散作用就越大，直到出现粘性耗散与波的追赶效应出现平衡为止。出现平衡时波形内部高梯度区所对应的厚度不再变化，为几个分子平均自由程的量级，宏观上可以看成间断，也称为激波。激波厚度为1/10个微米的量级。激波内部有真实气体效应。VI3:有间断的流动参数与坐标系激波左边参数和右边参数分别记为和，另外激波的运动速度记为实验室坐标系：在该坐标系中，管道（或其它固定物体）不动。在前面提到的例子中，激波一侧的气体是静止的。一般坐标系：在一般坐标系中，两侧气体都有速度。激波坐标系：坐标系原点固定于激波上，在该坐标系中，激波不动，即激波速度为0。VI3:有间断的流动参数关系式设一般坐标系中的流动参数为那么在激波坐标系中，相应的流动参数为VI3:有间断的流动守恒关系式坐标系：选择激波坐标系。守恒关系式：左边的（质量、动量与能量）通量，等于右边的通量,即VI3:有间断的流动兰金于戈尼奥(R-H)关系式利用和前面的守恒关系式得R-H关系式利用记号，则R-H(Rankine-Hugoniot)关系式(也称激波跳跃关系式)也可以写成VI3:有间断的流动R-H关系式的显示表达对于气体动力学方程，R-H关系式可以写为VI3:有间断的流动历史上的兰金于戈泥奥关系式由历史牛顿于1687年用等温假设研究声音的传播。Monge于1770年提出特征线理论，被Earnshaw(1858)和Riemann(1859)独立用于气体动力学问题。Poisson于1808年首次考虑有限振幅的波动（可看成激波理论的首作）。同时他定义了气体比热比的符号。Laplace于1816年改用等熵假设改进了牛顿的理论。Earnshaw(1851)发现雷电的超音速传播现象。激波一词首先由Toepler于1864年提出。DeLaval于1888年发明能产生超音速出流的收缩扩张管道（称为拉瓦尔喷管）。Vieille于1899年发明激波管，通过弄破分开高压气体和低压气体的膜片，获得强激波并产生高温高速气流。VI3:有间断的流动历史续接着Poisson的工作，Airy(1848),Stokes(1849),Rankine(1858,1870),Earnshaw(1958-1860),Riemann(1859),Hugoniot(1885-1887),Rayleigh(1910),Taylor(1910),争吵着发展和完善了经典激波理论。Riemann(1859)认为穿越激波熵不变。Rankine(1869)与Hugoniot(1887)分别独立地发现穿越激波质点熵增加。Rankine(1869)考虑热传导。Hugoniot(1887)发现，在p-v相平面上，激波参数不满足绝热曲线，而满足动态绝热曲线（被称为HugoniotCurve）。并且导出上式被称为（Rankine-）Hugoniot关系式。VI3:有间断的流动激波坐标系中R-H关系式的解考虑量热完全气体，R-H关系式为VI3:有间断的流动绝热线和Hugoniot线VI3:有间断的流动激波坐标系中普朗特公式由得另外由临界音速的定义有从而VI3:有间断的流动注意临界音速是这样定义的，气流到某一截面，其速度等于当地音速。用的是总焓守恒这一关系式由于跨越（定常）激波总焓守恒，所以左边的临界音速等于右边的临界音速。但推广到非定常激波，需要特别注意。VI3:有间断的流动普朗特公式意义可以写成。分两种情况：如果，因激波是压缩波，所以因此波前必定是超音速（），波后是亚音速（）。如果，因激波是压缩波，所以因此波前必定是超音速（），波后是亚音速（）。VI3:有间断的流动激波前后马赫数关系式由速度系数关系式和得VI3:有间断的流动激波前后密度与速度关系式由得VI3:有间断的流动激波前后压力与温度关系式由得VI3:有间断的流动激波前后总压关系式由得VI3:有间断的流动激波前后总温总密度关系由于并且所以即激波前后总温不变于是VI3:有间断的流动激波前后熵的变化由VI3:有间断的流动激波坐标系中定性分析“如果”，那么如果，则称为弱激波，因为VI3:有间断的流动一般坐标系将前面的速度用相对速度，马赫数用相对马赫数代替，则所有关系式维持不变。另外定义压力比关系不变普朗特关系VI3:有间断的流动一般坐标系马赫数关系速度关系VI3:有间断的流动一般坐标系压力密度关系总压关系VI3:有间断的流动速度压力关系式由和消去马赫数得右行和左行激波的速度压力关系式与某状态相连的激波的p-V关系式与状态相连的右行激波满足与状态相连的左行激波满足一般坐标系熵关系式VI3:有间断的流动多激波定理设左边马赫数和压力给定，经过若干道激波使气流压力增加到某一给定值。那么，经过一道激波比经过多道激波引起的总压损失更大，熵增更大。VI3:有间断的流动普朗特公式意义可以写成或已经知道，在亚音速区，马赫数小于速度系数；在超音速区，马赫数大于速度系数。分两种情况：VI3:有间断的流动普朗特公式意义续如果，因激波是压缩波，所以因此即因此，第I束特征线向左行激波会聚。VI3:有间断的流动普朗特公式意义续如果，因激波是压缩波，所以因此即因此，第III束特征线向右行激波会聚。VI3:有间断的流动同族激波性质第I族特征线会聚的激波为I族激波，第III族特征线会聚的激波为III族激波。同族激波也可以说成是同向激波。也可以说成是左行（第I族、反向）激波或右行(第III族、正向)激波。同族激波追赶定理：对于两个相邻的同族激波，后面的激波比前面的运动的快，必然赶上前面的。请用两种方法证明。VI3:有间断的流动接触间断前面提到的2简单波为接触间断，即速度和压力连续，密度（从而熵）有间断。间断以质点速度运动。物理上也称为滑移线（面）。两相流的交界面为接触间断。数学表达式为思考题：证明接触间断也满足R-H关系式VI3:有间断的流动基本解的特征线表示VI3:有间断的流动思考题：唯一性问题考虑一以速度运动的激波，左右状态分别为，请问该激波是否可以用多个波（两个激波、一个激波加一个稀疏波等）代替？（激波表达唯一性问题）。VI3:有间断的流动思考题：稳定性问题考虑一以速度运动的激波，左右状态分别为，请问该激波在什么条件下是稳定的，也就是说，给定激波一个轻微的扰动，激波是恢复原样？还是分解成其它波（或变得不成样子）？（激波稳定性问题）。VI3:有间断的流动黎曼问题定义前面介绍了激波与接触间断，间断前后参数都满足R-H关系式。如果给定一般间断，并且和不满足R-H关系式，即不存在，使得下式满足那么将出现说明情况？解是什么？这就是黎曼问题。黎曼问题在激波管（书p272-274-）等问题中有应用,另外大量用于计算流体力学。VI3:有间断的流动黎曼问题的思考给定间断左边的状态和右边的状态。如果是激波，那么激波速度是未知数。而激波跳跃关系式有三个关系式。因此，未知数个数(1)少于关系式(3)个数，即问题为超定问题，一般无解，即任意给定的间断一般不为激波。既然一般初始间断不是一个激波，那么有可能立即分解为几个基本解，如两个激波或更多、激波加接触间断或激波加中心稀疏波。这些解都是最简单的解。物质世界越简单越稳定。到底是何种组合，首先要考虑问题的确定性。VI3:有间断的流动黎曼问题的思考续假设初始间断分解为两个激波，它们中间的（均匀流）状态记为。两个激波各有一个速度。因此给定左边和右边流动参数后，有3个流动参数和2个激波速度。而两个激波的跳跃关系式加起来有6个。于是未知数个数（5个）少于于方程个数（6个）。于是问题变成超定的，一般无解。初始间断不可能分解为三个激波。这是因为，如果有三个，那么必然有两个属于同族激波，从而前者被后者赶上。刚分开又和好，自然界不存在。VI3:有间断的流动黎曼问题的思考续假设分解为三个波（激波、接触间断或中心稀疏波），则中间出现两段均匀流区，共6个未知数。如果中间波为接触间断，那么左右波无论是激波还是中心稀疏波，都可以使未知数个数等于方程个数（课堂解释）。VI3:有间断的流动黎曼问题的思考续因此，一般情况下，黎曼问题的解为一个向右的激波或中心稀疏波、中间一个接触间断、一个向左的激波或中心稀疏波。特殊情况下上述波的数量减少到两个，更特殊的情况减少到一个。VI3:有间断的流动黎曼问题思考续思考题：能否出现4个或以上的波的情况？（证明不能）具体出现何种情况，还得根据左右流动状态决定。利用穿越各波的关系式，消去参数，剩下一个，看是否有合理的解。对于给定的左右状态，合理的解一般是唯一的。目前有通用软件求解黎曼问题的精确解（判断出现何种情况和给出相应的解）。VI3:有间断的流动黎曼问题精确解给定初始间断的左右状态，确定t0时刻的波态、波的强度和波与波之间的流动特性的问题称为黎曼问题。VI3:有间断的流动黎曼问题的精确解续五种波态:RCS;SCR;SCS;RCR;RCVCRVI3:有间断的流动黎曼问题的精确解续五种波态边界(p-V图，N表示退化)具体边界见计算流体力学基本原理第111页.VI3:有间断的流动RCS情况左移稀疏波关系式右移激波关系式VI3:有间断的流动RCS情况续由得VI3:有间断的流动RCS情况续由得VI3:有间断的流动RCS情况续由得从而VI3:有间断的流动RCS情况续由两者相减得上式即为确定压力的代数方程（迭代法求）。VI3:有间断的流动RCS情况续得到压力后，由等熵关系式确定密度，再由确定速度，再由确定密度VI3:有间断的流动思考题用上述类似思路，考虑下面问题的解:SCR问题RCR问题SCS问题.VI3:有间断的流动黎曼问题得一般公式求和的一般公式为参考文献:E.F.Toro,RiemannSolversandNumericalMethodsforFluidFlows,Springer,2019.VI3:有间断的流动黎曼问题得一般公式这里，由下式定义VI3:有间断的流动黎曼问题得一般公式这里，由下式定义VI3:有间断的流动黎曼问题得一般公式系数由下面式子定义VI3:有间断的流动波的相互作用波(简单波、激波、接触间断)在物理平面上的反射。波的相互作用。激波的反射、衍射、折射（略）。激波的合并与分裂（两个同向激波可以合并成一个吗？一个激波可以分裂成两个吗？略）。VI3:有间断的流动膨胀波在固体壁面上的反射考虑右行膨胀波（），在静止气体中向右传播，在右边固体壁面上反射。反射时相当于静止气体突然向左以某速度运动，由于壁面处速度保持为零，所以仍有，即右行膨胀波在固壁上反射为左行膨胀波。如果是中心简单波，则反射规律见p256-257.VI3:有间断的流动压缩波在固体壁面上的反射考虑右行压缩波（），在静止气体中向右传播，在右边固体壁面上反射。反射时相当于静止气体突然向右以某速度运动，由于壁面处速度保持为零，所以仍有，即右行压缩波在固壁上反射为左行压缩波。VI3:有间断的流动激波在固体壁面上的反射考虑右行激波，在静止气体中向右传播，在右边固体壁面上反射。反射时相当于静止气体突然向右以某有限速度运动，由于壁面处速度保持为零，所以反射波必为激波。反射规律见p266-267。VI3:有间断的流动波在开口端的反射与固体壁面速度恒为0相反，在开口端（大气压）压力恒定。右行膨胀波（左边压力低）反射为左行压缩波（由于右端开口端压力固定，所以左边压力低的状态仍然保持，对于左行波，质点穿越波时被压缩）。同理，右行压缩波和激波都反射为膨胀波。VI3:有间断的流动等熵波的相互作用异族膨胀波相互碰撞，其透射波仍为膨胀波异族压缩波相互碰撞，其透射波仍为压缩波膨胀波和压缩波相互膨胀，压缩波的透射波仍为压缩波，膨胀波的透射波仍为膨胀波。原因（极限思考法）：考虑等熵波的极限情况即小扰动情况，此时各波的相互干扰近似为线性干扰，相互近似无干扰。严格说明请参阅p258-259.VI3:有间断的流动非等熵波与任意波的相互作用异族激波的相互碰撞。激波与膨胀波的相互碰撞。激波与接触间断的相互碰撞。膨胀波与接触间断的相互碰撞。同族激波的超越。激波超越膨胀波。膨胀波超越激波。VI3:有间断的流动一般求解法则从碰撞或接触开始，可以假定两个波之间的中间状态消失，以最左边的状态和最右边的状态定义一个黎曼问题，黎曼问题的解即为碰撞或接触后的解。因此，碰撞或接触后，一般出现一个向左的膨胀波或激波、一个处于中间的接触间断和一个向右的激波或膨胀波。但有些情况更特殊（更简单）。VI3:有间断的流动特殊规律一异族激波的相互碰撞，各自的透射波均为激波，中间出现接触间断(p268)。特殊规律二激波与膨胀波相互碰撞，激波的透射波为激波，膨胀波的透射波为膨胀波，中间出现接触间断。特殊规律三任何波与接触间断相互作用，透射波性质不变，接触间断的透射波仍为接触间断，同时产生反射波。因此接触间断有固体边界、开口边界和波的三重作用。激波与接触间断相互碰撞，激波的透射波为激波，激波的反射波可以是激波也可以是膨胀波，接触间断的透射波为接触间断。特殊规律四同性波相互追赶、或者强波追赶异性波，追赶波的透射波与入射波同类型。例：激波追赶激波，追赶激波的透射波仍为激波，但反射波可以是激波也可以是膨胀波。例：强激波追赶膨胀波，追赶激波的透射波为激波，反射波可以是膨胀波，也可以是激波。例：强膨胀波追赶激波，追赶膨胀波的透射波为膨胀波，反射波可以是膨胀波，也可以是激波。特殊规律五弱波追赶异性波，透射波的性质改变例：弱激波追赶膨胀波，其透射波变为膨胀波。例：弱膨胀波追赶激波，其透射波变为激波。讨论题证明：两个异族激波相互碰撞，各自的透射波必为激波证明碰撞后，两个激波的透射波可能是激波，也可能是膨胀波。中间状态还可能出现接触间断。证明续与状态1相连的右行激波满足与状态1相连的左行激波满足证明续与状态1相连的右行膨胀波满足与状态1相连的左行膨胀波满足证明续与状态2相连的右行激波满足与状态2相连的左行激波满足证明续与状态2相连的右行膨胀波满足与状态2相连的左行膨胀波满足在V-p相图上的碰撞显然，由于激波A向右移，激波B向左移，有在V-p相图上的透射情形1透射激波a向左移，透射激波b向右移在V-p相图上的透射情形2透射膨胀a向左移，透射激波b向右移在V-p相图上的透射情形3透射激波a向左移，透射膨胀波b向右移结论只有透射波均为激波时，才存在中间状态4，能与原来的状态1和状态2相连。因此，异性激波相互碰撞，各自的透射波均为激波。证明完毕。谢谢！
展开阅读全文