Physically consistent simulation of mesoscale chemical kinetics: The non-negative FIS-alpha method

Biochemical pathways involving chemical kinetics in medium concentrations (i.e., at mesoscale) of the reacting molecules can be approximated as chemical Langevin equations (CLE) systems. We address the physically consistent non-negative simulation of the CLE sample paths as well as the issue of non-Lipschitz diffusion coefficients when a species approaches depletion and any stiffness due to faster reactions. The non-negative Fully Implicit Stochastic alpha (FIS alpha) method in which stopped reaction channels due to depleted reactants are deleted until a reactant concentration rises again, for non-negativity preservation and in which a positive definite Jacobian is maintained to deal with possible stiffness, is proposed and analysed. The method is illustrated with the computation of active Protein Kinase C response in the Protein Kinase C pathway. (C) 2011 Elsevier Inc. All rights reserved.

Physically Based, Hydrologic Model Results Based on Three Precipitation Products

The main objective of the study is to examine the accuracy of and differences among simulated streamflows driven by rainfall estimates from a network of 22 rain gauges spread over a 2,170 km2 watershed, NEXRAD Stage III radar data, and Tropical Rainfall Measuring Mission (TRMM) 3B42 satellite data. The Gridded Surface Subsurface Hydrologic Analysis (GSSHA), a physically based, distributed parameter, grid-structured, hydrologic model, was used to simulate the June-2002 flooding event in the Upper Guadalupe River watershed in south central Texas. There were significant differences between the rainfall fields estimated by the three types of measurement technologies. These differences resulted in even larger differences in the simulated hydrologic response of the watershed. In general, simulations driven by radar rainfall yielded better results than those driven by satellite or rain-gauge estimates. This study also presents an overview of effects of land cover changes on runoff and stream discharge. The results demonstrate that, for major rainfall events similar to the 2002 event, the effect of urbanization on the watershed in the past two decades would not have made any significant effect on the hydrologic response. The effect of urbanization on the hydrologic response increases as the size of the rainfall event decreases.

Modeling of Channel Potential and Subthreshold Slope of Symmetric Double-Gate Transistor

A new physically based classical continuous potential distribution model, particularly considering the channel center, is proposed for a short-channel undoped body symmetrical double-gate transistor. It involves a novel technique for solving the 2-D nonlinear Poisson's equation in a rectangular coordinate system, which makes the model valid from weak to strong inversion regimes and from the channel center to the surface. We demonstrated, using the proposed model, that the channel potential versus gate voltage characteristics for the devices having equal channel lengths but different thicknesses pass through a single common point (termed ``crossover point''). Based on the potential model, a new compact model for the subthreshold swing is formulated. It is shown that for the devices having very high short-channel effects (SCE), the effective subthreshold slope factor is mainly dictated by the potential close to the channel center rather than the surface. SCEs and drain-induced barrier lowering are also assessed using the proposed model and validated against a professional numerical device simulator.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

A Jacobian-Based Algorithm for Planning Attitude Maneuvers Using Forward and Reverse Rotations

Algorithms for planning quasistatic attitude maneuvers based on the Jacobian of the forward kinematic mapping of fully-reversed (FR) sequences of rotations are proposed in this paper. An FR sequence of rotations is a series of finite rotations that consists of initial rotations about the axes of a body-fixed coordinate frame and subsequent rotations that undo these initial rotations. Unlike the Jacobian of conventional systems such as a robot manipulator, the Jacobian of the system manipulated through FR rotations is a null matrix at the identity, which leads to a total breakdown of the traditional Jacobian formulation. Therefore, the Jacobian algorithm is reformulated and implemented so as to synthesize an FR sequence for a desired rotational displacement. The Jacobian-based algorithm presented in this paper identifies particular six-rotation FR sequences that synthesize desired orientations. We developed the single-step and the multiple-step Jacobian methods to accomplish a given task using six-rotation FR sequences. The single-step Jacobian method identifies a specific FR sequence for a given desired orientation and the multiple-step Jacobian algorithm synthesizes physically feasible FR rotations on an optimal path. A comparison with existing algorithms verifies the fast convergence ability of the Jacobian-based algorithm. Unlike closed-form solutions to the inverse kinematics problem, the Jacobian-based algorithm determines the most efficient FR sequence that yields a desired rotational displacement through a simple and inexpensive numerical calculation. The procedure presented here is useful for those motion planning problems wherein the Jacobian is singular or null.

