文章快速检索  
  高级检索
跨介质运动物体的附加质量
张晓强1,2, 冯金富1, 吝科1, 胡俊华1    
1. 空军工程大学 航空航天工程学院, 西安 710038;
2. 中国人民解放军95889部队, 酒泉 735018
摘要:跨介质运动物体具备能够适应多种流体介质环境的突出优势,但由于受自由面及流体影响,力学环境复杂,其运动特性难以预测。附加质量法是工程上解决物体与流体耦合问题的常用方法,但现有的附加质量结论无法直接推广到带自由面问题的跨介质过程。针对跨介质过程的时变附加质量问题,提出了一种新的求解方法,研究了不同浸深时某柱体外形跨介质过程中出水阶段和入水阶段的附加质量变化规律。结论表明:出水阶段和入水阶段的附加质量随浸深的变化趋势一致,曲线相似;相同浸深状态时出水阶段和入水阶段的法向附加质量相近,线性特征明显;相同浸深状态时出水阶段和入水阶段的轴向附加质量、附加转动惯量和附加静距差距明显,2个过程不可逆。本文所提出的方法可为各类跨介质结构体设计提供技术支持,也可推广至其他带自由面问题的附加质量求解中。
关键词: 跨介质     附加质量     自由面流动     水动力     数值计算    
Added mass of trans-media moving object
ZHANG Xiaoqiang1,2, FENG Jinfu1 , LIN Ke1, HU Junhua1     
1. Aeronautics and Astronautics Engineering College, Air Force Engineering University, Xi'an 710038, China;
2. PLA 95889, Jiuquan 735018, China
Abstract: Trans-media moving object was equipped with outstanding advantages in multiple fluid media environment. However, its moving characteristics were hard to predict because of the free surface problem and complex stress. Added mass method was common to solve fluid-structure coupling problems in engineering, but the existing conclusions could not be extended to trans-media process with free surface. A new method was proposed aiming to solve the time-varying added mass problems in trans-media process. Changing laws of added mass of some cylinder were studied, which lied in different depths in both separation step and entry step of the trans-media process. It is shown that the two steps have same changing trend and similar depth-varying curve with the depth varying; in the same depth, normal added mass of separation step is near to that of the entry step and both depth-varying curves have linear characteristics; as for axial added mass, added inertia moment, and added static moment, the results of two steps in the same depth are different obviously, so the two steps are irreversible. The proposed method is helpful for design of trans-media structure, and can be extended to other situations of added mass computing with free surface.
Key words: trans-media     added mass     free surface flow     hydrodynamic force     numerical calculation    

跨介质是近年来提出的一种新概念,用于航行器、飞行器和潜射类兵器,其主要特点在于可以跨越界面从而实现在不同介质环境中的有效运动,能够隐蔽高效突防、执行多重任务[1]。此概念可追溯到20世纪中期苏联、美国先后都曾开展过的“飞行潜艇”项目,据报道,2008年,美国国防高级研究计划局(DARPA)重启此项计划,这一新式军事装置将能像潜艇一样在水中行进,又能像飞机一样在空中飞行,前景十分广阔,但技术上面临着诸如外形、结构、动力和材料等诸多难题,可检索到的文献不多,主要以探索性研究、弹道仿真为主[2, 3, 4]。出水阶段和入水阶段是物体跨介质运动的关键阶段,运动过程中,物体受自由面流动、介质变化及耦合作用等因素影响[5, 6],受力情况复杂,附加质量法虽然是工程上解决此类问题的成熟办法,在水下机器人控制、结构体参数辨识等方面应用较多[7, 8],但跨介质过程中附加质量的时变特性极大地增加了研究的难度,成为其模型构建和运动特性分析的关键制约因素。传统的理论推导法只适用于无穷域简单外形物体[9]的附加质量求解,水池试验方法周期长、花费高,而数值方法无比尺效应、流场信息丰富,国内外对于自由面流动的数值模拟已经做了大量的研究工作[10, 11, 12],这也为时变附加质量的求解奠定了基础。目前对于运动物体的附加质量研究主要针对水下运动体和船舶类物体[13, 14, 15, 16, 17],而对于跨介质运动物体附加质量的研究比较少,主要包括强迫振荡法及其数值化方法,不足之处在于振荡所引入的外界能量会导致自由液面以及流场的变化,与实际运动环境不符,增大了计算误差。本文针对运动物体跨越自由界面过程中的时变附加质量求解问题,采用约束数值试验手段,提出一种快速有效的计算策略,并以球头柱体为对象研究了其跨介质过程中的附加质量变化规律,可以为新概念飞行器/武器的总体设计及优化提供技术支撑。

