Modeling of the Effects of Process Variations on Circuit Delay at 65nm

A novel methodology for modeling the effects of process variations on circuit delay performance is proposed by relating the variations in process parameters to variations in delay metric of a complex digital circuit. The delay of a 2-input NAND gate with 65nm gate length transistors is extensively characterized by mixed-mode simulations which is then used as a library element. The variation in saturation current Ionat the device level, and the variation in rising/falling edge stage delay for the NAND gate at the circuit level, are taken as performance metrics. A 4-bit x 4-bit Wallace tree multiplier circuit is used as a representative combinational circuit to demonstrate the proposed methodology. The variation in the multiplier delay is characterized, to obtain delay distributions, by an extensive Monte Carlo analysis. An analytical model based on CV/I metric is proposed, to extend this methodology for a generic technology library with a variety of library elements.

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.

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.

Modeling of room temperature current-voltage measurements on homo-junction HgCdTe diodes exhibiting nonequilibrium effects

HgCdTe mid wave infrared (MWIR) n(+)/nu/p(+) homo-junction photodiodes with planar architecture are designed, fabricated, and measured at room temperature. An improved analytical I-V model is reported by incorporating trap assisted tunneling and electric field enhanced Shockley-Read-Hall generation recombination process due to dislocations. Tunneling currents are fitted before and after the Auger suppression of carriers with energy level of trap (E-t), trap density (N-t), and the doping concentrations of n(+) and nu regions as fitting parameters. Values of E-t and N-t are determined as 0.79 E-g and similar to 9 x 10(14) cm(-3), respectively, in all cases. Doping concentration of nu region was found to exhibit nonequilibrium depletion from a value of 2 x 10(16) to 4 x 10(15) cm(-3) for n(+) doping of 2 x 10(17) cm(-3). Pronounced negative differential resistance is observed in the homo-junction HgCdTe diodes. (C) 2012 American Institute of Physics. [doi:10.1063/1.3682483]

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.

Piecewise Linearization Technique for Compact Charge Modeling of Independent DG MOSFET

Charge linearization techniques have been used over the years in advanced compact models for bulk and double-gate MOSFETs in order to approximate the position along the channel as a quadratic function of the surface potential (or inversion charge densities) so that the terminal charges can be expressed as a compact closed-form function of source and drain end surface potentials (or inversion charge densities). In this paper, in case of the independent double-gate MOSFETs, we show that the same technique could be used to model the terminal charges quite accurately only when the 1-D Poisson solution along the channel is fully hyperbolic in nature or the effective gate voltages are same. However, for other bias conditions, it leads to significant error in terminal charge computation. We further demonstrate that the amount of nonlinearity that prevails between the surface potentials along the channel actually dictates if the conventional charge linearization technique could be applied for a particular bias condition or not. Taking into account this nonlinearity, we propose a compact charge model, which is based on a novel piecewise linearization technique and shows excellent agreement with numerical and Technology Computer-Aided Design (TCAD) simulations for all bias conditions and also preserves the source/drain symmetry which is essential for Radio Frequency (RF) circuit design. The model is implemented in a professional circuit simulator through Verilog-A, and simulation examples for different circuits verify good model convergence.

Modeling of Ascending/Descending Velocity of Metal Droplet Emulsified in Pb-Salt System

Metal-slag emulsion is an important process to enhance the reaction rate between the two phases; thus, it improves the heat and mass transfer of the process significantly. Various experimental studies have been carried out, and some system specific relations have been proposed by various investigators. A unified, theoretical study is lacking to model this complex phenomenon. Therefore, two simple models based on fundamental laws for metal droplet velocity (both ascending and descending) and bubble velocity, as well as its position at any instant of time, have been proposed. Analytical solutions have been obtained for the developed equations. Analytical solutions have been verified for the droplet velocity, traveling time, and size distribution in slag phase by performing high-temperature experiments in a Pb-salt system and comparing the obtained data with theory. The proposed model has also been verified with published experimental data for various liquid systems with a wide range of physical properties. A good agreement has been found between the analytical solution and the experimental and published data in all cases.

Optimal Control of Markov Processes with Age-Dependent Transition Rates

