A comparison of a semi-implicit with an explicit scheme in a three-dimensional hydrodynamic model

As the capacity of computers increases, the size and resolution of numerical models can be increased. In tidal models, however, using an explicit scheme together with high spatial resolution results in an unreasonably small time-step demanded by the CFL condition for stability. This condition is usually the most restrictive and applies to the propagation of the gravity wave.A fully nonlinear three-dimensional model has been developed, using the Galerkin method in the vertical dimension, in which the gravity wave terms are treated by an alternating-direction implicit scheme, while the friction, viscous and advective terms are treated explicitly. This permits a stable solution with a longer time-step than that required in an equivalent explicit scheme, while not needing as much computational effort as a fully implicit scheme. This semi-implicit model is compared with an explicit model in terms of efficiency, accuracy, and stability. Tidal and wind-driven flows and free oscillations in a rectangular sea model are examined, using various boundary conditions. The semi-implicit scheme takes approximately 1.1 times as long to run (per time-step) on the CRAY-1 computer as the explicit scheme.     

自适应多尺度第二代小波配点法探地雷达数值模拟

基于第二代小波变换的提升方案构造了插值小波,将雷达波场函数进行了二维小波变换,得到所有尺度上与计算网格相联系的小波系数和尺度系数.对所有尺度上的小波系数进行分析,根据解的局部性与小波系数阈值的控制,实现网格压缩和配点的自适应调节.保留大于给定阈值的小波系数及对应网格点,令小于给定阈值的小波系数为零,并舍弃其对应网格点.达到光滑区域采用较少的计算网格点,在奇异性较大的区域采用较多的计算网格点的目的.通过对自适应网格进行邻域校正、重构检查等附加修正,推导了场值更新的显式时间迭代方案.最后,以均匀、阶梯与复杂三个典型GPR模型为例,与常规数值计算结果对比表明：自适应小波配点法（AWCM）利用第二代小波的多尺度分解和快速变换的特点,可以使计算网格随着时间步适应解的移动和变化,允许计算资源更有效地使用,具有高压缩率,达到跟踪奇异性的目的,特别适合于探地雷达正演中波传问题的模拟.     

A semi-implicit three-dimensional numerical model for non-hydrostatic pressure free-surface flows on an unstructured,sigma grid

A semi-implicit 3-D numerical formulation for solving non-hydrostatic pressure free-surface flows on an unstructured,sigma grid is proposed.Pressure-splitting and 9 semi-implicit methods are inherited and reformed from Casulli’s z-coordinate formulation.The non-orthogonal sigma-coordinate transformation leads to additional terms. The resulting linear system for the non-hydrostatic correction is diagonally dominant but unsymmetric,and it is solved by the BiCGstab method.In contrast with z-coordinate non-hydrostatic models,the new model fits vertical boundaries much better,which is important for the long-time simulation of sediment transport and riverbed deformation.A lock-exchange density flow is computed to determine whether the new scheme is able to simulate non-hydrostatic free-surface flows.The new model is further verified using the field data of a natural river bend of the lower Yangtze River.Good agreement between simulations and earlier research results,field data is obtained, indicating that the new model is applicable to hydraulic projects in real rivers.     

Comparison of finite difference and element models of internal tides on the Malin–Hebrides shelf

A three-dimensional baroclinic finite element model with a coarse and fine (i.e. local refinement along the shelf edge) grid is used to examine the influence of shelf edge grid refinement upon the internal tide generation and propagation off the west coast of Scotland. Comparisons are made with observations in the region and with a published solution using a finite difference model. The calculations show that provided that the finite element grid is refined in the internal tide generation area and the adjacent region through which the internal tide propagates, then a numerically accurate solution is obtained. In the regions of strong internal tide generation with a local grid refinement, internal wave energy can accumulate at small scales and must be removed by a scale-selective filter.     

A method of accelerating transport simulation when groundwater pumping is simulated

