A Multidimensional Kinetic Upwind Method for Euler Equations

The study of directional derivative lead to the development of a rotationally invariant kinetic upwind method (KUMARI)3 which avoids dimension by dimension splitting. The method is upwind and rotationally invariant and hence truly multidimensional or multidirectional upwind scheme. The extension of KUMARI to second order is as well presented.

Characterization of Clean and Fouled Rail Track Ballast Subsurface Using Seismic Surface Survey Method: Model and Field Studies

The efficiency of track foundation material gradually decreases due to insufficient lateral confinement, ballast fouling, and loss of shear strength of the subsurface soil under cyclic loading. This paper presents characterization of rail track subsurface to identify ballast fouling and subsurface layers shear wave velocity using seismic survey. Seismic surface wave method of multi-channel analysis of surface wave (MASW) has been carried out in the model track and field track for finding out shear wave velocity of the clean and fouled ballast and track subsurface. The shear wave velocity (SWV) of fouled ballast increases with increase in fouling percentage, and reaches a maximum value and then decreases. This character is similar to typical compaction curve of soil, which is used to define optimum and critical fouling percentage (OFP and CFP). Critical fouling percentage of 15 % is noticed for Coal fouled ballast and 25 % is noticed for clayey sand fouled ballast. Coal fouled ballast reaches the OFP and CFP before clayey sand fouled ballast. Fouling of ballast reduces voids in ballast and there by decreases the drainage. Combined plot of permeability and SWV with percentage of fouling shows that after critical fouling point drainage condition of fouled ballast goes below acceptable limit. Shear wave velocities are measured in the selected location in the Wollongong field track by carrying out similar seismic survey. In-situ samples were collected and degrees of fouling were measured. Field SWV values are more than that of the model track SWV values for the same degree of fouling, which might be due to sleeper's confinement. This article also highlights the ballast gradation widely followed in different countries and presents the comparison of Indian ballast gradation with international gradation standards. Indian ballast contains a coarser particle size when compared to other countries. The upper limit of Indian gradation curve matches with lower limit of ballast gradation curves of America and Australia. The ballast gradation followed by Indian railways is poorly graded and more favorable for the drainage conditions. Indian ballast engineering needs extensive research to improve presents track conditions.

An improved method for calculating activities from distribution equilibria

The method of Gibbs-Duhem integration suggested by Speiser et al. has been modified to derive activities from distribution equilibria. It is shown that, in general, the activities of components in melts with a common anion can be calculated, without using their standard Gibbs energies of formation, from eqUilibrium ratios and the knowledge of activities in the metal phase. Moreover, if systems are so chosen that the concentration of one element in the metal phase lies in the Henry's law region (less than 1 %), information on activities in the metal phase is not required. Conversely, activities of elements in an alloy can be readily calculated from equilibrium distribution ratios alone, if the salt phase in equilibrium contains very small amounts of one element. Application of the method is illustrated using distribution ratios from the literature on AgCI-CuCI, AgBr-CuBr, and CuDo.5 -PbD systems. The results indicate that covalent bonding and van der Waals repulsive interactions in certain types of fused salt melts can significantly affect the thermodynamic properties of mixing.

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.

A simple numerical method for three-dimensional analysis of simple expansion chamber mufflers of rectangular as well as circular cross-section with a stationary medium

Three-dimensional effects are a primary source of discrepancy between the measured values of automotive muffler performance and those predicted by the plane wave theory at higher frequencies. The basically exact method of (truncated) eigenfunction expansions for simple expansion chambers involves very complicated algebra, and the numerical finite element method requires large computation time and core storage. A simple numerical method is presented in this paper. It makes use of compatibility conditions for acoustic pressure and particle velocity at a number of equally spaced points in the planes of the junctions (or area discontinuities) to generate the required number of algebraic equations for evaluation of the relative amplitudes of the various modes (eigenfunctions), the total number of which is proportional to the area ratio. The method is demonstrated for evaluation of the four-pole parameters of rigid-walled, simple expansion chambers of rectangular as well as circular cross-section for the case of a stationary medium. Computed values of transmission loss are compared with those computed by means of the plane wave theory, in order to highlight the onset (cutting-on) of various higher order modes and the effect thereof on transmission loss of the muffler. These are also compared with predictions of the finite element methods (FEM) and the exact methods involving eigenfunction expansions, in order to demonstrate the accuracy of the simple method presented here.

Thermodynamics of SmMnO3 and SmMn2O5 phases determined by the E.M.F. method

