Stochastic modelling of dispersion in shallow water

A random walk model to describe the dispersion of pollutants in shallow water is developed. By deriving the Fokker-Planck equation, the model is shown to be consistent with the two-dimensional advection-diffusion equation with space-varying dispersion coefficient and water depth. To improve the behaviour of the model shortly after the deployment of the pollutant, a random flight model is developed too. It is shown that over long simulation periods, this model is again consistent with the advection-diffusion equation. The various numerical aspects of the implementation of the stochastic models are discussed and finally a realistic application to predict the dispersion of a pollutant in the Eastern Scheldt estuary is described.     

The nonlinear spin-up of a stratified ocean

Abstract

The adjustment of a nonlinear, quasigeostrophic, stratified ocean to an impulsively applied wind stress is investigated under the assumption that barotropic advection of vortex tube length is the most important nonlinearity. The present study complements the steady state theories which have recently appeared, and extends earlier, dissipationless, linear models.

In terms of Sverdrup transport, the equation for baroclinic evolution is a forced advection-diffusion equation. Solutions of this equation subject to a “tilted disk” Ekman divergence are obtained analytically for the case of no diffusion and numerically otherwise. The similarity between the present equation and that of a forced barotropic fluid with bottom topography is shown.

Barotropic flow, which is assumed to mature instantly, can reverse the tendency for westward propagation, and thus produce regions of closed geostrophic contours. Inside these regions, dissipation, or equivalently the eddy field, plays a central role. We assume that eddy mixing effects a lateral, down-gradient diffusion of potential vorticity; hence, within the closed geostrophic contours, our model approaches a state of uniform potential vorticity. The solutions also extend the steady-state theories, which require weak diffusion, by demonstrating that homogenization occurs for moderately strong diffusion.

The evoiution of potential vorticity and the thermocline are examined, and it is shown that the adjustment time of the model is governed by dissipation, rather than baroclinic wave propagation as in linear theories. If dissipation is weak, spin-up of a nonlinear ocean may take several times that predicted by linear models, which agrees with analyses of eddy-resolving general circulation models. The inclusion of a western boundary current may accelerate this process, although dissipation will still play a central role.     

Analysis of key parameters in a diffusion type beach profile evolution model

Diffusion type formulations are commonly used in beach profile evolution models. The practical idea behind that is to map the behaviour of the beach profile onto a simple mathematical model that exhibits the same behaviour under defined operating conditions. The success of this approach is based on the accurate determination of key parameters in the diffusion model that govern its behaviour, using observed beach behaviour in the field. In order to determine these parameters, i.e. diffusion coefficient and a time and space varying source function, we used observations of historic beach profiles at Milford-on-Sea beach in Christchurch Bay, Dorset, United Kingdom. The relationship between the diffusion coefficient and Dean's equilibrium profile was investigated, leading to a new interpretation of the diffusion coefficient in terms of the sediment characteristics. The analysis also shows the significance of the diffusion process in the medium to long term evolution of the beach profile. A canonical correlation analysis (CCA) was undertaken in order to identify patterns of behaviour between wave conditions and source terms, and the possible correlations between them. The analysis provides strong evidence of a useful link between the source term in the simple dynamical equation and the distribution of wave steepness.     

A process-based model for sediment transport under various wave and current conditions

The purpose of this study is to investigate the capability of a newly developed process-based model for sediment transport under a wide variety of wave and current conditions.The model is based on the first-order boundary layer equation and the sediment advection-diffusion equation.In particular,a modified low Reynolds number k-e model is coupled to provide the turbulence closure.Detailed model verifications have been performed by simulating a number of laboratory experiments,covering a considerable range of hydrodynamic conditions such as sinusoidal waves,asymmetric waves and wave-current interactions.The model provides satisfactory numerical results which agree well with the measured results,including the time-averaged/dependent sediment concentration profiles and sediment flux profiles,as well as the time series of concentration at given elevations.The observed influences of wave orbital velocity amplitude,wave period and sediment grain size are correctly reproduced,indicating that the fundamental physical mechanisms of those processes are properly represented in the model.It is revealed that the present model is capable of predicting sediment transport under a wide range of wave and current conditions,and can be used to further study the morphodynamic processes in real coastal regions.     

Truncation errors of selected finite difference methods for two-dimensional advection-diffusion equation with mixed derivatives

