Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

Dielectric anomaly in Li-doped zinc oxide thin films grown by sol-gel route

Sol-gel route was employed to grow polycrystalline thin films of Li-doped ZnO thin films (Zn1-xLixO, x=0.15). Polycrystalline films were obtained at a growth temperature of 400-500 degrees C. Ferroelectricity in Zn0.85Li0.15O was verified by examining the temperature variation of the real and imaginary parts of dielectric constant, and from the C-V measurements. The phase transition temperature was found to be 330 K. The room-temperature dielectric constant and dissipation factor were 15.5 and 0.09 respectively, at a frequency of 100 kHz. The films exhibited well-defined hysteresis loop, and the values of spontaneous polarization (P-s) and coercive field were 0.15 mu C/cm(2) and 20 kV/cm, respectively, confirming the presence of ferroelectricity.

A dynamic renormalization group study of active nematics

We carry out a systematic construction of the coarse-grained dynamical equation of motion for the orientational order parameter for a two-dimensional active nematic, that is a nonequilibrium steady state with uniaxial, apolar orientational order. Using the dynamical renormalization group, we show that the leading nonlinearities in this equation are marginally irrelevant. We discover a special limit of parameters in which the equation of motion for the angle field bears a close relation to the 2d stochastic Burgers equation. We find nevertheless that, unlike for the Burgers problem, the nonlinearity is marginally irrelevant even in this special limit, as a result of a hidden fluctuation-dissipation relation. 2d active nematics therefore have quasi-long-range order, just like their equilibrium counterparts.

Estimation of thermal breakdown voltage of HVDC cables - A theoretical framework

The insulation in a dc cable is subjected to both thermal and electric stress at the same time. While the electric stress is generic to the cable, the temperature rise in the insulation is, by and large, due to the Ohmic losses in the conductor. The consequence of this synergic effect is to reduce the maximum operating voltage and causes a premature failure of the cable. The authors examine this subject in some detail and propose a comprehensive theoretical formulation relating the maximum thermal voltage (MTV) to the physical and geometrical parameters of the insulation. The heat flow patterns and boundary conditions considered by the authors here and those found in earlier literature are provided. The MTV of a dc cable is shown to be a function of the load current apart from the resistance of the insulation. The results obtained using the expressions, developed by the authors, are compared with relevant results published in the literature and found to be in close conformity.

Quadrature generation techniques for frequency multiplication based oscillators

Frequency multiplication (FM) can be used to design low power frequency synthesizers. This is achieved by running the VCO at a much reduced frequency, while employing a power efficient frequency multiplier, and also thereby eliminating the first few dividers. Quadrature signals can be generated by frequency- multiplying low frequency I/Q signals, however this also multiplies the quadrature error of these signals. Another way is generating additional edges from the low-frequency oscillator (LFO) and develop a quadrature FM. This makes the I-Q precision heavily dependent on process mismatches in the ring oscillator. In this paper we examine the use of fewer edges from LFO and a single stage polyphase filter to generate approximate quadrature signals, which is then followed by an injection-locked quadrature VCO to generate high- precision I/Q signals. Simulation comparisons with the existing approach shows that the proposed method offers very good phase accuracy of 0.5deg with only a modest increase in power dissipation for 2.4 GHz IEEE 802.15.4 standard using UMC 0.13 mum RFCMOS technology.

Rheology of dense sheared granular flows

Rapid granular flows are defined as flows in which the time scales for the particle interactions are small compared to the inverse of the strain rate, so that the particle interactions can be treated as instantaneous collisions. We first show, using Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.6 are rapid granular flows. Since collisions are instantaneous, a kinetic theory approach for the constitutive relations is most appropriate, and we present kinetic theory results for different microscopic models for particle interaction. The significant difference between granular flows and normal fluids is that energy is not conserved in a granular flow. The differences in the hydrodynamic modes caused by the non-conserved nature of energy are discussed. Going beyond the Boltzmann equation, the effect of correlations is studied using the ring kinetic approximation, and it is shown that the divergences in the viscometric coefficients, which are present for elastic fluids, are not present for granular flows because energy is not conserved. The hydrodynamic model is applied to the flow down an inclined plane. Since energy is not a conserved variable, the hydrodynamic fields in the bulk of a granular flow are obtained from the mass and momentum conservation equations alone. Energy becomes a relevant variable only in thin 'boundary layers' at the boundaries of the flow where there is a balance between the rates of conduction and dissipation. We show that such a hydrodynamic model can predict the salient features of a chute flow, including the flow initiation when the angle of inclination is increased above the 'friction angle', the striking lack of observable variation of the volume fraction with height, the observation of a steady flow only for certain restitution coefficients, and the density variations in the boundary layers.

Pressure and thermally induced stages of wear in dry sliding of a steel ball against an aluminium-silicon alloy flat

Wear of etched near-eutectic aluminium silicon alloy slid against a steel ball under ambient is explored. The sliding velocity is kept low (0.01 m/s) and the nominal contact pressure is varied in a 15-40 MPa range. Four stages of wear are identified; ultra mild wear, mild wear, severe wear and post severe oxidative wear. The first transition is controlled by the protrusions of silicon particles, projecting out of the aluminium alloy matrix. Once these protrusions disappear under pressure and sliding, oxidation and bulk energy dissipation mechanisms take over to institute transitions to other stages of wear. The phenomenological characteristics of wear stages are explored using a variety of techniques including nanoindentation, focused ion beam milling, electron microscopy, X-ray photoelectron spectroscopy (XPS), energy dispersive X-ray spectroscopy (EDS) and optical interferometry. (c) 2010 Elsevier B.V. All rights reserved.

Structural studies of decaying fluid turbulence: Effect of initial conditions

