不同对流参数化方案试验中凝结加热的特征及对暴雨中尺度模拟结果的影响

使用20 km分辨率的MM5模式,分别选用KUO,GRELL,KAIN-FRITSCH和BETTS-MILLER(以下简称KU,GR,KF和BM)等4种不同对流参数化方案,对1996年8月3~4日石家庄特大暴雨过程作数值模拟试验,分析比较了4个不同试验中网格尺度(显式方案)和次网格尺度(对流参数化方案)凝结加热的水平、垂直分布和时变特征;研究探讨了凝结加热分布及差异对暴雨中尺度模拟结果的可能影响。分析显示,暴雨过程中,4个不同对流参数化方案试验所得到的次网格尺度凝结加热基本都呈单峰特征、加热峰值在对流层中层,但加热层厚度和强度在不同试验间存在差别;4个试验的网格尺度凝结加热的垂直范围表现出较好的一致性,加热重心位于对流层低层,但加热强度仍有所不同;GR和KF及BM试验的总凝结加热率的垂直分布特征主要受其网格尺度凝结加热率特征的影响、加热重心在对流层低层,而KU试验的总凝结加热率的垂直分布特征由其次网格尺度凝结加热率特征所决定、加热重心在对流层中层。研究表明,尽管4个试验在暴雨期间总凝结加热的垂直分布差异并不显著,但对暴雨中尺度模拟的影响却不能忽视。凝结加热的分布特征及演变直接影响与暴雨发生发展密切关联的物理量场的中尺度结构和演变;凝结加热对暴雨中尺度的影响具有连锁性,由加热差异波及局部环流细致结构和强度及其变化的差异,进而影响暴雨发生发展的细致特征。在20 km或更高一些分辨率的条件下,对于描述温带/中纬度暴雨的发展和结构,选用KF方案得到的模拟结果可能更具物理合理性;而KU方案模拟结果容易出现格点气柱的水汽和温度被过量调整的不合理情况。要得到一个可信的中尺度模拟结果,对降水模拟结果进行细化特征的验证、特别是随时间演变特征的验证分析是非常重要的,因为降水的细致演变特征与凝结加热及与之相联系的物理量场的中尺度演变特征密切关联。     

太平洋SSTA对中国东部夏季降水的影响II--数值模拟

基于观测分析的结果,采用NCARCCM3模式,设计4组7个数值试验,研究太平洋海表温度异常对中国东部夏季降水影响及其可能的物理过程。结果表明：数值模拟与观测分析所得的结果一致,前期冬季西北太平洋黑潮海区海温异常,引起亚洲中高纬和东亚东部地区大气环流异常,导致长江流域夏季降水异常;前期春季赤道中东太平洋海温异常,西太平洋副热带高压异常;同期夏季北太平洋中纬中太平洋海区海温异常,激发夏季EUP遥相关型,影响东部夏季降水。     

Experimental and numerical investigations on the intact and damage survivability of a floating–moored oscillating water column device

This paper investigates the intact and damage survivability of a floating–moored Oscillating Water Column (OWC) device using physical model experiments and Computational Fluid Dynamics (CFD) simulations. Different extreme wave conditions have been tested using irregular and regular wave conditions. The device was moored to the tank floor via four vertical taut lines and the effect of the mooring line pre–tension on the device response was studied. It was found that the instantaneous position of the floating device was a key factor in the survivability analysis such that a certain irregular wave train that might not include the largest wave could induce the maximum response. Reducing the pre–tension minimized the maximum surge, but significantly increased the maximum tension due to mooring slack events causing snatch loads. A design regular wave with a period equal to the peak period and a height of 1.9–2.0 times the significant wave height could reasonably predict the same maximum line tension as the irregular sea state, but a smaller wave height was required to achieve the maximum surge. A single failure in the mooring system increased the maximum tension by 1.55 times the intact tension. For a damaged mooring system, using the same design regular wave condition derived from the survivability analysis with an intact mooring system could result in overestimating the maximum tension by more than 20% in comparison to the tension from the irregular sea state, but a smaller regular wave height or a different regular wave condition representing another sea state could lead to the same maximum tension. This highlighted the importance of investigating survival conditions with a damaged mooring system instead of simply using the same conditions derived for the intact mooring system.     

Investigation of drag crisis phenomenon using CFD methods

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 × 10⁴ and 7.57 × 10⁵. 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.     

The effects of blade twist and nacelle shape on the performance of horizontal axis tidal current turbines

The paper presents the effects of blade twist and nacelle shape on the performance of horizontal axis tidal current turbines using both analytical and numerical methods. Firstly, in the hydrodynamic design procedure, the optimal profiles of untwisted and twisted blades and their predicted theoretical turbine performance are obtained using the genetic algorithm method. Although both blade profiles produce desired rated rotational speed, the twisted blade achieves higher power and thrust performance. Secondly, numerical simulation is performed using sliding mesh technique to mimic rotating turbine in ANSYS FLUENT to validate the analytical results. The Reynolds-Averaged Navier-Stokes (RANS) approximation of the turbulence parameters is applied to obtain the flow field around the turbine. It is found that power and axial thrust force from BEMT (Blade Element Momentum Theory) method are under-predicted by 2% and 8% respectively, compared with numerical results. Afterwards, the downstream wake field of the turbine is investigated with two different nacelle shapes. It is found that the rotor performance is not significantly affected by the different nacelle shapes. However, the structural turbulence caused by the conventional nacelle is stronger than that by the NACA-profiled shape, and the former can cause detrimental effect on the performance of the downstream turbines in tidal farms.     

Experimental and numerical study on vortex-induced motions of a deep-draft semi-submersible

