全自动超声速流动计算架构及实现

1

作者 刘君 陈洁 +2 位作者 卢俊宇 汪骥 徐春光 《国防科技大学学报》 北大核心 2025年第1期83-93,共11页

实现CFD计算流程中网格生成及计算流程的自动化,能显著提高CFD仿真效率,具有重要应用价值。以扎染算法为基础的全自动CFD模拟技术架构,实现输入实体模型后无须几何清理,即可自动生成网格、快速开始无黏超声速流场计算的功能。通过二维... 实现CFD计算流程中网格生成及计算流程的自动化,能显著提高CFD仿真效率,具有重要应用价值。以扎染算法为基础的全自动CFD模拟技术架构,实现输入实体模型后无须几何清理,即可自动生成网格、快速开始无黏超声速流场计算的功能。通过二维数值算例验证了该技术的计算精度和非结构有限体积法相当,但计算效率明显提升;开发出基于手绘模型实时开始计算的Auto-CFD软件,理论上能对二维任意外形进行全自动模拟计算;将该技术从二维空间推广到三维,以激光扫描汽车得到的不规则点云作为外形演示了本技术对复杂外形问题的适应能力。基于扎染算法发展的Auto-CFD技术架构,兼容主流差分格式,具有良好的网格适应性,有望解决现有Auto-CFD软件存在的技术问题。 展开更多

关键词 非结构网格有限差分法 笛卡儿网格 网格生成 软件自动化 扎染算法

在线阅读 下载PDF 职称材料

水下拖曳系统弱流固耦合方法

2

作者 王志博 黄帅瑜 +2 位作者 顾津菁 赵淼 姚惠之 《哈尔滨工程大学学报》 北大核心 2025年第3期410-417,共8页

为了解决拖曳缆系统采用常系数估算流体载荷造成缆索动力学响应模拟不准确的问题,本文将有限差分法与雷诺平均方程相结合,依托初始缆形采用网格生成器脚本功能生成计算域,设计拓扑源面,并填充结构化网格,以高质量黏性网格为基础,采用计... 为了解决拖曳缆系统采用常系数估算流体载荷造成缆索动力学响应模拟不准确的问题,本文将有限差分法与雷诺平均方程相结合,依托初始缆形采用网格生成器脚本功能生成计算域,设计拓扑源面,并填充结构化网格,以高质量黏性网格为基础,采用计算流体力学求解器结合湍流模型计算得到沿缆分布的水动力,采用插值方法将水动力载荷插值到离散缆段几何中心处,利用有限差分方法更新拖曳系统动力学计算结果获得新的缆形,循环迭代地建立了单向耦合的弱流固耦合方法。计算结果表明:该耦合方法具有良好的收敛性和稳定性,计算结果与实测缆形更加吻合。该方法也获得了拖曳缆的绕流场流动分离特点和流体载荷沿着缆的分布的一般规律,证实了沿缆的阻力系数取固定的经验系数会产生较大误差。 展开更多

关键词 拖曳系统 流固耦合 计算流体力学 有限差分方法 雷诺平均方程 插值 粘性网格 阻力

在线阅读 下载PDF 职称材料

油藏组分模型的无网格扩展有限体积法

3

作者 饶翔 徐云峰 +2 位作者 刘伟 周玉辉 刘怡娜 《计算力学学报》 CAS CSCD 北大核心 2024年第4期709-717,748,共10页

本文开发了油藏组分模型的首个无网格数值求解器。该无网格求解器使用生成时受拓扑约束较少的点云和扩展有限体积法分别离散油藏计算域和控制方程,可以直接使用网格类模拟器中已有的非线性求解器计算全局非线性离散方程组而获取油藏计... 本文开发了油藏组分模型的首个无网格数值求解器。该无网格求解器使用生成时受拓扑约束较少的点云和扩展有限体积法分别离散油藏计算域和控制方程,可以直接使用网格类模拟器中已有的非线性求解器计算全局非线性离散方程组而获取油藏计算域压力、相饱和度和组分浓度的分布。采用两个典型的油藏组分模型算例(包括CO_(2)驱和蒸汽驱)用以检验该无网格模拟器的计算表现。算例结果表明,相较于网格类方法,该无网格模拟器能够显著降低对复杂油藏计算域的离散难度,且能避免笛卡尔网格有限体积法的网格取向效应,提高了计算精度和非线性求解的收敛效率。 展开更多

