磁流体力学数值计算处理方法

介绍电磁内爆实验中求解磁流体力学方程组的处理方法。方程组采用分裂法以克服非线性耦合可能引起的数值不稳定性；采用隐式差分格式避免过小时间步长限制。隐式格式采用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法叠代求解，使用比较简单的叠代初值的给定方法。对一些特殊边界处理作了说明。     

一种追踪多介质流体界面运动的NND数值模拟方法

把界面捕捉等效方程、Level-Set方程和欧拉方程组耦合,在Stiffened状态方程下,采用高分辨率NND格式求解流体力学方程组并用Level-Set函数捕捉界面的位置。对二维情况下激波和气泡相互作用的问题进行数值模拟,并与波传算法的模拟结果进行比较。计算结果表明该方法能有效的抑制间断附近的非物理振荡,有很强的捕捉界面的能力。     

拟线性粘弹性方程新H~1-Galerkin最低阶混合元格式的高精度分析

利用双线性元和零阶Raviart-Thomas元,针对拟线性粘弹性方程建立新的H~1-Galerkin混合元逼近格式.在半离散格式下,给出原始变量u的H~1模和应力=?ut的H(div;?)模的超逼近性和超收敛结果.同时,导出向后欧拉格式和Crank-Nicolson-Galerkin格式的最优误差估计.最后,通过数值算例表明逼近格式是有效的.     

《数学之友》投稿及查询启事

1.《数学之友》2017年选题范围涉及：数学教育基本理论研究、中国数学教育传统及其特征研究、数学核心素养问题研究、数学比较教育研究、中小学数学课程与教学改革的理论与实践研究、现代教育技术与数学课程及教学整合研究、解题的理论与实践研究、中学生数学学习的心理研究等。2.来稿论文务请格式规范、言之有物、有一定的独创性,作者信息要完备（姓名、单位、邮编）,来稿一般应有参考文献。     

非定常俯抑振荡下的横向喷流数值模拟

采用高精度格式数值求解RANS方程,研究了定常状态下横向喷流流场,压力分布计算结果与实验结果基本吻合,并捕捉到喷流干扰流场中多种流动结构.在非定常计算过程中,飞行器的振动引起了法向力和俯仰力矩系数的相位滞后,推力放大因子随俯仰角周期变化.飞行器振动过程中,喷流流场的动态气动特性与稳态喷流有明显的区别,因此在利用横向喷流对飞行器进行姿态控制时,应该考虑由于飞行器姿态的变化对横向喷流所产生的非定常影响问题.     

气体运动论数值算法在微槽道流中的应用研究

简要介绍基于Boltzmann模型方程的气体运动论数值算法基本思想及其对二维微槽道流动问题数值计算的推广，并阐述适用于微尺度流动问题的气体运动论边界条件数值处理方法。通过对压力驱动的二维微槽道流动问题进行数值模拟，将不同Knudsen数下的微槽道流计算结果分别与有关DSMC模拟值和经滑移流理论修正的N—S方程解进行比较分析，表明基于Boltzmann模型方程的气体运动论数值算法对微槽道气体流动问题具有很好的模拟能力。     

Level Set追踪等温非牛顿熔全前沿界面

应用Level Set方法追踪薄壁型腔内Hele-Shaw熔体流动前沿界面,采用5阶加权本质无振荡格式耦合中心差分格式实现了充填阶段的动态模拟.准确追踪到了不同时刻熔体前沿界面,并得到了对应的压力等值线分布,数值结果表明Level Set方法是准确追踪注塑成型熔体前沿界面的一种行之有效的方法.     

采空区流场非达西渗流一种新的迭代算法

对Bachmat非线性渗流方程求解问题，结合采空区问题的力学特点，提出基于变渗透率达西 (Darcy)渗流求解非线性渗流的迭代算法. 建立了冒落采空区非线性漏风流态的有限元数 值模型，通过对复杂边界采空区的漏风渗流的计算，得到与实际流态更接近的流动规律(风 压等值线和流函数线)和采空区漏风强度分布(速度场). 结果表明，迭代算法是振荡性收 敛的，收敛速度快；与Darcy渗流相比，采空区非线性渗流速度场趋向于平缓，在工作面附 近30\,m范围内最大风速降低到原来的0.61倍，其余大部分区域速度略有增加. 迭代方法满 足工程要求.     

球面聚心气相爆轰波传播过程的数值研究

应用频散可控耗散格式对球面聚心气相爆轰波的传播过程进行了数值模拟 研究. 通过跟踪波阵面上压力和温度变化，分析这些参数在爆轰波传播过程中的演变规律， 及其与几何尺度和初始条件之间的依赖关系. 研究结果表明，仅在远场波面压力的变化近似 只依赖于r/R，对称中心附近则需要考虑初始半径$R$的影响；波面压力与初始压力的变化呈 线性关系；汇聚过程中温度升高比压力慢得多.     

高阶多维半离散中心迎风格式及其应用

提出了求解多维对流-扩散方程的四阶半离散中心迎风格式。该格式以中心加权基本无振荡(CWENO)重构为基础,同时考虑到在Riemann扇内波传播的局部速度,从而更加准确地估计出了局部Riemann扇的宽度,最终既回避了网格的交错,又降低了格式的数值粘性,建立了介于迎风格式和中心格式之间的半离散中心迎风格式。本文还将该四阶半离散中心迎风格式与涡度-流函数方法相结合,有效地求解了二维不可压Euler方程组和Navier-Stokes方程组。