同位网格摄动有限体积格式求解浮力驱动方腔流

利用对流扩散方程的摄动有限体积格式，在Rayleigh数从10$^{3}$ 到10$^{8}$的范围内对浮力驱动方腔流动问题作了数值模拟. 对流扩散方程的摄动 有限体积格式具有一阶迎风格式的简洁形式，使用相同的基点，重构近似精度高，特别是两 相邻控制体中心到公共界面的距离相等或不相等，PFV格式公式相同等优点. 在数值模拟中， 无论均匀网格还是非均匀网格均获得与DSC方法、自适应有限元法、多重网格法等Benchmark 解相符较好的数值结果，证明UPFV格式对高Rayleigh数对流传热问题的适用性和有效性.     

关于交通流中扰动传播和发展的数值模拟

针对不同拥挤程度交通流中扰动的传播和发展，进行了系统的数值模拟研究. 讨论了PW模型、一维管流模型、单侧传播模型等几种典型交通流模型的不同特征；指出了 不同差分格式对数值结果的影响；也分析了数学模型中不同形式的平衡函数的作用. 结果发 现在适当的数学模型、平衡函数和离散格式下，能够对交通流中扰动的传播和发展，特别是 扰动波的传播速度，取得与实际测量数据相当接近的模拟效果. 建议了分别适用于不同密度 交通流的计算格式.     

随机波浪作用下近岸波流场的数值模拟

结合近岸波浪抛物型缓坡模型和近岸波流场模型,对近岸不规则波浪及其破碎后所产生的流场进行了数值模拟. 在不规则波浪场的模拟中,采用JONSWAP波浪谱对入射单向不规则波浪要素按等分频率法进行离散,应用考虑波浪不规则性和破碎效应的抛物型缓坡方程对波浪场进行数值模拟,并基于抛物型缓坡方程中的波浪势函数等参数计算波浪辐射应力,以波浪辐射应力为主要动力因素基于近岸流数学模型对近岸波浪破碎所产生的近岸流场进行数值模拟,并对数值模拟结果进行了验证. 模拟结果表明该模型可有效地用于研究波浪破碎产生的近岸波流场.      

基于VPM-THINC/QQ模型的波浪高保真模拟

为实现波浪传播的高保真数值模拟,采用包含单元均值和点值(volume-average/point-value method,VPM)的有限体积法求解纳维-斯托克斯方程和具有二次曲面性质和高斯积分的双曲正切函数(THINC method with quadratic surface representation and Gaussian quadrature,THINC/QQ)方法来重构自由面,建立以开源求解库OpenFOAM底层函数库为基础的VPM-THINC/QQ模型. 在本模型中添加推板造波法实现波浪的产生功能,采用松弛法实现消波功能,构建高精度黏性流数值波浪水槽. 分别采用VPM-THINC/QQ模型和InterFoam求解器(OpenFOAM软件包中广泛使用的多相流求解器)开展规则波的数值模拟,重点探究网格大小和时间步长等因素对波浪传播过程的影响,定量地分析波高衰减程度;为验证本模型的适应性,对长短波进行模拟. 结果表明,在相同网格大小或时间步长条件下,VPM-THINC/QQ模型的预测结果与参考值吻合较好,波高衰减较少,且无相位差,在波浪传播过程的模拟中呈现出良好的保真性. 本文工作 为波浪传播的模拟研究提供了一种高精度的黏性数值波浪水槽模型.      

基于VPM-THINC/QQ模型的波浪高保真模拟

为实现波浪传播的高保真数值模拟,采用包含单元均值和点值(volume-average/point-value method,VPM)的有限体积法求解纳维-斯托克斯方程和具有二次曲面性质和高斯积分的双曲正切函数(THINC method with quadratic surface representation and Gaussian quadrature,THINC/QQ)方法来重构自由面,建立以开源求解库OpenFOAM底层函数库为基础的VPM-THINC/QQ模型.在本模型中添加推板造波法实现波浪的产生功能,采用松弛法实现消波功能,构建高精度黏性流数值波浪水槽.分别采用VPM-THINC/QQ模型和InterFoam求解器(OpenFOAM软件包中广泛使用的多相流求解器)开展规则波的数值模拟,重点探究网格大小和时间步长等因素对波浪传播过程的影响,定量地分析波高衰减程度;为验证本模型的适应性,对长短波进行模拟.结果表明,在相同网格大小或时间步长条件下,VPM-THINC/QQ模型的预测结果与参考值吻合较好,波高衰减较少,且无相位差,在波浪传播过程的模拟中呈现出良好的保真性.本文工作为波浪传播的模拟研究提供了一种高精度的黏性数值波浪水槽模型.     

