2013年中国力学大会：基于几何形态为曲面的连续介质的有限变形理论研究一般固定曲面上的流动

2013年8月中国力学大会-2013 基于几何形态为曲面的连续介质的有限变形理论 研究一般固定曲面上的流动 复旦大学力学与工程科学系 谢锡麟陈瑜史倩 研究背景、研究思路 理论研究:一般固定曲面上的流动(连续介质作为二维 Riemann微分流形) ·理论研究:固定曲面上二维流动的涡量动力学理论框架 典型运动事例:凹凸螺旋面上的流动,固定曲面上圆柱绕流

2013年8月 中国力学大会-2013 基于几何形态为曲面的连续介质的有限变形理论 研究一般固定曲面上的流动 复旦大学 力学与工程科学系 谢锡麟 陈瑜 史倩 • 研究背景、研究思路 • 理论研究：一般固定曲面上的流动（连续介质作为二维Riemann微分流形） • 理论研究：固定曲面上二维流动的涡量动力学理论框架 • 典型运动事例：凹凸螺旋面上的流动，固定曲面上圆柱绕流

星体上流体运动 固 定曲面上二维流动

星体上流体运动 —— 固定曲面上二维流动

曲面上两类梯度算子 曲面梯 度算子 g1∞g 8'o-0og6g)(x)=g g1∞g 8, 0g+,(Ti8+b, n)og+ 8, 8(-rik8+bi'n v, a,(go-g0g'a b( o-n)og+a bi(go-g )o 曲面梯度算子定义基于一般张量赋范空间上的微分学 02Φ 02Φ 可以应用于在曲面上有定义的所有张量场 axa Levi-Civitavo→Φ=go-Va(Φg,g)(x) 梯度算子 d g +g,8g11s19P-④,g,③Ig=V(8°-8)8g Levi- Civita联络算子仅对切平面上的分量有效 VV2@-VV中=Rm①,+R0①, where rs=bbn分:Vb=Vb

曲面上两类梯度算子 Levi-Civita 梯度算子                 . . . . . . ( ) ( ) i j i j i j j i j i j i i l j i k j i j k j j i j l i l l i k li j i lk l i l j i j l j li j l i j l j i l l l l g g g g x g g x g g g g b n g g g b n x b g g x x g g g g n g g b g n                                                              曲面 梯 度算子     . . . . ( ) l i j j i i l j i k j i j k j l i j li k j i lk i l j l i l j x g g x g g g g g g g x g g g g                                   曲面梯度算子定义基于一般张量赋范空间上的微分学 可以应用于在曲面上有定义的所有张量场 Levi-Civita 联络算子仅对切平面上的分量有效     2 2 q p p q x x x x x x          , : - = i i i t s i s s s q p j p q j tqp j j qp t tqp q tp tq p q ps p qs R R where R b b b b b b                      ；

曲面梯度算子的应用之 内蕴形式广义sces公式」 z×n)。-①d=「Vo-①+Hno-①da t xn ∑ Φo-(x×n)dl=|o-V+HΦo-ndo 质量守恒 Xie et al. Science China g 2013 at (x,t)da+p:n·(V)dl at (a, t)do+/V(pv)do 动量守恒 apv) at (a, t)do+y, n(pv)Vdl=/ pa do=Ften+Fpre+Fvis+Fsur 表面张力Fm:=74r×m=7/Hno 内压力 Fpre:=-r×(pn)l= p- pHn do 内摩擦Fm:=4p(rxn)·(v+四av)m V·(vV+Vv)+Ht 0

曲面梯度算子的应用之一 —— 内蕴形式广义Stokes公式 表面张力 内压力 动量守恒 质量守恒 内摩擦 Xie et al. Science China G,2013

Mass conservation (x,t)+V·(p)、 at a(x,、、8(x,t)+pV=p+p0=0,:=Vv Novel NsEs Curvatures roles m(+7+7)+ka可+fm pan=Hr-p)+u[26v: i]+four (ay, t)+VV,Vi 0 (a, t)+u[9ViV, Vi+VI(vsVs)+Kgvi+fsi biivivj= H(y-p)+u[,+four Gauss曲率伴随着速度在切平面的分量直接参与至切平面上的动量平衡 而平均曲率伴随着表面张力系数、内压力参与至法向动量平衡

Novel NSEs Curvatures Roles Mass Conservation Gauss 曲率伴随着速度在切平面的分量直接参与至切平面上的动量平衡； 而平均曲率伴随着表面张力系数、内压力参与至法向动量平衡