首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
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.  相似文献   

2.
A 2-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear wave forces on a fixed box-shaped ship in a harbor.The domain is divided into an inner domain and an outer domain.The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations.The other area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions.Along the interface boundaries between the inner domain and the outer domain,the volume flux is assumed to be continuous and the wave pressures are equal.Relevant physical experiment is conducted to validate the present model.It is shown that the numerical results agree with the experimental data.Compared with the coupled model with the flow in the inner domain governed by the Laplace equation,the present coupled model is more efficient and its solution procedure is more simple,which is particularly useful for the study on the effect of the nonlinear wave forces on a fixed box-shaped ship in a large harbor.  相似文献   

3.
In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton's method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C~0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.  相似文献   

4.
A 3-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear waves acting on a box-shaped ship fixed in a harbor. The domain is divided into the inner domain and the outer domain. The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations. The remaining area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions. Along the interface boundaries between the inner domain and the outer domain, the volume flux is assumed to be continuous and the wave pressures are equal. Relevant physical experiment is conducted to validate the present model and it is shown that the numerical results agree with the experimental data. Compared the coupled model with the flow in the inner domain governed by the Laplace equation, the present coupled model is more efficient and its solution procedure is simpler, which is particularly useful for the study on the effect of the nonlinear waves acting on a fixed box-shaped ship in a large harbor.  相似文献   

5.
This paper presents a refined parabolic approximation model of the mild slope equation to simu-late the combination of water wave refraction and diffraction in the large coastal region.The bottom frictionand weakly nonlinear term are included in the model.The difference equation is established with the Crank-Nicolson scheme.The numerical test shows that some numerical prediction results will be inaccurate in com-plicated topography without considering weak nonlinearity;the bottom friction will make wave height damp-ing and it can not be neglected for calculation of wave field in large areas.  相似文献   

6.
Simulation of Fully Nonlinear 3-D Numerical Wave Tank   总被引:6,自引:0,他引:6  
A fully nonlinear numerical wave tank (NWT) has been simulated by use of a three-dimensional higher order bouodary element method (HOBEM) in the time domain. Within the frame of potential flow and the adoption of simply Rankine source, the resulting boundary integral equation is repeatedly solved at each time step and the fully nonlinear free surface boundary conditions are integrated with time to update its position and boundary values. A smooth technique is also adopted in order to eliminate the possible saw-tooth numerical instabilities. The incident wave at the uptank is given as theoretical wave in this paper. The outgoing waves are absorbed inside a damping zone by spatially varying artificial damping on the free surface at the wave tank end. The numerical results show that the NWT developed by these approaches has a high accuracy and good numerical stability.  相似文献   

7.
A depth-integrated, non-hydrostatic model with hybrid finite difference and finite volume numerical algorithm is proposed in this paper. By utilizing a fraction step method, the governing equations are decomposed into hydrostatic and non-hydrostatic parts. The first part is solved by using the finite volume conservative discretization method, whilst the latter is considered by solving discretized Poisson-type equations with the finite difference method. The second-order accuracy, both in time and space, of the finite volume scheme is achieved by using an explicit predictor-correction step and linear construction of variable state in cells. The fluxes across the cell faces are computed in a Godunov-based manner by using MUSTA scheme. Slope and flux limiting technique is used to equip the algorithm with total variation dimensioning property for shock capturing purpose. Wave breaking is treated as a shock by switching off the non-hydrostatic pressure in the steep wave front locally. The model deals with moving wet/dry front in a simple way. Numerical experiments are conducted to verify the proposed model.  相似文献   

8.
Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan.  相似文献   

9.
A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity.The Euler predictor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones,respectively,thus leading both the time and spatial derivatives to the second-order accuracy.The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift.Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach,respectively,and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor.From comparison of numerical results with theoretical or experimental ones,it is found that they are in reasonable agreement.  相似文献   

