Temporal Moments for Reactive Transport through Fractured Impermeable/Permeable Formations

The transport of reactive solutes through fractured porous formations has been analyzed. The transport through the porous block is represented by a general multiprocess nonequilibrium equation (MPNE), which, for the fracture, is represented by an advection-dispersion equation with linear equilibrium sorption and first-order transformation. An implicit finite-difference technique has been used to solve the two coupled equations. The transport characteristics have been analyzed in terms of zeroth, first, and second temporal moments of the solute in the fracture. The solute behavior for fractured impermeable and fractured permeable formations are first compared and the effects of various fracture and matrix transport parameters are analyzed. Subsequently, the transport through a fractured permeable formation is analyzed to ascertain the effect of equilibrium sorption, rate-limited sorption, and the multiprocess nonequilibrium transport process. It was found that the temporal moments were nearly identical for the fractured impermeable and permeable formations when both the diffusion coefficient and the first-order transformation coefficient were relatively large. The multiprocess nonequilibrium model resulted in a smaller mass recovery in the fracture and higher dispersion than the equilibrium and rate-limited sorption models. DOI: 10.1061/(ASCE)HE.19435584.0000586. (C) 2012 American Society of Civil Engineers.

Ion transport studies in lithium phospho-molybdate glasses containing Cl- ion

