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

2.
利用有限元方法离散椭圆型缓坡方程,能适用于复杂区域,并很好地拟合不规则边界;采用改进共轭梯度法求解离散方程组,可以大大降低计算内存要求,提高计算效率。利用结合上述两种方法的模式对规划的日照港区水域进行了波浪数值计算,并将计算结果与物理模型试验值进行比较,结果表明:该模式能适用于较大区域的波浪场计算,并可以得到较好的计算结果。  相似文献   

3.
径向点插值法在波浪传播数值模拟中的应用   总被引:1,自引:0,他引:1       下载免费PDF全文
针对波浪数值模拟中基于矩形网格的数值方法在深水到浅水的网格间距选择与复杂边界处理上的缺陷,以及基于正交曲线网格和无结构网格的数值方法前处理工作复杂的问题,引入最近在计算力学中发展起来的无网格法——径向点插值法,对经典的双曲型缓坡方程进行空间离散,并在时间上采用四阶Adams-Bashforth-Moulton格式求解建立近岸波浪传播数学模型,通过椭圆形浅滩地形和环形河道的波浪传播计算验证,表明该无网格方法可较为有效地模拟近岸波浪的传播变形,且在处理复杂边界时具有较高的精度.  相似文献   

4.
三棱柱形网格下自由表面流的三维数值模拟   总被引:1,自引:0,他引:1       下载免费PDF全文
为了模拟自由表面流动问题,建立了基于三棱柱形网格的三维Navier-Stokes方程的半隐式动压离散模型。出流边界采用在动量方程右端附加线性衰减项的海绵层处理方式来消减波能。引入k-ε双方程紊流模型求解涡黏性系数,使方程达到封闭。用规则周期波通过淹没障碍物的传播变形和复合明渠流动算例,验证了当遇到短波高频问题时,所述模型模拟结果与测量值吻合良好,体现了模型的精确性和有效性。  相似文献   

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

6.
复杂边界孔口出流的数值模拟   总被引:1,自引:0,他引:1       下载免费PDF全文
运用复势理论建立了求解明流泄洪洞泄流的一种新的计算模型,并提出了未知流量的求解方法。通过数值求解边界积分方程组,对流量、水面线和边壁压力等水力参数进行了数值模拟。计算表明,提出的方法计算速度较快、收敛性较好、精度较高,计算值与试验值符合良好。  相似文献   

7.
为实现颗粒沉降问题的精确数值模拟,探索稠密颗粒两相流中流体与颗粒间的影响机制,基于浸入边界法(IBM),建立了高解析度CFD-DEM-IBM流固耦合数值模拟方法。分别通过流体动力学方法(CFD)及离散单元法(DEM)描述连续流体及非连续颗粒,引入浸入边界法处理颗粒移动边界,并在Navier-Stokes动量方程中附加体力项以体现流固相互作用,采用交错迭代算法在同一时间步内多次迭代求解直至收敛以实现流固间的强耦合。通过单、双、群颗粒的沉降行为模拟,结果表明:与传统CFD-DEM方法相比,该方法能够准确考虑颗粒间及流固两相间的相互作用,获得高解析度的流场信息。模拟结果与前人数值计算结果吻合,验证了该方法的准确性与有效性及其对稠密颗粒两相流问题的适用性与优越性。  相似文献   

8.
非线性波浪变形计算的三维边界元方法   总被引:1,自引:0,他引:1       下载免费PDF全文
采用0-1混合型边界元剖分计算域边界面,以提高边界元方法模拟三维波场波浪变形问题的数值计算精度。采用Spongelayer阻尼消波和Sommerfeld放射条件相匹配的边界处理方式,以消除计算域内的反射波。借助时域内的波面位置追踪,实现了三阶Stokes波在数值波动水槽中的波浪变形计算,其结果吻合于理论值。  相似文献   

9.
墩柱结构是一种常见的工程建筑物,正确求解作用在墩柱上的波浪力具有重要的意义。主要考虑大尺度墩柱上的波浪力的计算,应用线性水波绕射理论来计算作用在大尺度墩柱上的波浪力。应用线性小振幅波理论把水波绕射问题转化为一个二维的Helmholtz方程来求解,再应用Nyst?m方法来求解二维的Helmholtz方程,通过求得的数值解与解析解的对比,表明给出的数值模型有很高的精度。  相似文献   

10.
墩柱结构是一种常见的工程建筑物,正确求解作用在墩柱上的波浪力具有重要的意义.主要考虑大尺度墩柱上的波浪力的计算,应用线性水波绕射理论来计算作用在大尺度墩柱上的波浪力.应用线性小振幅波理论把水波绕射问题转化为一个二维的Helmholtz方程来求解,再应用Nystm方法来求解二维的Helmholtz方程,通过求得的数值解与解析解的对比,表明给出的数值模型有很高的精度.  相似文献   