10.
This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equations. Numerical method for modelling the short waves is investigated in detail. The forces, such as Coriolis forces, wind stress, atmosphere and bottom friction, are considered. A two-dimensional implicit difference scheme of Boussinesq equations is proposed. The low-reflection outflow open boundary is suggested. By means of this model,both velocity fields of circulation current in a channel with step expansion and the wave diffraction behind a semi-infinite breakwater are computed, and the results are satisfactory.  相似文献   

11.
In order to assess the impact of deep-sea mining on the in situ benthic life, we measured the microbial standing stock and concentration of organic nutrients in the deep-sea sediments of the Central Indian Ocean Basin in the Indian pioneer area. Sediments were collected using box core and grab samples during September 1996. The total bacterial numbers ranged from 10 10 -10 11 cells per g -1 dry weight sediment. There was a marginal decrease in the number of bacteria from surface to 30 cm depth, though the subsurface section registered a higher number than did the surface. The highest numbers were encountered at depths of 4-8 cm. The retrievable number of bacteria were two orders less in comparison with the direct total counts of bacteria. An almost homogeneous distribution of bacteria, total organic carbon, living biomass, and lipids throughout the depth of cores indicates active microbial and benthic processes in the deep sea sediments. On the other hand, a uniform distribution of total counts of bacteria, carbohydrates, and total organic carbon in all the cores indicates their stable nature and suggests that they can serve as useful parameters for long-term monitoring of the area after the benthic disturbance. Further studies on temporal variability in this region would not only verify the observed norms of distribution of these variables but would also help to understand restabilization processes after the simulated benthic disturbance.  相似文献   

12.
An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-waveletconvolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L 2 norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.  相似文献   

13.
海上大直径钢管桩打桩过程中,桩周土体受到强烈扰动而发生强度弱化,掌握桩周土体强度弱化规律对于准确预测打桩过程、保证工程安全具有重要意义。为研究土体强度弱化规律,开展了环剪试验模拟打桩对桩周土体的扰动,测试土体强度随剪切速率的变化规律,建立了描述土体强度弱化规律的拟合公式,引入到打桩分析软件中。研究结果表明:土体的强度折减程度不仅与土体本身的性质有关还受到土体的埋深和剪切速率的影响,埋深越深土体强度折减程度越低,剪切速率越高土体强度折减越高,在打桩分析中可采用这里推荐的线性折减方法来模拟不同深度处土体强度的折减规律。  相似文献   

14.
The interdependence between the seismo-acoustic properties of a marine sediment and its geotechnical/physical parameters has been known for many years, and it has been postulated that this should allow the extraction of geotechnical information from seismic data. Though in the literature many correlations have been published for the surficial layer, there is a lack of information for greater sediment depths. In this article, a desktop study on a synthetic seafloor model illustrates how the application of published near-surface prediction equations to subsurface sediments (up to several tens of meters burial depth) can lead to spurious predictions. To test this further, acoustic and geotechnical properties were measured on a number of sediment core samples, some of which were subjected to loading in acoustically-equipped consolidation cells (oedometers) to simulate greater burial depth conditions. For low effective pressures (representing small burial depths extending to around 10 meters subsurface), the general applicability of established relationships was confirmed: the prediction of porosity, bulk density, and mean grain size from acoustic velocity and impedance appears generally possible for the investigated sedimentary environments. As effective pressure increases through, the observed relationships deviate more and more from the established ones for the near-surface area. For the samples tested in this study, in some instances increasing pressure even resulted in decreasing velocities. There are several possible explanations for this abnormal behavior, including the presence of gas, overconsolidation, or bimodal grain size distribution. The results indicate that an appropriate depth correction must be introduced into the published prediction equations in order to obtain reliable estimates of physical sediment properties for greater subsurface depths.  相似文献   

