The Wiener-Hopf solution of a class of mixed boundary value problems arising in surface water wave phenomena

Two mixed boundary value problems associated with two-dimensional Laplace equation, arising in the study of scattering of surface waves in deep water (or interface waves in two superposed fluids) in the linearised set up, by discontinuities in the surface (or interface) boundary conditions, are handled for solution by the aid of the Weiner-Hopf technique applied to a slightly more general differential equation to be solved under general boundary conditions and passing on to the limit in a manner so as to finally give rise to the solutions of the original problems. The first problem involves one discontinuity while the second problem involves two discontinuities. The reflection coefficient is obtained in closed form for the first problem and approximately for the second. The behaviour of the reflection coefficient for both the problems involving deep water against the incident wave number is depicted in a number of figures. It is observed that while the reflection coefficient for the first problem steadily increases with the wave number, that for the second problem exhibits oscillatory behaviour and vanishes at some discrete values of the wave number. Thus, there exist incident wave numbers for which total transmission takes place for the second problem. (C) 1999 Elsevier Science B.V. All rights reserved.

Automatic decomposition of unstructured meshes employing genetic algorithms for parallel FEM computations

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

Mechanistic investigations on the efficient catalytic decomposition of peroxynitrite by ebselen analogues

In this study, ebselen and its analogues are shown to be catalysts for the decomposition of peroxynitrite (PN). This study suggests that the PN-scavenging ability of selenenyl amides can be enhanced by a suitable substitution at the phenyl ring in ebselen. Detailed mechanistic studies on the reactivity of ebselen and its analogues towards PN reveal that these compounds react directly with PN to generate highly unstable selenoxides that undergo a rapid hydrolysis to produce the corresponding seleninic acids. The selenoxides interact with nitrite more effectively than the corresponding seleninic acids to produce nitrate with the regeneration of the selenenyl amides. Therefore, the amount of nitrate formed in the reactions mainly depends on the stability of the selenoxides. Interestingly, substitution of an oxazoline moiety on the phenyl ring stabilizes the selenoxide, and therefore, enhances the isomerization of PN to nitrate.

On the explicit determination of the polar decomposition in n-dimensional vector spaces

A method for the explicit determination of the polar decomposition (and the related problem of finding tensor square roots) when the underlying vector space dimension n is arbitrary (but finite), is proposed. The method uses the spectral resolution, and avoids the determination of eigenvectors when the tensor is invertible. For any given dimension n, an appropriately constructed van der Monde matrix is shown to play a key role in the construction of each of the component matrices (and their inverses) in the polar decomposition.

Oxidation and decomposition of NH3 over combustion synthesized Al2O3 and CeO2 supported Pt, Pd and Ag catalysts

The catalytic oxidation and decomposition of NH3 have been carried out over combustion synthesized Al2O3 and CeO2 supported Pt, Pd and Ag catalysts using temperature programmed reaction (TPR) technique in a packed bed tubular reactor. Metals are ionically dispersed over CeO2 and fine metal particles are found on Al2O3. NH3 oxidation occurs over 1% Pt/Al2O3, 1% Pd/Al2O3 and 1% Ag/Al2O3 at 175, 270 and 350 C respectively producing N-2, NO, N2O and H2O, whereas 1% Pt/CeO2, 1% Pd/CeO2 and 1% Ag/CeO2 give N-2 along with NO, N2O and H2O at 200, 225 and 250degreesC respectively. N-2 predominates over other nitrogen-containing products during the reaction on all catalysts. At less O-2 concentration, N-2 and H2O are the only products obtained during NH3 Oxidation. NH3 decomposition over all the catalysts occurs above 450degreesC.

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.

Unimolecular HCl elimination from 1,2-dichloroethane: A single pulse shock tube and ab initio study

Thermal decomposition of 1,2-dichloroethane (1,2-DCE) has been studied in the temperature range of 10501175 K behind reflected shock waves in a single pulse shock tube. The unimolecular elimination of HCl is found to be the major channel through which 1,2-DCE decomposes under these conditions. The rate constant for the unimolecular elimination of HCl from 1,2-dichloroethane is found to be 10(13.98+/-0.80) exp(-57.8+/-2.0/RT) s(-1), where the activation energy is given in kcal mol(-1) and is very close to that value for CH3CH2Cl (EC). Ab initio (HF and MP2) and DFT calculations have been carried out to find the activation barrier and the structure of the transition state for this reaction channel from both EC and 1,2-DCE. The preexponential factors calculated at various levels of theory (BF/6-311++G**, MP2/6-311++G**, and B3LYP/6-311++G**) are (approximate to10(15) s(-1)) significantly larger than the experimental results. If the torsional mode in the ground state is treated as free internal rotation the preexponential factors reduce significantly, giving excellent agreement with experimental values. The DFT results are in excellent (fortuitous?) agreement with the experimental value for activation energy for 1,2-DCE while the MP2 and HF results seem to overestimate the barrier. However, DFT results for EC is 4.5 kcal mol(-1) less than the previously reported experimental values. At all levels, theory predicts an increase in HCI elimination barrier on beta-Cl substitution on EC.

A new parallel overlapped domain decomposition method for nonlinear dynamic finite element analysis