11.
丁选明  陈育民  孔纲强 《岩土力学》2010,31(Z2):109-114
假定各物理量沿径向不变化,建立了低应变瞬态集中荷载作用下现浇大直径管桩振动响应的计算模型和波动方程。采用Laplace变换法,求得了波动方程的频域解析解,采用Fourier逆变换求得了时域响应。将文中解计算结果与三维频域解析解进行了对比分析,文中解的入射峰-反射峰时间差与三维频域解差别很小,说明对于PCC桩这种大直径薄壁管桩,采用径向不变假定对计算结果几乎没有影响。因此,在PCC桩低应变检测波形分析时,采用基于径向不变假定的二维解是完全合理的。将文中二维频域解的计算结果与二维时域解的结果进行了对比分析,结果表明,2种解在入射波、第一个桩底反射波、第二个桩底反射波峰值大小和到达时间非常吻合,这说明考虑环向位移与否对计算结果没有太大的影响,研究PCC桩低应变动力响应的问题时忽略水平方向的位移是完全可以的。  相似文献   

12.
通常地层介质都是粘滞性的,由于地层的吸收作用,地震波在向下传播过程中能量被衰减,特别是高频能量的衰减,使反射信号频带宽度变窄。为了补偿被吸收的能量,笔者采用粘滞声波方程波场延拓方法进行能量吸收和衰减补偿。首先,从粘滞性波动方程波场延拓函数理论出发,推导了波场延拓函数稳定性条件;其次,对其有效性进行了测试和分析;最后,通过设计数值模型进行正演模拟对比试验,证明了基于粘滞性声波方程的吸收补偿方法能够很好地对粘性介质吸收衰减进行能量补偿,改善处理效果。  相似文献   

13.
针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。  相似文献   

14.
付晓东  盛谦  张勇慧 《岩土力学》2015,36(7):2057-2062
定量研究了非连续变形分析(DDA)方法中的黏性阻尼与数值阻尼。首先,基于Newmark直接积分法,推导了块体系统的运动方程。其次,通过动力学中的黏性阻尼理论建立了DDA中动力系数、时间步长与黏性阻尼比的关系式,探讨了DDA中的常加速度积分方案的数值阻尼分区及阻尼比计算方法,进而得到两种阻尼共同作用时的阻尼比表达式,并分析了频域内阻尼比的分布情况。最后,以谐振激励下的块体振动为例,通过对比不同阻尼作用下块体位移的DDA计算值与理论解,验证了本文提出的阻尼比计算公式的正确性。研究表明:黏性阻尼对低频的衰减作用明显,数值阻尼则可以很快地消除高频干扰,而二者共同作用下可降低阻尼的频率相关性。该研究成果为DDA的振动、波动等动力计算的阻尼取值提供了理论依据。  相似文献   

15.
正演数值模拟是反演成像的基础。为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法。该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常。在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证。模型算例结果表明:该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2 501×2 501的模型模拟时间为4.18 s;适用于任意复杂地形模型。  相似文献   

16.
基于常Q模型的解耦分数阶拉普拉斯算子粘滞波动方程,可以分开模拟振幅衰减和相位错动。但该方程拉普拉斯算子的阶数是随空间变化的,因此数值求解存在一定困难。这里基于截断的泰勒展开,经过一系列近似,推导出拉普拉斯算子的阶数与空间无关的解耦分数阶粘滞弹性波动方程。采用中心差分计算时间导数,使用交错网格伪谱法计算空间导数。数值算例表明,新的方程在处理非均匀介质时具有精度高,计算简便的优点。  相似文献   

17.
林树海  王伟利 《现代地质》2012,26(6):1193-1198
为提高电磁波场数值计算精度,对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一种目前流行的电磁场时域数值计算方法,已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算,时间域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法,该方法使用快速傅里叶度换来计算麦克斯韦方程中的空间导数。由于傅里叶变换的准确性,时间域伪谱法使计算精度提高,数值计算时所需的格子数大大减少,这极大地节省了计算机内存空间,适合于大规模正反演问题的数值计算。  相似文献   

18.
季宪军  梁瑛  欧国强  杨顺  王钧  陆桂红 《岩土力学》2015,36(8):2402-2408
由无黏粗颗粒与黏性泥浆组成的黏性碎屑流,其运动过程会产生不连续变形,基于连续介质假定的流体理论无法描述。根据散体材料理论,在考虑黏性泥浆影响情况下,以PFC3D为平台,编写黏性碎屑坡面运动数值模型试验程序,根据泥浆(成都黏土,密度1.413 g/cm3)室内拉伸试验和旋转剪切试验结果,设置数值模型参数,开展黏性碎屑流坡面运动数值模型试验,再现黏性碎屑坡面运动过程及运动过程中不连续变形现象,并通过同尺寸黏性碎屑坡面运动物理模型试验进行验证。结果表明:基于散体材料理论的PFC3D离散单元法能很好地再现黏性碎屑坡面运动过程及运动过程中不连续变形现象,为深入分析黏性泥浆介质影响下黏性碎屑坡面运动过程提供新的途径。  相似文献   

19.
Viscous boundaries are widely used in numerical simulations of wave propagation problems in rock mechanics and rock engineering. By using such boundaries, reflected waves from artificial boundaries can be eliminated; therefore, an infinite domain can be modeled as a finite domain more effectively and with a much greater accuracy. Little progress has been made, thus far, with the implementation and verification of a viscous boundary in the numerical, discrete element, discontinuous deformation analysis (DDA) method. We present in this paper a new viscous boundary condition for DDA with a higher absorbing efficiency in comparison to previously published solutions. The theoretical derivation of the new viscous boundary condition for DDA is presented in detail, starting from first principles. The accuracy of the new boundary condition is verified using a series of numerical benchmark tests. We show that the new viscous boundary condition works well with both P waves as well as S waves.  相似文献   

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

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