期刊
  出版年
  关键词
结果中检索 Open Search
Please wait a minute...
选择: 显示/隐藏图片
定常不可压Navier-Stokes方程的两水平grad-div稳定化有限元方法
王雅莉, 郑波, 尚月强
计算物理    2024, 41 (4): 418-425.   DOI: 10.19596/j.cnki.1001-246x.8731
摘要122)   HTML4)    PDF (5605KB)(444)   

使用标准的混合有限元方法数值求解定常不可压Navier-Stokes方程所得速度解的精度常常受压力的影响。为了克服或减弱压力对速度精度的影响,本文将grad-div稳定化方法和两水平有限元方法相结合,提出数值求解定常不可压Navier-Stokes方程的两水平grad-div稳定化有限元方法。首先在粗网格上求解grad-div稳定化的非线性Navier-Stokes问题,然后在细网格上分别求解grad-div稳定化的Stokes型、Newton型和Oseen型的线性问题。最后给出数值算例验证两水平grad-div稳定化有限元方法的高效性。

图表 | 参考文献 | 相关文章 | 多维度评价
一氟二氯乙烷分子在外电场中的光谱特征和解离特性
刘启帆, 闫好奎, 布玛丽亚·阿布力米提, 向梅, 安桓, 郑敬严
计算物理    2023, 40 (4): 453-460.   DOI: 10.19596/j.cnki.1001-246x.8591
摘要150)   HTML5)    PDF (5761KB)(644)   

利用密度泛函理论在B3LYP/6-311+G(d, p)水平上研究不同外加电场下(-0.05~0.05 a.u.)一氟二氯乙烷分子的光谱特征和解离特性, 其中包括该分子基态结构、总能量、偶极矩、最高占据轨道能级、最低空轨道能级、能隙、红外光谱、拉曼光谱、紫外可见光吸收光谱及C原子与Cl原子间势能曲线等数据在电场下的变化趋势。在y轴方向上, 随着负向电场的增强, C原子与Cl原子核间距增长、最高占据轨道能级减小, 体系总能量、最低未占据轨道能级、能隙先增大后减小, 偶极矩先减小后增加。外加电场会影响一氟二氯乙烷分子红外光谱、拉曼光谱与紫外可见光吸收光谱的吸收强度和吸收峰频率, 红外光谱、拉曼光谱与紫外可见光吸收光谱随着电场变化出现红移或蓝移现象。C原子与Cl原子间势垒随着负向电场增加逐渐减小, 并在电场达到-0.05 a.u.时C原子与其中一个Cl原子发生断裂, 当分子中一个C—Cl键断裂后, 施加强度为-0.04 a.u.的电场时另外一个C—Cl键发生断裂, 分子在电场下发生逐步解离, 研究结果完善了一氟二氯乙烷分子受外电场影响的理论数据。

图表 | 参考文献 | 相关文章 | 多维度评价
外电场下CHBr3分子的光谱和解离特性
乃皮赛·吾买尔江, 闫好奎, 布玛丽亚·阿布力米提, 王丹琦, 向梅, 安桓
计算物理    2022, 39 (5): 624-630.   DOI: 10.19596/j.cnki.1001-246x.8502
摘要264)   HTML9)    PDF (1379KB)(1226)   

利用密度泛函理论在B3LYP/6-311G++(d, p)基组上优化和计算分析三溴甲烷分子处于外电场中的基态稳定构型, 包括分子结构、总能量、偶极矩、键长、电荷分布、轨道能级、红外光谱和拉曼光谱及解离势能曲线等。结果表明: 随着Y轴(平行于1C-4Br连线)方向外电场的增加(0a.u.~0.02a.u.), 分子总能量减小, 1C-4Br键键长增大, 分子偶极矩增大, 分子能隙逐渐减小, 分子红外光谱呈现蓝移。对解离势能曲线的计算分析发现, 强度为0.02a.u.的外电场使得1C-4Br键断裂而降解, 为三溴甲烷进行外电场降解提供参考依据。

图表 | 参考文献 | 相关文章 | 多维度评价
不可压缩流的并行两水平稳定有限元算法
朱家莉, 尚月强
计算物理    2022, 39 (3): 309-317.   DOI: 10.19596/j.cnki.1001-246x.8411
摘要273)   HTML19)    PDF (5666KB)(1402)   

