Coupled lubrication and Stokes flow finite elements

A method is developed for performing a local reduction of the governing physics for fluid problems with domains that contain a combination of narrow and non‐narrow regions, and the computational accuracy and performance of the method are measured. In the narrow regions of the domain, where the fluid is assumed to have no inertia and the domain height and curvature are assumed small, lubrication, or Reynolds, theory is used locally to reduce the two‐dimensional Navier–Stokes equations to the one‐dimensional Reynolds equation while retaining a high degree of accuracy in the overall solution. The Reynolds equation is coupled to the governing momentum and mass equations of the non‐narrow region with boundary conditions on the mass and momentum flux. The localized reduction technique, termed ‘stitching,’ is demonstrated on Stokes flow for various geometries of the hydrodynamic journal bearing—a non‐trivial test problem for which a known analytical solution is available. The computational advantage of the coupled Stokes–Reynolds method is illustrated on an industrially applicable fully‐flooded deformable‐roll coating example. The examples in this paper are limited to two‐dimensional Stokes flow, but extension to three‐dimensional and Navier–Stokes flow is possible. Copyright © 2003 John Wiley & Sons, Ltd.     

Kinetics of phenyl radical reactions with propane,n‐butane,n‐hexane,and n‐octane: Reactivity of C6H5 toward the secondary C?H bond of alkanes

The kinetics of C₆H₅ reactions with n‐C_nH_2n+2 (n = 3, 4, 6, 8) have been studied by the pulsed laser photolysis/mass spectrometric method using C₆H₅COCH₃ as the phenyl precursor at temperatures between 494 and 1051 K. The rate constants were determined by kinetic modeling of the absolute yields of C₆H₆ at each temperature. Another major product C₆H₅CH₃ formed by the recombination of C₆H₅ and CH₃ could also be quantitatively modeled using the known rate constant for the reaction. A weighted least‐squares analysis of the four sets of data gave k (C₃H₈) = (1.96 ± 0.15) × 10¹¹ exp[?(1938 ± 56)/T], and k (n‐C₄H₁₀) = (2.65 ± 0.23) × 10¹¹ exp[?(1950 ± 55)/T] k (n‐C₆H₁₄) = (4.56 ± 0.21) × 10¹¹ exp[?(1735 ± 55)/T], and k (n?C₈H₁₈) = (4.31 ± 0.39) × 10¹¹ exp[?(1415 ± 65)T] cm³ mol^?1 s^?1 for the temperature range studied. For the butane and hexane reactions, we have also applied the CRDS technique to extend our temperature range down to 297 K; the results obtained by the decay of C₆H₅ with CRDS agree fully with those determined by absolute product yield measurements with PLP/MS. Weighted least‐squares analyses of these two sets of data gave rise to k (n?C₄H₁₀) = (2.70 ± 0.15) × 10¹¹ exp[?(1880 ± 127)/T] and k (n?C₆H₁₄) = (4.81 ± 0.30) × 10¹¹ exp[?(1780 ± 133)/T] cm³ mol^?1 s^?1 for the temperature range 297‐‐1046 K. From the absolute rate constants for the two larger molecular reactions (C₆H₅ + n‐C₆H₁₄ and n‐C₈H₁₈), we derived the rate constant for H‐abstraction from a secondary C? H bond, k_s?CH = (4.19 ± 0.24) × 10¹⁰ exp[?(1770 ± 48)/T] cm³ mol^?1 s^?1. © 2003 Wiley Periodicals, Inc. Int J Chem Kinet 36: 49–56, 2004     

Kinetics and mechanism of the reactions of hexaaqua rhodium (III) with sulphur (IV) in aqueous medium

An O-bonded sulphito complex, Rh(OH₂)₅(OSO₂H)²⁺, is reversibly formed in the stoppedflow time scale when Rh(OH₂) ₆ ³⁺ and SO₂/HSO ₃ ⁻ buffer (1 <pH< 3) are allowed to react. For Rh(OH₂)₅OH₂₊+ SO₂ □ Rh(OH₂)₅(OSO₂H)²⁺ (k₁/k_-1), k₁ = (2.2 ±0.2) × 10³ dm³ mol⁻¹ s⁻¹, k₋1 = 0.58 ±0.16 s⁻¹ (25°C,I = 0.5 mol dm⁻³). The protonated O-sulphito complex is a moderate acid (K _d = 3 × 10⁻⁴ mol dm⁻³, 25°C, I= 0.5 mol dm⁻³). This complex undergoes (O, O) chelation by the bound bisulphite withk= 1.4 × 10⁻³ s⁻¹ (31°C) to Rh(OH₂)₄(O₂SO)⁺ and the chelated sulphito complex takes up another HSO ₃ ⁻ in a fast equilibrium step to yield Rh(OH₂)₃(O₂SO)(OSO₂H) which further undergoes intramolecular ligand isomerisation to the S-bonded sulphito complex: Rh(OH₂)₃(O₂SO)(OSO₂)^- → Rh(OH₂)₃(O₂SO)(SO₃)⁻ (k _iso = 3 × 10⁻⁴ s⁻¹, 31°C). A dinuclear (μ-O, O) sulphite-bridged complex, Na₄[Rh₂(μ-OH)₂(OH)₂(μ-OS(O)O)(O₂SO)(SO₃) (OH₂)]5H₂O with (O, O) chelated and S-bonded sulphites has been isolated and characterized. This complex is sparingly soluble in water and most organic solvents and very stable to acid-catalysed decomposition     

