适用于Mie-Grüneisen状态方程的一种通用激波限制器

为了求解含Mie-Grüneisen状态方程多介质可压缩流动问题,针对界面追踪方法,发展了一种通用的激波限制器。根据平均守恒误差及L1误差随参数的变化情况,确定了参数的范围。由于表达式不含有状态方程中的物性参数,所以它适用于任何一种形式的状态方程。通过具体算例分析了激波限制器的开关闭状态。数值结果表明文中构造的激波限制器可以用于计算含强激波的一般状态方程的多介质可压缩流动问题。

Mie-Grüneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method（PPM）方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grüneisen状态方程。采用双渡近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡．而且数值模拟结果表明界面数值耗散也被控制在2～3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激渡等物理问题。且具有耗散小精度高的特点。

基于Mie-Grüneisen状态方程的水下爆炸数值模拟

由于Mie-Grüneisen状态方程形式比较复杂,给界面的处理带来很大难度。本文通过对Euler方程做分离变量和引入质量分数,完成了Mie-Grüneisen状态方程下的多介质可压缩流动的数值模拟,使计算过程得到简化,并通过算例验证了该方法的可靠性。在多项式形式的Mie-Grüneisen状态方程下,利用该方法模拟了球形炸药在水中爆炸后,气体和水相互作用的近场情况。在计算模型中引入的气体质量分数,很好反映了流场中不同区域内气体、水及水气并存3种不同的状态。

Grüneisen状态方程的数值解及其应用

从Grüneisen状态方程、冲击绝热关系出发,给出求解相关微分方程组的Runge-Kutta法格式.研究钨材料的内部压力、Grüneisen系数、冷压、冷能与比体积之间的数值关系,得出Grüneisen系数、冷压、冷能及常温下压力随比体积的变化曲线.计算给出的压力-比体积关系同实验数据符合较好.在此基础上给出两个压缩状态(比体积分别为V1,V2)下压力随温度变化的曲线,说明温度的变化对压力的影响小于比体积变化的影响.

基于分子动力学方法的复杂合金Grüneisen状态方程参数计算及验证:以GH4169合金为例

针对复杂合金缺乏相应势函数的情况,以GH4169为例,提出了一种计算复杂合金Hugoniot参数C_(0)和λ的方法。采用分子动力学模拟技术,研究了GH4169合金的Hugoniot曲线、0 K等温线和Mie-Grüneisen状态方程。将合金中各相Hugoniot参数按质量百分比进行加权取值得到复杂合金Hugoniot参数。通过与现有文献实验数据对比,发现该方法计算得到的Hugoniot参数、Hugoniot曲线及0 K等温线均与已有数据吻合良好,证明了所提方法的可行性。

Grüneisen状态方程中所用假设的讨论

通过分析指出,从严格意义上说,Grüneisen状态方程成立的前提是等容条件下,格林乃森参数γ与温度无关.准确地说,Grüneisen状态方程中的γ应改写为γV.为了详细说明γ和γV之间的区别,该文从晶格动力学出发,给出了这两个量的物理意义.

CO_(2)状态方程适应性对比研究

气体状态方程被广泛应用于计算CO_(2)的各项物性参数,但CO_(2)物性参数对温度和压力较敏感,不同状态方程计算的结果有偏差。针对目前常用的状态方程开展适应性研究,基于不同温度和压力条件下CO_(2)物性参数实验值对比分析了不同状态方程的计算结果。结果表明:1)GERG对CO_(2)密度及压缩因子具有较高的计算精度。计算CO_(2)密度时蒸汽区平均相对偏差小于0.5%,气相区小于1%;2)SW方程计算CO_(2)定压比热的平均相对偏差为0.29%,计算焦汤系数的平均相对偏差为0.87%;3)Fenghour法计算的黏度最大相对偏差不超过0.5%,平均相对偏差为0.34%,计算结果误差较小;4)V-W法对CO_(2)导热系数的计算结果与实验值较吻合,其平均相对偏差仅为1.24%。研究结果为油气田开发过程中涉及的CO_(2)物性参数计算提供了一定的参考依据。

