Multi-scale analysis for environmental dispersion in wetland flow

Presented in this work is a multi-scale analysis for longitudinal evolution of contaminant concentration in a fully developed flow through a shallow wetland channel. An environmental dispersion model for the mean concentration is devised as an extension of Taylor’s classical formulation by a multi-scale analysis. Corresponding environmental dispersivity is found identical to that determined by the method of concentration moments. For typical contaminant constituents of chemical oxygen demand, biochemical oxygen demand, total phosphorus, total nitrogen and heavy metal, the evolution of contaminant cloud is illustrated with the critical length and duration of the contaminant cloud with constituent concentration beyond some given environmental standard level.     

Environmental dispersion in wetland flow

For the typical case of a pulsed contaminant emission into a shallow wetland channel, a theoretical analysis is presented in this paper for the decay of the width-averaged mean concentration under environmental dispersion. The velocity profile of a fully developed steady flow through the wetland channel is obtained with that for the well-known plane Poiseuille flow as a special case. An environmental dispersion model for the mean concentration is devised as an extension of Taylor’s classic analysis on dispersion, and corresponding environmental dispersivity is obtained by Aris’s method of moments and illustrated with an asymptotic time variation with stem dominated, transitional, and width-stem dominated stages. Analytical solution for the longitudinal decay of mean concentration due to environmental dispersion is rigorously derived and characterized with multiple time scales.     

Environmental dispersion in a tidal flow through a depth-dominated wetland

Presented in this paper is a theoretical analysis for longitudinal evolution of mean concentration of an environmental emission into a tidal wetland flow. The velocity distribution of the periodic flow through the wetland is derived, with that for a fully developed steady wetland flow included as a special case. The zero-th, first and the growth of the second order moments of the concentration are rigorously obtained by applying Aris’s method of concentration moments to derive the environmental dispersivity. The necessary time needed for the environmental dispersivity to attain a steady oscillating status is analyzed. The effects of some characteristic parameters, especially one representing the impact of vegetation in the wetland on both velocity profile and environmental dispersivity, and another identifying the effect of flow oscillation on the environmental dispersivity, are illustrated in detail. To reflect the dispersion enhancement by the flow oscillation, a typical example is given to characterize the critical length and duration of the contaminant cloud.     

含移动污染源的泥沙反应扩散模型及其解析解

该文研究了广义边界条件下含移动污染源的泥沙反应扩散模型 及其解析解，分析了解析解计算过程中的若干问题. 该文的结果可作为文[2]和[3]的推广.     

Large-strain numerical solution for coupled self-weight consolidation and contaminant transport considering nonlinear compressibility and permeability

Based on the one-dimensional (1D) consolidation equation and advection-dispersion transport equation, this paper presents a large-strain numerical solution for coupled self-weight consolidation and contaminant transport in saturated deforming porous media considering nonlinear compressibility and permeability relationships. The finite difference method is used to solve the governing equations for consolidation and transport. The proposed numerical solution for consolidation accounts for vertical strain, soil self-weight, and nonlinearly changing compressibility and hydraulic conductivity during consolidation. The solution for solute transport accounts for advection, diffusion, mechanical dispersion, linear and nonlinear equilibrium sorption, and porosity-dependent effective diffusion coefficient. The proposed numerical solution is verified against a self-weight consolidation field tank test, an analytical solution in the literature, and the CST1 numerical model. Using the verified solution, a series of parametric study is conducted to investigate the effect of several important parameters on the contaminant transport process for confined disposal of dredged contaminated sediments. The results indicate that the consolidation process and contaminant transport process induced by soil self-weight- can be very different from those induced by the more traditional external surcharge loading. Treating the self-weight loading as traditional external surcharge loading can underestimate the rate of contaminant outflow, especially in the early times. The compressibility and permeability relationships of sediment and the type of loading (i.e., self-weight loading versus external surcharge loading) can all significantly affect the contaminant transport process for confined disposal of dredged contaminated sediment.     

Low-amplitude instability as a premise for the spontaneous symmetry breaking in the new integrable semidiscrete nonlinear system

