首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 140 毫秒
1.
通量向量分裂格式的二维水流-水质模拟   总被引:16,自引:2,他引:16       下载免费PDF全文
赵棣华  姚琪  蒋艳  杨珏  逄勇 《水科学进展》2002,13(6):701-706
在有限体积法框架下应用通量向量分裂(FVS)格式进行平面二维水流 水质模拟。通过对偏微分方程进行有限体积的积分离散、利用通量的旋转不变性,把二维问题转化为一系列局部的一维问题进行求解,采用FVS格式计算各跨单元边界的水量、动量及污染物输运等通量。应用该格式计算了理想条件下的浓度输移,其结果与精确解拟合很好;模拟了长江江苏靖江段的水质及污染带,计算结果与水质监测值相当吻合,为长江江苏干流段水质评价提供了科学依据。  相似文献   

2.
基于黎曼近似解的溃堤洪水一维-二维耦合数学模型   总被引:2,自引:0,他引:2       下载免费PDF全文
采用一维-二维耦合数学模型对溃堤洪水进行模拟,可以发挥两者优势,提高计算效率。将基于黎曼近似解Godunov格式的一维、二维水流数学模型通过堰流公式进行耦合,实现交界面的水力要素交互;一维、二维模型均采用中心形式的有限体积法求解,分别采用HLL格式和Roe格式计算界面通量,具有一致的空间同步性。通过经典算例验证了模型的可靠性和稳定性,应用耦合模型对1998年松花江干流胖头泡溃堤洪水进行了模拟计算,较好地模拟了溃堤洪水的演进过程与淹没范围,对200年一遇洪水的分洪效果进行了分析,为此区域的分洪利用及松花江中上游防洪调度提供参考依据。  相似文献   

3.
复杂边界及实际地形上溃坝洪水流动过程模拟   总被引:7,自引:3,他引:4       下载免费PDF全文
建立了基于无结构三角网格下采用有限体积法求解的二维水动力学模型,用于模拟溃坝洪水在复杂边界及实际地形上的流动过程。该模型采用Roe格式的近似Riemann解计算界面水流通量,结合空间方向的TVD-MUSCL格式及时间方向的预测-校正格式,可使模型在时空方向具有二阶计算精度。模型中引入最小水深概念,提出了有效的干湿界面处理方法。模拟了理想条件下溃坝水流过程,研究不同最小水深取值对干河床上洪水演进的影响,并用两组简单溃坝水流的水槽试验资料对模型进行验证。采用该模型模拟了实际溃坝洪水的流动过程,所得计算结果与实测资料及已有模型计算结果较为符合。  相似文献   

4.
长江江苏段区域供水水源地水质可达性研究   总被引:7,自引:0,他引:7       下载免费PDF全文
采用平面水流、水质有限体积法及黎曼近似解模型,在无结构网格上对偏微分方程进行有限体积法的积分分离散,模型具备有限单元法及有限体积法两者的优点。采用同步遥感、水量、水质实测资料对模型进行了率定验证,在对长江水源地、沿江污染源调查以及岸线稳定段分析结果基础上,进行污染源概化以及确定主要研究江段。在对概化的污染源计算出在设计条件下的污染带分布图,并计算出污染带长、宽随排污量变化曲线;在上述研究基础上,分析出长江各水源地及岸线稳定段的水质变化状况。对长江江苏段各主要排污口的污染带进行了计算,并对长江江苏段区域供水水源地的水质可达性进行了研究。  相似文献   

5.
针对溃坝水流数值模拟面临的复杂地形和不规则边界等问题,基于结构网格建立了适应复杂地形和不规则边界的溃坝水流数值模拟有限体积模型(HydroM2D)。模型基于具有守恒特性的二维浅水方程,利用HLLC格式的近似Riemann解计算网格界面通量,利用MUSCL-Hancock法不断向前积分,使模型在时空上具有二阶精度;对源项进行离散处理确保模型的稳定性;模型引入有效干湿边界和不规则地形边界处理方法,准确模拟了干湿单元的动态交替和复杂边界上的水流特性。最后分别利用水槽试验、物理模型和实际算例对模型进行验证。结果表明,该模型对不同情景下的溃坝洪水模拟结果和实测资料以及现有模型模拟结果具有较高的一致性,模拟精度较高,稳定性较好,具有推广应用价值。  相似文献   

