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.     

The frequency scale of speech intonation.

In intonation research, prominence-lending pitch movements have either been described on a linear or on a logarithmic frequency scale. An experiment has been carried out to check whether pitch movements in speech intonation are perceived on one of these two scales or on a psychoacoustic scale representing the frequency selectivity of the auditory system. This last scale is intermediary between the other two scales. Subjects matched the excursion size of prominence-lending pitch movements in utterances resynthesized in different pitch registers. Their task was to adjust the excursion size in a comparison stimulus in such a way that it lent equal prominence to the corresponding syllable in a fixed test stimulus. The comparison stimulus and the test stimulus had pitches running parallel on either the logarithmic frequency scale, the psychoacoustic scale, or the linear frequency scale. In one-half of the experimental sessions, the test stimulus was presented in the low register, while the comparison stimulus was presented in the high register, and, conversely, for the other half of the sessions. The result is that, in all cases, stimuli are matched in such a way that the average excursion sizes in different registers are equal on the psychoacoustic scale.     

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.     

BRST cohomology and highest weight vectors. I

We initiate a program to study certain recent problems in non-compact coset CFT by the BRST approach. We derive a reduction formula for the BRST cohomology by making use of a twisting by highest weight modules. As illustrations, we apply the formula to the bosonic string model and a rank one non-compact coset model [DPL]. Our formula provides a completely new approach to non-compact coset construction.Partially supported by NSF Grant DMS-8703581     

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     

Projection of Markov Measures May Be Gibbsian

We study the induced measure obtained from a 1-step Markov measure, supported by a topological Markov chain, after the mapping of the original alphabet onto another one. We give sufficient conditions for the induced measure to be a Gibbs measure (in the sense of Bowen) when the factor system is again a topological Markov chain. This amounts to constructing, when it does exist, the induced potential and proving its Hölder continuity. This is achieved through a matrix method. We provide examples and counterexamples to illustrate our results.