Instability, intermittency, and multiscaling in discrete growth models of kinetic roughening

We show by numerical simulations that discretized versions of commonly studied continuum nonlinear growth equations (such as the Kardar-Parisi-Zhangequation and the Lai-Das Sarma-Villain equation) and related atomistic models of epitaxial growth have a generic instability in which isolated pillars (or grooves) on an otherwise flat interface grow in time when their height (or depth) exceeds a critical value. Depending on the details of the model, the instability found in the discretized version may or may not be present in the truly continuum growth equation, indicating that the behavior of discretized nonlinear growth equations may be very different from that of their continuum counterparts. This instability can be controlled either by the introduction of higher-order nonlinear terms with appropriate coefficients or by restricting the growth of pillars (or grooves) by other means. A number of such ''controlled instability'' models are studied by simulation. For appropriate choice of the parameters used for controlling the instability, these models exhibit intermittent behavior, characterized by multiexponent scaling of height fluctuations, over the time interval during which the instability is active. The behavior found in this regime is very similar to the ''turbulent'' behavior observed in recent simulations of several one- and two-dimensional atomistic models of epitaxial growth.

Transient analysis of mammalian cell growth in hollow fibre bioreactor

A mathematical model describing the dynamics of mammalian cell growth in hollow fibre bioreactor operated in closed shell mode is developed. Mammalian cells are assumed to grow as an expanding biofilm in the extra-capillary space surrounding the fibre. Diffusion is assumed to be the dominant process in the radial direction while axial convection dominates in the lumen of the bioreactor. The transient simulation results show that steep gradients in the cell number are possible under the condition of substrate limitation. The precise conditions which result in nonuniform growth of cells along the length of the bioreactor are delineated. The effect of various operating conditions, such as substrate feed rate, length of the bioreactor and diffusivity of substrate in different regions of the bioreactor, on the bioreactor performance are evaluated in terms of time required to attain the steady-state. The rime of growth is introduced as a measure of effectiveness factor for the bioreactor and is found to be dependent on two parameters, a modified Peclet number and a Thiele modulus. Diffusion, reaction and/or convection control regimes are identified based on these two parameters. The model is further extended to include dual substrate growth limitations, and the relative growth limiting characteristics of two substrates are evaluated. (C) 1997 Elsevier Science Ltd.

Insulin Growth Factor-2 Binding Protein 3 (IGF2BP3) Is a Glioblastoma-specific Marker That Activates Phosphatidylinositol 3-Kinase/Mitogen-activated Protein Kinase (PI3K/MAPK) Pathways by Modulating IGF-2

Glioblastoma is the most common and malignant form of primary astrocytoma. Upon investigation of the insulin-like growth factor (IGF) pathway, we found the IGF2BP3/IMP3 transcript and protein to be up-regulated in GBMs but not in lower grade astrocytomas (p<0.0001). IMP3 is an RNA binding protein known to bind to the 5'-untranslated region of IGF-2 mRNA, thereby activating its translation. Overexpression-and knockdown-based studies establish a role for IMP3 in promoting proliferation, anchorage-independent growth, invasion, and chemoresistance. IMP3 overexpressing B16F10 cells also showed increased tumor growth, angiogenesis, and metastasis, resulting in poor survival in a mouse model. Additionally, the infiltrating front, perivascular, and subpial regions in a majority of the GBMs stained positive for IMP3. Furthermore, two different murine glioma models were used to substantiate the above findings. In agreement with the translation activation functions of IMP3, we also found increased IGF-2 protein in the GBM tumor samples without a corresponding increase in its transcript levels. Also, in vitro IMP3 overexpression/knockdown modulated the IGF-2 protein levels without altering its transcript levels. Additionally, IGF-2 neutralization and supplementation studies established that the proproliferative effects of IMP3 were indeed mediated through IGF-2. Concordantly, PI3K and MAPK, the downstream effectors of IGF-2, are activated by IMP3 and are found to be essential for IMP3-induced cell proliferation. Thus, we have identified IMP3 as a GBM-specific proproliferative and proinvasive marker acting through IGF-2 resulting in the activation of oncogenic PI3K and MAPK pathways.

New Ventures’ influence on the Growth of SME Employment, Production and Exports in India: Does Institutional Credit matter?

New ventures are considered to be a major source of small firm growth. In Indian context the contribution of new ventures in terms of new employment, production and exports has largely remained unexplored. It is equally important and unexplored, the significance of the contribution of bank credit to the growth of new ventures in India. This paper is an attempt to throw light on these two aspects. The research is based on secondary data of the liberalized period provided by Ministry of Micro, Small and Medium Enterprises, Government of India and Reserve Bank of India. To analyze the influence of bank credit growth on new ventures and the influence of new ventures on growth of additional employment, additional production and additional exports, we used a Bi-Variate Vector Auto Regression. Based on the model generated, Granger causality tests are conducted to obtain the results. The study found that rate of growth of bank credit causes the number of new ventures, implying any increase in the rate of growth of bank credit will be beneficial to the growth of new ventures. The study also concluded that new ventures are not causing the growth of additional employment or additional production. However new ventures cause the growth of additional exports. This is reasonable as entrepreneurs start their new ventures with minimum possible employment and relatively low rate of capacity utilization and they come up to take advantage of the process of globalization by catering to the international market.