15.
16.
Oedometer tests have been carried out on 70 undisturbed surficial clays (at approximately 250 mm below the mudline), mostly collected by free-fall corers from sites widely scattered throughout the deep-sea North Atlantic. Acoustic measurements were also made, initially on contiguous samples and ultimately on the same sample using a geophysically instrumented oedometer which also collected electrical resistivity data. Apart from those quiescent areas below the carbonate compensation depth, such as north of the West Indies where very fine clays exist, most of the samples are silty clays whose geotechnical-geophysical properties are dependent on the type of clay minerals present (and their ability to take in moisture), the sand-size fraction, and the quantity of carbonate present. Thus the pure clays have high compressibilities which decrease on the addition of coarse particles, while the converse is true for the acoustic parameters, these increasing with the sand fraction. Using the notion of the intrinsic compression line for all samples, and comparison to it of the measured compression curves, it is clear that, contrary to some previously held ideas, most deep-sea clays are normally consolidated; the addition of carbonate has the effect of creating an open, stronger sediment skeleton. Interestingly, where information is available, the variation with depth of a sample's acoustic velocity follows the void ratio pressure relationship of the compression curve. This allows the construction of an in-situ sediment compression curve using the in-situ geophysical observations.  相似文献   

17.
A set of 27 marine planktonic bacteria isolated from the polar regions was characterized by 16S rDNA sequencing and physiological and biochemical testing. More than half of these bacteria were positive for caseinase, gelatinase and 13-glucosidase, and could utilize glucose, maltose or malic acid as carbon source for cell growth. Twelve isolates expressed nitrate reduction activities. Except for one antarctic isolate BSwlO175 belonging to Actinobacteria phylum, these isolates were classified as γ-Proteobacteria, suggesting that γ-Proteobacteria dominated in cultivable marine bacterioplankton at both poles. Genus Pseudoalteromonas was the predominant group in the Chukchi Sea and the Bering Sea, and genus ShewaneUa dominated in cultivable bacterioplankton in the Prydz Bay. With sequence similarities above 97%, genus Psychrobacter was found at both poles. These 27 isolates were psychrotolerant, and significant 16S rDNA sequence similarities were found not only between arctic and antarctic marine bacteria ( 〉 99% ), but also between polar marine bacteria and bacteria from other aquatic environments ( ≥ 98.8% ), including temperate ocean, deep sea, pond and lake, suggesting that in the polar oceans less temperature-sensitive bacteria may be cosmopolitan and have a bipolar, even global, distribution at the species level.  相似文献   

18.
19.
Specific properties of the interannual sea level variations and annual tides in the Northwestern Pacific were studied. Several tide stations were monitored. The monthly mean sea level for the year of 1995 was analyzed at each tide station. A seismic event in 1995, some tectonic activity around the subject area, and the Kuroshio (the oceanic western boundary current) may possibly contaminate results which would have occurred from the astronomical annual tide alone.  相似文献   

20.
Particle fluxes were measured 7 m above the sea bottom during the predisturbance, disturbance, and postdisturbance periods by using time series sediment traps attached to seven deep-sea moorings deployed in the INDEX experiment site in the Central Indian Basin. The predisturbance particle fluxes varied between 22.3 to 55.1 mg m -2 day -1 . Increased and variable particle fluxes were recorded by the sediment traps during the disturbance period. The increase observed was 0.5 to 4 times more than the background predisturbance fluxes. The increases in particle fluxes (~4 times) recorded by the sediment trap located in the southwestern direction (DMS-1) were the greatest, which could be the result of preferential movement of resuspended particles generated during the deep-sea benthic disturbance along the general current direction prevailing in this area during the experimental period. Also, the traps located closer to the disturbance area recorded greater fluxes than did the traps far away, across the Deep Sea Sediment Resuspension System path. This variability in recorded particle fluxes by the traps around the disturbance area clearly indicates that physical characteristics such as grain size and density of the resuspended particles produced during the disturbance had an important effect on particle movement. The postdisturbance measurements during ~5 days showed a reduction in particle fluxes of ~50%, indicating rapid particle settlement.  相似文献   

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

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