具有拓扑稳定测度的非自治动力系统的复杂性

针对具有拓扑稳定测度的非自治动力系统的初值敏感依赖性和传递属性进行研究,得到:若一个具有拓扑稳定测度的非自治动力系统能被一列具有初值敏感依赖性的非自治动力系统逼近,则其也具有初值敏感依赖性;若一个具有拓扑稳定测度的非自治动力系统能被一列具有（r,s）-传递性或弱混合性的非自治动力系统逼近,则其也具有（r,s）-传递性或弱混合性。     

《数值流形法》序

得知郑宏教授的《数值流形法》即将脱稿,激动之情难以言表.我完全同意郑教授将流形法的特点归结为“升阶”和“切割”.受郑教授的委托,我愿在此将建立流形法的历程做一简单的回忆.其实早在1985年———也就是在我完成不连续变形分析法(DDA)之前,“升阶”在我心中已完全成熟.“升阶”似乎能取得一致收敛,弥补有限元只能按能量收敛的缺憾.我没有急于发表是因为一直没能想到令我满意的解决“不连续”的办法.起初曾想过在DDA块体上引入高阶多项式,或在块体上再划网格,也就是用分片多项式来提高DDA块体的变形精度.但我对这些方案都不甚满意,因为实际计算时多项式的阶数会受到限制,引入网格又会招致有限元既有的毛病——解的精度受控于网格质量.     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

旋转射弹高速倾斜入水多相流场与弹道数值模拟

基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程，结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响，然后对不同入水角下倾斜入水过程进行分析，得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性，但会降低弹体侧向稳定性，使射弹受到的阻力系数、俯仰力矩系数变小；入水角越小，形成的空泡越不对称，由射弹运动状态的改变引起的空泡形态变化越明显，在超空泡航行阶段，弹体运动较稳定，不同角度下流体动力系数差别很小，当弹体下表面刺破空泡壁沾湿时，弹体运动状态发生较大变化，流体动力系数迅速增大，此时入水角度过小，弹体容易失稳；弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。     

细胞松弛素B阻断微核法鉴定男性个体辐射敏感性的方法初探

初步研究了男性个体辐射敏感性的鉴定方法及标准。采集50名男性志愿者的外周血，分别给予不同剂量X射线照射，采用细胞松弛素B阻断双核法测定微核率(MNF)，通过二阶多项拟合法，绘制微核剂量效应选项中心标准曲线，将个人微核剂量效应曲线与标准曲线比对后判断个体辐射敏感性。 $0.0\sim2.5 $ Gy剂量范围内，剂量效应二阶多项拟合的中心方程为(MNF=0.014 7+0.036 2D+0.023 1D 2， r=0.726)。50名志愿者中，辐射敏感的有13人，辐射抗性的有14人，基本符合正态分布。Spearman秩和相关分析结果显示，MNF在各个辐射剂量点与辐射敏感性均存在正相关,与辐射抗性呈负相关,MNF随剂量增加而增加。本研究初步建立了“以线代点”男性个体辐射敏感性鉴定方法,并发现男性外周血淋巴细胞的本底微核率与个体辐射敏感性呈正相关。     

考虑损伤的修正梯度弹性理论

梯度弹性理论在描述材料微结构起主导作用的力学行为时具有显著优势,将其与损伤理论相结合,可在材料破坏研究中考虑微结构的影响.基于修正梯度弹性理论,将应变张量、应变梯度张量和损伤变量作为Helmholtz自由能函数的状态变量,并在自然状态附近对自由能函数作Taylor展开,进而由热力学基本定律,推导出修正梯度弹性损伤理论本构方程的一般形式.编制有限元程序,模拟土样损伤局部化带的发展演化过程.结果表明,修正梯度弹性损伤理论消除了网格依赖性;损伤局部化带不是与损伤同时发生,而是在损伤发展到一定程度后再逐渐显现出来.     

h-ADAPTIVITY ANALYSIS BASED ON MULTIPLE SCALE REPRODUCING KERNEL PARTICLE METHOD

IntroductionAs one of the meshfree methods, multiple scale reproducing kernel particle method(RKPM)[1,2]bears, besides the common features of all meshfree methods, a uniqueinherent feature of multiresolution analysis by which the structural response can b…     

氢氧混合气体爆轰波的真实化学反应模型数值模拟

采用高精度的ENO格式和基于基元化学反应的真实化学反应模型求解氢氧混合气体一维爆轰波的精细结构。采用直接起爆方法得到稳定传播的爆轰波 ,计算的爆轰波阵面参数和实验相当符合。对爆轰波反应区化学反应的研究表明 ,参与反应的不同组分具有不同类型的变化特征。网格尺寸影响的研究表明 ,计算结果的精度随着网格尺寸的增加而增加 ,并能保持较好的收敛性。移动网格研究结果表明 ,网格运动速度和爆轰速度接近时 ,两者的相互作用对计算结果产生一定影响。     

建筑风压场数值模拟中压力波动的平抑

在以同位网格为基础的简单流场压力计算中,通常采用动量插值方法来平抑流场中的压力波动现象;但是对于建筑风场等复杂的钝体绕流问题,由该平抑方法得到的收敛风压场仍可能存在小幅波动。为彻底解决同位网格格式下的压力波动,除采用动量插值方法外,本文提出了在压力校正方程的界面流速中添加压力梯度差值项的方法。算例分析表明,该方法计算得到的建筑风压场完全避免了压力波动现象,风压解与试验结果吻合良好。