提出一种数值求解定常不可压缩Stokes方程的并行两水平Grad-div稳定有限元算法。首先在粗网格中求解Grad-div稳定化的全局解, 再在相互重叠的细网格子区域上并行纠正。通过对稳定化参数、粗细网格尺寸恰当的选取, 该方法可得到最优收敛率, 数值结果验证了算法的高效性。

图表 | 参考文献 | 相关文章 | 多维度评价
一种基于新型斜率限制器的理想磁流体方程的高分辨率熵相容格式
沈亚玲, 封建湖, 郑素佩, 李雅蓉
计算物理    2022, 39 (3): 297-308.   DOI: 10.19596/j.cnki.1001-246x.8405
摘要275)   HTML14)    PDF (10894KB)(917)   

针对磁流体动力学方程, 通过分析数据重建所需的条件, 构造一种基于MUSCL(Monotone Upstream-Centred Scheme for Conservation Laws)型重建方法的斜率限制器, 获得了一种求解理想磁流体动力学方程的高分辨率熵相容格式。该格式在解的光滑区域具有高精度; 在解的间断区域可以合理地控制耗散, 可有效避免非物理现象的产生。采用熵稳定格式、熵相容格式和新的高分辨率熵相容格式对一维、二维理想磁流体动力学方程进行数值模拟。结果表明: 新格式能准确地捕捉解的结构, 且具有无振荡、高分辨、鲁棒等特性。

图表 | 参考文献 | 相关文章 | 多维度评价
基于守恒高阶模型和支持向量机的多车道车辆换道模型
张立灿, 郭明旻, 林志阳, 张鹏, 段雅丽
计算物理    2022, 39 (1): 83-95.   DOI: 10.19596/j.cnki.1001-246x.8339
摘要179)   HTML3)    PDF (4453KB)(647)   

针对高速公路车辆换道问题, 提出一个多车道车辆换道模型。利用支持向量机(SVM)在多维特征下二分类问题的优势, 将SVM和Lagrange坐标下的高阶守恒模型(CHO)结合, 通过全离散跟车模型生成原始数据, 采用SMOTE(Synthetic Minority Oversampling Technique)算法对数据进行预处理, 采用双指标评估度SVM进行训练, 建立多车道车辆换道仿真模型。仿真结果表明: 基于支持向量机和CHO模型的换道模型, 驾驶车能够就当前的驾驶环境, 准确地作出决策, 有效地模拟高速公路上真实的多车道驾驶情况。

