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.     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

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     

Conformational analysis of dipeptide‐derived polyisocyanides

The conformational properties of polymers derived from isocyanodipeptides have been investigated with a combination of model calculations, X‐ray diffraction, and circular dichroism spectroscopy. Depending on the configuration of the side chains, defined arrays of hydrogen bonds along the polymeric backbone are formed. This leads to a well‐defined conformation as, for example, expressed in the formation of lyotropic liquid‐crystalline phases and increased helical stability. Upon the disruption of the hydrogen bonds by a strong acid, a less well‐defined macromolecular conformation is observed. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 1725–1736, 2003