首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 546 毫秒
1.
基于势流理论的数值水池可以快速计算波浪的传播及其对建筑物的作用,但是势流理论是基于波浪的无黏性假设的,而在工程中,通常需要在固体边界及波浪破碎的区域考虑黏性效应。针对基于求解Navier-Stokes方程的黏性水池计算量较大、速度较慢的缺点,采用耦合的方法模拟多向不规则波浪的传播,即在外域通过基于势流理论的数值水池产生多向不规则波浪,内域采用求解Navier-Stokes方程和流体体积法(VOF方法)对自由表面进行追踪,通过外域所提供的边界波浪,内域计算可以在较小区域进行计算,从而达到减少计算量、提高计算效率的目的。  相似文献   

2.
波生沿岸流数值模拟   总被引:4,自引:0,他引:4       下载免费PDF全文
为了更好地研究近岸海域波生沿岸流,建立了基于高阶Boussinesq水波方程的波生沿岸流时域数值模型。控制方程在中等水深范围内具有较好的色散性和变浅作用性能,同时具有二阶完全非线性特征,适合描述近岸区域波浪强非线性运动。通过采用松弛造波方法实现了非线性波浪的无反射入射,采用周期性侧边界条件模拟开敞边界。通过数值试验,讨论了模型中主要参数对数值结果的影响。利用率定后的参数模拟了均匀坡度海岸上产生的沿岸流,通过和实验数据的对比验证了模型的准确性和适用性。利用模型数值模拟了不同波浪入射条件(包括周期、波高和波浪入射角度)对波生沿岸流的影响。  相似文献   

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

4.
应用二阶完全非线性Boussinesq方程模拟破碎波浪   总被引:1,自引:0,他引:1       下载免费PDF全文
建立了基于高阶Boussinesq水波方程的一维波浪破碎数值模型。基于一组具有二阶完全非线性特征的Boussinesq水波方程,建立了交错网格下的高精度差分格式,推导了适用于该组方程的永形波解析解,其和松弛造波技术相结合实现了数值波浪水槽中(强)非线性波浪的无反射入射。通过模拟封闭容器内水体晃动问题对数值格式进行了验证,通过模拟孤立波在斜坡海岸上的浅化过程说明了将方程从弱非线性发展到完全非线性的必要性。采用涡粘方法处理波浪破碎,利用物理模型实验数据,分析了模型中各波浪破碎参数对数值结果的影响并对参数进行了率定。应用该模型对规则波在斜坡海岸上的传播、变浅以及破碎过程进行了数值模拟研究,数值结果同实验数据吻合良好,验证了模型的有效性。  相似文献   

5.
陈胖胖  王军 《冰川冻土》2012,34(6):1353-1357
数值模拟是河冰问题研究的重要手段之一, 但受限于天然河道的复杂地形条件, 研究工作开展中往往面临不规则边界的处理问题.基于河冰数值模拟工作现状, 总结和阐述了贴体网格技术对其的应用.针对二维数值模拟提出了采用滑动边界法和势流理论式相结合的方法来进一步提高边界处网格的正交性; 三维网格生成则宜采用σ竖向坐标变换法; 河冰模拟中水面位置确定通常采用标高函数法; 对于非稳态河冰数值模拟, 动边界处理技术可考虑改进并应用干湿法.  相似文献   

6.
波流共存场中多向随机波浪传播变形数学模型   总被引:1,自引:0,他引:1       下载免费PDF全文
基于波作用量守恒方程建立了波流共存场中多向随机波浪传播变形数学模型,模型中考虑了波浪绕射的影响和水流引起的波浪弥散多普勒效应,应用包含水流和地形影响的激破波模式计算波浪破碎的能量耗散,采用一阶上迎风有限差分格式离散控制方程。分别计算了有无近岸流情况下单向和多向随机波浪的波高分布,考虑水流影响的数值计算结果与物理模型实验数据吻合良好,比较分析表明,所建立的数学模型能够复演由于离岸流引起的波高增大,可用于波流共存场多向随机波浪传播变形的模拟和预报。  相似文献   

7.
复杂边界下的流场数值模拟常基于非结构化网格进行求解,建立一种在非结构混合网格上求解水深平均的二维浅水方程的模型,以便精确模拟复杂边界、提高计算效率。该模型时间项离散采用隐格式使得模型具有较好的稳定性,对流项和扩散项分别采用总变差减小(Total Variation Diminishing,TVD)格式和构造辅助点的方法来离散,同时采用水深平均的标准k-ε模型来封闭湍流模型。选用两个经典验证算例检验模型,计算结果表明,基于非结构混合网格开发的模型具有较高的精度,且收敛性能较好。  相似文献   

8.
利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。  相似文献   

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

10.
基于椭圆型缓坡方程和消息传递接口机制,提出了缓坡地形上非线性波浪变形的一种并行数值模拟方法,并在微机群集并行系统上进行了数值实现.对一种典型缓坡地形上的非线性波浪传播变形进行了数值模拟,数值结果与实验数据的比较表明,在椭圆型缓坡方程中应用非线性波色散关系可有效考虑波浪的非线性效应,同时从获得的并行效率和并行加速比来看,并行数值方法可有效提高数值模拟效率.  相似文献   

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

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

