SPH后处理研究

提出一套相对完整的光滑粒子流体动力学(SPH)后处理方法.凸区域和严重变形的非凸区域上的SPH计算结果都可用该方法处理.首先,对SPH粒子集进行Delaunay三角化,可得到由粒子作为节点的三角单元集;根据每个单元中三个节点是否彼此为粒子作用对,决定是否将该单元从单元集中删除.将保留下的单元作为有限单元并利用它们节点上的函数值,根据有限元插值方法即可得单元内部任何一点的函数值.根据该方法,可提取介质的自由表面.数值算例表明方法可行.对含固体壁面、严重粒子飞溅、多重介质互相作用的情况,提出相应对策;为SPH以至其它无网格方法形成较通用的后处理软件提供可行的途径.     

部分状粒子处理方法及其对云微物理参数测量的影响

作为云微物理过程测量的重要利器,机载云降水粒子成像仪在云降水物理与人工影响天气研究中具有重要的作用.从采样结果来看,机载云降水粒子成像仪所测粒子图像中含有大量的粒子图像仅是粒子的一部分而已,即部分状粒子.因其数量较多,对该类粒子所选处理方法不同,会引起测量结果的很大差异.本文介绍并分析了现有部分状粒子处理方法的优劣,通过对部分状粒子的再定义与粒子形状分类,提出了一个融合粒子形状识别技术、"粒径重构"和"中心在内"方法的新的部分状粒子处理方法;利用实测数据,对所提方法与现有方法进行了云微物理参量处理结果的对比,发现本文所提方法与"粒径重构"方法处理结果比较一致,能较好地克服"整体在内"与"中心在内"两种方法存在的缺陷;同时,在针柱状粒子占比较多情形下,本文所提方法要比"粒径重构"方法处理后的结果相对合理.因此本文所提方法对仪器所测粒子数据处理具有更好的适应性.     

任意不可压多介质流的粒子有限元方法模拟

基于粒子有限元方法(particle finite element method,PFEM),利用细分混合单元的界面识别思想,模拟种类任意多的不可压多介质流问题.对分步算法采用基于有限增量微积分理论的稳定措施,以适应流体特性差异;将混合单元细分为代表单一流体的小单元,进而得到流体间的边界;通过加密边界、控制粒子速度、自动检查穿透来防止粒子穿透外部边界.瑞利-泰勒不稳定性和水柱在空气中倒塌的模拟与已有结果的对比验证了PFEM及界面识别方法的可靠性和准确性.七种流体混合的模拟结果表明PFEM可有效处理任意多种类不相溶流体的混合流动问题.     

轴对称变形核中配对粒子态的构成

量子多粒子体系的配对粒子态的构成与体系对称性密切相关.轴对称变形核中的配对粒子态可有三种构成方式,即单极对,±Ω和Signature对,分别是(J2,J3),J3,和(J32,R1(π))的共同本征态.在只计及对力时,三种形式的对力得出的能谱完全相同,常被等价地混用.但在推转壳模型中,由于三种配对态在空间转动下的性质不尽相同,能谱及有关性质随ω的变化将有不同程度的差异.     

弹性力学的重构核粒子边界无单元法

将重构核粒子法(RKPM)和边界积分方程方法结合，提出了一种新的边界积分方程无网格方法——重构核粒子边界无单元法(RKP-BEFM).对弹性力学问题，推导了其重构核粒子边界无单元法的公式，研究其数值积分方案，建立了重构核粒子边界无单元法离散化边界积分方程，并推导了重构核粒子边界无单元法的内点位移和应力积分公式.重构核粒子法形成的形函数具有重构核函数的光滑性，且能再现多项式在插值点的精确值，所以本方法具有更高的精度.最后给出了数值算例，验证了本方法的有效性和正确性. 关键词： 重构核粒子法 弹性力学 边界无单元法     

全同多准粒子转动带的内部结构

