高能(450 kV)X射线参考辐射场附加过滤研究

利用MCNP程序对450kV X光机建立模型，计算各参考辐射质下的原级谱、平均光子能量。对ISO 4037-1中已有范围的管电压与附加过滤片厚度关系曲线进行二次拟合，确定模拟厚度范围。利用MCNP计算各附加过滤组合下的平均光子能量值，并与外推法得出的平均光子能量值进行对比，将相对偏差最小时的附加过滤作为最终确定值。从部分能量点的模拟与测量结果的对比可看出，模拟得出的平均光子能量在2%以内与PTB及中国计量院测量结果保持一致，附加过滤最大相对偏差为18%。     

火箭整流罩外气动噪声环境的大涡模拟研究

基于五阶加权本质无振荡（WENO）格式构造隐式大涡模拟方法（ILES）,对跨声速来流条件下（Ma=0.8）火箭整流罩外噪声环境进行数值模拟。通过与风洞试验结果及国外文献进行对比,ILES方法能够在较粗网格下准确预测壁面湍流脉动特性。跨声速流动在壁面折角处出现分离、再附、激波/边界层干扰现象,均方根脉动压力系数出现峰值,同时该区域噪声能量在全频段都较高,易引起结构抖振效应。最后,根据ILES模拟结果,指出工程常用的外噪声经验公式的不足,并提出改进措施。     

超声速复杂流动中湍流模型的性能评估

对8°,16°,20°,24°超声速二维压缩拐角进行了数值模拟,系统评估了SA(Spalart-Allmaras), k-ω, SST(Shear Stress Transport)3种工程常用湍流模型对激波/边界层干扰复杂流动的模拟适用能力.对于"弱"干扰、小分离,工程常用湍流模型能够较准确地预测其压力、摩擦力、热流分布,而对于"强"干扰、大分离则会造成较大的偏差.同时,随着分离区的增大,各湍流模型的适用性能差别更加明显,数值方法也有一定的影响.      

火星进入器作强迫震荡运动壁面脉动压力数值模拟

为研究超声速阶段进入器作强迫震荡运动对壁面脉动压力环境的影响规律，本文耦合进入器刚体运动方程与流体力学方程，采用动网格技术，对火星进入器模型开展非定常数值模拟，获取壁面不同位置处的脉动压力信息。研究表明：在超声速阶段，进入器作强迫震荡运动诱导的脉动压力远大于进入器保持相对静止时仅由非定常流动诱导的脉动压力。来流马赫数为1.2时，进入器作强迫震荡运动对脱体激波影响较小，脱体激波强度较弱且形态变化较小，攻角的震荡导致同一测点距离脱体激波的位置发生周期性改变，舱体迎风面及配平翼迎风面的脉动压力环境主要受攻角变化的影响；来流马赫数为3时，进入器作强迫震荡运动对脱体激波的影响较大，脱体激波震荡剧烈，诱导舱体迎风面及配平翼迎风面产生极其恶劣的脉动压力环境，功率谱分析表明激波震荡诱导的脉动压力能量主要集中在30 Hz左右。     

火星进入器壁面脉动压力环境数值模拟

为研究跨超声速阶段来流马赫数、来流攻角及配平翼展开角的变化对进入器壁面脉动压力环境的影响规律，本文采用脱体涡方法对火星进入器模型开展非定常数值模拟，获取壁面不同位置处的脉动压力信息。研究表明：在跨超声速阶段，进入器壁面脉动压力环境随马赫数的增加而趋于减缓。配平翼迎风面分离区受脱体激波影响明显，当来流马赫数较小时，可压缩效应较弱，分离区涡流运动剧烈，诱导的脉动压力环境较强；随着来流马赫数的增加，脱体激波对分离区抑制作用增强，分离区运动受到限制，诱导的脉动压力环境趋于平缓。此外，随着来流攻角增加，配平翼迎风面上再附点的位置向翼根方向转移，从而使翼根处的脉动压力环境趋于恶劣。当配平翼展开180°时，分离区再附点位置基本固定，配平翼迎风面脉动压力环境得到一定程度的减缓。功率谱分析表明，在配平翼迎风面上诱导的脉动压力能量主要集中在中低频区域。     

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA (Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

房地产开发规模的优化决策

针对我国房地产供求失衡的状况，从需求管理的角度，以控制房地产投资规模、保持供需平衡为目标，利用经济控制论的库存理论，建立了房地产开发规模优化决策的数学模型，推导出模型的最优解。经实例分析验证，该数学模型对于确定最优的房地产开发规模而言，具有简洁、实用的特点，对于提高房地产市场管理及投资经营决策水平具有广泛的应用价值。     

基于熵限制的Baldwin-Lomax湍流模型

针对Baldwin-Lomax(BL)模型在模拟超声速复杂流动中的不足,提出一种基于熵限制的新型修正方式.分别采用原始的BL模型(BL-origin)、修正后的BL模型(BL-entropy)以及Wilcox k-ω两方程模型对超声速平板、超声速二维压缩拐角以及膨胀-压缩拐角3个典型算例进行计算,结果表明BL-entropy 能够较好地获取长度尺度进而得到均匀合理的涡粘性分布.同时,该修正方式简单高效并且适用于其他与超声速边界层有关的模型,具有较大的发展潜力.     

比熵增概念及其在湍流模型中的应用

通过对不同马赫数下(Ma=0.7、2.25、6)平板边界层进行直接数值模拟研究,提出比熵增的概念.比熵增表征单位机械能的耗散,并且在数值上表现出较好的单调性和马赫数无关性,因而能够可靠地表征边界层的范围.应用比熵增概念修正BL零方程湍流模型长度尺度,构造BL-entropy,同时采用原始的BL模型,SA一方程模型对后台阶(Ma=0.128)以及高超声速锥柱裙组合体(Ma=7.05)流动进行数值模拟,并与实验数据进行比较,结果表明BL-entropy能够准确地获取长度尺度进而得到均匀合理的涡粘性分布,极大地提升了原始BL模型模拟复杂流动的能力.     

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-Stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA(Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.