Existence and Stability for Spherical Crystals Growing in a Supersaturated Solution

We consider the growth of a spherical crystal in a supersaturatedsolution. In the first part, existence and uniqueness resultsfor radially symmetric growth are obtained, provided that thesupersaturation is not too large; conversely, when the far-fieldsupersaturation exceeds a critical value, it is shown that theradially symmetric solution ceases to exist in finite time.In the second part, we examine the linear stability of a radiallysymmetric similarity solution (in which the radius grows ast?) to shape perturbations. The results are compared with previousquasi-static analyses, and, in particular, the critical radiusat which the crystal becomes unstable is found to be largerfor small supersaturations, but smaller for large supersaturations,than those predicted by the quasi-static analysis     

Dependence of the H2-H2 Interaction on the Monomer Bond Lengths: Steps Toward an Accurate ab initio Estimate

We compute the vibrational coupling between two H₂ molecules from ab initio quantum chemical calculations of the H₂-H₂ potential carried out at the full configuration interaction level of theory using the atom-centered aug-cc-pVTZ basis set for hydrogen. We compare the full configuration interaction results with those obtained using two variants of coupled cluster theory and find that a fully iterative treatment of connected triples may be required to estimate the H₂-H₂ vibrational coupling accurately using coupled cluster theory.     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Lateral interactions in the thermal desorption of H2 from clean Cu/Ni(110) alloy surfaces

The adsorption of hydrogen on a clean Cu10%/Ni90% (110) alloy single crystal was studied using flash desorption spectroscopy (FDS), Auger electron spectroscopy (AES), and work function measurements. Surface compositions were varied from 100% Ni to 35% Ni. The hydrogen chemisorption on a-surface of 100% nickel revealed strong attractive interactions between the hydrogen atoms in accordance with previous work on Ni(100). Three desorption states (β₁, β₂ and α) appeared in the desorption spectra. The highest temperature (α) state was occupied only after the initial population of the β₂-state. As the amount of copper was increased in the nickel substrate, desorption from the higher energy binding α-state was reduced, indicating a decrease in the attractive interactions among hydrogen atoms. The hydrogen coverage at saturation was not affected by the addition of copper to the nickel substrate until the copper concentration was greater than 25% at which a sharp reduction in saturation coverage occurred. This phenomenon was apparently due to the adsorption of hydrogen on Ni atoms followed by occupation of NiNi and CuNi bridged adsorption sites, while occupation of CuCu sites was restricted due to an energy barrier to migration.     

OPTIMAL INVESTMENT IN THE PROTECTION OF A VULNERABLE BIOLOGICAL RESOURCE

It is assumed that the probability of destruction of a biological asset by natural hazards can be reduced through investment in protection. Specifically a model, in which the hazard rate depends on both the age of the asset and the accumulated invested protection capital, is assumed. The protection capital depreciates through time and its effectiveness in reducing the hazard rate is subject to diminishing returns. It is shown how the investment schedule to maximize the expected net present value of the asset can be determined using the methods of deterministic optimal control, with the survival probability regarded as a state variable. The optimal investment pattern involves “bang-bang-singular” control. A numerical scheme for determining jointly the optimal investment policy and the optimal harvest (or replacement) age is outlined and a numerical example involving forest fire protection is given.     

Double-layer structural-acoustic coupling for cylindrical shell by using a combination of WDA and BIE

This work formulates the double-layer structural-acoustic coupling problem for cylindrical shell by using a combination of the wave-number domain approach (WDA) and the boundary integral equation (BIE). Expressions for the spectral radial velocity of the outer surface of a finite fluid-filled/submerged (FFS) cylindrical thin shell are formulated by means of the transfer matrix equation in wave-number domain. It is shown that the spectral variables on the inner surface of the shell are related to those on the outer surface of the shell. The far field sound radiation from this kind of shell is numerically evaluated for various fluid cases. An experimental verification is performed, and a good correlation between the theoretical results and the experimental results shows that the theoretical study work in this paper is correct.     

Miscibility of bisphenol-A polycarbonate with poly(methyl methacrylate)

The miscibility of bisphenol-A polycarbonate (PC) with poly(methyl methacrylate) (PMMA) has been reexamined using differential scanning calorimetry (DSC) and optical indications for phase separation on heating, i.e., lower critical solution temperature (LCST) behavior. Various methods have been used to prepare the blends including methylene chloride (CH₂Cl₂) and tetrahydrofuran (THF) solution casting, melt mixing, and precipitation of PC and PMMA simultaneously from THF solution by using the nonsolvents methanol and heptane. It is shown that the resulting phase behavior for PC/PMMA blends is strongly affected by the blend preparation method. However, these blends are miscible over the whole blend composition range (unambiguous single composition-dependent T_g's and LCST behavior) when prepared by precipitation from solution using heptane as the nonsolvent. To the contrary, solution-cast and melt-mixed PC/PMMA blends were all phase separated, which may be attributed to the “solvent” effect and LCST behavior, respectively, not discovered in previous reports. Methanol precipitation does not lead to fully mixed blends, which demonstrates the importance of the choice of nonsolvent when using the precipitation method.     

Semi‐analytical integration of the 8‐node plane element stiffness matrix using symbolic computation

The semi‐analytical integration of an 8‐node plane strain finite element stiffness matrix is presented in this work. The element is assumed to be super‐parametric, having straight sides. Before carrying out the integration, the integral expressions are classified into several groups, thus avoiding duplication of calculations. Symbolic manipulation and integration is used to obtain the basic formulae to evaluate the stiffness matrix. Then, the resulting expressions are postprocessed, optimized, and simplified in order to reduce the computation time. Maple symbolic‐manipulation software was used to generate the closed expressions and to develop the corresponding Fortran code. Comparisons between semi‐analytical integration and numerical integration were made. It was demonstrated that semi‐analytical integration required less CPU time than conventional numerical integration (using Gaussian‐Legendre quadrature) to obtain the stiffness matrix. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006