Ion conducting glasses in xLiCl-20Li(2)O-(80-x) 0.80P(2)O(5)-0.20MoO(3)] glass system have been prepared over a wide range of composition (X = 5, 10, 15, 20 and 25 mol%). The electrical conductivity and dielectric relaxation of these glasses were analyzed using impedance spectroscopy in the frequency range of 10 Hz-10 MHz and in the temperature range of 313-353 K. D.c. activation energies extracted from Arrhenius plots using regression analysis, decreases with increasing LiCl mol%. A.c. conductivity data has been fitted to both single and double power law equation with both fixed and variable parameters. The increased conductivity in the present glass system has been correlated with the volume increasing effect and the coordination changes that occur due to structural modification resulting in the creation of non-bridging oxygens (NBO's) of the type O-Mo-O- bonds in the glass network. Dielectric relaxation mechanism in these glasses is analyzed using Kohlrausch-Williams-Watts (KWW) stretched exponential function and stretched exponent (beta) is found to be insensitive to temperature.

An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

In this paper we present a massively parallel open source solver for Richards equation, named the RichardsFOAM solver. This solver has been developed in the framework of the open source generalist computational fluid dynamics tool box OpenFOAM (R) and is capable to deal with large scale problems in both space and time. The source code for RichardsFOAM may be downloaded from the CPC program library website. It exhibits good parallel performances (up to similar to 90% parallel efficiency with 1024 processors both in strong and weak scaling), and the conditions required for obtaining such performances are analysed and discussed. These performances enable the mechanistic modelling of water fluxes at the scale of experimental watersheds (up to few square kilometres of surface area), and on time scales of decades to a century. Such a solver can be useful in various applications, such as environmental engineering for long term transport of pollutants in soils, water engineering for assessing the impact of land settlement on water resources, or in the study of weathering processes on the watersheds. (C) 2014 Elsevier B.V. All rights reserved.

Recovery of the solitons using a lattice Boltzmann model

We formulate a lattice Boltzmann model which simulates Korteweg-de Vries equation by using a method of higher moments of lattice Boltzmann equation. Using a series of lattice Boltzmann equations in different time scales and the conservation law in time scale to, we obtain equilibrium distribution function. The numerical examples show that the method can be used to simulate soliton.

A New Model for Sediment Transport

The problem of predicting sediment transportation by water waves is treated analytically with the rate of wave energy dissipation or wave damping. With resorting to the theory of shallow water waves and the basis of Yamamoto’s Coulomb-damped poroelastic model, the Boussinesq-type equation has been derived over a variation depth bed. For convenience Cnoidal wave is just discussed, The Cnoidal wave with complex wave length and wave velocity, which are as a function of wave frequency, water depth, permeability, Poisson’s ratio and complex elastic moduli of bed soil, is applied to analyse the rate of sediment transportation. Considering the sediment transportation depended on the shear stress near-bed or the horizontal velocity, the conclusion of Yamamoto’s experiment in clay bed has been extended to general situation. It could be figured out that the model should provide a method to avoid the undistinguishable factors during sediment transport processes and relate mass transport with the sediment peculiarities.

A GIS-based model of soil erosion and transport

Soil erosion is a natural process that occurs when the force of wind, raindrops or running water on the soil surface exceeds the cohesive forces that bind the soil together. In general, vegetation cover protects the soil from the effects of these erosive forces. However, land management activities such as ploughing, burning or heavy grazing may disturb this protective layer, exposing the underlying soil. The decision making process in rural catchment management is often supported by the predictive modelling of soil erosion and sediment transport processes within the catchment, using established techniques such as the Universal Soil Loss Equation [USLE] and the Agricultural Nonpoint Source pollution model [AGNPS]. In this article, the authors examine the range of erosion models currently available and describe the application of one of these to the Burrishoole catchment on the north-west coast of Ireland, which has suffered heavy erosion of blanket peat in recent years.

On the transport properties of fluid-particle flow

The hydrodynamic forces acting on a solid particle in a viscous, incompressible fluid medium at low Reynolds number flow is investigated mathematically as a prerequisite to the understanding of transport processes in two-phase flow involving solid particles and fluid. Viscous interaction between a small number of spherical particles and continuous solid boundaries as well as fluid interface are analyzed under a “point-force” approximation. Non-spherical and elastic spherical particles in a simple shear flow area are then considered. Non-steady motion of a spherical particle is briefly touched upon to illustrate the transient effect of particle motion.

A macroscopic continuum description of particle-fluid flow is formulated in terms of spatial averages yielding a set of particle continuum and bulk fluid equations. Phenomenological formulas describing the transport processes in a fluid medium are extended to cases where the volume concentration of solid particles is sufficiently high to exert an important influence. Hydrodynamic forces acting on a spherical solid particle in such a system, e.g. drag, torque, rotational coupling force, and viscous collision force between streams of different sized particles moving relative to each other are obtained. Phenomenological constants, such as the shear viscosity coefficient, and the diffusion coefficient of the bulk fluid, are found as a function of the material properties of the constituents of the two-phase system and the volume concentration of solid. For transient heat conduction phenomena, it is found that the introduction of a complex conductivity for the bulk fluid permits a simple mathematical description of this otherwise complicated process. The rate of heat transfer between particle continuum and bulk fluid is also investigated by means of an Oseen-type approximation to the energy equation.

Validation of a lattice Boltzmann model for gas-solid reactions with experiments

A lattice Boltzmann method is used to model gas-solid reactions where the composition of both the gas and solid phase changes with time, while the boundary between phases remains fixed. The flow of the bulk gas phase is treated using a multiple relaxation time MRT D3Q19 model; the dilute reactant is treated as a passive scalar using a single relaxation time BGK D3Q7 model with distinct inter- and intraparticle diffusivities. A first-order reaction is incorporated by modifying the method of Sullivan et al. [13] to include the conversion of a solid reactant. The detailed computational model is able to capture the multiscale physics encountered in reactor systems. Specifically, the model reproduced steady state analytical solutions for the reaction of a porous catalyst sphere (pore scale) and empirical solutions for mass transfer to the surface of a sphere at Re=10 (particle scale). Excellent quantitative agreement between the model and experiments for the transient reduction of a single, porous sphere of Fe 2O 3 to Fe 3O 4 in CO at 1023K and 10 5Pa is demonstrated. Model solutions for the reduction of a packed bed of Fe 2O 3 (reactor scale) at identical conditions approached those of experiments after 25 s, but required prohibitively long processor times. The presented lattice Boltzmann model resolved successfully mass transport at the pore, particle and reactor scales and highlights the relevance of LB methods for modelling convection, diffusion and reaction physics. © 2012 Elsevier Inc.

Spin-dependent current noises in transport through coupled quantum dots

Current fluctuations can provide additional insight into quantum transport in mesoscopic systems. The present work is carried out for the fluctuation properties of transport through a pair of coupled quantum dots which are connected with ferromagnetic electrodes. Based on an efficient particle-number-resolved master equation approach, we are concerned with not only fluctuations of the total charge and spin currents, but also of each individual spin-dependent component. As a result of competition among the spin polarization, Coulomb interaction, and dot-dot tunnel coupling, rich behaviors are found for the self- and mutual-correlation functions of the spin-dependent currents.

Quantum transport from the perspective of quantum open systems

By viewing the non-equilibrium transport setup as a quantum open system, we propose a reduced-density-matrix based quantum transport formalism. At the level of self-consistent Born approximation, it can precisely account for the correlation between tunneling and the system internal many-body interaction, leading to certain novel behavior such as the non-equilibrium Kondo effect. It also opens a new way to construct time-dependent density functional theory for transport through large-scale complex systems. (c) 2006 Elsevier B.V. All rights reserved.

Electron transport through a double-quantum-dot structure with intradot and interdot Coulomb interactions

Electron transport through a double-quantum-dot structure with intradot and interdot Coulomb interactions is studied by a Green's function (GF) approach. The conductance is calculated by a Landauer-Buttiker formula for the interacting systems derived using the nonequilibrium Keldysh formalism and the GF's are solved by the equation-of-motion method. It is shown that the interdot-coupling dependence of the conductance peak splitting matches the recent experimental observations. Also, the breaking of the electron-hole symmetry is numerically demonstrated by the presence of the interdot repulsion. [S0163-1829(99)01640-9].

Numerical investigation of particle velocity distributions in aeolian sand transport

Particle velocity distribution in a blowing sand cloud is a reflection of saltation movement of many particles. Numerical analysis is performed for particle velocity distribution with a discrete particle model. The probability distributions of resultant particle velocity in the impact-entrainment process, particle horizontal and vertical velocities at different heights and the vertical velocity of ascending particles are analyzed. The probability distributions of resultant impact and lift-off velocities of saltating particles can be expressed by a log-normal function, and that of impact angle comply with an exponential function. The probability distribution of particle horizontal and vertical velocities at different heights shows a typical single-peak pattern. In the lower part of saltation layer, the particle horizontal velocity distribution is positively skewed. Further analysis shows that the probability density function of the vertical velocity of ascending particles is similar to the right-hand part of a normal distribution function, and a general equation is acquired for the probability density function of non-dimensional vertical velocity of ascending particles which is independent of diameter of saltating particles, wind strength and height. These distributions in the present numerical analysis are consistent with reported experimental results. The present investigation is important for understanding the saltation state in wind-blown sand movement. (C) 2009 Elsevier B.V. All rights reserved.

The equation of state study in UU collisions at CSR, Lanzhou

The cooling storage ring, to be built at Lanzhou, will be able to deliver heavy ion beams up to uranium up to 0.52 GeV/u. It is expected to make considerable contribution to nuclear EOS study in the high net baryon-density region. With a relativistic transport model, we performed simulations for U+U collisions with different orientations. It is shown that by combining the forward neutron multiplicity and an event-wise elliptic flow selection, it is possible to identify the tip - tip and body - body head-on collisions. The effective identification of these two extreme configurations will allow us to study the EOS at the highest baryon density in the U+U collisions.

Effects of elastic and inelastic NN scattering cross sections on pi(-)/pi(+) ratios in heavy-ion collisions at intermediate energies

Based on the isospin-dependent Boltzmann-Uehling-Uhlenbeck transport model and the scaling model according to nucleon effective mass, effects of elastic and inelastic NN scattering cross sections on pi(-)/pi(+) in the neutron-rich reaction of Ca-48 + Ca-48 at a beam energy of 400 MeV/nucleon are studied. It is found that cross-section effects of both NN elastic and inelastic scatterings affect Delta(1232), pi(-) and pi(+) production, as well as the value of pi(-)/pi(+).

Can the nuclear symmetry potential at supra-saturation densities be negative?

In the framework of an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model, for the central Au-197 + Au-197 reaction at an incident beam energy of 400 MeV/nucleon, the effect of nuclear symmetry potential at supra-saturation densities on the preequilibrium clusters emission is studied. It is found that for the positive symmetry potential at supra-saturation densities the neutron-to-proton ratio of lighter clusters with mass number A less than or similar to 3 [(n/p)(A less than or similar to 3)] is larger than that of the heavier clusters with mass number A > 3 [(n/p)(A>3)], whereas for the negative symmetry potential at supra-saturation densities the (n/p)(A less than or similar to 3) is smaller than the (n/p)(A>3). This may be considered as a probe of the negative symmetry potential at supra-saturation densities.