运用克拉珀龙方程解决理想气体状态转换和变质量问题

理想气体状态方程p1V1/T1=p2V2/T2只适用于定质量的理想气体,在解决变质量问题时,经常需要通过巧妙选取研究对象将“变质量”问题转化为“定质量”问题,而克拉珀龙方程pV=m MRT则可直接用于变质量的理想气体.根据这一方程式,可以得到如下三个推论。

不同状态方程对于氢热物性预测的对比研究

从状态方程的角度出发,针对不同状态方程在低温高压范围内对于氢热力学性质的预测展开研究,利用目前使用较多的立方型以及多参数状态方程对于氢的PVT数据、等容热容、等压热容、焓、熵等热力学性质进行计算。将计算结果与NIST数据库给出的标准值进行比较,并对方程参数进行修正,进一步探讨了各个压力段不同状态方程的优劣以及误差产生的原因,为之后状态方程的使用、改进以及液氢相关设备的设计提出了建议。

四次型状态方程计算若干物质的固-液-气三相物性

为考察四次型状态方程(CCOR-PT EOS)的性能,对苯酚、乙烷、正壬烷、丙酮、乙酸乙酯、二氧化碳、甲苯、正庚烷、乙炔、正己烷、辛烷11个物质的固-液-气三相压力、体积、温度(PVT)物性,使用CCOR-PT EOS进行计算,并与文献数据相比较。同时考察了物质的固有参数对计算结果的影响。结果表明,该方程能够较好地描述包括固体升华压力、固体密度等在内的三相PVT物性,计算结果对临界压缩因子ZC的敏感程度大于偏心因子ω。

乙酸乙烯与乙烯混合物亥姆霍兹自由能状态方程及应用

为向高压聚乙酸乙烯生产过程提供全面可靠的热力学性质数据及预测模型,在搜集文献中热力学性质实验数据及添加合理趋势限制条件的基础上,采用Powell直接拟合算法开发了乙酸乙烯及乙酸乙烯与乙烯混合物的多参数亥姆霍兹自由能经验状态方程。结果表明:乙酸乙烯状态方程适用范围为三相点温度到650 K,压力可外推至100 MPa;乙酸乙烯与乙烯混合物状态方程适用温度为104~520 K,压力为0.1~200 MPa。经与实验数据比较分析,乙酸乙烯状态方程计算密度的总体偏差为0.3%;混合物状态方程计算压力50~195 MPa内密度的偏差在±4%内,在高压区域内计算精确度表现良好;乙酸乙烯与乙烯混合物的压缩因子与实验数据吻合较好,偏差总体分布在±3%。分析表明所开发状态方程可满足工程设计需求。

基于立方型气体状态方程的激波风洞准一维流动数值研究

采用基于立方型气体状态方程的准一维流动数值模拟方法研究了反射式高焓激波风洞的真实气体流动,重点关注了高压真实气体效应对风洞全场流动时空结构和驻室区气流参数的影响,并以理论分析揭示了高压真实气体效应对激波管内流动的作用机理。研究表明:对于以冷高压气体驱动的激波风洞,使用考虑分子体积和分子间作用力的真实气体状态方程能够更准确地描述气体的状态和风洞内的流动状况。高压真实气体效应主要在冷驱动气体中发生作用,其作用效果主要是使当地声速增大,从而使得入射稀疏波和反射稀疏波的传播速度加快;另一方面,高压气体效应在高温气体效应较显著的被驱动气体中作用微弱,且对激波管产生激波的强度和激波后的流动状态影响甚微。稀疏波的加快传播改变了激波管波系的相干时空关系。提前抵达的稀疏波可在一定情况下侵蚀激波风洞的有效试验时间。对于所测试的激波风洞构型,在150 MPa氢气驱动110 kPa氮气的工况下,高压效应导致的有效试验时间缩短约38%。适当加长驱动段长度和采用高温气体驱动均可有效减弱高压真实气体效应的影响。

一种新爆轰产物状态方程及其在炸药爆轰性能预测上的应用

