The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

The influence of temperature gradient and lowering speed on the melt-solid interface shape of GaxIn1-xSb alloy crystals grown by vertical Bridgman technique

Radially homogeneous bulk alloys of GaxIn1-xSb in the range 0.7 < x < 0.8, have been grown by vertical Bridgman technique. The factors affecting the interface shape during the growth were optimised to achieve zero convexity. From a series of experiments, a critical ratio of the temperature gradient (G) of the furnace at the melting point of the melt composition to the ampoule lowering speed (v) was deduced for attaining the planarity of the melt-solid interface. The studies carried out on directional solidification of Ga0.77In0.23Sb mixed crystals employing planar melt-solid interface exhibited superior quality than those with nonplanar interfaces. The solutions to certain problems encountered during the synthesis and growth of the compound were discussed. (C) 1999 Elsevier Science B.V. All rights reserved.

Thermodynamic modeling of quaternary Al-Ga-In-As system

The present research describes the modeling of the thermodynamic properties of the liquid Al-Ga-In-As alloys at 1073 and 1173 K, and investigates the solid-liquid equilibria in the systems. The isothermal molar excess free energy function for the liquid alloys is represented in terms of 37 parameters pertaining to six of the constituent binaries, four ternaries and the quaternary interactions in the system. The corresponding solid alloys which consist of AlAs, GaAs and InAs are assumed to be quasi-regular ternary solutions. The solidus and liquidus compositions are calculated at 1073 and 1173 K using the derived values of the partial components for the solid and liquid alloys at equilibrium. They are in good agreement with those of the experimentally determined values available in the literature. (C) 1999 Elsevier Science S.A. All rights reserved.

Thermodynamics of surfaces and adsorption in the Fe-S, Fe-N, and Fe-S-N systems at 1823 K

The present investigation calculates the surface tension and adsorption functions of the Fe-S, Fe-N, and Fe-S-N melts at 1823 K using the modified form of Butler's equations and the derived values of the surface interaction parameters of the systems. The calculated values are found to be in good agreement with those of the experimental data of the systems. The present analysis indicates similar adsorption behavior of sulfur for the Fe-S and Fe-S-N systems at 1823 K. Although a linear adsorption behavior is observed in the Fe-N system, an inverse relationship in the adsorption behavior exists between sulur and nitrogen in the Fe-S-N system.

Transplanting the n-terminus from Kv1.4 to Kv1.1 generates an inwardly rectifying K+ channel

A chimeric channel, 4N/1, was generated from two outwardly rectifying K+ channels by linking the N-terminal cytoplasmic domain of hKv1.4 (N terminus ball and chain of hKv1.4) with the transmembrane body of hKvl.l (Delta 78N1 construct of hKvl.l). The recombinant channel has properties similar to the six transmembrane inward rectifiers and opens on hyperpolarization with a threshold of activation at -90 mV. Outward currents are seen on depolarization provided the channel is first exposed to a hyperpolarizing pulse of -100mV or more. Hyperpolarization at and beyond -130mV provides evidence of channel deactivation. Delta 78N1 does not show inward currents on hyperpolarization but does open on depolarizing from -80mV with characteristics similar to native hKvl.l. The outward currents seen in both Delta 78N1 and 4N/1 inactivate slowly at rates consistent with C-type inactivation. The inward rectification of the 4N/1 chimera is consistent with the inactivation gating mechanism. This implies that the addition of the N-terminus from hKv1.4 to hKvl.l shifts channel activation to hyperpolarizing potentials. These results suggest a mechanism involving the N-terminal cytoplasmic domain for conversion of outward rectifiers to inward rectifiers. (C) 1999 Lippincott Williams & Wilkins.

Microwave heating in multiphase systems: evaluation of series solutions

The 1D electric field and heat-conduction equations are solved for a slab where the dielectric properties vary spatially in the sample. Series solutions to the electric field are obtained for systems where the spatial variation in the dielectric properties can be expressed as polynomials. The series solution is used to obtain electric-field distributions for a binary oil-water system where the dielectric properties are assumed to vary linearly within the sample. Using the finite-element method temperature distributions are computed in a three-phase oil, water and rock system where the dielectric properties vary due to the changing oil saturation in the rock. Temperature distributions predicted using a linear variation in the dielectric properties are compared with those obtained using the exact nonlinear variation.