We present results from a systematic numerical study of structural properties of an unforced, incompressible, homogeneous, and isotropic three-dimensional turbulent fluid with an initial energy spectrum that develops a cascade of kinetic energy to large wave numbers. The results are compared with those from a recently studied set of power-law initial energy spectra [C. Kalelkar and R. Pandit, Phys. Rev. E 69, 046304 (2004)] which do not exhibit such a cascade. Differences are exhibited in plots of vorticity isosurfaces, the temporal evolution of the kinetic energy-dissipation rate, and the rates of production of the mean enstrophy along the principal axes of the strain-rate tensor. A crossover between "non-cascade-type" and "cascade-type" behavior is shown numerically for a specific set of initial energy spectra.

Fermi level depinning at the germanium Schottky interface through sulfur passivation

We demonstrate the depinning of Fermi level on both p- and n-type germanium after sulfur passivation by aqueous (NH4)(2)S treatment. Schottky contacts realized using metals with a wide range of work functions produce nearly ideal behavior confirming that the Fermi level is depinned. Examination of the passivated surface using x-ray photoelectron spectroscopy reveals bonding between Ge and sulfur.It is shown that good Ohmic contacts to n-type Ge and a hole barrier height (phi(Bp)) of 0.6 eV to p-type Ge can be achieved after this passivation treatment, with Zr Schottky contacts. This is the highest phi(Bp) reported so far.

Power Optimal Network-on-Chip Interconnect Design

A large part of today's multi-core chips is interconnect. Increasing communication complexity has made essential new strategies for interconnects, such as Network on Chip. Power dissipation in interconnects has become a substantial part of the total power dissipation. Techniques to reduce interconnect power have thus become a necessity. In this paper, we present a design methodology that gives values of bus width for interconnect links, frequency of operation for routers, in Network on Chip scenario that satisfy required throughput and dissipate minimal switching power. We develop closed form analytical expressions for the power dissipation, with bus width and frequency as variables and then use Lagrange multiplier method to arrive at the optimal values. We present a 4 port router in 90 nm technology library as case study. The results obtained from analysis are discussed.

Generic Routing Rules and a Scalable Access Enhancement for the Network-on-Chip RECONNECT

RECONNECT is a Network-on-Chip using a honeycomb topology. In this paper we focus on properties of general rules applicable to a variety of routing algorithms for the NoC which take into account the missing links of the honeycomb topology when compared to a mesh. We also extend the original proposal [5] and show a method to insert and extract data to and from the network. Access Routers at the boundary of the execution fabric establish connections to multiple periphery modules and create a torus to decrease the node distances. Our approach is scalable and ensures homogeneity among the compute elements in the NoC. We synthesized and evaluated the proposed enhancement in terms of power dissipation and area. Our results indicate that the impact of necessary alterations to the fabric is negligible and effects the data transfer between the fabric and the periphery only marginally.

Recombination and trapping in n-type cobalt-doped germanium

The role of cobalt centers in promoting the recombination and trapping processes in n-type germanium has been investigated. Data on lifetime measurements carried out by the steadystate photoconductivity and photo-magneto-electric methods in the temperature range 145 to 300°K on n-type germanium samples containing cobalt in the concentration range 1·1013 to 5.·014/cm3 are presented. The results are analysed on the basis of Sah-Shockley's multi-level formula to yield the capture cross-sections Sp= (hole capture cross-section at doubly negatively charged center) and Sn-(electron capture cross-section at singly negatively charged center) and temperature dependence. It is found that Sp= is (22 ± 6). 10-16 cm2 and Sn- is ∼ 0·1. 10-16 cm2 at 145°K. Sp= varies (n = 3·5 to 4·5) in the range 145-220°K; above 225°K the index 'n' tends to a smaller value. Sn- is practically temperature independent below 180°K and increases with increase of temperature above 180°K. The value of Sp= and its temperature variation lead one to the conclusion that during capture at attractive centers, the phonon cascade mechanism is responsible for the dissipation of the recombination energy.

A new method of producing local enhancement of buoyancy in liquid flows

We describe here a novel method of generating large volumetric heating in a liquid. The method uses the principle of ohmic heating of the liquid, rendered electrically conducting by suitable additives if necessary. Electrolysis is prevented by the use of high frequency alternating voltage and chemically treated electrodes. The technique is demonstrated by producing substantial heating in an initially neutral jet of water. Simple flow visualisation studies, made by adding dye to the jet, show marked changes in the growth and development of the jet with heat addition.

Probing top-gated field effect transistor of reduced graphene oxide monolayer made by dielectrophoresis

We demonstrate a top-gated field effect transistor made of a reduced graphene oxide (RGO) monolayer (graphene) by dielectrophoresis. The Raman spectrum of RGO flakes of typical size of 5 mu m x 5 mu m shows a single 2D band at 2687 cm(-1), characteristic of single-layer graphene.The two-probe current-voltage measurements of RGO flakes, deposited in between the patterned electrodes with a gap of 2.5 mu m using ac dielectrophoresis, show ohmic behavior with a resistance of similar to 37 k Omega. The temperature dependence of the resistance (R) of RGO measured between 305 K and 393 K yields a temperature coefficient of resistance [dR/dT]/R similar to -9.5 x 10(-4)/K, the same as that of mechanically exfoliated single-layer graphene. The field-effect transistor action was obtained by electrochemical top-gating using a solid polymer electrolyte (PEO + LiClO4) and Pt wire. The ambipolar nature of graphene flakes is observed up to a doping level of similar to 6 x 10(12)/cm(2) and carrier mobility of similar to 50 cm(2)/V s. The source-drain current characteristics show a tendency of current saturation at high source-drain voltage which is analyzed quantitatively by a diffusive transport model. (C) 2010 Elsevier Ltd. All rights reserved.

Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.