利用GPU实现基于物理模型的流体运动仿真

纳维-斯托克斯方程是描述流体运动的基础物理模型。首先简要介绍了在计算机上实现该模型的数值解算方法。在流体运动的计算机视觉仿真中,密度渲染才是最终结果,因此,又阐述了流体的密度场与速度场的关系,并且给出了流体仿真的基本流程。由于可编程图形处理器(GPU)具有并行性和速度快等特点,流体模拟的解算和渲染在图形处理器上运行更加高效,但是利用GPU编程实现流体仿真必须要实现3个核心转换,最后以支持3维纹理的Direct3D10为例,给出了利用GPU实现流体运动仿真两个简单实例。     

A finite differences formulation for the linear and nonlinear dynamics of 2D catenary risers

A finite differences (FD) solution method is proposed for the numerical treatment of the dynamic equilibrium problem of 2D catenary risers. The method is based on the so-called Box approximation, which in the scope of the present contribution is applied to the complete nonlinear model as well as to the reduced linearized formulation. The application of the Box method transforms the original governing systems into convenient sets of algebraic equations, which in turn are solved efficiently by the relaxation method. Extensive numerical calculations are presented that describe the dynamic behaviour of the structure and evaluate the amplification in loading due to the dynamic components. The effect of the geometric nonlinearities is assessed through comparative calculations that concern both mathematical formulations examined in the present, i.e. the complete nonlinear, and the reduced linearized model. Special attention is paid to the heave excitations as they amplify significantly the magnitudes of the loading components.     

Numerical calculation of free-surface potential flow around a ship using the modified Rankine source panel method

A modified Rankine source panel method is presented for solving a linearized free-surface flow problem with respect to the double-body potential. The method of solution is based on the distribution of Rankine sources on the hull as well as its image and on the free surface. An iterative algorithm is used for determining the free surface and wave resistance using upstream finite difference operator. A verification of numerical modeling is made using the Wigley hull and the validity of the computer program is examined by comparing the details of wave profiles and wave making resistance with Series 60 model.     

Verification of significant wave representation method

The significant wave representation method is the simplest method for computing the transformation of significant wave height across-shore. However, many engineers are reluctant to use this method because many researchers have pointed out that the method possibly contains a large estimation error. Nevertheless, Rattanapitikon et al. [Rattanapitikon, W., Karunchintadit, R., Shibayama, T., 2003. Irregular wave height transformation using representative wave approach. Coastal Engineering Journal, JSCE 45(3), 489–510.] showed that the wave representation method could be used to compute the transformation of root mean square wave heights. It may also be possible to use it for computing the significant wave height transformation. Therefore, this study was carried out to examine the possibility of simulating significant wave height transformation across-shore by using the significant wave representation method. Laboratory data from small- and large-scale wave flumes were used to calibrate and examine the models. Six regular wave models were applied directly to irregular waves by using the significant wave height and spectral peak period. The examination showed that three regular wave models (with new coefficients) could be used to compute the significant wave height transformation with very good accuracy. On the strength of both accuracy and simplicity of the three models, a suitable model is recommended for computing the significant wave height transformation. The suitable model was also modified for better predictions. The modified model (with different coefficients) can be used to compute either regular wave height or significant wave height transformation across-shore.     

Solution of the propagation of the waves in open channels by the transfer matrix method

Many problems in mechanics can be solved by the use of the transfer matrix method. The use of this method in hydraulics engineering is not widespread and only limited studies are available. In this study, linearized St. Venant equations were used and the use of transfer matrix in ocean engineering was investigated for long waves in open channels, and numerical application was carried out. The results obtained through the transfer matrix method, which is quite easy to use, program and comprehend, showed similar results obtained from the characteristics method and finite differences method.     

极限波浪运动特性的非线性数值模拟

利用时域高阶边界元方法建立了模拟极限波浪运动的完全非线性数值模型,其中自由水面满足完全非线性自由水面条件.采用半混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶Runge-Kutta方法更新下一时间步的波面和速度势,同时应用镜像格林函数消除水槽两个侧面和底面上的积分.研究中利用波浪聚焦的方法产生极限波浪,并且在水槽中开展了物理模型实验,将测点试验数据与数值结果进行了对比,两者吻合得很好.对极限波浪运动的非线性和流域内速度分布进行了研究.     

Simulation of nonlinear wave run-up with a high-order Boussinesq model

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions. As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability to accurately simulate a moving wet–dry boundary is of considerable practical importance within coastal engineering, and the extension described in this work significantly improves the nearshore versatility of the present high-order Boussinesq approach.     

Detection of nighttime sea fog/stratus over the Huanghai Sea using MTSAT-1R IR data

A dual channel difference (DCD) method is applied to detect nighttime sea fog/stratus over the Huanghai Sea using the infrared (IR) data of shortwave (3.5–4.0 μm) and longwave (10.3–11.3 μm) channels from the Multi-functional Transport Satellite (MTSAT)-1R, i.e., shortwave minus longwave brightness temperature difference (SLTD). Twenty-four sea fog events over the Huanghai Sea during March to July of 2006 and 2007 are chosen to determine a suitable value of SLTD for nighttime sea fog/stratus detection, and ...     

Modified Moving Particle Semi-implicit methods for the prediction of 2D wave impact pressure

As a gridless particle method, the MPS (Moving Particle Semi-implicit) method has proven useful in a wide variety of engineering applications including free-surface hydrodynamic flows. Despite its wide range of applicability, the MPS method suffers from some shortcomings such as non-conservation of momentum and spurious pressure fluctuation. By introducing new formulations for the pressure gradient and a new formulation of the source term of the Poisson Pressure Equation (PPE), and by allowing a slight compressibility, we have proposed modified MPS methods for the prediction of wave impact pressure on a coastal structure. The improved performance of the modified methods is shown through the simulation of numerous wave impact problems (including the impacts by a dam break flow, a flip-through and two cases of slightly-breaking waves) in comparison with the experimental data.     

Nonlinear random wave-induced drag force on a vegetation field

This paper provides a practical method by which the drag force on a vegetation field beneath nonlinear random waves can be estimated. This is achieved by using a simple drag formula together with an empirical drag coefficient given by Mendez et al. (Mendez, F.J., Losada, I.J., Losada, M.A., 1999. Hydrodynamics induced by wind waves in a vegetation field. J. Geophys. Res. 104 (C8), 18383–18396). Effects of nonlinear waves are included by using Stokes second order wave theory where the basic harmonic motion is assumed to be a stationary Gaussian narrow–band random process. An example of calculation is also presented.