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.     

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.     

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.     

Polynomial-time computation of the degree of a dominant morphism in zero characteristic. II

Consider a projective algebraic variety W which is an irreducible component of the set of all common zeros of a family of homogeneous polynomials of degrees less than d in n + 1 variables over a field of zero characteristic. Consider a dominant rational morphism from W to W′ given by homogeneous polynomials of degree d′. We suggest algorithms for constructing objects in general position related to this morphism. They generalize some algorithms from the first part of the paper to the case dim W > dim W′. These algorithms are deterministic and polynomial in (dd′)n and the size of the input. Bibliography: 12 titles. __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 325, 2005, pp. 181–224.     

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     

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.