A Dual Scale Model for Macrosegregation in Alloy Solidification

A numerical approach for coupling the temperature and concentration fields using a micro/macro dual scale model for a solidification problem is presented. The dual scale modeling framework is implemented on a hybrid explicit-implicit solidification scheme. The advantage of this model lies in more accurate consideration of microsegregation occurring at micro-scale using a subgrid model. The model is applied to the case of solidification of a Pb-40% Sn alloy in a rectangular cavity. The present simulation results are compared with the corresponding experimental results reported in the literature, showing improvement in macrosegregation predictions. Subsequently, a comparison of macrosegregation prediction between the results of the present method with those of a parameter model is performed, showing similar trends.

Identification of delamination in composite beams using spectral estimation and a genetic algorithm

An efficient strategy for identification of delamination in composite beams and connected structures is presented. A spectral finite-element model consisting of a damaged spectral element is used for model-based prediction of the damaged structural response in the frequency domain. A genetic algorithm (GA) specially tailored for damage identification is derived and is integrated with finite-element code for automation. For best application of the GA, sensitivities of various objective functions with respect to delamination parameters are studied and important conclusions are presented. Model-based simulations of increasing complexity illustrate some of the attractive features of the strategy in terms of accuracy as well as computational cost. This shows the possibility of using such strategies for the development of smart structural health monitoring softwares and systems.

Phase diagram of a two-species lattice model with a linear instability

We discuss the properties of a one-dimensional lattice model of a driven system with two species of particles in which the mobility of one species depends on the density of the other. This model was introduced by Lahiri and Ramaswamy (Phys. Rev. Lett., 79, 1150 (1997)) in the context of sedimenting colloidal crystals, and its continuum version was shown to exhibit an instability arising from linear gradient couplings. In this paper we review recent progress in understanding the full phase diagram of the model. There are three phases. In the first, the steady state can be determined exactly along a representative locus using the condition of detailed balance. The system shows phase separation of an exceptionally robust sort, termed strong phase separation, which survives at all temperatures. The second phase arises in the threshold case where the first species evolves independently of the second, but the fluctuations of the first influence the evolution of the second, as in the passive scalar problem. The second species then shows phase separation of a delicate sort, in which long-range order coexists with fluctuations which do not damp down in the large-size limit. This fluctuation-dominated phase ordering is associated with power law decays in cluster size distributions and a breakdown of the Porod law. The third phase is one with a uniform overall density, and along a representative locus the steady state is shown to have product measure form. Density fluctuations are transported by two kinematic waves, each involving both species and coupled at the nonlinear level. Their dissipation properties are governed by the symmetries of these couplings, which depend on the overall densities. In the most interesting case,, the dissipation of the two modes is characterized by different critical exponents, despite the nonlinear coupling.

An extended Cahn-Hilliard model for interfaces with cubic anisotropy

For studying systems with a cubic anisotropy in interfacial energy sigma, we extend the Cahn-Hilliard model by including in it a fourth-rank term, namely, gamma (ijlm) [partial derivative (2) c/(partial derivativex(i) partial derivativex(j))] [partial derivative (2) c/(partial derivativex(l) partial derivativex(m))]. This term leads to an additional linear term in the evolution equation for the composition parameter field. It also leads to an orientation-dependent effective fourth-rank coefficient gamma ([hkl]) in the governing equation for the one-dimensional composition profile across a planar interface. The main effect of a non-negative gamma ([hkl]) is to increase both sigma and interfacial width w, each of which, upon suitable scaling, is related to gamma ([hkl]) through a universal scaling function. In this model, sigma is a differentiable function of interface orientation (n) over cap, and does not exhibit cusps; therefore, the equilibrium particle shapes (Wulff shapes) do not contain planar facets. However, the anisotropy in the interfacial energy can be large enough to give rise to corners in the Wulff shapes in two dimensions. In particles of finite sizes, the corners become rounded, and their shapes tend towards the Wulff shape with increasing particle size.

NETGEM: Network Embedded Temporal GEnerative Model for gene expression data