1 附加质量计算策略

以物体水空跨越为例,建立如图 1所示坐标系。

图 1 倾斜出水过程的坐标系 Fig. 1 Coordinate system of aslant separation

参考文献[7],给出物体倾斜水空跨越动力学方程:

式中:ρw为水的密度;m、v、S、LJz分别为物体的质量、速度、特征截面积、特征长度和转动惯量;vx、vyvz分别为物体速度在3个坐标轴方向上的分量;Cx、Cy、mz、ωzθ分别为物体x轴方向的阻力系数、y轴方向的阻力系数、绕z轴的转动力矩系数、绕z轴的转动角速度和运动倾角;xByg分别为出水过程物体浮心、质心的坐标;T、GB分别为出水过程中物体所受的推力、重力和浮力;λ11λ22λ26λ66为出水过程的附加质量,即轴向附加质量、法向附加质量、附加静矩和附加转动惯量。

vy=0,ωz=0,物体沿x轴方向运动,方程组可简化为

vx=0,ωz=0,物体沿y轴方向运动,方程组可简化为

vx=0,vy=0,物体绕z轴转动,方程组可简化为

式(4)~式(6)右端后4项代表物体所受的水作用力/力矩,改变加速度与角加速度的方向,单独写出,如式(7)~式(12)所示:

式中:Fx+、Fx-、Fy+Fy-分别为物体沿x轴正向、x轴负向、y轴正向和y轴负向运动时所受的水作用力;My+、My-、Mω+Mω-分别为物体沿y轴正向运动、y轴负向运动、绕z轴顺时针转动和绕z轴逆时针转动时所受的水作用力矩。

于是可得

由流体力学的知识易知,若能获取物体表面的流体信息,那么通过对物体表面的压强进行面积分即可得出Fx+、Fx-、Fy+Fy-,通过对物体表面的压强同时进行面积分和距离积分即可得出My+、My-、Mω+Mω-,附加质量计算公式虽然由出水动力学方程引出,但从过程来看,只要能确定不同运动参数下的流场环境,物体是出水状态还是入水状态均不影响结论的成立,故式(13)~式(16)同样适合于入水过程。

综上可得出附加质量的计算策略:建立运动物体的流场数值模型,设定2组速度相同、加速度相反的运动规律(v,v′,ω,ω′)与(v,-v′,ω,-ω′)分别计算,之后对物体表面的压强进行积分,利用式(13)~式(16)即可得出。

2 流场数值模型

本文的数值模拟包含水和空气2种流体,流动速度远远小于声速,所以认为水和空气均为不可压缩流体,把流动当作两相流,自由面为水与空气的交界面,整个过程是一个非定常过程,为了能够精确地跟踪描述水空自由面,采用应用广泛的流体体积(Volume of Fluid,VOF)方法处理。

2.1 流动控制方程[18]

含相界面的气-液两相流可用如式(17)和式(18)所示的Navier-Stokes方程统一描述:

式中:U=(ux,uy,uz)为流体质点在x、yz轴3个方向的速度矢量;P为静压;ρμ分别为流体的密度和动力黏性系数;g为重力加速度;FSV为表面张力的等价体积力形式,借助于流体体积函数F获得。

2.2 VOF模型[19]

定义一个流体体积函数F,使F的值等于一个单元内流体体积与该单元体积之比。如果F=1,说明该单元内充满流体;如果F=0,说明该单元内不含流体;如果0<F<1,说明该单元内含有自由界面。F由如下守恒方程决定:

式中:下标q代表流体属性,q=1代表空气,q=2代表水;SFq代表与流体属性相关的源项。本文的计算中,方程右端的源项为零,于是方程可写为

约束方程为

2.3 湍流模型

考虑水的黏性,利用稳定性较好的湍流模型“标准k-ε模型”对控制方程进行封闭[18, 19],其模型方程如下。

湍流动能方程k:

扩散方程ε:

式中:k为紊流脉动动能;ui为流体i方向的速度;xixj为互为相交的2个方向;ε为紊流脉动动能的耗散率;μt为流体的湍动黏性系数;Gk为由层流速度梯度而产生的湍流动能;Gb为由浮力产生的湍流动能;C1εC2ε为经验常数;σkσε分别为k方程和ε方程的湍流Pr(Prandtl)数;SkSε为用户定义的源项。

