Coupled consolidation and contaminant transport model for simulating migration of contaminants through the sediment and a cap

Capping contaminated sediments in waterways is an alternate remediation technique to dredging and is typically much cheaper than dredging. When cap material is placed on top of contaminated sediment, it has both a short-term and long-term hydraulic impact on the underlying sediment. A numerical model of consolidation, based on a nonlinear finite strain theory for a consolidating fine-grained sediment bed was developed. The nonlinear equation of consolidation was solved in a material (or reduced) coordinate using an explicit finite difference numerical scheme. An one-dimensional advection–diffusion equation with sorption and decay was solved on a convective coordinate using a finite volume total variation diminishing (TVD) scheme for the contaminant concentration within the consolidating sediment. The contaminant transport model was coupled with the consolidation model. The time and space varying porosities, permeabilities, and advective velocities computed by the consolidation model were linked to the transport model at the same time level. A number of benchmark tests that are relevant to the consolidation of a fine-grained sediment were designed and tested. The relative contribution of consolidation-induced transient advective velocities on the migration of a contaminant during consolidation was also investigated. The coupled model performance was validated by simulating the transport of hazardous chemicals under consolidation in a confined aquatic disposal (CAD) site in the Lower Duwamish Waterway, Seattle.     

对流扩散方程在成品油顺序输送混油分析中的应用

本文研究了对流与扩散对成品油顺序输送混油过程的影响;推导了紊流条件下,描述混油过程的对流占优的扩散方程;将该方程分解为纯对流方程和纯扩散方程,分别应用特征线法和差分法求解,数值计算结果和实际操作经验相符,能很好地解释混油的形成和发展.     

Convergence of an operator splitting method on a bounded domain for a convection-diffusion-reaction system

We solve a convection-diffusion-sorption (reaction) system on a bounded domain with dominant convection using an operator splitting method. The model arises in contaminant transport in groundwater induced by a dual-well, or in controlled laboratory experiments. The operator splitting transforms the original problem to three subproblems: nonlinear convection, nonlinear diffusion, and a reaction problem, each with its own boundary conditions. The transport equation is solved by a Riemann solver, the diffusion one by a finite volume method, and the reaction equation by an approximation of an integral equation. This approach has proved to be very successful in solving the problem, but the convergence properties where not fully known. We show how the boundary conditions must be taken into account, and prove convergence in L1,loc of the fully discrete splitting procedure to the very weak solution of the original system based on compactness arguments via total variation estimates. Generally, this is the best convergence obtained for this type of approximation. The derivation indicates limitations of the approach, being able to consider only some types of boundary conditions. A sample numerical experiment of a problem with an analytical solution is given, showing the stated efficiency of the method.     

基于效益最大化的城际铁路属性设置及调整

为便于城际铁路企业及时、合理地调整城际铁路运输属性,使其能在区域交通走廊各种交通工具竞争激烈的条件下健康运营和可持续性发展,分析了城际间各交通工具的属性特征,并对其进行了定量化描述。以城际铁路企业利润最大化为优化目标,考虑了城际铁路应取得最小社会效益的约束,建立了城际铁路运输属性设置及调整的多元非线性模型。利用罚函数法把模型转化为无约束的非线性模型,运用步长加速法进行求解。最后,把所设计的模型及算法应用到长株潭城际铁路初期运输属性设置中,结果表明本文所建立的模型及采用的算法是可行的,得到的城际铁路的技术经济特征设置值能为企业决策提供参考。     

Modeling and analysis of reactive solute transport in deformable channels with wall adsorption–desorption

We show well posedness for a model of nonlinear reactive transport of chemical in a deformable channel. The channel walls deform due to fluid–structure interaction between an unsteady flow of an incompressible, viscous fluid inside the channel and elastic channel walls. Chemical solutes, which are dissolved in the viscous, incompressible fluid, satisfy a convection–diffusion equation in the bulk fluid, while on the deforming walls, the solutes undergo nonlinear adsorption–desorption physico‐chemical reactions. The problem addresses scenarios that arise, for example, in studies of drug transport in blood vessels. We show the existence of a unique weak solution with solute concentrations that are non‐negative for all times. The analysis of the problem is carried out in the context of semi‐linear parabolic PDEs on moving domains. The arbitrary Lagrangian–Eulerian approach is used to address the domain movement, and the Galerkin method with the Picard–Lindelöf theorem is used to prove existence and uniqueness of approximate solutions. Energy estimates combined with the compactness arguments based on the Aubin–Lions lemma are used to prove convergence of the approximating sequences to the unique weak solution of the problem. It is shown that the solution satisfies the positivity property, that is, that the density of the solute remains non‐negative at all times, as long as the prescribed fluid domain motion is ‘reasonable’. This is the first well‐posedness result for reactive transport problems defined on moving domains of this type. Copyright © 2015 John Wiley & Sons, Ltd.     

