Constitutive modeling of ice in the high strain rate regime

The objective of the present work is to propose a constitutive model for ice by considering the influence of important parameters such as strain rate dependence and pressure sensitivity on the response of the material. In this regard, the constitutive model proposed by Carney et al. (2006) is considered as a starting basis and subsequently modified to incorporate the effect of brittle cracking within a continuum damage mechanics framework. The damage is taken to occur in the form of distributed cracking within the material during impact which is consistent with experimental observations. At the point of failure, the material is assumed to be fluid-like with deviatoric stress almost dropping down to zero. The constitutive model is implemented in a general purpose finite element code using an explicit formulation. Several single element tests under uniaxial tension and compression, as well as biaxial loading are conducted in order to understand the performance of the model. Few large size simulations are also performed to understand the capability of the model to predict brittle damage evolution in un-notched and notched three point bend specimens. The proposed model predicts lower strength under tensile loading as compared to compressive loading which is in tune with experimental observations. Further the model also asserts the strain rate dependency of the strength behavior under both compressive as well as tensile loading, which also corroborates well with experimental results. (C) 2010 Elsevier Ltd. All rights reserved.

Calcium ionophore, A23187 and its amino acid complexes: spectroscopic and molecular modeling studies.

The circular dichroism, fluorescence, Nuclear Magnetic Resonance and BLM conductance studies indicate that A23187 forms a stable complex with amino acids at low ionophore concentrations (<10(-4)M). However, A23187 prefers to be in a dimeric structure with no significant binding to amino acids, at concentrations higher than 10(-4)M. It was also observed that at lower concentrations, at which the amino acids bind to the ionophore, the affinity for calcium ions was several orders of magnitude lower than that at higher ionophore concentrations. We have also conducted molecular modeling studies to examine the structure of the A23187 dimer and its amino acid complexes. The results of these modeling studies strongly support our experimental results and validate the formation of a hydrogen bonded and energetically stable A23187 dimer and its amino acid complexes.

Modes of binding of 2'-AMP to RNase T1. A computer modeling study

The modes of binding of adenosine 2'-monophosphate (2'-AMP) to the enzyme ribonuclease (RNase) T1 were determined by computer modelling studies. The phosphate moiety of 2'-AMP binds at the primary phosphate binding site. However, adenine can occupy two distinct sites--(1) The primary base binding site where the guanine of 2'-GMP binds and (2) The subsite close to the N1 subsite for the base on the 3'-side of guanine in a guanyl dinucleotide. The minimum energy conformers corresponding to the two modes of binding of 2'-AMP to RNase T1 were found to be of nearly the same energy implying that in solution 2'-AMP binds to the enzyme in both modes. The conformation of the inhibitor and the predicted hydrogen bonding scheme for the RNase T1-2'-AMP complex in the second binding mode (S) agrees well with the reported x-ray crystallographic study. The existence of the first mode of binding explains the experimental observations that RNase T1 catalyses the hydrolysis of phosphodiester bonds adjacent to adenosine at high enzyme concentrations. A comparison of the interactions of 2'-AMP and 2'-GMP with RNase T1 reveals that Glu58 and Asn98 at the phosphate binding site and Glu46 at the base binding site preferentially stabilise the enzyme-2'-GMP complex.

The structures of bacteriorhodopsin with different retinal-Schiff base orientations--computer modeling and energy minimization studies

Bacteriorhodopsin has been the subject of intense study in order to understand its photochemical function. The recent atomic model proposed by Henderson and coworkers based on electron cryo-microscopic studies has helped in understanding many of the structural and functional aspects of bacteriorhodopsin. However, the accuracy of the positions of the side chains is not very high since the model is based on low-resolution data. In this study, we have minimized the energy of this structure of bacteriorhodopsin and analyzed various types of interactions such as - intrahelical and interhelical hydrogen bonds and retinal environment. In order to understand the photochemical action, it is necessary to obtain information on the structures adopted at the intermediate states. In this direction, we have generated some intermediate structures taking into account certain experimental data, by computer modeling studies. Various isomers of retinal with 13-cis and/or 15-cis conformations and all possible staggered orientations of Lys-216 side chain were generated. The resultant structures were examined for the distance between Lys-216-schiff base nitrogen and the carboxylate oxygen atoms of Asp-96 - a residue which is known to reprotonate the schiff base at later stages of photocycle. Some of the structures were selected on the basis of suitable retinal orientation and the stability of these structures were tested by energy minimization studies. Further, the minimized structures are analyzed for the hydrogen bond interactions and retinal environment and the results are compared with those of the minimized rest state structure. The importance of functional groups in stabilizing the structure of bacteriorhodopsin and in participating dynamically during the photocycle have been discussed.