ESR Studies of Phase Transitions in Double Propionates: Dicalcium Barium Propionate Ca2Ba(C2H5COO)6

Single crystal [(111) and (100) planes], and powder ESR of Mn2+ (substituting for Ca2+) in Ca2Ba(C2H5COO)6 in the temperature range 220°C to -160°C shows (i) a doubling of both the physically and chemically inequivalent sites, and a change in the magnitude (150 G at -6°C to 170 G at -8°C) as well as the orientation of the D tensor across the -6°C transition and (ii) an inflection point in the D vs T plot across the -75°C transition. The oxygen octahedra around the Ca2+ sites are inferred to undergo alternate rotations, showing the participation of the carboxyl oxygens in the -6°C transition. A relation of the type D=D0(1+αT+βT2) seems to fit the D variation satisfactorily.

Classical Langevin dynamics of a charged particle moving on a sphere and diamagnetism: A surprise

It is generally known that the orbital diamagnetism of a classical system of charged particles in thermal equilibrium is identically zero —the Bohr-van Leeuwen theorem. Physically, this null result derives from the exact cancellation of the orbital diamagnetic moment associated with the complete cyclotron orbits of the charged particles by the paramagnetic moment subtended by the incomplete orbits skipping the boundary in the opposite sense. Motivated by this crucial but subtle role of the boundary, we have simulated here the case of a finite but unbounded system, namely that of a charged particle moving on the surface of a sphere in the presence of an externally applied uniform magnetic field. Following a real space-time approach based on the classical Langevin equation, we have computed the orbital magnetic moment that now indeed turns out to be non-zero and has the diamagnetic sign. To the best of our knowledge, this is the first report of the possibility of finite classical diamagnetism in principle, and it is due to the avoided cancellation.

Quantum-Ohmic Resistance Fluctuation In Disordered Conductors - An Invariant Imbedding Approach

It is now well known that in extreme quantum limit, dominated by the elastic impurity scattering and the concomitant quantum interference, the zero-temperature d.c. resistance of a strictly one-dimensional disordered system is non-additive and non-self-averaging. While these statistical fluctuations may persist in the case of a physically thin wire, they are implicitly and questionably ignored in higher dimensions. In this work, we have re-examined this question. Following an invariant imbedding formulation, we first derive a stochastic differential equation for the complex amplitude reflection coefficient and hence obtain a Fokker-Planck equation for the full probability distribution of resistance for a one-dimensional continuum with a Gaussian white-noise random potential. We then employ the Migdal-Kadanoff type bond moving procedure and derive the d-dimensional generalization of the above probability distribution, or rather the associated cumulant function –‘the free energy’. For d=3, our analysis shows that the dispersion dominates the mobilitly edge phenomena in that (i) a one-parameter B-function depending on the mean conductance only does not exist, (ii) an approximate treatment gives a diffusion-correction involving the second cumulant. It is, however, not clear whether the fluctuations can render the transition at the mobility edge ‘first-order’. We also report some analytical results for the case of the one dimensional system in the presence of a finite electric fiekl. We find a cross-over from the exponential to the power-low length dependence of resistance as the field increases from zero. Also, the distribution of resistance saturates asymptotically to a poissonian form. Most of our analytical results are supported by the recent numerical simulation work reported by some authors.

ESTIMATION OF EVAPOTRANSPIRATION IN AN EXPERIMENTAL WATERSHED IN A DECIDUOUS FOREST OF SOUTH INDIA

