爆炸波对高层建筑干扰流场的数值模拟.pdfVIP

爆炸波对高层建筑干扰流场的数值模拟 刘君 陈宏 刘伟 (嗣舫科技炎学航天技术蒹， 长抄市410073) 籀要：本文从三维非定常Ⅳ一S方程出发，采用时间和空问二阶精度的^矾理)22汀和新型肋差分 格式，对爆炸波扫过高层建筑时形成的干扰汽场进行了数值模拟和理论分析． 关键诃：爆炸波，人防工程，计算模拟 1 引言 现代战争有时要以城市为战场，即使在和平环境，城市也存在由于各种可燃物质导致 爆炸事件的可能性：因此研究爆炸波对高层建筑的作用规律具有现实意义和应用价值．对 这～问题的研究除了使用激波管等实验设餐进行试验外，还可以采用数值模拟的手段．本 文从兰雏非定常守恒型Ⅳ一S方程出发，采用时间和空向二阶精度、高分辨率、高效率的 舳拼差分格式，对可燃气云在高层建筑附近爆炸时形成的干扰流场进行了数值模拟， 研究了建筑物表面受到的压力载荷分布和整体冲击力的作用规律。为提到计算精度，本文 还在舳2M差分格式基础上发展了一种新型的科0差分格式。对上述物理模型进行 重复性计算：比较发现ENO差分格式可明显提高流场中物理量在极值点附近的分辨率， 但是降低计算效率． 本文的工作对城市内高层建筑的安全工程问题有一定的指导意义。 2控制方程 采用如下无量编形式： (x，弘力=F，y+，2+)／n+，(砧，v们=(“+，v’，1I’)，口帕+ ‘P=，，疋，P=p’，【疋(Z)2】，T=T’蠢／R’ 其中上标“·”表示有量纲物理量，下标“oD”表示环境静止大气物理量；尺’为建筑物 的底部半径，Z为静止大气的音速．这样，三维非定常守恒型N—s方程可以写成： 型+翌±堡!+丛生±墨! (2．1) 魂 融 +堑掣：o 砂 院 其中U=(A肛f，肌pp，e)7，E，F和G是对流项通量．E、只和G，是粘性项通量， 口 体表示形式较为通用，在此省略·为封闭方程组还需补充完全气体状态方程：p=pT． 计算时为了准确描述物体外形，需要将上式变化到任意系(f，厶，7，力-Fi为提高计算 效率，进行薄层近似处理：具体过程可参考文献【1】．采用时间相关法把兰维问题分解为 三个～维问题求解，写成算子形式： ￡是三维问题算子：k，岛和0是一维问题算子，粘性顶用中心差分格式离散，对流 项差分格式如下． 3差分格式简介 近年来国内外对差分格式研究进步很大，其中rpD类型的格式在应用于求解包含激 波等间断的流场计算中取得实效，受到计算流体力学工作者重视。在目前看到的时间和空 间二阶精度的r国格式中。作者刘君在张涵信院士指导下提出的NND2^，格式，由于投 有涉及到矩阵运算，在求解守恒型Ⅳ一S方程时具有较高的计算效率．在文献f2】中应用 肋2M格式效值模拟冲击波扫过兰维地下通道时的流动现象，取得较为满意的结果， 本文计算依然采用这一格式． 7 现以一维层U丝职方程为例，对舳2jl彳格式进行简单介绍： 型+竺：0 ’ (3．” A=A++彳一，F=F++F一，这样一来(3．1)式变为： aUaF+OF一 一 ——+～+一=0 (3．2) 研 良 缸 方程(3．2)用NND2M格式离散为： 咿=u；一尝(H—n；一日4’ (3-3) 式中上标刀为时问