In solving groundwater transport problems with numerical models, the computation time (CPU processing time) of transport simulation is approximately inversely proportional to the transport time-step size. Therefore, large time-step sizes are favorable for achieving short computation time. However, transport time-step size must be sufficiently small to avoid numerical instability if an explicit scheme is used (and to guarantee enough model accuracy if an implicit scheme is used). For a transport model involving groundwater pumping, a small transport time-step size is often required due to the high groundwater velocities near the pumping well. Small grid spacing often specified near the pumping well also limits the time-step size. This paper presents a method to increase transport time-step size in a transport model when groundwater pumping is simulated. The key to this approach is to numerically decrease the groundwater seepage velocities in grid cells near the pumping well by increasing the effective porosity so that the transport time-step size can be increased without violating stability constraints. Numerical tests reveal that by using the proposed method, the computation time of transport simulation can be reduced significantly, while the transport simulation results change very little.     

Assessment of Anthropogenic Impact on Marine Ecosystems and Biological Resources in the Process of Oil and Gas Field Development in the Shelf Area

Current methodology of environmental impact assessment in connection with the environmental consequences of hydrocarbons production in the shelf area is analyzed. Basing on the ecosystem approach, a scheme of environmental impact estimates is suggested, envisaging the use of a set of gradations (scales) to characterize spatial and temporal scope of impacts and their consequences, as well as criteria (thresholds) of impact permissibility, taking into account the natural variability in the population and ecosystem parameters. Estimates of the environmental and fishery-related consequences of hydrocarbon production in the sea shelf area at different stages of the process are compared to the available measured and calculated data. The suggested scheme of environmental impact assessment and criteria of permissible impacts are recommended for expert analysis, forecasting, and monitoring of the environmental situation in the sea shelf areas of Russia under the conditions of anthropogenic impact.     

Resolving frontal structures: on the payoff using a less diffusive but computationally more expensive advection scheme

This article presents some advantages using a shape-preserving total variation diminishing (TVD) advection scheme in an ecosystem model. The superbee flux-limiter has been used to the second-order Lax–Wendroff advection scheme to make it TVD. We performed simulations for three shelf sea regions with different characteristic time scales, namely, the North Sea, the Barents Sea, and the Baltic Sea. To explore the advantages, we also performed reference runs with the much simpler and computationally cheaper upwind advection scheme. Frontal structures are much better resolved with the TVD scheme. The bottom salinity in the Baltic Sea stays at realistic values throughout the 10-year simulation with the TVD scheme, while with the upwind scheme, it drifts towards lower values and the permanent haline stratification in the Baltic is almost completely eroded within one seasonal cycle. Only when applying TVD for both the vertical and horizontal advections the model succeeded to preserve haline stratification in the decadal simulation. Lower trophic level patterns are far better reproduced with the TVD scheme, and for the estimated cod larval survival, the advantages seem to be even stronger. Simulations using the TVD-derived prey fields identified distinct regions such as Dogger Bank to favor potential larvae survival (PLS), which did not appear as particularly favorable in the upstream simulations. The TVD scheme needs about 25?% more time on the central processing unit (CPU) in case of a pure hydrodynamic setup with only two scalar state variables (Barents Sea application). The additional CPU time cost increases for a coupled physical–biological model application with a large number of state variables. However, this is more than compensated by all the advantages found, and, hence, we conclude that it is worthwhile to use the TVD scheme in our ecosystem model.     

Resolving frontal structures: on the payoff using a less diffusive but computationally more expensive advection scheme