Forested areas play a dominant role in the global hydrological cycle. Evapotranspiration is a dominant component most of the time catching up with the rainfall. Though there are sophisticated methods which are available for its estimation, a simple reliable tool is needed so that a good budgeting could be made. Studies have established that evapotranspiration in forested areas is much higher than in agricultural areas. Latitude, type of forests, climate and geological characteristics also add to the complexity of its estimation. Few studies have compared different methods of evapotranspiration on forested watersheds in semi arid tropical forests. In this paper a comparative study of different methods of estimation of evapotranspiration is made with reference to the actual measurements made using all parameter climatological station data of a small deciduous forested watershed of Mulehole (area of 4.5 km2 ), South India. Potential evapotranspiration (ETo) was calculated using ten physically based and empirical methods. Actual evapotranspiration (AET) has been calculated through computation of water balance through SWAT model. The Penman-Montieth method has been used as a benchmark to compare the estimates arrived at using various methods. The AET calculated shows good agreement with the curve for evapotranspiration for forests worldwide. Error estimates have been made with respect to Penman-Montieth method. This study could give an idea of the errors involved whenever methods with limited data are used and also show the use indirect methods in estimation of Evapotranspiration which is more suitable for regional scale studies.

Evaluation of Charges for Power Transmission and Losses in Bilateral Power Contracts

The paper presents a method for transmission loss charge allocation in deregulated power systems based on Relative Electrical Distance (RED) concept. Based on RED between the generator and load nodes and the predefined bilateral power contracts, charge evaluation is carried out. Generally through some power exchange mechanism a set of bilateral contracts are determined that facilitate bilateral agreements between the generation and distribution entities. In this paper the possible charges incurred in meeting loads like generation charge, transmission charge and charge due to losses are evaluated. Case studies have been carried out on a few practical equivalent systems. Due to space limitation results for a sample 5 bus system are presented considering ideal load/generation power contracts and deviated load/generation power contracts. Extensive numerical testing indicates that the proposed allocation scheme produces loss allocations that are appropriate and that behave in a physically reasonable manner.

Rotating Beams and Nonrotating Beams With Shared Eigenpair

In this paper, we look for rotating beams whose eigenpair (frequency and mode-shape) is the same as that of uniform nonrotating beams for a particular mode. It is found that, for any given mode, there exist flexural stiffness functions (FSFs) for which the jth mode eigenpair of a rotating beam, with uniform mass distribution, is identical to that of a corresponding nonrotating uniform beam with the same length and mass distribution. By putting the derived FSF in the finite element analysis of a rotating cantilever beam, the frequencies and mode-shapes of a nonrotating cantilever beam are obtained. For the first mode, a physically feasible equivalent rotating beam exists, but for higher modes, the flexural stiffness has internal singularities. Strategies for addressing the singularities in the FSF for finite element analysis are provided. The proposed functions can be used as test-functions for rotating beam codes and for targeted destiffening of rotating beams.

Gibbs energies of formation of rare earth oxysulfides

The standard Gibbs energy change accompanying the conversion of rare earth oxides to oxysulfides by reaction of rare earth oxides with diatomic sulfur gas has been measured in the temperature range 870 to 1300 K using the solid state cell: Pt/Cu+Cu2S/R2O2S+R2O3‖(CaO)ZrO2‖Ni+NiO, Pt where R=La, Nd, Sm, Gd, Tb, and Dy. The partial pressure of diatomic sulfur over a mixture of rare earth oxide (R2O3) and oxysulfide (R2O2S) is fixed by the dissociation of Cu2S to Cu in a closed system. The buffer mixture of Cu+Cu2S is physically separated from the rare earth oxide and oxysulfide to avoid complications arising from interaction between them. The corresponding equilibrium oxygen partial pressure is measured with an oxide solid electrolyte cell. Gibbs energy change for the conversion of oxide to the corresponding oxysulfide increases monotonically with atomic number of the rare earth element. Second law enthalpy of formation also shows a similar trend. Based on this empirical trend Gibbs energies of formation of oxysulfides of Pr, Eu, Ho, and Er are estimated as a function of temperature.

