Internal wave generation in an improved two-dimensional Boussinesq model

A set of Boussinesq-type equations with improved linear frequency dispersion in deeper water is solved numerically using a fourth order accurate predictor-corrector method. The model can be used to simulate the evolution of relatively long, weakly nonlinear waves in water of constant or variable depth provided the bed slope is of the same order of magnitude as the frequency dispersion parameter. By performing a linearized stability analysis, the phase and amplitude portraits of the numerical schemes are quantified, providing important information on practical grid resolutions in time and space. In contrast to previous models of the same kind, the incident wave field is generated inside the fluid domain by considering the scattered wave field in one part of the fluid domain and the total wave field in the other. Consequently, waves leaving the fluid domain are absorbed almost perfectly in the boundary regions by employment of damping terms in the mass and momentum equations. Additionally, the form of the incident regular wave field is computed by a Fourier approximation method which satisfies the governing equations accurately in water of constant depth. Since the Fourier approximation method requires an Eulerian mean current below wave trough level or a net mass transport velocity to be specified, the method can be used to study the interaction of waves and currents in closed as well as open basins. Several computational examples are given. These illustrate the potential of the wave generation method and the capability of the developed model.     

粘性泥石流阻力和运动方程验证分析

现有的粘性泥石流阻力方程可以大致归纳为三类：一类是将泥石流视为固液两相流,通过理论分析建立的泥石流阻力运动方程;二类是认为粘性泥石流符合宾汉姆流体阻力方程;三类是依据泥石流体具有基本符合库伦公式的剪切强度,将泥石流视为固体颗粒散体重力流,在理论分析基础上建立的阻力和运动方程。用粘性泥石流动力学实验数据对三类阻力运动方程进行了验证,结果表明第三类阻力和运动方程与实际接近。     

系泊船非线性波浪力时域计算:二维模型

为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算.     

Propagation of a solitary wave over rigid porous beds

The unsteady two-dimensional Navier–Stokes equations and Navier–Stokes type model equations for porous flows were solved numerically to simulate the propagation of a solitary wave over porous beds. The free surface boundary conditions and the interfacial boundary conditions between the water region and the porous bed are in complete form. The incoming waves were generated using a piston type wavemaker set up in the computational domain. Accuracy of the numerical model was verified by comparing the numerical results with the theoretical solutions. The main characteristics of the flow fields in both the water region and the porous bed were discussed by specifying the velocity fields. Behaviors of boundary layer flows in both fluid and porous bed regions were also revealed. Effects of different parameters on the wave height attenuation were studied and discussed. The results of this numerical model indicate that for the investigated incident wave as the ratio of the porous bed depth to the fluid depth exceeds 10, any further increase of the porous bed depth has no effect on wave height attenuation.     

Essential properties of Boussinesq equations for internal and surface waves in a two-fluid system

Boussinesq equations describing motions of internal waves in a two-fluid system with the presence of free surface are theoretically derived, and the associated essential properties are examined in this study. Eliminating the dependence on the vertical coordinate from all variables, four equations constitute the Boussinesq model with two flexible parameters, z_u and z_l, which indicate the specific elevations, respectively, in the upper and lower fluids. Similar to the Boussinesq model for a single-layer fluid, z_u and z_l are determined by matching the linear dispersion relation with Lamb's solution. This determines the optimal model. In the analysis stage, this problem is classified into two cases, the thicker-upper-layer case and the thicker-lower-case case, to avoid the possible divergence of wave properties as the thickness ratio grows. Since there exist two modes of motions that may be excited, cases of both modes are separately analyzed. Linear characteristics including the amplitude ratios and normalized particle velocities are analyzed. Second-order harmonic waves are examined to validate nonlinear behaviors of present model. Results of linear and nonlinear investigations show that the present model indeed extends the applicable range of traditional Boussinesq equations.     

Two-dimensional scour simulations based on coupled model of shallow water equations and sediment transport on unstructured meshes

A two-dimensional scour model based on coupled system of shallow water equations (SWEs) and sediment transport on unstructured mesh is developed. The coupled system of hydrodynamic and morphodynamic equations is solved by finite volume method using Godunov scheme. Roe's approximate Riemann solver is used to calculate the inviscid fluxes. The use of unstructured mesh makes the model applicable to complex domains. However, it is difficult to evaluate the eigenvalues and eigenvectors of the Jacobian matrix in the global coordinate. The method proposed herein to deal with this difficulty is to transform the system into the local coordinate with one of the axes in the same direction as the interface outward normal vector. In the local coordinate system, the Jacobian matrix is simplified and the eigenvalues are analyzed using asymptotic method. Regular expansion breaks down when the flow is near critical. Uniformity of the expansion is achieved by changing the scales. Rotational invariance theorem is used to relate the interfacial fluxes in the global and local coordinate systems. Special treatment of the source term on unstructured grid makes the scheme stable and physically balanced (both mass and momentum). The method proposed in this paper for the eigen-system is very efficient comparing to iterative numerical methods. Results from the test cases show good agreement with the experiments.     

Numerical Simulation on Hydraulic Performances of Quarter Circular Breakwater

Quarter circular breakwater （QCB） is a new-type breakwater developed from senti-circular breakwater （SCB）. The superstructure of QCB is composed of a quarter circular front wall, a horizontal base slab and a vertical rear wall. The width of QCB＇ s base slab is about half that of SCB, which makes QCB suitable to be used on relatively finn soil foundation. The numerical wave flume based on the Reynolds averaged Navier-Stokes equations for impressible viscosity fluid is adopted in this paper to simulate the hydraulic performances of QCB. Since the geometry of both breakwaters is similar and SCB has been studied in depth, the hydraulic performances of QCB are given in comparison with those of SCB.     

A numerical study of swash flows generated by bores

The dynamic processes of bore propagation over a uniform slope are studied numerically using a 2-D Reynolds Averaged Navier–Stokes (RANS) solver, coupled to a non-linear k − ε turbulence closure and a volume of fluid (VOF) method. The dam-break mechanism is used to generate bores in a constant depth region. Present numerical results for the ensemble-averaged flow field are compared with existing experimental data as well as theoretical and numerical results based on non-linear shallow water (NSW) equations. Reasonable agreement between the present numerical solutions and experimental data is observed. Using the numerical results, small-scale bore behaviors and flow features, such as the bore collapse process near the still-water shoreline, the ‘mini-collapse’ during the runup phase and the ‘back-wash bore’ in the down-rush phase, are described. In the case of a strong bore, the evolution of the averaged turbulence kinetic energy (TKE) over the swash zone consists of two phases: in the region near the still-water shoreline, the production and the dissipation of TKE are roughly in balance; in the region farther landwards of the still-water shoreline, the TKE decay rate is very close to that of homogeneous grid turbulence. On the other hand, in the case of a weak bore, the bore collapse generated turbulence is confined near the bottom boundary layer and the TKE decays at a much slower rate.     

一类三阶非线性微分方程组的全局稳定性条件

应用类比法构造一类三阶非线性微分方程组ｘ＝ｐ（ｙ）,ｙ＝ｑ（ｘ）ｐ（ｙ）＋ｇ（ｚ）,ｚ＝－ｆ（ｘ）－ｈ（ｘ,ｙ）的Ｌｉａｐｕｎｏｖ函数,得到该系统全局稳定性的定理。