非线性发展方程的精确解

基于齐次平衡法的思想,利用G′G-展开法求得KdV-Burgers-Kuramoto方程的精确解,这些解更具有一般性,其中包括了双曲函数解、三角函数周期解、有理数解。同时,这些解对于研究湍流运动和不稳定现象有重要意义。为更好的描述这些解,给出了三角函数周期解及扭状解的数值模拟图。     

Bifurcation analysis of wind-driven flows with MOM4

In this paper, the methodology of bifurcation analysis is applied to the explicit time-stepping ocean model MOM4 using a Jacobian–Free Newton–Krylov (JFNK) approach. We in detail present the implementation of the JFNK method in MOM4 but restrict the preconditioning technique to the case for which the density distribution is prescribed. For a prescribed density field case, we present bifurcation diagrams, for the first time in MOM4, for the wind-driven ocean circulation. In addition, we show that the JFNK method can reduce the spin-up time to a steady equilibrium in MOM4 considerably if an accurate solution is required.     

A view on capsizing under direct and parametric wave excitation

The behavior of a ship that is simultaneously heaving, swaying and rolling in waves is examined by using an early developed model by Thompson et al. [Thompson, J. M. T., Rainey, R. C. T. and Soliman, M. S. (1992) Mechanics of ship capsize under direct and parametric wave excitation. Philosophical Transactions of the Royal Society London A 338, 471-490]. A second-order nonlinear differential equation has been used as a mathematical model which is analyzed by using numerical methods to obtain geometrical phase space techniques of nonlinear dynamics. Bifurcation diagrams and Poincaré maps are two of the geometrical techniques used in this study. The result of the calculation is given to allow the naval architect to make a self comparison.     

An Efficient Hydrodynamic Model for Surface Waves

In the present study,a semi-implicit finite difference model for non-bydrostatic,free-surface flows is analyzed and discussed.The governing equations are the three-dimensional free-surface Reynolds-averaged Navier-Stokes equations defined on a general,irregular domain of arbitrary scale.At outflow,a combination of a sponge layer technique and a radiation boundary condition is applied to minimize wave reflection.The equations are solved with the fractional step method where the hydrostatic pressure component is determined first,while the non-hydrostatic component of the pressure is computed from the pressure Poisson equation in which the coefficient matrix is positive definite and symmetric.The advectiou and horizontal viscosity terms are discretized by use of a semi-Lagrangian approach.The resulting model is computationally efficient and unrestricted to the CFL condition.The developed model is verified against analytical solutions and experimental data,with excellent agreement.     

适用于沙坝上Bragg反射的二阶Boussinesq方程数学模型及其数值验证

在二阶 Boussinesq 方程基础上,通过引入含水深导数项对该方程进行了理论上的改进,使得该方程在应用于无限沙坝 Bragg反射问题时与理论解析解在更大范围内符合.基于该改进的高阶 Boussinesq 方程,在非交错网格下建立了混合 4 阶的Adams-Bashforth- Moulton 格式的数学模型.将数值模型应用到有限个连续沙坝上波浪传播变形问题的数值模拟中,通过两点法给出数值波浪反射系数,将这些反射系数与已有的实验数据进行对比,对比表明改进后的模型计算出的反射系数与实验结果吻合更好,这验证了本文理论改进的有效性.     

植被斜坡岸滩海啸波消减数值模拟研究

An explicit one-dimensional model based on the shallow water equations(SWEs) was established in this work to simulate tsunami wave propagation on a vegetated beach. This model adopted the finite-volume method(FVM)for maintaining the mass balance of these equations. The resistance force caused by vegetation was taken into account as a source term in the momentum equation. The Harten–Lax–van Leer(HLL) approximate Riemann solver was applied to evaluate the interface fluxes for tracing the wet/dry transition boundary. This proposed model was used to simulate solitary wave run-up and long-periodic wave propagation on a sloping beach. The calibration process suitably compared the calculated results with the measured data. The tsunami waves were also simulated to discuss the water depth, tsunami force, as well as the current speed in absence of and in presence of forest domain. The results indicated that forest growth at the beach reduced wave energy loss caused by tsunamis. A series of sensitivity analyses were conducted with respect to variable parameters(such as vegetation densities, wave heights, wave periods, bed resistance, and beach slopes) to identify important influences on mitigating tsunami damage on coastal forest beach.     

