首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
运用 Mc Creary和 Yu(1992 )的非线性 2 12 层模式的积分流场 ,研究了混沌混合与输运对海水及要素浓度的输运和混合问题。模式结果显示混沌混合区主要分布于赤道东太平洋不稳定区及东、西边界附近 ,中西太平洋的大部分区域为非混沌区。释放在混沌区内的示踪团随时间的弥散过程要较非混沌区内的示踪团复杂的多。因多模态不稳定波动的共同作用 ,尤其是在赤道上经向运动最强的 Yanai波的加入 ,使得赤道不稳定区的流体以混沌混合的方式均匀地弥散开 ,趋向于粗粒 (coarse- grained)意义下的均匀化。并有大量释放于赤道东太平洋不稳定区的示踪质点越过赤道进入对面半球。背景辐散流的加入能将部分示踪点输送出混合区使其不加入均匀化过程 ,故最后混入另一半球的质点数少于半数。过赤道交换速率则主要由混沌混合的强度有关。  相似文献   

2.
Cyclone-anticyclone asymmetry in a rotating fluid results in vortices with cyclonic rotation being attenuated more rapidly than vortices with anticyclonic rotation due to the Ekman bottom friction. To explain this effect, some authors invoked rather complex integral (averaged along the vertical) models with the parametrization of nonlinear friction. A simple analytical model, free of the procedure of formal averaging and based on a separate consideration of the equations for external flow in the nonviscous region and internal flow in the boundary layer, is investigated in this work. The corresponding equations are written in the so-called geostrophic momentum approximation, which makes it possible to take into account the nonlinear advective mass transport in the boundary layer at small Rossby numbers. The nonlinear equation of the hyperbolic type for the tangential velocity, which describes the process of attenuation of an axisymmetric vortex, is obtained from the condition of total mass conservation. Based on the solutions to this equation, it was shown that distinctions in the character of vortex attenuation are caused by deviations from the geostrophic regime in the nonviscous region. It was established that the concentration (compression) of anticyclonic vortices and the extension of cyclonic ones take place in the process of attenuation.  相似文献   

3.
A method is suggested for simulating axisymmetric laminar or turbulent flows formed during the motion of a vortex-ring bunch of given geometry and circulation toward a plane screen. Earlier, similar problems were simulated with the numerical solution of the Navier-Stokes equations for laminar flows. Turbulent flows have remained unconsidered until now. When a vortex ring approaches the screen, the secondary nonstationary flow is induced near the screen’s surface and this secondary flow causes the formation of the radial boundary layer (provided that air viscosity is taken into account). First, the medium spreads out from the critical point at the screen’s center with the negative pressure gradient along the radial coordinate and then detaches in the region of the positive pressure gradient. This radial wall flow and the corresponding boundary layer are considered in the quasi-stationary approximation. When the boundary layer detaches at successive instances, the flow is replenished with the radially moving secondary vortex rings whose circulations have the sign opposite to that of the circulation of the primary vortex ring. It is the interaction of the primary and secondary vortices that governs process dynamics, which differs substantially from that in the case when the formation of secondary vortices is disregarded. The suggested method is based on the method of discrete vortices (a perfect liquid) and the boundary-layer (laminar or turbulent) theory. During the development of the flow under investigation, the nonstationary ascending flow in the direction perpendicular to the screen’s plane is formed and then this flow decays and dissipates. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer show that the velocity of ascending vortices in the plane of the initial vortex bunch is less than one-tenth of the initial velocity of the descending vortex ring. The boundary layer is introduced into calculations with the sole goal of determining the parameters of the secondary vortex rings formed during boundary-layer detachments. The interaction of the primary and secondary vortices is then considered within the framework of a perfect medium. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer on the screen were correlated with the available data obtained in laboratory experiments for small Reynolds numbers. Qualitative agreement between the simulations and experiments is fairly satisfactory. The simulation for one combination of the circulation and vortex-ring geometry takes at most 10–15 min with the use of an average PC.  相似文献   

