潜堤对破碎区至冲泻区水动力特性影响数值分析

破碎区至冲泻区水动力特性分析是研究近岸地区泥沙输运机理和岸滩演变的关键。考虑潜堤修建的影响,对破碎区至冲泻区水动力特性开展研究。运用基于RANS方程的波浪数学模型,选取非线性涡粘性κ-ε紊动传输模型,采用高精度PLIC-VOF方法追踪自由面,并用实验数据验证模型精度。计算分析结果表明,潜堤的修建改变了破波点的位置,破碎区至冲泻区内水动力特性随之变化,入射波波高H、波长L、堤顶水深R和堤顶宽度B是主要影响因素。在相同的入射波条件下,B/R增大,破碎区至冲泻区内出现的最大紊动动能值和最大紊动耗散值减小。用量纲一参数Re,Fr和St来描述破碎区至冲泻区水动力特性,随着RL/BH增大,波浪将在离岸较近处发生破碎,Re和Fr增大,St减小。     

长江口斜压诊断模式三维流场数值模拟

建立了σ坐标系下长江口斜压诊断模式三维流场数学模型,采用k-kl二方程紊流闭合模型求解垂向涡粘系数,计算域内恒定、非均匀盐度场反映了计算域密度斜压效应。验证结果表明,斜压诊断模式的模拟成果较切合实际。模型成功复演了长江口水域由往复流向旋转流过渡的流场特性;潮波从四条汊道传入后在分汊口附近相遇形成多个会潮点,其位置随着径流的大小、潮汐的强弱、汊道的形态等因素而变化。     

基于非静压模型的波浪破碎模拟

波浪破碎是海岸工程所关注的关键水动力学问题之一,而波浪破碎的数值模拟技术的研究方兴未艾。高效的浅水方程基于静压假定,而通过引入动压项建立的完全非静压模型,可成功应用于色散水波的模拟。自由表面的捕捉采用的垂向坐标变换技术,较之VOF（Volume of Fluid）模型,计算效率较高。但垂向坐标变换不能模拟大曲率自由表面变形,即波浪破碎过程。对于破碎波的模拟,一种高效的模式分裂法应用至非静压模型,即在波浪破碎局部水域将模型分裂为静压模型和非静压模型,破碎波波峰附近退化至静压模型,并持续至波浪破碎结束再恢复为动压模型。通过典型算例,验证了模式分裂法的适用性及对于波浪破碎过程的模拟精度,鉴于模式分裂法对于波浪破碎过程的模拟未引入新的概化模型,计算效率较高,可应用于大尺度的海岸带波浪的变形、破碎及越浪的数值模拟。     

Modelling flash flood propagation in urban areas using a two-dimensional numerical model

This paper reports on the numerical modelling of flash flood propagation in urban areas after an excessive rainfall event or dam/dyke break wave. A two-dimensional (2-D) depth-averaged shallow-water model is used, with a refined grid of quadrilaterals and triangles for representing the urban area topography. The 2-D shallow-water equations are solved using the explicit second-order scheme that is adapted from MUSCL approach. Four applications are described to demonstrate the potential benefits and limits of 2-D modelling: (i) laboratory experimental dam-break wave in the presence of an isolated building; (ii) flash flood over a physical model of the urbanized Toce river valley in Italy; (iii) flash flood in October 1988 at the city of Nîmes (France) and (iv) dam-break flood in October 1982 at the town of Sumacárcel (Spain). Computed flow depths and velocities compare well with recorded data, although for the experimental study on dam-break wave some discrepancies are observed around buildings, where the flow is strongly 3-D in character. The numerical simulations show that the flow depths and flood wave celerity are significantly affected by the presence of buildings in comparison with the original floodplain. Further, this study confirms the importance of topography and roughness coefficient for flood propagation simulation.     

Wave-turbulence instability in nonequilibrium hydrodynamics systems

The problem of interaction between a large-scale wave and the small-scale wave-induced turbulence in non-equilibrium systems of hydrodynamic type is considered. The results of numerical analysis of the model of the turbulence-wave interaction in the hydrodynamic medium (an interior gravity wave + turbulence in a stratified fluid with a shear flow) and in plasma (an ion-sound wave + the Lengmuir turbulence) are presented. It is shown that significant mutual amplification of the wave and turbulence is possible at some values of the parameters.     

斜坡上异重流的三维数值模拟

针对异重流的流动特征,建立了适用于具有各向异性浮力紊动特征的三维异重流运动的数学模型,并模拟了异重流在15°斜坡底面上的潜行过程。计算结果准确地模拟了异重流的运动特征和形态,其前锋的潜行速度与实验结果相当吻合。该模型采用非结构同位网格上的SIMPLEC算法能适应复杂边界和地形,可应用于自然界实际环境中异重流的演进计算。     

应用二阶完全非线性Boussinesq方程模拟破碎波浪