The new integrable semidiscrete multicomponent nonlinear system characterized by two coupling parameters is presented. Relying upon the lowest local conservation laws the concise form of the system is given and its selfconsistent symmetric parametrization in terms of four independent field variables is found. The comprehensive analysis of quartic dispersion equation for the system low-amplitude excitations is made. The criteria distinguishing the domains of stability and instability of low-amplitude excitations are formulated and a collection of qualitatively distinct realizations of a dispersion law are graphically presented. The loop-like structure of a low-amplitude dispersion law of reduced system emerging within certain windows of adjustable coupling parameter turns out to resemble the loop-like structure of a dispersion law typical of beam-plasma oscillations. Basing on the peculiarities of low-amplitude dispersion law as the function of adjustable coupling parameter it is possible to predict the windows of spontaneous symmetry breaking even without an explicit knowledge of the system Lagrangian function. Having been rewritten in terms of properly chosen modified field variables the reduced four wave integrable system can be qualified as consisting of two coupled nonlinear lattice subsystems, namely the self-dual ladder network and the vibrational ones.     

Taylor dispersion in a packed tube

Presented in this paper is a theoretical analysis for longitudinal scalar spread of mean concentration under a fully developed flow in a tube packed with porous media. A general form of momentum equation for superficial flow in porous media is introduced as a combination of the Navier–Stokes equation and Darcy’s law plus a superficial dispersion term due to phase discontinuity between the fluid flow and solid frame. The analytical solution presented for the fully developed superficial flow includes that for the Poiseulle flow in an evacuated tube as a limiting case. As an extension of Taylor’s classical work on dispersion of soluble matter in solvent flowing slowly through an evacuated tube, a one-dimensional dispersion equation valid for overall environmental assessment of contaminant is rigorously derived by cross-sectionally averaging the superficial mass equation and introducing a closure relation for a new unknown out of the averaging procedure, and corresponding Taylor dispersivity determined is shown to be a generalization of Taylor’s well-known result for the Poiseulle flow.     

二维振荡流动中污染云团的收缩*

当垂向扩散时间尺度与流动的周期相当时,在转流过程中,污染云团将会出现收缩.这时水平剪切分散导数将会出现负值奇性.本文根据作者两维延迟扩散方程[7]: 其中u(t),v (t)为深度平均水平速度.导出X(t,τ),Y(t,τ)坐标位移,D_ij(t,τ)为剪切扩散导数的方程.一般情况下,?D_ij(t,τ)/?τ是正的.不存在奇异性.但在转流的初期.记忆函数D_ij(t,τ)就有可能是负的.本文给出了D_ij和X、Y的解析表示式.     

Advection-Diffusion in Reversing and Oscillating Flows: 1. The Effect of a Single Reversal

Motivated by the need to understand effluent dispersion in shallowtidal waters, a two-dimensional analysis of advection and diffusionin a reversing flow has been carried out. The flow speed varieslinearly with time, passing through zero at time t=0. A pointsource discharges contaminant into this flow at a steady rate,so that water which is close to the source around the time offlow reversal will become highly contaminated. Thus a peak inthe contaminant concentration field will appear, moving downstreamafter the reversal at a speed close to that of the flow. Thisconcentration peak has certain characteristics similar to acloud of contaminant from an instantaneous discharge at t=0.The solution of the advection-diffusion equation is in the formof an integral of concentration fields due to instantaneousreleases of contaminant at all previous times. At large timesafter the flow reversal, asymptotic analysis yields good approximationsto this integral. The use of Laplace's method is equivalentto ignoring longitudinal diffusion (the boundary-layer approximation);however, the expansion obtained in this way is not uniformlyvalid near the concentration peak, indicating that longitudinaldiffusion plays an important role in the development of thispeak. Uniformly valid expansions are obtained for the concentrationaround the peak, and also around the source where the boundary-layerapproximation always breaks down. Numerical integration hasalso been carried out, the results being used to produce contourplots of concentration for various times either side of theflow reversal.     

A solution-scheme for the convective-diffusion equation

This paper presents a versatile solution-scheme for the convective-diffusion equation. A small-time, asymptotic, solution for an instantaneous point source of scalar contaminant is expressed as a three-dimensional, Hermite polynomial expansion and manipulated, using superposition, to generate the contaminant concentration field that results at larger times and for arbitrary, continuous or instantaneous, source contaminant distribution. This equation is commonly used to model contaminant dispersion in complex environmental flows so that the considerable degree of generality, flexibility and efficiency of this solution-scheme highly commends it to this application. The off-diagonal terms in the diffusivity tensor and the non-zero gradient of this term and the mean-velocity field are shown to make a significant contribution to the evolution of the contaminant concentration field resulting from the instantaneous release of contaminant from a point source.