2.4 动网格方法

数值计算中,整个计算域由网格组成,物体的运动意味着物面边界的移动,即运动条件需要定义在边界条件上,使用自定义函数(User-Defined Function,UDF)对其进行描述,计算过程中被求解器动态地加载,本文中设定的运动规律(v,v′,ω,ω′)即通过UDF实现。

随着边界的运动,部分网格被不断压缩,部分网格被不断拉伸,压缩和拉伸产生的畸变会影响计算的精度,畸变达到一定程度还会产生负体积网格而导致计算失败。物体的跨介质运动过程在数值环境中就是边界的运动,所以计算过程中需要用到动网格方法,包括smooth、layering和remesh等多种网格更新方法[20],下一时间步的运动情况是由当前时间步的计算结果确定的,各种时间步的体网格的更新是基于边界条件新的位置。

smooth根据边界点上的已知位移来光滑调整流域内节点的位置,不改变网格的拓扑结构,能保证网格质量,适用于小位移问题,对于大位移和强切变问题仅靠节点松弛不能保证网格质量,而且有可能出现网格相交的情况,网格的品质会发生改变,需要用到remesh,通过设置网格控制参数如网格的最小、最大长度,网格单元的偏斜度、尺度函数参数可以控制网格移动后的品质。

3 结果验证及分析

考虑到无穷域球体的附加质量存在理论解析解,也有比较可靠的试验资料,故选用球体作为参考外形。如图 2所示,整体计算域为20 m×10 m×10 m的长方体,球体直径D=1 m,非结构网格能较好地适应球体曲面的变化,故模型的面网格采用三角网格,计算域采用四面体、六面体和楔形体;为了兼顾模型曲面网格的精细度与总体计算网格数量,将整个计算域划分为静止区域和变化区域,划分网格时在质量优先的准则下对静止区域使用六面体、楔形体网格,对变化区域使用四面体网格,在计算过程中,静止区域的网格保持不变,变化区域的网格根据边界的运动进行实时更新,模型壁面网格尺寸为0.03 m,外围边界网格尺寸为0.3 m,最终网格数量149万,网格质量绝大部分超过0.7,网格质量及局部网格如图 3图 4所示。

图 2 计算域示意图 Fig. 2 Schematic diagram of computational region
图 3 网格质量分布 Fig. 3 Mesh quality distribution
图 4 局部网格示意图 Fig. 4 Schematic diagram of regional mesh

选用压力-速度耦合求解器和couple算法,计算域上表面为压力出口,其他表面为壁面,时间步长Δt=0.001 s,计算无穷域球体附加质量系数时浸深位置为自由面y=2 m处。4组运动规律下的求解结果如表 1所示,(v,v′,-v′)代表球体分别以(v,v′)和(v,-v′)进行运动。

表 1中可以看出,y轴方向与理论值的最大相对误差为0.700 02%,z轴方向最大相对误差为0.543 39%,均与理论值吻合较好,说明选定的速度与加速度不影响计算结果的正确性。

表 1 无穷域球体附加质量系数 Table 1 Added mass coefficient of sphere in infinite region
计算条件 (v,v′,-v′)理论值y轴方向y轴方向相对误差/%z轴方向z轴方向相对误差/%
(1,3,-3)0.50.496 520.696 240.497 330.533 84
(1,10,-10)0.50.496 500.700 020.497 280.543 39
(2,3,-3)0.50.496 640.672 060.497 360.528 48
(2,10,-10)0.50.496 760.647 240.497 440.511 22

选择运动参数(v,v′,-v′)=(1,10,-10),计算球体垂直出水过程中不同浸深时y轴方向、z轴方向的附加质量系数,结果如图 5所示。从图 5中可以发现,球体顶端在接触自由面时,yz轴方向附加质量系数已经分别下降了29.6%和16.2%,符合实际的运动过程,说明方法适合于带自由面的附加质量求解。

图 5 球体附加质量系数曲线 Fig. 5 Added mass coefficient curves of sphere
4 跨介质过程中的附加质量

物体在跨介质运动过程中,由于出水阶段和入水阶段受自由水面影响,所以附加质量是随时间变化的参量,为研究其变化规律,以某简化球头柱体外形为例,头部为半球(r=0.2 m),中间为圆柱(r=0.2 m,h=1 m),尾部为圆台(r1=0.2 m,r2=0.1 m,h=0.3 m),计算域与网格划分结果见图 6图 7,出水阶段和入水阶段使用同一套网格,在数值计算时只需改变重力加速度的方向即可。

