共查询到20条相似文献，搜索用时 754 毫秒
1.
The aim of this paper is to evaluate the accuracy, stability and efficiency of the overset grid approach coupled with the RANS (Reynolds Averaged Navier-Stokes) model via the benchmark computations of flows around a stationary smooth circular cylinder. Two dimensional numerical results are presented within a wide range of Reynolds numbers (6.31 × 104 ∼ 7.57 × 105) including the critical flow regime. All the simulations are carried out using the RANS solver pimpleFoam provided by OpenFOAM, an open source CFD (Computational Fluid Dynamics) toolkit. Firstly, a grid convergence study is performed. The results of the time-averaged drag and lift force coefficients, root-mean square value of lift force coefficient and Strouhal number (St number) are then compared with the experimental data. The velocity, vorticity fields and pressure distribution are also given. One main conclusion is that the numerical solutions in regard to a fixed cylinderare not deteriorated due to the implementation of the overset grid. Furthermore, it can be an appealing approach to facilitate simulations of Vortex Induced Vibrations (VIV), which involves grid deformation. The present study is a good start to implement the overset grid to solve VIV problems in the future. 相似文献
2.
When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 104 and 7.57 × 105. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem. 相似文献
3.
In the present study, locomotion of a real longfin inshore squid (Doryteuthis pealeii) was numerically investigated. Geometry of a real squid was obtained using computed tomography (CT) images. In addition to a two-dimensional axisymmetric squid model, a modified squid model with no cavities around her head and an ellipse shaped model were generated with a fineness ratio (the ratio of body length to maximum body diameter) of 7.56. These numerical models were exposed to an acceleration with two different velocity programs. Added mass and basset forces on bodies were calculated during acceleration of the squid models starting from rest. Pressure and viscous drag forces were also calculated due to pressure variation along the squid models and friction on the surface of the models. The effect of a nozzle diameter on jet velocities and propulsive efficiency at all bodies were evaluated when time dependent velocity profiles (from 0 to 10 m/s in 0.5 and 1 s time durations) were set for the inlet of computational domain. The modified squid model required least thrust force during acceleration phase of time dependent velocity profile compared to the other models while the 0.02 m nozzle diameter provided largest propulsive efficiency for all models. 相似文献
4.
Most international discussions on ballast water management onboard ships have been virtually on the prevention of the transport of biotic alien species. This study focuses on how to improve this limited knowledge regarding the transport of butyltin compounds via ballasting. Butyltin compounds in ballast water, extracted from three sea-going merchant ships with similar voyage routes and berth ports but with different ballasting practices, were determined. Only monobutyltin (8–57 ng l?1) and dibutyltin (10–32 ng l?1) onboard two general cargo ships were detected in the ballast water taken outside the ports. In contrast, significant tributyltin (23–93 ng l?1) was detected onboard the container ship in the ballast water pumped inside the loading/unloading ports. The comparison results from this study indicate that the risk of tributyltin transfer may be effectively prevented by ballasting outside the ports. 相似文献
5.
The ∼8.15 ka Storegga submarine slide was a large (∼3000 km3), tsunamigenic slide off the coast of Norway. The resulting tsunami had run-up heights of around 10–20 m on the Norwegian coast, over 12 m in Shetland, 3–6 m on the Scottish mainland coast and reached as far as Greenland. Accurate numerical simulations of Storegga require high spatial resolution near the coasts, particularly near tsunami run-up observations, and also in the slide region. However, as the computational domain must span the whole of the Norwegian-Greenland sea, employing uniformly high spatial resolution is computationally prohibitive. To overcome this problem, we present a multiscale numerical model of the Storegga slide-generated tsunami where spatial resolution varies from 500 m to 50 km across the entire Norwegian-Greenland sea domain to optimally resolve the slide region, important coastlines and bathymetric changes. We compare results from our multiscale model to previous results using constant-resolution models and show that accounting for changes in bathymetry since 8.15 ka, neglected in previous numerical studies of the Storegga slide-tsunami, improves the agreement between the model and inferred run-up heights in specific locations, especially in the Shetlands, where maximum run-up height increased from 8 m (modern bathymetry) to 13 m (palaeobathymetry). By tracking the Storegga tsunami as far south as the southern North sea, we also found that wave heights were high enough to inundate Doggerland, an island in the southern North Sea prior to sea level rise over the last 8 ka. 相似文献
6.
We describe experiments with multi-directional focused waves interacted with a vertical circular cylinder in a 3D wave basin. The focus of this study is on the run-up of multi-directional focused waves, wave forces, and wave pressures on the cylinder. Part I, the study on wave run-up, has already been presented by Li et al. (2012). In this paper, the analysis of the wave force on the vertical cylinder is presented.In this experiment, a cylinder with 0.25 m in diameter was adopted and different wave parameters, such as focused wave amplitude, peak frequency, frequency bandwidth and directional spreading index, are considered. The model scale kpa (kp is the wave number corresponding to peak frequency, a is the radium of the cylinder) varies from 0.32 to 0.65. The maximum forces of multi-directional focused wave on cylinder were measured and investigated. The results showed that the wave parameters have a significant influence on the wave force, and that the spatial profile of the surface of multi-directional focused wave can also affect its force on the cylinder, which is different from two-dimensional wave. In addition, the ‘secondary loading cycle’ phenomenon was also observed and discussed. In our experiments, the ‘secondary loading cycles’ occur when kA > 0.36 for all cases. While in some referred small scale experiments, the secondary load cycles are observed even for kA = 0.2, when the waves are longer enough. To larger model scale, the pronounced secondary load cycle occurs with larger wave steepness waves. 相似文献
7.
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%. 相似文献
8.
《Ocean Modelling》2010,31(4):241-255
This study demonstrates the sensitivity of the near-surface properties in the tropical Atlantic Ocean to the high-frequency of the winds in numerical simulations. At intra-seasonal timescales (2–50 days), two distinct period ranges dominate the variability in the upper ocean: periods between 2 and 20 days, which are essentially wind-forced and periods between 20 and 50 days, due mostly to Tropical Instability Waves (TIWs). Using a numerical model forced by different wind fields, it is shown that the characteristics of the intra-seasonal variability in the ocean surface mixed-layer are strongly dependent on the wind forcing. Submonthly winds are shown to force large variability in the upper ocean that can strikingly decrease the amplitude of the TIWs in the mixed-layer and their imprint on the horizontal distribution of sea surface temperatures. Wind products containing too much energy at submonthly periods thus prevent wind-forced simulations from reproducing a realistic surface signature of TIWs, when compared to satellite observations of sea surface temperature. In addition, submonthly wind variability may be responsible for part of the observed interannual variability of the TIW signature in the temperature. The impact of submonthly winds is strongest in the mixed-layer: beneath the mixed-layer, all simulations show similar characteristics of the TIWs. 相似文献
9.
The forced constant acceleration exit of two-dimensional bodies through a free-surface is computed for various 2D bodies (symmetric wedges, asymmetric wedges, truncated wedges and boxes). The calculations are based on the fully non-linear time-stepping complex-variable method of Vinje and Brevig (1981). The model was formulated as an initial boundary-value problem (IBVP) with boundary conditions specified on the boundaries (dynamic and kinematic free-surface boundary conditions) and initial conditions at time zero (initial velocity and position of the body and free-surface particles). The formulated problem was solved by means of a boundary-element method using collocation points on the boundary of the domain and stepped forward in time using Runge–Kutta and Hamming predictor–corrector methods. Numerical results for the deformed free-surface profile, pressure along the wetted region of the bodies and force experienced by the bodies are given for the exit. The analytical added-mass force is presented for the exit of symmetric wedges and boxes with constant acceleration using conformal mappings. To verify the numerical results, the added-mass force and the numerical force are compared and give good agreement for the exit of a symmetric wedge at a time zero (t = 0) as expected but only moderate agreement for the box. 相似文献
10.
Prey selection and knowledge of the amounts of water processed by the early stages of the common jellyfish Aurelia aurita may at certain times of the year be crucial for understanding the plankton dynamics in marine ecosystems with mass occurrences of this jellyfish. In the present study we used two different methods (“clearance method” and “ingestion-rate method”) to estimate the amount of water cleared per unit of time of different types and sizes of prey organisms offered to A. aurita ephyrae and small medusae. The mean clearance rates of medusae, estimated with Artemia sp. nauplii as prey by both methods, agreed well, namely 3.8 ± 1.4 l h? 1 by the clearance method and 3.2 ± 1.1 l h? 1 by the ingestion-rate method. Both methods showed that copepods (nauplii and adults) and mussel veligers are captured with considerably lower efficiency, 22 to 37% and 14 to 30%, respectively, than Artemia salina nauplii. By contrast, the water processing rates of ephyrae measured by the clearance method with A. salina nauplii as prey were 3 to 5 times lower than those measured by the ingestion-rate method. This indicates that the prerequisite of full mixing for using the clearance method may not have been fulfilled in the ephyrae experiments. The study demonstrates that the predation impact of the young stages of A. aurita is strongly dependent on its developmental stage (ephyra versus medusa), and the types and sizes of prey organisms. The estimated prey-digestion time of 1.3 h in a steady-state feeding experiment with constant prey concentration supports the reliability of the ingestion-rate method, which eliminates the negative “container effects” of the clearance method, and it seems to be useful in future jellyfish studies, especially on small species/younger stages in which both type and number of prey can be easily and precisely assessed. 相似文献
11.
《Coastal Engineering》2006,53(11):929-945
A finite difference model based on a recently derived highly-accurate Boussinesq-type formulation is presented. Up to the third-order space derivatives in terms of the velocity variables are retained, and the horizontal velocity variables are re-formulated in terms of a velocity potential. This decreases the total number of unknowns in two horizontal dimensions from seven to five, simplifying the implementation, and leading to increased computational efficiency. Analysis of the embedded properties demonstrates that the resulting model has applications with errors of 2 to 3% for (wavenumber times depth) kh ≤ 10 in terms of dispersion and kh ≤ 4 in terms of internal kinematics. The stability and accuracy of the discrete linearised systems are also analysed for both potential and velocity formulations and the advantages and disadvantages of each are discussed. The velocity potential model is then used to study physically demanding problems involving highly nonlinear wave run-up on a bottom-mounted (surface-piercing) plate. New cases involving oblique incidence are considered. In all cases, comparisons with recent physical experiments demonstrate good quantitative accuracy, even in the most demanding cases, where the local wave steepness can exceed (waveheight divided by wavelength) H / L = 0.20. The velocity potential model is additionally shown to have numerical advantages when dealing with wave–structure interactions, requiring less smoothing around exterior structural corners. 相似文献
12.
Burak Aydog˘an Berna Ayat Mehmet N. Öztürk Esin Özkan Çevik Yalçın Yüksel 《Ocean Engineering》2010,37(5-6):443-453
Owing to their complex character, modeling flow patterns of narrow straits has always been a challenge, even with the numerical techniques of today. This study was aimed at predicting vertical current profiles of a given point in a narrow strait, the Strait of Istanbul. On account of the speed and simplicity it offers, and of its remarkable success in solving complex problems, the feed forward back propagation (FFBP) artificial neural network (ANN) technique was chosen for this study. The model was built on 7039 hours of concurrent measurements of current profiles, meteorological conditions, and surface elevations. The model predicted 12 outputs of East and North velocity components at different depths in a given location. Various alternative models with different inputs and neuron numbers were evaluated attaining the best model by trial and error. Predictions from proposed ANN model were in accordance with the observations with average root mean square error of 0.16 m/s. The same input parameters were then used to build models that predicted current velocities 1–12 h into the future. Results of these predictions show good overall agreement with observations and that FFBP ANN can be used as a reliable tool for forecasting current profiles in straits. 相似文献
13.
The main purpose of this study is to establish a better understanding of the relationship between drag reduction and surface roughness. Experiments were conducted to measure the force and flow characteristics of a circular cylinder with different types of artificial surface roughness over the range 6 × 103 < Re < 8 × 104 (Re is based on the cylinder diameter D). The roughness cylinder was formed by covering the exterior surface of the cylinder with uniformly distributed (1) sandpaper, (2) netting, and (3) dimples. The roughness coefficient ranged from k/D = 0.0028 to 0.025 (k is the roughness height). A detailed quantitative measurement of the flow field around the cylinder using Particle Imaging Velocimetry (PIV) was carried out. The hydrodynamic force coefficients (drag and lift) of the rough cylinders are compared against those of a smooth cylinder measured under the same flow conditions. It is found that certain configuration of surface roughness significantly reduces the mean drag coefficient of the cylinder, particularly at large Reynolds numbers. In addition, the root-mean-square (r.m.s.) lift coefficient of the rough cylinders is considerably lower than that of a smooth cylinder. 相似文献
14.
A horizontal, circular cylinder fitted with one bilge keel is forced to rotate harmonically around its axis. The bilge keel load and hull pressure distribution are investigated. A fully submerged condition (infinite fluid), and three partly-submerged conditions are considered. A two-dimensional numerical study is performed, and the results are validated against recently published experimental data by van’t Veer et al. [30]. In addition, comparisons for mass and drag coefficients are also made with experimental data for plate in infinite fluid (Keulegan and Carpenter [8]), and wall-mounted plate (Sarpkaya and O’Keefe [9]) in oscillatory flow.A Navier–Stokes solver based on the Finite Volume Method is adopted for solving laminar flow of incompressible water. The free-surface condition is linearized by neglecting the nonlinear free-surface terms and the influence of viscous stresses in the free surface zone, while the body-boundary condition is exact. This simplified modeling of the problem required the mesh to be fine only around the bilge keels, leading to a total number of cells around N ≃ 1 ×104, which reduced computational cost significantly.The influence of draft and amplitude of oscillations on the bilge keel force and hull pressure distribution are considered. The bilge keel force is presented in terms of non-dimensional drag and mass coefficients including higher harmonic components. The numerical results are also compared with the industry standard empirical method for calculation of roll damping proposed by Ikeda et al. [4]. In general, a good agreement between the results of the present numerical method and the experimental data is obtained and the differences with those predicted by the empirical method are addressed. 相似文献
15.
Shandong province is located on the east coast of China and has a coastline of about 3100 km. There are only a few tsunami events recorded in the history of Shandong Province, but the tsunami hazard assessment is still necessary as the rapid economic development and increasing population of this area. The objective of this study was to evaluate the potential danger posed by tsunamis for Shandong Province. The numerical simulation method was adopted to assess the tsunami hazard for coastal areas of Shandong Province. The Cornell multi-grid coupled tsunami numerical model (COMCOT) was used and its efficacy was verified by comparison with three historical tsunami events. The simulated maximum tsunami wave height agreed well with the observational data. Based on previous studies and statistical analyses, multiple earthquake scenarios in eight seismic zones were designed, the magnitudes of which were set as the potential maximum values. Then, the tsunamis they induced were simulated using the COMCOT model to investigate their impact on the coastal areas of Shandong Province. The numerical results showed that the maximum tsunami wave height, which was caused by the earthquake scenario located in the sea area of the Mariana Islands, could reach up to 1.39 m off the eastern coast of Weihai city. The tsunamis from the seismic zones of the Bohai Sea, Okinawa Trough, and Manila Trench could also reach heights of >1 m in some areas, meaning that earthquakes in these zones should not be ignored. The inundation hazard was distributed primarily in some northern coastal areas near Yantai and southeastern coastal areas of Shandong Peninsula. When considering both the magnitude and arrival time of tsunamis, it is suggested that greater attention be paid to earthquakes that occur in the Bohai Sea. In conclusion, the tsunami hazard facing the coastal area of Shandong Province is not very serious; however, disasters could occur if such events coincided with spring tides or other extreme oceanic conditions. The results of this study will be useful for the design of coastal engineering projects and the establishment of a tsunami warning system for Shandong Province. 相似文献
16.
The purpose of this research work is to study the effect of specific surface s, the fluid–solid contact surface per volume unit, on the wave energy dissipation by porous structures consisting in dense arrays of emergent vertical cylinders. Experiments have been carried out in a 10 m long wave flume. Three cylinder diameters D are considered in order to study the effects of the specific surface while keeping the porosity constant. In a first series, the length of the porous zone is kept constant for the three cylinder diameters tested. The measurements, which include various wave steepness conditions, demonstrate the role of specific surface s on both wave attenuation and interference processes. The larger the specific surface is, the stronger the wave damping is. Damping is found to be almost proportional to 1/D when laminar, turbulent and inertial effects are of same order. Results are compared to numerical calculations based on either a constant rate of wave damping within the porous medium per unit wavelength or a quadratic damping developed using a force expression based on the work of [26]. This latter model, calibrated with drag and inertia coefficients, shows a good agreement with measurements. In a second series, both porous length and water depth are kept proportional to the cylinder diameter for the three diameters. Scale effects are then discussed and underline the importance of the flow regime within the porous medium. 相似文献
17.
Rocky shore ecosystems are considered sentinels of climate warming because they are in close contact with the atmosphere and their shallow waters present low thermal inertia. Concerns on the vulnerability of rocky shore species subject to climate warming make the investigation of their thermal tolerance an urgent topic. The aim of this study was to determine the upper thermal limits of species that are common in tidal pools of rocky shore ecosystems of the Northeast Atlantic. The method used was the Critical Thermal Maximum (CTMax), which allowed the ranking of species in terms of their upper thermal limits as follows: Coryphoblennius galerita (32.0 °C), Palaemon serratus (33.0 °C), Gobius paganellus (33.1 °C), Palaemon elegans (33.4 °C), Lipophrys pholis (33.9 °C) and Paralipophrys trigloides (35.0 °C). Intraspecific variability was always lower than 2%. 相似文献
18.
Large eddy simulations of the flow around a circular cylinder at high Reynolds numbers are reported. Five Reynolds numbers were chosen, such that the drag crisis was captured. A total of 18 cases were computed to investigate the effect of gridding strategy, turbulence modelling, numerical schemes and domain width on the results. It was found that unstructured grids provide better resolution of key flow features, when a ‘reasonable’ grid size is to be maintained.When using coarse grids for large eddy simulation, the effect of turbulence models and numerical schemes becomes more pronounced. The dynamic mixed Smagorinsky model was found to be superior to the Smagorinsky model, since the model coefficient is allowed to dynamically adjust based on the local flow and grid size. A blended upwind-central convection scheme was also found to provide the best accuracy, since a fully central scheme exhibits artificial wiggles, due to dispersion errors, which pollute the solution.Mean drag, fluctuating lift Strouhal number and base pressure are compared to experiments and empirical estimates for Reynolds numbers ranging from 6.31 × 104 to 5.06 × 105. In terms of the drag coefficient, the drag crisis is well captured by the present simulations, although the other integral quantities (rms lift and Strouhal number) show larger discrepancies. For the lowest Reynolds number, the drag is seen to be more sensitive to the domain width than the spanwise grid spacing, while at the higher Reynolds numbers the grid resolution plays a more important role, due to the larger extent of the turbulent boundary layer. 相似文献
19.
《Estuarine, Coastal and Shelf Science》1997,44(3):377-383
The numbers of benthic foraminiferans at four sites in the Clyde Sea area showed no consistent temporal variation throughout 1993. In the finest surface sediments, numbers ranged between 200 and 400 cells cm−3, compared to only 25–50 cells cm−3in the coarsest sediments. On two occasions, high populations of cells less than 63 μm were found in the surface layers. These were thought to represent recruitment peaks since these ‘ juvenile ’ cells grew rapidly when maintained in the laboratory. A total of 56 taxa were identified from the region, the greatest diversity being recorded in the finest sediments. Rose Bengal stained foraminiferans (i.e. presumed living) were found below the anoxic–oxic boundary. The fate of these cells was considered by examining their ability to migrate through fine sediments, and their capacity to survive (based on evidence of pseudopodial activity) periods of anoxia. This study has highlighted the numerical importance of foraminiferans, particularly in fine surface coastal sediments, but questions whether the high populations of ‘ stained ’ cells found in deeper sediments play a significant ecological role. 相似文献
20.
《Deep Sea Research Part I: Oceanographic Research Papers》2007,54(8):1390-1401
Heat transports estimated CTD data collected during the World Ocean Circulation Experiment (WOCE) along the January 1993 30°S hydrographic transect (A10) and the output from a numerical model show a mean heat transport of 0.40 and 0.55±0.24 PW (standard deviation), respectively. The model shows a large annual cycle in heat transport (more than 30% of the variance) with a maximum (minimum) heat transport in July (February) of 0.68 (0.41) PW. Using these data, a method is proposed and evaluated to calculate the heat transport from temperature data obtained from a trans-basin section of expendable bathythermographs (XBTs) profiles. In this method, salinity is estimated from Argo profiles and CTD casts for each XBT temperature observation using statistical relationships between temperature, latitude, longitude and salinity computed along constant-depth surfaces. Full-depth temperature/salinity profiles are obtained by extending the profiles to the bottom of the ocean using deep climatological data. The meridional transport is then determined by using the standard geostrophic method, applying NCEP-derived Ekman transports, and requiring that the salt flux through the Bering Straits be conserved. The results indicate that the methods described here can provide heat transport estimates with a maximum uncertainty of ±0.18 PW (1 PW=1015 W). Most of this uncertainty is due to the climatology used to estimate the deep structure and issues related to not knowing the absolute velocity field and most especially characterizing barotropic motions. Nevertheless, when the methodology is applied to temperatures collected along 30°S (A10) and direct model integrations, the results are very promising. Results from the numerical model suggest that ageostrophic non-Ekman motions can contribute less than 0.05 PW to heat transport estimates in the South Atlantic. 相似文献