4.
Methods from chaos theory are applied to the analysis of the circulation in the Southern Ocean, using velocity fields produced by a realistic global ocean model. We plot the intersections of individual trajectories encircling Antarctica with a vertical plane in the Drake passage. This so-called Poincaré section shows a drastic difference between regular trajectories in a core region of the Antarctic Circumpolar Current (ACC), and chaotic, mixing trajectories in the surrounding region. It also shows that there is a region with overturning circulation of approximately 3.5 Sv in the ACC, with downwelling on the northern side and upwelling on the southern side, which may be related to the Deacon cell.  相似文献   

5.
The effects of chaotic advection in a three-layer ocean model   总被引:1,自引:0,他引:1  
A three-layer model of an inviscid incompressible fluid is considered in a quasigeostrophic approximation within the concept of background currents. A singular topographic obstacle yields the formation of a vortical motion in the layers. It is always present in the lower layer, while in the upper and middle layers it can be observed only at certain velocities of the external flows. Three current functions describing a singular topographic flow in the lower layer and two regular flows in the upper layer are obtained. In case of a nonstationary harmonic perturbation of the external flow, chaotic particle transport is possible in these flows. This paper analyzes the chaotic properties of this model. Depending on the type of unperturbed frequency curves of the fluid particle revolution is determined by the model parameters (stratification) and the incident flow, the patterns of the chaotic transport will substantially differ. Two of the limiting dependences of the revolution frequency are determined, namely, the dependence for a singular vortex and that for a regular vortex with the smallest area. Other dependences are intermediate ones. Two limiting types of revolution frequencies determine the two different scenarios of chaotic advection.  相似文献   

6.
Propagation of a solitary wave over rigid porous beds   总被引:1,自引:0,他引:1  
The unsteady two-dimensional Navier–Stokes equations and Navier–Stokes type model equations for porous flows were solved numerically to simulate the propagation of a solitary wave over porous beds. The free surface boundary conditions and the interfacial boundary conditions between the water region and the porous bed are in complete form. The incoming waves were generated using a piston type wavemaker set up in the computational domain. Accuracy of the numerical model was verified by comparing the numerical results with the theoretical solutions. The main characteristics of the flow fields in both the water region and the porous bed were discussed by specifying the velocity fields. Behaviors of boundary layer flows in both fluid and porous bed regions were also revealed. Effects of different parameters on the wave height attenuation were studied and discussed. The results of this numerical model indicate that for the investigated incident wave as the ratio of the porous bed depth to the fluid depth exceeds 10, any further increase of the porous bed depth has no effect on wave height attenuation.  相似文献   

7.
It is suggested that the downward helicity flux (through the upper boundary of the viscous turbulent boundary layer) be treated as a measure of the intensity of atmospheric vortices, including tropical cyclones, tornadoes, and dust devils. As follows immediately from the general helicity balance equation known in the literature, this flux is determined by the product of the cubed maximum wind speed and the width of the strip swept by the maximum wind during vortex movement. For intense vortices in their steady-state, mature stage, this helicity flux can also serve as a measure of the rate of helicity destruction by the forces of viscous turbulent friction. Examples of applying the introduced notion to the diagnostics of tornadoes and their classification according to a destructive force are given. A comparative analysis (according to helicity flux values) of dust devils on the Earth and Mars, on the one hand, and tornadoes, on the other, is presented.  相似文献   

8.
《Coastal Engineering》2005,52(1):25-42
The performance of the standard kε, Wilcox high-Reynolds-number kω, Wilcox low-Reynolds-number kω and Smagorinsky's subgrid scale (SGS) turbulence models is examined against the flow around a circular cylinder 0.37 diameter above a rigid wall. The governing equations are solved using finite difference method in a non-orthogonal boundary-fitted curvilinear coordinate system. A mesh dependence study for the four turbulence models is carried out on computational meshes with different densities. In addition, the performance of the kω models with either wall function or no-slip boundary condition on the cylinder surface is examined on the finest mesh. It is found that the SGS model over-predicts the shedding of vortices from the cylinder and is sensitive to the computational mesh and the model constant Cs used. The standard kε and the Wilcox kω models predict the mean velocity field quite well but generally under-predict the velocity and hydrodynamic force oscillations using wall functions on the cylinder surface. It is also found that the Wilcox kω models with the no-slip boundary condition on the cylinder surface give better predictions on the shedding of vortices than their counterparts using the wall function boundary condition.  相似文献   