关键词 广义有限差分方法 无网格扩展有限体积法 油藏组分模型 油藏数值模拟

在线阅读 下载PDF 职称材料

先进中子学栅格程序KYLIN-2的开发与初步验证 被引量：10

4

作者 柴晓明 涂晓兰 +7 位作者 郭凤晨 尹强 黄世恩 芦韡 卢宗健 姚栋 李庆 吴文斌 《强激光与粒子束》 EI CAS CSCD 北大核心 2017年第1期90-96,共7页

针对先进核反应堆中结构复杂的燃料组件,中国核动力研究设计院开发了先进中子学栅格中子学程序KYLIN-2,该程序采用子群方法进行共振处理,采用特征线方法(MOC)进行复杂几何中子输运计算,并利用采用广义粗网格有限差分加速方法(GCMFD)来... 针对先进核反应堆中结构复杂的燃料组件,中国核动力研究设计院开发了先进中子学栅格中子学程序KYLIN-2,该程序采用子群方法进行共振处理,采用特征线方法(MOC)进行复杂几何中子输运计算,并利用采用广义粗网格有限差分加速方法(GCMFD)来加速中子输运求解流程,采用基于改进预估校正临界-燃耗迭代方法(PPC)的切比雪夫方法求解复杂燃耗链,同时,为了方便用户使用,开发形成了支持复杂组件图形化建模工具和后处理显示工具。通过初步的数值检验,针对典型压水堆燃料组件,KYLIN-2具有较高的计算精度,满足工程使用需求。 展开更多

关键词 中子输运方程 特征线方法 子群方法 广义粗网格有限差分加速方法 切比雪夫燃耗方法

在线阅读 下载PDF 职称材料

基于广义有限差分法的输流直管振动响应特性研究 被引量：10

5

作者 张挺 林震寰 +2 位作者 郭晓梅 张恒 范佳铭 《振动与冲击》 EI CSCD 北大核心 2019年第24期165-171,共7页

基于输流直管横向振动微分方程,采用广义有限差分法和Houbolt法分别对空间和时间上的偏微分项进行离散,建立高阶精度的无网格法数值模式。将该数值模式应用于两端支撑轴向运动梁模型和两端固支输流直管模型,所得的固有频率和振幅时程等... 基于输流直管横向振动微分方程,采用广义有限差分法和Houbolt法分别对空间和时间上的偏微分项进行离散,建立高阶精度的无网格法数值模式。将该数值模式应用于两端支撑轴向运动梁模型和两端固支输流直管模型,所得的固有频率和振幅时程等与前人研究成果和理论解均吻合良好;并对不同总点数N、时间步长△t、子区域选点数n s的数值结果进行对比,表明提出的数值模型在求解输流直管振动响应问题上具有良好的准确性和鲁棒性;对比分析了三种不同支撑条件(两端固支、两端简支和一端固支一端简支)下输流直管的振动响应特性。结果表明:两端简支时输流直管中点处的振幅最大,振动频率最小;两端固支时输流直管中点处的振幅最小,振动频率最大;且在端部约束限制条件不对称时,其振动幅值最大值出现位置会向弱约束端偏移。 展开更多

关键词 输流直管 无网格法 广义有限差分法(gfdm) Houbolt法 横向振动

在线阅读 下载PDF 职称材料

铸造CAE系统的网格剖分技术 被引量：5

6

作者 刘晶峰 李洪友 方建成 《华侨大学学报（自然科学版）》 CAS 北大核心 2008年第3期327-330,共4页

