首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
利用算子分裂迎风均衡格式解对流为主溶质运移问题   总被引:2,自引:0,他引:2  
水污染模拟问题是水流问题与溶质运移问题的耦合问题.各种常见的数值解法在以对流为主溶质运移问题的求解中都会遇到困难,如用有限单元法或有限差分法时,会产生数值弥散与过量这两类误差.引入算子分裂迎风均衡格式法求解对流为主的水污染模拟问题,较好地克服了数值弥散和数值解出现振荡问题,该格式具有良好的稳定性、单调性及守恒性特点.  相似文献   

2.
一维浅水流动方程的Godunov格式求解   总被引:16,自引:4,他引:16       下载免费PDF全文
以准确Riemann解为基础,建立了求解一维非平底浅水流动方程的Godunov格式,用"水位方程法(WaterLevel Formulation,WLF)"求解Riemann解,结合中心差分和Riemann解离散底坡项,保证了计算格式的和谐性.经算例验证,方法健全、通用,且分辨率高.  相似文献   

3.
基于有限体积方法和结构化网格,建立了海底滑坡引起的波浪传播数值模型。模型控制方程为考虑了海床随时间变化的二维浅水方程。采用中心迎风格式计算控制体界面数值通量,采用线性重构技术、局部海床高程处理技术和全隐式离散底摩阻项,保证了格式的和谐性、守恒性和水深非负性,有效处理了海岸动边界问题。时间积分采用具有强稳定性质的二阶龙格-库塔方法(ssp RK)。针对滑坡体兴波经典算例开展数值模拟,将计算结果与解析解、实验结果及其他模型计算结果进行比较和分析。结果表明,对于所考虑的计算工况,模型能较合理地模拟滑坡兴波的产生、传播和爬高过程。  相似文献   

4.
非结构网格上的三维浅水流动数值模型   总被引:7,自引:3,他引:7       下载免费PDF全文
针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。  相似文献   

5.
根据水质模型的具体特点,对不同的方程采用不同方法,水流问题用有限元法;对流弥散方程先用算子分裂的方法分解为两个方程,即对流方程和弥散方程,前者用高精度广义迎风格式求解,对弥散方程则采用多单元均衡格式法求解,最后合成为高精度广义迎风均衡格式求出溶质浓度。通过对数值实验例子的计算和实验溶质迁移的模拟,可以看出在求解对流弥散定解问题时,广义迎风均衡格式克服了有限元数值波动和浓度出现负值的问题,与有限元相比有较大改进。  相似文献   

6.
以三阶中心加权本质无振荡重构为基础,采用一维一维进行计算的方法,给出了求解二维浅水方程的高分辨率三阶半离散中心迎风格式.引入的重构方法既提高了格式的精度,又保证格式是无振荡的.时间的离散用最优的三阶SSP(Strong Stability Preserving)Runge-Kutta方法.源项的离散用辛普森公式.计算方法保持了中心差分格式简单的优点,即不需用黎曼解算器和进行特征分解过程.数值模拟结果与其它方法所得结果一致,表明了方法的有效性和稳定性.  相似文献   

7.
浅水方程数值计算方法的研究   总被引:1,自引:5,他引:1       下载免费PDF全文
求解以水位为变量的连续方程,并根据Navier-Stokes方程压力修正算法的基本思想,建立了浅水方程的水位修正算法,放宽了对离散时间步长的限制.通过对离散方程系数矩阵的重新构造,建立了高分辨率有限元格式,该格式既具有较高的离散精度又避免了数值解的伪振荡.对动量方程的阻力项做负坡线性化处理,提高了露滩计算的稳定性.数值模拟结果与解析解吻合良好,表明所建立的数值计算方法是正确的和可靠的.  相似文献   

8.
以Roe的近似Riemann解为基础,将源项按特征方向进行特征分解,建立了带源项浅水方程的通量平衡Godunov求解格式。此格式具有迎风的性质并保证了变宽、非平底坡浅水方程计算的和谐性和存在底摩擦时的收敛性。通过实例验证了此法具有和谐、健全、通用性好、分辨率高等优点。  相似文献   

9.
提高畦灌施肥地表水流与溶质运移数值模型的稳定性、收敛性及计算精度,有利于改善地面畦灌施肥系统的设计与评价工作.在利用隐-显混合时间格式对一维畦灌地表水流与溶质运移耦合模型中包含的各矢量项进行时间离散基础上,借助有限差分法、有限体积法和有限单元法分别对由隐时间格式生成的物理矢量线性近似式空间导数、物理矢量空间导数、溶质扩散矢量和地形矢量项进行空间离散,对最终形成的控制方程代数方程组进行数值求解,构建起基于混合数值解法的一维畦灌施肥地表水流与溶质运移耦合模型.  相似文献   

10.
苗雨  晏飞  万云冬  王元汉 《岩土力学》2008,29(7):1872-1876
杂交边界点法是一种边界类型的无网格法,它以修正变分原理和移动最小二乘近似为基础,同时具有边界元法和无网格法的优良特性。将该方法同双重互易法结合,将问题的解分为通解和特解两部分,通解使用杂交边界点方法求解,特解则利用局部径向基函数近似。该方法输入数据只是求解域上离散的点,不需要额外的方程来计算域内物理量,后处理十分简便。数值算例表明,该方法是一种计算精度很高的数值方法,适合于岩土工程的若干问题。  相似文献   