In this paper a new parallel algorithm for nonlinear transient dynamic analysis of large structures has been presented. An unconditionally stable Newmark-beta method (constant average acceleration technique) has been employed for time integration. The proposed parallel algorithm has been devised within the broad framework of domain decomposition techniques. However, unlike most of the existing parallel algorithms (devised for structural dynamic applications) which are basically derived using nonoverlapped domains, the proposed algorithm uses overlapped domains. The parallel overlapped domain decomposition algorithm proposed in this paper has been formulated by splitting the mass, damping and stiffness matrices arises out of finite element discretisation of a given structure. A predictor-corrector scheme has been formulated for iteratively improving the solution in each step. A computer program based on the proposed algorithm has been developed and implemented with message passing interface as software development environment. PARAM-10000 MIMD parallel computer has been used to evaluate the performances. Numerical experiments have been conducted to validate as well as to evaluate the performance of the proposed parallel algorithm. Comparisons have been made with the conventional nonoverlapped domain decomposition algorithms. Numerical studies indicate that the proposed algorithm is superior in performance to the conventional domain decomposition algorithms. (C) 2003 Elsevier Ltd. All rights reserved.

Decomposition temperatures of Cu2Ln2O5 (Ln = Tb, Dy, Ho, Er, Tm, Yb, and Lu) compounds

The equilibrium decomposition temperatures of Cu2Ln2O5 (Ln = Tb, Dy, Ho, Er, Tm, Yb, and Lu) compounds have been measured using a combined DTA-TGA apparatus under a flowing Ar + O2 gas mixture, in which the partial pressure of oxygen was controlled at 5.0 × 103 Pa. The Cu2Ln2O5 compounds yield Ln2O3 and Cu2O on decomposition. The decomposition temperature increases monotonically with the atomic number of the lanthanide element. This suggests that the stability of the Cu2Ln2O5 compounds with respect to the component binary oxides increases with decreasing radius of the Ln3+ ion.

Kinetic decomposition of Ni2SiO4 in oxygen potential gradients

Nickel orthosilicate (Ni2SiO4) has been found to decompose into its component binary oxides in oxygen potential gradients at 1373 K. Nickel oxide was formed at the high oxygen potential boundary, while silica was detected at the low oxygen potential side. Significant porosity and fissures were observed near the Ni2SiO4/SiO2 interface and the SiO2 layer. The critical oxygen partial pressure ratio required for decomposition varied from 1.63 to 2.15 as the oxygen pressures were altered from 1.01 ⊠ 105 to 2.7X 10−4 Pa, well above the dissociation pressure of Ni2SiO4. Platinum markers placed at the boundaries of the Ni2SiO4 sample indicated growth of NiO at the higher oxygen potential boundary, without any apparent transport of material to the low oxygen potential side. However, significant movement of the bulk Ni2SiO4 crystal with respect to the marker was not observed. The decomposition of the silicate occurs due to the unequal rates of transport of Ni and Si. The critical oxygen partial pressure ratio required for decomposition is related both to the thermodynamic stability of Ni2SiO4 with respect to component oxides and the ratio of diffusivities of nickel and silicon. Kinetic decomposition of multicomponent oxides, first discovered by Schmalzried, Laqua, and co-workers [H. Schmalzried, W. Laqua, and P. L. Lin, Z. Natur Forsch. Teil A 34, 192 (1979); H. Schmalzried and W. Laqua, Oxid. Met. 15, 339 (1981); W. Laqua and H. Schmalzried, Chemical Metallurgy—A Tribute to Carl Wagner (Metallurgical Society of the AIME, New York, 1981), p. 29] has important consequences for their use at high temperatures and in geochemistry.

Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

An investigation of carbon nanotubes obtained from the decomposition of methane over reduced Mg1− xM xAl2O4 spinel catalysts

Carbon nanotubes produced by the treatment of Mg1−xMxAl2O4 (M = Fe, Co, or Ni; x = 0.1, 0.2, 0.3, or 0.4) spinels with an H2–CH4 mixture at 1070 °C have been investigated systematically. The grains of the oxide-metal composite particles are uniformly covered by a weblike network of carbon nanotube bundles, several tens of micrometers long, made up of single-wall nanotubes with a diameter close to 4 nm. Only the smallest metal particles (<5 nm) are involved in the formation of the nanotubes. A macroscopic characterization method involving surface area measurements and chemical analysis has been developed in order to compare the different nanotube specimens. An increase in the transition metal content of the catalyst yields more carbon nanotubes (up to a metal content of 10.0 wt% or x = 0.3), but causes a decrease in carbon quality. The best compromise is to use 6.7 wt% of metal (x = 0.2) in the catalyst. Co gives superior results with respect to both the quantity and quality of the nanotubes. In the case of Fe, the quality is notably hampered by the formation of Fe3C particles.

Determining Top K Nodes in Social Networks using the Shapley Value

In this paper, we consider the problem of selecting, for any given positive integer k, the top-k nodes in a social network, based on a certain measure appropriate for the social network. This problem is relevant in many settings such as analysis of co-authorship networks, diffusion of information, viral marketing, etc. However, in most situations, this problem turns out to be NP-hard. The existing approaches for solving this problem are based on approximation algorithms and assume that the objective function is sub-modular. In this paper, we propose a novel and intuitive algorithm based on the Shapley value, for efficiently computing an approximate solution to this problem. Our proposed algorithm does not use the sub-modularity of the underlying objective function and hence it is a general approach. We demonstrate the efficacy of the algorithm using a co-authorship data set from e-print arXiv (www.arxiv.org), having 8361 authors.