Seismic passive earth pressure coefficients using the method of characteristics

The method of characteristics was used to generate passive earth pressure coefficients for an inclined wall retaining cohesionless backfill material in the presence of pseudostatic horizontal earthquake body forces. The variation of the passive earth pressure coefficients K-pq and K-pgamma with changes in horizontal earthquake acceleration coefficient due to the components of soil unit weight and surcharge pressure, respectively, has been obtained; a closed-form solution for K-pq is also provided. The passive earth resistance has been found to decrease sharply with an increase in the magnitude of horizontal earthquake acceleration. The computed passive earth pressure coefficients were found to be the lowest when compared to all of the previous solutions available in the literature.

The restriction mapping problem revisited

In computational molecular biology, the aim of restriction mapping is to locate the restriction sites of a given enzyme on a DNA molecule. Double digest and partial digest are two well-studied techniques for restriction mapping. While double digest is NP-complete, there is no known polynomial-time algorithm for partial digest. Another disadvantage of the above techniques is that there can be multiple solutions for reconstruction. In this paper, we study a simple technique called labeled partial digest for restriction mapping. We give a fast polynomial time (O(n(2) log n) worst-case) algorithm for finding all the n sites of a DNA molecule using this technique. An important advantage of the algorithm is the unique reconstruction of the DNA molecule from the digest. The technique is also robust in handling errors in fragment lengths which arises in the laboratory. We give a robust O(n(4)) worst-case algorithm that can provably tolerate an absolute error of O(Delta/n) (where Delta is the minimum inter-site distance), while giving a unique reconstruction. We test our theoretical results by simulating the performance of the algorithm on a real DNA molecule. Motivated by the similarity to the labeled partial digest problem, we address a related problem of interest-the de novo peptide sequencing problem (ACM-SIAM Symposium on Discrete Algorithms (SODA), 2000, pp. 389-398), which arises in the reconstruction of the peptide sequence of a protein molecule. We give a simple and efficient algorithm for the problem without using dynamic programming. The algorithm runs in time O(k log k), where k is the number of ions and is an improvement over the algorithm in Chen et al. (C) 2002 Elsevier Science (USA). All rights reserved.

Electrochemical Properties of Li3N Dissolved in Molten LiCl at 900 K

Thermodynamic properties of Li3N dissolved in the molten LiCl salt at 900 K were explored using electrochemical methods. It was difficult to determine precisely the decomposition voltage of Li3N dissolved in the molten salt by cyclic voltammetry. The oxidation wave of N3– ion could not be located with high accuracy. However, the lithium activity of the Pb-Li alloy in equilibrium with the molten salt containing dissolved Li3N under nitrogen atmosphere could be measured electrochemically with high accuracy using the Li/Li + reference electrode. Under the conditions used in this study, the potential of the Li-Pb electrode is equal to the decomposition voltage of Li3N. The activity of Li3N in molten LiCl was determined for anionic fractions of N3– ranging from xN3– = 10–4 to 0.028. The nitride ion concentration in the salt was determined by chemical titration. The activity coefficient of the Li3N at high dilution was found to be very low, around 10–4. The activity coefficient increases sharply with composition and has a value of 0.25 at xN3– = 0.028. ©2001 The Electrochemical Society. All rights reserved.

Novel phase-transition behavior near liquid/liquid critical points of aqueous solutions: Formation of a third phase at the interface

We report a novel phase behavior in aqueous solutions of simple organic solutes near their liquid/liquid critical points, where a solid-like third phase appears at the liquid/liquid interface. The phenomenon has been found in three different laboratories. It appears in many aqueous systems of organic solutes and becomes enhanced upon the addition of salt to these solutions.

Effect of lattice orientation on crack tip constraint in ductile single crystals

