Mathematical model of COREX melter gasifier: Part I. Steady-state model

The COREX melter gasifier is a countercurrent reactor to produce liquid iron. Directly reduced iron (DRI), noncoking coal, and other additives are charged to the melter gasifier at their respective temperatures, and O-2 is blown through the tuyeres. Functionally, a melter gasifier is divided into three zones: a moving bed, fluidized bed, and free board. A model has been developed for the moving bed, where the tuyere region is two-dimensional (2-D) and the rest is one-dimensional (1-D). It is based on multiphase conservation of mass, momentum, and heat. The fluidized bed has been treated as 1-D. Partial equilibrium is calculated for the free board. The calculated temperature of the hot metal, the top gas, and the chemistry of the top gas agree with the reported plant data. The model has been used to study the effects of bed height, injection of impure O-2, coal chemistry, and reactivity on the process performance.

Field emission properties of carbon nanotube arrays with defects and impurities

Field emission from carbon nanotubes (CNTs) in the form of arrays or thin films give rise to several strongly correlated process of electromechanical interaction and degradation. Such processes are mainly due to (1) electron-phonon interaction (2) electromechanical force field leading to stretching of CNTs (3) ballistic transport induced thermal spikes, coupled with high dynamic stress, leading to degradation of emission performance at the device scale. Fairly detailed physics based models of CNTs considering the aspects (1) and (2) above have already been developed by these authors, and numerical results indicate good agreement with experimental results. What is missing in such a system level modeling approach is the incorporation of structural defects and vacancies or charge impurities. This is a practical and important problem due to the fact that degradation of field emission performance is indeed observed in experimental I-V curves. What is not clear from these experiments is whether such degradation in the I-V response is due to dynamic reorientation of the CNTs or due to the defects or due to both of these effects combined. Non-equilibrium Green’s function based simulations using a tight-binding Hamiltonian for single CNT segment show up the localization of carrier density at various locations of the CNTs. About 11% decrease in the drive current with steady difference in the drain current in the range of 0.2-0.4V of the gate voltage was reported in literature when negative charge impurity was introduced at various locations of the CNT over a length of ~20nm. In the context of field emission from CNT tips, a simplistic estimate of defects have been introduced by a correction factor in the Fowler-Nordheim formulae. However, a more detailed physics based treatment is required, while at the same time the device-scale simulation is necessary. The novelty of our present approach is the following. We employ a concept of effective stiffness degradation for segments of CNTs, which is due to structural defects, and subsequently, we incorporate the vacancy defects and charge impurity effects in the Green’s function based approach. Field emission induced current-voltage characteristics of a vertically aligned CNT array on a Cu-Cr substrate is then simulated using a detailed nonlinear mechanistic model of CNTs coupled with quantum hydrodynamics. An array of 10 vertically aligned and each 12 m long CNTs is considered for the device scale analysis. Defect regions are introduced randomly over the CNT length. The result shows the decrease in the longitudinal strain due to defects. Contrary to the expected influence of purely mechanical degradation, this result indicates that the charge impurity and hence weaker transport can lead to a different electromechanical force field, which ultimately can reduce the strain. However, there could be significant fluctuation in such strain field due to electron-phonon coupling. The effect of such fluctuations (with defects) is clearly evident in the field emission current history. The average current also decreases significantly due to such defects.

Model track studies on fouled ballast using ground penetrating radar and multichannel analysis of surface wave

Ballast fouling is created by the breakdown of aggregates or outside contamination by coal dust from coal trains, or from soil intrusion beneath rail track. Due to ballast fouling, the conditions of rail track can be deteriorated considerably depending on the type of fouling material and the degree of fouling. So far there is no comprehensive guideline available to identify the critical degree of fouling for different types of fouling materials. This paper presents the identification of degree of fouling and types of fouling using non-destructive testing, namely seismic surface-wave and ground penetrating radar (GPR) survey. To understand this, a model rail track with different degree of fouling has been constructed in Civil engineering laboratory, University of Wollongong, Australia. Shear wave velocity obtained from seismic survey has been employed to identify the degree of fouling and types of fouling material. It is found that shear wave velocity of fouled ballast increases initially, reaches optimum fouling point (OFP), and decreases when the fouling increases. The degree of fouling corresponding after which the shear wave velocity of fouled ballast will be smaller than that of clean ballast is called the critical fouling point (CFP). Ground penetrating radar with four different ground coupled antennas (500 MHz, 800 MHz, 1.6 GHz and 2.3 GHz) was also used to identify the ballast fouling condition. It is found that the 800 MHz ground coupled antenna gives a better signal in assessing the ballast fouling condition. Seismic survey is relatively slow when compared to GPR survey however it gives quantifiable results. In contrast, GPR survey is faster and better in estimating the depth of fouling. (C) 2011 Elsevier B.V. All rights reserved.

