共查询到20条相似文献,搜索用时 31 毫秒
1.
暂态历程的精细计算方法 总被引:6,自引:0,他引:6
在2 ̄N类算法计算指数矩阵的基础上,提出了精细积分法计算暂态历程问题,其数值结果可以比拟于精确解的数值结果。数值例题表明了方法高度准确的特点。分析了算法的精度,指出了高度准确的条件。 相似文献
2.
数值求解非均匀介质中的输运问题广泛应用于科学计算和工程领域.介质的强非均匀性给相关问题的准确求解带来极大的困难.近年来,本课题组将有限分析法拓展到该领域,建立了非均匀介质中输运问题的有限分析法.该算法基于网格奇点邻域内类拉普拉斯方程局部解析解构建,算法具有很高的精度,且不依赖于介质的非均匀性强度.不管相邻网格传导率差异如何,仅需对原始网格进行很少地细分就可以获得非常准确的计算结果,因此与其他传统数值算法相比,可以大幅提高计算精度和效率.该算法可广泛应用于求解非均匀多孔介质中的渗流、复合材料中的热传导及电场分布等问题. 相似文献
3.
根据泥质夹层的低渗特性及空间分布,本文提出了一种含泥质夹层油藏网格渗透率的粗化计算方法,并在此基础上,将自适应网格算法应用于含泥质夹层油藏的数值模拟,提升其计算效率.在计算过程中,网格的动态划分仅依据流体物理量的变化,泥质夹层区域不全部采用细网格,仅针对流动锋面处的泥质夹层采用细网格,其余泥质夹层处采用不同程度的粗网格.相较于传统算法,网格数大幅下降.数值算例表明,自适应网格算法的计算结果精度与全精细网格一致,能够准确模拟出泥质夹层对于流体的阻碍作用,同时计算效率得到大幅提升,约为全精细网格算法的3~7 倍. 相似文献
4.
5.
非线性问题和分叉问题及其数值方法 总被引:5,自引:0,他引:5
本文给出了一个一般性的分叉定义,说明了伪弧长算法在分叉计算中的应用,概述了静分叉点定位、用单纯形算法准确确定静分叉后各分叉解枝初始方向的算法,以及Hopf分叉点定位和大范围连续追踪周期解轨道的数值方法。 相似文献
6.
在IPSAR算法基础上推导了适用于考虑液度变化影响的密相液固两相湍流流动数值计算的DIPSAR算法,并采用低雷诺数模型,对竖直上升管中密相液固两相湍流分别采用IPSAR算法和DIPSAR算法进行了数值计算,计算值与实测值较符合,计算结果的比较表明DIPSAR算法能更有效地预测密相液固两相湍流流动 相似文献
7.
三维编织复合材料模量的双尺度有限元计算 总被引:4,自引:1,他引:4
针对三维编织复合材料的力学性能进行了双尺度有限元(TSA)数值计算,给出了计算模型和算法过程,并将数值结果与文献中的实验数据进行了比较,验证了算法的物理准确性。编织复合材料的力学性能不仅依赖于材料的基本组份,也与细观构造相关。双尺度有限元计算可以数值模拟出三维编织复合材料的整体力学性能,从而为材料的研发提供指导。本文的双尺度有限元三维数值计算方法可以推广到其他增强/孔隙等多相复合材料的数值模拟。 相似文献
8.
随着磁头滑块的飞行高度不断降低,给气体润滑方程的数值求解带来了诸如计算时间过长、甚至计算发散等方面的问题。为了获得1Tbit/in2的存储密度,磁头滑块尾部的最小飞行高度接近1.5nm。本文基于作者提出的修正气膜润滑方程的线性流率(LFR)模型,考虑磁头滑块表面高度的不连续性,建立了基于有限体积法的气膜润滑方程离散格式,并把网格自适应技术与多重网格法应用到离散方程的迭代算法中,发展了可模拟最小飞行高度为0.5nm时磁头滑块压力分布的数值模拟方法与有效算法。文中以一个具有复杂表面形状的磁头滑块为例,检验了计算方法与算法的有效性。数值结果表明:在磁头滑块最小飞行高度较低时,必须要考虑滑块表面高度的不连续性,否则就得不到收敛的数值计算结果;与FK-Boltzmann模型相比,LFR模型具有较高的计算效率,采用网格自适应技术与多重网格法能有效地提高求解气膜润滑方程的计算效率。 相似文献
9.
分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于
某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于
其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳
分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于
某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于
其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳
定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法,它能在
保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实
了此特点,该方法被有效应用于充填流动过程的数值模拟. 相似文献
10.
11.
提出了一种基于节点分离概念的Lagrange有限元冲击破碎分析的算法,采用节点分离技术实现
了网格断裂,并用畸变侵蚀技术处理严重畸变单元。利用C++编程实现了节点分离计算模型的创建和畸变
侵蚀,结合LS-dyna的重启动分析功能,编程调用LS-dyna求解器以及畸变侵蚀程序实现了破碎分析。分别
用节点分离的Lagrange有限元算法、SPH 无网格算法对超高速碰撞问题进行了数值模拟,并与文献中的实
验结果和Euler有限元数值模拟结果进行了对比。结果表明,节点分离算法具有计算速度快、数值稳定、边界
明确等优点,能准确有效模拟超高速碰撞问题。 相似文献
12.
成形过程数值模拟的张量时空求解策略 总被引:1,自引:0,他引:1
成形过程的数值模拟和分析涉及几何、材料、接触摩擦等高度非线性的耦合作用,采用传统的增量算法会导致巨大计算量。本文建议采用张量时间函数的非增量时空算法,在整个时间和空间域上迭代求解。由于采用新的分离变量构思,以张量表达时间函数,可提高问题求解的速度和精度。但对问题的列式和数值求解方法提出了更高的要求。文中讨论其方法与数值实施。 相似文献
13.
<正> 3 新颖的算法 3.1 数值解的范围 要想使复杂流场的准确数值仿真成功,计算方面的处理方法应该是既强调对流体基本物理规律和力学规律的理解,也注意建立它们的准确模型.计算模型必须基于最少的一组过程和机理,人们认为这些模型能给出已知的分析解和观测结果的主要特点.对每一个主要过程 相似文献
14.
基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。 相似文献
15.
为了验证某型大展弦比电动飞机气动设计的合理性以及为飞行性能及品质计算提供数据,采用有限体积法离散求解三维可压雷诺平均Navier-Stokes方程,并选用Spalart-Allmaras湍流模型对该电动飞机流场进行CFD数值模拟。结果表明,该型电动飞机气动设计合理,巡航速度升阻比最高可达23,具有较高气动效率;通过CFD数值模拟得到了全机升力系数、阻力系数和升阻比。为了验证CFD计算结果,对该型电动飞机进行了缩比模型风洞实验,结果显示,CFD数值模拟法计算结果与风洞实验结果高度吻合,说明CFD计算结果准确。该方法可为大展弦比电动飞机气动设计提供指导。 相似文献
16.
17.
流体动力系统激波的数值模拟 总被引:1,自引:0,他引:1
本文为了对流体动力系统含有激波的跨音速流动进行数值模拟,给出了曲线坐标系下的LU-AUSMLW算法。该算法综合利用了AUSMPW格式与LU-SGS方法的优点。本文通过采用三阶MUSCL格式获得高阶精度。为了验证该算法的精度,对通道和叶栅中含有激波的跨音速流动进行了数值模拟。本文数值试验的计算结果与文献计算结果相符很好。 相似文献
18.
19.
涡量是流体运动的一个基本物理量,旋涡是流体运动中常见的一种基本形态,又是湍流的一种基本结构,它在理论上和工程实践中占有相当重要的位置。本文在吴镇远(J.C.Wu)方法的基础上,提出了一种新的数值计算方法——涡量修正方法。对炉内速度场及涡量场进行了计算,其速度场的计算结果与实验基本符合,在网格划分基本相同的条件下与SIMPLE算法所得结果相比,涡量修正算法更接近于实验结果,为流场的数值模拟提供了一条新的途径。 相似文献