9.
Open boundaries are important when simulating water waves. In this study, a transparent boundary condition at an open boundary was developed for simulating nonlinear water waves propagating to a distant area using the Moving Particle Semi-implicit method. The novelty of this study is that the technique of wave analysis used in the experiment was introduced into the particle simulation to absorb incident waves; the simulation cost was reduced by employing inflow and outflow regions instead of a long dissipation region. Incident waves in front of the boundary were evaluated using Fourier analysis, and the particles on the transparent boundary were forced to move at the velocity of the analytical solution for Stokes waves in order to absorb the incident waves. The analysis was restricted to periodic waves. Wave propagation was simulated for two wave periods using the developed transparent boundary condition. The results showed that this transparent boundary transmitted the incident waves with small reflection and the simulation cost was lower than that for wave damping by a conventional highly viscous region.  相似文献   

10.
Using an eddy-resolving two layer primitive-equation model forced by symmetric wind stress, we investigate the asymmetry of the recirculation of a double gyre (subpolar gyre and subtropical gyre). In the case where the layer thickness change is large, cyclonic recirculation becomes unstable and splits into meso-scale vortices more easily than anti-cyclonic recirculation in their developing stage. The subpolar gyre is, therefore, filled with more vortices than the subtropical gyre. Moreover, the effect of the upper layer potential vorticity on the lower layer potential vorticity in the subpolar gyre is stronger than that in the subtropical gyre. The characteristic of turbulence in the subpolar gyre is different from that found in the subtropical gyre and, therefore, the asymmetry of the recirculation of the double gyre is maintained by this difference. The asymmetry can not be produced in a quasi-geostrophic model because it ignores the nonlinearity associated with layer thickness change. Moreover, we investigate the effects of layer thickness and lateral viscosity on the asymmetry of the recirculation of the double gyre. In the case of realistic physical parameters, the asymmetry of the recirculation of the double gyre is noticeable from the view of the activities of the eddies. In the case with the shallowest upper layer, the position of separation of the western boundary current moves further southward.  相似文献   

11.
The rolling motion of a ship has been successfully modelled using a semi-empirical nonlinear differential equation by a number of researchers. Experimental data has been used to model nonlinear damping and righting lever characteristics and comparison with observed behaviour has been reasonably good. The present article describes a numerical, phenomenological approach to analyse this type of behaviour. The stability of the periodic motion, and in particular the possibility of capsize, is explored with reference to qualitative prediction techniques. The appearance of chaotic motions in regular beam seas is a new feature which should be of interest to designers. The inability of traditional quantitative methods, such as the perturbation technique, to detect chaos is a further justification for using numerical simulation guided by dynamical systems theory.  相似文献   

12.
杨继 《海洋技术学报》2008,27(1):96-100
混沌理论在水文系统中的应用主要集中在水文时间序列的混沌性识别和水文混沌预测的模型研究上。重点探讨了水文系统混沌分析应用领域的研究进展,并分析认为,提取海浪混沌特征,探索海浪预报的一种新方法,很值得尝试并深入研究。  相似文献   

13.
If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.  相似文献   

14.
The objective of the present work is to discuss the implementation of an active wave generating–absorbing boundary condition for a numerical model based on the Volume Of Fluid (VOF) method for tracking free surfaces. First an overview of the development of VOF type models with special emphasis in the field of coastal engineering is given. A new type of numerical boundary condition for combined wave generation and absorption in the numerical model VOFbreak2 is presented. The numerical boundary condition is based on an active wave absorption system that was first developed in the context of physical wave flume experiments, using a wave paddle. The method applies to regular and irregular waves. Velocities are measured at one location inside the computational domain. The reflected wave train is separated from the incident wave field in front of a structure by means of digital filtering and subsequent superposition of the measured velocity signals. The incident wave signal is corrected, so that the reflected wave is effectively absorbed at the boundary. The digital filters are derived theoretically and their practical design is discussed. The practical use of this numerical boundary condition is compared to the use of the absorption system in a physical wave flume. The effectiveness of the active wave generating–absorbing boundary condition finally is proved using analytical tests and numerical simulations with VOFbreak2.  相似文献   

