Modeling of seepage flow through layered solis

A two-dimensional finite difference model, which solves mixed type of Richards' equation, whose non-linearity is dealt with modified Picard's iteration and strongly implicit procedure to solve the resulting equations, is presented. Modeling of seepage flow through heterogeneous soils, which is common in the field is addressed in the present study. The present model can be applied to both unsaturated and saturated soils and can handle very dry initial condition and steep wetting fronts. The model is validated by comparing experimental results reported in the literature. Newness of this two dimensional model is its application on layered soils with transient seepage face development, which has not been reported in the literature. Application of the two dimensional model for studying unconfined drainage due to sudden drop of water table at seepage face in layered soils is demonstrated. In the present work different sizes of rectangular flow domain with different types of layering are chosen. Sensitivity of seepage height due to problem dimension of layered system is studied. The effect of aspect ratio on seepage face development in case of the flow through layered soil media is demonstrated. The model is also applied to random heterogeneous soils in which the randomness of the model parameters is generated using the turning band technique. The results are discussed in terms of phreatic surface and seepage height development and also flux across the seepage face. Such accurate modeling of seepage face development and quantification of flux moving across the seepage face becomes important while modeling transport problems in variably saturated media.

Effect of convection and viscoelasticity on spinodal decomposition in fluids

Current analytical work on the effect of convection and viscoelasticity on the early and late stages of spinodal decomposition is briefly described. In the early stages, the effect of viscoelastic stresses was analysed using a simple Maxwell model for the stress, which was incorporated in the Langevin equation for the momentum field. The viscoelastic stresses are found to enhance the rate of decomposition. In the late stages, the pattern formed depends on the relative composition of the two species. Droplet spinodal decomposition occurs when the concentration of one of the species is small. Convective transport does not have a significant effect on the growth of a single droplet, but it does result in an attractive interaction between non - Brownian droplets which could lead to coalescence. The effect of convective transport for the growth of random interfaces in a near symmetric quench was analysed using an 'area distribution function', which gives the distribution of surface area of the interface in curvature space. It was found that the curvature of the interface decreases proportional to t in the late stages of spinodal decomposition, and the surface area also decreases proportional to t.

Application of kinetic schemes to all types of meshes

Kinetic schemes as pursued in CFD Centre are obtained by taking suitable moments of upwind schemes for Boltzmann equation without collision term. The primary ones among these are KFVS, LSKUM, KFMG and these have been applied successfully to a variety of flow problems using various meshes. These schemes have been found to be very robust.

Universalities in vortex transport at the melting transition

We consider the jump in resistance at the melting transition, which is experimentally observed to be constant, independent of magnetic field (vortex density). We present an explanation of this effect based on vortex cuttings, and universalities of the structure factor at the freezing transition (the Hansen-Verlet criterion).

Transport through an electrostatically defined quantum dot lattice in a two-dimensional electron gas

Quantum dot lattices (QDLs) have the potential to allow for the tailoring of optical, magnetic, and electronic properties of a user-defined artificial solid. We use a dual gated device structure to controllably tune the potential landscape in a GaAs/AlGaAs two-dimensional electron gas, thereby enabling the formation of a periodic QDL. The current-voltage characteristics, I (V), follow a power law, as expected for a QDL. In addition, a systematic study of the scaling behavior of I (V) allows us to probe the effects of background disorder on transport through the QDL. Our results are particularly important for semiconductor-based QDL architectures which aim to probe collective phenomena.

An inverse problem for Helmholtz equation

This article is concerned with subsurface material identification for the 2-D Helmholtz equation. The algorithm is iterative in nature. It assumes an initial guess for the unknown function and obtains corrections to the guessed value. It linearizes the otherwise nonlinear problem around the background field. The background field is the field variable generated using the guessed value of the unknown function at each iteration. Numerical results indicate that the algorithm can recover a close estimate of the unknown function based on the measurements collected at the boundary.

Ion Transport in Liquid Salt Solutions with Oxide Dispersions:``Soggy Sand'' Electrolytes