---

Zusammenfassung Die vorliegende Arbeit gibt eine vielseitig brauchbare Lösungsmethode für die konvektive Diffusionsgleichung. Die für kurze Zeiten gültige asymptotische Lösung für eine plötzlich wirkende Punktquelle der diffundierenden Skalargröße wird als dreidimensionale Entwicklung in Hermite'schen Polynomen dargestellt. Dies wird benützt um das Feld für größere Zeiten darzustellen, für beliebige, kontinuierliche oder diskontinuierliche Quellverteilungen. Diese Gleichung wird allgemein benützt als Modell für Dispersion von Schadstoffen in komplizierten Strömungen bei natürlichen Umgebungsverhältnissen; die Allgemeingültigkeit, Anpassungsfähigkeit und Wirksamkeit der Lösung machen die Methode für diese Anwendungen besonders geeignet. Es wird gezeigt, daß die Nichtdiagonalterme des Diffusionstensors sowie die nichtverschwindenden Gradienten dieses Termes und der mittleren Geschwindigkeit wesentliche Beiträge zur Entwicklung des Feldes leisten, die von der Punktquelle erzeugt wurde.

     

多孔介质热弥散系数的分形模型北大核心CSCD

热弥散系数是与流体的物性和多孔介质结构有关的,表征多孔介质传热传质强弱的重要参数.该文建立了分形多孔介质的孔喉结构模型,研究了在孔喉结构处流体由湍流状态变为层流状态的局部水头损失和速度弥散效应,在考虑微观孔喉结构和速度弥散效应的影响下,推导了热弥散系数关系式.研究表明,热弥散系数与孔喉比、孔喉结构个数和迂曲分形维数成正比,与孔隙率和面积分形维数成反比.进一步研究发现,孔喉比在1~150范围内对速度弥散效应有显著影响,流体在孔喉结构处存在局部水头损失,导致速度弥散效应增强,热弥散系数增大.     

The backward group preserving scheme for 1D backward in time advection‐dispersion equation

In this article, we propose a backward group preserving scheme (BGPS) on the advection‐dispersion equation (ADE) for tackling of the contamination problems. The BGPS has been successfully applied on the backward heat conduction problems as well as the backward in time Burgers equation, but it has never been applied to solve the ADE. The BGPS is able to recover the spatial distribution of groundwater contaminant concentration in this work. Several numerical examples are worked out, and we show, based on those numerical examples, that the BGPS is applicable to the ADE and the method can also handle the ADE with piecewise constant dispersion coefficients. When a steep gradient is appeared in the solution, several steps of the BGPS can be used to retrieve the desired initial data and its result is better than the marching‐jury backward beam equation (MJBBE) method as far as our examples are concerned. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Dispersion relations of elastic waves in three-dimensional cubical piezoelectric phononic crystal with initial stresses and mechanically and dielectrically imperfect interfaces

A three-dimensional cubical piezoelectric phononic crystal is theoretically studied in this paper, formed by pre-stressed piezoelectric rectangular blocks and imperfect interfaces. Firstly, transfer matrices of pre-stressed piezoelectric rectangular blocks and imperfect interfaces are obtained based on the structural characteristics. The Bloch waves consist of three groups of sub-coupled elastic waves corresponding with three orthogonal periodical directions in the three-dimensional periodical structure. Furthermore, based on the transfer matrices of typical single cell and the Bloch theorem, it is established that the theoretical model of above-mentioned three-dimensional cubical piezoelectric phononic crystal to obtain the dispersion relations of Bloch waves. Finally, the influences of non-dimensional geometrical parameters (structural and modal parameters) and physical parameters (initial stress and mechanically and dielectrically imperfect interface parameters) on the dispersion relations are discussed based on the graphically numerical results. Numerical calculation results show that the existences of initial stresses and mechanically and dielectrically imperfect interfaces equivalently lead to the alterations of structural flexibilities or rigidities. The theoretical model and numerical discussions will provide a direct guidance of multi-material additive manufacturing for pre-stressed and inhomogeneous periodic structures with partial and global dispersion properties.     

高分辨率差分格式的数值分析与限制因子的构造

