共查询到20条相似文献，搜索用时 31 毫秒
1.
《Coastal Engineering》2005,52(1):43-62
A vertical two-dimensional (2D) numerical model for time dependent local scour below offshore pipelines subject to unidirectional steady flow is developed. The governing equations for the flow and sediment transport are solved by using finite difference method in a general curvilinear coordinate system. The performance of two turbulence models, the standard k–ɛ model and Smagorinsky subgrid scale (SGS) model, on modeling time dependent scour processes is examined. Both suspended load and bed load are considered in the scour model. The suspended-load model is verified against two channel sediment transport cases. The change of bed level is calculated from the continuity equation of total sediment transport. A new time marching scheme and a sand slide scheme are proposed for the scour calculation. It is found that the proposed time marching scheme and sand slide model work well for both clear-water and live-bed scour situations and the standard k–ɛ turbulence closure is more preferable than the SGS model in the 2D scour model developed in this study. 相似文献
2.
The flow structure of a swash event over a uniform slope is studied using a RANS-VOF numerical model coupled with a v2–f turbulence closure. The model is compared with experimental data of recent laboratory experiments. The ability of the turbulence modelling for simulating swash flow and the evolution of the computed bed shear stress during run-up and run-down are investigated. The agreement between numerical results and measured data, such as water depth, depth-averaged velocity and bed shear stress is very good during run-up. Main discrepancies are found during run-down. The paper also examines the aeration of the water layer in the swash flow, taking advantage of the PLIC method for computing the air–water interfaces. Air is continuously entrapped in the swash front and released at its rear during run-up. A detailed analysis indicates that the flow reversal is initiated near the bottom at the outer boundary of the swash zone and progresses landward. The study highlights the asymmetry between run-up and run-down. During run-up, the swash front propagation determines the turbulence properties and the bed shear stress profile on the beach, whereas the flow properties are more homogeneously distributed in the swash area during run-down. 相似文献
3.
The submerged 3D turbulent jet flow behavior around a pile on a rigid bed and on a scoured bed was studied experimentally and numerically. ADV was used to obtain the jet velocity distributions and Realizable k–ε turbulence model was used for the prediction of flow field around a pile. The jet flow area was three-dimensional and thus numerical model was a three-dimensional model. The numerical results were used to predict the velocity close to the pile and bed shear stress on the bed. In general, the results indicated that the flow field was also in agreement with the findings of previous experiments in literature and the related principles in the subject area. The experimental results demonstrated that Acoustic Doppler Velocimeter (ADV) measurements were almost identical with the Realizable k–ε turbulence model results for turbulence intensity I=10%. 相似文献
4.
A numerical model is developed to predict the onset of local scour below offshore pipelines in steady currents and waves. The scour is assumed to start when the pressure gradient underneath the pipeline exceeds the floatation gradient of the sediments. In this model, the water flow field above the bed is determined by solving the two-dimensional (2-D) Reynolds-averaged Navier–Stokes equations with a k-ω turbulence closure. The seepage flow below the seabed is calculated by solving the Darcy's law (Laplace's equation) with known pressure distribution along the common boundaries of the flow domains-seabed. The numerical method used for both the turbulent flow around the pipeline and Darcy's flow in the seabed is a fractional finite element method. The average pressure gradient along the buried pipe surface is employed in the evaluation of onset condition with a calibration coefficient. The numerical model is validated against experimental data available in literature. A unified onset condition for steady currents and waves is proposed. Influences of flow parameters, including water depth, embedment depth, boundary layer thickness, Reynolds number (Re) and Keuleagan–Carpenter (KC) number, on the pressure drop coefficient over the pipeline are studied systematically. 相似文献
5.
Wave–current flow is a phenomenon that is present in many practical engineering situations. Over the past several decades, this type of flow has been increasingly investigated under controlled laboratory conditions. This paper presents a numerical study of wave–current flow in the ocean basin of the LabOceano (COPPE/UFRJ). A homogeneous multiphase model based on the RANS equations and the k–ɛ turbulence model implemented in ANSYS-CFX code were used. A cross section of the ocean basin was represented. A regular wave with a height of 0.08 m and a period of 1.80 s (i.e., a wave steepness of H/L = 0.016), propagating on favourable currents, was simulated. The behaviour of the free surface elevation over time and the streamlines along the basin for wave and wave–current flows were presented. The numerical results were compared to the non-viscous theory given by the Rayleigh equation applied to the problem of wave–current interaction. Good agreement was found between the wave length estimated by the numerical results and the analytical solutions, with a deviation of less than 2%. 相似文献
6.
A new coupling model of wave interaction with porous medium is established in which the wave field solver is based on the two dimensional Reynolds Averaged Navier-Stokes (RANS) equations with a kε closure. Incident waves, which could be linear waves, cnoidal waves or solitary waves, are produced by a piston-type wave maker in the computational domain and the free surface is traced through the Piecewise Linear Interface Construction-Volume of Fluid (PLIC-VOF) method. Nonlinear Forchheimer equations are adopted to calculate the flow field within the porous media. By introducing a velocity–pressure correction equation, the wave field and the porous flow field are highly and efficiently coupled. The two fields are solved simultaneously and no boundary condition is needed at the interface of the internal porous flow and the external wave. The newly developed numerical model is used to simulate wave interaction with porous seabed and the numerical results agree well with the experimental data. The additional numerical tests are also conducted to study the effects of seabed thickness, porosity and permeability coefficient on wave damping and the pore water pressure responses. 相似文献
7.
Liquid sloshing in storage tank is a fundamental problem of great engineering importance. Sloshing motion can be laminar or turbulent. However, the necessity for inclusion of turbulence in CFD simulation of sloshing flows has not yet been established. In this paper, three roll–induced sloshing cases are studied to assess the merits and shortcomings of the laminar model and three most–commonly used turbulence models (RANS k–ε, LES and Very LES). To overcome the deficiencies in the RANS and LES, the new Very LES (VLES) model, which combines the RANS k–ε and LES, is developed in this paper. The free surface profiles are reconstructed by a coupled Level–Set and Volume–of–Fluid (CLSVOF) method. To the authors’ knowledge, the comprehensive and systematical assessment of the effect of turbulence on sloshing simulation has not been reported in the literature. The numerical results are evaluated using experimental measurements from Delorme and Souto−Iglesias. The present study indicates that the inclusion of an appropriate turbulence model has a profound influence on the simulations of violent and non–violent sloshing flows. The VLES and LES models can provide accurate predictions of free surface profiles and impact pressures, whereas the laminar flow assumption and the RANS model cannot adequately capture the energy dissipation in the sloshing simulation and lead to the inaccurate flow predictions. 相似文献
8.
《Coastal Engineering》2001,44(1):13-36
Interactions between a solitary wave and a submerged rectangular obstacle are investigated both experimentally and numerically. The Particle Image Velocimetry (PIV) technique is used to measure the velocity field in the vicinity of the obstacle. The generation and evolution of vortices due to flow separation at the corners of the obstacle are recorded and analyzed. It is found that although the size of the vortex at the weatherside of the obstacle is smaller than that at the leeside, the turbulence intensity is, however, stronger. A numerical model, based on the Reynolds Averaged Navier–Stokes (RANS) equations with a k–ϵ turbulence model, is first verified with the measurements. Overall, the agreement between the numerical results and laboratory velocity measurements is good. Using the RANS model, a series of additional numerical experiments with different wave heights and different heights of the rectangular obstacle are then performed to test the importance of the energy dissipation due to the generation of vortices. The corresponding wave transmission coefficient, the wave reflection coefficient and the energy dissipation coefficient are calculated and compared with solutions based on the potential flow theory. As the height of the obstacle increases to D/h=0.7, the energy dissipation inside the vortices can reach nearly 15% of the incoming wave energy. 相似文献
9.
Chunyan Li Nan Walker Aixin Hou Ioannis Georgiou Harry Roberts Ed Laws J. Alex McCorquodale Eddie Weeks Xiaofeng Li Jessica Crochet 《Estuarine, Coastal and Shelf Science》2008
Circular shaped density plumes of low turbidity, low fecal indicator (Escherichia coli and enterococci) concentrations, and high salinity have been observed near the Industrial Canal in Lake Pontchartrain, north of the City of New Orleans. A conceptual model in polar coordinates and a numerical model are developed, together with data analysis, to illustrate the dense plume. It is demonstrated that the northward expansion of the plume occurs under northerly winds. The northward expansion of the plume occurs under northerly winds that drive downwind flow at the surface and upwind radial flow at the bottom. Northerly wind-induced straining, similar to tidal straining, promotes vertical stratification. As a result, the water becomes stratified near a thin bottom layer (<1 m), within which density currents are facilitated. The stability of the stratified plume suppresses wind-induced turbulent mixing inside the plume. The bottom water outside of the plume is more effectively stirred by the wind, the result being that the suspended sediment concentration outside of the plume area is much higher than inside. This contrast in mixing makes the plume visible from the surface by satellites even though the stratification is at the bottom. Laterally, wind stress produces a torque (vorticity) in areas of non-uniform depth such that upwind flow is developed in deep water and downwind flow in shallow water. The continuity requirement produces an upwind flow along the axis of the Industrial Canal (IC). The upwind flow is balanced by the downwind flow over the shallower peripheral areas along the coast. 相似文献
10.
Circulation in Tasman Bay 总被引：4，自引：4，他引：0
R. A. Heath 《新西兰海洋与淡水研究杂志》2013,47(3):469-483
Direct current measurements at four locations in Tasman Bay and numerical model results are used to analyse the mean flow in Tasman Bay. The mean circulation conforms to that previously found from drift card experiments: a clockwise circulation in Golden Bay, and an anti‐clockwise flow in Tasman Bay, with a return south‐westerly flow on the coast near Nelson. Typical mean speeds are 0.02–0.05 m.s‐1. The circular flow appears asymmetrical in both bays, with a stronger outflow along Farewell Spit in Golden Bay and near D'Urville Island in Tasman Bay. An analytical tidal solution is used to exhibit the influence of Cook Strait in producing smaller tidal amplitudes in eastern Tasman Bay. Tidal speeds of 0.15–0.30 nus‐1 are typical, with tidal ellipses having degenerated into north‐east, south‐west lines. 相似文献
11.
A computational fluid dynamics (CFD) model was developed to simulate the turbulent flow and species transport of deep-sea high temperature hydrothermal plumes. The model solves numerically the density weighted unsteady Reynolds-averaged Navier–Stokes equations and energy equation and the species transport equation. Turbulent entrainment and mixing is modeled by a k–ε turbulence closure model. The CFD model explicitly considers realistic vent chimney geometry, vent exit fluid temperature and velocity, and background stratification. The model uses field measurements as model inputs and has been validated by field data. These measurements and data, including vent temperature and plume physical structure, were made in the ABE hydrothermal field of the Eastern Lau Spreading Center. A parametric sensitivity study based on this CFD model was conducted to determine the relative importance of vent exit velocity, background stratification, and chimney height on the mixing of vent fluid and seawater. The CFD model was also used to derive several important scalings that are relevant to understanding plume impact on the ocean. These scalings include maximum plume rise height, neutrally buoyant plume height, maximum plume induced turbulent diffusivity, and total plume vertically transported water mass flux. These scaling relationships can be used for constructing simplified 1-dimensional models of geochemistry and microbial activity in hydrothermal plumes. Simulation results show that the classical entrainment assumptions, typically invoked to describe hydrothermal plume transport, only apply up to the vertical level of ~0.6 times the maximum plume rise height. Below that level, the entrainment coefficient remains relatively constant (~0.15). Above that level, the plume flow consists of a pronounced lateral spreading flow, two branches of inward flow immediately above and below the lateral spreading, and recirculation flanking the plume cap region. Both turbulent kinetic energy and turbulence dissipation rate reach their maximum near the vent; however, turbulent viscosity attains its maximum near the plume top, indicating strong turbulent mixing in that region. The parametric study shows that near vent physical conditions, including chimney height and fluid exit velocity, influence plume mixing from the vent orifice to a distance of ~10 times the vent orifice diameter. Thus, physical parameters place a strong kinetic constraint on the chemical reactions occurring in the initial particle-forming zone of hydrothermal plumes. 相似文献
12.
The results of joint analysis of temperature variations near mesopause from long-term measurements at the Zvenigorod Scientific Station of the Obukhov Institute of Atmospheric Physics RAS in 1960–2015 and variations of surface air temperature characterizing global climate change. Together with variations of temperature at the mesopause T ms from measurements of the hydroxyl emissions we analyzed the temperature variations near mesopause T m reduced to the same level of solar activity. The observed strong decrease in temperature near mesopause during last decades, particularly in winter, with its tendency to slow down since the 1980’s is was detected against the background of general increase in the surface air temperature of the Northern Hemisphere T NHs and the Earth as a whole. It was revealed a sharp drop in winter temperature near mesopause in 1970s. and its synchronicity with the shift in climatic features at the surface associated with changes in formation of El Nino events and their impact on the global climate. The general significant negative correlation of temperature variations near mesopause and T NHs detected from 56-year observational data was not accompanied by any significant coherence between the most long-period temperature variations from the cross-wavelet analysis. To assess the possible manifestation of this coherence the results of numerical simulations with a global climate model were used. According to model simulations for the 20–21 centuries taking into account anthropogenic forcings for significant coherence between long-term variations T m and T NHs the prolonged observations are required for temperature near mesopause–about a century or more. 相似文献
13.
14.
In this paper, a well-developed numerical model based on the immersed boundary (IB) method is used to study oscillatory flows over a bed with large-amplitude ripples in a systematic manner. The work shows that the complex flow over the rippled bed can be numerically dealt with in Cartesian coordinate by the IB method and that the IB method is able to provide main features of the flows near the ripples. An accurate simulation of vortices generation as a result of flow separation at the rippled bed is obtained. It is found that the oscillatory flows start to separate during the flow deceleration when the Keulegan–Carpenter (KC) number is small. The steady streaming for various ripple steepness is simulated and the criterion for separating the single and double structure streaming is also discussed. Moreover, a new type of steady streaming which consists of a pair of embedded recirculations in the vicinity of the ripple trough is obtained for relatively steep ripples in this work. The numerical results, including the steady streaming in particular, may be helpful to improve the understanding of the sediment transport and the seabed evolution with natural ripples under sea waves. 相似文献
15.
16.
To obtain a better understanding of the oscillatory soil liquefaction around an offshore pipeline, a three-dimensional integrated model for the wave–seabed–pipeline interaction (WSPI) is proposed by combining the Reynolds-Averaged Navier–Stokes equations for flow simulations and the dynamic Biot’s equation (“u-p” approximation) for the poro-elastic seabed model. Compared with previous investigations, the wave–current interaction is included in the present WSPI system. At a given time step, the wave pressure extracted from the flow model is applied on the seabed surface to determine the corresponding oscillatory seabed response around an offshore pipeline. The integrated numerical model is first validated using previous laboratory experiments. Then, a parametric study is conducted to examine the effects of flow obliquity and pipeline burial depth on the soil response around an offshore pipeline. Numerical results indicate that the soil under the pipeline is more susceptible to liquefaction at a reduced flow obliquity and pipeline burial depth. Moreover, the liquefaction depth in the case where the wave travels along the current can increase by 10%–30% compared to that in the case where the wave travels against the current, when the magnitude of the current velocity is 1 m/s. 相似文献
17.
A smoothed particle hydrodynamic (SPH) model is developed to simulate wave interaction with porous structures. The mean flow outside the porous structures is obtained by solving Reynolds Averaged Navier–Stokes (RANS) equations and the turbulence field is calculated by a large eddy simulation (LES) model. The porous flow is described by the spatially averaged Navier–Stokes type equations with the resistance effect of the porous media being represented by an empirical frictional source term. The interface boundaries between the porous flow and the outside flow are modeled by means of specifying a transition zone along the interface. The model is validated against other available numerical results and experimental data for wave damping over porous seabed with different levels of permeability. The validated model is then employed to investigate wave breaking over a submerged porous breakwater and good agreements between the SPH model results and the experimental data are obtained in terms of free surface displacement. In addition the predicted velocity, vorticity and pressure fields near the porous breakwater and in the breaking wave zone are also analyzed. 相似文献
18.
A numerical wave tank is first established using the Navier–Stokes equations and the VOF method assuming laminar flow. The standard k–ε, realizable k–ε and RNG k–ε turbulent models are then incorporated to the numerical tank. An effective numerical method for wave absorption utilizing the energy-dissipating property of porous media is also included. To validate the accuracy of the proposed models, the propagation of a solitary wave, where analytical solution is available for comparison, is first simulated. This is followed by the simulation of irregular wave runup on a composite seawall, wave propagation over submerged bars and wave refraction and diffraction over an elliptic shoal, where experimental data are available for comparison. All computed results agree well with either the analytical solution or the experimental data. 相似文献
19.
20.
Ping Zeng Huiquan Chen Baichuan Ao Ping Ji Xiaojian Wang Zhiliang Ou 《Coastal Engineering》2002,44(4)
A numerical and experimental modeling is presented for studying the transport of waste heat from a nuclear power plant into coastal water by using a full-field physical model with scale distortion, a local physical model with normal scale and a depth-averaged k− turbulence model with a modified second-order upwind scheme. Field investigations are also used to provide the calibration and validation of the modeling. A case study simulating the turbulent tidal flow and waste heat transport in the coastal water near Daya Bay Nuclear Power Plant in Southern China was conducted. The experimental result of the case study shows that the water temperature in coastal water was a little oversimulated near the surface and was a little undersimulated near the bottom of heated-water layer by the full-field physical model. The numerical study shows that the depth-averaged k− turbulence model presented a satisfied prediction of turbulent tidal flow and far-field temperature distribution in coastal water, although the near-field stratification due to the heated water effluent was not accounted for. The result of the effect of scale distortion on physical model shows that a full-field physical model with a scale distortion of 10 produced a satisfied result of temperature distribution in the present case study. 相似文献