Background: Temporal analysis of gene expression data has been limited to identifying genes whose expression varies with time and/or correlation between genes that have similar temporal profiles. Often, the methods do not consider the underlying network constraints that connect the genes. It is becoming increasingly evident that interactions change substantially with time. Thus far, there is no systematic method to relate the temporal changes in gene expression to the dynamics of interactions between them. Information on interaction dynamics would open up possibilities for discovering new mechanisms of regulation by providing valuable insight into identifying time-sensitive interactions as well as permit studies on the effect of a genetic perturbation. Results: We present NETGEM, a tractable model rooted in Markov dynamics, for analyzing the dynamics of the interactions between proteins based on the dynamics of the expression changes of the genes that encode them. The model treats the interaction strengths as random variables which are modulated by suitable priors. This approach is necessitated by the extremely small sample size of the datasets, relative to the number of interactions. The model is amenable to a linear time algorithm for efficient inference. Using temporal gene expression data, NETGEM was successful in identifying (i) temporal interactions and determining their strength, (ii) functional categories of the actively interacting partners and (iii) dynamics of interactions in perturbed networks. Conclusions: NETGEM represents an optimal trade-off between model complexity and data requirement. It was able to deduce actively interacting genes and functional categories from temporal gene expression data. It permits inference by incorporating the information available in perturbed networks. Given that the inputs to NETGEM are only the network and the temporal variation of the nodes, this algorithm promises to have widespread applications, beyond biological systems. The source code for NETGEM is available from https://github.com/vjethava/NETGEM

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]

Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Real-time density matrix renormalization group dynamics of spin and charge transport in push-pull polyenes and related systems

In this paper we investigate the effect of terminal substituents on the dynamics of spin and charge transport in donor-acceptor substituted polyenes [D-(CH)(x)-A] chains, also known as push-pull polyenes. We employ a long-range correlated model Hamiltonian for the D-(CH)(x)-A system, and time-dependent density matrix renormalization group technique for time propagating the wave packet obtained by injecting a hole at a terminal site, in the ground state of the system. Our studies reveal that the end groups do not affect spin and charge velocities in any significant way, but change the amount of charge transported. We have compared these push-pull systems with donor-acceptor substituted polymethine imine (PMI), D-(CHN)(x)-A, systems in which besides electron affinities, the nature of p(z) orbitals in conjugation also alternate from site to site. We note that spin and charge dynamics in the PMIs are very different from that observed in the case of push-pull polyenes, and within the time scale of our studies, transport of spin and charge leads to the formation of a ``quasi-static'' state.

Genetic Algorithm for Damage Assessment

Genetic Algorithms (GAs) are recognized as an alternative class of computational model, which mimic natural evolution to solve problems in a wide domain including machine learning, music generation, genetic synthesis etc. In the present study Genetic Algorithm has been employed to obtain damage assessment of composite structural elements. It is considered that a state of damage can be modeled as reduction in stiffness. The task is to determine the magnitude and location of damage. In a composite plate that is discretized into a set of finite elements, if a jth element is damaged, the GA based technique will predict the reduction in Ex and Ey and the location j. The fact that the natural frequency decreases with decrease in stiffness is made use of in the method. The natural frequency of any two modes of the damaged plates for the assumed damage parameters is facilitated by the use of Eigen sensitivity analysis. The Eigen value sensitivities are the derivatives of the Eigen values with respect to certain design parameters. If ωiu is the natural frequency of the ith mode of the undamaged plate and ωid is that of the damaged plate, with δωi as the difference between the two, while δωk is a similar difference in the kth mode, R is defined as the ratio of the two. For a random selection of Ex,Ey and j, a ratio Ri is obtained. A proper combination of Ex,Ey and j which makes Ri−R=0 is obtained by Genetic Algorithm.

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.

An EDC-based turbulent premixed combustion model

In the present study, a new turbulent premixed combustion model is proposed by integrating the Coherent Flame Model with the modified eddy dissipation concept, and relating the fine structure mass fraction to the flame surface density. First, experimental results of turbulent flame speed available from literature are compared with the predicted results at different turbulence intensities to validate the flame surface density model. It is observed that the model is able to predict the turbulent burning speeds accurately. Then, a comprehensive validation is carried out utilizing data on a turbulent lifted methane flame issuing into a vitiated co-flow. Detailed comparison of temperature and species concentrations between experiment and simulation is performed at different heights of the flame. Overall, the model is found to predict both the spatial variation and peak values of the scalars at various heights satisfactorily.

