基于车身绕流的低雷诺数湍流模型改进研究

本文将汽车绕流模块化为各典型局部流动,通过常用湍流模型对各典型局部流动进行数值模拟,结果验证了湍流模型对转捩的捕捉能力是准确模拟汽车绕流的关键. 在分析汽车绕流分离及转捩机理的基础上,优化了稳态和瞬态求解方法,改进了湍流模型对转捩的预测能力,进而提高了湍流模型在汽车流场模拟上的精度. 针对汽车绕流的稳态问题,将流线曲率因子及 响应阈值引入 LRN $k$-$\varepsilon $ (low Reynolds number $k$-$\varepsilon $) 模型,获得了一种能够更准确预 测转捩的改进低雷诺数湍流模型 (modified LRN $k$-$\varepsilon $),改善了原模型对湍流耗散率的过强依赖性及全应力发展预测不足等问题;针对汽车绕流瞬态求解,通过分析 RANS/LES 混合湍流模型的构造思想及特点,引入约束大涡模拟方法,结合本文提出的改进的 LRN $k$-$\varepsilon $ 湍流模型,提出了一种能准确捕捉转捩现象 的转捩 LRN CLES 模型. 分别将改进的模型用于某实车外流场和风振噪声仿真中,通过 Ansys Fluent 求解器计算,并将计算结果与常用湍流模型的仿真结果、HD-2 风洞试验结果和实车道路实验结果进行对比,表明改进后的湍流模型能够更准确模拟复杂实车的稳态和瞬态特性,为汽车气动特性的研究提供了可靠理论依据及有效数值解决方法.      

土体的本构模型和超重力物理模拟

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.      

高熵合金的力学性能及变形行为研究进展

高熵合金是近年来提出的一种新的合金设计理念,打破了一般合金中以1种或2种元素为主,辅以极少量其他元素来改善合金性能的传统思想,由多种元素以等原子或近似等原子比混合后形成具有独特原子结构特征的单一固溶体合金.高熵合金的多主元特性使其在变形过程中表现出多重机制(包括位错机制、形变孪生、相变等)的协同,因而高熵合金已经展示了优异的力学性能,如高强、高硬、高塑性、抗高温软化、抗辐照、耐磨等,被认为是最具有应用潜力的新型高性能金属结构材料,已经成为国际固体力学和材料科学领域研究的热点.本文首先介绍了高熵合金独特的结构特征,即具有短程有序结构和严重的晶格畸变;随后对近年来针对不同类型高熵合金(包括具有面心立方相、体心立方相、密排六方相、多相以及亚稳态高熵合金)力学性能、变形行为方面的研究成果,特别是强韧化机制以及相关的原子尺度模拟,进行了较为系统的综述;最后强调了高熵合金未来研究中所面临的一些主要问题和挑战,并对其研究进行了展望.     

基于LBM方法的风力发电机尾流结构仿真

针对传统CFD数值计算方法难以实现风力机动态旋转及其旋转状态下的流固耦合计算,本文结合格子玻尔兹曼(LBM)方法易于处理动态复杂边界的特点及大涡模拟(LES)方法在非稳态涡流结构捕捉上的优势,采用LBM-LES联合方法进行三维风力发电机整机气动性能及尾流结构仿真研究,同时采用尺度自适应方法对尾涡结构进行跟踪和精细化计算。针对NREL PhaseⅥ型试验机进行模拟,得到了与实验结果吻合的流动形态及尾流结构演变规律,分析了尾流区速度演变规律并对比了不同亚格子湍流模型对计算结果的影响.     

质子辐照硼硅酸盐玻璃盖片的物理效应分析

作为航天器电源系统的重要组成部分,太阳电池需要更高的转换效率和可靠性以及更长的使用寿命。通过在太阳电池表面覆盖抗辐照玻璃盖片,可以增强太阳电池对粒子辐射的防护,延长太阳电池的服役寿命,使航天器获得可靠的能源供应。硼硅酸盐玻璃就是一种理想的太阳电池玻璃盖片材料。采用蒙特卡罗方法,结合SRIM软件模拟研究质子辐照硼硅酸盐玻璃的损伤物理机理。基于粒子与物质相互作用的理论以及基本公式,通过分析不同入射能量的质子在硼硅酸盐玻璃中的阻止本领、电离能损、位移能损、空位的产生情况,对辐照损伤的物理机制进行研究。结果表明:能量为30~120 keV的质子辐照损伤主要发生在硼硅酸盐玻璃表面;质子沉积、空位分布等均为Bragg峰型分布;电离能损是能量损失的主要部分,随入射能量的增加而增大,导致电子的电离和激发;位移能损在玻璃内部随能量降低而增大,导致硼、氧和硅等空位缺陷的产生;电离效应和缺陷的产生是硼硅酸盐玻璃色心形成的重要原因。     