建立了基于高阶Boussinesq水波方程的一维波浪破碎数值模型。基于一组具有二阶完全非线性特征的Boussinesq水波方程,建立了交错网格下的高精度差分格式,推导了适用于该组方程的永形波解析解,其和松弛造波技术相结合实现了数值波浪水槽中(强)非线性波浪的无反射入射。通过模拟封闭容器内水体晃动问题对数值格式进行了验证,通过模拟孤立波在斜坡海岸上的浅化过程说明了将方程从弱非线性发展到完全非线性的必要性。采用涡粘方法处理波浪破碎,利用物理模型实验数据,分析了模型中各波浪破碎参数对数值结果的影响并对参数进行了率定。应用该模型对规则波在斜坡海岸上的传播、变浅以及破碎过程进行了数值模拟研究,数值结果同实验数据吻合良好,验证了模型的有效性。     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

Simulated annealing in 3-D seismic modelling and elastic structure of the mediterranean basin

3-D elastic structures images of the Mediterranean basin are obtained by solving a tomographic problem by means of a simulated annealing algorithm. The database consists of local Rayleigh wave group velocities determined for the fundamental mode in the 10–90s period range. These local velocities were obtained from surface wavetrains recorded at very broad band stations of MEDNET network. From a methodological point of view, the characteristics of the annealing algorithm and the amount of time required for computation suggest the implementation of parallel computational procedures to shorten the computational time. According to the thermodynamic analogies used in the annealing algorithm, we propose an alternative procedure to estimate uncertainties and correlations for variables defining the elastic model (shear velocities and layer thicknesses) without necessity of computing partial derivatives. From an applied point of view, the 3-D elastic model obtained shows, as main results, relevant differences for the eastern and western pans of the Mediterranean basin, specially down to 60 km. The changes in the position of the top and bottom asthenosphere surfaces can also be followed throughout the Mediterranean basin; remarkably, these changes are in agreement with some tectonic features.     

波流环境中垂向圆管射流三维运动和稀释过程模拟

采用大涡数值模拟方法建立了波流环境中垂向圆管射流三维物质输运数学模型,在此基础上分析了波流共同作用下射流运动和稀释过程,并探讨了波浪对射流三维时均浓度特征的影响。数模结果与实测数据的对比表明,该模型能够正确复演波流环境中射流速度和浓度沿水深分布的规律。引入了射流的三维浓度特征指标,即断面最小稀释度及其垂向位置和断面可视范围面积,发现随着波高或波周期的增大,波流环境中射流的断面可视范围面积逐渐增加,这说明波浪对射流的稀释起到一定的促进作用。     

非结构网格上的三维浅水流动数值模型

针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。     

洞塞式消能工的数值模拟

针对有压洞内洞塞式消能工附近局部流动,分析研究了水流紊动特点,应用三维RNG k-ε紊流模型,利用物理模型试验资料对其计算结果进行了验证,利用模型对洞塞式消能工引起的水头损失进行计算分析,分析结果表明,当洞塞段相对长度从很小开始增加时,水头损失系数呈现先急剧降低而后明显增加,最后缓慢线性增加的趋势;在有压洞半径一定的情况下,洞塞段过流半径的变化比洞塞段过流长度的变化对水头损失系数影响更为明显,可为洞塞式消能工的设计提供相关参考依据。     

基于非结构混合网格求解二维浅水方程的一种数值方法

复杂边界下的流场数值模拟常基于非结构化网格进行求解,建立一种在非结构混合网格上求解水深平均的二维浅水方程的模型,以便精确模拟复杂边界、提高计算效率。该模型时间项离散采用隐格式使得模型具有较好的稳定性,对流项和扩散项分别采用总变差减小(Total Variation Diminishing,TVD)格式和构造辅助点的方法来离散,同时采用水深平均的标准k-ε模型来封闭湍流模型。选用两个经典验证算例检验模型,计算结果表明,基于非结构混合网格开发的模型具有较高的精度,且收敛性能较好。     

植物消波机制的实验与理论解析研究进展

从以物理模型实验及理论解析为代表的植物消波研究手段与宏观和微观尺度下植物消波分析方法两个方面总结了国内外的研究现状和最新研究进展。对物理模型实验中的植物因子、水动力因子、植物模型制作和实验布置等要素以及理论解析模型进行归纳,分别阐述了以波高衰减和波能耗散为主要分析对象的植物消波分析方法,指出现有物理模型实验中植物几何比尺与波要素比尺难以统一和理论解析研究局限于二维情况等问题。综合考虑季节更替、植物摆动、波浪破碎、波浪随机性及非线性、潮流等因素是今后植物消波研究的重点突破方向。     

远破波作用数值模拟的CLEAR-VOF模型

在非结构化网格基础上,采用三阶精度的三步有限元方法离散N-S方程,CLEAR-VOF方法追踪运动水体自由表面,建立了数值波浪水槽模型.模拟了远破波冲击位于倾斜底坡上直立堤的过程.对直立堤上波压力和堤前流场的数值结果与模型实验结果做了比较.结果表明,远破波波浪力和流场的数值结果同PIV实验结果吻合较好,该数值波浪水槽模型可较好地模拟远破波的破碎冲击过程.     

Hybrid differential evolution and particle swarm optimization for optimal well placement

