非线性波消波及波场分布研究   总被引：1，自引：1，他引：0  

王巧莎  李明海 《海洋科学》2018,42(6):30-38

本文基于雷诺平均N-S方程,并结合RNG k-e方程建立了粘性数值波浪水槽,对不同波陡、不同相对水深、不同相对波高的非线性规则波的阻尼消波问题和波场分布进行研究。文中提出了两种描述消波区内部阻尼变化的阻尼函数,分别适用于小波陡情形和高波陡情形。研究结果表明,小波陡组消波区可设为一个波长,阻尼系数取10~4～10~5即可满足消波要求,计算结果与实验结果及造波理论吻合良好;高波陡组消波区可设为两个波长,阻尼系数取10~4～10~5亦可满足消波要求,计算结果与实验结果吻合良好。此外,当波陡较小时,波场内反射情况的小幅改变即可对整个波场造成影响,特别是当水深较浅时这种影响极为明显,需谨慎考虑。当波陡较大时,水波能量较高,整个波场沿水波传播方向可观测到明显的衰减现象,在具体试验中需进行考虑。  相似文献   

Three dimensional simulation of fluid flow in X‐ray CT images of porous media     

A. Al‐Omari  E. Masad 《国际地质力学数值与分析法杂志》2004,28(13):1327-1360

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction     

M. Zeghal  U. El Shamy 《国际地质力学数值与分析法杂志》2004,28(14):1361-1383

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

P‐wave transmission across fractures with nonlinear deformational behaviour     

X. B. Zhao  J. Zhao  J. G. Cai 《国际地质力学数值与分析法杂志》2006,30(11):1097-1112

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity     

Jean-Francois Guéguen  Pho Hoang-Trong  Hilaire Legros 《Geophysical Journal International》1987,89(3):839-855

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.  相似文献   

PSO-RBFNN模型及其在岩土工程非线性时间序列预测中的应用   总被引：1，自引：0，他引：1  

梁桂兰  徐卫亚  何育智  赵延喜 《岩土力学》2008,29(4):995-1000

岩土工程受力变形演化是一个典型的非线性问题,其演化的高度非线性和复杂性,很难用简单的力学、数学模型描述,但可用粒子群优化径向基神经网络对岩土工程应力、位移非线性时间序列进行动态实时预测。网络径向基层的单元数通过均值聚类法确定后,所有其它参数：中心位置、形状参数、网络权值,均通过粒子群优化算法在全局空间优化确定。工程实例应用表明,该模型预测结果准确、精度高,有良好的应用前景。  相似文献   

影响低渗透油藏低速非线性渗流的实验研究     

王正波  岳湘安  韩冬  王金峰 《矿物学报》2008,28(1):48-54

以膨胀性的粘土为主要研究对象,通过粘土浆体的流变性实验和填砂管渗流实验,并使用一种新研制的微流量测量工具,研究了低渗岩心在较低流速下的渗流特性;利用微尺度流动理论,Einstein粘度定律等进行了理论分析,对引起低渗透油藏低速非线性渗流的影响因素进行了初步的探索.结论认为,粘土的膨胀性以及其水化后形成浆体的流变性是引起低速非线性渗流的主要因素之一.  相似文献   

海洋石油管道单点提升分析   总被引：3，自引：0，他引：3  

邢静忠  柳春图等 《海洋工程》2002,20(3):29-33

基于管道微单元体平衡建立了海管单点提升的非线性力学模型的控制微分方程组 ,使用变弧长的无量纲代换将动边界问题化为固定边界的两点边值问题 ,利用Maple环境下编制的两点边值问题的打靶法程序得到了该问题在各个提升阶段的数值解答和在单点提升过程中管道的极限弯矩约为 0 .71q1 3(EI) 2 3。  相似文献   

Internal generation of waves: Delta source function method and source term addition method     

Gunwoo Kim  Changhoon Lee  Kyung-Duck Suh   《Ocean Engineering》2007,34(17-18):2251-2264

In this study, we investigate two internal wave generation methods in numerical modeling of time-dependent equations for water wave propagation, i.e., delta source function method and source term addition method, the latter of which has been called the line source method in literatures. We derive delta source functions for the Boussinesq-type equations and extended mild-slope equations. By applying the fractional step splitting method, we show that the delta source function method is equivalent to the source term addition method employing the energy velocity. This suggests that the energy velocity should be used rather than the phase velocity for the transport of incident wave energy in the source term addition method. Finally, the performance of the delta source function method is verified by accurately generating nonlinear cnoidal waves as well as linear waves for horizontally one-dimensional cases.  相似文献   

A New Pursuit Mode of Nonlinear Wave Surface     

孙大鹏李玉成  郭海滨 《中国海洋工程》2004,18(2):257-266

The numerical mode of nonlinear wave transformation based on both the Laplace equation for water field and the Bemoulli equation for water surface is a kind of time-domain boundary problem with initial conditions. And the basis for establishing the numerical mode of nonlinear wave in time domain is to trace the position of wave free surface and to calculale the instantaneous surface height and surface potential function. This paper firstly utilizes the ‘0-1‘ combined BEM to separate the boundary by means of discretization of Green‘ s integral equation based on the Laplace equation, then separates the free surface of wave with FEM and derives the FEM equation of wave surface that satisfies the nonlinear boundary conditions. By jointly solving the above BEM and FEM equations, the wave potential and surface height could be obtained with iteration in time domain. Thus a new kind of nonlinear numerical mode is established for calculating wave transformation. The wave test in the numerical wave tank shows that the numerical simulation with this mode is of high accuracy.  相似文献