基于有限差分法和标准模板库(STL),研究铸造计算机辅助工程(CAE)系统中的网格剖分关键技术,具体包括STL实体截面轮廓快速生成、容错处理,以及截面网格信息的快速提取等技术.由于实现了容错处理技术,对于既定的铸件/铸型系统,其对应的ST... 基于有限差分法和标准模板库(STL),研究铸造计算机辅助工程(CAE)系统中的网格剖分关键技术,具体包括STL实体截面轮廓快速生成、容错处理,以及截面网格信息的快速提取等技术.由于实现了容错处理技术,对于既定的铸件/铸型系统,其对应的STL数据文件即使存在一定的错误也可以顺利实施网格剖分,这就充分保证网格剖分的安全性和稳定性.特别地,由于采用了较为简单的数据结构及简洁的算法,所以确保网格剖分程序的高效率运行.最后,利用上述网格剖分技术成功实现三维铸件/铸型系统的网格剖分. 展开更多

关键词 网格剖分 截面轮廓 容错处理 铸造CAE 有限差分法

在线阅读 下载PDF 职称材料

广义α算法在拖曳阵列动态仿真中的应用 被引量：4

7

作者 张广磊 张维竞 刘涛 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2012年第5期574-579,共6页

针对有限差分Box method时间域上变量离散求解非线性拖缆动力学方程存在不稳定的性缺点,在求解拖缆方程时在时间域上采用广义α算法,并保留在空间域上有限差分Box method的离散方法.将新算法应用于三维非线性拖缆动态方程仿真求解,通过... 针对有限差分Box method时间域上变量离散求解非线性拖缆动力学方程存在不稳定的性缺点,在求解拖缆方程时在时间域上采用广义α算法,并保留在空间域上有限差分Box method的离散方法.将新算法应用于三维非线性拖缆动态方程仿真求解,通过对目前几种常用算法的仿真实验对比分析,表明新算法不仅具有较高的计算精度,而且计算时间短、稳定性好,为求解非线性拖缆方程提供了一个较好的求解方法,可以应用到实际应用中,对准确预报拖曳阵列在各种状态下的运动规律具有重要意义. 展开更多

关键词 拖缆 非线性 有限差分 时域 广义α算法

在线阅读 下载PDF 职称材料

基于广义有限差分法求解非线性自由液面的液体晃动问题 被引量：4

8

作者 张挺 任聿飞 +1 位作者 杨志强 范佳铭 《四川大学学报（工程科学版）》 EI CAS CSCD 北大核心 2016年第1期8-14,共7页

广义有限差分法属新发展的区域型无网格法,在每个演算时间层内可避免重建网格及数值积分等工作,可用于模拟非线性动态边界的物理问题。针对2维矩形水槽的自由液面晃动问题,以Lagrangia法为基础,采用广义有限差分法配合蛙跳法对于此移动... 广义有限差分法属新发展的区域型无网格法,在每个演算时间层内可避免重建网格及数值积分等工作,可用于模拟非线性动态边界的物理问题。针对2维矩形水槽的自由液面晃动问题,以Lagrangia法为基础,采用广义有限差分法配合蛙跳法对于此移动边界问题分别在空间上和时间上进行离散。通过水槽液体自由晃动和受迫振动的数值案例对比,结果表明提出的数值模式在处理自由液面的液体晃动问题的准确性、有效性以及稳定性。 展开更多

关键词 液体晃动 广义有限差分法 无网格法 拉格朗日法 蛙跳法

在线阅读 下载PDF 职称材料

共享内存系统中高效并行FDTD计算方案 被引量：3

9

作者 赖生建 王秉中 黄廷祝 《电子科技大学学报》 EI CAS CSCD 北大核心 2010年第5期680-683,共4页

针对电大问题的FDTD仿真,在共享内存系统中提出一种不交换数据的并行FDTD计算方案。该方案利用共享内存系统的通用多核处理器架构,直接读取并行场域边界面的数据实现并行计算,采用自主开发的多线程技术实现FDTD的并行计算。通过数值仿... 针对电大问题的FDTD仿真,在共享内存系统中提出一种不交换数据的并行FDTD计算方案。该方案利用共享内存系统的通用多核处理器架构,直接读取并行场域边界面的数据实现并行计算,采用自主开发的多线程技术实现FDTD的并行计算。通过数值仿真测试,提出了最大并行计算效率的区域划分方式,并验证了所提出的方案是一种理想的基于共享内存系统的并行FDTD计算方案。 展开更多