Linear stability, transient growth and the role of viscosity stratification in compressible Couette flow

Linear stability and the nonmodal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the uniform shear flow with constant viscosity, and (b) the nonuniform shear flow with stratified viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers (M). For a given M, the critical Reynolds number (Re) is significantly smaller for the uniform shear flow than its nonuniform shear counterpart; for a given Re, the dominant instability (over all streamwise wave numbers, α) of each mean flow belongs to different modes for a range of supersonic M. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean flow to perturbations. It is shown that the energy transfer from mean flow occurs close to the moving top wall for “mode I” instability, whereas it occurs in the bulk of the flow domain for “mode II.” For the nonmodal transient growth analysis, it is shown that the maximum temporal amplification of perturbation energy, Gmax, and the corresponding time scale are significantly larger for the uniform shear case compared to those for its nonuniform counterpart. For α=0, the linear stability operator can be partitioned into L∼L̅ +Re2 Lp, and the Re-dependent operator Lp is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: G(t∕Re)∼Re2. In contrast, the dominance of Lp is responsible for the invalidity of this scaling law in nonuniform shear flow. An inviscid reduced model, based on Ellingsen-Palm-type solution, has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and nonmodal instability, it is shown that the viscosity stratification of the underlying mean flow would lead to a delayed transition in compressible Couette flow.

Modeling of growth and motion of equiaxed dendrites in a convecting melt

A numerical micro-scale model is developed to study the behavior of dendrite growth in presence of melt convection. In this method, an explicit, coupled enthalpy model is used to simulate the growth of an equiaxed dendrite, while a Volume of Fluid (VOF) method is used to track the movement of the dendrite in the convecting melt in a two-dimensional Eulerian framework. Numerical results demonstrate the effectiveness of the enthalpy model in simulating the dendritic growth involving complex shape, and the accuracy of VOF method in conserving mass and preserving the complex dendritic shape during motion. Simulations are performed in presence of uniform melt flow for both fixed and moving dendrites, and the difference in dendrite morphology is shown.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

Model projections for household energy use in India

Energy use in developing countries is heterogeneous across households. Present day global energy models are mostly too aggregate to account for this heterogeneity. Here, a bottom-up model for residential energy use that starts from key dynamic concepts on energy use in developing countries is presented and applied to India. Energy use and fuel choice is determined for five end-use functions (cooking, water heating, space heating, lighting and appliances) and for five different income quintiles in rural and urban areas. The paper specifically explores the consequences of different assumptions for income distribution and rural electrification on residential sector energy use and CO(2) emissions, finding that results are clearly sensitive to variations in these parameters. As a result of population and economic growth, total Indian residential energy use is expected to increase by around 65-75% in 2050 compared to 2005, but residential carbon emissions may increase by up to 9-10 times the 2005 level. While a more equal income distribution and rural electrification enhance the transition to commercial fuels and reduce poverty, there is a trade-off in terms of higher CO(2) emissions via increased electricity use. (C) 2011 Elsevier Ltd. All rights reserved.

A model investigation of vegetation-atmosphere interactions on a millennial timescale

A terrestrial biosphere model with dynamic vegetation capability, Integrated Biosphere Simulator (IBIS2), coupled to the NCAR Community Atmosphere Model (CAM2) is used to investigate the multiple climate-forest equilibrium states of the climate system. A 1000-year control simulation and another 1000-year land cover change simulation that consisted of global deforestation for 100 years followed by re-growth of forests for the subsequent 900 years were performed. After several centuries of interactive climate-vegetation dynamics, the land cover change simulation converged to essentially the same climate state as the control simulation. However, the climate system takes about a millennium to reach the control forest state. In the absence of deep ocean feedbacks in our model, the millennial time scale for converging to the original climate state is dictated by long time scales of the vegetation dynamics in the northern high latitudes. Our idealized modeling study suggests that the equilibrium state reached after complete global deforestation followed by re-growth of forests is unlikely to be distinguishable from the control climate. The real world, however, could have multiple climate-forest states since our modeling study is unlikely to have represented all the essential ecological processes (e. g. altered fire regimes, seed sources and seedling establishment dynamics) for the reestablishment of major biomes.

A coupled VOF-IBM-enthalpy approach for modeling motion and growth of equiaxed dendrites in a solidifying melt

