基于GPU加速的Boussinesq类波浪传播变形数值模型     

孙家文  朱桐  房克照  刘忠波 《海洋工程》2020,38(2):111-119

Boussinesq波浪模型是一类相位解析模型,在时域内求解需要较高的空间和时间分辨率以保证计算精度。为提高计算效率,有必要针对该类模型开展并行算法的研究。与传统的中央处理器(CPU)相比,图形处理器(GPU)有大量的运算器,可显著提高计算效率。基于统一计算设备架构CUDA C语言和图形处理器,实现了Boussinesq模型的并行运算。将本模型的计算结果同CPU数值模拟结果和解析解相比较,发现得到的结果基本一致。同时也比较了CPU端与GPU端的计算效率,结果表明,GPU数值模型的计算效率有明显提升,并且伴随数值网格的增多,提升效果更为明显。  相似文献   

URANS investigation of the interaction between the free surface and a shallowly submerged underwater vehicle at steady drift     

《Applied Ocean Research》2019

An axisymmetric underwater vehicle (UV) at a steady drift angle experiences the complex three-dimensional crossflow separation. This separation arises from the unfavorable circumferential pressure gradient developed from the windward side toward the leeward side. As is well known, the separated flow in the leeward side gives rise to the formation of a pair of vortices, which affects considerably the forces and moments acting on the UV. In this regard, the main purpose of the present study is to evaluate the role of the leeward vortical flow structure in the hydrodynamic behavior of a shallowly submerged UV at a moderate drift angle traveling beneath the free surface. Accordingly, the static drift tests are performed on the SUBOFF UV model using URANS equations coupled with a Reynolds stress turbulence model. The simulations are carried out in the commercial code STARCCM+ at a constant advance velocity based on Froude number equal to F_n = 0.512 over submergence depths and drift angles ranging from h = 1.1D to h = ∞ and from β = 0 to β = 18.11°, respectively. The validation of the numerical model is partially conducted by using the existing experimental data of the forces and moment acting on the totally submerged bare hull model. Significant interaction between the low-pressure region created by the leeward vortical flow structure and the free surface is observed. As a result of this interaction, the leeward vortical flow structure appears to be largely responsible for the behavior of the forces and moments exerted on a shallowly submerged UV at steady drift.  相似文献   

Modeling of iceberg drift in the marginal ice zone of the Barents Sea     

《Applied Ocean Research》2019

The paper is focused on the analysis of the drift of tabular iceberg observed in 2009 in the marginal ice zone of the North-West Barents Sea. Momentum balance equations are derived from the Kirchhoff equations describing plane motion of solid body in an ideal fluid. Field works performed on the drift iceberg and on the drift icenear the iceberg are described. Results of the field works and numerical simulations of the iceberg drift and rotation are performed and discussed. It is shown that acceleration of water flow around the iceberg has visible influence on the iceberg drift. Kinetic energy balance of drift iceberg is used to estimate the forces applied to the iceberg by the drift ice.  相似文献   

附有限制条件的间接平差秩亏时解法初探          下载免费PDF全文

秦永宽  黄声享  张书毕 《海洋测绘》2009,29(1):7-9

在附有限制条件的间接平差的应用中,若所列观测方程中实际用到的未知参数的个数少于所选参数的总数u时,就会出现法方程系数阵秩亏的问题。为了解决此问题,在理论公式推导的基础上,给出了一种新的解决方法,并结合算例进行了验证。  相似文献   

地心坐标反解大地坐标的非线性方程组牛顿迭代法分析     

李琦  王小辉 《测绘科学》2009,34(6)

目前的一些地心坐标向大地坐标的转换模型在计算速度,稳定性或精度方面存在一定的局限性。本文探讨了非线性方程组数值迭代的求解方法及其在MATLAB 7.0中的编程实现,并将结果与基于一元三次方程求解的严密方法做了比较分析,结果证明该方法是可以在实际中参考应用的。  相似文献   

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

胡香  张威巍  周玉柱  刘芳 《测绘学院学报》2009,(3)

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

Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: Application to skipjack tuna   总被引：2，自引：1，他引：1  

Inna Senina  John Sibert  Patrick Lehodey 《Progress in Oceanography》2008,78(4):319-335

A Spatial Ecosystem and Population Dynamic Model (SEAPODYM) is used in a data assimilation study aiming to estimate model parameters that describe dynamics of Pacific skipjack tuna population on ocean-based scale. The model based on advection–diffusion–reaction equations explicitly predicts spatial dynamics of large pelagic predators, while taking into account data on several mid-trophic level components, oceanic primary productivity and physical environment. In order to improve its quantitative ability, the model was parameterized through assimilation with commercial fisheries data, and optimization was carried out using maximum likelihood estimation approach. To address the optimization task we implemented an adjoint technique to obtain an exact, analytical evaluation of the likelihood gradient. We conducted a series of computer experiments in order to (i) determine model sensitivity with respect to variable parameters and, hence, investigate their observability; (ii) estimate observable parameters and their errors; and (iii) justify the reliability of the computed solution. Parameters describing recruitment, movement, habitat preferences, natural and fishing mortality of skipjack population were analysed and estimated. Results of the study suggest that SEAPODYM with achieved parameterization scheme can help to investigate the impact of fishing under various management scenarios, and also conduct forecasts of a given species stock and spatial dynamics in a context of environmental and climate changes.  相似文献   

A two-dimensional finite volume hydrodynamic model for coastal areas: Model development and validation     

A.H. Daoud  K.A. Rakha  A.G. Abul-azm   《Ocean Engineering》2008,35(1):150-164

A two dimensional implicit finite volume scheme for solving the shallow-water equations is developed. The effects of the Coriolis force, surface wind stress, and waves are included. A non-uniform rectilinear forward staggered grid is used with Cartesian coordinates. The time integration is performed using the Euler implicit technique. The convective flux is treated using the deferred correction method. The viscous terms are discretized using a second order central difference approximation. The SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm is used for coupling the velocity components and the water elevation gradient for the water level correction. The system of equations is solved sequentially using the Strongly Implicit Procedure (SIP). To simulate wave driven current, a phase averaged wave model is used first to simulate wave transformation and calculate radiation stresses. The performance of the developed model is validated for different sources of external forces and different combinations of boundary conditions. The validation cases include tidal circulation in a harbor and wave induced currents behind a breakwater parallel to the coastline. The model is finally applied to simulate the flow pattern in a closed artificial lagoon and along the coastline near Damietta Port located along the Northern coast of Egypt. Results of the developed model agree well with the published results for the considered cases.  相似文献   

Simulation of nonlinear wave run-up with a high-order Boussinesq model   总被引：2，自引：0，他引：2  

David R. Fuhrman  Per A. Madsen 《Coastal Engineering》2008

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.  相似文献   

Velocity potential formulations of highly accurate Boussinesq-type models     

Harry B. Bingham  Per A. Madsen  David R. Fuhrman 《Coastal Engineering》2009

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.  相似文献