Ultraviolet photoemission spectroscopy of hydrogen complex deactivation on InP:Zn(1 0 0) surfaces

Ultraviolet photoemission spectroscopy is used to study the kinetics of the H-Zn complex deactivation in Zn doped InP(1 0 0). Hydrogen injected into the material electronically passivates the local carrier concentration. Reverse-biased anneals of the InP under ultra-high vacuum show a dramatic change in the work function of the material with increasing temperature. Spectral features are also shown to be sensitive to sample temperature. To our knowledge, we show the first view of hydrogen retrapping at the surface using photoemission spectroscopy. A simple photoelectron threshold energy analysis shows the state of charge compensation of the material.     

High pressure elastic and plastic deformations of silica: In situ diamond anvil cell Raman experiments

Normal silica glass is usually referred to as low density amorphous silica as it can be converted to high density amorphous silica by a hydrostatic pressure (polyamorphic transition). In this work in situ Raman experiments are performed in a diamond anvil cell up to 18 GPa. The pressure effects on the structure of silica after successive compression decompression experiments are analyzed. The mode Grüneisen parameters corresponding to the elastic compression of high density amorphous silica are obtained and compared with those of normal silica. A reorganization of the high density amorphous silica below 3 GPa is evidenced.     

Fluorine incorporation in silica glass by the MCVD process: Study of fluorine incorporation zone,evaluation of optical properties and structure of the glass

A theoretical approach was made to find out a complete fluorine incorporation zone on a ternary diagram which serves as a useful graphical representation to select the flows of the supplied reagents for incorporation of the suitable amount of fluorine into cladding glass of optical fiber preform made by the MCVD process using CCl₂F₂ as a source of fluorine under oxygen abundance, oxygen deficiency and intermediate oxygen state conditions. The possible mechanism for incorporation of fluorine into cladding glass of optical fiber is also evaluated on the basis of the thermodynamical data. The fluorine incorporation mechanism in silica glass by the MCVD process is found to be dependent on the CCl₂F₂/SiCl₄ ratio in the input gas mixture. Fluorine doping is found to be effective for removing the strained Si–O–Si bonds, which govern the optical transparency in deep ultra-violet (DUV) and vacuum ultra-violet (VUV) regions. The maximum refractive index depression of ?0.5 × 10^?3 is obtained with incorporation of fluorine into silica cladding glass by the MCVD process using CCl₂F₂ as a dopant precursor with suitable flow of SiCl₄ vapor along with O₂ through backward deposition pass. The structure of fluorine doped silica glass preform samples containing 1.70–1.79 mol% fluorine incorporated by the MCVD process based on the analyses of ¹⁹F MAS spectra done by high-resolution ¹⁹F NMR spectroscopy reveal the presence of two distinct types of fluorine environments. The majority of the fluorine environments are formed in SiO_1.5F polyhedral and less abundant species is observed to be highly unusual, yielding a fivefold coordinated silicon of the type SiO₂F polyhedral which become increased with increasing the fluorine content.     

A semi-analytic meshfree method for Almansi–Michell problems of piezoelectric cylinders

Saint-Venant’s Problem, Almansi–Michell Problems, Meshfree Methods, Piezoelectricity. We present a semi-analytical method for analyzing prismatic nonhomogeneous piezoelectric cylinders with arbitrary cross-sectional geometry. The prescribed loads considered in this study include axial forces, torques, moments, and voltage resultants prescribed at the cylinder’s ends, as well as body forces, lateral surface shears, voltages, and pressures as long as they can be represented by a power series in the axial coordinate. This problem can be considered as an extension of Saint-Venant and Almansi–Michell problems for elastic bodies to piezoelectric bodies. In this computationally efficient method, the cross-sectional plane is discretized with a meshfree approach, and the solution is obtained analytically with respect to the axial coordinate. A number of examples are provided to demonstrate the veracity and utility of the proposed method.     

Ductile–brittle transitions in the fracture of plastically deforming,adhesively bonded structures. Part II: Numerical studies

To enable the effective and reliable use of structural adhesive bonding in automotive applications, the cohesive properties of a joint need to be determined over a wide range of loading rates. In this paper, a strategy for determining these properties has been described and used to analyze a set of experimental results presented in a companion paper. In the particular system studied, a crack growing in a toughened quasi-static mode could make a catastrophic transition to a brittle mode of fracture. The cohesive parameters for both the toughened and brittle modes of crack growth were determined by comparing numerical predictions from cohesive-zone simulations to the results of experimental tests performed using double-cantilever beam specimens and tensile tests. The cohesive parameters were found to be essentially rate-independent for the toughened mode, but the toughness dropped by a factor of four upon a transition to the brittle mode. The results of wedge tests were used as an independent verification of the cohesive parameters, and to verify that the quasi-static properties remained rate-independent to very high crack velocities corresponding to conditions of low-velocity impact. The effects of friction, and the use of the wedge test to determine cohesive parameters, were also explored.     

Atomistic phenomena in dense fluid shock waves

The shock structure problem is one of the classical problems of fluid mechanics and at least for non-reacting dilute gases it has been considered essentially solved. Here we present a few recent findings, to show that this is not the case. There are still new physical effects to be discovered provided that the numerical technique is general enough to not rule them out a priori. While the results have been obtained for dense fluids, some of the effects might also be observable for shocks in dilute gases.