关键词 计算效率 时域有限差分法 通用多核处理器 并行算法 共享内存

在线阅读 下载PDF 职称材料

水下滑坡产生涌浪波特性的数值模拟研究 被引量：5

10

作者 荆海晓 陈国鼎 李国栋 《应用力学学报》 CAS CSCD 北大核心 2018年第3期503-509,共7页

通过在浅水方程中引入滑块运动项来考虑滑坡涌浪波的产生过程,采用四阶预测矫正格式和高阶有限差分格式对扩展浅水方程进行了数值求解,建立了一种基于浅水方程的滑坡涌浪有限差分模型。采用已有的实验结果与所建立模型的数值结果进行对... 通过在浅水方程中引入滑块运动项来考虑滑坡涌浪波的产生过程,采用四阶预测矫正格式和高阶有限差分格式对扩展浅水方程进行了数值求解,建立了一种基于浅水方程的滑坡涌浪有限差分模型。采用已有的实验结果与所建立模型的数值结果进行对比验证。在此基础上对平底水下滑坡运动情况,不同滑块长度、最大厚度、运动速度下产生的涌浪波首波特性进行了数值计算分析。结果表明:本文模型模拟结果与实验结果符合程度较好,其决定系数值均达到0.89以上,表明本文模型具有较高精度;涌浪波波长的主要影响因素为滑块长度,随着滑块长度的增加,涌浪波波长随之增加,而滑块运动速度和滑块厚度对生成的涌浪波波长无明显影响;涌浪波波高主要受滑块运动速度(滑块Fr)和滑块厚度影响,随着滑块厚度的增加,波高随之增加。当滑块Fr小于1时,首波波高随滑块Fr的增大而增大,而当滑块Fr大于1时,首波波高随滑块Fr的增大而减小,首波波高在滑块Fr等于1附近达到最大值;涌浪波波面变形主要受滑块厚度和滑块速度影响,滑块厚度越大、滑块Fr越接近1,涌浪波波面变形越明显。 展开更多

关键词 滑坡涌浪 浅水方程 有限差分法 滑块Fr 水下滑坡

在线阅读 下载PDF 职称材料

基于广义有限差分法的多孔介质墙体热湿耦合迁移模型求解 被引量：1

11

作者 郭兴国 罗银辉 刘向伟 《南昌大学学报（工科版）》 CAS 2023年第3期242-247,共6页

针对多孔介质墙体的热湿耦合迁移控制方程的求解问题,提出了一种基于广义有限差分法的新型无网格区域法数值模型。该模型在空间域上采用广义有限差分法离散,在时间域上采用隐式欧拉法离散,得到递归形式的线性方程组并用Matlab软件进行... 针对多孔介质墙体的热湿耦合迁移控制方程的求解问题,提出了一种基于广义有限差分法的新型无网格区域法数值模型。该模型在空间域上采用广义有限差分法离散,在时间域上采用隐式欧拉法离散,得到递归形式的线性方程组并用Matlab软件进行程序编写和运算。将计算结果与采用解析解、有限差分法、TDMA法和实验的结果进行比较,验证解法的准确性。分析了区域总点数、子区域选点数以及时间步长对数值计算精度的影响,结果表明该数值方法在求解热湿耦合迁移微分方程上具有良好的稳定性。 展开更多

关键词 热湿耦合迁移 广义有限差分法 无网格法

在线阅读 下载PDF 职称材料

微扰Landau-Ginzburg-Higgs方程的保结构数值分析 被引量：1

12

作者 胡伟鹏 张宇 邓子辰 《西北工业大学学报》 EI CAS CSCD 北大核心 2012年第6期957-960,共4页