Modeling of precipitation in reverse micellar systems

A model of the precipitation process in reverse micelles has been developed to calculate the size of fine particles obtained therein. While the method shares several features of particle nucleation and growth common to precipitation in large systems, complexities arise in describing the processes of nucleation, due to the extremely small size of a micelle and of particle growth caused by fusion among the micelles. Occupancy of micelles by solubilized molecules is governed by Poisson statistics, implying most of them are empty and cannot nucleate of its own. The model therefore specifies the minimum number of solubilized molecules required to form a nucleus which is used to calculate the homogeneous nucleation rate. Simultaneously, interaction between micelles is assumed to occur by Brownian collision and instantaneous fusion. Analysis of time scales of various events shows growth of particles to be very fast compared to other phenomena occurring. This implies that nonempty micelles either are supersaturated or contain a single precipitated particle and allows application of deterministic population balance equations to describe the evolution of the system with time. The model successfully predicts the experimental measurements of Kandori ct al.(3) on the size of precipitated CaCO3 particles, obtained by carbonation of reverse micelles containing aqueous Ca(OH)(2) solution.

Analysis, Design, Modeling, and Implementation of an Active Clamp HF Link Converter

This paper deals with the system oriented analysis, design, modeling, and implementation of active clamp HF link three phase converter. The main advantage of the topology is reduced size, weight, and cost of the isolation transformer. However, violation of basic power conversion rules due to presence of the leakage inductance in the HF transformer causes over voltage stresses across the cycloconverter devices. It makes use of the snubber circuit necessary in such topologies. The conventional RCD snubbers are dissipative in nature and hence inefficient. The efficiency of the system is greatly improved by using regenerative snubber or active clamp circuit. It consists of an active switching device with an anti-parallel diode and one capacitor to absorb the energy stored in the leakage inductance of the isolation transformer and to regenerate the same without affecting circuit performance. The turn on instant and duration of the active device are selected such that it requires simple commutation requirements. The time domain expressions for circuit dynamics, design criteria of the snubber capacitor with two conflicting constrains (over voltage stress across the devices and the resonating current duration), the simulation results based on generalized circuit model and the experimental results based on laboratory prototype are presented.

Modeling re-ignition and chuffing in solid rocket motors

A nonlinear model is developed to numerically simulate dynamic combustion inside a solid rocket motor chamber. Using this model, the phenomena of re-ignition and chuffing are investigated under low-L* conditions. The model consists of two separate submodels (coupled to each other), one for unsteady burning of propellant and the other for unsteady conservation of mass and energy within the chamber. The latter yields instantaneous pressure and temperature within the chamber. The instantaneous burning rate is calculated using a one-dimensional, nonlinear, transient gas-phase model previously developed by the authors. The results presented in this paper show that the model predicts not only the critical L*, but also the various regimes of L*-instabihty. Specifically, the results exhibit (1) amplifying pressure oscillations leading to extinction, and (2) re-ignition after a dormant period following extinction. The re-ignition could be observed only when a radiation heat flux (from the combustion chamber to the propellant surface) was included. Certain high-frequency oscillations, possibly due to intrinsic instability, are observed when the pressure overshoots during re-ignition. At very low values of initial L*, successive cycles of extinction/reignition displaying typical characteristics of chuffing are predicted. Variations of the chuffing frequency and the thickness of propellant burned off during a chuff with L* are found to be qualitatively the same as that reported from experimental observations.

Modeling of nonreacting and reacting turbulent spray jets using a fully stochastic separated flow approach