This article presents some advantages using a shape-preserving total variation diminishing (TVD) advection scheme in an ecosystem model. The superbee flux-limiter has been used to the second-order Lax–Wendroff advection scheme to make it TVD. We performed simulations for three shelf sea regions with different characteristic time scales, namely, the North Sea, the Barents Sea, and the Baltic Sea. To explore the advantages, we also performed reference runs with the much simpler and computationally cheaper upwind advection scheme. Frontal structures are much better resolved with the TVD scheme. The bottom salinity in the Baltic Sea stays at realistic values throughout the 10-year simulation with the TVD scheme, while with the upwind scheme, it drifts towards lower values and the permanent haline stratification in the Baltic is almost completely eroded within one seasonal cycle. Only when applying TVD for both the vertical and horizontal advections the model succeeded to preserve haline stratification in the decadal simulation. Lower trophic level patterns are far better reproduced with the TVD scheme, and for the estimated cod larval survival, the advantages seem to be even stronger. Simulations using the TVD-derived prey fields identified distinct regions such as Dogger Bank to favor potential larvae survival (PLS), which did not appear as particularly favorable in the upstream simulations. The TVD scheme needs about 25 % more time on the central processing unit (CPU) in case of a pure hydrodynamic setup with only two scalar state variables (Barents Sea application). The additional CPU time cost increases for a coupled physical–biological model application with a large number of state variables. However, this is more than compensated by all the advantages found, and, hence, we conclude that it is worthwhile to use the TVD scheme in our ecosystem model.

     

Experiments with a numerical model of coastal currents and tidal mixing fronts

A three-dimensional numerical model is used to simulate the development of disturbances on shelf-sea coastal currents and fronts. The model, which has a free surface, uses a finite difference grid ☐ scheme based on sigma coordinates. It has a semi-implicit scheme for the barotropic flow and a hydrid advection scheme to retain sharp fronts. The results demonstrate that (i) eddy formation follows changes at the inflow of a coastal current, (ii) a simple radiation boundary condition at the outflow produces nearly identical results for different outflow boundary positions, (iii) eddy growth, with matching behaviour of surface and bottom fronts, follows a small displacement on a tidal mixing front and (iv) effects of friction and mixing can significantly alter the behaviour of the front and the relative strength of the cyclonic and anticyclonic eddies formed.     

A 3D unstructured non-hydrostatic ocean model for internal waves

A 3D non-hydrostatic model is developed to compute internal waves. A novel grid arrangement is incorporated in the model. This not only ensures the homogenous Dirichlet boundary condition for the non-hydrostatic pressure can be precisely and easily imposed but also renders the model relatively simple in its discretized form. The Perot scheme is employed to discretize horizontal advection terms in the horizontal momentum equations, which is based on staggered grids and has the conservative property. Based on previous water wave models, the main works of the present paper are to (1) utilize a semi-implicit, fractional step algorithm to solve the Navier-Stokes equations (NSE); (2) develop a second-order flux-limiter method satisfying the max–min property; (3) incorporate a density equation, which is solved by a high-resolution finite volume method ensuring mass conservation and max–min property based on a vertical boundary-fitted coordinate system; and (4) validate the developed model by using four tests including two internal seiche waves, lock-exchange flow, and internal solitary wave breaking. Comparisons of numerical results with analytical solutions or experimental data or other model results show reasonably good agreement, demonstrating the model’s capability to resolve internal waves relating to complex non-hydrostatic phenomena.     

Numerical studies of large-amplitude internal waves shoaling and breaking at shelf slopes

Hydro carbon fields beyond the shelf break are presently being explored and developed, which has increased the scientific focus in this area. Measurements from the slopes reveal large variability in temperature and velocity, and some of the observed events are due to interactions between large-amplitude oscillations of the thermocline and the topography. The present study focuses on the strong currents that are generated near the seabed during shoaling and breaking of internal waves along shelf slopes. The parameter regime used is similar to the one for the Nordic Seas. The results show that, during shoaling of large internal waves along (gentle) slopes, the energy is transferred towards smaller scales and strong velocities (over 1 m s^− 1) can be generated. To resolve all scales involved is still not feasible, and therefore, the model results are sensitive to the grid size and the subgrid scale closure.     

数值差分格式及格点设置对土壤温度模拟结果的影响