13.
岩体材料物理网格对流形元覆盖系统形成的影响   总被引:5,自引:1,他引:4  
朱爱军  邓安福  颜昌武  邓卫东 《岩土力学》2004,25(12):1933-1936
流形元的数学覆盖与物理网格相对独立,使流形元覆盖系统的形成比有限元单元网格的形成有更大的灵活性和方便性,但在形成流形元覆盖系统时对物理网格也有一定的依赖。笔者讨论了使用规则网格作为数学覆盖时,岩体材料普遍存在的裂隙和岩层界面这两种物理网格对流形元覆盖系统形成的影响,指出数学覆盖的构造要与岩体裂隙和岩层界面相适应,并提出了相应的流形元覆盖系统的形成方法,有利于流形元法应用于大型岩体工程问题。  相似文献   

14.
数据百分比网格法适用于区域化探资料数据的处理和研究,由于不受图幅比例尺、采样方法和分析精度的限制,也不影响图幅的拼接,因此是一种简单、易行的数学地质方法。通过对宁芜及溧水地区化探资料的进一步处理,取得了较为满意的效果。  相似文献   

15.
16.
宋梅  张石峰 《新疆地质》2004,22(3):311-314
通常情况下,含水层的参数可分别借助于Dupuit和Theis井流公式,在局部地段利用稳定流和非稳定流抽水试验获得.在新井设计阶段,一般是依据现有参照井数据进行线性或非线性插值计算,但插值受插值表的约束,效率不高.本次研究采用2次优化方法,用网格法确定所求参数最优估计值作为下一步的叠代初值.用共轭梯度法逐步逼近所求参数的精确解,高效地利用参照井的基本数据得到了精度较高的新井参数,为准确地设计新井提供可靠的依据,实地应用效果良好.  相似文献   

17.
梅发国  高彦伟  郭华 《世界地质》2002,21(4):385-389
针对差分法求解过程中细网格计算所需时间和空间要求高以及在粗网格上计算又达不到所要求的精度等问题,利用粗网格上的残差校正特性消除迭代误差的低频分量,同时利用细网格上的松弛光滑特性消除迭代误差的高频分量,采用不同疏密的网格距消除不同频率范围内的误差分量,将多重网格方法成功应用于频域波动方程的求解,并给出了实际算例。算例中,分别用Gauss-Seidel迭代法和多重网格法求其数值解,并绘出效果图;比较得出多重网格方法具有精度高、收敛速度快和易于实现的特性,更适合于高维波动问题的求解。  相似文献   

18.
新疆天山某煤炭勘查区,需查明一水平内火烧煤层的下边界及其富水性。在综合分析该区地质资料及物性参数的基础上,采用了以磁法和瞬变电磁法为主的物探方法。磁法勘探测网设计为20m×80m,并在后期对西部测线进行了加长。根据测区实测磁异常分布图、低通滤波后磁异常分布图、向上延拓30m磁异常分布图及垂向二阶导数磁异常分布图,圈定火烧区主要分布在测区的北部。瞬变电磁法勘探测网设计为40m×80m,依据下1、下5及下10煤层底板富水异常分布图,解释下1、下5煤底板的富水性主要集中在测区的北部,与火烧区范围高度一致,而下10煤层底板富水性极不均匀,仅在测区北部煤层露头附近,受砂岩裂隙的影响,接受降水补给形成零星砂岩裂隙富水区。该解释成果经矿方钻孔验证,结论基本可靠。  相似文献   

19.
为了安全回采受底板承压水害威胁的采煤工作面,可采取注浆改造的防治水方法,而注浆改造设计是注浆改造的基础和前提。由于传统设计方法较为繁琐且效率较低,本次研究,首先对方格网布孔法的设计原理与基本思路进行了介绍,然后以新桥煤矿北二采区2207工作面为例,对根据工作面大小和底板水文地质情况设计方格网,输入基本参数,得出需要施工的钻孔参数,并根据物探资料及实际揭露的底板灰岩岩溶裂隙发育情况的变化,及时调整钻孔数量与位置进行了详细说明。应用效果表明,该方法具有设计简便,钻孔布置均匀,重点突出,清除了设计盲区,底板注浆改造效果好的特点,是一种操作简便、安全高效的注浆设计方法。  相似文献   

20.
交错网格下的浅水方程高分辨Gauss型格式   总被引:1,自引:0,他引:1       下载免费PDF全文
在交错网格上,基于高精度的Gauss积分公式,针对浅水波方程设计了对模拟涌波具有高分辨率的完全二阶精度的数值计算格式。由于采用了交错网格,差分格式不需要解Riemann问题,因此本文格式具有计算简单、工作量少、编程简便等特点。另外,在一维单个方程时,本文格式在CFL(Courant Friedrich Lewy)条件限制下为TVD(Total Variation Diminishing)格式,在二维和三维情况下格式具有MmB(Maximum and Minimum Bounds Preserving)性质。利用国家高性能计算中心(合肥)的曙光1000型分布存储大规模并行机,对在交错网格下所构造的求解浅水方程的高分辨差分格式进行了并行实现,几个算例的计算结果令人满意。  相似文献   

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

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