基于Hamilton变分原理,构造了微扰Landau-Ginzburg-Higgs方程的一阶广义多辛对称形式,随后对该形式采用多辛差分离散构造其保结构离散格式,最后通过计算机模拟,研究了微扰对Lan-dau-Ginzburg-Higgs方程孤子解的影响,为微扰动力学系统的... 基于Hamilton变分原理,构造了微扰Landau-Ginzburg-Higgs方程的一阶广义多辛对称形式,随后对该形式采用多辛差分离散构造其保结构离散格式,最后通过计算机模拟,研究了微扰对Lan-dau-Ginzburg-Higgs方程孤子解的影响,为微扰动力学系统的数值研究提供了新的途径。 展开更多

关键词 有限差分方法 哈密尔顿 孤子解 广义多辛 微扰Landau—Ginzburg-Higgs方程 保结构

在线阅读 下载PDF 职称材料

多模光纤光栅中新型的光波传输分析法 被引量：1

13

作者 王云明 戴劲草 +1 位作者 张明德 孙小菡 《应用科学学报》 CAS CSCD 2004年第1期25-29,共5页

提出一种新型的分析多模光纤光栅的准三维时域有限差分束传播法(Q-3D-FDTD-BPM).在圆柱坐标系下将场在角向分量的解析式分离出来,使三维方程变为若干个二维方程求解,从而简化了三维离散化方程.同时使用综合道格拉斯(GD)法和交替隐向迭代... 提出一种新型的分析多模光纤光栅的准三维时域有限差分束传播法(Q-3D-FDTD-BPM).在圆柱坐标系下将场在角向分量的解析式分离出来,使三维方程变为若干个二维方程求解,从而简化了三维离散化方程.同时使用综合道格拉斯(GD)法和交替隐向迭代法(ADIM),获得比中心差分算子高两个数量级的空间离散精度.然后根据场强加权转换方案,加权综合各子方程解,从而得到多模光纤中三维场解.基于这种方法,模拟多模光纤光栅中场分布,得到与实验结果多模式特性基本一致的光栅传输光谱图. 展开更多

关键词 多模光纤光栅 光波传输 时域有限差分束传播法 综合道格拉斯法 交替隐向迭代法

在线阅读 下载PDF 职称材料

探针馈电圆柱共形微带天线阻抗特性的FDTD法分析 被引量：2

14

作者 何芒 徐晓文 《微波学报》 CSCD 北大核心 2004年第3期1-5,共5页

将圆柱坐标系下非分裂式理想匹配层吸收边界条件 (UPML)引入到圆柱共形微带天线的全波分析中 ,并给出了其中的场量迭代方程 ,明显地减小了计算空间 ,简化了编程 ;提出了圆柱坐标系下考虑有限尺寸探针半径的新算法 ,使得计算结果更加精... 将圆柱坐标系下非分裂式理想匹配层吸收边界条件 (UPML)引入到圆柱共形微带天线的全波分析中 ,并给出了其中的场量迭代方程 ,明显地减小了计算空间 ,简化了编程 ;提出了圆柱坐标系下考虑有限尺寸探针半径的新算法 ,使得计算结果更加精确 ;应用GPOF方法预估了时域晚时响应 ,使计算时间减少了 70 %～ 90 %。本文给出的模型能够对任意多层圆柱共形结构的微带天线进行计算 ,因此对此类天线的CAD设计具有实际意义。 展开更多

关键词 微带天线 FDTD法 馈电 全波分析 阻抗特性 理想匹配层 吸收边界条件 圆柱坐标系 迭代 探针

在线阅读 下载PDF 职称材料

基于第二代小波的PDEs数值新算法 被引量：2

15

作者 胡俊华 田锦会 《信息技术》 2013年第9期88-90,95,共4页

根据第二代小波变换计算效率高和易实现自适应的特性,将微分方程的数值解映射到小波系数空间,参考小波系数选择合适的阈值,构造自适应网格。在大于阈值的小波系数所对应的区域内,适当增加小波函数的尺度,提高基函数的分辨率,因而能捕捉... 根据第二代小波变换计算效率高和易实现自适应的特性,将微分方程的数值解映射到小波系数空间,参考小波系数选择合适的阈值,构造自适应网格。在大于阈值的小波系数所对应的区域内,适当增加小波函数的尺度,提高基函数的分辨率,因而能捕捉到数值解的突变结构;同时,将小波逆变换引入到方程数值解的逼近过程中,避免了直接在小波基上的运算,提高了算法的执行效率。在自适应网格上关于空间变量的导数计算采用变步长的有限差分法,在时间方向上用向前差分进行迭代求解。 展开更多

