基于WENO格式的一维溃坝波的数值计算

应用高精度加权本质无振荡(WENO)格式构造了溃坝水流计算的数学模型,模拟了坝体瞬间全溃时洪水波的运动过程,并将计算值与理论解进行了比较。结果表明,WENO格式精度高,稳定性能好,能够很好地模拟溃坝波的演进过程。     

激波和爆震波在U形管中传播过程模拟

为缩短爆震室轴向尺寸,并满足爆燃波向爆震波转变(DDT)或激波向爆震波转变(SDT)距离的需求,采用曲管爆震室代替直管爆震室来研究曲管中的爆震特性。应用数值方法对激波在U形管中传播、加速过程,SDT过程以及爆震波在U形管中的传播过程进行了模拟。U形管的曲率直径和入射激波强度对SDT有很大影响,较小的曲率直径对SDT有较大的促进作用;在入射激波速度大于960 m/s,曲率直径小于78 mm时,激波可以成功向爆震波转变;爆震波在弯管段形成后,向上游和下游两个方向传播,传播的速度和压力均不同;向上游传播的爆震波可以显示出"三波点"的运动轨迹;弯曲光管对爆震波在短距离内形成是有利的,可以用来作为预爆管。     

一维溃坝波的数值模拟——基于WENO重构的解耦格式

采用五阶精度的WENO格式对界面变量进行数值重构,形成高精度的数值离散格式,对其捕捉溃坝波的能力进行分析。用roe平均对一维浅水方程组进行局部线性化,通过直接和解耦2种格式进行求解,数值模拟溃坝洪水的推进过程,得到沿程的水位和流速分布。与理论解的比较表明基于WENO5重构的解耦格式,具有很高的分辨率,展现了溃坝涌波的运动特性。     

一维WENO格式及其在标量守恒方程数值计算中的应用

研究高阶精度加权本质无振荡(WENO)格式及其在标量守恒律方程中的应用.应用五阶WENO空间离散格式和三阶TVD Runger-kutta时间离散格式对一维标量守恒方程以及二维标量守恒方程进行了数值模拟.数值结果表明该格式具有高精度性和本质无振荡性.     

一种基于WENO格式的一维溃坝水流数值模拟研究

采用一维Saint-Venant方程组,应用WENO格式和Runge-Kutta时间离散的思想,进行溃坝水流的数值模拟,得出了水位和流速的沿程分布,并与理论解比较,发现数值解在间断波附近没有出现数值振荡,水位和流速数值解与理论解吻合较好,表明WENO格式是一种进行溃坝水流模拟的非常理想的差分格式.     

改进型WENO格式的一维爆轰波数值模拟

基于3种改进型加权本质无振荡（WENO）格式,以5阶的WENO-Z格式为参考,采用与经典WENO格式不同形式的非线性权重计算方法,对稳定、中度不稳定和高度不稳定状态下的1维爆轰波系统进行了数值模拟,克服了经典的WENO格式在计算临界点时精度降低的问题.数值实验表明,对于3种状态下的爆轰系统,基于Lagrange多项式重构的WENO-Zη格式计算稳定性好,模拟效果与WENO-Z格式较为一致,比较适合于爆轰波的数值模拟.对于稳定和中度不稳定的爆轰系统,WENO-NS格式和WENO-P格式在数值模拟中会出现小的数值振荡,但计算稳定;而对于高度不稳定爆轰系统,WENO-NS格式计算不稳定,其改进形式WENO-P格式计算稳定性较好.      

基于AUSMDV格式的高超声速钝体绕流的数值模拟

围绕高超声速高温流场的高温效应问题,采用有限体积法,以及熵修正的AUSMDV格式,利用时间分裂法求解了完全气体和化学非平衡流场。求解了SOD2激波管问题,并与精确解进行了对比,分析和验证了AUSMDV格式捕获间断的能力。化学非平衡流求解中采用Park化学动力模型,化学反应源项采用准定常法(The New Quasi-Steady-State-Approach)处理,克服了化学反应源项引起的刚性问题,结果与CEA程序的计算结果进行了比较,研究了高超声速圆柱绕流流场,分析了化学反应对流场的影响。结果证明了计算方法与程序的有效性和实用性。     

叶栅激波流场的数值模拟

给出了ＡＵＳＭ 格式在任意曲线坐标下的形式,将ＡＵＳＭ 格式和三阶精度ＭＵＳＣＬ格式融合,并结合ＬＵ－ＳＧＳ隐式时间推进法,对典型双圆弧叶栅无粘跨音速和超音速流动进行了数值模拟．本文计算结果与文献计算结果相符     

改进型加权本质无振荡格式及其对典型激波问题的高精度数值模拟

为了提高加权本质无振荡(WENO)格式对流场结构的分辨率,提出了性能更优的WENOEX格式。在经典WENO格式的基础上,考虑各子模板光滑因子和整体模板光滑性之间的相互影响,引入表征整体模板光滑度的相关项,构造了新的光滑因子,优化了权分配。理论分析表明,WENO-EX格式保持了原有的设计精度,对间断模板的非线性权配置更大,具有更低的耗散,更高的分辨率。为验证改进格式的计算性能,采用不同的WENO格式来求解线性迁移方程控制的光滑和间断问题,并且求解了欧拉方程控制的双马赫反射和激波旋涡相互作用等几种典型激波问题。计算结果表明:WENO-EX格式在计算光滑问题时,在相同网格下可以取得更小的计算误差;在计算间断问题时,对流场结构有更高的分辨率。     