为准确描述爆轰产物在高温高压下的热力学关系,拓展应用范围,实现CHNO单质/混合炸药与含能金属盐爆轰性能的可靠预测,基于Exp⁃6势的维里系数理论值建立了新气态爆轰产物状态方程Virial⁃Peng⁃Long(VPL),并基于更准确的“冷压”项与考虑金属高压热运动变化的“晶格振动”项建立了新凝聚态金属产物三项式状态方程Wu⁃Chen⁃Peng(WCP)。应用VPL计算了几种典型CHNO单质/混合炸药的爆轰参数,应用VPL与WCP评价了几种典型含能金属盐的爆轰CJ参数与驱动做功能力。结果表明,相比VLW和BKW,基于VPL能够更准确评价CHNO单质/混合炸药的爆轰性能,对太安(PETN)爆速预测偏差在±2.1%以内,最高不超过2.5%;对RHT⁃901驱动圆筒稳定速度预测误差绝对值<1%。而VPL与WCP能够准确预测含能金属盐爆轰性能,其中对叠氮化铅爆轰CJ参数计算相对误差不超过±4%;对铜叠氮化物驱动飞片速度预测相对误差绝对值小于1%。

山区输电线路跳线风偏状态方程及临界风速研究

为研究山区微地形中的竖向风速对跳线风偏的影响,基于刚性直棒法对导线、跳线及绝缘子串增加竖向风荷载分析,确定跳线处最大风偏位移值作为风偏指标,建立风荷载作用前后跳线静力平衡平面的几何变换关系,提出了考虑竖向风速影响的跳线风偏状态方程。建立半径5 km的地形模型,包含跳线风偏事故地点与微气象站点,采用数值模拟方法将微气象站点风速推演到事故地点风速,获得了考虑微地形影响的跳线风偏闪络事故风速。将事故线路参数及竖向风速代入跳线风偏状态方程,得到跳线风偏临界水平风速,通过与事故风速进行比较验证了跳线风偏状态方程的可靠性。最后根据风偏状态方程分析竖向风速的变化对跳线风偏临界水平风速的影响。结果表明:随着竖向风速的增大,跳线风偏临界水平风速线性减小,当风攻角为20°时临界水平风速可降低16%。考虑竖向风速影响的跳线风偏临界水平风速取值较为安全,建议在山区输电线路跳线设计过程中加以考虑。

一种基于Virial理论的爆轰产物状态方程

爆轰参数理论计算热化学代码的核心在于选择合适的状态方程用于描述物质的热力学状态。为扩大状态方程适用范围、减少状态方程的经验参数,基于维里(Virial)理论建立了能够准确计算到第i个Virial系数的VLI状态方程。采用VLI状态方程研究了第2~9维里系数随温度的变化规律,对比了不同阶数维里系数的相对大小及其对氮气热力学状态曲线的影响程度;计算了氮气的高压Hugoniot曲线和低压等温线,并与VLW和JCZ两种状态方程的计算值以及相应的试验数据进行对比,进一步计算了典型炸药的爆轰性能,验证了VLI状态方程的可行性和适用性。

状态方程参数化耦合制冷工质热导率性质预测模型

针对准确的制冷工质热导率性质有益于制冷热泵系统的设计和优化,引入标准化参数化方案,得到扰动链统计缔合流体理论状态方程(PC-SAFT)的输入参数,结合剩余熵标度理论(RES),预测了包括氢氟烯烃(HFOs)、氢氟烃(HFCs)、烷烃(HCs)和CO_(2)(R744)在内的11种常见制冷工质的热导率性质,使用自助法,获得了标准化参数化方案与常规参数化方案输入参数的不确定度,并展开不确定度分析。通过与文献实验结果及已报道热导率性质模型比对,所提热导率性质模型与实验数据的总体平均相对偏差为3.32%。结果表明:所提热导率性质预测模型的整体表现优于常规参数化热导率模型和CE参考项热导率模型,且标准化参数化方案有潜力为混合工质体系的研发和推广提供可靠的基础物性数据;两种参数化方案输入参数在不确定度范围内波动时,对气相迁移性质的预测影响程度相似,但标准化参数化方案在液相迁移性质的预测上表现更优;此外,基于标准化参数化方案预测迁移性质需特别关注偏心因子的精度。

