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.     

Higher-dimensional homogeneous cosmological models as autonomous dynamical systems

The field equations for homogeneous models in an arbitrary number of dimensions form a Hamiltonian system with constraint forces. Due to the monotonic behavior of the determinant of the induced metric, the evolution of the system can be interpreted as the motion of a particle in an explicitly time-dependent potential. Considering vacuum models, we show that this explicit time-dependence can be eliminated. Using the scaling properties of the Ricci tensor we obtain an autonomous system, for which we can also find a Liapunov function in terms of the n-dimensional Ricci curvature.     

Gluon jets from quarkonia

A detailed analysis is made within perturbative QCD of the decays into gluon jets of P-wave and S-wave quarkonia (J^PC = 0⁺⁺, 1⁺⁺, 2⁺⁺, 1^+? and 0^?+, 1^??). For all those states we compute to leading non-trivial order in α_s the conventional fraction f(?, δ) of the jet-like hadronic final states, characterized by the back-to-back cone of half angle δ, in which all but a small fraction ? of the total energy is emitted. To the same order we give also the average values of thrust momenta and spherocity distributions.     

Determination of trace levels of diosmin in a pharmaceutical preparation by adsorptive stripping voltammetry at a glassy carbon electrode.

A systematic study on the electrochemical behavior of diosmin in Britton-Robinson buffer (pH 2.0-10.0) at a glassy carbon electrode (GCE) was made. The oxidation process of the drug was found to be quasi-reversible with an adsorption-controlled step. The adsorption stripping response was evaluated with respect to various experimental conditions, such as the pH of the supporting electrolyte, the accumulation potential and the accumulation time. The observed anodic peak current at +0.73 V vs. Ag/AgCl reference electrode increased linearly over two orders of magnitude from 5.0x10(-8) M to 9.0x10(-6) M. A limit of detection down to 3.5x10(-8) M of diosmin at the GCE was achieved with a mean recovery of 97+/-2.1%. Based on the electrochemical data, an open-circuit accumulation step in a stirred sample solution of BR at pH 3.0 was developed. The proposed method was successfully applied to the determination of the drug in pharmaceutical formulations. The results compared favorably with the data obtained via spectrophotometric and HPLC methods.     

A MODEL OF TROPICAL MARINE RESERVE‐FISHERY LINKAGES

ABSTRACT. The excessive and unsustainable exploitation of our marine resources has led to the promotion of marine reserves as a fisheries management tool. Marine reserves, areas in which fishing is restricted or prohibited, can offer opportunities for the recovery of exploited stock and fishery enhancement. In this paper we examine the contribution of fully protected tropical marine reserves to fishery enhancement by modeling marine reserve‐fishery linkages. The consequences of reserve establishment on the long‐run equilibrium fish biomass and fishery catch levels are evaluated. In contrast to earlier models this study highlights the roles of both adult (and juvenile) fish migration and larval dispersal between the reserve and fishing grounds by employing a spawner‐recruit model. Uniform larval dispersal, uniform larval retention and complete larval retention combined with zero, moderate and high fish migration scenarios are analyzed in turn. The numerical simulations are based on Mombasa Marine National Park, Kenya, a fully protected coral reef marine reserve comprising approximately 30% of former fishing grounds. Simulation results suggest that the establishment of a fully protected marine reserve will always lead to an increase in total fish biomass. If the fishery is moderately to heavily exploited, total fishery catch will be greater with the reserve in all scenarios of fish and larval movement. If the fishery faces low levels of exploitation, catches can be optimized without a reserve but with controlled fishing effort. With high fish migration from the reserve, catches are optimized with the reserve. The optimal area of the marine reserve depends on the exploitation rate in the neighboring fishing grounds. For example, if exploitation is maintained at 40%, the ‘optimal’ reserve size would be 10%. If the rate increases to 50%, then the reserve needs to be 30% of the management area in order to maximize catches. However, even in lower exploitation fisheries (below 40%), a small reserve (up to 20%) provides significantly higher gains in fish biomass than losses in catch. Marine reserves are a valuable fisheries management tool. To achieve maximum fishery benefits they should be complemented by fishing effort controls.     

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     

Nonlinear Theory of Spectra of Thermally Stimulated Currents in Complex Crystals with Hydrogen Bonds

The nonlinear theory of thermally stimulated depolarization currents is developed. The theory explains the processes of hetero- and homocharge relaxation in complex crystals with hydrogen bonds and allows the relaxation oscillator parameters to be calculated using the quadratic approximation for the external electric field.