用粒子数守恒方法研究了全同带的内部结构.几乎全同转动惯量的出现是下列三种效应竞争的结果:即壳效应(包括形变的改变)、配对(反顺排)效应和堵塞(反配对)效应.粒子数守恒方法计算出的全同多准粒子带的转动惯量与实验值符合很好,计算中无自由参数.     

中能重离子反应中的多单元轻粒子阵列探测器

本文描述了一个能在中能(10—100MeV/u)重离子核反应测量中鉴别轻粒子、测量粒子多重性,并能给出轻粒子能量和位置信息的大立体角(θ:5°—20°,φ0°—360°)、多单元(36单元)阵列探测器.初步调试结果表明该探测器对轻粒子(p、α等)有较好的分辨.     

SPH方法进出口数值边界条件研究

基于虚粒子概念,提出进出口数值边界条件处理方法.在进口边界外设置进口区域,并赋予进口虚粒子相应的物理量;出口边界外设置出口和缓冲区域,出口虚粒子的物理量由计算获得,缓冲区的虚粒子初始物理量是指定的;根据每个时间步的流动情况,改变粒子的区域属性,添加/删除相应的粒子.利用拉格朗日形式的SPH方法,通过管内流动问题验证数值进出口边界方法的适用性,研究进出口边界条件在激波管、绕流问题中的应用.提出的进出口边界处理方法,避免了边界附近流体粒子积分截断问题,保证流体粒子能够流出边界,激波能够透射边界.     

基于中国散裂中子源的商用静态随机存取存储器中子单粒子效应实验研究

利用中国散裂中子源反角白光中子束线开展13款商用静态随机存取存储器的中子单粒子效应实验.研究了测试图形、特征尺寸和版图工艺差异对单粒子效应的影响.结果表明测试图形对器件的单粒子翻转截面影响不大,但对部分器件的多单元翻转占比有较大的影响;特征尺寸对器件单粒子翻转截面的影响没有明显的规律,但对多单元翻转的影响规律明显,多单元翻转占比和最大位数都随着特征尺寸的降低而增大;器件版图工艺差异对器件的单粒子翻转截面和多单元翻转占比都有较大的影响.此外,通过与高原辐照实验结果对比,发现在反角白光中子源获得的多单元翻转占比小于高原辐照实验的结果,其原因是反角白光中子源实验中,中子的最高能量和高能成分占比偏小,且中子束流只有垂直入射.因此,利用反角白光中子源评估器件的大气中子单粒子效应时可能会低估多单元翻转情况.本文的结果可为研究者利用反角白光中子源开展相关研究提供参考.     

一种有限元-光滑粒子流体动力学耦合算法

介绍一种新的有限元-光滑粒子流体动力学耦合算法.通过将SPH法中的粒子固连到单元边界上实现有限元计算与SPH计算的耦合,粒子可以连接到单元边界上的任意点,且一条单元边界可以与多个粒子连接,耦合界面处粒子分布与单元分布无相互限制.应用该算法模拟应力波传播和高速冲击问题.结果表明,算法在应力波传播计算中具有良好的精度,在高速冲击问题中计算精度与SPH法相近的,而计算效率大大提高.     

基于RANSAC方法的星表自动匹配方法

本文把星表匹配问题转化成计算机视觉中的对应点匹配问题,接着使用二维Delaunay三角化方法对星表进行三角化,然后用RANSAC方法对星表的自动匹配问题进行了研究。研究结果表明,我们的方法能够快速有效地完成星表匹配。在总共960个样本中,除了两幅不符合匹配的星表之外,全部获得了正确的匹配结果。     

水珠滴落的最小二乘粒子有限元方法模拟