Borax Mediated Layer-by-Layer Self-Assembly of Neutral Poly(vinyl alcohol) and Chitosan

We report a multilayer film of poly(vinyl alcohol) (PVA)-borate complex and chitosan by using a layer-by-layer approach. PVA is an uncharged polymer, but hydroxyl functional groups of PVA can be crosslinked by using borax as a cross-linking agent. As a result electrostatic charges and intra- and interchain cross-links are introduced in the PVA chain and provide physically cross-linked networks. The PVA-borate was then deposited on a flat Substrate as well as on colloidal particles with chitosan as an oppositely charged polyelectrolyte. Quartz crystal microbalance. scanning electron microscopy, and atomic force microscopy were used to follow the growth of thin film oil flat substrate. Analogous experiments were performed on melamine formaldehyde colloidal particles (3-3.5 mu m) to quantify the process for the preparation of hollow rnicrocapsules. Removal of the core in 0.1 N HCI results in hollow microcapsules. Characterization of microcapsules by transmission electron microscopy revealed formation of stable microcapsules. Further, self-assembly of PVA-borate/chitosan was loaded with the anticancer drug doxorubicin, and release rates were determined at different pH Values to highlight the drug delivery potential of this system.

Electron- Spin-Resonance studies of phase transitions in double propionates

Detailed ESR investigations of Mn2+ substituting for Ca2+ in Ca2Sr(C2H5COO)6, (DSP) and Ca2Pb(C2H5COO)6, (DLP) and Ca2Ba(C2H5COO)6, (DBP), in single crystals and powders, over the temperature range from 300°C to -180°C have been carried out to study the successive phase transitions in these compounds. Spectra have been analyzed in terms of axial spin Hamiltonians and the temperature dependences of the parameters studied. Across the I-II transition, new physically and chemically inequivalent sites appear indicating the disappearance of the diad axes on which the propionate groups are located, bringing out the connection between the motional states of the propionate groups and the occurrence of ferroelectricity. The II-III transition also causes chemically inequivalent sites to develop, indicating that the transitions may not be isomorphous as believed previously. Similarities and dissimilarities of the ESR spectra of DLP, DSP and DBP are discussed in relation to the phase transitions.

ESR studies of phase transitions in double propionates

Detailed ESR investigations of Mn2+ substituting for Ca2+ in Ca2Sr(C2H5COO)6 (DSP), Ca2Pb(C2H5COO)6 (DLP) and Ca2Ba(C2H5COO)6 (DBP), in single crystals and powders, over the temperature range from 200°C to -180°C have been carried out to study the successive phase transitions in these compounds. (DSP: [Tetragonal] ← 8.5°C → [tetragonal, ferroelectric] [tetragonal] ← -169°C → [monoclinic, ferroelectric]; DLP : [tetragonal] ← 60°C → [tetragonal, ferroelectric] ← -71.5°C → [monoclinic, ferroelectric]; [Cubic] ← -6°C → [orthorhombic] ← -75°C → [?]). Spectra have been analysed in terms of axial spin Hamiltonians and the temperature dependences of the parameters studied. In DSP and DLP across the I ↔ II transition, new physically and chemically inequivalent sites appear indicating the disappearance of the diad axes on which the propionate groups are located, bringing out the connection between the motional states of the propionate groups and the occurence of ferroelectricity. The II ↔ III transition also causes chemically inequivalent sites to develop, indicating that the transitions may not be isomorphous as believed previously. In DBP, the -6°C transition leads to (i) a doubling of both physically and chemically inequivalent sites (ii) a small (150 G at -6°C to 170 G at -8°C), but abrupt change in the magnitude of the zero-field splitting tensor D, and (iii) displacements of the orientations of the D tensors. Results are interpreted in terms of alternate rotations of the oxygen octahedra, showing participation of the carboxyl oxygens in the transition. No drastic changes in the parameters occur across the -75°C transition consistent with its second order nature. Similarities and dissimilarities of the ESR spectra of the three compounds in relation to the phase transitions, are discussed.