A coupled methodology for simulating the simultaneous growth and motion of equiaxed dendrites in solidifying melts is presented. The model uses the volume-averaging principles and combines the features of the enthalpy method for modeling growth, immersed boundary method for handling the rigid solid-liquid interfaces, and the volume of fluid method for tracking the advection of the dendrite. The algorithm also performs explicit-implicit coupling between the techniques used. A two-dimensional framework with incompressible and Newtonian fluid is considered. Validation with available literature is performed and dendrite growth in the presence of rotational and buoyancy driven flow fields is studied. It is seen that the flow fields significantly alter the position and morphology of the dendrites. (C) 2012 Elsevier Inc. All rights reserved.

Growth Kinetics of Nanoclusters in Solution

We study the growth kinetics of nanoclusters in solution. There are two generic factors that drive growth: (a) reactions that produce the nanomaterial; and (b) diffusion of the nanomaterial due to chemical-potential gradients. We model the growth kinetics of ZnO nanoparticles via coupled dynamical equations for the relevant order parameters, We study this model both analytically and numerically. We find that there is a crossover in thenanocluster growth law: from L(t) similar to t(1/2) in the reaction-controlled regime to L(t) t(1/3) in the diffusion-controlled regime.

Fatigue crack growth due to overloads in plain concrete using scaling laws

Scaling laws are represented in power law form and can be utilized to extract the characteristic properties of a new phenomenon with the help of self-similar solutions. In this work, an attempt has been made to propose a scaling law analytically, for plain concrete when subjected to variable amplitude loading. Due to the application of overload on concrete structures, acceleration in the crack growth process takes place. A closed form expression has been developed to capture the acceleration in crack growth rate in conjunction with the principles of dimensional analysis and self-similarity. The proposed model accounts for parameters such as, the tensile strength, fracture toughness, overload effect and the structural size. Knowing the governed and the governing parameters of the physical problem and by using the concepts of self-similarity, a relationship is obtained between the different parameters involved. The predicted results are compared with experimental crack growth data for variable amplitude loading and are found to capture the overload effect with sufficient accuracy. Through a sensitivity analysis, fracture toughness is found to be the most dominant parameter in accelerating the crack length due to application of overload.

Distinct mechanisms of DNA repair in mycobacteria and their implications in attenuation of the pathogen growth

About a third of the human population is estimated to be infected with Mycobacterium tuberculosis. Emergence of drug resistant strains and the protracted treatment strategies have compelled the scientific community to identify newer drug targets, and to develop newer vaccines. In the host macrophages, the bacterium survives within an environment rich in reactive nitrogen and oxygen species capable of damaging its genome. Therefore, for its successful persistence in the host, the pathogen must need robust DNA repair mechanisms. Analysis of M. tuberculosis genome sequence revealed that it lacks mismatch repair pathway suggesting a greater role for other DNA repair pathways such as the nucleotide excision repair, and base excision repair pathways. In this article, we summarize the outcome of research involving these two repair pathways in mycobacteria focusing primarily on our own efforts. Our findings, using Mycobacterium smegmatis model, suggest that deficiency of various DNA repair functions in single or in combinations severely compromises their DNA repair capacity and attenuates their growth under conditions typically encountered in macrophages. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Vapor-liquid-solid growth of Si nanowires: A kinetic analysis

A steady state kinetic model has been developed for the vapor-liquid-solid growth of Si whiskers or nanowires from liquid catalyst droplets. The steady state is defined as one in which the net injection rate of Si into the droplet is equal to the ejection rate due to wire growth. Expressions that represent specific mechanisms of injection and ejection of Si atoms from the liquid catalyst droplet have been used and their relative importance has been discussed. The analysis shows that evaporation and reverse reaction rates need to be invoked, apart from just surface cracking of the precursor, in order to make the growth rate radius dependent. When these pathways can be neglected, the growth rate become radius independent and can be used to determine the activation energies for the rate limiting step of heterogeneous precursor decomposition. The ejection rates depend on the mechanism of wire growth at the liquid-solid interface or the liquid-solid-vapor triple phase boundary. It is shown that when wire growth is by nucleation and motion of ledges, a radius dependence of growth rate does not just come from the Gibbs-Thompson effect on supersaturation in the liquid, but also from the dependence of the actual area or length available for nucleation. Growth rates have been calculated using the framework of equations developed and compared with experimental results. The agreement in trends is found to be excellent. The same framework of equations has also been used to account for the diverse pressure and temperature dependence of growth rates reported in the literature. © 2012 American Institute of Physics.

Analysis of fatigue crack growth in RC beams

In this work, a fatigue crack propagation model developed using dimensional analysis for plain concrete is used in conjunction with the steel closing force to predict the crack growth behavior of reinforced concrete beams. A numerical procedure is followed using the proposed model to compute the fatigue life of RC beams and the dissipated energy in the steel reinforcement due to shake down behavior. Through a sensitivity study, it is found that the structural size is the most sensitive parameter on which the crack growth rate is dependent. Furthermore, the moment carrying capacity of an RC beam is computed as function of crack size by considering the effect of bond slip.