均匀流中双色波传播特性的模拟研究

针对双色波浪与均匀流相互作用问题,采用时域高阶边界元方法建立自由水面满足完全非线性边界条件的数学模型。求解中采用混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过与已发表试验结果对比,验证了本模型的准确性。通过数值计算研究了水流参数对各组成波及衍生的高阶波幅值、波浪和水流间能量交换的影响规律。     

非结构化网格下大范围波浪的折绕射计算

从含流缓坡方程出发,推导出光程函数方程。采用基于非结构化网格下的有限体积法对光程函数方程和波作用守恒方程进行数值离散和联合求解,从而构建了一个考虑绕射的计算近岸大范围波浪传播过程的数值模型。模型的空间步长不再受制于波长的限制,同时非结构化网格可以很好地拟合复杂岸线变形。模型分别通过了施奈尔定律、直立防波堤后的波浪绕射和圆形浅滩上的波浪变形验证,结果表明,该数值模式能有效模拟复杂边界条件和大范围水域下近岸波浪的传播过程中的折射和绕射等变形。     

波浪场中物质输移离散系数的理论分析

海岸水域波浪引起的物质输移扩散存在着不同于水流引起的物质输移扩散的特征. 通过解析的方法对垂向扩散方程进行求解, 推导出波浪和潮流共同作用下的海岸水域物质输移离散系数的理论公式. 分析中将总的离散系数分为潮流、质量输移流、纯波浪波动作用以及潮流与质量输移流、潮流与波浪波动、质量输移流与波浪波动的相互作用6部分, 通过分析各组成部分特征, 得出海岸水域物质输移中波浪所产生的离散效应主要是质量输移流贡献, 波浪水质点运动及波流相互作用对时间平均离散系数的贡献较小的结论.仅波浪的作用(不包括潮流)的有关结果与他人数值模拟和实验结果的吻合良好, 证明了该理论推导的正确性. 给出了时间平均离散系数及离散系数波动幅值随波浪周期和波高的变化特征,同时将结果应用于垂向环流, 给出了垂向环流对应于质量输移流所产生的离散系数的特征.      

弯道水流二维浅水模型的修正研究

借鉴有关弯道水流流速分布的研究成果,计入深度平均流速与真实流流速分布差值引起的扩散效应,在正交曲线坐标系下建立了平面二维浅水模型.采用以标准κ-ε模型为基础的曲率效应修正紊流模型模拟紊动应力项,在一定程度上考虑了流线弯曲水流紊动应力的各向异性.应用控制体积法和交错网格法离散方程,并用SIMPLEC算法求解离散方程;同时采用修正后的模型对90°弯道水流进行了数值模拟,并与原模型的计算结果及实测资料进行了比较,结果表明该模型能够有效地模拟流线弯曲水流的水力特性.     

模拟畸形波的聚焦波浪模型

利用改进的高阶谱方法建立了模拟极限波的二维聚焦模型,通过与Baldock(1996)的实验结果和理论值的比较,验证了模型的正确性,并分析了波浪非线性的相互作用对聚焦结果的影响. 通过改进Longuet-Higgins海浪模型,给出了4种实验室聚焦模拟畸形波的波浪模型:极限波聚焦模型+随机波模型;极限波聚焦模型+规则波模型;相位角分布范围调制聚焦模型;相同相位角组成波个数调制聚焦模型. 基于上述完全非线性数值波浪模型,采用不同的能量分配方式,在有限模拟长度和时间内得到了具有不同$H_{\max}/H_{s}$值的畸形波.      

非均匀水流水域波浪的传播变形