In this work, the effect of lattice orientation on the fields prevailing near a notch tip is investigated pertaining to various constraint levels in FCC single crystals. A modified boundary layer formulation is employed and numerical solutions under mode I, plane strain conditions are generated by assuming an elastic-perfectly plastic FCC single crystal. The analysis is carried out corresponding to different lattice orientations with respect to the notch line. It is found that the near-tip deformation field, especially the development of kink or slip shear bands is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the notch tip are also strongly influenced by the level of T-stress. The present results clearly establish that ductile single crystal fracture geometries would progressively lose crack tip constraint as the T-stress becomes more negative irrespective of lattice orientation. Also, the near-tip field for a range of constraint levels can be characterized by two-parameters such as K-T or J-Q as in isotropic plastic solids.

Roy equation analysis of ππ scattering

We analyse the Roy equations for the lowest partial waves of elastic ππ scattering. In the first part of the paper, we review the mathematical properties of these equations as well as their phenomenological applications. In particular, the experimental situation concerning the contributions from intermediate energies and the evaluation of the driving terms are discussed in detail. We then demonstrate that the two S-wave scattering lengths a00 and a02 are the essential parameters in the low energy region: Once these are known, the available experimental information determines the behaviour near threshold to within remarkably small uncertainties. An explicit numerical representation for the energy dependence of the S- and P-waves is given and it is shown that the threshold parameters of the D- and F-waves are also fixed very sharply in terms of a00 and a20. In agreement with earlier work, which is reviewed in some detail, we find that the Roy equations admit physically acceptable solutions only within a band of the (a00,a02) plane. We show that the data on the reactions e+e−→ππ and τ→ππν reduce the width of this band quite significantly. Furthermore, we discuss the relevance of the decay K→ππeν in restricting the allowed range of a00, preparing the grounds for an analysis of the forthcoming precision data on this decay and on pionic atoms. We expect these to reduce the uncertainties in the two basic low energy parameters very substantially, so that a meaningful test of the chiral perturbation theory predictions will become possible.

Analytical solutions for heat transfer during cyclic melting and freezing of a phase change material used in electronic or electrical packaging

In this paper we develop an analytical heat transfer model, which is capable of analyzing cyclic melting and solidification processes of a phase change material used in the context of electronics cooling systems. The model is essentially based on conduction heat transfer, with treatments for convection and radiation embedded inside. The whole solution domain is first divided into two main sub-domains, namely, the melting sub-domain and the solidification sub-domain. Each sub-domain is then analyzed for a number of temporal regimes. Accordingly, analytical solutions for temperature distribution within each subdomain are formulated either using a semi-infinity consideration, or employing a method of quasi-steady state, depending on the applicability. The solution modules are subsequently united, leading to a closed-form solution for the entire problem. The analytical solutions are then compared with experimental and numerical solutions for a benchmark problem quoted in the literature, and excellent agreements can be observed.

K(2)O(2): The most stable oxide of K

We have analyzed the stability of various oxides of K and find that K(2)O(2) is the most stable one. The additional stability is traced to the presence of oxygen dimers in K(2)O(2) which interact to form molecular orbitals. Other oxides such as KO(2) and KO(3) which also have dimers/trimers of oxygens are found to be less stable. This is traced to the shorter O-O bonds that one finds in them which gives rise to a significant coulomb repulsion between the electrons on the oxygen atoms making up the dimer/trimer, making them less stable.

Minimum error solutions of the Boltzmann equation for shock structure

‘Best’ solutions for the shock-structure problem are obtained by solving the Boltzmann equation for a rigid sphere gas by applying minimum error criteria on the Mott-Smith ansatz. The use of two such criteria minimizing respectively the local and total errors, as well as independent computations of the remaining error, establish the high accuracy of the solutions, although it is shown that the Mott-Smith distribution is not an exact solution of the Boltzmann equation even at infinite Mach number. The minimum local error method is found to be particularly simple and efficient. Adopting the present solutions as the standard of comparison, it is found that the widely used v2x-moment solutions can be as much as a third in error, but that results based on Rosen's method provide good approximations. Finally, it is shown that if the Maxwell mean free path on the hot side of the shock is chosen as the scaling length, the value of the density-slope shock thickness is relatively insensitive to the intermolecular potential. A comparison is made on this basis of present results with experiment, and very satisfactory quantitative agreement is obtained.