Streszczenie angielskie: Using solid oxide galvanic cells of the type: MnO + Sm2O3 + SmMnO3 / O-2/ Ni + NiO and Mn3O4 + SmMnO3 + SmMn2O5 / O-2 / air the equilibrium oxygen pressure for three-phase equilibria described by the following reactions of formation of ternary phases: MnO + 1/2Sm2O3 + 1/4O2 = SmMnO3 1/3Mn3O4 + SmMnO3 + 1/3O2 = SmMn2O5 was determined in the temperature range from 1173 to 1450 K. From the obtained experimental data the corresponding Gibbs free energy change for above reactions of phases formation was derived: ΔG0f,SmMnO3(+/ - 250J) = -131321(+/ - 2000) + 48.02(+/ - 0:35)T / K ΔG0f,SmMn2O5(+/ - 2000 J) = -107085(+/ - 2200) + 69.74(+/ - 1:70)T / K Using obtained results and available literature data, thermodynamic data tables for the two ternary phases have been compiled from 298.15 to 1400 K. Streszczenie polskie: W pracy przedstawiono wyniki badań dotyczące własności termodynamicznych manganinów samaru, wyznaczone metodą pomiaru SEM ogniw ze stałym elektrolitem: MnO + Sm2O3 + SmMnO3 / O-2/ Ni + NiO ogniwo I Mn3O4 + SmMnO3 + SmMn2O5 / O-2 / powietrze ogniwo II oraz określono równowagowe ciśnienie parcjalne tlenu dla reakcji tworzenia SmMnO3 i SmMn2O5 w zakresie temperatur 1173�1450 K: MnO + 1/2Sm2O3 + 1/4O2 = SmMnO3 1/3Mn3O4 + SmMnO3 + 1/3O2 = SmMn2O5 Z tych danych doświadczalnych wyznaczono zależności temperaturowe energii swobodnych tworzenia powyższych manganinów samaru: ΔG0f,SmMnO3(+/ - 250J) = -131321(+/ - 2000) + 48.02(+/ - 0:35)T / K ΔG0f,SmMn2O5(+/ - 2000 J) = -107085(+/ - 2200) + 69.74(+/ - 1:70)T / K W tablicach I i II zamieszczono dane termodynamiczne dla dwóch potrójnych faz otrzymane poprzez kompilacje własnych danych doświadczalnych z danymi literaturowymi.

Structural, optical, and electrical properties of bulk single crystals of InAsxSb(1-x) grown by rotatory Bridgman method

Radially-homogeneous and single-phase InAsxSb(1−x) crystals, up to 5.0 at. % As concentration, have been grown using the rotatory Bridgman method. Single crystallinity has been confirmed by x-ray and electron diffraction studies. Infrared transmission spectra show a continuous decrease in optical energy gap with the increase of arsenic content in InSb. The measured values of mobility and carrier density at room temperature (for x = .05) are 5.6×104 cm2/V s and 2.04×1016 cm−3, respectively.

A Method of Evaluating Electrical Response of Solion Diodes

A method of evaluating the transient electrical response of a solion diode when excited by different current stimuli is given. This method is extended to obtain the transient response of the solion when connected in a circuit. To illustrate the utility of this method a circuit incorporating a solion diode has been analyzed.

An operator-splitting finite element method for the efficient parallel solution of multidimensional population balance systems

A finite element method for solving multidimensional population balance systems is proposed where the balance of fluid velocity, temperature and solute partial density is considered as a two-dimensional system and the balance of particle size distribution as a three-dimensional one. The method is based on a dimensional splitting into physical space and internal property variables. In addition, the operator splitting allows to decouple the equations for temperature, solute partial density and particle size distribution. Further, a nodal point based parallel finite element algorithm for multi-dimensional population balance systems is presented. The method is applied to study a crystallization process assuming, for simplicity, a size independent growth rate and neglecting agglomeration and breakage of particles. Simulations for different wall temperatures are performed to show the effect of cooling on the crystal growth. Although the method is described in detail only for the case of d=2 space and s=1 internal property variables it has the potential to be extendable to d+s variables, d=2, 3 and s >= 1. (C) 2011 Elsevier Ltd. All rights reserved.

A density matrix renormalization group method study of optical properties of porphines and metalloporphines

The symmetrized density matrix renormalization group method is used to study linear and nonlinear optical properties of free base porphine and metalloporphine. Long-range interacting model, namely, Pariser-Parr-Pople model is employed to capture the quantum many-body effect in these systems. The nonlinear optical coefficients are computed within the correction vector method. The computed singlet and triplet low-lying excited state energies and their charge densities are in excellent agreement with experimental as well as many other theoretical results. The rearrangement of the charge density at carbon and nitrogen sites, on excitation, is discussed. From our bond order calculation, we conclude that porphine is well described by the 18-annulenic structure in the ground state and the molecule expands upon excitation. We have modeled the regular metalloporphine by taking an effective electric field due to the metal ion and computed the excitation spectrum. Metalloporphines have D(4h) symmetry and hence have more degenerate excited states. The ground state of metalloporphines shows 20-annulenic structure, as the charge on the metal ion increases. The linear polarizability seems to increase with the charge initially and then saturates. The same trend is observed in third order polarizability coefficients. (C) 2012 American Institute of Physics. [doi: 10.1063/1.3671946]