There is no gainsaying that determining the optimal number, type, and location of hydrocarbon reservoir wells is a very important aspect of field development planning. The reason behind this fact is not farfetched—the objective of any field development exercise is to maximize the total hydrocarbon recovery, which for all intents and purposes, can be measured by an economic criterion such as the net present value of the reservoir during its estimated operational life-cycle. Since the cost of drilling and completion of wells can be significantly high (millions of dollars), there is need for some form of operational and economic justification of potential well configuration, so that the ultimate purpose of maximizing production and asset value is not defeated in the long run. The problem, however, is that well optimization problems are by no means trivial. Inherent drawbacks include the associated computational cost of evaluating the objective function, the high dimensionality of the search space, and the effects of a continuous range of geological uncertainty. In this paper, the differential evolution (DE) and the particle swarm optimization (PSO) algorithms are applied to well placement problems. The results emanating from both algorithms are compared with results obtained by applying a third algorithm called hybrid particle swarm differential evolution (HPSDE)—a product of the hybridization of DE and PSO algorithms. Three cases involving the placement of vertical wells in 2-D and 3-D reservoir models are considered. In two of the three cases, a max-mean objective robust optimization was performed to address geological uncertainty arising from the mismatch between real physical reservoir and the reservoir model. We demonstrate that the performance of DE and PSO algorithms is dependent on the total number of function evaluations performed; importantly, we show that in all cases, HPSDE algorithm outperforms both DE and PSO algorithms. Based on the evidence of these findings, we hold the view that hybridized metaheuristic optimization algorithms (such as HPSDE) are applicable in this problem domain and could be potentially useful in other reservoir engineering problems.     

一种改进的近岸波浪破碎数值模型

提出一种基于完全非线性Boussinesq方程的改进波浪破碎模型,用于模拟近岸浅水波浪破碎.模型借助水滚的概念,提出一种确定k方程紊动模式中紊动源项的计算公式,通过解k方程计算Boussinesq方程中的涡黏系数,实现对破碎波的模拟.岸边界采用窄缝法,使得模型可用于波浪爬坡的计算.用实验室实测波高和增减水资料对模型进行了验证,得到了一致的结果.紊动源、紊动动能以及紊动耗散率的计算结果表明:①在破波点处紊动源项值最大,随着波浪向岸边传播,逐渐减小;②破波点处,水平方向的对流和扩散在紊动能量输移中发挥重要作用;③岸边附近紊动源与紊动耗散接近平衡.     

A numerical study on 2-D compressive wave propagation in rock masses with a set of joints along the radial direction normal to the joints

An explosion in a borehole or an accidental explosion in a tunnel will generate a two-dimensional (2-D) compressive wave that travels through the surrounding rock mass. For the problem of 2-D compressive wave propagation in a rock mass with parallel joints in the radial direction normal to the joints, parametric studies on the transmission ratio and the maximum rebound ratio are performed in universal distinct element code. The variation of the transmission ratio with the ratio of joint spacing to wavelength is generalized into a characteristic curve, which can be used as a prediction model for estimating the transmission ratio. The relationship between the maximum rebound ratio and the influence factors is obtained in charts. The charts can be used as a prediction model for estimating the maximum rebound ratio. The proposed prediction models for estimating the transmission ratio and the maximum rebound ratio are applied to a field explosion test, Mandai test in Singapore. The minimum possible values of peak particle velocity (PPV) at the monitoring points are estimated by using the proposed prediction model for estimating the transmission ratio along the radial direction normal to the joints. On the other hand, the maximum possible values of PPV are estimated by using the proposed prediction model for estimating the maximum rebound ratio along the same radial direction. The comparison shows a good agreement between the field-recorded PPVs and the estimated range of PPVs given by the minimum possible PPVs and the maximum possible PPVs at the monitoring points in Mandai test. The good agreement between the estimated and field-recorded values validates the proposed prediction models for estimating peak particle velocity in a rock mass with a set of joints due to application of a compressive wave at the boundary of a tunnel or a borehole.     

基于离散颗粒模型的岩体结构面动-静态刚度系数对比的数值模拟研究

岩体工程计算分析中结构面刚度系数是至关重要的力学参数,计算分析的精度和可靠程度与这个参数密不可分,结构面刚度系数取值仍然是一个难点.岩体中应力波传播至结构面处将会发生反射和透射现象,利用应力波透射系数可反演结构面动态刚度系数.本文从细观力学角度运用颗粒离散元方法,开发分段线性接触模型及应力波吸收边界模型,开展宏观岩体中...     

Mesh size influence on dynamic modeling

Mesh size influence on dynamic modelling in 1-D and 2-D models is investigated by using discrete element code UDEC. It is shown that the modelling with smaller mesh ratio obtains higher accuracy, which agrees well with the common acceptance. Furthermore, it is found that the error, induced by the mesh size, increases with increasing distance from the wave source and also with increasing wave frequency. Therefore, it is suggested that the mesh distortion may resulted by a combination effect of mesh size, wave frequency and distance from the wave source. Comparing the results between 1-D and 2-D models indicates that the 2-D modelling produces higher error than 1-D model, which is explained as due to the error superposition in x and y directions in the 2-D model.