6.
长江中下游河湖洪水演进的数值模拟   总被引:13,自引:1,他引:13       下载免费PDF全文
以长江中下游(宜昌-大通)防洪系统为对象,在水文、河道及湖盆地形现状条件下,建立了一个能适应各种复杂条件的一二维非恒定流模型,来进行长江干流、河网、湖泊、分蓄洪区垸及水库的洪水演进和调度仿真.所建模型的洞庭湖部分采用无结构网格二维非恒定有限体积格式,以适应湖区复杂的边界形状和保持水量平衡.河网区部分采用一维非恒定流显隐结合的分块三级算法,以准确实现河网汊点流量的自动分配和往复流动.为了提高模拟精度和扩展模拟功能,在水流数值模拟的范围内侧重讨论了内外动边界处理、分蓄洪运用及阻力项计算等环节,提出了合理可行的数值处理方法.采用20世纪80年代至90年代共6年汛期洪水资料对所建模型进行了严格的率定和检验,高精度的模拟结果证实了模型的合理性和有效性.  相似文献   

7.
计算域内地形的概化精度对模拟结果有较大影响,针对实际地形中同时存在狭长河谷与广阔泛洪区的问题,在层次自适应网格模型的基础上,研究了网格快速加密与合并的步骤,并以水位梯度与局部弗劳德数为基础设计了网格自适应准则,当水流运动情势变化时网格密度自动调整,实现计算精度与效率的平衡。在此网格模型基础上,采用有限体积法求解二维浅水方程,利用梯度限制器技术及龙格-库塔法提高模型的空间、时间计算精度。算例表明,层次自适应网格模型既能实现随水流运动动态变化并捕获水位计算敏感区,也能对局部区域进行静态固定加密,自适应性良好,具有较好的推广应用价值。  相似文献   

8.
应用多尺度有限单元法模拟非均质多孔介质中的三维地下水流问题。与传统有限单元法相比,多尺度有限单元法的基函数具有能反映单元内参数变化的优点,所以这种方法能在大尺度上抓住解的小尺度特征获得较精确的解。在介绍多尺度有限单元法求解非均质多孔介质中三维地下水流问题的基本原理之后,对参数水平方向渐变垂直方向突变的非均质多孔介质中的三维地下水流和Borden实验场的三维地下水流分别用多尺度有限单元法和传统等参有限单元法进行了计算,结果表明在模拟高度非均质多孔介质中的三维地下水流问题时,多尺度有限单元法比传统有限单元法有效,既节省计算量又有较高的精度;在模拟非均质性弱的多孔介质中的三维地下水流问题时,多尺度有限单元法虽然也能在大尺度上获得较为精确的解,但效果不明显。  相似文献   

9.
城市地区存在密集的建筑,采用传统二维数值模拟方法模拟城市洪水,需要提取城市复杂的建筑物作为计算边界和加密计算网格,难以应用于大尺度城市洪水模拟。为此,引入容积率系数修改原始二维浅水方程,以反映建筑物对城市洪水演进的影响,并采用人工逆风通量向量分裂有限体积法构建数值计算模型,求解修改后的二维浅水方程。为了检验模型,采用两个算例进行验证,并对计算结果对比分析。研究表明,该模型不需要进行建筑边界提取和精细计算网格,精度上,完全能够满足大尺度城市洪水模拟精度要求。  相似文献   

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

11.
水体交换年龄模型研究   总被引:2,自引:0,他引:2       下载免费PDF全文
采用剖开算子法,把水体年龄控制方程分成几个连续的初值问题。在任意三角形网格中,分别对不同性质的算子采用各自适合的算法,即采用特征线法求解对流分步,采用半隐式有限元法求解扩散分步和传播分步。利用水体交换模型试验结果和水体年龄对称特性解析解对所建立的年龄模型进行了验证计算,结果表明,建立的水体年龄模型计算结果与试验结果、解析解结果吻合,可较好地预测水体交换年龄。  相似文献   

12.
Appropriate thermal design of an artificial soil freezing system should include the linking between the performance of the refrigeration system and the thermal regime in the soil, as well as the influence of variable thermal properties and complicated geometries.

The paper presents design systems, based on computer programs, where these features are demonstrated. The refrigeration capacity is given as a boundary condition at the freezing pipe, either as temperature or flux.