Numerical modeling of several turbulent nonreacting and reacting spray jets is carried out using a fully stochastic separated flow (FSSF) approach. As is widely used, the carrier-phase is considered in an Eulerian framework, while the dispersed phase is tracked in a Lagrangian framework following the stochastic separated flow (SSF) model. Various interactions between the two phases are taken into account by means of two-way coupling. Spray evaporation is described using a thermal model with an infinite conductivity in the liquid phase. The gas-phase turbulence terms are closed using the k-epsilon model. A novel mixture fraction based approach is used to stochastically model the fluctuating temperature and composition in the gas phase and these are then used to refine the estimates of the heat and mass transfer rates between the droplets and the surrounding gas-phase. In classical SSF (CSSF) methods, stochastic fluctuations of only the gas-phase velocity are modeled. Successful implementation of the FSSF approach to turbulent nonreacting and reacting spray jets is demonstrated. Results are compared against experimental measurements as well as with predictions using the CSSF approach for both nonreacting and reacting spray jets. The FSSF approach shows little difference from the CSSF predictions for nonreacting spray jets but differences are significant for reacting spray jets. In general, the FSSF approach gives good predictions of the flame length and structure but further improvements in modeling may be needed to improve the accuracy of some details of the Predictions. (C) 2011 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Numerical modeling of the contaminant Transport through stratified soil

The transport processes of the dissolved chemicals in stratified or layered soils have been studied for several decades. In case of the solute transport through stratified layers, interface condition plays an important role in determining appropriate transport parameters. First‐ type and third‐ type interface conditions are generally used in the literature. A first‐type interface condition will result in a continuous concentration profile across the interface at the expense of solute mass balance. On the other hand, a discontinuity in concentration develops when a third‐ type interface condition is used. To overcome this problem, a combined first‐ and third‐ type condition at the interface has been widely employed which yields second‐ type condition. This results in a similar break‐through curve irrespective of the layering order, which is non‐physical. In this work, an interface condition is proposed which satisfies the mass balance implicitly and brings the distinction between the breakthrough curves for different layering sequence corroborating with the experimental observations. This is in disagreement with the earlier work by H. M. Selim and co‐workers but, well agreement with the hypothetical result by Bosma and van der Zee; and Van der Zee.

The Mathematical modeling of inhomogeneites in ventricular tissue

Cardiac arrhythmias such as ventricular tachycardia (VT) or ventricular fibrillation (VF) are the leading cause of death in the industrialised world. There is a growing consensus that these arrhythmias arise because of the formation of spiral waves of electrical activation in cardiac tissue; unbroken spiral waves are associated with VT and broken ones with VF. Several experimental studies have been carried out to determine the effects of inhomogeneities in cardiac tissue on such arrhythmias. We give a brief overview of such experiments, and then an introduction to partial-differential-equation models for ventricular tissue. We show how different types of inhomogeneities can be included in such models, and then discuss various numerical studies, including our own, of the effects of these inhomogeneities on spiral-wave dynamics. The most remarkable qualitative conclusion of our studies is that the spiral-wave dynamics in such systems depends very sensitively on the positions of these inhomogeneities.

Multi-Mode Phonon-Controlled Field Emission from Carbon Nanotubes: Modeling and Experiments

The main idea proposed in this paper is that in a vertically aligned array of short carbon nanotubes (CNTs) grown on a metal substrate, we consider a frequency dependent electric field, so that the mode-specific propagation of phonons, in correspondence with the strained band structure and the dispersion curves, take place. We perform theoretical calculations to validate this idea with a view of optimizing the field emission behavior of the CNT array. This is the first approach of its kind, and is in contrast to the the conventional approach where a DC bias voltage is applied in order to observe field emission. A first set of experimental results presented in this paper gives a clear indication that phonon-assisted control of field emission current in CNT based thin film diode is possible.

Measurement and modeling of Alloy-Spinel-Corundum equilibrium in the Ni-Mn-Al-0 system at 1873 K

