非线性波消波及波场分布研究   总被引：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.  相似文献   

关于Diophantine方程的解     

乐茂华 《广东海洋大学学报》2004,24(4):47-48

对于正整数a ,设S(a)是Smarandache函数。证明了 :方程S(1·2 ) +S(2·3) +… +S(x(x +1) ) =S(x(x +1) (x +2 ) /3)仅有正整数解x =1。  相似文献   

Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

About Darcy's law in non-Galilean frame     

C. Geindreau  E. Sawicki  J.-L. Auriault  P. Royer 《国际地质力学数值与分析法杂志》2004,28(3):229-249

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

用随机介质模型方法描述孔洞型油气储层   总被引：5，自引：0，他引：5  

朱生旺  魏修成  曲寿利  刘春园  吴开龙 《地质学报》2008,82(3):420-427

孔洞型储集层具有强烈的非均质性,在空间分布上变化剧烈,规律性差,其地球物理参数(如速度、密度、弹性参数等)在空间上的变化很难用传统的层状介质模型进行刻画.本文借助随机模型方法,研究孔洞型储层中不同孔洞尺度、分布密度与模型参数的关系,建立反映实际孔洞油气储层介质地球物理参数空间分布统计特征的随机介质模型,以更方便地对孔洞型油气储层进行描述,并为碳酸盐岩孔洞型储层地震反射特征的研究提供有效的方法途径.  相似文献   

Robust adaptive control of underactuated ships on a linear course with comfort     

K. D. Do  J. Pan  Z. P. Jiang   《Ocean Engineering》2003,30(17):2201-2225

This paper addresses an important problem in ship control application—the robust stabilization of underactuated ships on a linear course with comfort. Specifically, we develop a multivariable controller to stabilize ocean surface ships without a sway actuator on a linear course and to reduce roll and pitch simultaneously. The controller adapts to unknown parameters of the ship and constant environmental disturbances induced by wave, ocean current and wind. It is also robust to time-varying environmental disturbances, time-varying change in ship parameters and other motions of the ship such as surge and heave. The roll and pitch can be made arbitrarily small while the heading angle and sway are kept to be in reasonably small bounds. The controller development is based on Lyapunov’s direct method and backstepping technique. A Lipschitz continuous projection algorithm is used to update the estimate of the unknown parameters to avoid the parameters’ drift due to time-varying environmental disturbances. Simulations on a full-scale catamaran illustrate the effectiveness of our proposed controller.  相似文献   

Vertical variations of wave-induced radiation stress tensor   总被引：3，自引：0，他引：3  

Zheng Jinhai  Yan Yixin 《海洋学报(英文版)》2001,20(4):597-605

INTRODUcrIONThe concept of radiation stress was deve1oPed by tonguet--Higgins and Stewart (1964 ),who intreduced the definition of radiation stress as the excess mornentum due to the presence ofwaves, on the basis of time-averaged laws of Newtonian fluid mechanics and the assmption ofa unifOrm velocity distribution over depth. Subequently, the theory has been applied success-fully in the investigation of phenomena such as wave set-up and set--down (Bowen et al.,l968), longshore currents …  相似文献   

声波在水-多孔介质海底界面上的反射与透射     

彭临慧  赵燕鹏  郁高坤 《中国海洋大学学报(自然科学版)》2007,37(4):671-675

声波在海底界面的反射和透射是海底散射、海底混响、海底目标探测的重要问题。利用Biot多孔介质声传播理论对声波在水-多孔海底界面上的反射和透射进行了分析,具体给出反射声波的反射系数,3种透射声波的透射系数以及声能透射系数随入射波入射角和频率(10~40 kHz)的变化关系,分析了各种透射波对透射声能的贡献。多孔海底介质参数分别采用Stoll和Chotiros给出的2组参数进行理论计算。  相似文献   

A linear Boltzmann equation to model wave scattering in the marginal ice zone   总被引：1，自引：0，他引：1  

Michael H. Meylan  Diane Masson   《Ocean Modelling》2006,11(3-4):417-427

We present a linear Boltzmann equation to model wave scattering in the Marginal Ice Zone (the region of ocean which consists of broken ice floes). The equation is derived by two methods, the first based on Meylan et al. [Meylan, M.H., Squire, V.A., Fox, C., 1997. Towards realism in modeling ocean wave behavior in marginal ice zones. J. Geophys. Res. 102 (C10), 22981–22991] and second based on Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136]. This linear Boltzmann equation, we believe, is more suitable than the equation presented in Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136] because of its simpler form, because it is a differential rather than difference equation and because it does not depend on any assumptions about the ice floe geometry. However, the linear Boltzmann equation presented here is equivalent to the equation in Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136] since it is derived from their equation. Furthermore, the linear Boltzmann equation is also derived independently using the argument in Meylan et al. [Meylan, M.H., Squire, V.A., Fox, C., 1997. Towards realism in modeling ocean wave behavior in marginal ice zones. J. Geophys. Res. 102 (C10), 22981–22991]. We also present details of how the scattering kernel in the linear Boltzmann equation is found from the scattering by an individual ice floe and show how the linear Boltzmann equation can be solved straightforwardly in certain cases.  相似文献