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

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

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

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model     

《Applied Ocean Research》2018

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.  相似文献   

礁坪上波浪传播的数值模拟     

何栋彬  房克照  孙家文  刘忠波 《海洋通报》2018,(2):129-138

本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。  相似文献   

Boussinesq相位解析的海岸水动力学数学模型研究进展     

张尧  谢欣  陶爱峰  王岗  陈新平  郑金海 《海洋通报》2018,(5):481-493

由于在平衡计算效率和精度上具有优势,Boussinesq相位解析数学模型研究不断取得突破,已成为波浪和水流精细化模拟的较优解析方式,为海岸工程、环境、地质等问题提供了实用和高效的研究手段。本文对已有Boussinesq类模型的研究进行了评述,深入探讨其重要发展、实际应用和理论瓶颈,从高阶非静压修正、GPU准三维高性能算法编译、波浪破碎和泥沙运移沉积等4个方面提出未来可能的科学突破方向。  相似文献   

Generation of oceanic internal gravity waves by a cyclonic surface stress disturbance     

《Dynamics of Atmospheres and Oceans》2019

Atmospheric cyclones with strong winds significantly impact ocean circulation, regional sea surface temperature, and deep water formation across the global oceans. Thus they are expected to play a key role in a variety of energy transport mechanisms. Even though wind-generated internal gravity waves are thought to contribute significantly to the energy balance of the deep ocean, their excitation mechanisms are only partly understood.The present study investigates the generation of internal gravity waves during a geostrophic adjustment process in a Boussinesq model with axisymmetric geometry. The atmospheric disturbance is set by an idealized pulse of cyclonic wind stress with a Rankine vortex structure. Strength, radius and duration of the forcing are varied. The effect upon wave generation of stratification with variable mixed-layer depth is also examined.Results indicate that internal gravity waves are generated after approximately one inertial period. The outward radial energy flux is dominated by waves having structure close to vertical mode-1 and with frequency close to the inertial frequency. Less energetic higher mode waves are observed to be generated close to the sea floor underneath the storm. The total radiated energy corresponds to approximately 0.02% of the wind input. Deeper mixed-layer conditions as well as weaker stratification reduce this fraction.The low energy transfer rates suggest that other processes that drive vertical motion like surface heat fluxes, turbulent motion, mixed region collapse and storm translation are essential for significant energy extraction by internal gravity waves to occur.  相似文献   

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments     

Soumyaranjan Sahoo  Bhabagrahi Sahoo 《水文研究》2019,33(9):1362-1377

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.  相似文献   

Lagged rejuvenation of groundwater indicates internal flow structures and hydrological connectivity     

Tamara Kolbe  Jean Marçais  Jean-Raynald de Dreuzy  Thierry Labasque  Kevin Bishop 《水文研究》2020,34(10):2176-2189

Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km² forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.  相似文献   

An alternative Boussinesq equation considering the effect of hysteresis on coastal groundwater waves          下载免费PDF全文

Jun Kong  Zhaoyang Luo  Chengji Shen  Guofen Hua  Hongjun Zhao 《水文研究》2016,30(15):2657-2670

This paper presents an alternative Boussinesq equation considering hysteresis effect via a third‐order derivative term. By introducing an improved moisture–pressure retention function, this equation describes, with reasonable precision, groundwater propagation in coastal aquifers subject to Dirichlet boundary condition of different oscillation frequencies. Test results confirmed that it is necessary to consider horizontal and vertical flows in unsaturated zone, because of their variable influences on hysteresis. Hysteresis in unsaturated zone can affect the water table wave number of groundwater wave motion, such as wave damping rate and phase lag. Oscillations with different periods exert different hysteresis effect on wave propagation. Truncation/shrinkage of unsaturated zones also affects the strength of hysteresis. These impacts can be reflected in the alternative Boussinesq equation by adjusting the parameter representing the variation rate of moisture associated with pressure change, as opposed to traditional computationally expensive hysteresis algorithms. The present Boussinesq equation is simple to use and can provide feasible basis for future coupling of groundwater and surface water models. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Study of irregular wave run-up over fringing reefs based on a shock-capturing Boussinesq model     

《Applied Ocean Research》2019

This paper presents the results of a parametric study of irregular wave run-up over fringing reefs using the shock-capturing Boussinesq wave model Funwave-TVD to better understand the role of fringing reefs in the mitigation of wave-driven flooding. Laboratory experiments were newly performed with a typical fringing reef profile and typical hydrodynamic conditions to validate the model. Experimental data shows irregular wave run-ups are dominated by the low-frequency motions and confirms the run-up resonant phenomenon over the back-reef slope, which has been revealed in previous numerical studies. It is demonstrated that irregular wave evolution and run-up over fringing reefs are reasonably reproduced by the present model with a proper grid size. However, the infragravity run-up height and highest 2% run-up height over the back-reef slope are under-predicted due to the underestimation of the infragravity wave height over the reef flat. The validated model was then utilized to model irregular wave transformations and run-ups under different conditions. Through a series of numerical experiments, the effects of key hydrodynamic and reef geometry parameters, including the reef flat width, water depth over the reef flat, fore-reef slope angle and back-reef slope angle, on the irregular wave run-up were investigated. Variations of spectral components of irregular wave run-ups were examined to better understand the physical process underlying the effect of each parameter.  相似文献   

Soil saturation index of salt marshes subjected to spring‐neap tides: a new variable for describing marsh soil aeration condition     

Pei Xin  Badin Gibbes  Ling Li  Zhiyao Song  David Lockington 《水文研究》2010,24(18):2564-2577

Organized spatial distribution of plants (plant zonation) in salt marshes has been linked to the soil aeration condition in the rhizosphere through simplistic tidal inundation parameters. Here, a soil saturation index (ratio of saturation period to tidal period at a soil depth) is introduced to describe the soil aeration condition. This new index captures the effects of not only the tidal inundation period and frequency but also the flow dynamics of groundwater in the marsh soil. One‐dimensional numerical models based on saturated flow with the Boussinesq approximations and a two‐dimensional variably saturated flow model were developed to explore the behaviour of this new soil aeration variable under the influence of spring‐neap tides. Simulations revealed two characteristic zones of soil aeration across the salt marsh: a relatively well aerated near‐creek zone and a poorly aerated interior zone. In the near‐creek zone, soils undergo periodic wetting and drying as the groundwater table fluctuates throughout the spring‐neap cycle. In the interior, the soil remains largely water saturated except for neap tide periods when limited drainage occurs. Although such a change of soil aeration condition has been observed in previous numerical simulations, the soil saturation index provides a clear delineation of the zones that are separated by an ‘inflexion point’ on the averaged index curve. The results show how the saturation index represents the effects of soil properties, tidal parameters and marsh platform elevation on marsh soil aeration. Simulations of these combined effects have not been possible with traditional tidal inundation parameters. The saturation index can be easily derived using relatively simple models based on five non‐dimensional variables. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献