《复杂流体讲义》由会员分享，可在线阅读，更多相关《复杂流体讲义（79页珍藏版）》请在金锄头文库上搜索。
1、现代力学与热科学进展复杂流体及其应用朱克勤清华大学航天航空学院2011年3月2021/3/111力学面临的机遇和挑战中国科学技术协会主编，力学学科发展报告，中国科学技术出版社，2007年北京引言引言2021/3/112引言引言力学与天文学是最早形成的两门自然科学。从牛顿时代开始，到十九世纪末，力学以质点、质点系、刚体、理想弹性体和理想流体为模型，运用微积分等数学工具形成了自己完整的理论体系。进入二十世纪后，力学开始以自然界和工程技术中遇到的复杂介质和复杂系统为研究对象，力学研究领域的不断开拓，一方面导致力学新分支学科不断出现，另一方面，使得力学成为现代工程技术（比如：航空航天工程、船舶工程、土
2、木工程、机械工程、热能工程和兵器工程等）的重要基础。 。2021/3/113引言引言2000年底，美国工程院评出20世纪对人类社会影响最大的20项技术，许多关键技术的进展与力学相关。以排在前3位的技术为例：1)电力系统技术。叶轮机、发电机以及输电线路的设计都离不开力学。二十世纪后50年，从力学设计导致叶轮机效率提高约1/3，经济效益达5000亿美元，而力学设计导致锅炉燃烧效率提高的经济效益也非常可观。2)汽车制造技术。力学设计使汽车发动机的效率近50年提高约1/3，仅小轿车节省的燃料费就达2000亿美元，排气污染减少90以上。3)航空技术。几乎每一阶段的重大进步均与力学家的贡献密不可分。科学技
3、术的进步永无止境。再过100年，20世纪引以为豪的技术成就只是人类现代文明的一个新的起点。2021/3/114引言引言图-144(M=2.2)协和号(M=2.02)协和号(1976-2003)由法英联合研制，它和图-144(1975-1987)同为世界上仅有的商业超音速客机。1996年2月7日，协和飞机从伦敦飞抵纽约仅耗时2小时52分59秒。2021/3/115科学大师谈力学“尽管我们今天确实知道古典力学不能用来作为统治全部物理学的基础，可是它在物理学中仍然占领着我们全部思想的中心。” A.Einstein物理学的进化“自然的一切现象，完全可以根据力学的原理用相似的推理一一演示出来。” 牛顿自
4、然哲学的数学原理1643-17271879-19552021/3/116力学与现代工程的关系“力学是航天航空的基石” 王永志 “力学搭起了基础科学与工程技术之间的桥梁” 黄克智“力学能为缓解能源短缺,提高能源利用率做出重要贡献”过增元“宇宙之大,基本粒子之小,力无所不在”杨卫“机械科学技术中的关键问题依赖力学的发展”温诗铸2021/3/117问渠那得清如许，为有源头活水来 宋朱熹观书有感流体力学的源头活水：研究对象的拓展和新研究方法的探寻引言引言2021/3/118混沌：少了一颗钉子，.丢了一个国家。Bernard对流在二十世纪初发现引言引言2021/3/119孤立波首先由S.Russell(
5、1834)在运河中发现引言引言2021/3/1110原地重现孤立波的实验（1995）引言引言2021/3/1111经典流体力学主要研究牛顿流体的运动规律和应用，二十世纪以来，近代流体流体力学迅速发展，其主要标志之一是研究对象开始从牛顿流体拓展到复杂流体。引言引言2021/3/1112问题：为什么要关注复杂流体？ICTAM2012大会将在北京举行2021/3/11132021/3/1114II.FluidPhysicsResearchThefluidphysicsprogramencompassesawiderangeofresearchinphysicsandengineeringscienc
6、e,includingstudiesofheatandmasstransferprocesses,fluiddynamics,andthephysicsofcomplexfluids.A.Complexfluids1)Colloidsandsuspensions2)Nanoscalefabricationinthefluidphase3)Granularmechanics4)Non-NewtonianfluidB.InterfacialphenomenaC.MultiphaseflowandphasechangeD.BiofluidsNASA ResearchAnnouncement2021/
7、3/11151.1复杂流体的例子泥浆火山熔岩钢水2021/3/1116血液牙膏生活中的：稀饭、果酱、酸奶、沥青、油漆、黏合剂等复杂流体有许多不同于牛顿流体的独特性质1.1复杂流体的例子同学发言：请再举出几个复杂流体的例子2021/3/1117电流变液1.2复杂流体的流动特性2021/3/1118Newtonian fluidViscoelastic fluidSpraysoffluids1.2复杂流体的流动特性2021/3/1119AsuspensionsedimentinginafluidIn a Newtonian fluidIn a viscoelastic fluid1.2复杂流体的流
8、动特性2021/3/1120DropimpactoffluidsNewtonian fluidViscoelastic fluid 1.2复杂流体的流动特性2021/3/1121T.CubaudandT.G.Mason,Folding of viscous threads in diverging microchannels,Phys.Rev.Lett.96,114501(2006).1.2复杂流体的流动特性2021/3/11221.2复杂流体的流动特性Many complex materials can not be described by simple models!GroismanA,
9、SteinbergV.EfficientmixingatlowReynoldsnumbersusingpolymeradditives.Nature,2001,410:905-8 2021/3/1123Weissenberg效应1.2复杂流体的流动特性本讲座将集中讨论复杂的粘弹性流体2021/3/11241.2复杂流体的流动特性定义：粘弹性流体是一种既具有粘性又具有弹性的介质。首先介绍粘弹性流体的几个经典模型：它们的构造方法非常简单。弹性体是固体力学的理想化模型（弹簧），粘性流体是流体力学的理想化模型（粘性阻尼器）。粘弹性流体是两者组合而成的体系。spring(Hooklaw)dashpot(
10、Newtonianfrictionlaw)2021/3/1125KelvinmodelMaxwellmodelOldroyd-Bmodel1.3粘弹性流体的经典模型问题：如何导出以上系统的应力应变关系（本构关系）？基本原则：并联:应力相加，应变相同；串联:应力相同，应变相加2021/3/11261.3粘弹性流体的经典模型2021/3/11271.3粘弹性流体的经典模型实验表明，以上经典模型过于简单，无法描述某些真实粘弹性材料的行为模式，需要探寻新的开拓方法和新的模型（源头活水）。在进一步开拓复杂粘弹流体本构关系的各种探索中，最大胆的设想由G.W.ScottBlair在1947年提出G.W.Sc
11、ottBlair,Theroleofpsychophysicsinrheology,JournalofColloidScience,1947,Vol.2,pp.21-32G.W.ScottBlair,Psychoreology:linkbetweenthepastandthepresent,JournalofTextureStudies,1974,Vol.5,pp.3-122021/3/11282.ScottBlair模型G.W.ScottBlair在他的经典论文“心理物理学在流变学中的作用”中指出，弹性体的应力与应变的零阶时间导数成正比，牛顿流体的应力与应变的一阶时间导数成正比，进一步的研究
12、则需要考虑应力与应变的分数阶导数成正比的复杂粘弹性流体。G.W.ScottBlair,Theroleofpsychophysicsinrheology,JournalofColloidScience,1947,Vol.2,pp.21-322021/3/1129Thisisathree-parametermodelandintroducedbyScottBlair.Spring(1676)Dashpot(1686)Fractionalelement (1947)分数阶导数在描述许多粘弹性材料的流变学行为中十分有效。2.ScottBlair模型2021/3/1130FractionalMaxwel
13、lfluidA.Hernndez-Jimnez,etal,Relaxation modulus in PMMA and PTFE fitting by fractional Maxwell model,Polym.Test.21(2002)325331.PolymerMethylmethacrylateMaxwellfluidPolytetrafluorethyleneThisisafour-parametermodelofviscoelasticfluids.Conclusion: fractionalelementplaysavitalroleinthedescriptionofcompl
14、exviscoelasticfluids!2.ScottBlair模型2021/3/1131Can the meaning of a derivative of integer order dny/dxn have meaning when n is 1/2?(LHospital1695)LHospital1661-1704以上内容，欢迎提问以上内容，欢迎提问2.ScottBlair模型2021/3/1132一些著名的数学大师都曾着迷于Hospital问题，比如：Euler1707-1783Fourier1768-1830Laplace1749-1827Abel1802-1829Liouvil
15、le1809-1882Riemann1826-18662.1ScottBlair模型的数学基础Riemanndevelopedadifferenttheoryoffractionaloperationsduringhisstudentdays,butitwaspublishedonlyposthumouslyin1876.ThefirstuseoffractionaloperationwasAbelin1823(21岁).2021/3/1133In1819startingwithy =xm,S.F.Lacroixpresentedhisexpressionof-orderderivativei
16、ntermsofLegendressymbolwhichdefinitionofa-orderwasintroducedbyLaplacein1812.2.1ScottBlair模型的数学基础与Laplace定义的对比Notation:then-foldintegral;then-orderderivative.2021/3/11342.1ScottBlair模型的数学基础Then-folditeratedintegraloff (t)isgivenbytheCauchysformulaForexampleTheRiemann-Liouvilleoperatoroffractionalinte
17、grationisdefinedasForexample2021/3/11352.1ScottBlair模型的数学基础ThenwegettheRiemann-LiouvilleoperatoroffractionalderivativewhichcoincideswithLaplacesdefinitionof-orderderivative.Taking =1/2andm =1intheRiemann-LiouvilleoperatoryieldsTakingm-order(misinteger)derivativegives2021/3/1136Recentlymathematically
18、fractionalcalculushasobtainedmuchsuccessinthestudyofphysicsincludingcomplexviscoelasticfluids.R.Hilfer,Applications of Fractional Calculus in Physics,WorldScientific,2000讨论：除了数学定义和运算，针对ScottBlair模型下一步应该研究的关键问题是什么？目前在聚合材料中，分数阶微积分已经成为分析应力松弛现象的一种极为重要的工具。2.1ScottBlair模型的数学基础关键问题：ScottBlair模型的力学机理和基础。202
19、1/3/11372.2ScottBlair模型的力学基础H.Schiessel&A.Blumen(1993)firstlyconstructedfractionalrheologicalconstitutiveequationsonthebasisofwellknownmechanicalmodels.H.Schiessel&A.Blumen,Hierarchical analogues to fractional relaxation equations,J.Phys.A:Math.Gen.1993,Vol.26,pp.5057-5069Spring-dashpotladder2021/3/
20、11382.2ScottBlair模型的力学基础Schiessel&Blumen利用拉氏变换，证明了系统各级弹簧和阻尼器参数满足一定递推关系时，其应变拉氏变换与应力拉氏变换服从以下关系再利用逆变换得到具体的过程将用另一个我们所提出的更加简单的例子来说明。2021/3/11392.2ScottBlair模型的力学基础HerewepresentanovelmechanicalmodeloffractionalelementQuestion: how do we obtain the constitutive equation of the tree?Spring-dashpottreewhichw
21、asenlightenedfromaresistor-capacitorself-similarstructure.I.Podlubny,Fractional Differential Equations,AcademicPress,1999.P280,Fig.10.42021/3/11402.2ScottBlair模型的力学基础Schiessel&Blumen使用的拉氏变换法，对一层的树形结构：对两层的树形结构：2021/3/11412.2ScottBlair模型的力学基础对三层的树形结构：递推求解，得到该系统的应变拉氏变换与应力拉氏变换服从以下关系2021/3/1142令右边为A，利用结构
22、层次为无穷的特点所产生的自相似性，可得到2021/3/11432.2ScottBlair模型的力学基础解得用另一种方法：Heaviside算子法逆变换得到2021/3/1144KelvinmodelMaxwellmodelOldroyd-Bmodel2.2ScottBlair模型的力学基础2021/3/1145弹簧假定系统的本构关系2.2ScottBlair模型的力学基础阻尼器总应变自相似总应力2021/3/1146Heavisideoperatorpisdisposedasaparameterduringthealgebraicoperation.2.2ScottBlair模型的力学基础20
23、21/3/11472.2ScottBlair模型的力学基础Heavisidedevelopedoperationalcalculusbetween1880and1887,whichisoneofthethreemostimportantmathematicaldiscoveriesofthelate19thCenturyandcausedmuchcontroversy.R.Courant&D.Hilbert,Methods of Mathematical Physics,Vol.IIPartialdifferentialequations,IntersciencePublishers,John
24、Wiley&Sons,1962,P507Heaviside(1850-1925)HeavisidesoperationalcalculuswasplacedonarigorousmathematicalbasisbyJanMikusinski,whoconstructedanalgebraicsettingfortheoperationalmethods.J.Mikusinski,Operational Calculus,Pergamon Press,NewYork,1983在运算微积中，算子p作为参数进行代数运算，有依据吗？2021/3/1148瞬态问题或混合问题有重要应用背景(比如:机电工
25、程)，讨论这一问题的文献很多，其中重点是Heaviside符号算子法。该方法处理问题直捷惊人，往往能给出不能以其它方法同样简单地获得的明确解答。原先发表这一方法时对于符号运算步骤并无严格道理可讲；事实上，Heaviside对职业数学家的疑虑甚至颇表不屑。然而Heaviside方法的成就压倒一切，使人们非得从数学上去弄清它的道理不可，结果完全证明这种方法有理论依据，而终于大大促进符号方法的发展。引自R.Courant和D.Hilbert的经典名著“数学物理方法”第五章附录二“瞬态问题和Heaviside运算微积”2.2ScottBlair模型的力学基础2021/3/11492.2ScottBla
26、ir模型的力学基础I.Podlubny,Fractional Differential Equations,AcademicPress,1999.(P274,Fig.10.3)homeworkKelvinmodel2021/3/11502.2ScottBlair模型的力学基础homework22021/3/11513.1圆管起动流动量方程初边值条件是速度分解为定常和非定常两部分之和为Heaviside阶跃函数 3.复杂粘弹性流体的流动2021/3/1152定常部分解的非定常部分应满足齐次方程与定常部分迭加在一起，得到3.1圆管起动流2021/3/11533.1圆管起动流2021/3/1154问
27、题：Scott-Blair模型的圆管起动流是否会有不同的特性？3.1圆管起动流2021/3/1155分数元模型动量方程其中得到3.1圆管起动流2021/3/1156用Heaviside运算微积和分数阶微积分得到其中称为Mittag-Leffler函数,是指数函数的推广,为指数函数3.1圆管起动流2021/3/11572021/3/11582021/3/11592021/3/11602021/3/1161小结小结lThe constitutive equations of spring-dashpot systems can be easily derived by operational me
28、thods. l lA new mechanical system of fractional element is presented.l lThe exact solution of starting flow of fractional element in a pipe is obtained.l lThe starting flow of fractional element in a pipe will stop finally except .2021/3/1162讨论讨论牛顿流体在静止时不能承受剪切力，为何分数元流体会出现类固体的特性？Kelvin模型Maxwell模型Oldr
29、oyd-B模型2021/3/1163讨论讨论2021/3/11643.2圆管振荡流牛顿流体3.复杂粘弹性流体的流动2021/3/1165牛顿流体3.2圆管振荡流2021/3/11662021/3/1167 Fractional Maxwell modelMaxwell model when ab 1分数阶Maxwell流体的圆管振荡流3.2圆管振荡流2021/3/1168Exact solution3.2圆管振荡流2021/3/1169Maxwellfluid PTFE3.2圆管振荡流2021/3/11703.3复杂粘弹性流体的Couette流Tan,W.C.&Xu,M.Y.,Planesur
30、facesuddenlysetinmotioninaviscoelasticfluidwithfractionalMaxwellmodel.ActaMech.Sinica(2002)18:342349W.Shaowei&X.Mingyu,ExactsolutiononunsteadyCouetteflowofgeneralizedMaxwellfluidwithfractionalderivative,ActaMechanica(2006)187:103112HaitaoQi&HuiJin,Unsteadyrotatingflowsofaviscoelasticfluidwiththefrac
31、tionalMaxwellmodelbetweencoaxialcylinders,ActaMech.Sinica(2006)22:3013052021/3/11714. 分数阶微积分在流体力学中的其它应用Sugimoto,N,PropagationofnonlinearacousticwavesinatunnelwithanarrayofHelmholtzresonators.JFM,1992,244:55-784.1在排列有一组Helmholtz共振器的隧道中，非线性声波的传播2021/3/11724.2油岩层中温度场的延伸Control Equation: fractional diff
32、usion-wave equation (Mainardi 1994)Classical diffusion equation when 2b 1Classical wave equation when 2b 2Ultraslow diffusion process when 02b 1Intermediate diffusion process when 12b 24. 分数阶微积分在流体力学中的其它应用2021/3/11732021/3/1174ChenW.Aspeculativestudyof2/3-orderfractionalLaplacianmodelingofturbulence
33、,Chaos,2006,16:023126陈文，复杂流体和复杂流动（湍流）：一个比较研究.中国力学大会2007，会议论文摘要集（下）P431它同时包含拉氏算子描述的小尺度分子粘性和2/3阶拉氏算子表达的大尺度涡粘性，并清晰地将雷诺方程和标度率联系起来4. 分数阶微积分在流体力学中的其它应用4.3湍流2021/3/11755. 结束语Although the potential of fractional constitutive equations for viscoelastic fluids remains to be fully evaluated, it is our belief t
34、hat this will be a primary direction for future developments in the quest for a predictive branch of fluid mechanics of some complex fluids.2021/3/1176分数阶导数及其应用的研究是一项具有挑战性的工作分数阶微积分类似的应用在其它科学和技术领域不断被发现，包括电子技术、信号分析、控制理论、分形与混沌、热力学、生物物理和流变学等。问题：常数 c 的 阶导数是否应该等于零?5. 结束语2021/3/1177TheCaputosfractionalderivativesisdefinedasComparewithRiemann-Liouvilleoperatoroffractionalderivative5. 结束语2021/3/1178本讲座结束！2021/3/1179