RDX基含铝炸药JWL-Miller状态方程标定及对舰船结构的毁伤特性

为了精确计算某RDX基含铝炸药近场水下爆炸对舰船结构的毁伤特性,提出了一种基于水下爆炸试验和OPTIMUS的炸药JWL-Miller状态方程参数标定优化方法,建立了该型炸药水下爆炸载荷计算模型及对舰船结构毁伤的有限元模型。通过任意拉格朗日-欧拉方法(ALE)和引入JWL-Miller模型数值模拟研究不同爆距、不同攻角工况对舰船结构的毁伤特性,分析冲击波载荷和气泡载荷传播过程,计算得到爆距、攻角对局部结构毁伤影响规律,给出了断裂部位破口尺寸和外板最大凹陷距离。结果表明,运用优化得到的炸药状态方程仿真计算的冲击波压力峰值、冲击波能量结果与试验结果吻合较好;船体结构随着气泡的膨胀和收缩会呈现先中拱弯曲、后中垂弯曲的运动变形过程,随着近距爆距的改变,结构的破坏模式逐渐从冲塞形破口、四周花瓣开裂向长条形撕裂状破口再向外板向内凹陷改变;随着近距攻角的增加,舰船舷侧外板边界处剪应力在冲击波和气泡脉动载荷作用下达到应力衡准,破口面积进一步增加,气泡的溃灭速度随药包与结构之间攻角的增大而加快。

基于VLW状态方程的氢氟爆轰参数理论分析及数值模拟研究

氟气是具有强氧化性的气体,与氢气会发生剧烈的化学反应,为了得到氢氟爆轰的爆轰参数,基于VLW状态方程与爆轰理论对氟气与氢气在标准状况预混条件下爆轰性能进行了理论计算,当氢气与氟气的摩尔比为1∶1时,计算得到爆压为15.31 MPa,爆温为6 479 K。利用materials studio与CP/2K分子动力学模拟计算了氢与氟的反应,对现有氢氟反应机理进行了完善。最后通过FLUENT流体动力学仿真软件对摩尔比为1∶1时的氢氟爆轰进行了仿真计算,因爆轰波在管道内多次叠加,所以计算出超压可达46 MPa,高于理论爆压值的15.31 MPa;管道内计算出的温度可达6 800 K,与理论爆温值相差4.6%,对研究氢氟爆轰条件与能量释放规律具有一定借鉴意义。

聚合物pVT关系状态方程

聚合物的比体积(V)随压力(p)和温度(T)的变化过程主要通过聚合物的压力-比体积-温度(pVT)关系状态方程(EOS)进行描述。聚合物pVT关系状态方程可用作计算一定温度和压力条件下,聚合物的比体积,是研究聚合物相态演变及物质掺混分离规律的重要依据,也是高分子物理和聚合物加工模拟计算必不可少的模型之一。采用理论/半理论和经验方法建立了不同类型及特点的聚合物pVT关系状态方程,根据应用的不同需求而存在显著差别。随着高分子领域数字化智能化的发展,研究并建立了具有适用性和高预测精度的聚合物pVT关系状态方程,对新聚合物材料的开发和新型聚合物加工方法的建立具有重要意义。

范德瓦尔斯气体的状态方程与多方过程

范德瓦尔斯气体模型是对理想气体模型的一种修正,很多普通物理教材和文献对范德瓦尔斯气体的状态方程进行了讨论,但对过程方程很少讨论,尤其是结合物理图像的讨论。本文结合数学表达式和相关物理图像首先讨论了范德瓦尔斯气体的状态方程,并代入实际气体进行分析,再探究特殊的热力学过程,进而讨论一般的热力学过程,并与理想气体进行比较。讨论在热力学过程中功、内能、热量、热容量以及多方指数的变化,并结合相关实际应用,可以加深对范德瓦尔斯气体相关物理图像的理解。