Nonequilibrium effects in models of three-phase flow in porous media

In this paper we extend to three-phase flow the nonequilibrium formalism proposed by Barenblatt and co-workers for two-phase porous media flow. The underlying idea is to include nonequilibrium effects by introducing a pair of effective water and gas saturations, which are linked to the actual saturations by a local evolution equation. We illustrate and analyze how nonequilibrium effects lead to qualitative and quantitative differences in the solution of the three-phase flow equations.     

An accurate time integration method for simplified overland flow models

An accurate time integration method for the diffusion-wave and kinematic-wave approximated models for the overland flow is proposed. The discretization of the first- and second-order spatial derivatives in the basic equation is obtained by using the second-order Lax–Wendroff and the three-point centred finite difference schemes, respectively. For the solution in time, the system of ordinary differential equations, obtained by the linearization of the celerity and of the hydraulic diffusivity by Taylor series expansions, is integrated analytically. The stability and the numerical dissipation and dispersion are investigated by the Fourier analysis. A proper Courant number, and the corresponding time step for the numerical simulations can be established. In addition, the proposed diffusion-wave and kinematic-wave models are straightforwardly extended to the two-dimensional flow. Test cases for both one- and two-dimensional problems, compare the solutions of the diffusion-wave and kinematic-wave models with analytical solutions, with experimental results and with numerical solutions obtained by the Saint–Venant equations. These simulations show that the proposed numerical–analytical models accurately predict the overland flow for several situations, in particular for unsteady rainfall rate and for spatial variations of the surface roughness.     

Comparison of alternative upscaled model formulations for simulating DNAPL source dissolution and biodecay

This paper compares the performance of analytical and numerical approaches for modeling DNAPL dissolution with biodecay. A solution derived from a 1-D advective transport formulation (“Parker” model) is shown to agree very closely with high resolution numerical solutions. A simple lumped source mass balance solution in which with decay is assumed proportional to DNAPL mass (“Falta1” model) over- or underpredicts aqueous phase biodecay depending on the magnitude of the exponential factor governing the relationship between dissolution rate and DNAPL mass. A modification of the Falta model that assumes decay proportional to the source exit concentration is capable of accurately simulating source behavior with strong aqueous phase biodecay if model parameters are appropriately selected or calibrated (“Falta2” model). However, parameters in the lumped models exhibit complex interdependencies that cannot be quantified without consideration of transport processes within the source zone. Combining the Falta2 solution with relationships derived from the Parker model was found to resolve these limitations and track the numerical model results. A method is presented to generalize the analytical solutions to enable simulation of partial mass removal with changes in source parameters over time due to various remedial actions. The algorithm is verified by comparison with numerical simulation results. An example application is presented that demonstrates the interactions of partial mass removal, enhanced biodecay, enhanced mass transfer and source zone flow reduction applied at various time periods on contaminant flux reduction. Increasing errors that arise in numerical solutions with coarse discretization and high decay rates are shown to be controlled by using an adjusted decay coefficient derived from the Parker analytical solution.     

贵州遵义下寒武统黑色页岩镍、钼多金属矿床的形态硫特征及成矿模式

黄家湾矿是贵州遵义地区寒武系底部新发现的黑色页岩镍钼多金属元素矿床。根据形态硫测试结果,文章分析了黄家湾下寒武统牛蹄塘组黑色页岩剖面和矿石中不同形态硫（硫酸盐硫、黄铁矿硫、有机硫）的特征,推测镍、钼金属富集层处还原性硫的来源。     

Mathematical Analysis of a Model of River Channel Formation

The study of overland flow of water over an erodible sediment leads to a coupled model describing the evolution of the topographic elevation and the depth of the overland water film. The spatially uniform solution of this model is unstable, and this instability corresponds to the formation of rills, which in reality then grow and coalesce to form large-scale river channels. In this paper we consider the deduction and mathematical analysis of a deterministic model describing river channel formation and the evolution of its depth. The model involves a degenerate nonlinear parabolic equation (satisfied on the interior of the support of the solution) with a super-linear source term and a prescribed constant mass. We propose here a global formulation of the problem (formulated in the whole space, beyond the support of the solution) which allows us to show the existence of a solution and leads to a suitable numerical scheme for its approximation. A particular novelty of the model is that the evolving channel self-determines its own width, without the need to pose any extra conditions at the channel margin.     

A stochastic ground motion model with separable temporal and spectral nonstationarities

A fully nonstationary stochastic model for strong earthquake ground motion is developed. The model employs filtering of a discretized white‐noise process. Nonstationarity is achieved by modulating the intensity and varying the filter properties in time. The formulation has the important advantage of separating the temporal and spectral nonstationary characteristics of the process, thereby allowing flexibility and ease in modeling and parameter estimation. The model is fitted to target ground motions by matching a set of statistical characteristics, including the mean‐square intensity, the cumulative mean number of zero‐level up‐crossings and a measure of the bandwidth, all expressed as functions of time. Post‐processing by a second filter assures zero residual velocity and displacement, and improves the match to response spectral ordinates for long periods. Copyright © 2008 John Wiley & Sons, Ltd.     

