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.     

Three-dimensional Schwarzschild spacetimes

In the following we find all of the three-dimensional flat spacetimes which are static and rotationally symmetric.     

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.     

Rb saturation spectroscopy with an optically narrowed FM diode laser

A fast-frequency modulated (FM) diode laser has been optically narrowed using the technique of resonant optical feedback, to provide linewidths in each FM mode of ≈ 200 kHz peak-to-peak. With a drive frequency of 50 MHz and modulation index of 0.2, the FM laser has been used for the first time to obtain Doppler-free FM spectra of the Rb D₁ line at 795 nm. The potential use of this system for laser frequency stabilisation is discussed.     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

Hydromagnetic instability of streaming fluids in porous medium

The instability of the plane interface between two uniform, superposed, electrically conducting and counter-streaming fluids through a porous medium is considered in the presence of a horizontal magnetic field. In the absence of surface tension, perturbations transverse to the direction of streaming are found to be unaffected by the presence of streaming if perturbations in the direction of streaming are ignored. For perturbations in all other directions there exists instability for a certain wavenumber range. The instability of this system is postponed by the presence of magnetic field. The magnetic field and surface tension are able to suppress this Kelvin-Helmholtz instability for small wavelength perturbations and the medium porosity reduces the stability range given in terms of a difference between the streaming velocities and the Alfvén velocity.This research forms a part of the research project awarded to the first author (R.C.S.) by the University Grants Commission.     

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.     

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