假设检验与CFD不确定度量化分析

针对CFD验证与确认中不确定度量化研究,提出了假设检验量化技术,采用了Grubs数据处理技术并具有好的鲁棒性,应用了假设检验原理分析了计算结果,可对CFD计算具指导作用,数值例子也说明了这一点。     

面向流场计算的专用超级计算机体系结构研究

 本文在对流体力学等一类大型计算问题的并行计算模型进行分析的基础上,提出了一种可以实现大规模并行处理的专用超级计算机体系结构方案。该系统由一台宿主机和一个多计算机阵列组成一台个人超级计算机。利用专用性,减小复杂性,提高并行性是系统的主要特色。系统采用二维超环结构作为计算机结点之间的互连结构,反映了计算中数据交换的局部性特点,具有好的通信性能。每个计算机结点本身又采用流水线向量处理结构,从而将高层次的多机并行与低层次的向量处理相结合,实现了大规模的并行计算。本文讨论了系统的软、硬件结构并研究和评价了通信性能。     

微型飞行器的仿生流体力学——昆虫前飞时的气动力和能耗

 用数值模拟方法研究了昆虫前飞时的气动力和需用功率。由N S方程的数值解提供速度场和压力场,从而得到涡量、气动力和力矩 (惯性力矩用解析方法计算 )。基于流场结构,解释了非定常气动力产生的原因;基于气动力和力矩,得到需用功率。悬停飞行中揭示出的 3个非定常高升力机制 (不失速机制,拍动初期的快速加速运动,拍动后期的快速上仰运动 )在前飞时仍然适用 (即使在快速前飞时,V∞ =2～ 2.5m/s,失速涡也不脱落 )。在低速飞行时 (V∞ ≈ 0.5m/s)平衡重量的升力既来自于翅膀的下拍运动也来自于上挥运动,并主要由翅膀的升力贡献;克服身体阻力的推力主要来自于翅膀的上挥运动,由翅膀的阻力贡献。在中等速度下 (V∞ ≈ 1.0m/s),升力主要来自于下拍运动,其中一半由翅膀升力贡献,一半由翅膀阻力贡献;推力主要来自上挥运动,也是一半由翅膀升力贡献,一半由翅膀阻力贡献。在快速飞行时 (V∞ ≈ 2.0m/s),升力主要来自于下拍运动,主要由翅膀阻力贡献;推力来自上挥运动,主要由翅膀升力贡献。悬停时,下拍和上挥做功同样大;前飞时,下拍做功较上挥大得多 :V∞ =0.5,1.0和 2.0m/s时,下拍做的功分别是上挥的 1.6,2.6和 3.5倍。     

带裂纹板在冲击载荷作用下动态应力强度因子的数值计算

 本文使用八节点等参元和六节点畸变等参元的有限元方法,计算了中心穿透裂纹板在冲击载荷作用下的动态应力强度因子,并与Chen和Aberson的计算结果作了比较,本文的结果具有较高的精度和较少的计算量。文中还给出了动载影响系数k₁（t）和静态应力强度因子的关系式。并考虑了材料特性和板的几何尺寸对于动态应力强度因子的影响。     

DYNAMIC BEHAVIOR OF A CRACKED FLEXIBLE ROTOR SUPPORTED ON JOURNAL BEARINGS