土壤温度是反映气候系统和生态系统能量循环的重要地球物理学参量,土壤温度的模拟精度直接影响着气候系统模式以及陆面物理过程模式的模拟结果.为了提高模式对土壤温度的模拟能力,本文利用土壤热扩散方程的傅里叶解析解定量研究了差分方案、格点设置以及时间步长对土壤温度模拟结果的影响;提出了一种优化的格点设置方案,并利用巴丹吉林沙漠观测数据检验了该方案的性能.研究结果表明:三种差分方案中,显式方案的模拟误差最小,Crank-Nicolson方案其次,隐式方案的模拟误差最大;每一种格点设置方案均存在一个使模拟结果误差最小的最优化时间步长;常用格点设置方案的最优化时间步长为5358 s,最小标准差为0.156 K,优化方案的最优化时间步长为1694 s,最小标准差为0.0465 K;取时间步长为1800 s时,采用常用格点设置方案,巴丹吉林沙漠10 cm深度土壤温度模拟结果的标准差为1.61 K,而采用优化方案,模拟结果的标准差降至0.21 K,改进效果明显.     

波动方程的应力-速度有限元解耦交叠格式

本基于有限差分交叠格式和解耦有限元方法的基本概念,以应力－速度为变量,提出了求解波动的应力－速度有限元解耦交叠格式,这一格式不仅时空解耦,而且为显式,它适合于线性及非线性波动问题的数值模拟,已有的应力－速度有限元交叠格式（即格子法）为本的特例。通过解析解数值检验表明,本建议的方法具有较高的精度,而格子法计算精度较低。     

Implicit time integration for pseudodynamic tests

The use of unconditionally stable implicit time integration techniques for pseudodynamic tests has been recently proposed and advanced by several researchers. Inspired by such developments, a pseudodynamic test scheme based on an unconditionally stable implicit time integration algorithm and dual displacement control is presented in this paper. The accuracy of the proposed scheme is proved with error-propagation analysis. It is shown by numerical examples and verification tests that the error-correction method incorporated can eliminate the spurious higher-mode response, which can often be excited by experimental errors. The practicality of the proposed scheme lies in the fact that the implementation is as easy as that of explicit schemes and that the convergence criteria required are compatible with the accuracy limits of ordinary test apparatus.     

Application of modified Patankar schemes to stiff biogeochemical models for the water column

In this paper, we apply recently developed positivity preserving and conservative Modified Patankar-type solvers for ordinary differential equations to a simple stiff biogeochemical model for the water column. The performance of this scheme is compared to schemes which are not unconditionally positivity preserving (the first-order Euler and the second- and fourth-order Runge–Kutta schemes) and to schemes which are not conservative (the first- and second-order Patankar schemes). The biogeochemical model chosen as a test ground is a standard nutrient–phytoplankton–zooplankton–detritus (NPZD) model, which has been made stiff by substantially decreasing the half saturation concentration for nutrients. For evaluating the stiffness of the biogeochemical model, so-called numerical time scales are defined which are obtained empirically by applying high-resolution numerical schemes. For all ODE solvers under investigation, the temporal error is analysed for a simple exponential decay law. The performance of all schemes is compared to a high-resolution high-order reference solution. As a result, the second-order modified Patankar–Runge–Kutta scheme gives a good agreement with the reference solution even for time steps 10 times longer than the shortest numerical time scale of the problem. Other schemes do either compute negative values for non-negative state variables (fully explicit schemes), violate conservation (the Patankar schemes) or show low accuracy (all first-order schemes).     

An externally forced barotropic circulation model for the New York Bight