The generalized Onsager model for the secondary flow in a high-speed rotating cylinder

The generalizations of the Onsager model for the radial boundary layer and the Carrier-Maslen model for the end-cap axial boundary layer in a high-speed rotating cylinder are formulated for studying the secondary gas flow due to wall heating and due to insertion of mass, momentum and energy into the cylinder. The generalizations have wider applicability than the original Onsager and Carrier-Maslen models, because they are not restricted to the limit A >> 1, though they are restricted to the limit R e >> 1 and a high-aspect-ratio cylinder whose length/diameter ratio is large. Here, the stratification parameter A = root m Omega(2)R(2)/2k(B)T). This parameter A is the ratio of the peripheral speed, Omega R, to the most probable molecular speed, root 2k(B)T/m, the Reynolds number Re = rho w Omega R(2)/mu, where m is the molecular mass, Omega and R are the rotational speed and radius of the cylinder, k(B) is the Boltzmann constant, T is the gas temperature, rho(w) is the gas density at wall, and mu is the gas viscosity. In the case of wall forcing, analytical solutions are obtained for the sixth-order generalized Onsager equations for the master potential, and for the fourth-order generalized Carrier-Maslen equation for the velocity potential. For the case of mass/momentum/energy insertion into the flow, the separation-of-variables procedure is used, and the appropriate homogeneous boundary conditions are specified so that the linear operators in the axial and radial directions are self-adjoint. The discrete eigenvalues and eigenfunctions of the linear operators (sixth-order and second-order in the radial and axial directions for the Onsager equation, and fourth-order and second-order in the axial and radial directions for the Carrier-Maslen equation) are determined. These solutions are compared with direct simulation Monte Carlo (DSMC) simulations. The comparison reveals that the boundary conditions in the simulations and analysis have to be matched with care. The commonly used `diffuse reflection' boundary conditions at solid walls in DSMC simulations result in a non-zero slip velocity as well as a `temperature slip' (gas temperature at the wall is different from wall temperature). These have to be incorporated in the analysis in order to make quantitative predictions. In the case of mass/momentum/energy sources within the flow, it is necessary to ensure that the homogeneous boundary conditions are accurately satisfied in the simulations. When these precautions are taken, there is excellent agreement between analysis and simulations, to within 10 %, even when the stratification parameter is as low as 0.707, the Reynolds number is as low as 100 and the aspect ratio (length/diameter) of the cylinder is as low as 2, and the secondary flow velocity is as high as 0.2 times the maximum base flow velocity. The predictions of the generalized models are also significantly better than those of the original Onsager and Carrier-Maslen models, which are restricted to thin boundary layers in the limit of high stratification parameter.

Simulation and experimental study on compositional evolution of Li-Co in LiCoO2 thin films during sputter deposition

The compositional evolution in sputter deposited LiCoO2 thin films is influenced by process parameters involved during deposition. The electrochemical performance of these films strongly depends on their microstructure, preferential orientation and stoichiometry. The transport process of sputtered Li and Co atoms from the LiCoO2 target to the substrate, through Ar plasma in a planar magnetron configuration, was investigated based on the Monte Carlo technique. The effect of sputtering gas pressure and the substrate-target distance (dst) on Li/Co ratio, as well as, energy and angular distribution of sputtered atoms on the substrate were examined. Stable Li/Co ratios have been obtained at 5 Pa pressure and dst in the range 5−11 cm. The kinetic energy and incident angular distribution of Li and Co atoms reaching the substrate have been found to be dependent on sputtering pressure. Simulations were extended to predict compositional variations in films prepared at various process conditions. These results were compared with the composition of films determined experimentally using x-ray photoelectron spectroscopy (XPS). Li/Co ratio calculated using XPS was in moderate agreement with that of the simulated value. The measured film thickness followed the same trend as predicted by simulation. These studies are shown to be useful in understanding the complexities in multicomponent sputtering.

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 breakdown of carbon nanotube devices and the predictability of breakdown position