``Soggy sand'' electrolyte, which essentially consists of oxide dispersions in nonaqueous liquid salt solutions, comprises an important class of soft matter electrolytes. The ion transport mechanism of soggy sand electrolyte is complex. The configuration of particles in the liquid solution has been observed to depend in a nontrivial manner on various parameters related to the oxide (concentration, size, surface chemistry) and solvent (dielectric constant, viscosity) as well as time. The state of the particles in solution not only affects ionic conductivity but also effectively the mechanical and electrochemical properties of the solid liquid composite. Apart from comprehensive understanding of the underlying phenomena that govern ion transport, which will benefit design of better electrolytes, the problem has far-reaching implications in diverse fields such as catalysis, colloid chemistry, and biotechnology.

On the Presumed Kinetic Consequences of Pre-equilibrium. Implications for the Michaelis-Menten Equation

The overall rate equation for a reaction sequence consisting of a pre-equilibrium and rate-determining steps should not be derived on the basis of the concentration of the intermediate product (X). This is apparently indicated by transition state theory (as the path followed to reach the highest energy transition state is irrelevant), but also proved by a straight-forward mathematical approach. The thesis is further supported by the equations of concurrent reactions as applied to the partitioning of X between the two competing routes (reversal of the pre-equilibrium and formation of product). The rate equation may only be derived rigorously on the basis of the law of mass action. It is proposed that the reactants acquire the overall activation energy prior to the pre-equilibrium, thus forming X in a high-energy state en route to the rate-determining transition state. (It is argued that conventional energy profile diagrams are misleading and need to be reinterpreted.) Also, these arguments invalidate the Michaelis-Menten equation of enzyme kinetics, and necessitate a fundamental revision of our present understanding of enzyme catalysis. (The observed ``saturation kinetics'' possibly arises from weak binding of a second molecule of substrate at the active site; analogous conclusions apply to reactions at surfaces).

Carrier-transport studies of III-nitride/Si3N4/Si isotype heterojunctions

GaN/Si3N4/n-Si and InN/Si3N4/n-Si heterojunctions (HJs) were fabricated using plasma-assisted molecular beam epitaxy for a comparison study. Single-crystalline wurtzite structures of GaN and InN epilayers were confirmed by high-resolution X-ray diffraction and thickness of ultrathin Si3N4 layer was measured by transmission electron microscopy. n-GaN/Si3N4/n-Si HJs show diode-like rectifying current-voltage (I-V) characteristic, while n-InN/Si3N4/n-Si HJs show symmetric nonlinear I-V behavior. The I-V characteristics of both HJs were discussed in terms of the band diagram of HJs and the carrier transport mechanism. The activation energies of carrier conduction were estimated to be similar to 29 meV for GaN/Si3N4/Si and similar to 95 meV for InN/Si3N4/Si HJs. (C) 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Quenching of Meridional Circulation in Flux Transport Dynamo Models

Guided by the recent observational result that the meridional circulation of the Sun becomes weaker at the time of the sunspot maximum, we have included a parametric quenching of the meridional circulation in solar dynamo models such that the meridional circulation becomes weaker when the magnetic field at the base of the convection zone is stronger. We find that a flux transport solar dynamo tends to become unstable on including this quenching of meridional circulation if the diffusivity in the convection zone is less than about 2x10(11) cm(2) s(-1). The quenching of alpha, however, has a stabilizing effect and it is possible to stabilize a dynamo with low diffusivity with sufficiently strong alpha-quenching. For dynamo models with high diffusivity, the quenching of meridional circulation does not produce a large effect and the dynamo remains stable. We present a solar-like solution from a dynamo model with diffusivity 2.8x10(12) cm(2) s(-1) in which the quenching of meridional circulation makes the meridional circulation vary periodically with solar cycle as observed and does not have any other significant effect on the dynamo.

Transport in nanoporous zeolites: Relationships between sorbate size, entropy, and diffusivity