图表 | 参考文献 | 相关文章 | 多维度评价
金、银纳米结构中的非局域吸收特性
刘建晓, 刘啸岚, 陆华丽, 高英杰, 姜严, 唐万春
计算物理    2021, 38 (2): 206-214.   DOI: 10.19596/j.cnki.1001-246x.8211
摘要230)   HTML2)    PDF (7628KB)(948)   
采用时域有限差分方法(FDTD)研究金、银纳米结构的非局域吸收。对金属材料的Drude-Lorentz模型进行离散差分,详细推导各项迭代系数。研究一维、二维、三维金、银纳米结构的吸收特性。讨论一维分层纳米薄膜吸收率与材料的关系。结果表明:非局域分层介质板仍然符合等效介质理论。在一维、二维情况下,金、银纳米结构的非局域吸收峰只与纳米尺度相关而与材料无关。三维纳米结构的非局域吸收特性不仅与纳米尺度相关还与材料特性相关。三维纳米结构尺度更小,非局域效应表现也更为明显,吸收峰的蓝移也更大。因此,三维纳米颗粒有望在纳米材料以及纳米器件的研究中发挥更大的作用。
参考文献 | 相关文章 | 多维度评价
基于守恒高阶模型的动态交通分配建模与仿真
李浩宇, 林志阳, 张鹏, 段雅丽
计算物理    2020, 37 (6): 687-699.   DOI: 10.19596/j.cnki.1001-246x.8143
摘要483)   HTML1)    PDF (3978KB)(1144)   
针对路网中的动态交通分配问题,采用高阶守恒模型(CHO)进行建模与数值研究,并推广高阶守恒模型二进二出路口的Riemann问题;同时将高阶守恒模型与动态网络加载(DNL)模型相结合,通过变分不等式对动态网络加载模型进行分析.数值模拟采用一阶有限体积法求解高阶守恒模型,同时采用梯度下降方法迭代求解动态网络加载模型的变分不等式问题,最终以动态用户最优条件为目标实现分配均衡.数值结果表明CHO模型与DNL模型结合解决动态交通分配问题是可行的,对传统模型的研究有一定的指导意义.
参考文献 | 相关文章 | 多维度评价
掺杂菱形BN片的石墨烯纳米带的电子特性
马瑞, 张华林
计算物理    2019, 36 (1): 99-105.   DOI: 10.19596/j.cnki.1001-246x.7803
摘要399)   HTML2)    PDF (10443KB)(1255)   
采用基于密度泛函理论的第一性原理方法,系统研究掺杂菱形BN片的石墨烯纳米带的电子特性.掺杂使扶手椅型石墨烯纳米带(AGNRs)的带隙增大,不同位置掺杂AGNRs的带隙大小略有差异.在无磁性态,无论是否掺杂,锯齿型石墨烯纳米带(ZGNRs)都为金属.在铁磁态,掺杂使ZGNRs由金属转变为半导体.而处于反铁磁态时,无论是否掺杂,ZGNRs都为半导体,掺杂使其带隙发生改变.掺杂的AGNRs和ZGNRs的结构稳定,掺杂ZGNRs的基态为反铁磁态.掺杂菱形BN片可以有效调控GNRs的电子特性.
参考文献 | 相关文章 | 多维度评价
二维三角格点系统中的拓扑陈数
郁华玲, 高雨, 翟章印
计算物理    2018, 35 (5): 606-612.   DOI: 10.19596/j.cnki.1001-246x.7696
摘要623)   HTML1)    PDF (7847KB)(1360)   
利用紧束缚模型对二维三角周期格点中各能带的陈数分布进行研究.通过严格对角化方法得到体系能量本征值和对应的本征态,再利用Kubo公式计算出量子化的霍尔电导、态密度及各扩展态对应的陈数.在傅里叶变换下将哈密顿量转换到k空间从而得到体系的能谱分布.研究表明:次近邻格点之间的跳跃积分t'的不同取值影响体系各能带对应的陈数分布,计算得到当t'=1/2时体系三个能带从低到高对应的陈数分布为{-4,5,-1},t'=-1/2时其对应陈数分布变化为{2,-4,2},而t'=±1/4时对应的陈数分布都为{2,-1,-1}.同时发现:能谱帯隙的宽度和对应霍尔平台的宽度一致,并且k空间的能带越平坦,其对应的在霍尔电导跳跃处的态密度峰就越高越尖锐,而该处霍尔电导跳跃就越陡峭.
参考文献 | 相关文章 | 多维度评价
基于数值模拟的飞行员高速气流吹袭损伤分析
张朋, 田佳林, 韩鲁佳, 刘洁
计算物理    2016, 33 (2): 170-176.  
摘要308)   HTML0)    PDF (948KB)(1191)   
研究飞行员高速气流吹袭的损伤机理,应用尺度适应模型对人椅系统外流场进行数值模拟.在马赫数为0.6,迎角为-90°~90°条件下的数值计算结果与风洞试验结果吻合较好.对不同马赫数、不同迎角的无防护状态下的人椅系统进行数值模拟,研究发现高速气流吹袭下飞行员面颈部、胸腹部、上臂内侧、小腿迎风侧压力较大;增大人椅系统的迎角可以改善飞行员身体各部分的受力情况.
参考文献 | 相关文章 | 多维度评价
Burgers-Korteweg-de Vries复合方程的格子Boltzmann方法模拟
段雅丽, 陈先进, 孔令华
计算物理    2015, 32 (6): 639-648.  
摘要482)      PDF (1536KB)(1241)   
针对Burgers-Korteweg-de Vries(cBKdV)复合方程提出一种格子Boltzmann模型.通过恰当地处理色散项uxxx并运用Chapman-Enskog展开从格子Boltzmann方程推导出宏观方程,从而得到联系微观量与宏观量的局部平衡分布函数.对不同微分方程进行数值实验,数值解与解析解非常吻合,相比于其它数值结果,该格子Boltzmann模型的数值结果更精确,说明该数值模型的高效性.
相关文章 | 多维度评价
激光相变烧蚀二维结构网格程序的数值计算
陶应学, 赵强, 陈发良
计算物理    2014, 31 (2): 165-172.  
摘要293)      PDF (1321KB)(1208)   
研究激光相变烧蚀二维结构网格程序,对液-气相变、固液汽三相物质进行整体模拟.激光相变烧蚀二维结构网格程序采用气化相变开始后在迎激光方向局部加密的重分方法,动量守恒方程采用Wilkins三角形回路积分法,能量守恒方程采用基于节点重构的扩散方程有限体积格式.
相关文章 | 多维度评价
Klein-Gordon-Schrödinger方程的辛Fourier拟谱格式
王兰, 马院萍, 孔令华, 段雅丽
计算物理    2011, 28 (2): 275-282.  
摘要292)      PDF (293KB)(977)   
主要讨论Klein-Gordon-Schrödinger方程的Fourier拟谱辛格式,包括中点公式和Störmer/Verlet格式.首先构造一个哈密尔顿方程,针对此哈密尔顿方程,在空间方向用Fourier拟谱离散得到一个有限维的哈密尔顿系统,对此有限维系统在时间方向用Störmer/Verlet方法离散得到KGS方程的完全显式的辛格式.中点格式虽然是隐式的但效率也很高,且具有质量守恒律.数值实验表明,辛格式能够在长时间内很好地模拟各类孤立波.
相关文章 | 多维度评价
煤油超燃冲压发动机三维大规模并行数值模拟
杨顺华, 乐嘉陵, 赵慧勇, 郑忠华
计算物理    2009, 26 (4): 534-540.  
摘要314)      PDF (418KB)(1416)   
在我国巨型计算机上,采用1024个CPU对煤油燃料超燃冲压发动机燃烧流场进行三维大规模并行数值模拟.计算软件采用自主研发的并行软件AHL3D,控制方程采用雷诺平均的N-S方程,无粘项计算采用Steger-Warning矢通量分裂格式,湍流模型采用k-ω双方程模型,煤油分子式采用正癸烷代替.计算给出的发动机壁面压强分布与试验测量结果有较好的一致性.结果表明,凹槽是发动机主要的着火区和火焰稳定区,同时,由于燃料的喷注形成的回流区也起到一定的稳焰作用.计算结果验证了AHL3D程序和采用的理论模型可以用于模拟煤油燃料超燃冲压发动机内部的复杂流动.
相关文章 | 多维度评价
径向基函数及其耦合方法在电磁场数值计算中的应用
张淮清, 俞集辉, 郑亚利
计算物理    2009, 26 (2): 299-310.  
摘要261)      PDF (660KB)(1400)   
针对带电导体附近急剧变化的位函数和场函数这一难于处理的边界条件,将小波函数的紧支撑特性和全域径向基函数(RBF)的高精度逼近能力相结合,提出电磁场边值问题求解的耦合方法并应用于接地金属槽/箱的数值计算中;将径向基函数无网格方法引入波导本征值的计算中,给出其求解本征问题的思路,建立相应的离散方程,分析矩形、圆形和脊形波导的本征值并与有限元方法进行比较.数值仿真实验表明,径向基函数及其耦合方法在分析电磁场边值和本征值问题时是有效的且具有实现简单、节点少和精度高的优势.
相关文章 | 多维度评价
纳米尺度绕流现象中流体分子运动行为研究
李印实, 孙杰, 何雅玲
计算物理    2008, 25 (6): 711-717.  
摘要267)      PDF (480KB)(1026)   
采用分子动力学模拟方法,研究纳米尺度低Re数(Re=20)条件下,氩流体流过铂(FCC(100))金属圆柱的绕流现象.从流线和速度矢量分布两方面刻画涡的周期性产生、发展和脱落过程.结果显示,一个涡脱落周期大约为0.3 ns,斯特劳哈尔数St约为0.2;时均对称涡的长度是圆柱直径的1.4倍,涡径与圆柱直径相当.从速度场和密度场的角度刻画回流区,并进一步研究典型区域(圆柱中心线附近区域及下游邻近区域)的速度分布特性.
相关文章 | 多维度评价
吸气式高超声速飞行器气动力气动热的数值模拟方法及应用
贺旭照, 赵慧勇, 乐嘉陵
计算物理    2008, 25 (5): 555-560.  
摘要344)      PDF (325KB)(1305)   
对吸气式高超声速飞行器而言,物面热流和摩阻的准确预测对飞行器设计及安全十分关键.介绍采用CFD准确预测气动力和气动热的方法,包括流动的控制方程、湍流模型及湍流的先进壁面函数边界条件,介绍流动的数值求解方法.对典型超声速层流和湍流流动的摩擦阻力和热流进行详细的验证与确认,考察CFD工具在使用先进壁面函数边界条件后,湍流计算的法向网格无关性能力.对设计的一种吸气式高超声速飞行器的气动力和气动热进行数值模拟,为飞行器的气动设计及热防护提供了可靠的数据.
相关文章 | 多维度评价
一种新的壁面距离计算方法——循环盒子法
赵慧勇, 贺旭照, 乐嘉陵
计算物理    2008, 25 (4): 427-430.  
摘要351)      PDF (304KB)(1197)   
很多湍流模型都需要内点到壁面的距离,当网格量特别大时,计算壁面距离不仅计算量很大,而且费时.通过对当前壁面距离计算方法的分析,提出一种新的计算方法——循环盒子方法.循环盒子法计算效率高,通用性和健壮性较好,可以得到和直接法相同的结果.通过对7个二维和三维算例的计算表明,循环盒子法的计算时间比直接法缩小一个量级.而且,对于同一个问题,随着网格的加密,循环盒子法效率更高.
相关文章 | 多维度评价
双时间步方法的应用分析
赵慧勇, 乐嘉陵
计算物理    2008, 25 (3): 253-258.  
摘要519)      PDF (286KB)(1239)   
对双时间步方法隐式迭代的稳定性、启动问题、子迭代初值和收敛性准则进行分析.一般认为,双时间步的真实时间步长是基于精度,而不是稳定性的基础上给定,因此真实时间步长可以取得很大,但对VonNeumann稳定性分析表明,对于子迭代采用隐式的算法,稳定性的要求对真实时间步有限制.并通过对Sod激波管问题的计算,验证该分析.
相关文章 | 多维度评价
用遗传算法拟合Eu-151穆斯堡尔谱
孙雪峰, 林俊, 黄红波, Ali, 何云, 戴耀东, 夏元复
计算物理    2003, 20 (3): 264-266.  
摘要313)      PDF (195KB)(1016)   
以Eu2O3(立方晶系)为例,论述遗传算法(GA,Genetic Algorithm)求解穆斯堡尔谱的过程,并将解谱结果与文献数据和Moss Win解谱结果进行比较.
相关文章 | 多维度评价
一维波动方程反问题求解的正则化方法
吴建成, 张大力, 刘家琦
计算物理    1995, 12 (3): 415-420.  
摘要314)      PDF (316KB)(1115)   
讨论了一维波动方程utt-∂x(μ(x)ux)=0在一般的初、边值条件和附加条件下系数μ(x)的求解方法.把反问题归结为一不适定的非线性积分方程组,利用正则化方法克服了反问题的不适定性.
相关文章 | 多维度评价
强激光照射下受预内压作用的球壳的热力响应
陈发良
计算物理    1994, 11 (3): 269-277.  
摘要237)      PDF (510KB)(936)   
分析受强激光照射和预内压联合作用的球壳的热力效应。在光强相对较低的情形,激光照射处理为壳体表面热流输入。推导了瞬态热传导解析解;并选取一组典型的几何及材料参数,针对一系列激光光强情形用大型有限元程序对热传导和热弹塑性变形作了数值计算,得到了热传导的清晰的物理图象,讨论了壳体结构变形的空间分布和随时间的演变规律以及对激光光强和照射时间的依赖关系。
相关文章 | 多维度评价
TVD数值技术在弱电离化学非平衡斜激波反射流场计算中的应用
胡光初, 乐嘉陵, 曹文祥
计算物理    1994, 11 (2): 195-202.  
摘要323)      PDF (445KB)(904)   
本文成功地将二阶迎风TVD数值格式运用到弱电离、化学非平衡斜激波反射流场的数值计算中。气体力学方程与组元质量守恒方程组全耦合求解是目前非平衡流场数值计算中的新方向,文中详细地给出了气体力学方程和组元质量守恒方程组全耦合求解过程,最后给出了三种不同情况斜激波反射流场数值结果并进行了详细讨论,本文所得结果与文献[4]由I.I.Glass提供的实验结果和文献[1]由H.M.Glaz提供的数值结果相符合。
相关文章 | 多维度评价