第十五届全国分离流、旋涡和流动控制会议 2014年10月24-26日北京 可变形壁面上涡量动力学相关理论结果及应用 —变形率、切应力、涡量法向梯度 复旦大学力学与工程科学系 谢锡麟陈瑜
第十五届全国分离流、旋涡和流动控制会议 2014年10月24-26日 北京 可变形壁面上涡量动力学相关理论结果及应用 —— 变形率、切应力、涡量法向梯度 复旦大学 力学与工程科学系 谢锡麟 陈瑜
几身时化映 何诱间的照 X(5,n,5,)= t +h 5 同导变物观 力的化理点 学局的区 之部参域 间基数微 7 的展区分将 关开域同几 叶片间流动区域 参数区域 (边界始终固定) 系张‘胚何 旋转轴(边界可作有限变形运动) 量⌒至形 控几态 =0(51)g8g8g(x,1) 制基何不 参数区域 方于形规 x-曲线、8() =(x,1)g8g8g() 程曲态则 x"k曲线 g2(x2,1) 线规且 曲线 曲线 曲线 8.(=,n) 以坐则随 g,(xa x”-曲线 x2-曲线 期标且时物理区域 x∠曲线 建系不间 x2-曲线 立自随变 曲线坐标X=X(x,)
映 照 观 点 ( 1 ) 将 几 何 形 态 不 规 则 且 随 时 间 变 化 的 物 理 区 域 微 分 同 胚 至 几 何 形 态 规 则 且 不 随 时 间 变 化 的 参 数 区 域 ; ( 2 ) 基 于 曲 线 坐 标 系 自 身 诱 导 的 局 部 基 展 开 张 量 控 制 方 程 , 以 期 建 立 几 何 同 力 学 之 间 的 关 系
X(x)DxR3(=→()=x1()+=(:8小+n(起小少 于 号m人 curve line 曲 面 X(rr,5 g2 x x-line 的 K-s 1 curve 维曲面梯度算子 正 VoΦg )op=g0(g,8g)( 交 系 Φ,g⊙g⑧g+Φ,bg⊙n⑧g+Φ,b!g^⊙g,⑧n 全空间梯度算子vs; V+n aX a(x25=0) Lei算]va(g,8g) g1∞g
基 于 曲 面 的 三 维 半 正 交 系 1 1 2 2 1 2 1 2 3 X , : , , X , , , , , , , , , x x X x t D x t x t x t t X x t t x x t n x x t t X i , 0, X n x t , , k i j k j i i k j i k k j i j k k j i j ki j k i k g x t g g g x x t x g g g b g n g b g g n k i j i j j i k j C x i g g g g 曲面梯度算子 全空间梯度算子 Levi-Civita算子
三维不可压缩流动解决策略—边界条件处理 (基于 E Weinan方法以及本组关于曲面上二类梯度算子的相关理论) 基于曲面的半正交基{g1,82,m Helmholtz- Stokes分解: 定义:=vg1+yg2∈ VV=0 VV=O V=Vxv,满足 V×=Vn,BC|V×n=Vnar+H=0 按曲面切平面中的 helmholtz- Stokes分解ie2013) △s=V.b b=V5+V×(mn),此处 Vb∈T∑ △m=(V×b)·m y=V×(mn)其中△n=Vn (r×n)Vmd=/v.(vndo V. ndo
1 2 1 2 g g T 基于曲面的半正交基g n 1 2 , g , · 0 · B.C. Helmholtz Stok V V n es 分解: ,满足 , · 0;· · C C n V n 3 3 3 H 0 x 定义: 按曲面切平面中的Helmholtz-Stokes 分解(Xie 2013) V n· C n 其中 三维不可压缩流动解决策略 ——边界条件 处理 (基于 E.Weinan 方法 以及 本组关于曲面上二类梯度算子的相关理论)
固定壁面情形 可变形壁面上 变形率张量 「Casw公式适用 (吴介之理论 切平面 涡线 D (V③+VV) O×n 2 D=0-vpn8n+3(O×n)n+n8(O×n) 「Casw公式 ∑ +(F×n)8n+n8(Wxn)-n1nV,此处W=1V87,n,n 纯壁面变形影响 注:1.相关理论研究为复杂流场的数值研究提供了“自我检查”的理论依据。 2.从方法论而言,在揭示物理意义方面, Cartesian典则基是最简单的基,但 并非一定就是最合适的基
可变形壁面上 变形率张量 (吴介之理论) Caswell 公式 纯壁面变形影响 注:1. 相关理论研究为复杂流场的数值研究提供了“自我检查”的理论依据。 2. 从方法论而言,在揭示物理意义方面,Cartesian典则基是最简单的基,但 并非一定就是最合适的基。 1 : 2 D V V Caswell 公式适用