关键词 第二代小波变换 自适应网格 有限差分法 非线性偏微分方程

在线阅读 下载PDF 职称材料

多重网格方法在风电机组流场数值模拟中的应用

16

作者 王金铭 包常利 +1 位作者 卢奭瑄 霍丹 《太阳能学报》 EI CAS CSCD 北大核心 2012年第11期2005-2009,共5页

针对风电机组流场的数学模型为非定常、非线性偏微分方程组的定解问题将多重网格方法与有限差分法、牛顿迭代法相结合,应用到风电机组流场的数值模拟问题中,构造了多重网格方法的限制算子、插值算子及光滑迭代,给出一种具有快速收敛的... 针对风电机组流场的数学模型为非定常、非线性偏微分方程组的定解问题将多重网格方法与有限差分法、牛顿迭代法相结合,应用到风电机组流场的数值模拟问题中,构造了多重网格方法的限制算子、插值算子及光滑迭代,给出一种具有快速收敛的有限差分方法、牛顿迭代法与多重网格方法相结合的数值算法。数值算例表明该文所给数值方法的计算时间比传统迭代法减少50%以上,是风电机组流场数值模拟的一种高效、实用的数值计算方法。 展开更多

关键词 风电机组 流场 多重网格方法 有限差分法 牛顿迭代法 数值模拟

在线阅读 下载PDF 职称材料

一般正交曲线坐标系中的降维时域有限差分法

17

作者 王建永 赵长青 张敏 《微波学报》 CSCD 北大核心 2007年第3期17-19,共3页

相对于常规的时域有限差分法(FDTD),降维时域有限差分法(R-FDTD)减少了约1/3的存储量。本文推导得出了一般正交曲线坐标系(GOC)中的降维时域有限差分法的基本公式,给出了其计算步骤,为在一般正交曲线坐标系中应用降维时域有限差分法奠... 相对于常规的时域有限差分法(FDTD),降维时域有限差分法(R-FDTD)减少了约1/3的存储量。本文推导得出了一般正交曲线坐标系(GOC)中的降维时域有限差分法的基本公式,给出了其计算步骤,为在一般正交曲线坐标系中应用降维时域有限差分法奠定了基础。数值实验证实了这一方法的有效性。 展开更多

关键词 降维时域有限差分法(R-FDTD) 一般正交曲线坐标系(GOC) 存储量

在线阅读 下载PDF 职称材料

基于超表面的宽波带光束聚焦研究 被引量：4

18

作者 胡晨曦 王吉明 +3 位作者 吴彤 赫崇君 顾晓蓉 刘友文 《激光技术》 CAS CSCD 北大核心 2018年第5期681-686,共6页

超颖表面是一种基于亚波长结构的光学平板膜层,可在亚波长传输范围内调控入射光束的相位、振幅和偏振。为了代替传统的曲面光学元件,采用传输相位调控理论和广义折反射定律,设计了一种新型的超颖表面,并进行了程序模拟,取得了此亚波长... 超颖表面是一种基于亚波长结构的光学平板膜层,可在亚波长传输范围内调控入射光束的相位、振幅和偏振。为了代替传统的曲面光学元件,采用传输相位调控理论和广义折反射定律,设计了一种新型的超颖表面,并进行了程序模拟,取得了此亚波长结构对光束聚焦调控的数据。结果表明,当增加超颖表面的椭圆基元的长短轴长度时,材料的等效折射率增加,并且适用的波长范围增加到0.7μm^1.2μm;通过优化超表面结构参量,可实现在宽波带范围内的相位调控,进而获得聚焦光场的优化,在一定程度上可以代替传统光学元件实现光学聚焦。该研究结果在超分辨率成像及光刻等方面有一定参考价值,在一些特殊的需要亚波长结构调控光束的情况下可以使光路简单化,并且比传统的光学元件有着厚度方面的优势。 展开更多