We have investigated electrical transport properties of long (>10 mu m) multiwalled carbon nanotubes (NTs) by dividing individuals into several segments of identical length. Each segment has different resistance because of the random distribution of defect density in an NT and is corroborated by Raman studies. Higher is the resistance, lower is the current required to break the segments indicating that breakdown occurs at the highly resistive segment/site and not necessarily at the middle. This is consistent with the one-dimensional thermal transport model. We have demonstrated the healing of defects by annealing at moderate temperatures or by current annealing. To strengthen our mechanism, we have carried out electrical breakdown of nitrogen doped NTs (NNTs) with diameter variation from one end to the other. It reveals that the electrical breakdown occurs selectively at the narrower diameter region. Overall, we believe that our results will help to predict the breakdown position of both semiconducting and metallic NTs. Copyright 2012 Author(s). This article is distributed under a Creative Commons Attribution 3.0 Unported License. http://dx.doi.org/10.1063/1.4720426]

Nonequilibrium arrhythmic states and transitions in a mathematical model for diffuse fibrosis in human cardiac tissue

We present a comprehensive numerical study of spiral-and scroll-wave dynamics in a state-of-the-art mathematical model for human ventricular tissue with fiber rotation, transmural heterogeneity, myocytes, and fibroblasts. Our mathematical model introduces fibroblasts randomly, to mimic diffuse fibrosis, in the ten Tusscher-Noble-Noble-Panfilov (TNNP) model for human ventricular tissue; the passive fibroblasts in our model do not exhibit an action potential in the absence of coupling with myocytes; and we allow for a coupling between nearby myocytes and fibroblasts. Our study of a single myocyte-fibroblast (MF) composite, with a single myocyte coupled to N-f fibroblasts via a gap-junctional conductance G(gap), reveals five qualitatively different responses for this composite. Our investigations of two-dimensional domains with a random distribution of fibroblasts in a myocyte background reveal that, as the percentage P-f of fibroblasts increases, the conduction velocity of a plane wave decreases until there is conduction failure. If we consider spiral-wave dynamics in such a medium we find, in two dimensions, a variety of nonequilibrium states, temporally periodic, quasiperiodic, chaotic, and quiescent, and an intricate sequence of transitions between them; we also study the analogous sequence of transitions for three-dimensional scroll waves in a three-dimensional version of our mathematical model that includes both fiber rotation and transmural heterogeneity. We thus elucidate random-fibrosis-induced nonequilibrium transitions, which lead to conduction block for spiral waves in two dimensions and scroll waves in three dimensions. We explore possible experimental implications of our mathematical and numerical studies for plane-, spiral-, and scroll-wave dynamics in cardiac tissue with fibrosis.

Application of Artificial Neural Networks and Particle Swarm Optimization for the Management of Groundwater Resources

Ground management problems are typically solved by the simulation-optimization approach where complex numerical models are used to simulate the groundwater flow and/or contamination transport. These numerical models take a lot of time to solve the management problems and hence become computationally expensive. In this study, Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) models were developed and coupled for the management of groundwater of Dore river basin in France. The Analytic Element Method (AEM) based flow model was developed and used to generate the dataset for the training and testing of the ANN model. This developed ANN-PSO model was applied to minimize the pumping cost of the wells, including cost of the pipe line. The discharge and location of the pumping wells were taken as the decision variable and the ANN-PSO model was applied to find out the optimal location of the wells. The results of the ANN-PSO model are found similar to the results obtained by AEM-PSO model. The results show that the ANN model can reduce the computational burden significantly as it is able to analyze different scenarios, and the ANN-PSO model is capable of identifying the optimal location of wells efficiently.

Force and ATP hydrolysis dependent regulation of RecA nucleoprotein filament by single-stranded DNA binding protein

