Method of lines with boundary elements for 1‐D transient diffusion‐reaction problems

Time‐dependent differential equations can be solved using the concept of method of lines (MOL) together with the boundary element (BE) representation for the spatial linear part of the equation. The BE method alleviates the need for spatial discretization and casts the problem in an integral format. Hence errors associated with the numerical approximation of the spatial derivatives are totally eliminated. An element level local cubic approximation is used for the variable at each time step to facilitate the time marching and the nonlinear terms are represented in a semi‐implicit manner by a local linearization at each time step. The accuracy of the method has been illustrated on a number of test problems of engineering significance. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

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.     

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.     

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.