Are There Differences in the Response of Natural Stand and Plantation Biomass to Changes in Temperature and Precipitation? A Case for Two-needled Pines in Eurasia

A comparative discussion of the advantages and disadvantages of natural stands and plantations, including in terms of their productivity and stability, began from the moment of the first forest plantings and continues to this day. In the context of the progressive replacement of natural forests by plantations due to deforestation, the question of how will change the carbon storage capacity of forest cover when replacing natural forests with artificial ones in a changing climate becomes extremely relevant. This article presents the first attempt to answer this question at the transcontinental level on a special case for two-needled pine trees (subgenus Pinus L.). The research was carried out using the database compiled by the authors on the single-tree biomass structure of forest-forming species of Eurasia, in particular, data of 1880 and 1967 of natural and plantation trees, respectively. Multi-factor regression models are calculated after combining the matrix of initial data on the structure of tree biomass with the mean January temperature and mean annual precipitation, and their adequacy indices allow us to consider them reproducible. It is found that the aboveground and stem biomass of equal-sized and equal-aged natural and plantation trees increases as the January temperature and precipitation rise. This pattern is only partially valid for the branches biomass, and it has a specific character for the foliage one. The biomass of all components of plantation trees is higher than that of natural trees, but the percent of this excess varies among different components and depends on the level of January temperatures, but does not depend at all on the level of annual precipitation. A number of uncertainties that arose during the modeling process, as well as the preliminary nature of the obtained regularities, are noted.     

Scenarios of tsunami amplitudes in the north eastern coast of Sea of Marmara generated by submarine mass failure

In this paper the tsunamis resulting from a submarine mass failure such as slides and slumps triggered by earthquakes or other environmental effects, which is settled at the bottom of the north eastern Sea of Marmara are examined in one sample region. As the solution method, one hybrid method is developed. The main objective of this method is to combine an analytical solution presenting near-field tsunami amplitudes above the submarine mass failure with a numerical solution indicating the tsunami amplitudes in the coastal regions. For this purpose, one common linear boundary between analytical and numerical solution domains is defined. Movements of Submarine Mass Failures (SMF) are modeled using one simple kinematics source model and the amplitudes of the tsunamis at the region that are closer to the landslide are computed by using the analytical method. SMF is modeled approximately from the bottom geometry, and an average depth is used. Scenarios of SMF are established depending on the velocities and thicknesses of the failure, and near-field tsunami amplitudes are obtained in the open sea during the source time. After the source times, the solutions are found in the numerical region using TELEMAC-2D software system with the mentioned boundary above. In this boundary, the output of the analytical solutions is taken as the boundary conditions or the disturbances for the numerical method. With these disturbances, the numerical method is performed and the amplitudes are calculated in the coastal area. The generation, propagation and coastal amplifications of the tsunamis are illustrated at some certain points and regions both in the open sea and near the coast line. The results have been visualized and discussed.     

Numerical modeling of nonlinear resonance of semi-enclosed water bodies: Description and experimental validation

The purpose of this work is to a present a numerical model to solve a set of modified Boussinesq equations to analyse nonlinear resonance of semi-enclosed water bodies. The equations are solved on a finite element unstructured grid in order to achieve an optimal mesh resolution with the local geometry. The model is able to simulate long time lapses and realistic forcing in real bathymetries with a reasonable computational cost. To validate the numerical results, a set of experiments was carried out in a physical model of two adjacent elongated basins. Comparisons between numerical and experimental results for different geometries and nonlinear conditions show that the model is able to simulate with an excellent agreement the transient nonlinear resonant process.     

Determination of fractional energy loss of waves in nearshore waters using an improved high-order Boussinesq-type model

Fractional energy losses of waves due to wave breaking when passing over a submerged bar are studied systematically using a modified numerical code that is based on the high-order Boussinesq-type equations.The model is first tested by the additional experimental data,and the model's capability of simulating the wave transformation over both gentle slope and steep slope is demonstrated.Then,the model's breaking index is replaced and tested.The new breaking index,which is optimized from the several breaking i...