爆炸冲击载荷作用下板壳结构数值仿真分析

主要针对爆炸冲击载荷作用下板壳结构的试验破坏问题,利用LS-DYNA有限元分析软件,采用非线性动力学分析计算方法,考虑材料非线性和结构非线性等因素,模拟分析了板壳结构在接触爆炸冲击载荷作用下的动态响应。计算分析结果与试验结果相吻合,利用有限元分析方法能很方便解释试验过程和现象,为试验分析提供有效依据。     

变轨期间对接机构对接锁系温度场数值模拟

采用节点热网络法,建立了目标飞行器尾部组合系统(发动机组和对接机构)温度场数值模拟的物理、数学模型,明确了关键参数的计算方法,应用隐式R-K法,耦合求解非线性温度微分方程组.着重研究变轨期间对接锁系受高温变轨发动机的热影响.分析表明,发动机高温热辐射对邻近对接锁系的热影响很大,发动机热影响范围大致在θ±Δθ=(90±30)°和θ±Δθ=(270±30)°区间,发动机相对安装距离R、变轨起始时刻和工作时间是影响对接锁系温度的重要因素,而发动机安装倾角对对接锁系温度的影响非常小.     

水滴多尺寸分布对水滴撞击特性和结冰增长的影响

水滴直径是影响水滴撞击特性的重要参数,而撞击特性在很大程度上影响结冰冰形,因此,在结冰数值模拟中,如何合理选取水滴分布模型,对于保证计算精度有重要意义.本文基于求解雷诺时均N-S方程获得的流场结果,采用拉格朗日法计算水滴撞击特性,并采用Messinger模型进行结冰增长模拟,研究了水滴直径、水滴多尺寸分布对水滴撞击特性和结冰冰形的影响.结果表明,水滴的多尺寸分布对于结冰数值模拟结果具有重要影响,采用Langmuir D分布能够满足数值模拟精度要求.     

建筑物对污染物扩散影响的数值与风洞模拟研究

采用k-ε(RNG)与RSM湍流模型对处于方形建筑物不同位置污染源所排放污染物的扩散规律以及建筑物对流场的影响进行模拟,并且与相应的风洞试验结果进行了比较.流场分析结果表明:数值模拟能够较好地模拟建筑物前方迎风侧停滞回流、顶部回流以及后方空腔区等.浓度场分析结果表明:建筑物前方迎风侧以及顶部回流区污染物扩散的数值模拟结果与风洞试验结果基本相等,而在建筑物后方空腔区污染物的数值模拟结果略高于风洞试验结果.综合分析并与风洞试验结果相比较,RSM模型能够较好地模拟污染物浓度场以及建筑物周围流场的变化规律.     

红外成像非接触转捩测量低速风洞试验技术研究

发展红外成像非接触转捩测量低速风洞试验技术,旨在解决特殊气动布局外形及金属材料模型转捩位置测量问题.通过在模型表面产生热壁面、现场测试模型表面发射率、使用遮蔽板、在金属模型表面喷涂隔热氟碳漆等措施,解决了环境条件、发射率、辐射传递干扰、金属模型材料特性等阻碍红外成像技术应用的关键问题;通过数值计算及试验测试得到模型热壁面与环境温差在20℃范围内,热壁面背景温度对转捩位置基本没有影响,解决了热壁面对转捩位置影响问题;通过试验原理、试验方法、关键参数测试、转捩判据、准度考核等研究工作,构建了红外成像非接触转捩测量低速风洞试验技术;通过引导试验考核了试验系统.结果表明:该技术实用可靠,值得推广.     

基于CBS有限元的流动-传热-变形耦合计算方法

研究了基于CBS有限元法和常规有限元法相结合的流动-传热-变形耦合计算方法。在该方法中,流体流动和传热采用CBS有限元方法计算,固体变形采用常规的有限元方法计算,实现了流体域和固体域统一的有限元网格划分,简化了变形过程中的网格生成和不同网格间的数据交换。然后,依据此方法发展了计算程序,并通过算例分析,校验了计算方法的可行性和程序的计算能力。     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

高焓电弧风洞试验热化学非平衡流场数值模拟

针对高焓电弧风洞内部流动的热化学非平衡效应及气体组分和振动能量冻结效应导致的试验数据外推困难问题,基于高焓风洞喷管/试验段/试验模型一体化数值模拟的思路,通过数值求解三维热化学非平衡Navier-Stokes方程,开展了FD-15高焓电弧风洞典型运行状态下流场的数值模拟,与典型试验状态的气动热数据进行了对比验证,研究了试验数据外推飞行条件的方法及有效性问题,分析了提高驻室总压对试验数据外推的影响。研究表明:（1）风洞试验段来流离解度高,热化学非平衡效应及其冻结现象严重;（2）热流校核试验测量数据位于一体化数值模拟的完全催化热流和非催化热流之间,分布合理,验证了计算方法和程序的正确性;（3）试验模型安放位置对模型表面压力和热流存在影响,模型与喷管出口的距离越大,模型表面压力和热流越低;（4）当驻室总压较低时,通过双尺度模拟准则（模拟飞行条件总焓和双尺度参数ρ_∞L）外推热流失效,使用部分模拟准则（模拟飞行条件总焓和驻点压力）外推热流也会出现较大差异,在非催化条件下这一现象更加明显;（5）当驻室总压较高时,使用双尺度模拟准则或部分模拟准则外推飞行条件,产生的热流差异明显减小。     

氮气辐射强度的激波管测量与验证

超高速飞行器的热防护设计必须考虑激波层内高温气体发射与吸收的辐射能量,需要有效的辐射加热评估手段。相应飞行条件下的光谱辐射强度地面实验测量是验证数值模型和方法、理解高焓流动的重要手段。基于燃烧驱动激波管,发展辐射强度标定技术,针对富氮气环境,开展高温气体光谱辐射强度的高分辨定量化测试,掌握辐射特征,为数值验证提供基础数据。实验获得了激波速度5.70和6.20km/s条件下的气体光谱辐射强度精细结构,数据表明激波波后的非平衡过程对辐射强度存在很大影响。通过求解耦合化学反应动力学模型的Navier-Stokes方程和辐射特性模型,得到对应实验条件下的流场参数和辐射强度,计算结果和实验数据符合很好,验证了数值模拟方法。     

湍流燃烧模拟中化学反应的加速算法研究进展

为研究湍流燃烧数值模拟中化学反应机理计算的加速方法,讨论了动态自适应化学（Dynamic Adaptive Chemistry,DAC）方法和Krylov子空间近似的指数格式的应用情况。在湍流火焰大涡模拟中,使用DAC简化可以加速化学反应计算。然而,在并行燃烧数值模拟中,处理器核心的负载极度不平衡,加速效果有限。而Krylov子空间近似的指数格式的加速效果可以作用于每个处理器核心,更有利于整体计算效率的提高。在同等精度下,相比于隐式格式耦合DAC和MTS加速方法,Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著。