将两个不同的、考虑波流相互作用和能量耗散项的、依赖时间变化的双曲型缓坡方程分别化 为一组等价的控制方程组，具体分析了底摩阻项对相对频率和波数矢的影响，从而选择了合 适的数学模型. 将所选择的缓坡方程化为依赖时间变化的抛物型方程，并用ADI法进 行数值求解，建立了缓变水深水域非均匀水流中波浪传播的数值模拟模型. 通过和波流共线 的解析解的比较，说明数值解和解析解相一致. 结合Arthur(1950)水流这一经典算例，定 量地讨论了考虑联合折射-绕射作用后的波数和仅考虑折射作用的波数的差别及其对波高分 布的影响. 在基本同样的条件下, 本文的数值解与他人的计算结果一致.     

连续分层海洋中内波传播的一种数值模式

基于Euler方程,使用有限体积法建立了一种密度为连续分层情况下、适应水深变化的水域中内波传播的数值模式.为了使计算格式能够达到二阶精度,对流项的处理使用了TVD (total variation diminishing)格式.将SIMPLE算法引入连续分层海洋中内波的数值计算,为了简化计算并方便地适应多种TVD格式,在计算预估速度场时采用了显式格式,而没有采用传统的隐式格式;鉴于在原始的SIMPLE算法中没有涉及到由于密度扰动而引起的静水压力场的改变问题,给出了该问题的计算方法.因此改进了SIMPLE算法.出流边界的处理采用阻尼消波和Sommerfeld辐射条件相结合的方式,以使内波得到有效的衰减和释放.将等水深水域的数值解和理论解进行了比较,两者吻合较好;并对存在潜堤时数值计算的不同时刻密度变化的空间分布进行了详细的定性分析.计算结果表明,所建立的数值模式能有效地模拟内波的传播和变形.      

Incompressible SPH simulation of wave breaking and overtopping with turbulence modelling

In this paper a truly incompressible version of the smoothed particle hydrodynamics (SPH) method is presented to investigate the surface wave overtopping. SPH is a pure Lagrangian approach which can handle large deformations of the free surface with high accuracy. The governing equations are solved based on the SPH particle interaction models and the incompressible algorithm of pressure projection is implemented by enforcing the constant particle density. The two‐equation k–ε model is an effective way of dealing with the turbulence and vortices during wave breaking and overtopping and it is coupled with the incompressible SPH numerical scheme. The SPH model is employed to reproduce the experiment and computations of wave overtopping of a sloping sea wall. The computations are validated against the experimental and numerical data found in the literatures and good agreement is observed. Besides, the convergence behaviour of the numerical scheme and the effects of particle spacing refinement and turbulence modelling on the simulation results are also investigated in further detail. The sensitivity of the computed wave breaking and overtopping on these issues is discussed and clarified. Copyright © 2005 John Wiley & Sons, Ltd.     

A horizontally curvilinear non‐hydrostatic model for simulating nonlinear wave motion in curved boundaries

A horizontally curvilinear non‐hydrostatic free surface model that embeds the second‐order projection method, the so‐called θ scheme, in fractional time stepping is developed to simulate nonlinear wave motion in curved boundaries. The model solves the unsteady, Navier–Stokes equations in a three‐dimensional curvilinear domain by incorporating the kinematic free surface boundary condition with a top‐layer boundary condition, which has been developed to improve the numerical accuracy and efficiency of the non‐hydrostatic model in the standard staggered grid layout. The second‐order Adams–Bashforth scheme with the third‐order spatial upwind method is implemented in discretizing advection terms. Numerical accuracy in terms of nonlinear phase speed and amplitude is verified against the nonlinear Stokes wave theory over varying wave steepness in a two‐dimensional numerical wave tank. The model is then applied to investigate the nonlinear wave characteristics in the presence of dispersion caused by reflection and diffraction in a semicircular channel. The model results agree quantitatively with superimposed analytical solutions. Finally, the model is applied to simulate nonlinear wave run‐ups caused by wave‐body interaction around a bottom‐mounted cylinder. The numerical results exhibit good agreement with experimental data and the second‐order diffraction theory. Overall, it is shown that the developed model, with only three vertical layers, is capable of accurately simulating nonlinear waves interacting within curved boundaries. Copyright © 2011 John Wiley & Sons, Ltd.     

A new numerical model for simulations of wave transformation,breaking and long‐shore currents in complex coastal regions