We study optimal control of Markov processes with age-dependent transition rates. The control policy is chosen continuously over time based on the state of the process and its age. We study infinite horizon discounted cost and infinite horizon average cost problems. Our approach is via the construction of an equivalent semi-Markov decision process. We characterise the value function and optimal controls for both discounted and average cost cases.

Markov chain modeling of evolution of strains in reinforced concrete flexural beams

From the analysis of experimentally observed variations in surface strains with loading in reinforced concrete beams, it is noted that there is a need to consider the evolution of strains (with loading) as a stochastic process. Use of Markov Chains for modeling stochastic evolution of strains with loading in reinforced concrete flexural beams is studied in this paper. A simple, yet practically useful, bi-level homogeneous Gaussian Markov Chain (BLHGMC) model is proposed for determining the state of strain in reinforced concrete beams. The BLHGMC model will be useful for predicting behavior/response of reinforced concrete beams leading to more rational design.

Modeling of agglomeration inside a droplet with nanosuspensions in an acoustic field

This short communication reports results of particle agglomeration details of an acoustically levitated nanosilica droplet. The droplet undergoes thermo-physical and morphological changes under external heating load (convective or radiative) forming different solid structures due to particle agglomeration. We report an agglomeration model based on population balance approach coupled with species and energy conservation equations in the liquid phase and compare it with the experimentally observed structure formations using high speed photography. The analysis is able to predict similar spherical bowl shaped morphologies as observed experimentally using scanning electron microscopy and laser induced fluorescence. (C) 2012 Elsevier Ltd. All rights reserved.

Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

Modeling of terahertz heating effects in realistic tissues

Theterahertz (THz) propagation in real tissues causes heating as with any other electromagnetic radiation propagation. A finite-element (FE) model that provides numerical solutions to the heat conduction equation coupled with realistic models of tissues is employed in this study to compute the temperature raise due to THz propagation. The results indicate that the temperature raise is dependent on the tissue type and is highly localized. The developed FE model was validated through obtaining solutions for the steady-state case and showing that they were in good agreement with the analytical solutions. These types of models can also enable computation of specific absorption rates, which are very critical in planning/setting up experiments involving biological tissues.

Modeling of temperature and field-dependent electron mobility in a single-layer graphene sheet

In this paper, we address a physics-based analytical model of electric-field-dependent electron mobility (mu) in a single-layer graphene sheet using the formulation of Landauer and Mc Kelvey's carrier flux approach under finite temperature and quasi-ballistic regime. The energy-dependent, near-elastic scattering rate of in-plane and out-of-plane (flexural) phonons with the electrons are considered to estimate mu over a wide range of temperature. We also demonstrate the variation of mu with carrier concentration as well as the longitudinal electric field. We find that at high electric field (>10(6) Vm(-1)), the mobility falls sharply, exhibiting the scattering between the electrons and flexural phonons. We also note here that under quasi-ballistic transport, the mobility tends to a constant value at low temperature, rather than in between T-2 and T-1 in strongly diffusive regime. Our analytical results agree well with the available experimental data, while the methodologies are put forward to estimate the other carrier-transmission-dependent transport properties.

Rack Level Modeling of Air Flow Through Perforated Tile in a Data Center

Effective air flow distribution through perforated tiles is required to efficiently cool servers in a raised floor data center. We present detailed computational fluid dynamics (CFD) modeling of air flow through a perforated tile and its entrance to the adjacent server rack. The realistic geometrical details of the perforated tile, as well as of the rack are included in the model. Generally, models for air flow through perforated tiles specify a step pressure loss across the tile surface, or porous jump model based on the tile porosity. An improvement to this includes a momentum source specification above the tile to simulate the acceleration of the air flow through the pores, or body force model. In both of these models, geometrical details of tile such as pore locations and shapes are not included. More details increase the grid size as well as the computational time. However, the grid refinement can be controlled to achieve balance between the accuracy and computational time. We compared the results from CFD using geometrical resolution with the porous jump and body force model solution as well as with the measured flow field using particle image velocimetry (PIV) experiments. We observe that including tile geometrical details gives better results as compared to elimination of tile geometrical details and specifying physical models across and above the tile surface. A modification to the body force model is also suggested and improved results were achieved.