11.
泥石流是山区多发的一种地质灾害,它的发生和发展威胁着人们的生命和财产安全,影响着人们的正常生活,因而需要加强对其发生和发展过程的研究。结合泥石流的动力模型方程采用数值模拟方法再现泥石流发生和发展的过程,是研究和预测模拟泥石流灾害的有效手段。目前的动力模型方程大多只关注动力过程,却忽视了静动力过程的统一,这将导致在一些情况下产生错误的结果。本文研究了一维泥石流的静动力阻力特征,通过修正泥石流动力学方程的阻力项,得到了具有静动力统一特征的模型方程。并以Roe格式的近似Riemann解为基础,采用MUSCL线性重构方法建立了具有较高精度和分辨率的有限体积数值求解。具体算例的数值验证表明,方程阻力项的修正是合理的,所建立的数值求解也是稳定和有效的。  相似文献   

12.
格点法是在计算流体力学中首先发展起来的数值模拟新方法,共根本思想是对问题重新建模,建立直接模拟流体运动的离散格点模型。本文以一阶拟线性双典型方程及Kdv方程为例,推广这种求解一般偏微分方程的。它是运用多尺度分析方法,构造出格点模型的演化方程的局部平衡分布函数。数值试验表明,该方法程序实现简单,求解速度快,数值结果令人满意。  相似文献   

13.
突扩流动主流偏转特性数值模拟研究   总被引:2,自引:0,他引:2       下载免费PDF全文
以流函数方程作为控制方程,与N-S方程相比,不出现压强项和速度项,与涡量-流函数型N-S方程相比,没有涡量项和速度项,只有流函数项,避免了涡量边界条件问题,将控制方程中未知数的个数变为一个,通过解得流函数而得到其他各项。考虑对流项迎风性质,采用迎风差分格式,对突扩层流进行数值模拟,与DurstF流动试验及混合有限分析法模拟结果比较,基本吻合。通过一系列计算,得出主流经突扩断面偏转后相关特征值与雷诺数Re的关系,扩展了突扩流动的研究内容。  相似文献   

14.
为高效求解扩展型Boussinesq水波方程,建立了基于有限差分和有限体积方法的混合数值格式。将一维控制方程写为守恒形式,方程中通量部分采用有限体积方法求解,剩余部分采用有限差分方法求解。其中,有限体积方法采用Godunov类高分辨率格式,并结合HLL(Harten-Lax and van Leer)式黎曼问题近似解求界面数值通量,黎曼问题界面左右变量通过高精度状态插值方法(MUSCL)构筑。有限差分方法则采用具有二阶精度的中心差分公式进行。采用具有TVD(Total Variation Diminishing)性质的三阶龙格-库塔多步积分法进行时间积分。对数值模式进行了验证,数值结果同解析解或实验数据吻合良好。  相似文献   

15.
基于均匀网格,建立了沿水深积分的非静压波浪传播数值模型,模型的求解分静压步和动压步两部分。静压步的控制方程为全非线性浅水方程,采用有限体积格式求解,通过采用线性重构技术和全隐式离散底摩阻项,保证了格式的和谐性、守恒性和水深非负性,有效处理了海岸动边界问题。动压步通过应用有限差分方法求解泊松方程考虑动压力,使得模型具备模拟色散性波浪传播的能力。引入波浪破碎指标,波浪破碎后模型退化为静压模型,破碎波自动捕捉为间断。通过算例对所建立模型进行了验证。  相似文献   

16.
本文给出了一种在不同规则四边形网格上求解地下水非稳定流方程的混合拉氏变换局地标有限分析格式。由于将时间t直接包含在有限分析数值计算格式中,克服了传统数值法求解非稳定过程按时段迭代的缺陷,实现了非稳定过程计算的一步完成,提高了计算效率,这对多时段地下水资源管理和长期预报非常有用。  相似文献   

17.
对对流占优的三维溶质运移问题提出了分步广义迎风解法,首先利用N.N.Ya-nenko对水动力弥散方程分步求解的思想,将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题,对对流定解问题采用广义迎风对偶单元均衡法求解,对扩散定解问题采用一般的Galerkin有限元法求解,不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题,而且避免了求节点速度这一步,节省运算步骤,对井点的浓度变化给出了更合适的求解方法。  相似文献   

18.
三维溶质运移问题的分步广义迎风解法   总被引:1,自引:0,他引:1  
对对流占优的三维溶质运移问题提出了分步广义的迎风解法,首先利用N,N,Yanenko对水动力弥散方程分步求解的思想,将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题,对对流定解问题采用广义迎风对偶单元均衡法求解,对扩散定解问题采用一般的Galerkin有限元法求解,不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题,而且避免了求节眯速度  相似文献   

19.
承压含水层中地下水向井非稳定流动的LT有限分析法   总被引:1,自引:0,他引:1  
作者给出了一种求解承压非稳定井流方程的LT有限分析法.其特点是:将时间t直接包含在有限分析数值计算格式中,克服了传统数值法求解非稳定过程按时段迭代的缺陷,实现了非稳定过程计算的一步完成,给多时段地下水资源管理和长期预报带来了方便.用本文提出的方法分别计算了两个算例,结果表明,该方法具有较高的精度,与解析解吻合甚好,可用于实际地下水流的计算之中.  相似文献   

20.
陈亿军  薛强  孙可明  赵颖  万勇 《岩土力学》2012,33(5):1579-1584
建立了土质陡坡降雨侵蚀的数学模型,模型包括陡坡细沟间坡面流控制方程、变沟宽陡坡细沟流控制方程、陡坡细沟间坡面侵蚀泥沙连续方程及陡坡变沟宽细沟侵蚀泥沙连续方程。通过紊流冲击分布概率确定土壤剥蚀方程,确定了泥沙源/汇项的表达形式及细沟高程与宽度的表达式。通过有限元和有限差分相结合的方法求解该数学模型,利用有限单元法对数学模型的水动力方程、运动方程及泥沙连续方程进行空间上的离散,利用有限差分法解决时间域的问题,得出了顺序求解的数值计算公式及模型数值求解的具体方法和步骤。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号