复杂流体讲义PPT课件.ppt

现代力学与热科学进展复杂流体及其应用朱克勤清华大学航天航空学院2011年3月2024/1/27周六1.力学面临的机遇和挑战中国科学技术协会主编，力学学科发展报告，中国科学技术出版社，2007年北京引言引言2024/1/27周六2.引言引言力学与天文学是最早形成的两门自然科学。从牛顿时代开始，到十九世纪末，力学以质点、质点系、刚体、理想弹性体和理想流体为模型，运用微积分等数学工具形成了自己完整的理论体系。进入二十世纪后，力学开始以自然界和工程技术中遇到的复杂介质和复杂系统为研究对象，力学研究领域的不断开拓，一方面导致力学新分支学科不断出现，另一方面，使得力学成为现代工程技术（比如：航空航天工程、船舶工程、土木工程、机械工程、热能工程和兵器工程等）的重要基础。2024/1/27周六3.引言引言2000年底，美国工程院评出20世纪对人类社会影响最大的20项技术，许多关键技术的进展与力学相关。以排在前3位的技术为例：1)电力系统技术。叶轮机、发电机以及输电线路的设计都离不开力学。二十世纪后50年，从力学设计导致叶轮机效率提高约1/3，经济效益达5000亿美元，而力学设计导致锅炉燃烧效率提高的经济效益也非常可观。2)汽车制造技术。力学设计使汽车发动机的效率近50年提高约1/3，仅小轿车节省的燃料费就达2000亿美元，排气污染减少90以上。3)航空技术。几乎每一阶段的重大进步均与力学家的贡献密不可分。科学技术的进步永无止境。再过100年，20世纪引以为豪的技术成就只是人类现代文明的一个新的起点。2024/1/27周六4.引言引言图-144(M=2.2)协和号(M=2.02)协和号(1976-2003)由法英联合研制，它和图-144(1975-1987)同为世界上仅有的商业超音速客机。1996年2月7日，协和飞机从伦敦飞抵纽约仅耗时2小时52分59秒。2024/1/27周六5.科学大师谈力学“尽管我们今天确实知道古典力学不能用来作为统治全部物理学的基础，可是它在物理学中仍然占领着我们全部思想的中心。”A.Einstein物理学的进化“自然的一切现象，完全可以根据力学的原理用相似的推理一一演示出来。”牛顿自然哲学的数学原理1643-17271879-19552024/1/27周六6.力学与现代工程的关系“力学是航天航空的基石”王永志 “力学搭起了基础科学与工程技术之间的桥梁”黄克智“力学能为缓解能源短缺,提高能源利用率做出重要贡献”过增元“宇宙之大,基本粒子之小,力无所不在”杨卫“机械科学技术中的关键问题依赖力学的发展”温诗铸2024/1/27周六7.问渠那得清如许，为有源头活水来 宋朱熹观书有感流体力学的源头活水：研究对象的拓展和新研究方法的探寻引言引言2024/1/27周六8.混沌：少了一颗钉子，.丢了一个国家。Bernard对流在二十世纪初发现引言引言2024/1/27周六9.孤立波首先由S.Russell(1834)在运河中发现引言引言2024/1/27周六10.原地重现孤立波的实验（1995）引言引言2024/1/27周六11.经典流体力学主要研究牛顿流体的运动规律和应用，二十世纪以来，近代流体流体力学迅速发展，其主要标志之一是研究对象开始从牛顿流体拓展到复杂流体。引言引言2024/1/27周六12.问题：为什么要关注复杂流体？ICTAM2012大会将在北京举行2024/1/27周六13.2024/1/27周六14.II.FluidPhysicsResearchThefluidphysicsprogramencompassesawiderangeofresearchinphysicsandengineeringscience,includingstudiesofheatandmasstransferprocesses,fluiddynamics,andthephysicsofcomplexfluids.A.Complexfluids1)Colloidsandsuspensions2)Nanoscalefabricationinthefluidphase3)Granularmechanics4)Non-NewtonianfluidB.InterfacialphenomenaC.MultiphaseflowandphasechangeD.BiofluidsNASA ResearchAnnouncement2024/1/27周六15.1.1复杂流体的例子泥浆火山熔岩钢水2024/1/27周六16.血液牙膏生活中的：稀饭、果酱、酸奶、沥青、油漆、黏合剂等复杂流体有许多不同于牛顿流体的独特性质1.1复杂流体的例子同学发言：请再举出几个复杂流体的例子2024/1/27周六17.电流变液1.2复杂流体的流动特性2024/1/27周六18.Newtonian fluidViscoelastic fluidSpraysoffluids1.2复杂流体的流动特性2024/1/27周六19.AsuspensionsedimentinginafluidIn a Newtonian fluidIn a viscoelastic fluid1.2复杂流体的流动特性2024/1/27周六20.DropimpactoffluidsNewtonian fluidViscoelastic fluid 1.2复杂流体的流动特性2024/1/27周六21.T.CubaudandT.G.Mason,Folding of viscous threads in diverging microchannels,Phys.Rev.Lett.96,114501(2006).1.2复杂流体的流动特性2024/1/27周六22.1.2复杂流体的流动特性Many complex materials can not be described by simple models!GroismanA,SteinbergV.EfficientmixingatlowReynoldsnumbersusingpolymeradditives.Nature,2001,410:905-8 2024/1/27周六23.Weissenberg效应1.2复杂流体的流动特性本讲座将集中讨论复杂的粘弹性流体2024/1/27周六24.1.2复杂流体的流动特性定义：粘弹性流体是一种既具有粘性又具有弹性的介质。首先介绍粘弹性流体的几个经典模型：它们的构造方法非常简单。弹性体是固体力学的理想化模型（弹簧），粘性流体是流体力学的理想化模型（粘性阻尼器）。粘弹性流体是两者组合而成的体系。spring(Hooklaw)dashpot(Newtonianfrictionlaw)2024/1/27周六25.KelvinmodelMaxwellmodelOldroyd-Bmodel1.3粘弹性流体的经典模型问题：如何导出以上系统的应力应变关系（本构关系）？基本原则：并联:应力相加，应变相同；串联:应力相同，应变相加2024/1/27周六26.1.3粘弹性流体的经典模型2024/1/27周六27.1.3粘弹性流体的经典模型实验表明，以上经典模型过于简单，无法描述某些真实粘弹性材料的行为模式，需要探寻新的开拓方法和新的模型（源头活水）。在进一步开拓复杂粘弹流体本构关系的各种探索中，最大胆的设想由G.W.ScottBlair在1947年提出G.W.ScottBlair,Theroleofpsychophysicsinrheology,JournalofColloidScience,1947,Vol.2,pp.21-32G.W.ScottBlair,Psychoreology:linkbetweenthepastandthepresent,JournalofTextureStudies,1974,Vol.5,pp.3-122024/1/27周六28.2.ScottBlair模型G.W.ScottBlair在他的经典论文“心理物理学在流变学中的作用”中指出，弹性体的应力与应变的零阶时间导数成正比，牛顿流体的应力与应变的一阶时间导数成正比，进一步的研究则需要考虑应力与应变的分数阶导数成正比的复杂粘弹性流体。G.W.ScottBlair,Theroleofpsychophysicsinrheology,JournalofColloidScience,1947,Vol.2,pp.21-322024/1/27周六29.Thisisathree-parametermodelandintroducedbyScottBlair.Spring(1676)Dashpot(1686)Fractionalelement(1947)分数阶导数在描述许多粘弹性材料的流变学行为中十分有效。2.ScottBlair模型2024/1/27周六30.FractionalMaxwellfluidA.Hernndez-Jimnez,etal,Relaxation modulus in PMMA and PTFE fitting by fractional Maxwell model,Polym.Test.21(2002)325331.PolymerMethylmethacrylateMaxwellfluidPolytetrafluorethyleneThisisafour-parametermodelofviscoelasticfluids.Conclusion:fractionalelementplaysavitalroleinthedescriptionofcomplexviscoelasticfluids!2.ScottBlair模型2024/1/27周六31.Can the meaning of a derivative of integer order dny/dxn have meaning when n is 1/2?(LHospital1695)LHospital1661-1704以上内容，以上内容，欢迎提迎提问2.ScottBlair模型2024/1/27周六32.一些著名的数学大师都曾着迷于Hospital问题，比如：Euler1707-1783Fourier1768-1830Laplace1749-1827Abel1802-1829Liouville1809-1882Riemann1826-18662.1ScottBlair模型的数学基础Riemanndevelopedadifferenttheoryoffractionaloperationsduringhisstudentdays,butitwaspublishedonlyposthumouslyin1876.ThefirstuseoffractionaloperationwasAbelin1823(21岁).2024/1/27周六33.In1819startingwithy=xm,S.F.Lacroixpresentedhisexpressionof-orderderivativeintermsofLegendressymbolwhichdefinitionofa-orderwasintroducedbyLaplacein1812.2.1ScottBlair模型的数学基础与Laplace定义的对比Notation:then-foldintegral;then-orderderivative.2024/1/27周六34.2.1ScottBlair模型的数学基础Then-folditeratedintegraloff(t)isgivenbytheCauchysformulaForexampleTheRiemann-LiouvilleoperatoroffractionalintegrationisdefinedasForexample2024/1/27周六35.2.1ScottBlair模型的数学基础ThenwegettheRiemann-LiouvilleoperatoroffractionalderivativewhichcoincideswithLaplacesdefinitionof-orderderivative.Taking=1/2andm=1intheRiemann-LiouvilleoperatoryieldsTakingm-order(misinteger)derivativegives2024/1/27周六36.Recentlymathematicallyfractionalcalculushasobtainedmuchsuccessinthestudyofphysicsincludingcomplexviscoelasticfluids.R.Hilfer,Applications of Fractional Calculus in Physics,WorldScientific,2000讨论：除了数学定义和运算，针对ScottBlair模型下一步应该研究的关键问题是什么？目前在聚合材料中，分数阶微积分已经成为分析应力松弛现象的一种极为重要的工具。2.1ScottBlair模型的数学基础关键问题：ScottBlair模型的力学机理和基础。2024/1/27周六37.2.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-dashpotladder2024/1/27周六38.2.2ScottBlair模型的力学基础Schiessel&Blumen利用拉氏变换，证明了系统各级弹簧和阻尼器参数满足一定递推关系时，其应变拉氏变换与应力拉氏变换服从以下关系再利用逆变换得到具体的过程将用另一个我们所提出的更加简单的例子来说明。2024/1/27周六39.2.2ScottBlair模型的力学基础HerewepresentanovelmechanicalmodeloffractionalelementQuestion:how do we obtain the constitutive equation of the tree?Spring-dashpottreewhichwasenlightenedfromaresistor-capacitorself-similarstructure.I.Podlubny,Fractional Differential Equations,AcademicPress,1999.P280,Fig.10.42024/1/27周六40.2.2ScottBlair模型的力学基础Schiessel&Blumen使用的拉氏变换法，对一层的树形结构：对两层的树形结构：2024/1/27周六41.2.2ScottBlair模型的力学基础对三层的树形结构：递推求解，得到该系统的应变拉氏变换与应力拉氏变换服从以下关系2024/1/27周六42.令右边为A，利用结构层次为无穷的特点所产生的自相似性，可得到2024/1/27周六43.2.2ScottBlair模型的力学基础解得用另一种方法：Heaviside算子法逆变换得到2024/1/27周六44.KelvinmodelMaxwellmodelOldroyd-Bmodel2.2ScottBlair模型的力学基础2024/1/27周六45.弹簧假定系统的本构关系2.2ScottBlair模型的力学基础阻尼器总应变自相似总应力2024/1/27周六46.Heavisideoperatorpisdisposedasaparameterduringthealgebraicoperation.2.2ScottBlair模型的力学基础2024/1/27周六47.2.2ScottBlair模型的力学基础Heavisidedevelopedoperationalcalculusbetween1880and1887,whichisoneofthethreemostimportantmathematicaldiscoveriesofthelate19thCenturyandcausedmuchcontroversy.R.Courant&D.Hilbert,Methods of Mathematical Physics,Vol.IIPartialdifferentialequations,IntersciencePublishers,JohnWiley&Sons,1962,P507Heaviside(1850-1925)HeavisidesoperationalcalculuswasplacedonarigorousmathematicalbasisbyJanMikusinski,whoconstructedanalgebraicsettingfortheoperationalmethods.J.Mikusinski,Operational Calculus,Pergamon Press,NewYork,1983在运算微积中，算子p作为参数进行代数运算，有依据吗？2024/1/27周六48.瞬态问题或混合问题有重要应用背景(比如:机电工程)，讨论这一问题的文献很多，其中重点是Heaviside符号算子法。该方法处理问题直捷惊人，往往能给出不能以其它方法同样简单地获得的明确解答。原先发表这一方法时对于符号运算步骤并无严格道理可讲；事实上，Heaviside对职业数学家的疑虑甚至颇表不屑。然而Heaviside方法的成就压倒一切，使人们非得从数学上去弄清它的道理不可，结果完全证明这种方法有理论依据，而终于大大促进符号方法的发展。引自R.Courant和D.Hilbert的经典名著“数学物理方法”第五章附录二“瞬态问题和Heaviside运算微积”2.2ScottBlair模型的力学基础2024/1/27周六49.2.2ScottBlair模型的力学基础I.Podlubny,Fractional Differential Equations,AcademicPress,1999.(P274,Fig.10.3)homeworkKelvinmodel2024/1/27周六50.2.2ScottBlair模型的力学基础homework22024/1/27周六51.3.1圆管起动流动量方程初边值条件是速度分解为定常和非定常两部分之和为Heaviside阶跃函数 3.复杂粘弹性流体的流动2024/1/27周六52.定常部分解的非定常部分应满足齐次方程与定常部分迭加在一起，得到3.1圆管起动流2024/1/27周六53.3.1圆管起动流2024/1/27周六54.问题：Scott-Blair模型的圆管起动流是否会有不同的特性？3.1圆管起动流2024/1/27周六55.分数元模型动量方程其中得到3.1圆管起动流2024/1/27周六56.用Heaviside运算微积和分数阶微积分得到其中称为Mittag-Leffler函数,是指数函数的推广,为指数函数3.1圆管起动流2024/1/27周六57.2024/1/27周六58.2024/1/27周六59.2024/1/27周六60.2024/1/27周六61.小小结lThe constitutive equations of spring-dashpot systems can be easily derived by operational methods.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 .2024/1/27周六62.讨论牛顿流体在静止时不能承受剪切力，为何分数元流体会出现类固体的特性？Kelvin模型Maxwell模型Oldroyd-B模型2024/1/27周六63.讨论2024/1/27周六64.3.2圆管振荡流牛顿流体3.复杂粘弹性流体的流动2024/1/27周六65.牛顿流体3.2圆管振荡流2024/1/27周六66.2024/1/27周六67.Fractional Maxwell modelMaxwell model when ab 1分数阶Maxwell流体的圆管振荡流3.2圆管振荡流2024/1/27周六68.Exact solution3.2圆管振荡流2024/1/27周六69.Maxwellfluid PTFE3.2圆管振荡流2024/1/27周六70.3.3复杂粘弹性流体的Couette流Tan,W.C.&Xu,M.Y.,PlanesurfacesuddenlysetinmotioninaviscoelasticfluidwithfractionalMaxwellmodel.ActaMech.Sinica(2002)18:342349W.Shaowei&X.Mingyu,ExactsolutiononunsteadyCouetteflowofgeneralizedMaxwellfluidwithfractionalderivative,ActaMechanica(2006)187:103112HaitaoQi&HuiJin,UnsteadyrotatingflowsofaviscoelasticfluidwiththefractionalMaxwellmodelbetweencoaxialcylinders,ActaMech.Sinica(2006)22:3013052024/1/27周六71.4.分数阶微积分在流体力学中的其它应用Sugimoto,N,PropagationofnonlinearacousticwavesinatunnelwithanarrayofHelmholtzresonators.JFM,1992,244:55-784.1在排列有一组Helmholtz共振器的隧道中，非线性声波的传播2024/1/27周六72.4.2油岩层中温度场的延伸Control Equation:fractional diffusion-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.分数阶微积分在流体力学中的其它应用2024/1/27周六73.2024/1/27周六74.ChenW.Aspeculativestudyof2/3-orderfractionalLaplacianmodelingofturbulence,Chaos,2006,16:023126陈文，复杂流体和复杂流动（湍流）：一个比较研究.中国力学大会2007，会议论文摘要集（下）P431它同时包含拉氏算子描述的小尺度分子粘性和2/3阶拉氏算子表达的大尺度涡粘性，并清晰地将雷诺方程和标度率联系起来4.分数阶微积分在流体力学中的其它应用4.3湍流2024/1/27周六75.5.结束语Although the potential of fractional constitutive equations for viscoelastic fluids remains to be fully evaluated,it is our belief that this will be a primary direction for future developments in the quest for a predictive branch of fluid mechanics of some complex fluids.2024/1/27周六76.分数阶导数及其应用的研究是一项具有挑战性的工作分数阶微积分类似的应用在其它科学和技术领域不断被发现，包括电子技术、信号分析、控制理论、分形与混沌、热力学、生物物理和流变学等。问题：常数 c 的 阶导数是否应该等于零?5.结束语2024/1/27周六77.TheCaputosfractionalderivativesisdefinedasComparewithRiemann-Liouvilleoperatoroffractionalderivative5.结束语2024/1/27周六78.本讲座结束！2024/1/27周六79.