复旦学报(自然科学版) VoL. 52 No 4 2013年8月 Journal of Fudan University(Natural Science) Aug.2013 文章编号:0427-7104(2013)04-0547-11 有限变形理论的若干进展及其在 流体力学中的相关应用 谢锡麟,陈瑜,史倩 (复旦大学力学与工程科学系,上海200433) 摘要:概要性地叙述了作者新近提出的“当前物理构型对应之曲线坐标系显含时间的有限变形理论”、“几何 形态为曲面的连续介质的有限变形理论”,前者针对介质几何形态为 Euclid流形(体积形态),后者针对 Riemann 流形(曲面形态).类比于一般有限变形理论,上述理论均包括物理及参数构型构造·变形梯度定义及其基本性 质,变形刻画,输运定理以及守恒律方程。基于上述理论提出对应曲线坐标系显含时间的流函数-涡量解法,固定 曲面上二维不可压缩流动的流函数-涡量解法以及海面油污扩散控制方程,并给出了相关数值研究结果 关键词:有限变形理论;流函数-涡量解法;曲线坐标系显含时间;固定曲面上二维流动;可变形边界钝体绕 流;海面油污扩散 中图分类号:O33;035 文献标志码:A 1当前物理构型对应之曲线坐标系显含时间的有限变形理论 1.1研究背景 众所周知,当代飞机最为基础的飞行原理是按动量定理通过发动机提供向前推力,从而按 Bernoulli 定理由机翼产生升力——利用了流体力学中最为基础的规律.然而,自然界中各种各样的鱼儿游动与鸟儿 飞翔却没有“自身携带的发动机”,它们的自由运动往往仅依赖于躯体的摆动或者翅膀的挥舞——利用边 界的有限变形运动同周围流场间的相互作用而获得所需的各种动力.然而,至今我们对此方面的机制仍知 之甚少,我们当前所发展的航空和航海器的效率也远未及鱼儿与鸟儿 随着现代航空航海业的深度发展,国内外流体力学界开始广泛关注边界的有限变形运动同周围流场 间的相互作用.WuCJ等(2009,2010)3,主要基于CFD( Computational Fluid Dynamics,计算流体动力 学)以及控制理论等硏究3D三维)仿生鱼及鱼群的游动规律,发现类金枪鱼躯体上的旋涡脱落可同尾鳍 摆动造成的同向旋涡发生归并从而提供鱼的推进力;基于其所发展的实时控制理论,通过数值研究发现经 实时优化的二维仿生鱼的形态可显著提高鱼的推进力,流场反映为反向 Karman涡街的增强. Dong g j 等(2007)[研究了并列二维仿生鱼做同相位及反相位流向行波型运动对推进力等水动力学行为的影响, 发现了涡对涡街、单涡街和同相位及反相位同步涡街等4种大尺度旋涡结构.类似结构也在经典并列双圆 柱绕流以及横向振动单圆柱绕流中发现[.DuG等(2008)-0三维研究发现昆虫翅膀的可变形拍动(弯 拱及扭转,真实情况如此)较之刚性翅膀拍动可显著提高各项气动性能.LuXY等(2005)研究了壁面 可作流向行波运动的二维槽道流动,以模拟鱼游动时背脊的波动,发现优化的行波运动可以有效抑制分 离、提高推进及减小输入.WuCJ等(2003)-1研究二维机翼上表面的流向行波运动对流动的影响,基于 旋涡的空间特征提出“流动滚柱轴承”( fluid roller bearing)的概念,发现适当的表面行波运动可有效抑制 分离,显著提高机翼的空气动力性能.WuCJ等(2007)1进一步研究了二维圆柱后半部表面引入行波运 收稿日期:2013-05-08 基金项目:国家自然科学基金面上项目(1112069)资助;上海市教委2011年上海高等本科重点教学改革项目“‘现代 连续介质力学理论及实践’课程体系” 作者简介:谢锡麟(1974-),男,副教授,博士,E-mail:xiexilin@fudan.edu.cn. C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
文章编号：０４２７－７１０４（２０１３）０４－０５４７－１１ 收稿日期：２０１３－０５－０８ 基金项目：国家自然科学基金面上项目（１１１７２０６９）资助；上海市教委２０１１年上海高等本科重点教学改革项目“‘现代 连续介质力学理论及实践’课程体系” 作者简介：谢锡麟（１９７４—），男，副教授，博士，Ｅ－ｍａｉｌ：ｘｉｅｘｉｌｉｎ＠ｆｕｄａｎ．ｅｄｕ．ｃｎ． 有限变形理论的若干进展及其在 流体力学中的相关应用 谢锡麟，陈 瑜，史 倩 （复旦大学 力学与工程科学系，上海 ２００４３３） 摘 要：概要性地叙述了作者新近提出的“当前物理构型对应之曲线坐标系显含时间的有限变形理论”、“几 何 形态为曲面的连续介质的有限变形理论”，前者针对介质几何形态为 Ｅｕｃｌｉｄ流形（体积形态），后者针对 Ｒｉｅｍａｎｎ 流形（曲面形态）．类比于一般有限变形理论，上述理论均包括物理及参数构型构造，变形梯度定义及其基本性 质，变形刻画，输运定理以及守恒律方程．基于上述理论提出对应曲线坐标系显含时间的流函数 涡量解法，固定 曲面上二维不可压缩流动的流函数 涡量解法以及海面油污扩散控制方程，并给出了相关数值研究结果． 关键词：有限变形理论；流函数 涡量 解 法；曲线坐标系显含时间；固 定 曲 面 上 二 维 流 动；可变形边界钝体绕 流；海面油污扩散 中图分类号：Ｏ３３；Ｏ３５ 文献标志码：Ａ １ 当前物理构型对应之曲线坐标系显含时间的有限变形理论 １．１ 研究背景 众所周知，当代飞机最为基础的飞行原理是按动量定理通过发动机提供向前推力，从而按 Ｂｅｒｎｏｕｌｌｉ 定理由机翼产生升力———利用了流体力学中最为基础的规律．然而，自然界中各种各样的鱼儿游动与鸟儿 飞翔却没有“自身携带的发动机”，它们的自由运动往往仅依赖于躯体的摆动或者翅膀的挥舞———利用边 界的有限变形运动同周围流场间的相互作用而获得所需的各种动力．然而，至今我们对此方面的机制仍知 之甚少，我们当前所发展的航空和航海器的效率也远未及鱼儿与鸟儿［１－４］ ． 随着现代航空航海业的深度发展，国内外流体力学界开始广泛关注边界的有限变形运动同周围流场 间的相互作用．ＷｕＣＪ等（２００９，２０１０）［３，５，６］主要基于 ＣＦＤ（ＣｏｍｐｕｔａｔｉｏｎａｌＦｌｕｉｄＤｙｎａｍｉｃｓ，计算流体动力 学）以及控制理论等研究３Ｄ（三维）仿生鱼及鱼群的游动规律，发现类金枪鱼躯体上的旋涡脱落可同尾鳍 摆动造成的同向旋涡发生归并从而提供鱼的推进力；基于其所发展的实时控制理论，通过数值研究发现经 实时优化的二维仿生鱼的形态可显著提高鱼的推进力，流场反映为反向 Ｋａｒｍａｎ涡街的增强．ＤｏｎｇＧＪ 等（２００７）［７］研究了并列二维仿生鱼做同相位及反相位流向行波型运动对推进力等水动力学行为的影响， 发现了涡对涡街、单涡街和同相位及反相位同步涡街等４种大尺度旋涡结构．类似结构也在经典并列双圆 柱绕流［８］以及横向振动单圆柱绕流中发现［９］ ．ＤｕＧ 等（２００８）［１０］三维研究发现昆虫翅膀的可变形拍动（弯 拱及扭转，真实情况如此）较之刚性翅膀拍动可显著提高各项气动性能．ＬｕＸ Ｙ 等（２００５）［１１］研究了壁面 可作流向行波运动的二维槽道流动，以模拟鱼游动时背脊的波动，发现优化的行波运动可以有效抑制分 离、提高推进及减小输入．ＷｕＣＪ等（２００３）［１２］研究二维机翼上表面的流向行波运动对流动的影响，基于 旋涡的空间特征提出“流动滚柱轴承”（ｆｌｕｉｄｒｏｌｌｅｒｂｅａｒｉｎｇ）的概念，发现适当的表面行波运动可有效抑制 分离，显著提高机翼的空气动力性能．ＷｕＣＪ等（２００７）［１３］进一步研究了二维圆柱后半部表面引入行波运 第５２卷 第４期 ２０１３年８月 复 旦 学 报 （自然科学版） ＪｏｕｒｎａｌｏｆＦｕｄａｎＵｎｉｖｅｒｓｉｔｙ（ＮａｔｕｒａｌＳｃｉｅｎｃｅ） Ｖｏｌ．５２Ｎｏ．４ Ａｕｇ．２０１３
548 复旦学报(自然科学版) 第52卷 动的作用,发现适当的壁面行波运动可抑制 Karman涡街的形成,大幅度降低阻力 综上所述,绕流体边界的形态及其有限变形运动能够显著地改变流场中的主导旋涡结构及其空间演化 特性,从而大幅度地改变绕流体的气动或水动性能,甚至边界变形运动同流场的相互作用可直接提供动力 等.对此方面机制的研究与掌握必然对现代航空航海业具有极其重要的意义,故引起业界的广泛关注 1.2映照及张量分析观点 如图1所示,当研究可变形机翼的绕流时,我们可引入以下显含时间的曲线坐标系: X=(r(y,t)+·(R(y5)-r(y3,1)·cosn, X=5. X=(r(y,1)+·(R(y3)-r(n,t))·sin 此处(n,5)为刻画机翼表面(曲面)的双参数;r(n,3,t)和R(y,5)分别刻画机翼表面以及流场外界面,二者 之间所围区域可定义为所需数值模拟的物理流场.利用上述曲线坐标系,物理流场所对应的参数区域为相 对于时间固定的方块 当前物理构型 当前参数构型 (边界可变形) (边界始终固定)\ rn. 5, n) 图1适用于表面可做有限变形运动机翼绕流的显含时间的曲线坐标系示意图 Fig 1 Sketch of curvilinear coordinates including time explicitly that is suitable to flows around an airfoil with deformable boundaries 进一步,基于一般曲线坐标系的张量场分析,可将任意张量场及其场论微分运算在曲线坐标所诱导的 局部基下展开,如图2所示.如对流项的展开: V. vOV=V V, g =V(.D).(ar+r; Ve )(a D8 (, D) V=V (x, Dg (r, D) 其中r(x,1)为 Christoffel符号.籍此,我们可以获得各种控制方程在参数域上的分量控制方程 =(xng⑧g⑧g(x =d(xn)g②gg(x20 参数区域 x-曲线8( g(xn) x-曲线 gi(r,) x-曲线 x-曲线 物理区域 x-曲线 曲线坐标X=x,D 图2一般曲线坐标系及其所诱导的局部基示意图 Fig 2 Sketch of general curvilinear coordinates with its local induced basis C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
动的作用，发现适当的壁面行波运动可抑制 Ｋａｒｍａｎ涡街的形成，大幅度降低阻力． 综上所述，绕流体边界的形态及其有限变形运动能够显著地改变流场中的主导旋涡结构及其空间演化 特性，从而大幅度地改变绕流体的气动或水动性能，甚至边界变形运动同流场的相互作用可直接提供动力 等．对此方面机制的研究与掌握必然对现代航空航海业具有极其重要的意义，故引起业界的广泛关注． １．２ 映照及张量分析观点 如图１所示，当研究可变形机翼的绕流时，我们可引入以下显含时间的曲线坐标系： Ｘ１＝（ｒ（η，ζ，ｔ）＋ξ·（Ｒ（η，ζ）－ｒ（η，ζ，ｔ）））·ｃｏｓη， Ｘ２＝ζ， Ｘ３＝（ｒ（η，ζ，ｔ）＋ξ·（Ｒ（η，ζ）－ｒ（η，ζ，ｔ）））·ｓｉｎη 烅 烄 烆 ． 此处（η，ζ）为刻画机翼表面（曲面）的双参数；ｒ（η，ζ，ｔ）和Ｒ（η，ζ）分别刻画机翼表面以及流场外界面，二者 之间所围区域可定义为所需数值模拟的物理流场．利用上述曲线坐标系，物理流场所对应的参数区域为相 对于时间固定的方块． 图１ 适用于表面可做有限变形运动机翼绕流的显含时间的曲线坐标系示意图 Ｆｉｇ．１ Ｓｋｅｔｃｈｏｆｃｕｒｖｉｌｉｎｅａｒｃｏｏｒｄｉｎａｔｅｓｉｎｃｌｕｄｉｎｇｔｉｍｅｅｘｐｌｉｃｉｔｌｙｔｈａｔｉｓｓｕｉｔａｂｌｅｔｏ ｆｌｏｗｓａｒｏｕｎｄａｎａｉｒｆｏｉｌｗｉｔｈｄｅｆｏｒｍａｂｌｅｂｏｕｎｄａｒｉｅｓ 进一步，基于一般曲线坐标系的张量场分析，可将任意张量场及其场论微分运算在曲线坐标所诱导的 局部基下展开，如图２所示．如对流项的展开： Ｖ珝· Δ Ｖ珝＝Ｖｊ Δ ｊＶｉ ｇ珝ｉ＝Ｖｊ（ｘ，ｔ）· Ｖｉ ｘｊ＋Γｉ （ ） ｊｋＶｋ （ｘ，ｔ）ｇ珝ｉ（ｘ，ｔ），Ｖ＝Ｖｉ（ｘ，ｔ）ｇ珝ｉ（ｘ，ｔ）， 其中Γｉ ｊｋ（ｘ，ｔ）为 Ｃｈｒｉｓｔｏｆｆｅｌ符号．籍此，我们可以获得各种控制方程在参数域上的分量控制方程． 图２ 一般曲线坐标系及其所诱导的局部基示意图 Ｆｉｇ．２ Ｓｋｅｔｃｈｏｆｇｅｎｅｒａｌｃｕｒｖｉｌｉｎｅａｒｃｏｏｒｄｉｎａｔｅｓｗｉｔｈｉｔｓｌｏｃａｌｉｎｄｕｃｅｄｂａｓｉｓ ８４５ 复 旦 学 报（自然科学版） 第５２卷
第4期 谢锡麟等:有限变形理论的若干进展及其在流体力学中的相关应用 549 按上述映照及张量分析观点,所得控制方程的分量形式有如下特点:①分量形式为建立在参数区域 上的偏微分方程( Partial Differential Equations,PDE),计算域形态随时间始终保持方体而边界条件可随 时间变化,由此大大简化了差分格式等构造上的复杂性;②现分量方程较之直角坐标系的分量方程,可能 引入混合偏导数项、度量张量分量、 Christoffel符号等而显得更为复杂,但并不引入更高阶导数项,需指出 方程中显含的度量张量分量、 Christoffel符号等直接由曲线坐标系确定,事先就可获得解析表达式或通过 数值方法确定;③在壁面上,物理量按相对于反映曲面几何性质的局部基展开,往往便于边界条件的 处理 国际上较早基于张量分析获得一般曲线坐标系下的流动控制方程(物理量基于局部基展开)并进 行数值研究的理论成果包括国际著名透平机械专家吴仲华先生首创的关于S1及S2流面的研究;李 开泰、黄艾香在张量分析及其应用方面的较为系统的研究[15,他们的研究主要应用于旋转机械叶片附 近的流动.这些研究一般仅涉及静止或刚性的曲面壁面而未涉及壁面的有限变形运动,而我们将当前 物理构型对应的曲线坐标系(微分同胚)推广到显含时间情形,籍此可获得既“规则”又“固定”的参数区 域并基于张量分析获得其上的控制方程.谢锡麟(2012)16叙述了现代张量分析在连续介质力学中的 若干应用 1.3有限变形理论 我们按照郭仲衡所述的有限变形理论[,平行发展了当前物理构型对应之曲线坐标系显含时间的有 限变形理论[.有限变形理论的基本内容,可以依次归结为4部分:(1)物理构型及参数构型构造; (2)变形梯度及其基本性质;(3)变形刻画及输运定理(包括:①初始物理构型以及当前物理构型中有向 线元、面元以及体元之间的关系;②初始物理构型以及当前物理构型中有向线元、面元以及体元模之间的 关系;③当前物理构型中有向线元、面元以及体元的物质导数同其自身之间的关系;④当前物理构型中 有向线元、面元以及体元模的物质导数同其自身之间的关系),基于关系式③、④易于获得最一般的第二类 以及第一类线、面、体的输运定理;(4)守恒律控制方程 需指出,当前物理构型对应之曲线坐标系显含时间的有限变形理论同一般理论的主要差别在于速度 以及物质导数的表达式: △(,D)=X(x,)x X (,t)g:(x,t)+(x,t) a x (x,1)=9x 相对于一般情形,速度表达式出现增加项:a(x,D) ),(x,)=()(x,)(x,),完全 源于当前物理构型对应之曲线坐标系显含时间.进一步,对任意张量场的物质导数,有 9()D(x,)+8(,)(,)=(x,)+·-(x,D)·V②虫 相对于一般形式,增加了项一(x,1)·8,8小=()②单会②(x,1代表相对于 Euclid坐 标的全梯度算子对于变形梯度FA∞x(,1)g,(x,n)⑧G(),其基本意义及基本性质保持不变,即有: 正=((x,D)·F=:1L·F,adtF=tF,此处 det F4vgdet(en),0:=.(x,).由 d 此,变形刻画所有关系式同一般理论;当保留物质导数的整体形式,输运定理所有关系式同一般理论.需指 出,现研究的连续介质的几何形态仍为 Euclid流形,故在张量场分析上并无差异 1.4应用事例 近期我们基于当前物理构型对应之曲线坐标系显含时间的有限变形理论发展了对应显含时间曲线坐 标系的流函数-涡量解法[81.对已有报道的通过在圆柱边界上引入行波可消除涡街的研究[1,我们利用 现有方法实现了相关结果,并将相关效果推广至椭圆柱绕流,如图3(见第550页)所示 C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
按上述映照及张量分析观点，所得控制方程的分量形式有如下特点：① 分量形式为建立在参数区域 上的偏微分方程（ＰａｒｔｉａｌＤｉｆｆｅｒｅｎｔｉａｌＥｑｕａｔｉｏｎｓ，ＰＤＥ），计算域形态随时间始终保持方体而边界条件可随 时间变化，由此大大简化了差分格式等构造上的复杂性；② 现分量方程较之直角坐标系的分量方程，可能 引入混合偏导数项、度量张量分量、Ｃｈｒｉｓｔｏｆｆｅｌ符号等而显得更为复杂，但并不引入更高阶导数项，需指出 方程中显含的度量张量分量、Ｃｈｒｉｓｔｏｆｆｅｌ符号等直接由曲线坐标系确定，事先就可获得解析表达式或通过 数值方法确定；③ 在 壁 面 上，物理量按相对于反映曲面几何性质的局部基展开，往 往 便 于 边 界 条 件 的 处理． 国际上较早基于张量分析获得一 般 曲 线 坐 标 系 下 的 流 动 控 制 方 程（物 理 量 基 于 局 部 基 展 开）并 进 行数值研究的理论成果包括国际著名透平机械专家吴仲华先生首创的关于Ｓ１及Ｓ２流面的研究［１４］；李 开泰、黄艾香在张量分析及其应用方面的较为系统 的 研 究［１５］，他们的研究主要应用于旋转机械叶片附 近的流动．这些研究一般仅涉及静止或刚性的曲面壁面而未涉及壁面的有限变形运动，而 我 们 将 当 前 物理构型对应的曲线坐标系（微分同胚）推广到显含时间情形，籍此可获得既“规 则”又“固 定”的 参 数 区 域并基于张量分析获得其上的控制方程．谢 锡 麟（２０１２）［１６］叙 述 了 现 代 张 量 分 析 在 连 续 介 质 力 学 中 的 若干应用． １．３ 有限变形理论 我们按照郭仲衡所述的有限变形理论［１７］，平行发展了当前物理构型对应之曲线坐标系显含时间的有 限变形 理 论［１８］ ．有 限 变 形 理 论 的 基 本 内 容，可 以 依 次 归 结 为 ４ 部 分：（１）物 理 构 型 及 参 数 构 型 构 造； （２）变形梯度及其基本性质；（３）变形刻画及输运定理（包括：① 初始物理构型以及当前物理构型中有向 线元、面元以及体元之间的关系；② 初始物理构型以及当前物理构型中有向线元、面元以及体元模之间的 关系；③ 当前物理构型中有向线元、面元以及体元的物质导数同其自身之间的关系；④ 当前物理构型中 有向线元、面元以及体元模的物质导数同其自身之间的关系），基于关系式③、④易于获得最一般的第二类 以及第一类线、面、体的输运定理；（４）守恒律控制方程． 需指出，当前物理构型对应之曲线坐标系显含时间的有限变形理论同一般理论的主要差别在于速度 以及物质导数的表达式： Ｖ珝 Ｘ珤 ｔ （ξ，ｔ）＝ Ｘ珤 ｘｉ（ｘ，ｔ）ｘｉ ｔ （ξ，ｔ）＋Ｘ珤 ｔ （ｘ，ｔ）＝ ｘｉ ｔ （ξ，ｔ）ｇ珝ｉ（ｘ，ｔ）＋Ｘ珤 ｔ （ｘ，ｔ）． 相对于一般情形，速度表达式出现增加项：Ｘ珤 ｔ （ｘ，ｔ）＝ Ｘ珤 ｔ（ ） ，ｇ珝ｉ 瓗３ｇ珝ｉ（ｘ，ｔ）＝ Ｘ珤 （ ） ｔ ｉ （ｘ，ｔ）ｇ珝ｉ（ｘ，ｔ），完全 源于当前物理构型对应之曲线坐标系显含时间．进一步，对任意张量场的物质导数，有： Φ Φ ｔ （ξ，ｔ）＝ Φ ｔ （ｘ，ｔ）＋Φ ｘｉ（ｘ，ｔ）ｘｉ ｔ （ξ，ｔ）＝ Φ ｔ （ｘ，ｔ）＋Ｖ珝· Δ Φ－Ｘ珤 ｔ （ｘ，ｔ）· Δ Φ． 相对于一般形式，增加了项－Ｘ ｔ （ｘ，ｔ）· Δ Φ， Δ Φ＝ ｇ珗ｌ （ ） ｘｌ Φｇ珗ｌ Φ ｘｌ（ｘ，ｔ）代表Φ 相对于 Ｅｕｃｌｉｄ坐 标的全梯度算子．对于变形梯度 Ｆｘｉ ξ Ａ （ξ，ｔ）ｇｉ（ｘ，ｔ）ＧＡ （ξ），其基本意义及基本性质保持不变，即 有： ｄＦ ｄｔ＝ Ｖ珝 Δ （ ） （ｘ，ｔ） ·Ｆ＝∶Ｌ·Ｆ，ｄ ｄｔ ｄｅｔＦ＝θｄｅｔＦ，此处ｄｅｔＦ 槡ｇ 槡Ｇ ｄｅｔ ｘｉ ξ（ ） Ａ（ξ，ｔ） ，θ∶＝Ｖ珝· Δ （ｘ，ｔ）．由 此，变形刻画所有关系式同一般理论；当保留物质导数的整体形式，输运定理所有关系式同一般理论．需指 出，现研究的连续介质的几何形态仍为 Ｅｕｃｌｉｄ流形，故在张量场分析上并无差异． １．４ 应用事例 近期我们基于当前物理构型对应之曲线坐标系显含时间的有限变形理论发展了对应显含时间曲线坐 标系的流函数 涡量解法［１８，１９］ ．对已有报道的通过在圆柱边界上引入行波可消除涡街的研究［１３］，我们利用 现有方法实现了相关结果，并将相关效果推广至椭圆柱绕流，如图３（见第５５０页）所示． 第４期 谢锡麟等：有限变形理论的若干进展及其在流体力学中的相关应用 ９４５
复旦学报(自然科学版) 第52卷 0 P150 r150 6420 图3圆柱以及椭圆柱(长短轴比为1.5)表面引入行波以消除涡街发展:(左图)涡量分布, (右图)流函数分布.Re=400 Fig 3 Suppression of vortex street by traveling wave generated on the surfaces of circular and elliptical cylinders (the ratio of the long and short axes is 1. 5):(left subplot) vorticity distribution, (right subplot)stream function distribution, Re=400 2几何形态为曲面的连续介质的有限变形理论 2.1研究背景 考虑星体表面的大气运动,海面上油污扩散以及洪水蔓延过平原、洼地以及山丘等,流动的法向尺度 (流层厚度)远远小于流动的展向(流向)尺度.由此,可将此类流动模型化为几何形态为曲面的连续介质的 流动.由于三维空间中的曲面为双参数向量值映照,故上述流动又可称为二维流动( two dimension flows).目前,业内对二维流动的研究多限于平面形态的皂膜流动{,然而真正的代表性的二维流动应为 几何形态为一般曲面的流动,其典型形式可为固定曲面上的二维流动,或自身运动曲面上的二维流动(如 细胞膜等流体膜上的流动) Aris(1964)[21研究了曲面上的流动,包括运动刻画,质量守恒以及动量守恒等微分学控制方程,但其 就质量守恒等的分析结果有误[23,且关于运动刻画并未具有一般有限变形理论的系统性论述 2.2有限变形理论 几何形态为曲面的有限变形理论,其运动学与几何学仍可以参照一般有限变形理论,定义初始、当前 物理构型及其对应的初始、当前曲线坐标系,此处曲线坐标系在曲面的参数域中定义,实际也为曲面作为 流形的局部参数化(对应微分流形中的坐标卡);按微分学概念引入变形梯度并获得其基本性质;基于变形 梯度获得变形的4类刻画;基于变形刻画获得所有形式的输运方程[.动力学方面,主要基于我们发展的 内蕴形式第二类广义 Stokes公式l2)获得质量、动量、动量矩以及能量守恒等守恒律微分方程.内蕴形式 第二类广义 Stokes公式提供了将沿切平面并正交于边界曲线切向量方向作用的张量场转化至边界曲线 所围的曲面上积分的一般方法,如图4所示 (×)°-重=(。-更+H(-)d 此处φ为定义在曲面上的任意张量场,H:=b=gb为平均曲率,g=(g1,g)2(x)为曲面第一基本 量,4=(82(x3),n),为曲面第二基本量,(x2)为曲面自身诱导的局部协变基;一代表任意合法的张 量代数运算 C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
图３ 圆柱以及椭圆柱（长短轴比为１．５）表面引入行波以消除涡街发展：（左图）涡量分布， （右图）流函数分布．Ｒｅ＝４００． Ｆｉｇ．３ Ｓｕｐｐｒｅｓｓｉｏｎｏｆｖｏｒｔｅｘｓｔｒｅｅｔｂｙｔｒａｖｅｌｉｎｇｗａｖｅｇｅｎｅｒａｔｅｄｏｎｔｈｅｓｕｒｆａｃｅｓｏｆｃｉｒｃｕｌａｒａｎｄｅｌｌｉｐｔｉｃａｌｃｙｌｉｎｄｅｒｓ （ｔｈｅｒａｔｉｏｏｆｔｈｅｌｏｎｇａｎｄｓｈｏｒｔａｘｅｓｉｓ１．５）：（ｌｅｆｔｓｕｂｐｌｏｔ）ｖｏｒｔｉｃｉｔｙｄｉｓｔｒｉｂｕｔｉｏｎ，（ｒｉｇｈｔｓｕｂｐｌｏｔ）ｓｔｒｅａｍ ｆｕｎｃｔｉｏｎｄｉｓｔｒｉｂｕｔｉｏｎ．Ｒｅ＝４００． ２ 几何形态为曲面的连续介质的有限变形理论 ２．１ 研究背景 考虑星体表面的大气运动，海面上油污扩散以及洪水蔓延过平原、洼地以及山丘等，流动的法向尺度 （流层厚度）远远小于流动的展向（流向）尺度．由此，可将此类流动模型化为几何形态为曲面的连续介质的 流动．由于三维空间中的曲面为双参数向量值映照，故 上 述 流 动 又 可 称 为 二 维 流 动 （ｔｗｏｄｉｍｅｎｓｉｏｎａｌ ｆｌｏｗｓ）．目前，业内对二维流动的研究多限于平面形态的皂膜流动［２０］，然而真正的代表性的二维流动应为 几何形态为一般曲面的流动，其典型形式可为固定曲面上的二维流动，或自身运动曲面上的二维流动（如 细胞膜等流体膜上的流动）． Ａｒｉｓ（１９６４）［２１］研究了曲面上的流动，包括运动刻画，质量守恒以及动量守恒等微分学控制方程，但其 就质量守恒等的分析结果有误［２２］，且关于运动刻画并未具有一般有限变形理论的系统性论述． ２．２ 有限变形理论 几何形态为曲面的有限变形理论，其运动学与几何学仍可以参照一般有限变形理论，定义初始、当前 物理构型及其对应的初始、当前曲线坐标系，此处曲线坐标系在曲面的参数域中定义，实际也为曲面作为 流形的局部参数化（对应微分流形中的坐标卡）；按微分学概念引入变形梯度并获得其基本性质；基于变形 梯度获得变形的４类刻画；基于变形刻画获得所有形式的输运方程［１８］ ．动力学方面，主要基于我们发展的 内蕴形式第二类广义Ｓｔｏｋｅｓ公式［１８，２２］获得质量、动量、动量矩以及能量守恒等守恒律微分方程．内蕴形式 第二类广义Ｓｔｏｋｅｓ公式提供了将沿切平面并正交于边界曲线切向量方向作用的张量场转化至边界曲线 所围的曲面上积分的一般方法，如图４所示． ∮Σ ｔ （ ） τ珒×ｎ珗°－Φ ＝∫Σ ｔ （Σ Δ °－Φ＋Ｈ（ｎ珗°－Φ））ｄσ． 此处Φ 为定义在曲面上的任意张量场，Ｈ∶＝ｂｓ ｓ＝ｇｓｔｂｔｓ为平均曲率，ｇｉｊ＝（ｇ珝ｉ，ｇ珝ｊ）瓗３ （ｘΣ）为曲面第一基本 量，ｂｉｊ＝ 珗ｇｊ ｘｉ Σ （ ） （ｘΣ），ｎ 瓗３ 为曲面第二基本量，ｇ珝ｊ（ｘΣ）为曲面自身诱导的局部协变基；°－代表任意合法的张 量代数运算． ０５５ 复 旦 学 报（自然科学版） 第５２卷
第4期 谢锡麟等:有限变形理论的若干进展及其在流体力学中的相关应用 曲面方程:E(x:D 图4内蕴形式的第二类广义 Stokes公式示意图 Fig 4 Sketch of the intrinsic generalized Stokes formula of the second kind 考虑质量守恒、动量守恒、动量矩守恒以及能量守恒在几何形态为曲面的连续介质的应用,一般需要 考虑的各种作用表现为介质边界上的积分,且作用方向为切平面中的×n方向.对各种形式沿×n作用 的边界上曲线积分可以利用内蕴形式的第二类广义 Stokes公式转化为面积分.结合运动学中的各种输运 方程,我们就可获得各种守恒律控制方程的微分形式 设想切平面内的应力t=t,g⑧∈T(T),本文称为曲面应力.利用内蕴形式的第二类广义 Stokes公式,可几近平凡地获得其整体作用: (×).t=(.t+H(t)d=「 式中÷,=71+的,,=表示曲面上的梯度算子;表示曲面(作为 Riemann流形)上的协 变导数.结合第一类面输运定理以及质量守恒,则可得动量守恒的微分方程: pan=(+P)x;+(,+f)示, 此处子:=f十为面力分布,其整体作用形式直接表示为面积分fd;面力分布可以为重力、电磁 力等.Aris研究了曲面上的流动21,获得了上述动量方程(未考虑面力作用).本文基于内蕴形式的第二类 广义 Stokes公式获得此方程,相对而言显得更为直接且简便 需指出,我们可按实际研究情形考虑各种形式的作用,对于沿×n方向的作用可通过内蕴形式的第 二类广义 Stokes公式将边界上的曲线积分转化为面积分,结合输运方程就可获得守恒律控制方程 2.3应用事例 2.3.1固定曲面上不可压缩流动的流函数-涡量解法 按变形分析曲面上流动的祸量定义为,一CV1一示,此处剧wm量-云 vg=(g1,g2,n).一般连续性方程为:p+p=0,此处p为密度,:=V vax(vgV)(x2,1).不可 压缩流动的连续性方程为:0=0.由此,可引入流函数:V=c"(x,1).曲面应力可以考虑为: =(-p+7)1+(V+5V,)⑧∈T(T),P为沿切平面作用的内压力,y为表面张力,为流 体自身之间的黏性系数,I=g⑧g;此种形式的曲面应力自然满足无面力偶分布的动量矩守恒2.现 流函数-涡量控制方程具有如下形式: C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
图４ 内蕴形式的第二类广义Ｓｔｏｋｅｓ公式示意图 Ｆｉｇ．４ ＳｋｅｔｃｈｏｆｔｈｅｉｎｔｒｉｎｓｉｃｇｅｎｅｒａｌｉｚｅｄＳｔｏｋｅｓｆｏｒｍｕｌａｏｆｔｈｅｓｅｃｏｎｄｋｉｎｄ 考虑质量守恒、动量守恒、动量矩守恒以及能量守恒在几何形态为曲面的连续介质的应用，一般需要 考虑的各种作用表现为介质边界上的积分，且作用方向为切平面中的τ珒×ｎ珗方向．对各种形式沿τ珒×ｎ珗作用 的边界上曲线积分可以利用内蕴形式的第二类广义Ｓｔｏｋｅｓ公式转化为面积分．结合运动学中的各种输运 方程，我们就可获得各种守恒律控制方程的微分形式． 设想切平面内的 应 力ｔ＝ｔｉ ·ｊｇ珝ｉ ｇ珝ｊ ∈Ｔ２ （ＴΣ），本文称为曲面 应 力．利用内蕴形式的第二类 广 义 Ｓｔｏｋｅｓ公式，可几近平凡地获得其整体作用： ∮Σ ｔ （τ珒×ｎ珗）·ｔ＝∫Σ ｔ （Σ Δ ·ｔ＋Ｈ（ｎ珗·ｔ））ｄσ＝∫Σ ｔ Σ Δ ·ｔｄσ， 式中Σ Δ ·ｔ＝ Δ ｉ Σ ｔｉ ·ｊｇ珝ｊ ＋ｂｊ ｉｔｉ ·ｊｎ珗， Σ Δ ＝ｇ珝ｌ ｘｌ Σ 表示曲面上的梯度算子； Δ ｉ Σ 表示曲面（作为 Ｒｉｅｍａｎｎ流形）上的协 变导数．结合第一类面输运定理以及质量守恒，则可得动量守恒的微分方程： ρ ｄＶ珝 ｄｔ ＝ （Σ Δ ｉｔｉｊ ＋ｆ ）ｊ ｇ珝ｊ ＋ （ｂｊ ｉｔｉ ·ｊ ＋ｆ ）３ ｎ珗， 此处ｆ → ∶＝ｆｊ ｇ珝ｊ＋ｆ３ ｎ珗为面力分布，其整体作用形式直接表示为面积分∫Σ ｔ ｆ → ｄσ；面力分布可以为重力、电磁 力等．Ａｒｉｓ研究了曲面上的流动［２１］，获得了上述动量方程（未考虑面力作用）．本文基于内蕴形式的第二类 广义Ｓｔｏｋｅｓ公式获得此方程，相对而言显得更为直接且简便． 需指出，我们可按实际研究情形考虑各种形式的作用，对于沿τ珒×ｎ珗方向的作用可通过内蕴形式的第 二类广义Ｓｔｏｋｅｓ公式将边界上的曲线积分转化为面积分，结合输运方程就可获得守恒律控制方程． ２．３ 应用事例 ２．３．１ 固定曲面上不可压缩流动的流函数 涡量解法 按变形分析，曲面上流动的涡量定义为：ω珗 ＝ε ｓｔ３ Σ Δ ｓＶｔｎ珗∶＝ω３ ｎ珗，此 处 Ｅｄｄｉｎｇｔｏｎ张 量ε ｓｔ３ ＝ ｅｓｔ３ 槡ｇ ， 槡ｇ ＝（ｇ珝１，ｇ珝２，ｎ珗）．一般连续性方程为：ρ＋ρθ＝０，此处ρ为密度，θ∶＝ Δ ｉ Σ Ｖｉ ＝ １ 槡ｇ ｘｉ （ ） 槡ｇＶｉ （ｘΣ，ｔ）．不可 压缩流动 的 连 续 性 方 程 为：θ ＝ ０．由 此，可 引 入 流 函 数：Ｖｓ ＝ε ｓｔ３ ψ ｘｔ（ｘ，ｔ）．曲面应力可以考虑为： ｔ＝ －（ ） ｐ＋γＩ＋μ Σ Δ ｊＶｉ ＋ Σ Δ （ ） ｉＶｊ ｇ珝ｉ ｇ珝ｊ ∈Ｔ２ （ ） ＴΣ ，ｐ为沿切平面作用的内压力，γ为表面张力，μ为流 体自身之间的黏性系数，Ｉ＝ｇｉｊｇ珝ｉ ｇ珝ｊ；此种形式的曲面应力自然满足无面力偶分布的动量矩守恒［２３］ ．现 流函数 涡量控制方程具有如下形式： 第４期 谢锡麟等：有限变形理论的若干进展及其在流体力学中的相关应用 １５５
552 复旦学报(自然科学版) 第52卷 arax at (I,t)+vs a3 ar(x,)=(D3+2-5(KV1)+15fm 相对于平面不可压缩流动情形,涡量控制方程中增加了 Gauss曲率同速度的耦合作用项 图5显示了固定曲面上不可压缩流体绕过圆柱的流场,此处适当考虑了流体同承载曲面之间的摩 擦[∞.对现情形,我们仍观察到涡街;但相对一般二维情形,由于摩擦作用使得旋涡结构有所“模糊”,且曲 面的弯曲亦对旋涡结构有所“扭曲 r500 图5固定曲面上圆柱绕流流动:(左图)流函数分布三维视图,(右图)涡量分布平面投影.Re=500 Fig. 5 The flow around a circular cylinder on a fixed surface: (left subplot)three dimensional view of stream function distribution: (right subplot) planform of vorticity distribution. Re=500 图6显示沿凹凸螺旋面的不可压缩二维内流的涡量三维分布,在螺旋面的边界上满足黏性边界条件 (无滑移及无贯穿条件), Reynolds数为300(以螺旋面宽度作为特征尺度),并考虑了作为面力形式的摩擦 力.我们已从理论上获得,对任意固定曲面上的不可压缩黏性流动,若在固定边界上满足黏性边界条件(无 滑移及无贯穿,则在固定边界上变形率张量的特征值具有相同的绝对值2,其所对应的方向同固体边 界上切向量呈45°或者135t23.图7为内边界上变形率张量的特征问题分析,包括最大特征值所对应的特 征向量同内边界切向量夹角沿内边界的分布;涡量值与最大特征值之比沿内边界的分布.数值研究结果很 好地符合理论分析结果 图6(左图)固定凹凸螺旋曲面上的二维不可压缩内流涡量分布三维视图,内边界上的2个线段之间的区间 具有负的涡量值;(右图)固定凹凸螺旋曲面 Gauss曲率分布的投影视图.Re=300 Fig 6(left subplot) Three dimensional view of a two dimensional incompressible flow on a fixed undulated helicoidal surface: spatial distribution of the vorticity. Two line segments make the interval of the inner boundary with negative vorticity. (right subplot) Projective views of the spatial distributions of the Gaussian curvature of the undulated helicoidal surface, Re=300 C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
Δ Σ ψ ＝ｇｉｊ ２ ψ ｘｉ ｘｊ（ｘ，ｔ）－Γｋ ｉｊ ψ （ ） ｘｋ（ｘ，ｔ） ＝－ω３， ω３ ＝ ω３ ｔ（ ） ｘ，ｔ ＋Ｖｓ ω３ ｘｓ （ ） ｘ，ｔ ＝ μ ρ Σｓ ΔΣ Δ ｓω３ ＋２εｋｌ３ Σ Δ （ ） ｋ（ＫＧＶｌ）＋ １ ρ εｋｌ３ Σ Δ ｋｆｓｕｒ，ｌ 烅 烄 烆 ． 相对于平面不可压缩流动情形，涡量控制方程中增加了 Ｇａｕｓｓ曲率同速度的耦合作用项． 图５显示了固定曲面上不可压缩流体绕过圆柱的流场，此处适当考虑了流体同承载曲面之间的摩 擦［２０］ ．对现情形，我们仍观察到涡街；但相对一般二维情形，由于摩擦作用使得旋涡结构有所“模糊”，且曲 面的弯曲亦对旋涡结构有所“扭曲”． 图５ 固定曲面上圆柱绕流流动：（左图）流函数分布三维视图，（右图）涡量分布平面投影．Ｒｅ＝５００． Ｆｉｇ．５ Ｔｈｅｆｌｏｗａｒｏｕｎｄａｃｉｒｃｕｌａｒｃｙｌｉｎｄｅｒｏｎａｆｉｘｅｄｓｕｒｆａｃｅ：（ｌｅｆｔｓｕｂｐｌｏｔ）ｔｈｒｅｅｄｉｍｅｎｓｉｏｎａｌｖｉｅｗ ｏｆｓｔｒｅａｍｆｕｎｃｔｉｏｎｄｉｓｔｒｉｂｕｔｉｏｎ；（ｒｉｇｈｔｓｕｂｐｌｏｔ）ｐｌａｎｆｏｒｍｏｆｖｏｒｔｉｃｉｔｙｄｉｓｔｒｉｂｕｔｉｏｎ．Ｒｅ＝５００． 图６显示沿凹凸螺旋面的不可压缩二维内流的涡量三维分布，在螺旋面的边界上满足黏性边界条件 （无滑移及无贯穿条件），Ｒｅｙｎｏｌｄｓ数为３００（以螺旋面宽度作为特征尺度），并考虑了作为面力形式的摩擦 力．我们已从理论上获得，对任意固定曲面上的不可压缩黏性流动，若在固定边界上满足黏性边界条件（无 滑移及无贯穿），则在固定边界上变形率张量的特征值具有相同的绝对值 ω３ ２ ，其所对应的方向同固体边 界上切向量呈４５°或者１３５°［２３］ ．图７为内边界上变形率张量的特征问题分析，包括最大特征值所对应的特 征向量同内边界切向量夹角沿内边界的分布；涡量值与最大特征值之比沿内边界的分布．数值研究结果很 好地符合理论分析结果． 图６ （左图）固定凹凸螺旋曲面上的二维不可压缩内流涡量分布三维视图，内边界上的２个线段之间的区间 具有负的涡量值；（右图）固定凹凸螺旋曲面 Ｇａｕｓｓ曲率分布的投影视图．Ｒｅ＝３００． Ｆｉｇ．６ （ｌｅｆｔｓｕｂｐｌｏｔ）Ｔｈｒｅｅｄｉｍｅｎｓｉｏｎａｌｖｉｅｗｏｆａｔｗｏｄｉｍｅｎｓｉｏｎａｌｉｎｃｏｍｐｒｅｓｓｉｂｌｅｆｌｏｗｏｎａｆｉｘｅｄｕｎｄｕｌａｔｅｄ ｈｅｌｉｃｏｉｄａｌｓｕｒｆａｃｅ：ｓｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅｖｏｒｔｉｃｉｔｙ．Ｔｗｏｌｉｎｅｓｅｇｍｅｎｔｓｍａｋｅｔｈｅｉｎｔｅｒｖａｌｏｆｔｈｅ ｉｎｎｅｒｂｏｕｎｄａｒｙｗｉｔｈｎｅｇａｔｉｖｅｖｏｒｔｉｃｉｔｙ．（ｒｉｇｈｔｓｕｂｐｌｏｔ）Ｐｒｏｊｅｃｔｉｖｅｖｉｅｗｓｏｆｔｈｅｓｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｓｏｆ ｔｈｅＧａｕｓｓｉａｎｃｕｒｖａｔｕｒｅｏｆｔｈｅｕｎｄｕｌａｔｅｄｈｅｌｉｃｏｉｄａｌｓｕｒｆａｃｅ．Ｒｅ＝３００． ２５５ 复 旦 学 报（自然科学版） 第５２卷
第4期 谢锡麟等:有限变形理论的若干进展及其在流体力学中的相关应用 553 -2 60120180240300360 60120180240300360 图7内边界上变形率张量特征问题分析:(左图)涡量与最大特征值之比沿内边界的分布,注:理论预测值为2或-2. (右图)最大特征值所对应的特征向量同内边界切向量的夹角沿内边界的分布.注:理论预测值为45°或135° Fig 7 Eigenproblem analysis of the strain tensor on the inner boundary: (left subplot Distribution of the ratio of the orticity and the maximum eigenvalue of the strain tensor along the inner boundary of the undulated helicoidal surface. Remark: The theoretical indicated value of the ratio is 2 or -2. (right subplot Distribution of the angle between the eigenvector with respect to the maximum eigenvalue of the strain tensor and the tangent vector of the boundary along the inner boundary of the undulated helicoidal surface. Remark: The theoretical indicated value of the angle is 45 or 135 degree. 固定曲面上的二维流动可以对应星体表面上的大气运动,或者地球表面上的洪水蔓延过洼地、山丘或 者盘山公路等,这些流动的法向尺度(流层厚度)远远小于流动的流向尺度,故可将流动视为二维.近期,作 者类比于一般涡量动力学[2,发展了固定曲面上二维流动的涡量动力学理论框架[23,希望能为研究此类 流动铺垫了一定的理论基础 2海面上油污扩散 海面上油污扩散表现为几何形态为曲面的油污在海面上的运动过程,对应的物理构型以及参数构型 构造如图8所示,现运动刻画的参数表示为x=xx(x,t)∈R2.可得油污质点的速度表示为:V (x2,1)+,:(,1),质量守恒具有形式 2(x21)+(2,1)+p(V-HV)=0 为法向速度分量 初始物理构型 初始参数构型 曲面方程:E(x,D) 当前物理构型吃 曲面方程:E(x,4) 当前参数构型V 动刻画x=x(l 图8海面上油污扩散问题对应的构型构造:(左图)物理构型,(右图)参数构型 Fig 8 Sketch of the constructions of configurations with respect to oil diffusions on sea surfaces (left subplot) Physical configurations, (right subplot parametric configurations 动量守恒具有形式 切平面:P((,)+F+ ags(I C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
图７ 内边界上变形率张量特征问题分析：（左图）涡量与最大特征值之比沿内边界的分布．注：理论预测值为２或－２． （右图）最大特征值所对应的特征向量同内边界切向量的夹角沿内边界的分布．注：理论预测值为４５°或１３５°． Ｆｉｇ．７ Ｅｉｇｅｎ－ｐｒｏｂｌｅｍａｎａｌｙｓｉｓｏｆｔｈｅｓｔｒａｉｎｔｅｎｓｏｒｏｎｔｈｅｉｎｎｅｒｂｏｕｎｄａｒｙ：（ｌｅｆｔｓｕｂｐｌｏｔ）Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅｒａｔｉｏｏｆｔｈｅ ｖｏｒｔｉｃｉｔｙａｎｄｔｈｅｍａｘｉｍｕｍｅｉｇｅｎｖａｌｕｅｏｆｔｈｅｓｔｒａｉｎｔｅｎｓｏｒａｌｏｎｇｔｈｅｉｎｎｅｒｂｏｕｎｄａｒｙｏｆｔｈｅｕｎｄｕｌａｔｅｄｈｅｌｉｃｏｉｄａｌ ｓｕｒｆａｃｅ．Ｒｅｍａｒｋ：Ｔｈｅｔｈｅｏｒｅｔｉｃａｌｉｎｄｉｃａｔｅｄｖａｌｕｅｏｆｔｈｅｒａｔｉｏｉｓ２ｏｒ－２．（ｒｉｇｈｔｓｕｂｐｌｏｔ）Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅ ａｎｇｌｅｂｅｔｗｅｅｎｔｈｅｅｉｇｅｎｖｅｃｔｏｒｗｉｔｈｒｅｓｐｅｃｔｔｏｔｈｅｍａｘｉｍｕｍｅｉｇｅｎｖａｌｕｅｏｆｔｈｅｓｔｒａｉｎｔｅｎｓｏｒａｎｄｔｈｅｔａｎｇｅｎｔｖｅｃｔｏｒ ｏｆｔｈｅｂｏｕｎｄａｒｙａｌｏｎｇｔｈｅｉｎｎｅｒｂｏｕｎｄａｒｙｏｆｔｈｅｕｎｄｕｌａｔｅｄｈｅｌｉｃｏｉｄａｌｓｕｒｆａｃｅ．Ｒｅｍａｒｋ：Ｔｈｅｔｈｅｏｒｅｔｉｃａｌｉｎｄｉｃａｔｅｄ ｖａｌｕｅｏｆｔｈｅａｎｇｌｅｉｓ４５ｏｒ１３５ｄｅｇｒｅｅ． 固定曲面上的二维流动可以对应星体表面上的大气运动，或者地球表面上的洪水蔓延过洼地、山丘或 者盘山公路等，这些流动的法向尺度（流层厚度）远远小于流动的流向尺度，故可将流动视为二维．近期，作 者类比于一般涡量动力学［２４］，发展了固定曲面上二维流动的涡量动力学理论框架［２３］，希望能为研究此类 流动铺垫了一定的理论基础． ２．３．２ 海面上油污扩散 海面上油污扩散表现为几何形态为曲面的油污在海面上的运动过程，对应的物理构型以及参数构型 构造如图８所 示，现 运 动 刻 画 的 参 数 表 示 为 ｘΣ ＝ｘΣ（ξΣ，ｔ）∈ 瓗２．可得油污质点的速度表示为：Ｖ珝 ＝ Σ珝 ｔ （ｘΣ，ｔ）＋ｘｓ Σｇ珝ｓ，ｘｓ Σ∶＝ ｘｓ ｔ （ξΣ，ｔ）．质量守恒具有形式： ρ ｔ （ｘΣ，ｔ）＋ｘｓ Σ ρ ｘｓ Σ （ｘΣ，ｔ）＋ρ Δ ｓ Σ Ｖｓ （ ） －ＨＶ３ ＝０， Ｖ３ 为法向速度分量． 图８ 海面上油污扩散问题对应的构型构造：（左图）物理构型，（右图）参数构型 Ｆｉｇ．８ Ｓｋｅｔｃｈｏｆｔｈｅｃｏｎｓｔｒｕｃｔｉｏｎｓｏｆｃｏｎｆｉｇｕｒａｔｉｏｎｓｗｉｔｈｒｅｓｐｅｃｔｔｏｏｉｌｄｉｆｆｕｓｉｏｎｓｏｎｓｅａｓｕｒｆａｃｅｓ： （ｌｅｆｔｓｕｂｐｌｏｔ）Ｐｈｙｓｉｃａｌｃｏｎｆｉｇｕｒａｔｉｏｎｓ，（ｒｉｇｈｔｓｕｂｐｌｏｔ）ｐａｒａｍｅｔｒｉｃｃｏｎｆｉｇｕｒａｔｉｏｎｓ 动量守恒具有形式： 切平面：ρｘｌ Σ（ξΣ，ｔ）＋Γｌ ｐｑｘｐ Σｘｑ Σ ＋ ２ Σ珝 （ ） ｔ２ （ｘΣ，ｔ），ｇ珝ｌ 瓗３ ＋２ｘｓ Σ ｇ珝ｓ ｔ（ ） （ｘΣ，ｔ），ｇ珝ｌ （ ） 瓗３ ＝ 第４期 谢锡麟等：有限变形理论的若干进展及其在流体力学中的相关应用 ３５５
复旦学报(自然科学版) 第52卷 -2(,0)+m1+(5V++V-32-2()v-bv 法方向 berapit +2 Crv. t, n (y-p)H++(v+3+(:)V-2bvy 上述连续性方程及动量守恒可构成 Lagrange形式的控制方程组,可联立求解变形的参数刻画x= x(,t)∈R2,面密度分布以及内压力分布.实际进行数值研究时,可独立考虑各种作用的行为,包括油污 自身间的内压力以及内摩擦作用、油污同海面之间的摩擦作用以及油污重力等. 图9与图10显示了海面做行波运动时油污的扩散特征,设定海面运动形式为:z(x,y,t 0.1sin(-x-at),其中A=1.0,a=40 433 2.001.00 0.40 0.20 0.20 -0.500.00 0.00 0.20 0y0R0999o0r0+ -1.00-0.20 -0.80 3.50 O606002020w0660100 -100-1.00 图9海面做行波运动时的几何量空间分布(t=0.30):(左图)平均曲率H,(右图)HV Fig. 9 Spatial distributions of geometrical quantities in the case that the sea surface does traveling wave motion (r=0. 30):(left subplot )mean curvature H, (right subplot ) HV 0.20 0.20 00 010660a 1.00-100 图10海面做行波运动时的油污物理量的空间分布(t=0.40):(左图)密度,(右图)速率 Fig. 10 Spatial distributions of physical quantities of the oil in the case that the sea surface does traveling wave motion(t=0. 40):(left subplot )density,(right subplot)velocity magnitude 图11与图12显示了海面做驻波运动时油污的扩散特征,设定海面运动形式为:运动曲面方程 z(r, y,t)=Q Isin/ 2, )…sn(xy)sna,其中=.0,m=.0 上述数值研究中考虑了油污的重力作用,但尚未考虑油污自身间的内压力以及内摩擦作用,未考虑油 污同海面之间的摩擦作用.现有研究均表明初始厚度均匀分布的油污,随着时间演化会发生集聚现象 C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
－ ｐ ｘｌ Σ （ｘΣ，ｔ）＋ρｆｌ Σ ＋μ Σｌ ΔΔ ｓ Σ Ｖｓ ＋ Σｓ ΔΣ Δ ｓＶｌ ＋ＫＧＶｌ －３ｂｌｓ Ｖ３ ｘｓ Σ －２（Σｓ Δ ｂ ）ｌ ｓ Ｖ３ －ｂｌ ｓｂｓ （ ） ｔＶｔ ， 法方向： ρｂｐｑｘｐ Σｘｑ Σ ＋ ２ Σ珝 （ ） ｔ２ （ｘΣ，ｔ），ｎ 瓗３ ＋２ｘｓ Σ ｇ珝ｓ ｔ（ ） （ｘΣ，ｔ），ｎ （ ） 瓗３ ＝ （γ－ｐ）Ｈ ＋ρｆ３ Σ ＋μ（Σｓ ΔΔ ｓ Σ Ｖ３ ＋３ｂｓ ｔ Σ Δ ｓＶｔ ＋ （Σ Δ ｑｂｑ ｓ ）Ｖｓ －２ｂｓｔｂｓｔＶ ）３ ． 上述连续性方程及动量守恒可 构 成 Ｌａｇｒａｎｇｅ形 式 的 控 制 方 程 组，可 联 立 求 解 变 形 的 参 数 刻 画 ｘ＝ ｘ（ξ，ｔ）∈瓗２，面密度分布以及内压力分布．实际进行数值研究时，可独立考虑各种作用的行为，包括油污 自身间的内压力以及内摩擦作用、油污同海面之间的摩擦作用以及油污重力等． 图９与 图 １０ 显示了海面做行波 运动时油污的扩散特征，设定海面运动 形式为：ｚ（ｘ，ｙ，ｔ）＝ ０．１ｓｉｎ（ ２π λｘ－ωｔ），其中λ＝１．０，ω＝４．０． 图１１与图１２显示 了 海 面 做 驻 波 运 动 时 油 污 的 扩 散 特 征，设 定 海 面 运 动 形 式 为：运 动 曲 面 方 程： ｚ（ｘ，ｙ，ｔ）＝０．１ｓｉｎ（ ２π λ１ ｘ）·ｓｉｎ（ ２π λ２ ｙ）·ｓｉｎωｔ，其中λ＝１．０，ω＝４．０． 上述数值研究中考虑了油污的重力作用，但尚未考虑油污自身间的内压力以及内摩擦作用，未考虑油 污同海面之间的摩擦作用．现有研究均表明初始厚度均匀分布的油污，随着时间演化会发生集聚现象． ４５５ 复 旦 学 报（自然科学版） 第５２卷
第4期 谢锡麟等:有限变形理论的若干进展及其在流体力学中的相关应用 555 曰250 0.40 -0.50 0.20 0 -6.50 7.50 图11海面做驻波运动时的几何量空间分布(t=0.30):(左图)平均曲率,(右图)HV g. 11 Spatial distributions of geometrical quantities in the case that the sea surface does stationary wave motion (t=0. 30):(left subplot )mean curvature, (right subplot) HV 0099040 40 (m06 0.60 图12海面做驻波运动时油污物理量的空间分布(t=0.20):(左图)密度,(右图)速率 Fig. 12 Spatial distributions of physical quantities of the oil in the case that the sea surface does stationary wave motion(t=0. 20): (left subplot )density, (right subplot )velocity magnitude 3总结及讨论 本文阐述了作者新近发展的两套有限变形理论:当前物理构型对应之曲线坐标系显含时间的有限变 形理论,几何形态为曲面的连续介质的有限变形理论,分别对应连续介质几何形态为 Euclid流形和 Riemann流形,且给出了相关应用事例 当前物理构型对应之曲线坐标系显含时间的有限变形理论在一般连续介质力学的理论体系中,引 入初始物理构型以及当前物理构型,对二者可再分别引入初始参数构型以及当前参数构型,物理构型与参 数构型之间的关系即为一般曲线坐标系.为研究边界的有限变形运动对介质运动的影响,我们对当前物理 构型引入显含时间的曲线坐标系,表现为时空空间中的微分同胚.通过构造适当的曲线坐标系可将物理空 中几何形态不规则且随时间变化的运动区域微分同胚至参数空间中的几何形态规则且不随时间变化的 参数区域.进一步将连续介质的控制方程按曲线坐标系的局部基展开就可获得定义于参数区域上的控制 方程.将控制方程按局部基展开有益于研究边界或流场几何特性与流动行为之间的关系.就此理论,我们 发展了对应曲线坐标系显含时间的流函数-涡量解法,并应用于研究边界可作有限变形运动的钝体绕流 几何形态为曲面的连续介质的有限变形理论诸如细胞膜变形运动、平面或曲面皂膜上的流动、海面 上的油污扩散以及星体表面的大范围大气运动等,均表现为介质的法向尺度(厚度)远远小于其展向尺度, C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
３ 总结及讨论 本文阐述了作者新近发展的两套有限变形理论：当前物理构型对应之曲线坐标系显含时间的有限变 形理论，几何形态为曲面的连 续介质的有限变形理论，分别对应连续介质几何形态为 Ｅｕｃｌｉｄ 流 形 和 Ｒｉｅｍａｎｎ流形，且给出了相关应用事例． 当前物理构型对应之曲线坐标系显含时间的有限变形理论 在一般连续介质力学的理论体系中，引 入初始物理构型以及当前物理构型，对二者可再分别引入初始参数构型以及当前参数构型，物理构型与参 数构型之间的关系即为一般曲线坐标系．为研究边界的有限变形运动对介质运动的影响，我们对当前物理 构型引入显含时间的曲线坐标系，表现为时空空间中的微分同胚．通过构造适当的曲线坐标系可将物理空 间中几何形态不规则且随时间变化的运动区域微分同胚至参数空间中的几何形态规则且不随时间变化的 参数区域．进一步将连续介质的控制方程按曲线坐标系的局部基展开就可获得定义于参数区域上的控制 方程．将控制方程按局部基展开有益于研究边界或流场几何特性与流动行为之间的关系．就此理论，我们 发展了对应曲线坐标系显含时间的流函数 涡量解法，并应用于研究边界可作有限变形运动的钝体绕流． 几何形态为曲面的连续介质的有限变形理论 诸如细胞膜变形运动、平面或曲面皂膜上的流动、海面 上的油污扩散以及星体表面的大范围大气运动等，均表现为介质的法向尺度（厚度）远远小于其展向尺度， 第４期 谢锡麟等：有限变形理论的若干进展及其在流体力学中的相关应用 ５５５
556 复旦学报(自然科学版) 第52卷 由此可将此类介质的几何形态模型化为三维 Euclid空间中的光滑曲面,亦即二维 Riemann光滑流形.按 微分几何的观点, Euclid流形(平坦空间)同 Riemann流形(弯曲空间)在几何性质上有本质区别,对此在 场论上亦有显著差异.按现有理论,几何形态为曲面的连续介质其守恒律控制方程就直接体现为曲面几何 量同物理量之间的耦合作用.就此理论,我们发展了固定曲面上二维不可压缩流动的流函数-涡量解法,并 进行了固定曲面上圆柱绕流以及内流的数值研究;给出了海面上油污扩散的 Lagrange型控制方程及其数 值研究 本文所述有限变形理论的相关思想及方法的适用性及有效性需待更为广泛的数值实验及真实实验的 检验 参考文献: [1] Triantafyllou M S, Triantafyllou G S, Yue D K P. Hydrodynamics of fishlike swimming[J]. Annu Rev [2] Fish F E, Lauder G V. Passive and active flow control by swimming fishes and mammals[J]. Annu Rev Fluid mech,2006,38:193-224. [3] WuC. Wang L. Adaptive optimal control of the flapping rule of a fixed flapping plate[J]. Advances in Applied Mathematics and Mechanics, 2009, 1(3): 402-414 [4] LuX Y, Yin X Z, Yang J M, et al. Studies of hydrodynamics in fishlike swimming propulsion[J] [5] WuCJ, Wang L. Numerical simulations of self-propelled swimming of 3D bionic fish school[J]. Science [6] WuC J. Wang L. Where is the rudder of a fish?-The mechanism of swimming and control of self- propelled fish school[J]. Acta Mechanica Sinica, 2010, 26(1):45-65. [7] Dong G J, Lu X Y. Characteristics of flow over traveling wavy foils in a side-by-side arrangement[J] Physics of Fluids, 2007, 19: 057107 [8] Williamson C H K. Evolution of a single wake behind a pair of bluff bodies[J]. J Fluid Mech, 1985, 159: 1-18 [9] Govardhan R, Williamson C K. Modes of vortex formation and frequency response of a freely vibrating ylinder[J]. J Fluid Mech, 2000,420(1):85-130 [10] Du G, Sum M. Effects of unsteady deformation of flapping wing on its aerodynamics forces[J]. Applied Mathematics and Mechanics, 2008, 29(6):731-743. [11] Lu X Y, Yin XZ. Propulsive performance of a fishlike traveling wavy wall[J]. Acta Mechanica, 2005 175(1-4):197-215 [12] WuCJ, Xie Y Q, Wu J Z "Fluid roller bearing"effect and flow control[J]. Acta Mechanica Sinca 03,19(5):476-484. [13] Wu CJ, Wang L, Wu J Z. Suppression of the von Karman vortex street behind a circular cylinder by a travelling wave generated by a flexible surface[J]. J Fluid Mech,2007,574:365-391 [14]力学情报.动力叶轮机械气动力学专集[M].北京:中国科学院北京力学研究所,197 [15]李开泰,黄艾香.张量分析及其应用[M].北京:科学出版社,2004 [16]谢锡麟.现代张量分析在连续介质力学中的若干应用[C]∥吴有生,周如苹,颜开,等,第十一届全国水 动力学学术会议暨第二十四届全国水动力学研讨会并周培源教授诞辰110周年纪念大会文集,北京:海 洋出版社,2012:224-236 [17]郭仲衡,非线性弹性理论[M].北京:科学出版社,1980 [18] Xie X L, Chen Y, Shi Q. Some studies on mechanics of continuous mediums viewed as differential manifolds[J]. Sci China-Phys Mech Astron, 2013, 56(2):432-456 [19]陈瑜,谢锡麟,麻伟巍,二维类圆柱边界的有限变形运动对其尾迹空间动力学行为的影响[C]∥吴有 生,周如苹,颜开,等.第十一届全国水动力学学术会议暨第二十四届全国水动力学研讨会并周培源教 授诞辰110周年纪念大会文集.北京:海洋出版社,2012:212-223 [20] Boffetta G, Ecke R E. Two-dimensional turbulence[J]. Ammu Rev Fluid Mech,2012,44:427-451 C1994-2013ChinaAcademicJOurnalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
由此可将此类介质的几何形态模型化为三维 Ｅｕｃｌｉｄ空间中的光滑曲面，亦即二维 Ｒｉｅｍａｎｎ光滑流形．按 微分几何的观点，Ｅｕｃｌｉｄ流形（平坦空间）同 Ｒｉｅｍａｎｎ流形（弯曲空间）在几何性质上有本质区别，对此在 场论上亦有显著差异．按现有理论，几何形态为曲面的连续介质其守恒律控制方程就直接体现为曲面几何 量同物理量之间的耦合作用．就此理论，我们发展了固定曲面上二维不可压缩流动的流函数 涡量解法，并 进行了固定曲面上圆柱绕流以及内流的数值研究；给出了海面上油污扩散的 Ｌａｇｒａｎｇｅ型控制方程及其数 值研究． 本文所述有限变形理论的相关思想及方法的适用性及有效性需待更为广泛的数值实验及真实实验的 检验． 参考文献： ［１］ ＴｒｉａｎｔａｆｙｌｌｏｕＭＳ，ＴｒｉａｎｔａｆｙｌｌｏｕＧＳ，ＹｕｅＤ ＫＰ．Ｈｙｄｒｏｄｙｎａｍｉｃｓｏｆｆｉｓｈｌｉｋｅｓｗｉｍｍｉｎｇ［Ｊ］．ＡｎｎｕＲｅｖ ＦｌｕｉｄＭｅｃｈ，２０００，３２：３３－５３． ［２］ ＦｉｓｈＦＥ，ＬａｕｄｅｒＧＶ．Ｐａｓｓｉｖｅａｎｄａｃｔｉｖｅｆｌｏｗｃｏｎｔｒｏｌｂｙｓｗｉｍｍｉｎｇｆｉｓｈｅｓａｎｄｍａｍｍａｌｓ［Ｊ］．ＡｎｎｕＲｅｖ ＦｌｕｉｄＭｅｃｈ，２００６，３８：１９３－２２４． ［３］ ＷｕＣＪ，ＷａｎｇＬ．Ａｄａｐｔｉｖｅｏｐｔｉｍａｌｃｏｎｔｒｏｌｏｆｔｈｅｆｌａｐｐｉｎｇｒｕｌｅｏｆａｆｉｘｅｄｆｌａｐｐｉｎｇｐｌａｔｅ［Ｊ］．Ａｄｖａｎｃｅｓｉｎ ＡｐｐｌｉｅｄＭａｔｈｅｍａｔｉｃｓａｎｄＭｅｃｈａｎｉｃｓ，２００９，１（３）：４０２－４１４． ［４］ ＬｕＸ Ｙ，ＹｉｎＸ Ｚ，ＹａｎｇＪ Ｍ，ｅｔａｌ．Ｓｔｕｄｉｅｓｏｆｈｙｄｒｏｄｙｎａｍｉｃｓｉｎｆｉｓｈｌｉｋｅｓｗｉｍｍｉｎｇｐｒｏｐｕｌｓｉｏｎ［Ｊ］． Ｊｏｕｒｎａｌｏｆ Ｈｙｄｒｏｄｙｎａｍｉｃｓ，２０１０，２２（５）：１２－２２． ［５］ ＷｕＣＪ，ＷａｎｇＬ．Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｓｏｆｓｅｌｆ－ｐｒｏｐｅｌｌｅｄｓｗｉｍｍｉｎｇｏｆ３Ｄｂｉｏｎｉｃｆｉｓｈｓｃｈｏｏｌ［Ｊ］．Ｓｃｉｅｎｃｅ ｉｎＣｈｉｎａ（Ｅ），２００９，５２（３）：６５８－６６９． ［６］ ＷｕＣＪ，ＷａｎｇＬ．Ｗｈｅｒｅｉｓｔｈｅｒｕｄｄｅｒｏｆａｆｉｓｈ？－Ｔｈｅ ｍｅｃｈａｎｉｓｍ ｏｆｓｗｉｍｍｉｎｇａｎｄｃｏｎｔｒｏｌｏｆｓｅｌｆ－ ｐｒｏｐｅｌｌｅｄｆｉｓｈｓｃｈｏｏｌ［Ｊ］．ＡｃｔａＭｅｃｈａｎｉｃａＳｉｎｉｃａ，２０１０，２６（１）：４５－６５． ［７］ ＤｏｎｇＧＪ，ＬｕＸ Ｙ．Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｆｌｏｗｏｖｅｒｔｒａｖｅｌｉｎｇｗａｖｙｆｏｉｌｓｉｎａｓｉｄｅ－ｂｙ－ｓｉｄｅａｒｒａｎｇｅｍｅｎｔ［Ｊ］． ＰｈｙｓｉｃｓｏｆＦｌｕｉｄｓ，２００７，１９：０５７１０７． ［８］ ＷｉｌｌｉａｍｓｏｎＣＨＫ．Ｅｖｏｌｕｔｉｏｎｏｆａｓｉｎｇｌｅｗａｋｅｂｅｈｉｎｄａｐａｉｒｏｆｂｌｕｆｆｂｏｄｉｅｓ［Ｊ］．ＪＦｌｕｉｄＭｅｃｈ，１９８５，１５９：１－１８． ［９］ ＧｏｖａｒｄｈａｎＲ，ＷｉｌｌｉａｍｓｏｎＣ Ｈ Ｋ．Ｍｏｄｅｓｏｆｖｏｒｔｅｘｆｏｒｍａｔｉｏｎａｎｄｆｒｅｑｕｅｎｃｙｒｅｓｐｏｎｓｅｏｆａｆｒｅｅｌｙｖｉｂｒａｔｉｎｇ ｃｙｌｉｎｄｅｒ［Ｊ］．ＪＦｌｕｉｄＭｅｃｈ，２０００，４２０（１）：８５－１３０． ［１０］ ＤｕＧ，Ｓｕｍ Ｍ．Ｅｆｆｅｃｔｓｏｆｕｎｓｔｅａｄｙｄｅｆｏｒｍａｔｉｏｎｏｆｆｌａｐｐｉｎｇｗｉｎｇｏｎｉｔｓａｅｒｏｄｙｎａｍｉｃｓｆｏｒｃｅｓ［Ｊ］．Ａｐｐｌｉｅｄ ＭａｔｈｅｍａｔｉｃｓａｎｄＭｅｃｈａｎｉｃｓ，２００８，２９（６）：７３１－７４３． ［１１］ ＬｕＸＹ，ＹｉｎＸＺ．Ｐｒｏｐｕｌｓｉｖｅｐｅｒｆｏｒｍａｎｃｅｏｆａｆｉｓｈ－ｌｉｋｅｔｒａｖｅｌｉｎｇｗａｖｙｗａｌｌ［Ｊ］．ＡｃｔａＭｅｃｈａｎｉｃａ，２００５， １７５（１－４）：１９７－２１５． ［１２］ ＷｕＣＪ，ＸｉｅＹ Ｑ，ＷｕＪＺ．“Ｆｌｕｉｄｒｏｌｌｅｒｂｅａｒｉｎｇ”ｅｆｆｅｃｔａｎｄｆｌｏｗｃｏｎｔｒｏｌ［Ｊ］．ＡｃｔａＭｅｃｈａｎｉｃａＳｉｎｃａ， ２００３，１９（５）：４７６－４８４． ［１３］ ＷｕＣＪ，ＷａｎｇＬ，ＷｕＪＺ．ＳｕｐｐｒｅｓｓｉｏｎｏｆｔｈｅｖｏｎＫａｒｍａｎｖｏｒｔｅｘｓｔｒｅｅｔｂｅｈｉｎｄａｃｉｒｃｕｌａｒｃｙｌｉｎｄｅｒｂｙａ ｔｒａｖｅｌｌｉｎｇｗａｖｅｇｅｎｅｒａｔｅｄｂｙａｆｌｅｘｉｂｌｅｓｕｒｆａｃｅ［Ｊ］．ＪＦｌｕｉｄＭｅｃｈ，２００７，５７４：３６５－３９１． ［１４］ 力学情报．动力叶轮机械气动力学专集［Ｍ］．北京：中国科学院北京力学研究所，１９７６． ［１５］ 李开泰，黄艾香．张量分析及其应用［Ｍ］．北京：科学出版社，２００４． ［１６］ 谢锡麟．现代张量分析在连续介质力学中的若干应用［Ｃ］∥吴有生，周如苹，颜 开，等．第十一届全国水 动力学学术会议暨第二十四届全国水动力学研讨会并周培源教授诞辰１１０周年纪念大会文集．北京：海 洋出版社，２０１２：２２４－２３６． ［１７］ 郭仲衡．非线性弹性理论［Ｍ］．北京：科学出版社，１９８０． ［１８］ ＸｉｅＸ Ｌ，Ｃｈｅｎ Ｙ，ＳｈｉＱ．Ｓｏｍｅｓｔｕｄｉｅｓｏｎ ｍｅｃｈａｎｉｃｓｏｆｃｏｎｔｉｎｕｏｕｓ ｍｅｄｉｕｍｓｖｉｅｗｅｄａｓｄｉｆｆｅｒｅｎｔｉａｌ ｍａｎｉｆｏｌｄｓ［Ｊ］．ＳｃｉＣｈｉｎａ－ＰｈｙｓＭｅｃｈＡｓｔｒｏｎ，２０１３，５６（２）：４３２－４５６． ［１９］ 陈 瑜，谢锡麟，麻伟巍．二维类圆柱边界的有限变形运动对其尾迹空间动力学行为的影 响［Ｃ］∥吴 有 生，周如苹，颜 开，等．第十一届全国水动力学学术会议暨第二十四届全国水动力学研讨会并周培源教 授诞辰１１０周年纪念大会文集．北京：海洋出版社，２０１２：２１２－２２３． ［２０］ ＢｏｆｆｅｔｔａＧ，ＥｃｋｅＲＥ．Ｔｗｏ－ｄｉｍｅｎｓｉｏｎａｌｔｕｒｂｕｌｅｎｃｅ［Ｊ］．ＡｎｎｕＲｅｖＦｌｕｉｄＭｅｃｈ，２０１２，４４：４２７－４５１． ６５５ 复 旦 学 报（自然科学版） 第５２卷