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.     

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.     

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     

Effects of sodium salicylate on the febrile response and increased levels of cyclic AMP in cerebrospinal fluid during endogenous pyrogen-induced fever in rabbits

To further evaluate the causality between endogenous pyrogen (EP)-induced fever and cyclic adenosine-3',5'-monophosphate (cyclic AMP) level, the effects of sodium salicylate (SS) on the febrile response and increased levels of cyclic AMP in both cerebrospinal fluid (c.s.f.) and plasma during EP-induced fever in rabbits were observed. The results suggest that cyclic AMP is probably involved in the central mediation of EP-induced fever and that increased concentration of cyclic AMP in c.s.f. associated with EP-induced fever is not the result of temperature elevation but appears to be caused by the increased synthesis in the central nervous system. In addition it is confirmed that blood is impossibly a contributory source of increased cyclic AMP in c.s.f. during EP fever, and that SS may act subsequent to the increase in cyclic AMP.     

Surface-field induced interband tunneling in InAs

Metal/insulator/semiconductor junctions are prepared on degeneratep-type InAs substrates with hole concentrations ranging from 2.3×10¹⁷ cm^–3 to 2.7×10¹⁸ cm^–3. The low work function of the top metal Yb, Al, or Au and charged interface states influence a two-dimensional (2D) electron inversion layer at the InAs surface. The insulator barrier that is formed by thermal oxidation is designed sufficiently thin, so that the bias voltage applied at the metal electrode mainly drops across the depletion layer separating the electron channel from the bulk. The current-voltage (I–V) characteristics exhibit strong negative differential conductance due to interband, tunneling from the 2D subband into the 3D valence band with peak-to-valley current ratios up to 3.1, 18, and 32 at 300 K, 77 K, and 4.2 K, respectively. In agreement with a theoretical model based on coherentelastic tunneling, the form of the I–V curves resembles those of double-barrier resonant tunnel devices rather than those of 3D Esaki diodes. The series resistance is obtained from the saturation of the differential conductance dI/dV at high forward bias and from the shift of structures in d² I/dV ² arising from phonon assisted tunneling.Dedicated to G. Lautz on the occasion of his 65th birthday