A review of research on seismic behaviour of irregular building structures since 2002

The paper presents an overview of the progress in research regarding seismic response of plan and vertically irregular building structures. Three areas of research are surveyed. The first is the study of the effects of plan-irregularity by means of single-storey and multi-storey building models. The second area encompasses passive control as a strategy to mitigate torsional effects, by means of base isolation and other types of devices. Lastly, the third area concerns vertically irregular structures and setback buildings. Although fewer papers have been published in this last area with respect to the former ones, this state-of-the-art reports extensively on research efforts and progress into the seismic behaviour of irregular buildings in elevation to show the growing interest among specialists in the field.     

Adsorption kinetics, isotherm, and thermodynamic studies of adsorption of pollutant from aqueous solutions onto humic acid

In the present study,humic acid was used as an adsorbent for the investigation of the adsorption kinetics,isotherms,and thermo-dynamic parameters of hexavalent chromium from aqueous solution at varying pH,temperatures,and concentrations.Adsorption isotherms and equilibrium adsorption capacities were determined by the fittings of the experimental data to three well-known iso-therm models:Langmuir,Freundlich,and Redlich-Peterson.The results showed that the Langmuir and Redlich-Peterson models appear to fit the adsorption better than did the Freundlich adsorption model for the adsorption of chromium onto humic acid.The equilibrium constants were used to calculate thermodynamic parameters such as the change of free energy,enthalpy,and entropy.The derived adsorption constants (logaL) and their temperature dependencies from Langmuir isotherm have been used to calculate the corresponding thermodynamic quantities such as the free energy of adsorption,heat,and entropy of adsorption.The thermo-dynamic data indicate that Cr (VI) adsorption onto humic acid is entropically driven and characterized by physical adsorption.     

Hydrodynamics and sediment dynamics in the German Bight. A focus on observations and numerical modelling in the East Frisian Wadden Sea

This work deals with analysis of hydrographic observations and results of numerical simulations. The data base includes acoustic Doppler current profilers (ADCP) observations, continuous measurements on data stations and satellite data originating from the medium resolution imaging spectrometer (MERIS) onboard the European Space Agency (ESA) satellite ENVISAT with a spatial resolution of 300 m. Numerical simulations use nested models with horizontal resolutions ranging from 1 km in the German Bight to 200 m in the East Frisian Wadden Sea coupled with a suspended matter transport model. Modern satellite observations have now a comparable horizontal resolution with high-resolution numerical model of the entire area of the East Frisian Wadden Sea allowing to describe and validate new and so far unknown patterns of sediment distribution. The two data sets are consistent and reveal an oscillatory behaviour of sediment pools to the north of the back-barrier basins and clear propagation patterns of tidally driven suspended particulate matter outflow into the North Sea. The good agreement between observations and simulations is convincing evidence that the model simulates the basic dynamics and sediment transport processes, which motivates its further use in hindcasting, as well as in the initial steps towards forecasting circulation and sediment dynamics in the coastal zone.     

Spatial analysis of instream nitrogen loads and factors controlling nitrogen delivery to streams in the southeastern United States using spatially referenced regression on watershed attributes (SPARROW) and regional classification frameworks

Understanding how nitrogen transport across the landscape varies with landscape characteristics is important for developing sound nitrogen management policies. We used a spatially referenced regression analysis (SPARROW) to examine landscape characteristics influencing delivery of nitrogen from sources in a watershed to stream channels. Modelled landscape delivery ratio varies widely (by a factor of 4) among watersheds in the southeastern United States—higher in the western part (Tennessee, Alabama, and Mississippi) than in the eastern part, and the average value for the region is lower compared to other parts of the nation. When we model landscape delivery ratio as a continuous function of local‐scale landscape characteristics, we estimate a spatial pattern that varies as a function of soil and climate characteristics but exhibits spatial structure in residuals (observed load minus predicted load). The spatial pattern of modelled landscape delivery ratio and the spatial pattern of residuals coincide spatially with Level III ecoregions and also with hydrologic landscape regions. Subsequent incorporation into the model of these frameworks as regional scale variables improves estimation of landscape delivery ratio, evidenced by reduced spatial bias in residuals, and suggests that cross‐scale processes affect nitrogen attenuation on the landscape. The model‐fitted coefficient values are logically consistent with the hypothesis that broad‐scale classifications of hydrologic response help to explain differential rates of nitrogen attenuation, controlling for local‐scale landscape characteristics. Negative model coefficients for hydrologic landscape regions where the primary flow path is shallow ground water suggest that a lower fraction of nitrogen mass will be delivered to streams; this relation is reversed for regions where the primary flow path is overland flow. Published in 2009 by John Wiley & Sons, Ltd.