关键词 光学器件 超颖表面 时域有限差分法 广义折射和反射定律

在线阅读 下载PDF 职称材料

含肿瘤皮肤组织传热分析的广义有限差分法 被引量：9

19

作者 李艾伦 傅卓佳 +1 位作者 李柏纬 陈文 《力学学报》 EI CSCD 北大核心 2018年第5期1198-1205,共8页

生物传热分析在低温外科手术、肿瘤热疗、病热诊断等临床医学治疗和诊断中有着广泛的应用.由于健康皮肤组织内肿瘤的存在使得肿瘤附近区域的温度会明显升高,这一特性常被用于检测皮肤组织内的肿瘤生长,因此有必要开展生物传热数值分析... 生物传热分析在低温外科手术、肿瘤热疗、病热诊断等临床医学治疗和诊断中有着广泛的应用.由于健康皮肤组织内肿瘤的存在使得肿瘤附近区域的温度会明显升高,这一特性常被用于检测皮肤组织内的肿瘤生长,因此有必要开展生物传热数值分析的研究.本文以含肿瘤的皮肤组织为研究对象,将一种新型区域型无网格配点法——广义有限差分法应用于能描述含肿瘤皮肤组织传热过程的Pennes方程求解.广义有限差分法利用泰勒展开式与移动最小二乘法将计算区域内的每个离散点上的物理量导数表示成其与邻近点物理量及权重系数的线性组合,进而构建得到仅含各离散点未知物理量的线性方程组.该方法不仅具有无需划分网格、避免数值积分等无网格配点法的优点,同时还克服了大多数无网格配点法中插值矩阵高度病态的问题,为此类方法在大规模工程数值计算中的应用提供了可能性.本文首先介绍了模拟含肿瘤皮肤组织传热过程的广义有限差分法离散模型,随后通过不含肿瘤与含规则形状肿瘤的基准算例,检验广义有限差分法的计算精度与收敛性;在此基础上,通过数值模拟研究不同肿瘤形状及肿瘤位置分布对皮肤组织内温度分布的影响. 展开更多

关键词 广义有限差分法 配点法 移动最小二乘法 泰勒展开式 Pennes生物传热方程

在线阅读 下载PDF 职称材料

基于广义Kelvin模型的非定常盐岩蠕变模型 被引量：10

20

作者 韩伟民 闫怡飞 闫相祯 《中南大学学报（自然科学版）》 EI CAS CSCD 北大核心 2020年第5期1337-1349,共13页

针对川东北油气田盐岩地层普遍出现的套损现象,对该地层段盐岩的蠕变特性进行研究。用非定常黏壶元件替换广义Kelvin模型中的线性黏壶元件,建立非定常广义Kelvin模型;将其与能够描述盐岩稳态蠕变特征的Heard模型进行串联,构建新的四元件... 针对川东北油气田盐岩地层普遍出现的套损现象,对该地层段盐岩的蠕变特性进行研究。用非定常黏壶元件替换广义Kelvin模型中的线性黏壶元件,建立非定常广义Kelvin模型;将其与能够描述盐岩稳态蠕变特征的Heard模型进行串联,构建新的四元件NGKH盐岩非线性蠕变本构模型,可以描述盐岩的衰减蠕变阶段和稳态蠕变阶段特征。基于室内三轴蠕变试验结果,利用Levenberg-Marquardt算法对NGKH模型蠕变参数进行识别。基于FLAC3D的UDM接口程序,采用C++编程语言对其进行二次开发;并通过三轴蠕变试验的数值模拟对开发NGKH模型进行验证。研究结果表明:NGKH模型的拟合结果与室内试验结果吻合良好,该蠕变模型能够准确描述盐岩的非线性蠕变特性;通过黏弹性与非线性计算结果,验证了开发NGKH模型程序的正确性和合理性。 展开更多

关键词 蠕变本构模型 非定常广义Kelvin模型 Heard稳态蠕变模型 二次开发 有限差分法

在线阅读 下载PDF 职称材料