提出了一种最小二乘粒子有限元方法,用其模拟了二维水珠滴落水面并飞溅散开的过程.该法基于拉格朗日描述,在每个时间步上使用扩展的Delaunay划分更新计算网格,并应用α形方法识别自由面形状;用最小二乘有限元方法离散流体运动的Navier-Stokes方程,并推导了一种自适应时间步长方案以提高计算效率和鲁棒性;引入网格拉伸技术修正减小流体质量误差.对水滴飞溅进行仿真,得到了与商用软件Flow-3d比较符合的结果,且具有更清晰锐利的自由面. 关键词： 滴落 网格划分 α形')" href="#">α形 最小二乘有限元     

双参数曲面的泡泡网格化方法

为将双参数曲面离散成高质量的网格,首先在参数域内利用各向异性的非均匀泡泡布点方法优化布点,然后用各向异性Delaunay三角化方法将参数域网格化,最后用映射法得到双参数曲面的离散网格.参数域中的节点由二阶黎曼度量矩阵控制,该度量矩阵由三维曲面的网格度量矩阵和曲面参数方程的梯度计算得到.数值算例表明,泡泡布点法在参数域上能生成满足度量矩阵要求的节点集,将节点连接成网格并投影回曲面,所得曲面网格具有很高的质量.     

结合前沿推进的Delaunay三角化网格生成及应用

采用一种新的混合网格生成方法,生成复杂区域的非结构化网格.结合前沿推进法和Delaunay三角化两种非结构网格生成方法的特点,在边界处采用前沿推进法进行三角形初始网格的生成,在边界区域内部采用Delaunay三角化方法自动生成内部节点.分析表明,该算法简化网格生成过程,能够快速有效地生成非结构化网格.在计算时间以及网格的均匀性方面与其他方法相比具有一定的优势.最后,用混合网格生成方法生成方柱绕流的计算域网格,并运用基于特征线方程的分离算法进行流场计算.     

内部带特征约束的任意平面域的三角形网格生成方法

基于Ｄｅｌａｕｎａｙ三角化方法 ,给出了内部带有任意复杂特征约束的任意平面域的三角形网格生成方法。区域内部特征约束可以是点、直线、圆、圆弧、平面样条曲线以及由它们组合构成的封闭环或者不封闭环。     

用非结构网格与欧拉方程计算复杂区域的二维流动

提出用Delaunay三角化方法生成非结构网格的一种过程。所生成的网格可用于复杂多连通域内的可压流计算。采用Euler方程和格心有限体积法,研制出程序,给出了算例。     

一种新的非结构四边形网格生成方法及软件

提出一种间接生成非结构四边形网格的方法.采用Delaunay三角形网格生成的算法生成三角形单元,在三角形单元中定义生成四边形网格的前沿,根据三角形的特点(边、角大小)以及生成四边形前沿的方向,将两个三角形合并为一个四边形单元.此方法的特点是不增加节点,仅改变点与点之间的拓扑连接,使单元基本变为四边形占优的单元.然后研究四边形占优的网格转换为100%的四边形网格的转换模版.并应用于实际物理问题.数值实验显示了该方法生成非结构四边形网格的特点.     

三维非结构网格自动生成

采用各向异性平面非结构网格生成技术对曲面的参数平面进行三角剖分,从而得到曲面的非结构网格,作为三维非结构网格的边界网格.应用推进面法生成网格内点,增量法将生成的内点逐点插入现有网格进行网格细化,得到三维计算域的Delaunay非结构网格.讨论了非结构网格质量优化方法.给出几个算例说明方法的应用.     

基于约束Delaunay三角化的二维非结构网格生成方法

给出基于局部重构和边交换技术的两种约束Delaunay三角剖分方法并证明其收敛性.采用边界指示法恢复流场形状;在预设尺度的指导下融合流场边界曲率、中轴线、梯度限制等信息修正流场尺度;运用Spring方法布置边界点,通过符号面积函数和概率筛选法布置计算区域节点;运用Spring-Laplace方法优化节点位置,伴同边交换和边吞噬技术优化网格结构.该方法可自由进行局部自适应加密或稀疏,并应用于映射曲面网格生成和移动网格技术.