Two programs are discussed; a one-dimensional (1-D), finite difference program for the analysis of a single pipe, and a two-dimensional (2-D), finite element code, with seepage flow as an option, more suited for real problems.

The programs are used to display the influence of various design parameters, such as refrigeration capacities, material properties and geometries.

A comparison is made between the 2-D program and a simplified method developed by K.R. Khakimov, which shows that the simplified method gives a very conservative estimate of the frozen volume vs. elapsed time.

Finally, the 2-D program computing combined heat and seepage flow is applied to a laboratory model of a soil freezing system, and the agreement between measured and computed values of temperature and water flow is shown to be acceptable.  相似文献   


13.
河网海湾水动力联网数学模型   总被引:15,自引:3,他引:12       下载免费PDF全文
在一维显式、二维隐式水动力联网数学模型的基础上,建立了一种一维、二维全隐河网海湾水动力联网数学模型。该一维河网模型采用Preissmann四点隐式格式,用节点水位控制法进行数值计算,二维海湾模型采用改进型双向隐式(DSI)法进行数值求解。在河口一维、二维连接处,水力因子通过接口断面法传递,无需重叠一个一维河段来传递水力因子,避免了口门处二维网格需取较小尺度的问题。模型在珠江三角洲河网及横门、洪奇沥口门海域做了检验,验证结果表明,建立的河网海湾水动力联网数学模型是可行的,可以用于河网、河口治理工程的数值研究。  相似文献   

14.
Saltwater intrusion into coastal freshwater aquifers is an ongoing problem that will continue to impact coastal freshwater resources as coastal populations increase. To effectively model saltwater intrusion, the impacts of increased salt content on fluid density must be accounted for to properly model saltwater/freshwater transition zones and sharp interfaces. We present a model for variable density fluid flow and solute transport where a conforming finite element method discretization with a locally conservative velocity post-processing method is used for the flow model and the transport equation is discretized using a variational multiscale stabilized conforming finite element method. This formulation provides a consistent velocity and performs well even in advection-dominated problems that can occur in saltwater intrusion modeling. The physical model is presented as well as the formulation of the numerical model and solution methods. The model is tested against several 2-D and 3-D numerical and experimental benchmark problems, and the results are presented to verify the code.  相似文献   

15.
范留明 《岩土力学》2013,34(9):2715-2720
平面波的传播问题通常可以归结为一维波动方程的定解问题。在非均匀介质中,即使简单的一维波动方程也需要借助于数值方法获得近似解。3层5点古典差分格式是计算偏微分方程一种常用算法,作为一种显式迭代格式,需要满足稳定性条件 ,其中 为波速, 为空间采样间隔, 为时间采样间隔。当 时, ,古典差分格式达到临界稳定状态。在这种情况下,平面波在 时间内的传播距离恰好等于空间采样间隔,差分格式真实地反映了平面波的传播原理,因而可以得到一维波动方程的精确解。但是,由于在非均匀介质中存在不连续的波阻抗界面,此方法不适于计算非均匀介质的波场。为了将临界稳定情况下的古典差分格式推广应用至非均匀层状介质,提出了一种能够处理波阻抗界面的有限差分格式,并应用傅里叶分析法得到其稳定性条件。模型算例验证了此算法的正确性。  相似文献   

16.
深基坑降水中不同的地质模型有限元应用研究   总被引:2,自引:0,他引:2  
王翠英  王家阳 《岩土力学》2007,28(6):1259-1264
分析水-土耦合场相互作用时,将岩土视作弹塑性材料。在非线性本构关系下,考虑渗透作用,得出岩土体有限元方程,建立了基坑降水的三维水-土耦合的数学模型。对基坑降水中不同的地质模型潜水、承压水及有越流补给和实际工程条件井管、过滤管、止水帷幕分别进行了计算处理,并以武汉国际证券大厦降水工程为实例,对其降水过程进行了三维有限元计算。应用三维数值模拟方法,能很好地模拟基坑防渗帷幕、非均质和各向异性、承压-无压等一些解析法难以处理的实际工程条件和地质模型。  相似文献   

17.
An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

18.
给出电阻率、激发极化法数据的二维有限元反演方法和原理 ,并将此方法用于高龙金矿电法勘探资料解释中 ,解释结果与已知地质资料基本吻合。由此证明 ,利用双线性模型进行二维反演的有效性及其独特的优点。  相似文献   

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

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