15.
A three-dimensional (3D) numerical model of fixed Oscillating Water Column system (OWC) is presented and validated. The steady-state potential flow boundary value problem due to regular wave interaction with the OWC is solved by a first order mixed distribution panel method. Ocean response predictions are derived using a deterministic statistical model based on a spectral analysis method. The model validation focusses on diffraction predictions and involves convergence tests and numerical comparisons with independent potential flow computations. Predictions of both regular and irregular wave responses are also compared against experimental results. Sample results including the yearly-averaged power conversion efficiency are presented in the final section to illustrate the method’s suitability to a 3D hydrodynamic design optimisation.  相似文献   

16.
The method of slowly varying amplitudes is used to determine the reflection of surface gravity-capillary waves from regions of regular surface convection in the approximation of an ideal and homogeneous viscous fluid. A model of a quasi-steady distribution of temperature is employed. The conditions of both weak (nonresonant) interaction, when the wave passes almost completely through the convective zone, and strong resonant interaction are revealed. When interaction is strong, standing waves form in front of the convective region and the spectral composition of the wave packet changes. Viscosity and the accompanying boundary layers lead to the broadening of forbidden zones.  相似文献   

17.
任冰  高睿  金钊  王国玉  王永学 《海洋学报》2012,34(1):163-177
基于光滑粒子流体动力学(SPH)方法的二维数值波浪水槽模拟了规则波对透空式结构物的冲击作用。通过黎曼解和CSPM相结合的方法对连续方程和动量方程进行了修正。对造波边界采用虚粒子法模拟。提出了一种耦合计算方法来消除造波边界附近的压力波动现象,在结构物边界处设置了适合的冲击边界条件。应用修正的SPH模型模拟了规则波对浪溅区结构物的冲击作用,给出了结构物附近流场和压力场的变化特征,并应用物理模型试验结果对数值模型进行了验证。  相似文献   

18.
A combined method for the hydrodynamic characteristics of planing crafts   总被引:1,自引:0,他引:1  
The prediction of the total resistance of planing crafts at high speeds is very important. In this paper, a combined method is investigated for determining the hydrodynamic characteristics of planing crafts in the calm water. The study consists of a potential-based boundary element method (BEM) for the induced pressure resistance, the boundary layer theory for the frictional resistance and practical method for the spray resistance. The planing surface is represented by a number of elements with constant velocity potential at each element. The unknown-induced pressure is obtained by using the free surface elevation condition and the Kutta condition at the transom stern. Hydrodynamic-induced resistance and lift are determined by the calculated dynamic pressure distributions. The boundary layer analysis method is based on calculations of the momentum integral equation applied to obtain the frictional resistance. A particular practical approach is introduced to present the region of the upwash geometry for the spray. A numerical program has been developed for the present research and applied to the hull form of the craft. Four different hull forms of Series 62 model 4666 planing craft are presented. It is shown that the present combined method is efficient and the results are in good agreement with the experimental measurements over a wide range of volumetric Froude numbers.  相似文献   

19.
刘卓  郭冬建  朱江  曾庆存 《海洋学报》1999,21(4):96-105
随着沿海经济带的开发利用,河口海岸动力学问题,如航道淤积、河口演变、港湾建设及污染物的排放和扩散等,日益引起人们的重视.数值模拟方法是研究这些海岸环境问题的理想工具之一,但由于河口区域的海岸线和河岸线的形状非常复杂,加上众多的岛屿,给模式网格的划分造成很大的困难.为此,我们将自适应网格的生成方法应用于河口港湾的数值模拟研究,以长江口和杭州湾的水动力学模拟为例,该法较好地克服了均匀网格拟合边界的困难,得到了较合理的计算结果.  相似文献   

20.
利用运动学模型 ,研究了弯曲急流控制参数为定常和发生周期扰动时的流体质点运动轨迹以及混沌现象的产生机制。发现弯曲急流控制参数为定常时 ,流体质点的运动为规则的周期运动 ;而控制参数发生扰动时 ,可以发生流体质点在不同特征区的穿行 ,流体质点的运动具有混沌特性。弯曲急流的相速、振幅发生周期扰动均可使流体质点的运动出现混沌特性。作者从运动学上对弯曲急流区的漂流浮标混沌特性给出了解释。  相似文献   

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

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