本文利用差分方法余项效应理论,分析比较了一些典型的限制因子.对不同的限制因子,格式的表现明显差异主要是由其数值耗散性、色散性强弱不同所致.在分析比较格式的数值耗散性、色散性之后,本文提出了一种新的限制因子,得到的格式在解的剧烈变化区具有更高的分辨率,在光滑区避免了由于数值色散性较强导致的失真.数值试验表明该格式具有较好的性质.     

On Galerkin models for transport in groundwater

The Galerkin finite element model (GFEM) may provide oscillatory results when employed to predict contaminant transport in groundwater unless a very fine mesh is used. Adaptation of a very fine mesh may make the application of the GFEM impractical to field problems. The Petrov—Galerkin finite element models (PGFEMs) can provide oscillation free results for relatively coarser mesh. However, the PGFEM violates the Galerkin principle and introduces large “numerical” dispersion. The objective of this paper has been to develop accurate criteria to improve the applicability of the GFEM to obtain oscillation free accurate results for coarser mesh and compare its performance with that of the PGFEM. It has been shown that the GFEM provides oscillation free accurate results for coarser mesh with Peclet number Pe 20. Further, the GFEM prediction has always been more accurate than the PGFEM for a variety of source configurations and flow fields.     

指数型弥散过程的解析解

对具有指数型弥散系数的弥散过程建立了数学模型,应用积分变换把变系数的偏微分方程变为变系数的常微分方程,应用超几何函数方法和反演技术得到了两类边界条件下的解析解.利用解析解的表达式和计算结果,分析了指数型弥散过程和经典线性弥散过程的差异.     

污染物在非饱和带内运移的流固耦合数学模型及其渐近解

污染物在非饱和带中运移过程是多组分多相渗流问题.在考虑气相的存在对水相影响的前提下,基于流固耦合力学理论,建立了污染物在非饱和带内运移的流固耦合数学模型.对该强非线性数学模型采用摄动法及积分变换法进行拟解析求解,得出了解析表达式.对非饱和带内的孔隙压力分布、孔隙水流速以及污染物的浓度在耦合与非耦合气相条件下的分布规律进行解析计算.对该渐近解与Faust模型的计算结果进行了对比分析,结果表明:该模型解与Faust解基本吻合,且气相作用以及介质的变形对溶质的输运过程产生较大的影响,从而验证了解析表达式的正确性和实用性.这为定量化预报预测污染物在非饱和带中迁移转化和实验室确定压力-饱和度-渗透率三者之间的关系提供了可靠的理论依据.     

液固两相圆柱绕流尾迹内颗粒扩散分布的离散涡数值研究

基于离散涡方法求得的非定常水流场和颗粒的Lagrange运动方程,数值模拟了稀疏液固两相圆柱绕流尾迹内颗粒的扩散分布．获得了流动的涡谱与3种不同St数颗粒（St=0.25,1.0,40）在流场中的分布．通过引入扩散函数来定量表示颗粒在流场中的纵向扩散强度,并计算得到了不同St数颗粒的扩散函数随时间的变化．数值结果揭示出了液固两相圆柱绕流尾迹中的颗粒扩散分布与颗粒的St数和尾涡结构密切相关:1) 中小St数(St=0.25～4.0)颗粒在运动过程中不能进入涡核区,而在旋涡结构的外沿聚集,且颗粒的St数愈大,其越远离涡核区域;2） 在圆柱绕流尾迹区域内,中小St数(St=0.25～4.0)颗粒的纵向扩散强度随其St数的增大而减小．     

Special discontinuities in nonlinearly elastic media

Solutions of a nonlinear hyperbolic system of equations describing weakly nonlinear quasitransverse waves in a weakly anisotropic elastic medium are studied. The influence of small-scale processes of dissipation and dispersion is investigated. The small-scale processes determine the structure of discontinuities (shocks) and a set of discontinuities with a stationary structure. Among the discontinuities with a stationary structure, there are special ones that, in addition to relations following from conservation laws, satisfy additional relations required for the existence of their structure. In the phase plane, the structure of such discontinuities is represented by an integral curve joining two saddles. Special discontinuities lead to nonunique self-similar solutions of the Riemann problem. Asymptotics of non-self-similar problems for equations with dissipation and dispersion are found numerically. These asymptotics correspond to self-similar solutions of the problems.