非线性抛物型方程反问题的一种新算法

利用反问题新算法—时间域正演反演法研究非线性抛物型方程的逆时反问题。该方法处理反问题的主要思路是先求解相对应的正问题,获得解在特定时间网格处的近似值;然后构造一个合适的全纯映射,在象空间中获得解在时间网格对应点处的近似值;最后利用解析函数的唯一延拓性质实现反问题的逆时间反演。数值模拟例子说明了方法的有效性。     

波浪和风生浪流作用下的近岸海底泥沙运移模型的推广和应用

本文利用数值计算和经验公式相结合的方法,获得了在波浪和风生流和波浪共同作用下海底泥沙输移的二维模型.在此基础上,讨论了海底泥沙输移在风速、风向和泥沙粒径不同等三种情况下的变化特征;指出了经验公式对泥沙输移规律的适用范围以及风速、风向对于近岸处泥沙输移重要性.这一研究结果对近岸结构物的选址设计具有实际应用价值.     

Practical modified scheme of linear shallow-water equations for distant propagation of tsunamis

A simple but practical numerical model describing a distant propagation of tsunamis is newly proposed by introducing an additional term to the existing modified scheme. The numerical dispersion of the proposed model is manipulated to replace the physical dispersion of the linear Boussinesq equations without any limitation. The new model developed in this study is applied to propagation of a Gaussian hump over a constant water depth and the predicted free surface displacements are compared with available analytical solutions. A very reasonable agreement is observed.     

A Gibbs–Pitzer function for high-salinity seawater thermodynamics

A high-salinity Gibbs function for seawater is derived from Pitzer equations of the sea salt components, in conjunction with the 2003 Gibbs function of seawater for low salinities. Various properties, computed from both formulations by thermodynamic rules, are compared with each other, and with high-salinity measurements. The new Gibbs–Pitzer function presented in this paper is valid in the range 0–110 g kg⁻¹ in absolute salinity, −7 to +25 °C in temperature, and 0–100 MPa in applied pressure. The formulation is expressed in the International Temperature Scale 1990 (ITS-90), and is consistent with the International Standard for Fluid Water (IAPWS-95), and with the 2005/2006 equations of state of ice Ih.     

An application of Boussinesq equations to Bragg reflection of irregular waves

A numerical model based on the second-order fully nonlinear Boussinesq equations of Wei et al. [1995. Journal of Waterway, Port, Coastal and Ocean Engineering 121 (5), 251-263] is developed to simulate the Bragg reflection of both regular and irregular surface waves scattered by submerged bars. Particularly for incident regular waves, the computed results are observed to agree very well with the existing experimental data as presented by Davies and Heathershaw [1984. Journal of Fluid Mechanics 144, 419-446] and Kirby and Anton [1990. Proceedings of the 22nd International Conference on Coastal Engineering, ASCE, New York, pp. 757–768). In the case of incident irregular waves, the simulated results reveal that the distribution of Bragg reflection from irregular waves becomes more flat than that of regular waves. Due to lack of experimental data, the numerical results for incident irregular waves are compared with those of the evolution equation of the mild-slope equation [Hsu et al., 2002 Proceedings of the 24th Ocean Engineering Conference in Taiwan, pp. 70–77 (in Chinese)]. In addition, several parameters such as the number of bars, the relative height of bars and the spacing of bars affecting Bragg reflection are also discussed.     

海浪感应电磁场的理论计算

根据麦克斯韦电磁理论,在地磁场中运动的海水将产生感应电磁场。利用一个简单的数学物理模型对感应电磁场在海水内部的分布进行了计算,结果表明海浪产生的电磁场明显依赖于海浪波动的周期及浪高。在100 m的海水深度内,海浪产生的磁感应强度的大小为纳特数量级,而电场强度的大小为几个微伏每米。在同一海水深度处,磁感应强度随海水波动的周期呈现近线性变化,而电场强度的大小有一个极值,该极值随海水深度的增加向长周期方向移动。海浪产生的电磁场是影响海洋电磁探测数据精度的主要噪声之一。     

Numerical simulation of irregular wave overtopping against a smooth sea dike

Based on the filtered Navier-Stokes equations and Smagorinsky turbulence model,a numerical wave flume is developed to investigate the overtopping process of irregular waves over smooth sea dikes.Simulations of fully nonlinear standing wave and regular wave’s run-up on a sea dike are carried out to validate the implementation of the numerical wave flume with wave generation and absorbing modules.To model stationary ergodic stochastic processes,several cases with different random seeds are computed for each specified irregular wave spectrum.It turns out that the statistical mean overtopping discharge shows good agreement with empirical formulas,other numerical results and experimental data.     

Wave height transformation and set-up between a submerged permeable breakwater and a seawall

In this study,we investigated wave transformation and wave set-up between a submerged permeable breakwater and a seawall.Modified time-dependent mild-slope equations,which involve parameters of the porous medium,were used to calculate the wave height transformation and the mean water level change around a submerged breakwater.The numerical solution is verified with experimental data.The simulated results show that modulations of the wave profile and wave set-up are clearly observed between the submerged breakwater and the seawall.In contrast to cases without a seawall,the node or pseudo-node of wave height evolution can be found between the submerged breakwater and the seawall.Higher wave set-up occurs if the nodal or pseudo-nodal point appears near the submerged breakwater.We also examined the influence of the porosity and friction factor of the submerged permeable breakwater on wave transformation and set-up.     

Numerical solution of the asymmetric water impact of a wedge in three degrees of freedom

Impact problems associated with water entry have important applications in various aspects of naval architecture and ocean engineering. Estimation of hydrodynamic impact forces especially during the first instances after the impact is very important and is of interest. Since the estimation of hydrodynamic impact load plays an important role in safe design and also in evaluation of structural weight and costs, it is better to use a reliable and accurate prediction method instead of a simple estimation resulted by analyzing methods. In landing of flying boats, some phenomena such as weather conditions and strong winds can cause asymmetric instead of symmetric descent. In this paper, a numerical simulation of the asymmetric impact of a wedge, as the step of a flying boat, considering dynamic equations in two-phase flow is taken into account. The dynamic motion of the wedge in two-phase flow is solved based on finite volume method with volume of fluid (VOF) scheme considering dynamic equations. Then the effects of different angles of impact and water depth on the velocity change and slamming forces in an asymmetric impact are investigated. The comparison between the simulation results and experimental data verifies the accuracy of the method applied in the present study.     

Alternative forms of the higher-order Boussinesq equations: Derivations and validations

An alternative form of the Boussinesq equations is developed, creating a model which is fully nonlinear up to O(μ⁴) (μ is the ratio of water depth to wavelength) and has dispersion accurate to the Padé [4,4] approximation. No limitation is imposed on the bottom slope; the variable distance between free surface and sea bottom is accounted for by a σ-transformation. Two reduced forms of the model are also presented, which simplify O(μ⁴) terms using the assumption ε = O(μ^2/3) (ε is the ratio of wave height to water depth). These can be seen as extensions of Serre's equations, with dispersions given by the Padé [2,2] and Padé [4,4] approximations. The third-order nonlinear characteristics of these three models are discussed using Fourier analysis, and compared to other high-order formulations of the Boussinesq equations. The models are validated against experimental measurements of wave propagation over a submerged breakwater. Finally, the nonlinear evolution of wave groups along a horizontal flume is simulated and compared to experimental data in order to investigate the effects of the amplitude dispersion and the four-wave resonant interaction.