In Escherichia coli, the filament of RecA formed on single-stranded DNA (ssDNA) is essential for recombinational DNA repair. Although ssDNA-binding protein (SSB) plays a complicated role in RecA reactions in vivo, much of our understanding of the mechanism is based on RecA binding directly to ssDNA. Here we investigate the role of SSB in the regulation of RecA polymerization on ssDNA, based on the differential force responses of a single 576-nucleotide-long ssDNA associated with RecA and SSB. We find that SSB outcompetes higher concentrations of RecA, resulting in inhibition of RecA nucleation. In addition, we find that pre-formed RecA filaments de-polymerize at low force in an ATP hydrolysis- and SSB-dependent manner. At higher forces, re-polymerization takes place, which displaces SSB from ssDNA. These findings provide a physical picture of the competition between RecA and SSB under tension on the scale of the entire nucleoprotein SSB array, which have broad biological implications particularly with regard to competitive molecular binding.

Classification of Vehicles using Magnetic Dipole Model

This paper presents an efficient approach to the modeling and classification of vehicles using the magnetic signature of the vehicle. A database was created using the magnetic signature collected over a wide range of vehicles(cars). A vehicle is modeled as an array of magnetic dipoles. The strength of the magnetic dipole and the separation between the magnetic dipoles varies for different vehicles and is dependent on the metallic composition and configuration of the vehicle. Based on the magnetic dipole data model, we present a novel method to extract a feature vector from the magnetic signature. In the classification of vehicles, a linear support vector machine configuration is used to classify the vehicles based on the obtained feature vectors.

Rice LHS1/OsMADS1 Controls Floret Meristem Specification by Coordinated Regulation of Transcription Factors and Hormone Signaling Pathways

SEPALLATA (SEP) MADS box transcription factors mediate floral development in association with other regulators. Mutants in five rice (Oryza sativa) SEP genes suggest both redundant and unique functions in panicle branching and floret development. LEAFY HULL STERILE1/OsMADS1, from a grass-specific subgroup of LOFSEP genes, is required for specifying a single floret on the spikelet meristem and for floret organ development, but its downstream mechanisms are unknown. Here, key pathways and directly modulated targets of OsMADS1 were deduced from expression analysis after its knockdown and induction in developing florets and by studying its chromatin occupancy at downstream genes. The negative regulation of OsMADS34, another LOFSEP gene, and activation of OsMADS55, a SHORT VEGETATIVE PHASE-like floret meristem identity gene, show its role in facilitating the spikelet-to-floret meristem transition. Direct regulation of other transcription factor genes like OsHB4 (a class III homeodomain Leu zipper member), OsBLH1 (a BEL1-like homeodomain member), OsKANADI2, OsKANADI4, and OsETTIN2 show its role in meristem maintenance, determinacy, and lateral organ development. We found that the OsMADS1 targets OsETTIN1 and OsETTIN2 redundantly ensure carpel differentiation. The multiple effects of OsMADS1 in promoting auxin transport, signaling, and auxin-dependent expression and its direct repression of three cytokinin A-type response regulators show its role in balancing meristem growth, lateral organ differentiation, and determinacy. Overall, we show that OsMADS1 integrates transcriptional and signaling pathways to promote rice floret specification and development.

Assessing future rainfall projections using multiple GCMs and a multi-site stochastic downscaling model

Impact of global warming on daily rainfall is examined using atmospheric variables from five General Circulation Models (GCMs) and a stochastic downscaling model. Daily rainfall at eleven raingauges over Malaprabha catchment of India and National Center for Environmental Prediction (NCEP) reanalysis data at grid points over the catchment for a continuous time period 1971-2000 (current climate) are used to calibrate the downscaling model. The downscaled rainfall simulations obtained using GCM atmospheric variables corresponding to the IPCC-SRES (Intergovernmental Panel for Climate Change - Special Report on Emission Scenarios) A2 emission scenario for the same period are used to validate the results. Following this, future downscaled rainfall projections are constructed and examined for two 20 year time slices viz. 2055 (i.e. 2046-2065) and 2090 (i.e. 2081-2100). The model results show reasonable skill in simulating the rainfall over the study region for the current climate. The downscaled rainfall projections indicate no significant changes in the rainfall regime in this catchment in the future. More specifically, 2% decrease by 2055 and 5% decrease by 2090 in monsoon (HAS) rainfall compared to the current climate (1971-2000) under global warming conditions are noticed. Also, pre-monsoon (JFMAM) and post-monsoon (OND) rainfall is projected to increase respectively, by 2% in 2055 and 6% in 2090 and, 2% in 2055 and 12% in 2090, over the region. On annual basis slight decreases of 1% and 2% are noted for 2055 and 2090, respectively.