The oxygen content of liquid Ni-Mn alloy equilibrated with spinel solid solution, (Ni,Mn)O. (1 +x)A12O3, and α-Al2O3 has been measured by suction sampling and inert gas fusion analysis. The corresponding oxygen potential of the three-phase system has been determined with a solid state cell incorporating (Y2O3)ThO2 as the solid electrolyte and Cr + Cr2O3 as the reference electrode. The equilibrium composition of the spinel phase formed at the interface of the alloy and alumina crucible was obtained using EPMA. The experimental data are compared with a thermodynamic model based on the free energies of formation of end-member spinels, free energy of solution of oxygen in liquid nickel, interaction parameters, and the activities in liquid Ni-Mn alloy and spinel solid solution. Mixing properties of the spinel solid solution are derived from a cation distribution model. The computational results agree with the experimental data on oxygen concentration, potential, and composition of the spinel phase.

Modeling of seepage flow through layered solis

A two-dimensional finite difference model, which solves mixed type of Richards' equation, whose non-linearity is dealt with modified Picard's iteration and strongly implicit procedure to solve the resulting equations, is presented. Modeling of seepage flow through heterogeneous soils, which is common in the field is addressed in the present study. The present model can be applied to both unsaturated and saturated soils and can handle very dry initial condition and steep wetting fronts. The model is validated by comparing experimental results reported in the literature. Newness of this two dimensional model is its application on layered soils with transient seepage face development, which has not been reported in the literature. Application of the two dimensional model for studying unconfined drainage due to sudden drop of water table at seepage face in layered soils is demonstrated. In the present work different sizes of rectangular flow domain with different types of layering are chosen. Sensitivity of seepage height due to problem dimension of layered system is studied. The effect of aspect ratio on seepage face development in case of the flow through layered soil media is demonstrated. The model is also applied to random heterogeneous soils in which the randomness of the model parameters is generated using the turning band technique. The results are discussed in terms of phreatic surface and seepage height development and also flux across the seepage face. Such accurate modeling of seepage face development and quantification of flux moving across the seepage face becomes important while modeling transport problems in variably saturated media.

Characterization and finite element modeling of montmorillonite/polypropylene nanocomposites

Finite element modeling can be a useful tool for predicting the behavior of composite materials and arriving at desirable filler contents for maximizing mechanical performance. In the present study, to corroborate finite element analysis results, quantitative information on the effect of reinforcing polypropylene (PP) with various proportions of nanoclay (in the range of 3-9% by weight) is obtained through experiments; in particular, attention is paid to the Young's modulus, tensile strength and failure strain. Micromechanical finite element analysis combined with Monte Carlo simulation have been carried out to establish the validity of the modeling procedure and accuracy of prediction by comparing against experimentally determined stiffness moduli of nanocomposites. In the same context, predictions of Young's modulus yielded by theoretical micromechanics-based models are compared with experimental results. Macromechanical modeling was done to capture the non-linear stress-strain behavior including failure observed in experiments as this is deemed to be a more viable tool for analyzing products made of nanocomposites including applications of dynamics. (C) 2011 Elsevier Ltd. All rights reserved.

Spectroscopic Characterization and Modeling of Quadrupolar Charge-Transfer Dyes with Bulky Substituents

Joint experimental and theoretical work is presented on two quadrupolar D-pi-A-pi-D chromophores characterized by the same bulky donor (D) group and two different central cores. The first chromophore, a newly synthesized species with a malononitrile-based acceptor (A) group, has a V-shaped structure that makes its absorption spectrum very broad, covering most of the visible region. The second chromophore has a squaraine-based core and therefore a linear structure, as also evinced from its absorption spectra. Both chromophores show an anomalous red shift of the absorption band upon increasing solvent polarity, a feature that is ascribed to the large, bulky structure of the moleCules. For these molecules, the basic description of polar solvation in terms of a uniform reaction field fails. Indeed, a simple extension of the model to account for two independent reaction fields associated with the two molecular arms quantitatively reproduces the observed linear absorption and fluorescence as well as fluorescence anisotropy spectra, fully rationalizing their nontrivial dependence on solvent polarity. The model derived from the analysis of linear spectra is adopted to predict nonlinear spectra and specifically hyper-Rayleigh scattering and two-photon absorption spectra. In polar solvents, the V-shaped chromophore is predicted to have a large HRS response in a wide spectral region (approximately 600-1300 nm). Anomalously large and largely solvent-dependent HRS responses for the linear chromophores are ascribed to symmetry lowering induced by polar solvation and amplified in this bulky system by the presence of two reaction fields.