用人工神经网络方法测量生物组织光学参数γ

提出了一种利用亚扩散空间分辨漫反射预测生物组织约化散射系数μ_s'和相函数参量γ的人工神经网络方法 .采用蒙特卡罗方法得到光子经生物组织漫反射的数据样本,利用这些数据样本训练反向传播神经网络,用于从亚扩散散射光中预测γ的信息.为了解决同时预测μ_s'和γ两个参数时会产生较大误差的问题,分段数据训练两个BP网络,依次识别μ_s'和γ.研究发现3.64l_(th)(l_(th)表示平均输运自由程)是对γ的不敏感点,可用该点附近的数据样本训练网络用于预测μ_s',用2l_(th)范围内的数据样本训练网络用于预测γ.蒙特卡罗仿真结果表明,在1.3≤γ≤1.9范围内,预测结果与真实值的相对均方根误差在1%以内.与现有的测量方法相比,所提的人工神经网络方法更加简单,且提高了预测精度.     

碳离子在不同材料叶片的多叶光栅上产生的次级粒子研究

在碳离子放射治疗中，碳离子束与治疗头设备和患者身体相互作用产生的次级粒子可以到达患者体内的许多区域，在产生的次级粒子中以中子和$\gamma $射线的产额为最大。在不影响束流配送功能的情况下，减少碳离子放疗中产生的次级中子和$\gamma $射线对于降低放疗后出现的正常组织并发症及二次肿瘤风险有着非常重要的意义。本文利用蒙特卡罗(Monte Carlo)方法模拟计算了被动式束流配送系统下，400 MeV/u碳离子束照射到由不同材料叶片构成的多叶光栅(MLC)形成典型的10 cm×10 cm方形射野时，在水模体中产生的次级中子和$\gamma $射线所沉积的剂量及空间分布等。模拟结果显示:碳离子束通过MLC形成射野后在水模体中产生的次级中子主要分布在水模体的入射端，次级$\gamma $射线较为均匀的分布在整个水模体内，且较多分布在具有展宽Bragg峰(SOBP)射野在水模体中贯穿时的坪区。对于MLC叶片材料的选择，则需根据实际情况对叶片厚度以及次级粒子当量剂量的要求来确定。本文通过模拟研究不同MLC叶片材料产生次级粒子的情况，为被动式束流配送系统中MLC叶片及其他元件的材料选择提供了科学依据。     

三维格子Boltzmann方法研究上升双气泡的运动特性

应用格子Boltzmann三维模型，对双气泡在静水中的运动进行数值研究.采用八点差分和十八点差分格式分别求解一阶▽φ和二阶▽²φ可以有效避免气液密度比过大造成的数值不稳定问题.结果表明：当两个相同直径的气泡在上升时，位置靠上的气泡形状变化像单气泡上升一样，而位置靠下的气泡会受到前一个气泡尾迹的影响，并有很明显的形状变化.当两个气泡直径不同时，不管初始位置如何，大气泡总会对小气泡造成强烈的影响.     

磁化套筒惯性聚变初步探索研究与实验可行性分析北大核心CSCD

磁化套筒惯性聚变(MagLIF)构型可充分利用现有大型脉冲功率驱动装置,如聚龙一号等。基于磁流体力学方程组和1∶1比例氘氚(DT)混合燃料聚变模型,开发了零维MagLIF数值模拟程序并进行了初步探索研究。计算结果表明初始负载参数(如轴向磁场强度,预加热温度、时刻,负载半径等)与聚变产额之间有着密切的联系,在给定条件下,可依据计算给出的定性关系进行负载优化设计。值得注意的是,根据计算结果,即使在理想条件下,氘氚燃料要实现能量收支平衡,则驱动器的电流必须大于21.2 MA。这意味着聚龙一号装置(10 MA)无法开展集成化的MagLIF实验,进一步的校验计算验证了上述观点,并在此基础上提出铝套筒分解实验的建议和负载设计参数。所取得的计算结果有利于加深对MagLIF套筒压缩阶段物理过程的认知和理解。     

真空条件下LiCl溶液蒸发量随浓度变化规律研究

除湿溶液再生是除湿空调系统的核心过程,为了提高除湿溶液再生效率,以蒸发量作为评价标准,在不同的真空度下对不同浓度LiCl溶液进行再生。对不同再生条件下,溶液的再生情况进行数学模拟及实验研究,并将模拟结果与实验结果进行对比分析,结果表明,实验结果同模拟结果具有相同的变化趋势。在6 kPa和8 kPa的压强下,LiCl溶液的蒸发量随浓度的增大先上升后下降,而在1 kPa的压强下,LiCl溶液的蒸发量随浓度的增大呈现出先下降再升高接着又下降的趋势。