球形容器泄爆过程影响因素的数值模拟

为了研究球形容器内气体泄爆过程的影响因素及压力传播的基本规律,从能量守恒方程、质量守恒方程、动量守恒方程及化学组分平衡方程出发,建立了泄爆的数学模型,采用k-ε湍流模型和涡耗散概念模型即EDC燃烧模型,利用Fluent软件对该过程进行数值模拟,得到泄爆过程中不同时刻监测点的压力。用响应面曲线的方法讨论泄孔比、泄爆压力、容器体积对球形容器泄爆过程的影响,结果表明:泄孔比越大,危险系数越大;破膜压力过大或过小都将增大危险性;容器体积越大,越容易产生危险。因此,实际球形容器防爆泄压中,合理控制泄爆孔径、破膜压力和容器体积将对泄爆过程中产生的影响及泄爆效果至关重要。     

气相爆轰波贯穿胞格变化研究

根据气相爆轰波胞格结构的规则特则特性建立一种以子单元分解胞格结构的新方法,根据该子单元的性质和斜冲击波关系,求解胞格结构中的三波点对撞问题,并推导贯穿胞格时,爆轰速度的相对波动幅度仅依赖于胞格几何形状的关系,然后,采用爆炸波衰减模型,计算气相爆轰波贯穿胞格时的衰减过程,计算结果与实验符合得较好。     

轴对称二维氢氧爆轰驱动激波管的模拟

对由前向爆轰产生驱动气体并具有局部锥形收缩截面的轴对称激波管利用两相常比热完全气体的欧拉方程组和有限体积TVD格式在局部正交的网格上进行了二维数值模拟.将消除两种介质界面处数值振荡的新方法成功地推广到二维情况.分析了该激波管的力学特性.     

爆炸喷涂爆炸波参数数值计算

对爆炸喷涂燃烧室爆轰波的形成及爆轰产物热力学及动力学参数进行了数值计算．结果表明，爆炸喷涂气相爆轰系统爆轰特性主要取决于气相系统的选择和初始混合气体配比．C2H2 O2气相爆轰系统更适合于喷涂高熔点的陶瓷、金属陶瓷及纳米复合陶瓷粉体，其最优化工艺参数应该选择在零氧平衡附近．而H2 O2气相爆轰系统则适用于喷涂低熔点的金属及其合金，金属粉体的熔点越高，最优化参数越接近于零氧平衡点．     

浓度梯度对瓦斯爆炸影响的数值模拟

针对煤矿巷道中存在瓦斯浓度梯度的问题,基于时间上的TVD Runge-Kutta格式,空间上的5阶加权本质无振荡(WENO)格式离散控制方程组,自主研发了高精度大规模的并行计算程序. 利用该程序模拟了煤矿巷道中的爆轰波的传播过程. 研究结果得出了纵向和横向浓度梯度对瓦斯爆炸的影响规律,模拟结果与理论分析结果基本吻合.      

混合气压力对爆震波特征参数的影响规律

在混合气当量比一定的条件下，采用平衡常数法对烃-空气系统的爆震燃烧产物的平衡温度、平衡压力及爆震波速进行了计算及比较。在此基础上研究了混合气压力对爆震波特征参数的影响规律。结果表明：在同一当量比下，平衡温度和爆震波速随混合气压力变化规律一致，随着初始压力的增大而增大，逐渐趋于一最大值。平衡压力与初始压力呈线性增大。同时，对于不同燃料，分子量越大，其平衡温度、平衡压力就越大。     

气相爆轰波与空气冲击波在斜劈上Mach反射的特征参数比较

以传统的烟迹技术记录胞格结构,证实了气相爆轰波 在管道中与斜劈要互作用发生Mach反射,得到了系统的实验结果,在此基础上详细地分析了爆轰波反射后流场的特征参数（三波咪这线与楔面的夹角X与斜劈楔角X与斜劈楔角及气体初压P0的关系等）,并与空气冲击波在斜劈上Mach反射后流场的特征参数进行了比较,初步探讨了参数偏差的原因,这一研究有助于比较深入,全面地了解爆轰波的前导冲击波和横波在Mach反射中所起的作用。     

无结构三角形网格下运动界面追踪的WENO有限体积法

在处理运动界面追踪问题的流体体积函数(VOF)法的基础上,给出了一种无结构三角形网格下的高分辨率的运动界面捕捉方法.该方法采用高精度的加权本质无振荡(WENO)有限体积格式离散VOF函数的空间导数,采用三阶TVD Runge-kutta方法离散时间导数,采用Lax-Friedrichs通量作为数值流通量.数值试验结果表明,用该方法来进行旋转速度场和剪切速度场的运动界面追踪,可以得到与理论解非常一致的追踪结果.     

气体爆轰波前驱激波变化极限分析

根据气体爆轰波横波强度的实验观测和理论分析结果 ,通过分析横波碰撞前后前驱激波后的压力变化 ,提出了胞格周期中前驱激波速度上限和下限值的估算方法 .并将理论计算的结果与实测上、下限值进行了比较 ,证实了估算方法的可靠性