SupportedbytheNationalNaturalScienceFoundationandtheShaanxiNaturalScienceFoundationofChinaNomencIaturescd-externaldampingCP-machinedclearanceofbearingD.-externaldampingratio(=c,/(2m,eq))e.-theunbalanceeccentricityfx,f,-fluidfilmforcesk,-complexstiffnessofcrackedshaftk.-stiffnessofuncrackedshaftke,k,-stiffnessofcrackshaftin$,Vdirec-tlonsAse,As,-thelargeststiffnesschangemb,md-masslumpedatbearingandrotormid-spanS,-fiXedSommerfeldnumber(=S/n,SistheSommerfeldnumber)U-unbalanceparameter(=e./…     

The History and Dynamics of Comet 9P/Tempel 1

Since its discovery in 1867, periodic comet 9P/Tempel 1 has been observed at 10 returns to perihelion, including all its returns since 1967. The observations for the seven apparitions beginning in 1967 have been fit with an orbit that includes only radial and transverse nongravitational accelerations that model the rocket-like thrusting introduced by the outgassing of the cometary nucleus. The successful nongravitational acceleration model did not assume any change in the comet’s ability to outgas from one apparition to the next and the outgassing was assumed to reach a maximum at perihelion. The success of this model over the 1967–2003 interval suggests that the comet’s spin axis is currently stable. Rough calculations suggest that the collision of the impactor released by the Deep Impact spacecraft will not provide a noticeable perturbation on the comet’s orbit nor will any new vent that is opened as a result of the impact provide a noticeable change in the comet’s nongravitational acceleration history. The observing geometries prior to, and during, the impact will allow extensive Earth based observations to complement the in situ observations from the impactor and flyby spacecraft.     

Navier—Stokes（NS）方程组差分计算中的物理和网格尺度效应及NS方程组的简化

对于高Re数流动计算，在通常二阶精度NS差分格式和网格数条件下，存在某些粘性项落入修正微分方程截断误差项的问题。这类NS方程组计算实际是计算某种简化NS方程组，而且重复计算误差物理粘性项既浪费机时和内存，误差积累又会对数值解产生不可预测的影响，避免外述缺陷的办法一个提高NS差分格式的精度，另一个是丢掉可能落入截断误差项的物理粘性项，把NS方程组简化为广义NS方程组，广义NS计算避免了误差物理粘性项误差积累对数值解的不可知影响，又可节省内存和机时，对高Re数流体工程计算很有好处。利用广义NS方程组计算超声速绕前向和后向台阶流动的结果表明：广义NS方程组与NS方程组的数值结果很好相符符。     

制导炸弹大空域下的动力学模型

针对惯性导航/卫星定位方式制导炸弹在大空域下的工作过程和特点,建立了空中的运动数学模型,分析了大空域下制导炸弹运动的非线性因素,并在线性化模型的基础上,发展了一种时变的非线性模型来描述制导炸弹大空域下飞行动力学系统的方法。最后进行了仿真,结果表明,该方法是合理的,为后续制导炸弹控制系统的自适应设计与仿真提供了参考。     

基于CFD的电磁散射数值模拟

计算流体力学(CFD)技术中的时域有限体积法(FVTD)被推广到计算电磁学中,直接数值求解时变麦克斯韦方程组.FVTD使用电磁场矢量共置于网格单元中心、散射场积分形式的麦克斯韦方程组守恒形,空间离散使用基于特征值的近似黎曼解构建网格单元边界通量,时间推进采用四阶龙格-库塔法.TM和TE波极化下几种二维完全导电体的表面诱导电流密度和雷达散射截面(RCS)计算验证表明,时域有限体积法是一种高精度有效的时域方法.     

跪式起落架在武装直升机坠毁过程中能量吸收能力研究(Ⅰ)——数值仿真计算

 在简述直升机抗坠毁原理的基础上,模拟直升机机体的等效质量与跪式起落架构成的组合系统在 6m/ s硬着陆垂直撞击地面时的塑性动力响应和能量吸收过程。所用的模型为 :1基于真实几何构型和材料特性的起落架 FEM动力学计算模型;2简化的弹簧 -刚性杆系统模型。通过 Lagrange方程解出了直升机以6 m/ s的速度垂直撞击地面这一过程的动态响应,近似给出了起落架吸收的能量 (塑性功 )占初始动能的百分比;机体的动能变化曲线以及主缓冲器的载荷谱曲线。两者结果进行了比较,为直升机抗坠毁设计提供理论指导。