The sea level and the barotropic, frictional circulation response for the New York Bight are used to demonstrate the effects of external sea-level forcing, bathymetry, and variable friction. The governing equation is the steady, integrated vorticity equation and is computed by finite differencing over a curvilinear grid conforming to the 10- and 100-m isobaths and extending for 250 km alongshore. The boundary conditions are based on the hypothesis that the dynamics of the shelf are driven by the external sea-level gradient and the coastal no-flux condition; and consequently the conditions at the lateral boundaries are dependent thereon. Therefore, the external sea-level slope must be independently specified, and the lateral boundary conditions must be dependently generated. The diabathic component of the external sea slope forces the calm wind circulation by its effect on the transport through the upstream boundary; and the parabathic component has also an important modifying effect by forcing a shelf convergent transport. The parabathic sea slope at the coast is independent of its offshore value, being instead a direct product of the coastal boundary condition.The bottom friction is expressed as related to the sea level through a bottom length parameter and a veer angle, both of which are taken to increase shoreward. An additional bottom stress component, related to the surface stress, is determined for bottom depths less than the Ekman depth. Such bottom stress variability produces significant alterations in the nearshore flow field, over the constant bottom stress formulation, by reducing it and causing it to veer downgradient and downwind in the nearshore.The model is forced by different wind directions and the results are discussed. The circulations generally conform to the observed mean flow patterns, but with several smaller-scale features. The strong bathymetric feature of the Hudson Shelf Valley causes a polarized up- and downvalley flow for winds with an eastward or westward component, respectively. Under mean westerly winds, there is a divergence in the shelf valley flow at about the 60-m isobath. The Apex gyre existing off the western tip of Long Island becomes more extensive for winds changing from northeast to southwest. Mean flow reversals (to the northeast) occur off both Long Island and New Jersey for wind directions changing counterclockwise from northwest to southeast and from west to east, respectively. Southeastward transport over the outer New Jersey shelf tends to be enhanced by wind and external sea-level conditions; and the transport over the New Jersey midshelf, particularly in the lee of the shelf valley, tends to be weak and variable also under these mean conditions.     

瞬变电磁三维FDTD正演多分辨网格方法

瞬变电磁三维时域有限差分（FDTD）正演的网格剖分受最小网格尺寸、时间步长、边界条件、目标尺寸、模型尺寸等的影响,结构化网格一直存在最小网格尺寸受限于异常目标尺寸的矛盾;尽管非均匀网格能够在保证模型尺寸的前提下尽可能的降低网格数量,但由于Yee网格结构的限制,非均匀网格不能无限制的扩大单一方向的尺寸,这是为了避免边界网格区域出现长宽比过大的畸形网格,影响计算精度甚至导致结果发散.在非均匀网格剖分的基础上,本文提出了瞬变电磁三维FDTD正演的多尺度网格方法,即首先使用较大尺寸的粗网格进行第一次剖分,然后在希望加密的区域进行二次剖分,使计算域中包含粗、细两套网格.尽管细网格包含在粗网格内部,但其具有Yee网格的全部属性,因而可以在网格中设置不同的电性参数模拟不同形状的目标.基于Maxwell方程组推导了细网格内电场和磁场的迭代公式,基于泰勒展开给出了设置粗、细网格后产生的内部边界条件,使电磁场的传播在粗、细网格和时间步进上得到统一.采用均匀半空间中包含三维低阻异常的经典模型和三维接触带复杂模型进行精度验证,发现多分辨网格方法计算结果满足精度要求.使用"L"型异常模型计算采用多分辨网格方法和不采用多分辨网格的传统FDTD方法对比计算效率,发现多分辨网格算法能够显著提高计算效率,并能够保证计算精度.     

A hybrid absorbing boundary condition for elastic staggered‐grid modelling

We recently proposed an efficient hybrid scheme to absorb boundary reflections for acoustic wave modelling that could attain nearly perfect absorptions. This scheme uses weighted averaging of wavefields in a transition area, between the inner area and the model boundaries. In this paper we report on the extension of this scheme to 2D elastic wave modelling with displacement‐stress formulations on staggered grids using explicit finite‐difference, pseudo‐implicit finite‐difference and pseudo‐spectral methods. Numerical modelling results of elastic wave equations with hybrid absorbing boundary conditions show great improvement for modelling stability and significant absorption for boundary reflections, compared with the conventional Higdon absorbing boundary conditions, demonstrating the effectiveness of this scheme for elastic wave modelling. The modelling results also show that the hybrid scheme works well in 2D rotated staggered‐grid modelling for isotropic medium, 2D staggered‐grid modelling for vertically transversely isotropic medium and 2D rotated staggered‐grid modelling for tilted transversely isotropic medium.