Linearity-preserving flux correction and convergence acceleration for constrained Galerkin schemes

This paper is concerned with the development of general-purpose algebraic flux correction schemes for continuous (linear and multilinear) finite elements. In order to enforce the discrete maximum principle (DMP), we modify the standard Galerkin discretization of a scalar transport equation by adding diffusive and antidiffusive fluxes. The result is a nonlinear algebraic system satisfying the DMP constraint. An estimate based on variational gradient recovery leads to a linearity-preserving limiter for the difference between the function values at two neighboring nodes. A fully multidimensional version of this scheme is obtained by taking the sum of local bounds and constraining the total flux. This new approach to algebraic flux correction provides a unified treatment of stationary and time-dependent problems. Moreover, the same algorithm is used to limit convective fluxes, anisotropic diffusion operators, and the antidiffusive part of the consistent mass matrix.The nonlinear algebraic system associated with the constrained Galerkin scheme is solved using fixed-point defect correction or a nonlinear SSOR method. A dramatic improvement of nonlinear convergence rates is achieved with the technique known as Anderson acceleration (or Anderson mixing). It blends a number of last iterates in a GMRES fashion, which results in a Broyden-like quasi-Newton update. The numerical behavior of the proposed algorithms is illustrated by a grid convergence study for convection-dominated transport problems and anisotropic diffusion equations in 2D.     

An adaptive characteristics method for advective-diffusive transport

An adaptive characteristics method is presented for the solution of advective-diffusive groundwater transport problems. The method decomposes the transport processes into advective and diffusive transport components. Advective flows are defined by using a streamtube incrementing procedure, based on the method of characteristics, to define the position of advective front. Diffusive transport orthogonal to the front is represented by an array of propagating streamtube elements, with dimension determined from analytical solution of the one-dimensional diffusion equation. Adaptive time scaling is used to moderate the dimensions and aspect ratios of the advective and diffusive streamtube elements as appropriate to the dominant transport mechanism. Finite differences are used to solve for transport ahead of the advancing front. The distribution of streamtubes are predetermined from a direct boundary element algorithm. Comparison with analytical results for a one-dimensional transport geometry indicates agreement for Peclet numbers between zero and infinity. Solution for transport in two-dimensional domains illustrates excellent agreement for Peclet numbers from zero to 25.     

Solution of solute transport in unsaturated porous media by the method of characteristics

A degenerate convection‐diffusion problem is approximated using the scheme that is based on the relaxation method and also the method of characteristics. A mathematical model for solute transport in unsaturated porous media is included. Moreover, multiple site adsorption is considered. Convergence of the scheme is proved and numerical experiments in 1D and 2D are presented. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 732–761, 2003.     

Solution of the linear and nonlinear advection–diffusion problems on a sphere

Various linear advection–diffusion problems and nonlinear diffusion problems on a sphere are considered and solved using the direct, implicit and unconditionally stable finite-volume method of second-order approximation in space and time. In the absence of external forcing and dissipation, the method preserves the total mass of the substance and the norm of the solution. The component wise operator splitting allows us to develop the direct (noniterative) and fast numerical algorithm. The split problems in the longitudinal direction are solved using the Sherman-Morrison formula and Thomas algorithm. The direct solution of the split problems in the latitudinal direction requires the use of the bordering method for a block matrix, and the preliminary determination of the solution at the poles. The resulting systems with tridiagonal matrices are solved by the Thomas algorithm. The numerical experiments demonstrate that the method correctly describes the local advection–diffusion processes on the sphere, in particular, through the poles, and accurately simulate blow-up regimes (unlimited growing solutions) of nonlinear combustion, the propagation of nonlinear temperature and spiral waves, and solutions to Gray-Scott reaction–diffusion model.     

Differential quadrature solution of heat- and mass-transfer equations

The heat- and mass-transfer equations have an important role in various thermal and diffusion processes. These equations are nonlinear, due to the solution dependent diffusion coefficient and the source term. In this study, one- and two-dimensional nonlinear heat- and mass-transfer equations are solved numerically. To this end, the differential quadrature method is used to discretize the problem spatially and the resulting nonlinear system of ordinary differential equations in time are solved using the Runge–Kutta method. The solution is improved in time iteratively by solving considerably small sized linear system of resulting equations. To demonstrate its usefulness and accuracy, the proposed method is applied to four test problems, involving different nonlinearities.     

半导体问题的交替方向有限元方法