An experimental study on vortex-induced motions (VIM) of a deep-draft semi-submersible (DDS) was carried out in a towing tank, with the aim to investigate the VIM effects on the overall hydrodynamics of the structure. In order to study the fluid physics associated with VIM of the DDS, a comprehensive numerical simulation was conducted to examine the characteristics of vortex shedding processes and their interactions due to multiple cylindrical columns. The experimental measurements were obtained for horizontal plane motions including transverse, in-line and yaw motions as well as drag and lift forces on the structure. Spectral analysis was further carried out based on the recorded force time history. These data were subsequently used to validate the numerical model. Detailed numerical results on the vortex flow characteristics revealed that during the “lock-in”, the vortex shedding processes of the upstream columns enhance the vortex shedding processes of the downstream columns leading to the rapid increase of the magnitude of VIM. In addition to the experimental measurements, for the two uniform flow incidences (0° and 45°) investigated, comprehensive numerical data of the parametric study on the VIM characteristics at a wide range of current strength will also serve as quality benchmarks for future study and provide guidance for practical design.     

Numerical investigation on the suppression of VIV for a circular cylinder by three small control rods

Laminar flow past a circular cylinder with 3 small control rods is investigated by numerical simulation. This study is concerned with the suppression efficacy of vortex induced vibration by small control rods located around a main cylinder. The effects of the attack angle and rod-to-cylinder gap ratio on the hydrodynamics and vibration responses of the main cylinder are investigated. The attack angle of α = 45° is performed as the critical angle for VIV suppression of 3 control rods. The 3 control rods have no effect on VIV suppression when the attack angle is less than the critical angle. The 3 control rods have an excellent VIV suppression efficacy when the attack angle is larger than the critical angle. The transverse vibration frequency of the cylinder with 3 control rods is less than that for an isolated cylinder for all the configurations. The numerical results for the configurations of α = 45° & 60°, G/D = 0.6–1.2 show excellent suppression efficient among the cases investigated in this study. The best suppression efficient is found at α = 45°, G/D = 0.9 for 3 control rods. 2 rods in behind of the main cylinder perform more efficient than that of 1 rod in front for VIV suppression as the gap ratio of G/D less than 1.0.     

Numerical simulations of a highway bridge structure employing passive negative stiffness device for seismic protection

A new passive seismic response control device has been developed, fabricated, and tested by the authors and shown to be capable of producing negative stiffness via a purely mechanical mechanism, thus representing a new generation of seismic protection devices. Although the concept of negative stiffness may appear to be a reversal on the desired relationship between the force and displacement in structures (the desired relationship being that the product of restoring force and displacement is nonnegative), when implemented in parallel with a structure having positive stiffness, the combined system appears to have substantially reduced stiffness while remaining stable. Thus, there is an ‘apparent weakening and softening’ of the structure that results in reduced forces and increased displacements (where the weakening and softening is of a non‐damaging nature in that it occurs in a seismic protection device rather than within the structural framing system). Any excessive displacement response can then be limited by incorporating a damping device in parallel with the negative stiffness device. The combination of negative stiffness and passive damping provides a large degree of control over the expected performance of the structure. In this paper, a numerical study is presented on the performance of a seismically isolated highway bridge model that is subjected to various strong earthquake ground motions. The Negative Stiffness Devices (NSDs) are described along with their hysteretic behavior as obtained from a series of cyclic tests wherein the tests were conducted using a modified design of the NSDs (modified for testing within the bridge model). Using the results from the cyclic tests, numerical simulations of the seismic response of the isolated bridge model were conducted for various configurations (with/without negative stiffness devices and/or viscous dampers). The results demonstrate that the addition of negative stiffness devices reduces the base shear substantially, while the deck displacement is limited to acceptable values. This assessment was conducted as part of a NEES (Network for Earthquake Engineering Simulation) project which included shaking table tests of a quarter‐scale highway bridge model. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical modelling of aeolian erosion over a surface with non‐uniformly distributed roughness elements

The present study is focused on the analysis of the mean wall friction velocity on a surface including roughness elements exposed to a turbulent boundary layer. These roughness elements represent non‐erodible particles over an erodible surface of an agglomeration of granular material on industrial sites. A first study has proposed a formulation that describes the evolution of the friction velocity as a function of geometrical parameters and cover rate with different uniform roughness distributions. The present simulations deal with non‐uniform distributions of particles with a random sampling of diameters, heights, positions and arrangements. The evolution (relative to geometrical parameters of the roughness elements) of the friction velocity for several non‐uniform distributions of roughness elements was analysed by the equation proposed in the literature and compared to the results obtained with the numerical simulations. This comparison showed very good agreement. Thus, the formulation developed for uniform particles was found also to be valid for a larger spectrum of particles noted on industrial sites. The present work aims also to investigate in detail the fluid mechanics over several roughness particles. Copyright © 2013 John Wiley & Sons, Ltd.     

3D Modeling of sediment movement by ships-generated wakes in confined shipping channel

Ship-generated waves and return currents are capable of re-suspending significant quantities of bottom and bank sediments.However,most of the previous studies done on the subject do not show how and where sediment is re-suspended by the wakes and the directions of net transport.In this paper,a 3D numerical model based on hydro-sedimentary coupling is presented to search the relationship between the sediment movement,and the pattern of ship-generated waves around and far away from the vessel and the return currents around the ships.The hydrodynamic model is based on 3D Navier-Stokes equations including the standard k-ε model for turbulence processes,and the sediment transport model is based on a 3D equation for the re-suspended sediment transport.The computation results show that the areas of sediment concentration and transport(whether by resuspension or by the bedload) depend mainly on the position,the speed of the ship in the waterways,the kinematics of ship-generated waves and on the return flows.Thus,a map of sediment distribution and the modes of sediment transport generated by the passage of the ship are presented.