In this paper, we propose a model based on a new contravariant integral form of the fully nonlinear Boussinesq equations in order to simulate wave transformation phenomena, wave breaking, and nearshore currents in computational domains representing the complex morphology of real coastal regions. The aforementioned contravariant integral form, in which Christoffel symbols are absent, is characterized by the fact that the continuity equation does not include any dispersive term. A procedure developed in order to correct errors related to the difficulties of numerically satisfying the metric identities in the numerical integration of fully nonlinear Boussinesq equation on generalized boundary‐conforming grids is presented. The Boussinesq equation system is numerically solved by a hybrid finite volume–finite difference scheme. The proposed high‐order upwind weighted essentially non‐oscillatory finite volume scheme involves an exact Riemann solver and is based on a genuinely two‐dimensional reconstruction procedure, which uses a convex combination of biquadratic polynomials. The wave breaking is represented by discontinuities of the weak solution of the integral form of the nonlinear shallow water equations. The capacity of the proposed model to correctly represent wave propagation, wave breaking, and wave‐induced currents is verified against test cases present in the literature. The results obtained are compared with experimental measures, analytical solutions, or alternative numerical solutions. Copyright © 2015 John Wiley & Sons, Ltd.     

Interaction of a shock wave with a loose dusty bulk layer

The interaction of a planar shock wave with a loose dusty bulk layer has been investigated both experimentally and numerically. Experiments were conducted in a shock tube. The incident shock wave velocity and particle diameters were measured with the use of pressure transducers and a Malvern particle sizer, respectively. The flow fields, induced by shock waves, of both gas and granular phase were visualized by means of shadowgraphs and pulsed X-ray radiography with trace particles added. In addition, a two-phase model for granular flow presented by Gidaspow is introduced and is extended to describe such a complex phenomenon. Based on the kinetic theory, such a two-phase model has the advantage of being able to clarify many physical concepts, like particulate viscosity, granular conductivity and solid pressure, and deduce the correlative constitutive equations of the solid phase. The AUSM scheme was employed for the numerical calculation. The flow field behind the shock wave was displayed numerically and agrees well with our corresponding experimental results.      

基于Boussinesq方程的波浪模型

从欧拉方程出发，提供了另一种推导完全非线性Boussinesq方程的方法，并对方程的 线性色散关系和线性变浅率进行了改进. 改进后方程的线性色散关系达到了一阶Stokes波 色散关系的Pad\'{e}[4,4]近似，在相对水深达1.0的强色散波浪时仍保持较高的准确性，并且方程的非线性和线性 变浅率都得到了不同程度的改善. 方程的水平一维形式用预估-校正的有限差分格式求解， 建立了一个适合较强非线性波浪的Boussinesq波浪数值模型. 作为验证，模拟了波浪在潜 堤上的传播变形，计算结果和实验数据的比较发现两者符合良好.     

Selective lumping finite element method for nearshore current

A finite element method for the analysis of nearshore current, which is one of the principal currents in coastal seas, is presented in this paper. Because the nearshore current is induced by the variable distribution of the surface waves, it is necessary to analyse two main characteristics of the wave, i.e. direction and height. The current can be computed using the resulting wave characteristics. The present method makes it possible to employ procedures for which the same methods of solution are applicable for all basic equations of wave direction, height and current flow. The linear interpolation function is used for the discretization of spatial variables and a selective lumping two step explicit scheme is employed for the numerical integration in time. The numerical solutions obtained are compared with analytical, experimental and observed ones. From these comparative studies, it is concluded that the present finite element method provide a useful tool for the analysis of nearshore current.     

On the analysis of 2D nonlinear gravity waves with a fully nonlinear numerical model

A fully nonlinear numerical method, developed on the basis of Euler equations, is used to study the dynamics of nonlinear gravity waves, mainly in the aspects of the propagation of Stokes wave with disturbed sidebands, the evolution of one wave packet and the interaction of two wave groups. These cases have previously been studied with the higher order spectral method, which will be an approximately fully nonlinear scheme if the order of nonlinearity is not large enough, while the present method in the case of the 2D model has an integration scheme that is exact to the computer precision. As expected, in most cases the results are consistent between these two numerical models and it is confirmed again that this fully nonlinear numerical model is also capable of maintaining a high accuracy and good convergence, particularly in the long-term evolutionary process.