图 6 柱体外形计算域 Fig. 6 Computational region of cylinder shape
图 7 柱体外形及其局部网格 Fig. 7 Cylinder shape and its regional mesh

为便于比较,将浸深和附加质量无量纲化,即

式中:H为物体浸没于水下部分的长度;λ为附加质量;V为物体的体积。浸深为1代表刚好完全浸没水下,浸深为0代表刚好完全露出水面,浸深由0到1就代表了跨介质过程的入水阶段,由1到0就代表了跨介质过程的出水阶段。

选用组合(v,v′,ω,ω′)=(1,5,0.1,1)及(v,-v′,ω,-ω′)=(1,-5,0.1,-1),参照球体附加质量的设置条件进行计算,球头柱体出入水的附加质量系数求解结果如表 2表 3所示,图 8~图 11为结果对比图。

表 2 入水附加质量系数 Table 2 Added mass coefficient of entry step
水面/m无量纲 浸深附加质量系数
λ11λ66λ22λ26
y=-0.70.933 330.084 930.046 380.859 720.007 84
y=-0.50.800 000.075 050.036 760.724 750.042 74
y=-0.20.600 000.073 840.027 870.482 600.074 58
y=00.466 670.072 560.026 860.329 700.071 30
y=0.20.333 330.070 340.022 270.189 180.051 89
y=0.50.133 330.059 000.005 510.033 980.010 88
y=0.60.066 670.038 430.002 320.009 050.002 58
表 3 出水附加质量系数 Table 3 Added mass coefficient of separation step
水面/m无量纲 浸深附加质量系数
λ11λ66λ22λ26
y=-0.70.066 670.018 640.006 960.010 260.004 27
y=-0.50.200 000.042 000.016 370.070 410.027 28
y=-0.20.400 000.054 000.035 170.235 440.072 06
y=00.533 330.056 150.040 150.378 700.092 40
y=0.20.666 670.057 400.041 270.533 490.096 16
y=0.50.866 670.058 590.050 140.776 900.064 62
y=0.60.933 330.060 380.055 650.852 790.046 32

在本文的计算条件下,忽略物体小位移的影响,进行了6步计算,后续处理可得出2步、3步、4步、5步和6步的平均结果,不同步数平均结果的差异不超过3%,以6步平均结果为代表,直接给出其对比结果。从图 8~图 11中容易看出,无论出水还是入水,1步计算结果与6步计算的平均结果吻合较好,从曲线上看比较一致,原因是选择的外形相对简单,网格质量较高,计算结果比较理想。就实际数据而言,1步计算结果与6步计算的平均结果是有细微差别的,这主要源自于数值模型和多步计算。

图 8 随浸深变化的轴向附加质量系数 Fig. 8 Depth-varying coefficients of axial added mass
图 9 随浸深变化的附加转动惯量系数 Fig. 9 Depth-varying coefficients of added inertia moment
图 10 随浸深变化的法向附加质量系数 Fig. 10 Depth-varying coefficients of normal added mass
图 11 随浸深变化的附加静矩系数 Fig. 11 Depth-varying coefficients of added static moment

1步计算时物体的浸深反映了真实的自由面,所以严格意义上讲,只有1步计算的结果最能符合物理实际,但1步计算时流场发展不够充分,可能会导致结果偏离真实值,就工程应用而言,还需要考虑外形复杂程度、网格质量以及参数选取等因素对结果的影响。多步计算能够平抑结果的波动,但随着计算步数的增多,物体的浸深位置会发生变化,变化量的大小取决于设定的时间步长和运动速度,从而导致计算误差增大。只要选取适当的计算条件,计算结果就能够反映出实际的运动情况,而且能够降低结果对影响因素的敏感度,所以应用时可以根据情况选择使用。

跨介质的出水和入水过程中,轴向附加质量系数、法向附加质量系数、附加转动惯量系数和附加静距系数随浸深的变化曲线相似,除法向附加质量系数曲线外,其他3个均在浸深0.5附近存在趋缓区域,这与球头柱体的细长外形有关,在出水阶段和入水阶段,浸没于自由水面下的部分随球头柱体运动的规律比较接近,附加质量受此影响从而表现出趋势的一致性。从量级来看,轴向附加质量系数、附加静距系数最大均不超过0.1,附加转动惯量系数最大不超过0.06,相对法向附加质量系数而言比较小,这是因为球头柱体在全浸没状态下时法向的迎水面积远大于轴向的迎水面积,也说明其法向速度在水流体的作用下相对更难改变。

