关于非线性浅水波理论中的孤立波解

本文在[1]的基础上,详尽地得到了Boussinesq方程和KdV方程的孤立波解,并对波高和波形进行了细致的分析。为了更好地比较,本文还给出了高阶摄动的孤立波解。     

Boussinesq方程的永形波的级数解

对于用层平均速度势(或速度)与波高给出的非线性等水深Boussinesq方程,与KdV方程一样,存在两类永形波—孤立波和周期波。不同的是,对于Boussinesq方程,孤立波和周期波一样,是以反函数的积分形式给出的,这在使用上是不方便的。本文对这两类波,分别给出了相应的级数形式的正函数解。对于孤立波,级数的收敛速度是十分快的,且文中还给出了一个具有高精度的十分简单的近似解,对于周期波,只要有关参数选得合适,解的收敛速度也是相当快的。这就有效地克服了反函数形式解在使用上的不方便。     

Water table fluctuation with a random initial condition

The nonlinear Boussinesq equation is used to understand water table fluctuations in various ditch drainage problems. An approximate solution of this equation with a random initial condition and deterministic boundary conditions, recharge rate and aquifer parameters has been developed to predict a transient water table in a ditch-drainage system. The effects of uncertainty in the initial condition on the water table are illustrated with the help of a synthetic example. These results would find applications in ditch-drainage design.Notation A / tanh t - a lower value of the random variable representing the initial water table height at the mid point - a+b Upper value of the random variable representing the initial water table height at the midpoint - B tanh t - C 4/ - h variable water table height - h mean of the variable water table height - h _m variable water table height at the mid point - h _m mean of the variable water table height at the mid point - K hydraulic conductivity - L half spacing between the ditches - m ₀ initial water table height at the mid point - N Uniform rate of recharge - S specific yield - t time of observation - x distance measured from the ditch boundary - (4/SL)(NK)^1/2 - (L/4)(N/K)^1/2 - dummy integral variable     

改进型Boussinesq方程高精度紧致差分显格式

采用一种高精度的紧致差分显格式对改进型Boussinesq方程进行数值求解;采用具有TVD性质的三阶Runge-Kutta方法进行预报,用三次样条函数进行校正,时间精度可达到四阶;在空间离散上采用六阶精度的三点紧致显格式进行计算;运用以上数值格式对Beji和Nadaoka改进型Boussinesq方程进行了求解,求解证明:高精度的数值结果和已知的试验结果吻合良好.作为验证算例,同时对波浪在台阶上的传播进行了模拟,从效果对比上可以看出,所得结果明显比Kittitanasuan的计算结果更靠近试验值.     

Numerical solution of the ‘classical’ Boussinesq system

We consider the ‘classical’ Boussinesq system of water wave theory, which belongs to the class of Boussinesq systems modelling two-way propagation of long waves of small amplitude on the surface of water in a horizontal channel. (We also consider its completely symmetric analog.) We discretize the initial-boundary-value problem for these systems, corresponding to homogeneous Dirichlet boundary conditions on the velocity variable at the endpoints of a finite interval, using fully discrete Galerkin-finite element methods of high accuracy. We use the numerical schemes as exploratory tools to study the propagation and interactions of solitary-wave solutions of these systems, as well as other properties of their solutions.     

Numerical implementation of the asymptotic boundary conditions for steadily propagating 2D solitons of Boussinesq type equations

In the present paper, a difference scheme on a non-uniform grid is constructed for the stationary propagating localized waves of the 2D Boussinesq equation in an infinite region. Using an argument stemming form a perturbation expansion for small wave phase speeds, the asymptotic decay of the wave profile is identified as second-order algebraic. For algebraically decaying solution a new kind of nonlocal boundary condition is derived, which allows to rigorously project the asymptotic boundary condition at the boundary of a finite-size computational box. The difference approximation of this condition together with the bifurcation condition complete the algorithm. Numerous numerical validations are performed and it is shown that the results comply with the second-order estimate for the truncation error even at the boundary lines of the grid. Results are obtained for different values of the so-called ‘rotational inertia’ and for different subcritical phase speeds. It is found that the limits of existence of the 2D solution roughly correspond to the similar limits on the phase speed that ensure the existence of subcritical 1D stationary propagating waves of the Boussinesq equation.     

BOUSSINESQ MODELLING OF NEARSHORE WAVES UNDER BODY FITTED COORDINATE

A set of nonlinear Boussinesq equations with fully nonlinearity property is solved numerically in generalized coordinates,to develop a Boussinesq-type wave model in dealing with irregular computation boundaries in complex nearshore regions and to facilitate the grid refinements in simulations.The governing equations expressed in contravariant components of velocity vectors under curvilinear coordinates are derived and a high order finite difference scheme on a staggered grid is employed for the numerical implementation.The developed model is used to simulate nearshore wave propagations under curvilinear coordinates,the numerical results are compared against analytical or experimental data with a good agreement.     

Numerical Solution of Boussinesq Equations to Simulate Dam-Break Flows

To investigate the effect of nonhydrostatic pressure distribution, dam-break flows are simulated by numerically solving the one-dimensional Boussinesq equations by using a fourth-order explicit finite-difference scheme. The computed water surface profiles for different depth ratios have undulations near the bore front for depth ratios greater than 0.4. The results obtained by using the Saint Venant equations and the Boussinesq equations are compared to determine the contribution of individual Boussinesq terms in the simulation of dam-break flow. It is found that, for typical engineering applications, the Saint Venant equations give sufficiently accurate results for the maximum flow depth and the time to reach this value at a location downstream of the dam.     

当均布垂直荷载作用时半无限弹性体内侧向变位的计算

文献［１］提出了计算地基表面侧向变位的方法，本文则就均布垂直荷载作用的情况下，提出了计算半无限弹性体内侧向变位的方法，可作为文献［１］的补充。     

Heat/mass transfer from a neutrally buoyant sphere by mixed natural and forced convection in a simple shear flow

Building on the work of Yang et al. in 2011, the finite difference method and the Boussinesq approximation were applied to solve the time‐dependent Navier‐Stokes, convection diffusion and continuity equations in spherical coordinates. An idealized condition, the mass transfer from a neutrally buoyant sphere in a horizontal simple shear flow with natural convection was numerically simulated for the first time in this work. In the hybrid transfer case, the outwardly spiraling streamlines enhanced the transfer process, but the counter‐gravity spiraling streamlines near the sphere hindered the natural convection and the spatial dilution action weakened the natural convection transfer process. These competing effects led to nonmonotonic behavior of the Nusselt number with Reynolds number. Results from these previously undocumented cases were summarized into correlations for predicting Nusselt numbers at finite Reynolds numbers for various Grashof and Prandtl numbers. © 2018 American Institute of Chemical Engineers AIChE J, 64: 2816–2827, 2018