The spread of a passive contaminant in an open-channel reach is considered with use of a two-dimensional advection-diffusion equation with the included off-diagonal dispersion coefficients. This paper presents the calculation of truncation errors, namely numerical diffusion and numerical dispersion for various finite difference schemes. The accuracy of the considered finite-difference approximations is analysed by deriving and studying the relevant modified partial differential equation.     

An integral equation for the homogeneous Neumann condition1

Modelling the theoretical response of several important geophysical systems involves the solution of Poisson's equation with homogeneous Neumann boundary conditions (i.e. a zero normal gradient) imposed over either open or closed surfaces. A simple integral equation solution to this problem is derived from first principles. It is applicable to both types of surface and in this respect represents an improvement on existing integral equation techniques. However, the present surface integral equation displays a strong singularity of order 1/R³ which requires an appropriate interpretation for its implementation. A comparison of some numerical results with analytical data taken from the literature demonstrates that the proposed integral equation technique is suitably robust, accurate and efficient for practical application in geophysical interpretation.     

A displacement finite element modelling for convection-diffusion problems

The development of a displacement finite element formulation and its application to convective transport problems is presented. The formulation is based on the introduction of a generalized quantity defined as transport displacement. The governing equation is expressed in terms of this quantity and by using generalized coordinates a variational form of the governing equation is obtained. This equation may be solved by any numerical method, though it is of particular interest for application of the finite element method. Two finite element models are derived for the solution of convection-diffusion boundary value problems. The performance of the two element models is discussed and numerical results are given for different cases of convection and diffusion with two types of boundary conditions. The numerical results obtained show not only the efficiency of the numerical models in handling pure convection, pure diffusion and mixed convection-diffusion problems, but also good stability and accuracy. The applications of the developed numerical models are not limited to diffusion-convection problems but can also be applied to other types of problems such as mass transfer, hydrodynamics and wave propagation.     

The boundary layer of the residence time field

The residence time of a tracer in a control domain is usually computed by releasing tracer parcels and registering the time when each of these tracer parcels cross the boundary of the control domain. In this Lagrangian procedure, the particles are discarded or omitted as soon as they leave the control domain. In a Eulerian approach, the same approach can be implemented by integrating forward in time the advection–diffusion equation for a tracer. So far, the conditions to be applied at the boundary of the control domain were uncertain. We show here that it is necessary to prescribe that the tracer concentration vanishes at the boundary of the control domain to ensure the compatibility between the Lagrangian and Eulerian approaches. When we use the Constituent oriented Age and Residence time Theory (CART), this amounts to solving the differential equation for the residence time with boundary conditions forcing the residence time to vanish at the open boundaries of the control domain. Such boundary conditions are likely to induce the development of boundary layers (at outflow boundaries for the tracer concentration and at inflow boundaries for the residence time). The thickness of these boundary layers is of the order of the ratio of the diffusivity to the velocity. They can however be partly smoothed by tidal and other oscillating flows.     

沉积物-水界面污染物迁移扩散的研究进展

污染物在沉积物-水界面的迁移扩散对研究其环境生物地球化学循环过程和评估水生态系统质量具有重要意义.本文回顾了关于沉积物-水界面的基本研究历程,重点介绍沉积物-水界面的垂向结构以及扩散边界层（DBL）的作用,展示污染物在沉积物-水界面的多维度分布（一维垂向、二维平面和三维立体）及其在沉积物-水界面的扩散过程,详细总结影响污染物在沉积物-水界面迁移扩散的环境因素（包括温度、溶解氧或氧化还原条件、pH值、离子强度或盐度、沉积物组成、共存污染物、溶解性有机质、水动力条件、生物扰动、微生物以及其他因素）,讨论当前关于沉积物-水界面污染物扩散通量估算几种方法的优缺点,最后,对沉积物-水界面污染物扩散研究在未来发展需要关注的几个方面进行了展望.     

Adaptive hierarchical grid model of water-borne pollutant dispersion

Water pollution by industrial and agricultural waste is an increasingly major public health issue. It is therefore important for water engineers and managers to be able to predict accurately the local behaviour of water-borne pollutants. This paper describes the novel and efficient coupling of dynamically adaptive hierarchical grids with standard solvers of the advection–diffusion equation. Adaptive quadtree grids are able to focus on regions of interest such as pollutant fronts, while retaining economy in the total number of grid elements through selective grid refinement. Advection is treated using Lagrangian particle tracking. Diffusion is solved separately using two grid-based methods; one is by explicit finite differences, the other a diffusion-velocity approach. Results are given in two dimensions for pure diffusion of an initially Gaussian plume, advection–diffusion of the Gaussian plume in the rotating flow field of a forced vortex, and the transport of species in a rectangular channel with side wall boundary layers. Close agreement is achieved with analytical solutions of the advection–diffusion equation and simulations from a Lagrangian random walk model. An application to Sepetiba Bay, Brazil is included to demonstrate the method with complex flows and topography.     