Heat and mass transfer and chemical transformation in a cerium nitrate droplet

This paper deals with the thermo-physical changes that a droplet undergoes when it is radiatively heated in a levitated environment. The heat and mass transport model has been developed along with chemical kinetics within a cerium nitrate droplet. The chemical transformation of cerium nitrate to ceria during the process is predicted using Kramers' reaction mechanism which justifies the formation of ceria at a very low temperature as observed in experiments. The rate equation modeled by Kramers is modified suitably to be applicable within the framework of a droplet, and predicts experimental results well in both bulk form of cerium nitrate and in aqueous cerium nitrate droplet. The dependence of dissociation reaction rate on droplet size is determined and the transient mass concentration of unreacted cerium nitrate is reported. The model is validated with experiments both for liquid phase vaporization and chemical reaction. Vaporization and chemical conversion are simulated for different ambient conditions. The competitive effects of sensible heating rate and the rate of vaporization with diffusion of cerium nitrate is seen to play a key role in determining the mass fraction of ceria formed within the droplet. Spatially resolved modeling of the droplet enables the understanding of the conversion of chemical species in more detail.

A Novel Anticancer Agent, 8-Methoxypyrimido4 `,5 `: 4,5]thieno(2,3-b) Quinoline-4(3H)-One Induces Neuro 2a Neuroblastoma Cell Death through p53-Dependent, Caspase-Dependent and Independent Apoptotic Pathways

Neuroblastoma is the most common cancer in infants and fourth most common cancer in children. Despite recent advances in cancer treatments, the prognosis of stage-IV neuroblastoma patients continues to be dismal which warrant new pharmacotherapy. A novel tetracyclic condensed quinoline compound, 8-methoxypyrimido 4 `,5 `: 4,5] thieno(2,3-b) quinoline-4(3H)-one (MPTQ) is a structural analogue of an anticancer drug ellipticine and has been reported to posses anticancer property. Study on MPTQ on neuroblastoma cells is very limited and mechanisms related to its cytotoxicity on neuroblastoma cells are completely unknown. Here, we evaluated the anticancer property of MPTQ on mouse neuro 2a and human SH-SY5Y neuroblastoma cells and investigated the mechanisms underlying MPTQ-mediated neuro 2a cell death. MPTQ-mediated neuro 2a and SH-SY5Y cell deaths were found to be dose and time dependent. Moreover, MPTQ induced cell death reached approximately 99.8% and 90% in neuro 2a and SH-SY5Y cells respectively. Nuclear oligonucleosomal DNA fragmentation and Terminal dUTP Nick End Labelling assays indicated MPTQ-mediated neuro 2a cell death involved apoptosis. MPTQ-mediated apoptosis is associated with increased phosphorylation of p53 at Ser15 and Ser20 which correlates with the hyperphosphorylation of Ataxia-Telangiectasia mutated protein (ATM). Immunocytochemical analysis demonstrated the increased level of Bax protein in MPTQ treated neuro 2a cells. MPTQ-mediated apoptosis is also associated with increased activation of caspase-9, -3 and -7 but not caspase-2 and -8. Furthermore, increased level of caspase-3 and cleaved Poly ( ADP Ribose) polymerase were observed in the nucleus of MPTQ treated neuro 2a cells, suggesting the involvement of caspase-dependent intrinsic but not extrinsic apoptotic pathway. Increased nuclear translocation of apoptosis inducing factor suggests additional involvement of caspase-independent apoptosis pathway in MPTQ treated neuro 2a cells. Collectively, MPTQ-induced neuro 2a cell death is mediated by ATM and p53 activation, and Bax-mediated activation of caspase-dependent and caspase-independent mitochondrial apoptosis pathways.