对比出水和入水过程,同样的浸深下,入水时柱体的轴向附加质量系数高于出水,最小超出29%以上,入水时的法向附加质量系数略高于出水时的数据;出水时的附加转动惯量系数与附加静距系数均高于入水时的数据。虽然曲线形状相似,随浸深变化的趋势一致,但数据差别较大,这也从附加质量的角度证明了出水不是入水的逆过程。

5 结 论

本文综合运用动网格、两相流等数值方法,针对跨介质运动物体出水过程和入水过程的时变附加质量问题,提出了一种快速有效的计算策略,推导了附加质量的计算公式并进行了验证分析,以球头柱体外形为例,对其垂直跨越自由面过程中的附加质量随浸深变化的规律进行了研究,结果表明:

1) 新的计算策略能够对跨介质运动物体的时变附加质量进行计算,精度符合工程要求。

2) 跨介质运动物体的出水过程和入水过程的附加质量变化趋势相同,规律相近。

研究结论可应用于跨介质运动体的设计、性能分析预报,对其他带自由面问题如水下发射、空投入水兵器的研究也具有一定的参考价值。限于篇幅,本文未考虑波浪因素,而且运动限定在三自由度的纵向平面内,波浪、多自由度耦合对跨介质运动物体的附加质量产生何种影响是后续的重点研究方向。

