公共文化服务平台

孙宇涛: 作品数：21 被引量：23H指数：3; 供职机构：北京应用物理与计算数学研究所更多>>; 发文基金：国家自然科学基金中国工程物理研究院科技发展基金国家重点基础研究发展计划更多>>; 相关领域：理学航空宇航科学技术机械工程一般工业技术更多>>

合作作者

基于特征理论的二维无粘Lagrange流体力学有限体积法被引量：4: 2011年; 提出一个求解二维无粘Lagrange流体力学方程的中心型有限体积方法.采用特征理论求解网格节点处的速度及压力,并利用这些物理量更新节点位置及计算网格界面通量.方法适用于结构网格与非结构网格.典型数值实验的结果表明,格式具有较好的收敛性、对称性、能量守恒性及鲁棒性,且能自然地求解多物质流动问题.; 孙宇涛任玉新于明张树道

一种基于MOF界面重构的二维中心型MMALE方法被引量：2: 2016年; 发展了一种基于MOF(Moment of Fluid)界面重构的二维中心型MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法.其中,流体力学方程组采用中心型拉氏方法进行离散求解.混合网格的热力学封闭采用Tipton压力松弛模型.混合网格内的界面重构采用MOF方法,并对MOF方法作了简化和改进.重映步采用一种基于多边形剪裁算法的精确积分守恒重映方法.计算了若干数值例子,包括二维漩涡发展问题、Sedov问题、激波与氦气泡相互作用问题、水中强激波与空气泡相互作用问题、二维RT不稳定性问题等.数值算例表明,该方法具有二阶精度,能够计算界面两侧密度比和压力比很大的问题,并且其健壮性优于交错型MMALE方法,适合计算多介质复杂流体动力学问题.; 贾祖朋孙宇涛

冲击加载下“V”形界面的失稳与湍流混合被引量：3: 2016年; 采用多组分混合物质量分数模型和最小色散可控耗散格式的高分辨率有限体积方法,数值模拟了弱激波冲击不同角度的"V"形空气/SF_6界面的Richtmyer-Meshkov不稳定性问题.激波冲击界面后,在界面附近沉积涡量,形成沿界面规则排列的旋涡结构,同时界面扰动发展形成气泡和尖钉结构.本文统计了界面左端移动速度和界面混合宽度增长率等特征量的演化规律,并与已有的实验结果进行了对比,两者符合较好.讨论了物质界面处的流体向湍流混合发展的过程,随着界面旋涡结构的演化,涡结构之间开始发生相互诱导、并对等现象,并逐渐聚集在几个区域,而多尺度结构也因旋涡的诱导作用在这些区域中产生.通过对由雷诺数定义的惯性尺度进行分析,发现了具有上下边界的惯性尺度区域的形成,对动能能谱的分析发现了-5/3对数率的出现,这同样说明了惯性尺度区域的形成.由于湍流混合转捩与惯性尺度区域的形成是一致的,界面附近流场将发展为湍流.; 李俊涛孙宇涛潘建华任玉新; 关键词：RICHTMYER-MESHKOV不稳定性涡结构湍流混合

求解Euler和Navier-Stokes的多维迎风格式: 本文提出了一种求解Euler方程和Navier-Stokes(NS)方程的多维迎风型二阶有限体积格式.不同于'网格相关'(Grid-aligned)的有限体积方法或者维数分裂的有限差分方法,我们通过在沿控制体表面的法向量...; 任玉新孙宇涛; 关键词：EULER方程 NAVIER-STOKES方程; 文献传递

求解多维欧拉方程的二阶旋转输运格式被引量：4: 2005年; 提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法。不同于“网格相关”(Grid-aligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向)。我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓“红斑”(carbuncle)现象。为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响。这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式。; 孙宇涛任玉新; 关键词：数值模拟欧拉方程激波捕捉

基于特征理论的二维流体力学的单元中心型Lagrange计算方法: 多物质流体力学采用Lagrange方法进行数值计算具有很大的优势。二维Lagrange方法计算流体力学的一个关键点是如何确定离散网格节点的速度。目前广泛应用的交错网格型(staggered-grid) Lagrange方...; 于明孙宇涛

金属约束下定常非理想爆轰的理论研究: 2012年; 为了对金属约束条件下的定常非理想爆轰进行理论研究,对未反应炸药和爆轰产物采用JWL形式状态方程,对金属采用p(ρ,T)形式状态方程。采用过爆轰前导冲击波的流线偏转角在影响域内沿高度线性变化的假设,并且由未反应炸药和金属的冲击波极曲线的交点确定爆轰冲击边缘角,则可从未反应炸药斜冲击波极曲线关系式求出爆轰前导冲击波阵面的形状。采用爆轰前导冲击波阵面之后的流线是直线的假设,则爆轰流动控制方程由偏微分方程变为沿流线的常微分方程,沿着所有流线求解便给出爆轰化学反应区内声速线与化学反应结束线的位置。理论分析同时给出约束金属折射冲击波后面的流体的流动状态。理论结果给出的爆轰化学反应区结构特征和约束金属内的流动状态特征与高精度数值模拟的结果符合良好,说明本文中给出的理论方法具有良好的合理性和适用性。; 于明孙宇涛张文宏; 关键词：流线

爆轰系统建模与模拟不确定度量化的技术框架: 基于可靠性认证中建模与模拟不确定度量化的概念和要求,以及有利于对建模和模拟中存在的问题进行诊断和分离,并能够综合利用不同层级不确定度信息的原则,分析了建模和模拟确定度的来源和组成结构,根据爆轰系统在动作过程中所涉及的不同...; 马智博孙宇涛殷建伟

数值模拟离散误差的后验估计: 验证与确认是评估数值模拟精确性和可信度的主要技术手段。本文着眼于验证环节中的解验证，探讨了如何对数值模拟中离散误差造成的性能参数不确定度进行后验估计。基于Richardson 外推公式，本文研究了几种收敛解的估算方法，尝...; 殷建伟马智博孙宇涛; 关键词：数值模拟离散误差性能参数后验估计; 文献传递