Nodal domain integration model of two-dimensional advection-diffusion processes

The nodal domain integration method is applied to a two-dimensional advection-diffusion process in an anisotropic inhomogeneous medium. The domain is discretised into the union of irregular triangle finite elements with vertex-located nodal points and a linear trial function is used to approximate the governing flow equation's state variable in each element. Non-linear parameters are assumed quasi-constant for small durations in time in each element. The resulting numerical model represents the Galerkin and subdomain integration weighted residual methods and the integrated finite difference method as special cases. Both Dirichlet and Neumann boundary conditions are accommodated in a manner similar to the Galerkin finite element approach.     

Suspended load in flows on erodible bed

Steady state suspended-load of sediment transported in flow over erodible beds usually is treated by the advection-diffusion approach, though in recent years, it is being treated as a two-phase flow phenomenon incorporating kinetics of sediment particles. Among the advection-diffusion approaches, Rouse＇s equation is the well-known, although a number of researchers in later periods have attempted to improve it by modifying the mixing length concept taking into account other aspects. In this paper, the advection-diffusion approach and associated logarithmic law of flow velocity are revisited. It is concluded from the logarithmic law that the Reynolds shear stress is a linear function of height above the bed, which reduces to bed shear stress in the case of a long horizontal channel. As a consequence, it is shown that the volumetric concentration of sediment is best approximated by the sum of two power laws of height above the bed. An equation is derived for the suspended-load transport rate in terms of elementary functions.     

Nonlinear sloshing in a floating‐roofed oil storage tank under long‐period seismic ground motion

A hybrid analytical and FEM is proposed to investigate the nonlinear sloshing in a floating‐roofed oil storage tank under long‐period seismic ground motion. The tank is composed of a rigid cylindrical wall and a flat bottom, whereas the floating roof is treated as an elastic plate undergoing large deflection. The contained liquid is assumed to be inviscid and incompressible, and the flow is assumed to be irrotational. The method of analysis is based on representation of the liquid motion by superposing the analytical modes that satisfy the Laplace equation and the rigid wall and bottom boundary conditions. The FEM is then applied to solve the remaining kinematic and dynamic boundary conditions at the moving liquid surface coupled with the nonlinear equation of motion of the floating roof. This requires only the discretization of the liquid surface and the floating roof into finite elements, thus leading to a computationally efficient and accurate method compared with full numerical analysis. As numerical examples to illustrate the applicability of the proposed method, two oil storage tanks with single‐deck type floating roofs damaged during the 2003 Tokachioki earthquake are studied. It is shown that the nonlinear oscillation modes with the circumferential wave numbers 0, 2 and 3 caused by the finite liquid surface elevation as well as the membrane action due to large deflection of the deck produce excessively large stresses in the pontoon, which may cause the catastrophic failure of pontoon followed by the submergence of the roof. Copyright © 2012 John Wiley & Sons, Ltd.     

Motion of wetting fronts moving into partially pre-wet soil

We study the motion of wetting fronts for vertical infiltration problems as modeled by Richards’ equation. Parlange and others have shown that wetting fronts in infiltration flows can be described by traveling wave solutions. If the soil layer is not initially dry, but has an initial distribution of water content then the motion of the wetting front will change due to the interaction of the infiltrating flow with the pre-existing soil conditions. Using traveling wave profiles, we construct simple approximate solutions of initial-boundary value problems for Richards’ equation that accurately describe the position and moisture distribution of the wetting front. We show that the influences of surface boundary conditions and initial conditions produce shifts to the position of the wetting front. The shifts can be calculated by examining the cumulative infiltration, and are validated numerically for several problems for Richards’ equation and the linear advection–diffusion equation.     

Closed-form solutions for unsaturated flow under variable flux boundary conditions

It is shown that from any solution of the linear diffusion equation, we may construct a solution of a realistic form of the Richards equation for unsaturated flow. Compared to the usual direct linearization method, our inverse approach involves a quite different sequence of transformations. This opens the possibility of exact solutions with a wider variety of continuously varying flux boundary conditions. Closed-form solutions are presented for two examples. In these, the varying water flux boundary conditions resemble (i) the passage of a peaking storm and (ii) the continuous opening of a valve preceding a steady water supply. Unlike earlier more systematic approaches to this problem, our method does not require the numerical solution of an integral equation.     

