时一空守恒元解元（CE／SE）方法综述

时一空守恒元解元方法是近年来兴起的一种全新的高分辨率守恒型方程计算方法．它具有物理概念清晰、计算精度高和格式构造简单等优点,是一种具有广阔发展前景的计算方法．本文详细地介绍了cE／ss方法的基本原理、发展历史、应用情况和最新进展、并指出了当前研究的不足和发展方向．     

Gradient-augmented hybrid interface capturing method for incompressible two-phase flow

Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives(normals and curvatures) and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.     

用于测量材料高温动态力学性能的SHPB技术

讨论了高温下分离式Hopkinson压杆(SHPB)实验的两种方案:单独加热试样并快速对杆与同时加热试样和杆再修正温度梯度的影响。前者避免了杆被加热但是操作复杂,后者简单稳定。为了分析后者温度梯度的影响设计了一个简化模型,采用数值计算进行修正,提出了精度适当的假设,并且进行了实验验证。测量了试件和杆一起加热时杆上的温度分布,证实计算采用的温度分布是合理的。据此得到的专用于材料高温动态力学性能测试的SHPB数据处理公式简单实用,便于推广使用。     

An improved two-dimensional unstructured CE/SE scheme for capturing shock waves

In this paper, the accuracy of Chang’s unstructured space–time conservation element and solution element (CE/SE) scheme is analysed for the first time. Based on a redefinition of conservation elements and solution elements, an improved two-dimensional (2D) unstructured CE/SE scheme with an adjustable parameter β is proposed to accurately capture shock waves. The new scheme can be applied to any type of grid without special treatment. Compared with Chang’s original parameter α, larger β dose not cost extra computational resources. Numerical tests reveal that the new scheme is not only clear in physical concept, compact and highly accurate but also more capable of capturing shock waves than the popular fifth-order accurate weighted essentially non-oscillatory scheme.     

Simulation of shock-induced instability using an essentially conservative adaptive CE/SE method

An essentially conservative adaptive space time conservation element and solution element(CE/SE)method is proposed for the effective simulation of shock-induced instability with low computational cost.Its implementation is based on redefined conservation elements(CEs)and solution elements(SEs),optimized interpolations and a Courant number insensitive CE/SE scheme.This approach is used in two applications,the Woodward double Mach reflection and a twocomponent Richtmyer–Meshkov instability experiment.This scheme reveals the essential features of the investigated cases,captures small unstable structures,and yields a solution that is consistent with the results from experiments or other high order methods.     

爆轰驱动金属飞层对碰凸起和微射流形成的数值模拟研究

本文对柱面两极点起爆情况下滑移爆轰波驱动两层金属飞层对碰凸起和微射流形成进行了模拟研究. 铅飞层内界面走时计算结果与实验结果能够较好符合. 在两极位置铅飞层内部出现断裂并形成空腔, 内壁面则形成鼓包型凸起; 在赤道位置飞层内壁面凸起后断裂产生大尺度金属颗粒, 其和微喷射形成的小尺度颗粒叠加构成了对碰区凸起现象. 在铅飞层内表面微喷射现象的研究中发现, 两极附近的微喷物质最大速度逐渐下降, 而对碰区附近的微喷颗粒最大速度反而随时间逐渐增高. 之后, 通过设计沟槽型微喷计算模型, 验证了在两极和赤道上铅飞层内表面产生的初始微喷射最大速度能够由同一均匀缺陷表面所产生. 最后, 通过数值模拟分析研究初步给出了该问题中抑制金属飞层对碰凸起和微喷现象的方法.     

计算域可变的CEL方法

为了更好地兼顾拉氏方法和欧拉方法各自的特长,提出一种可将拉氏介质映射到欧拉计算域的耦合欧拉-拉格朗日(CEL)方法。通过这种映射,将欧拉-拉格朗日重叠区域的接触面协调问题转换为欧拉区域内的多介质计算问题,简化了CEL方法的构造过程。通过与侵彻实验和结构对爆炸冲击波响应实验的比较,验证了新算法,计算结果与实验数据符合较好。