Density matrix renormalization group algorithm for Bethe lattices of spin-1/2 or spin-1 sites with Heisenberg antiferromagnetic exchange

A density matrix renormalization group (DMRG) algorithm is presented for the Bethe lattice with connectivity Z = 3 and antiferromagnetic exchange between nearest-neighbor spins s = 1/2 or 1 sites in successive generations g. The algorithm is accurate for s = 1 sites. The ground states are magnetic with spin S(g) = 2(g)s, staggered magnetization that persists for large g > 20, and short-range spin correlation functions that decrease exponentially. A finite energy gap to S > S(g) leads to a magnetization plateau in the extended lattice. Closely similar DMRG results for s = 1/2 and 1 are interpreted in terms of an analytical three-site model.

Bose-Hubbard model in a strong effective magnetic field: Emergence of a chiral Mott insulator ground state

Motivated by experiments on Josephson junction arrays, and cold atoms in an optical lattice in a synthetic magnetic field, we study the ``fully frustrated'' Bose-Hubbard model with half a magnetic flux quantum per plaquette. We obtain the phase diagram of this model on a two-leg ladder at integer filling via the density matrix renormalization group approach, complemented by Monte Carlo simulations on an effective classical XY model. The ground state at intermediate correlations is consistently shown to be a chiral Mott insulator (CMI) with a gap to all excitations and staggered loop currents which spontaneously break time-reversal symmetry. We characterize the CMI state as a vortex supersolid or an indirect exciton condensate, and discuss various experimental implications.

Introduction of SV40ER and hTERT into mammospheres generates breast cancer cells with stem cell properties

Emerging evidence suggests that cancers arise in stem/progenitor cells. Yet, the requirements for transformation of these primitive cells remains poorly understood. In this study, we have exploited the `mammosphere' system that selects for primitive mammary stem/progenitor cells to explore their potential and requirements for transformation. Introduction of Simian Virus 40 Early Region and hTERT into mammosphere-derived cells led to the generation of NBLE, an immortalized mammary epithelial cell line. The NBLEs largely comprised of bi-potent progenitors with long-term self-renewal and multi-lineage differentiation potential. Clonal and karyotype analyses revealed the existence of heterogeneous population within NBLEs with varied proliferation, differentiation and sphere-forming potential. Significantly, injection of NBLEs into immunocompromised mice resulted in the generation of invasive ductal adenocarcinomas. Further, these cells harbored a sub-population of CD44(+)/CD24(-) fraction that alone had sphere- and tumor-initiating potential and resembled the breast cancer stem cell gene signature. Interestingly, prolonged in vitro culturing led to their further enrichment. The NBLE cells also showed increased expression of stemness and epithelial to mesenchymal transition markers, deregulated self-renewal pathways, activated DNA-damage response and cancer-associated chromosomal aberrations-all of which are likely to have contributed to their tumorigenic transformation. Thus, unlike previous in vitro transformation studies that used adherent, more differentiated human mammary epithelial cells our study demonstrates that the mammosphere-derived, less-differentiated cells undergo tumorigenic conversion with only two genetic elements, without requiring oncogenic Ras. Moreover, the striking phenotypic and molecular resemblance of the NBLE-generated tumors with naturally arising breast adenocarcinomas supports the notion of a primitive breast cell as the origin for this subtype of breast cancer. Finally, the NBLEs represent a heterogeneous population of cells with striking plasticity, capable of differentiation, self-renewal and tumorigenicity, thus offering a unique model system to study the molecular mechanisms involved with these processes. Oncogene (2012) 31, 1896-1909; doi:10.1038/onc.2011.378; published online 29 August 2011

Spectral-envelope and group-delay models for transient signals—applications to castanets and stop consonants

We present a novel approach to represent transients using spectral-domain amplitude-modulated/frequency -modulated (AM-FM) functions. The model is applied to the real and imaginary parts of the Fourier transform (FT) of the transient. The suitability of the model lies in the observation that since transients are well-localized in time, the real and imaginary parts of the Fourier spectrum have a modulation structure. The spectral AM is the envelope and the spectral FM is the group delay function. The group delay is estimated using spectral zero-crossings and the spectral envelope is estimated using a coherent demodulator. We show that the proposed technique is robust to additive noise. We present applications of the proposed technique to castanets and stop-consonants in speech.