参考文献
[1] 裴譞,张宇文, 李闻白,等.跨介质飞行器气/水两相弹道仿真研究[J].工程力学,2010,27(8):223-228. PEI X,ZHANG Y W,LI W B,et al.Simulation and analysis on the gas/water two-phase ballistics of trans-media aircraft[J].Engineering Mechanics,2010,27(8):223-228(in Chinese).
Cited By in Cnki (3) | Click to display the text
[2] 刘伟. 潜水飞机总体设计与气动外形结构设计分析[D].南昌:南昌航空大学,2012:1-5. LIU W.The diving plane design and the aerodynamic shape structural design and analysis[D].Nanchang:Nanchang Hangkong University,2012:1-5(in Chinese).
Cited By in Cnki (2)
[3] 李金洪,杨安强, 粟凌云.跨介质UAV水面滑跳转向特性建模与仿真[J].鱼雷技术,2012,20(6):401-406. LI J H,YANG A Q,SU L Y.Modeling and simulation of dynamic characteristics of slide jump steering trajectory for trans-media UAV[J].Torpedo Technology,2012,20(6):401-406(in Chinese).
Cited By in Cnki (1)
[4] 王伟,张宇文, 朱灼.跨介质飞行器弹道仿真分析[J].计算机仿真,2011,28(12):1-4. WANG W,ZHANG Y W,ZHU Z.Simulation and analysis on ballistic trajectory of trans-media aircraft[J].Computer Simulation,2011,28(12):1-4(in Chinese).
Cited By in Cnki (3) | Click to display the text
[5] 王永虎,石秀华, 李文哲,等.斜入水高速冲击的理论建模及缓冲分析[J].机械科学与技术,2008,27(6):766-769. WANG Y H,SHI X H,LI W Z,et al.Modeling and cushioning analysis of oblique water entry with high velocity[J].Mechanical Science and Technology for Aero-space Engineering,2008,27(6):766-769(in Chinese).
Cited By in Cnki (7) | Click to display the text
[6] 廖剑晖,由小川, 吕海波,等.发展时变附加质量方法模拟飞行器出水过程[J].工程力学,2012,29(4):202-209. LIAO J H,YOU X C,LV H B,et al.Development of a time-varying added mass method in the simulations of the water-exit process of underwater vehicles[J].Engineering Mechanics,2012,29(4):202-209(in Chinese).
Cited By in Cnki (3) | Click to display the text
[7] 陈玮琪. 智能辨识技术及其在物体出水水动力参数辨识中的应用研究[D].无锡:江南大学,2008:1-23. CHEN W Q.Study on intelligent identification for hydrodynamic parameters of water-exit body[D].Wuxi:Jiangnan University,2008:1-23(in Chinese).
Cited By in Cnki (6) | Click to display the text
[8] 徐亮,边宇枢, 宗光华.水下机器人路径控制与仿真[J].北京航空航天大学学报,2005,31(2):162-166. XU L,BIAN Y S,ZONG G H.Path tracking control and simulation of underwater vehicle[J].Journal of Beijing University of Aeronautics and Astronautics,2005,31(2):162-166(in Chinese).
Cited By in Cnki (3) | Click to display the text
[9] 王恭义. 物体在带自由表面的流体中作任意运动的附加质量[D].大连:大连理工大学,2007:1-29. WANG G Y.Added mass for bodies moving arbitrarily in fluid with free surface[D].Dalian:Dalian University of Technology,2007:1-29(in Chinese).
Cited By in Cnki (5)
[10] SCHOEFFMANN K, TASCHWER M, BOESZOERMENYI L.Numerical simulation of transient free surface flows using a moving mesh technique[J].Journal of Applied Mechanics,2006,73(6):1017-1025.
Click to display the text
[11] YU G X,AVITAL E, WILLIAMS J,et al.Large eddy simulation of flow past free surface piercing circular cylinders[J].Journal of Fluids Engineering,2008,130(10):1135-1150.
[12] WANG K, JIN S,LIU G.Numerical modeling of free-surface flows with bottom and surface-layer pressure treatment[J].Journal of Hydrodynamics,2009,21(3):352-359.
Click to display the text
[13] VALENTIN D, PRESAS A,EGUSQUIZA E,et al.Experimental study on the added mass and damping of a disk submerged in a partially fluid-filled tank with small radial confinement[J].Journal of Fluids and Structures,2014,50(1):1-17.
Click to display the text
[14] JAIMAN R K. Added mass and aeroelastic stability of a flexible plate interacting with mean flow in a confined channel[J].Journal of Applied Mechanics,2014,81(4):1154-1159.
Click to display the text
[15] 朱仁传,郭海强, 缪国平,等.一种基于CFD理论船舶附加质量与阻尼的计算方法[J].上海交通大学学报,2009,43(2):198-203. ZHU R C,GUO H Q,MIAO G P,et al.A Computational method for evaluation of added mass and damping of ship based on CFD theory[J].Journal of Shanghai Jiao Tong University,2009,43(2):198-203(in Chinese).
Cited By in Cnki (40) | Click to display the text
[16] 李刚,段文洋, 郭志彬.复杂构型潜器附加质量的研究[J].哈尔滨工业大学学报,2010,42(7):1145-1148. LI G,DUAN W Y,GUO Z B.Added mass of submerged vehicles with complex shape[J].Journal of Harbin Institute of Technology,2010,42(7):1145-1148(in Chinese).
Cited By in Cnki (9)
[17] 傅慧萍,李杰. 附加质量CFD计算方法研究[J].哈尔滨工程大学学报,2011,32(2):148-152. FU H P,LI J.Numerical studies of added mass based on the CFD method[J].Journal of Harbin Engineering University,2011,32(2):148-152(in Chinese).
Cited By in Cnki (6) | Click to display the text
[18] 李会雄,邓晟, 赵建福,等.LEVEL SET输运方程的求解方法及其对气-液两相流运动界面数值模拟的影响[J].核动力工程,2005,26(3):242-248. LI H X,DENG S,ZHAO J F,et al.Numerical simulation of interface movement in gas-liquid two-phase flows with level set method[J].Nuclear Power Engineering,2005,26(3):242-248(in Chinese).
Cited By in Cnki
[19] CHEN X F, LI Y C,WANG Y X,et al.Numerical simulation of wave interaction with perforated caisson breakwaters[J].China Ocean Engineering,2003,17(1):33-43.
Click to display the text
[20] 褚江. 非结构动网格生成方法研究[D].南京:南京理工大学,2006:1-5. CHU J.Research of the generation of dynamic unstructured meshes[D].Nanjing:Nanjing University of Science and Technology,2006:1-5(in Chinese).
Cited By in Cnki (27)
http://dx.doi.org/10.13700/j.bh.1001-5965.2015.0478
北京航空航天大学主办。
0

文章信息

张晓强, 冯金富, 吝科, 胡俊华
ZHANG Xiaoqiang, FENG Jinfu, LIN Ke, HU Junhua
跨介质运动物体的附加质量
Added mass of trans-media moving object
北京航空航天大学学报, 2016, 42(4): 821-828
Journal of Beijing University of Aeronautics and Astronsutics, 2016, 42(4): 821-828.
http://dx.doi.org/10.13700/j.bh.1001-5965.2015.0478

文章历史

收稿日期: 2015-07-16
录用日期: 2015-11-25
网络出版时间: 2015-12-24 09:47

相关文章

工作空间