共查询到20条相似文献,搜索用时 31 毫秒
1.
利用算子分裂迎风均衡格式解对流为主溶质运移问题 总被引:2,自引:0,他引:2
水污染模拟问题是水流问题与溶质运移问题的耦合问题.各种常见的数值解法在以对流为主溶质运移问题的求解中都会遇到困难,如用有限单元法或有限差分法时,会产生数值弥散与过量这两类误差.引入算子分裂迎风均衡格式法求解对流为主的水污染模拟问题,较好地克服了数值弥散和数值解出现振荡问题,该格式具有良好的稳定性、单调性及守恒性特点. 相似文献
2.
3.
基于有限体积方法和结构化网格,建立了海底滑坡引起的波浪传播数值模型。模型控制方程为考虑了海床随时间变化的二维浅水方程。采用中心迎风格式计算控制体界面数值通量,采用线性重构技术、局部海床高程处理技术和全隐式离散底摩阻项,保证了格式的和谐性、守恒性和水深非负性,有效处理了海岸动边界问题。时间积分采用具有强稳定性质的二阶龙格-库塔方法(ssp RK)。针对滑坡体兴波经典算例开展数值模拟,将计算结果与解析解、实验结果及其他模型计算结果进行比较和分析。结果表明,对于所考虑的计算工况,模型能较合理地模拟滑坡兴波的产生、传播和爬高过程。 相似文献
4.
针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。 相似文献
5.
根据水质模型的具体特点,对不同的方程采用不同方法,水流问题用有限元法;对流弥散方程先用算子分裂的方法分解为两个方程,即对流方程和弥散方程,前者用高精度广义迎风格式求解,对弥散方程则采用多单元均衡格式法求解,最后合成为高精度广义迎风均衡格式求出溶质浓度。通过对数值实验例子的计算和实验溶质迁移的模拟,可以看出在求解对流弥散定解问题时,广义迎风均衡格式克服了有限元数值波动和浓度出现负值的问题,与有限元相比有较大改进。 相似文献
6.
7.
8.
9.
提高畦灌施肥地表水流与溶质运移数值模型的稳定性、收敛性及计算精度,有利于改善地面畦灌施肥系统的设计与评价工作.在利用隐-显混合时间格式对一维畦灌地表水流与溶质运移耦合模型中包含的各矢量项进行时间离散基础上,借助有限差分法、有限体积法和有限单元法分别对由隐时间格式生成的物理矢量线性近似式空间导数、物理矢量空间导数、溶质扩散矢量和地形矢量项进行空间离散,对最终形成的控制方程代数方程组进行数值求解,构建起基于混合数值解法的一维畦灌施肥地表水流与溶质运移耦合模型. 相似文献
10.
11.
泥石流是山区多发的一种地质灾害,它的发生和发展威胁着人们的生命和财产安全,影响着人们的正常生活,因而需要加强对其发生和发展过程的研究。结合泥石流的动力模型方程采用数值模拟方法再现泥石流发生和发展的过程,是研究和预测模拟泥石流灾害的有效手段。目前的动力模型方程大多只关注动力过程,却忽视了静动力过程的统一,这将导致在一些情况下产生错误的结果。本文研究了一维泥石流的静动力阻力特征,通过修正泥石流动力学方程的阻力项,得到了具有静动力统一特征的模型方程。并以Roe格式的近似Riemann解为基础,采用MUSCL线性重构方法建立了具有较高精度和分辨率的有限体积数值求解。具体算例的数值验证表明,方程阻力项的修正是合理的,所建立的数值求解也是稳定和有效的。 相似文献
12.
格点法是在计算流体力学中首先发展起来的数值模拟新方法,共根本思想是对问题重新建模,建立直接模拟流体运动的离散格点模型。本文以一阶拟线性双典型方程及Kdv方程为例,推广这种求解一般偏微分方程的。它是运用多尺度分析方法,构造出格点模型的演化方程的局部平衡分布函数。数值试验表明,该方法程序实现简单,求解速度快,数值结果令人满意。 相似文献
13.
以流函数方程作为控制方程,与N-S方程相比,不出现压强项和速度项,与涡量-流函数型N-S方程相比,没有涡量项和速度项,只有流函数项,避免了涡量边界条件问题,将控制方程中未知数的个数变为一个,通过解得流函数而得到其他各项。考虑对流项迎风性质,采用迎风差分格式,对突扩层流进行数值模拟,与DurstF流动试验及混合有限分析法模拟结果比较,基本吻合。通过一系列计算,得出主流经突扩断面偏转后相关特征值与雷诺数Re的关系,扩展了突扩流动的研究内容。 相似文献
14.
为高效求解扩展型Boussinesq水波方程,建立了基于有限差分和有限体积方法的混合数值格式。将一维控制方程写为守恒形式,方程中通量部分采用有限体积方法求解,剩余部分采用有限差分方法求解。其中,有限体积方法采用Godunov类高分辨率格式,并结合HLL(Harten-Lax and van Leer)式黎曼问题近似解求界面数值通量,黎曼问题界面左右变量通过高精度状态插值方法(MUSCL)构筑。有限差分方法则采用具有二阶精度的中心差分公式进行。采用具有TVD(Total Variation Diminishing)性质的三阶龙格-库塔多步积分法进行时间积分。对数值模式进行了验证,数值结果同解析解或实验数据吻合良好。 相似文献
15.
16.
王文科 《西安地质学院学报》1996,18(3):56-62
本文给出了一种在不同规则四边形网格上求解地下水非稳定流方程的混合拉氏变换局地标有限分析格式。由于将时间t直接包含在有限分析数值计算格式中,克服了传统数值法求解非稳定过程按时段迭代的缺陷,实现了非稳定过程计算的一步完成,提高了计算效率,这对多时段地下水资源管理和长期预报非常有用。 相似文献
17.
姚磊华 《吉林大学学报(地球科学版)》1997,(4)
对对流占优的三维溶质运移问题提出了分步广义迎风解法,首先利用N.N.Ya-nenko对水动力弥散方程分步求解的思想,将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题,对对流定解问题采用广义迎风对偶单元均衡法求解,对扩散定解问题采用一般的Galerkin有限元法求解,不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题,而且避免了求节点速度这一步,节省运算步骤,对井点的浓度变化给出了更合适的求解方法。 相似文献
18.
三维溶质运移问题的分步广义迎风解法 总被引:1,自引:0,他引:1
姚磊华 《长春地质学院学报》1997,27(4):424-428
对对流占优的三维溶质运移问题提出了分步广义的迎风解法,首先利用N,N,Yanenko对水动力弥散方程分步求解的思想,将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题,对对流定解问题采用广义迎风对偶单元均衡法求解,对扩散定解问题采用一般的Galerkin有限元法求解,不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题,而且避免了求节眯速度 相似文献
19.
承压含水层中地下水向井非稳定流动的LT有限分析法 总被引:1,自引:0,他引:1
作者给出了一种求解承压非稳定井流方程的LT有限分析法.其特点是:将时间t直接包含在有限分析数值计算格式中,克服了传统数值法求解非稳定过程按时段迭代的缺陷,实现了非稳定过程计算的一步完成,给多时段地下水资源管理和长期预报带来了方便.用本文提出的方法分别计算了两个算例,结果表明,该方法具有较高的精度,与解析解吻合甚好,可用于实际地下水流的计算之中. 相似文献
20.
建立了土质陡坡降雨侵蚀的数学模型,模型包括陡坡细沟间坡面流控制方程、变沟宽陡坡细沟流控制方程、陡坡细沟间坡面侵蚀泥沙连续方程及陡坡变沟宽细沟侵蚀泥沙连续方程。通过紊流冲击分布概率确定土壤剥蚀方程,确定了泥沙源/汇项的表达形式及细沟高程与宽度的表达式。通过有限元和有限差分相结合的方法求解该数学模型,利用有限单元法对数学模型的水动力方程、运动方程及泥沙连续方程进行空间上的离散,利用有限差分法解决时间域的问题,得出了顺序求解的数值计算公式及模型数值求解的具体方法和步骤。 相似文献