The surface boundary condition in nonhydrostatic ocean models

With increasing resolution in numerical ocean models, nonhydrostatic pressure effects have to be accounted for. In sigma-coordinate mode split ocean models, this pressure may be regarded as a pressure correction. An elliptic equation must be solved for the nonhydrostatic pressure, and the gradients are used to correct the provisional hydrostatic velocity components in each time step. The focus in the present work is on the surface boundary condition for the elliptic equation. In the literature, both Dirichlet and Neumann boundary conditions are suggested and applied. To investigate the sensitivity of the numerical results to the choice of boundary condition, three numerical experiments are performed. The first and second experiments are studies of the propagation and steepening of nonlinear internal waves. The first study is on tank scale and the second experiment is on ocean scale. In the tank-scale experiment, the density and the flow fields are very robust to the choice of boundary condition. In the ocean-scale experiment, the waves produced with a Dirichlet boundary condition become more damped than the waves produced with a Neumann boundary condition. The third study involves a surface buoyant jet. It is shown that well-known characteristics of the plume front are reproduced with a Neumann boundary condition, but the rotating turbulent core of this front is lost with a Dirichlet condition. It is accordingly argued that the appropriate surface boundary condition in mode split nonhydrostatic ocean models is the Neumann condition.     

从瞬变电磁扩散场到拟地震波场的全时域反变换算法

将瞬变电磁满足的扩散方程转变为波动方程,然后利用地震类成像方法实现瞬变电磁虚拟波场成像,是实现瞬变电磁三维反演的有效手段之一.为了实现由扩散场到虚拟波场的转换,文中采用预条件正则化共轭梯度法求解波场反变换问题.首先,对几种离散方式进行比较,采用条件数最小的离散方式进行离散;然后选择最优的正则化参数,并利用超松弛预条件技术对系数矩阵进行预条件处理;最后,利用共轭梯度法进行迭代求解.超松弛预条件有效降低了系数矩阵的条件数,正则化方法使得反变换得到的波场稳定、可靠,共轭梯度法能够保证计算快速收敛.将反变换结果与已知虚拟波场函数对比,证明算法稳定、可信.将文中算法结果与前人研究结果进行对比,说明方法效果.通过实测数据的波场变换处理给出了文中方法的实际应用效果.结合反变换算法,对不同参数模型进行分析,总结了虚拟波场在色散介质中的传播规律.     

波动方程三维吸收边界分解与拟合深度偏移方法

在偏移问题中引入吸收边界条件，既可以消除由人工边界激发的虚假反射，从而提高剖面质量。又可以减少计算工作量．本文讨论了三维吸收边界条件方程，提出了求解具有吸收边界条件的三维波动方程偏移定解问题的分解与拟合方法。理论分析与合成记录及野外实际地震资料处理结果表明，本文方法为一有效的三维吸收边界深度偏移方法。     

Solution of the transport equations using a moving coordinate system

A convection-diffusion equation arises from the conservation equations in miscible and immiscible flooding, thermal recovery, and water movement through desiccated soil. When the convection term dominates the diffusion term, the equations are very difficult to solve numerically. Owing to the hyperbolic character assumed for dominating convection, inaccurate, oscillating solutions result. A new solution technique minimizes the oscillations. The differential equation is transformed into a moving coordinate system which eliminates the convection term but makes the boundary location change in time. We illustrate the new method on two one-dimensional problems: the linear convection-diffusion equation and a non-linear diffusion type equation governing water movement through desiccated soil. Transforming the linear convection diffusion equation into a moving coordinate system gives a diffusion equation with time dependent boundary conditions. We apply orthogonal collocation on finite elements with a Crank-Nicholson time discretization. Comparisons are made to schemes using fixed coordinate systems. The equation describing movement of water in dry soil is a highly non-linear diffusion-type equation with coefficients varying over six orders of magnitude. We solve the equation in a coordinate system moving with a time-dependent velocity, which is determined by the location of the largest gradient of the solution. The finite difference technique with a variable grid size is applied, and a modified Crank-Nicholson technique is used for the temporal discretization. Comparisons are made to an exact solution obtained by similarity transformation, and with an ordinary finite difference scheme on a fixed coordinate system.