Molecular dynamics simulations have been performed on monatomic sorbates confined within zeolite NaY to obtain the dependence of entropy and self-diffusivity on the sorbate diameter. Previously, molecular dynamics simulations by Santikary and Yashonath J. Phys. Chem. 98, 6368 (1994)], theoretical analysis by Derouane J. Catal. 110, 58 (1988)] as well as experiments by Kemball Adv. Catal. 2, 233 (1950)] found that certain sorbates in certain adsorbents exhibit unusually high self-diffusivity. Experiments showed that the loss of entropy for certain sorbates in specific adsorbents was minimum. Kemball suggested that such sorbates will have high self-diffusivity in these adsorbents. Entropy of the adsorbed phase has been evaluated from the trajectory information by two alternative methods: two-phase and multiparticle expansion. The results show that anomalous maximum in entropy is also seen as a function of the sorbate diameter. Further, the experimental observation of Kemball that minimum loss of entropy is associated with maximum in self-diffusivity is found to be true for the system studied here. A suitably scaled dimensionless self-diffusivity shows an exponential dependence on the excess entropy of the adsorbed phase, analogous to excess entropy scaling rules seen in many bulk and confined fluids. The two trajectory-based estimators for the entropy show good semiquantitative agreement and provide some interesting microscopic insights into entropy changes associated with confinement.

Low temperature charge transport and microwave absorption of carbon coated iron nanoparticles-polymer composite films

In this paper, the low temperature electrical conductivity and microwave absorption properties of carbon coated iron nanoparticles-polyvinyl chloride composite films are investigated for different filler fractions. The filler particles are prepared by the pyrolysis of ferrocene at 980 degrees C and embedded in polyvinyl chloride matrix. The high resolution transmission electron micrographs of the filler material have shown a 5 nm thin layer graphitic carbon covering over iron particles. The room temperature electrical conductivity of the composite film changes by 10 orders of magnitude with the increase of filler concentration. A percolation threshold of 2.2 and an electromagnetic interference shielding efficiency (EMI SE) of similar to 18.6 dB in 26.5-40 GHz range are observed for 50 wt% loading. The charge transport follows three dimensional variable range hopping conduction. (C) 2012 Elsevier Ltd. All rights reserved.

Macroscopic Transport of Mega-ampere Electron Currents in Aligned Carbon-Nanotube Arrays

We demonstrate that aligned carbon-nanotube arrays are efficient transporters of laser-generated megaampere electron currents over distances as large as a millimeter. A direct polarimetric measurement of the temporal and the spatial evolution of the megagauss magnetic fields (as high as 120 MG) at the target rear at an intensity of (10(18)-10(19)) W/cm(2) was corroborated by the rear-side hot electron spectra. Simulations show that such high magnetic flux densities can only be generated by a very well collimated fast electron bunch.

Linearization and Reconstruction of Non-linear Diffuse Optical Tomographic Image

Diffuse optical tomography (DOT) is one of the ways to probe highly scattering media such as tissue using low-energy near infra-red light (NIR) to reconstruct a map of the optical property distribution. The interaction of the photons in biological tissue is a non-linear process and the phton transport through the tissue is modelled using diffusion theory. The inversion problem is often solved through iterative methods based on nonlinear optimization for the minimization of a data-model misfit function. The solution of the non-linear problem can be improved by modeling and optimizing the cost functional. The cost functional is f(x) = x(T)Ax - b(T)x + c and after minimization, the cost functional reduces to Ax = b. The spatial distribution of optical parameter can be obtained by solving the above equation iteratively for x. As the problem is non-linear, ill-posed and ill-conditioned, there will be an error or correction term for x at each iteration. A linearization strategy is proposed for the solution of the nonlinear ill-posed inverse problem by linear combination of system matrix and error in solution. By propagating the error (e) information (obtained from previous iteration) to the minimization function f(x), we can rewrite the minimization function as f(x; e) = (x + e)(T) A(x + e) - b(T)(x + e) + c. The revised cost functional is f(x; e) = f(x) + e(T)Ae. The self guided spatial weighted prior (e(T)Ae) error (e, error in estimating x) information along the principal nodes facilitates a well resolved dominant solution over the region of interest. The local minimization reduces the spreading of inclusion and removes the side lobes, thereby improving the contrast, localization and resolution of reconstructed image which has not been possible with conventional linear and regularization algorithm.

Electrical Transport Properties of Buckled Carbon Nanotube Arrays

We report on the electrical transport properties of buckled carbon nanotube arrays synthesized by pyrolysis. Analyzing the experimental data based on the general theory of semiconductors, the arrays are predicted to be semiconducting and the band gap can be evaluated. The band gap of different arrays is in 25-50 meV range.