该文用交替方向有限元方法求解半导体问题的Energy Trans port (ET)模型。对模型中椭圆型的电子位势方程采用交替方向迭代法，对流占优扩散的电子浓度和空穴浓度方程采用特征交替方向有限元方法，热传导方程利用Patch逼近采用交替方向有限元方法求解。利用微分方程的先验估计理论和技巧，分别得到了椭圆型方程和抛物型方程的最优H+1和L+2误差估计。     

Phase space transport in a map of asteroid motion

IntroductionIt is well known that a complex phase space of a Hamiltonian system containing largemeasures of both regular and chaotic orbits is often partitioned by such partial obstructionas canton or Arnold web, ac.hich although not serving as absolute barriers, can significantlyimpede the motion of a chaotic orbit through a connected phase space region. This "stickiness"effect makes the phase space transport complicated. In fact, the chaotic transport or diffusionphenomenon can be met in ma…     

“运价”为负的运输模型计算机求解问题的处理

为解决一些计算机软件求解"运价"既有正值又有负值运输模型时"不可求解"的问题,本文采用"运价同额增减法"决策模型转换的方法,将原模型的"运价"全部转换为正值后再用计算机软件求解,并分别编写了EXCEL求解模板和求解程序对该方法的计算加以印证。结果表明,采用该方法求解得出的最优解(最优决策方案)与原模型求得的最优解完全一样,而最优值(最优决策效果)减去虚增(或加上虚减)的部分就是原模型的最优值。采用这种方法能成功地解决一些计算机软件"不可求解"的问题。     

A numerical scheme based on a solution of nonlinear advection–diffusion equations

A mathematical formulation of the two-dimensional Cole–Hopf transformation is investigated in detail. By making use of the Cole–Hopf transformation, a nonlinear two-dimensional unsteady advection–diffusion equation is transformed into a linear equation, and the transformed equation is solved by the spectral method previously proposed by one of the authors. Thus a solution to initial value problems of nonlinear two-dimensional unsteady advection–diffusion equations is derived. On the base of the solution, a numerical scheme explicit with respect to time is presented for nonlinear advection–diffusion equations. Numerical experiments show that the present scheme possesses the total variation diminishing properties and gives solutions with good quality.     

A degenerate elliptic-parabolic system arising in competitive contaminant transport

The objective of this work is to study a coupled system of degenerate and nonlinear partial differential equations governing the transport of reactive solutes in groundwater. We show that this system admits a unique weak solution provided the nonlinear adsorption isotherm associated with the reaction process satisfies certain physically reasonable structural conditions, by addressing a more general problem. In addition, we conclude, that the solute concentrations stay non-negative if the source term is componentwise non-negative and investigate numerically the finite speed of propagation of compactly supported initial concentrations, in a two-component test case.     

Eulerian–Lagrangian localized adjoint methods for reactive transport with biodegradation

The microbial degradation of organic contaminants in the subsurface holds significant potential as a mechanism for in-situ remediation strategies. The mathematical models that describe contaminant transport with biodegradation involve a set of advective–diffusive–reactive transport equations. These equations are coupled through the nonlinear reaction terms, which may involve reactions with all of the species and are themselves coupled to growth equations for the subsurface bacterial populations. In this article, we develop Eulerian–Lagrangian localized adjoint methods (ELLAM) to solve these transport equations. ELLAM are formulated to systematically adapt to the changing features of governing partial differential equations. The relative importance of retardation, advection, diffusion, and reaction is directly incorporated into the numerical method by judicious choice of the test functions that appear in the weak form of the governing equation. Different ELLAM schemes for linear variable–coefficient advective–diffusive–reactive transport equations are developed based on different operator splittings. Specific linearization techniques are discussed and are combined with the ELLAM schemes to solve the nonlinear, multispecies transport equations. © 1995 John Wiley & Sons, Inc.     

An unconditionally stable finite difference scheme for solving a 3D heat transport equation in a sub-microscale thin film

Heat transport at the microscale is of vital importance in microtechnology applications. The heat transport equation is different from the traditional heat diffusion equation since a second-order derivative of temperature with respect to time and a third-order mixed derivative of temperature with respect to space and time are introduced. In this study, we develop a finite difference scheme with two levels in time for the 3D heat transport equation in a sub-microscale thin film. It is shown by the discrete energy method that the scheme is unconditionally stable. The 3D implicit scheme is then solved by using a preconditioned Richardson iteration, so that only a tridiagonal linear system is solved for each iteration. The numerical procedure is employed to obtain the temperature rise in a gold sub-microscale thin film.     

基于动态BP神经网络的非线性预测及在港口货物吞吐量中的应用

文章给出了非线性预测的一般理论方法，运用动态ＢＰ神经网络对实际问题的历史数据进行学习，然后根据学习后获得的非线性机理来进行预测。并将此方法应用于港口货物吞吐量及出口量预报。仿真表明此方法是有效的。