Thermoelastoplastic deformation of a thick-walled cylinder with a radial crack

The thermoelastoplastic fracture mechanics problem of a thick-walled cylinder subjected to internal pressure and a nonuniform temperature field is solved by the method of elastic solutions combined with the finite-element method. The correctness of the solution is provided by using the Barenblatt crack model, in which the stress and strain fields are regular. The elastoplastic problem of a cracked cylinder subjected to internal pressure and a nonuniform temperature field are solved. The calculation results are compared with available data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 173–183, May–June, 2008.     

Axisymmetric frictionless contact of functionally graded materials

The main interest of this study is a new method to solve the axisymmetric frictionless contact problem of functionally graded materials (FGMs). Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into a series of sub-layers with shear modulus varying linearly in each sub-layer and continuous at the sub-interfaces. With this model, the axisymmetric frictionless contact problem of a functionally graded coated half-space is investigated. By using the transfer matrix method and Hankel integral transform technique, the problem is reduced to a Cauchy singular integral equation. The contact pressure, contact region and indentation are calculated for various indenters by solving the equations numerically. An erratum to this article can be found at     

Water wave radiation by a sphere submerged in water with an ice-cover

Using the multipoles method, we formulate the problems of radiation (both heave and sway) of water waves by a submerged sphere in deep as well as in uniform finite depth water with an ice-cover, with the ice-cover being modelled as an elastic plate of very small thickness. In each case this leads to an infinite system of linear equations which are solved numerically by standard techniques. The added-mass and damping coefficients for a heaving and swaying sphere are obtained and depicted graphically against the wave number for various values of the radius of the submerged sphere and flexural rigidity of the ice-cover to show the effect of the presence of ice-cover on these quantities. When the flexural rigidity is taken to be zero, the numerical results for the added-mass and damping coefficient for water with a free surface are recovered.     

Numerical method and models for anti-plane strain of a system with a thin elastic wedge

The paper presents a general method to find asymptotics for a (multi-)wedge system containing a thin wedge. It employs separation of the symmetric and anti-symmetric parts of the boundary displacements and tractions of the wedge. The method is applicable when the angle of the thin wedge turns to zero. A physical interpretation of the derived equations is obtained by using power expansions of non-polynomial functions, which appear after the Mellin transform. We establish that the first term in the expansion of the symmetric part corresponds to shear, while the first term of the anti-symmetric part describes deflection of the wedge axis. Numerical experiments, performed by using a code developed on the basis of the theory, show that using only the first terms of the expansions insignificantly influence accuracy: the approximate results coincide with the exact values of roots to the third significant digit even for the wedge angle of 30°.     

The comparison of detailed chemical kinetic mechanisms: Application to the combustion of methane

Detailed chemical kinetic mechanisms of complex chemical phenomena may be composed of hundreds of species and thousands of individual elementary reactions. It can be an extremely laborious and error‐prone procedure to compare two of these mechanisms, particularly if they come from different sources. We have created software tools which help to highlight the differences between mechanisms written in a Chemkin format and demonstrate their applicability to five literature mechanisms describing the high temperature oxidation of methane. © 2004 Wiley Periodicals, Inc. Int J Chem Kinet 36:467–471, 2004     

Large-scale computer simulation of fully developed turbulent channel flow with heat transfer

Recently, with the advent of supercomputers, there has been considerable interest in the use of direct numerical simulation to obtain information about turbulent shear flow at low Reynolds number. This paper presents a pseudospectral technique to solve the full three-dimensional time-dependent Navier-Stokes and advection-diffusion equations without the use of subgrid-scale modelling. The technique has not been previously used for fully developed turbulent channel flow simulation and is based on methods applied in other contexts. The emphasis of this paper is to provide a reasonably detailed account of how the simulation is done rather than to present new calculations of turbulence. The details of an algorithm for turbulent channel flow simulation and the grid and time step sizes needed to integrate through transient behaviour to steady state turbulence have not been published before and are presented here. Results from a Cray-2 simulation of fully developed turbulent flow in a channel with heat transfer are presented along with a critical comparison between experiment and computation. The first- and second-order moments agree well with experimental measurements; the agreement is poor for higher-order moments such as the skewness and flatness near the walls of the channel. Detailed information given about the effects of spatial grid resolution on a computed results is important for estimating the size of the computation required to study various aspects of a turbulent flow.     

The completely positive and doubly nonnegative completion problems

We prove the following. Let G be an undirected graph. Every partially specified symmetric matrix, the graph of whose specified entries is G and each of whose fully specified submatrices is completely positive (equal to BB^T for some entrywise nonnegative matrix B), may be completed to a completely positive matrix if and only if G is a block-clique graph (a chordal graph in which distinct maximal cliques overlap in at most one vertex). The same result holds for matrices that are doubly nonnegative (entrywise nonnegative and positive semidefinite).     

Fiber formation via solution spinning of the cellulose/ammonia/ammonium thiocyanate system

Fiber formation via the cellulose/ammonia/ammonium thiocyanate system by wet spinning has been investigated. This report presents a characterization of the structure and tensile properties of fibers spun under various coagulation conditions. Microscopic observations showed that the molecular size of coagulant was the dominant factor governing the crosssectional shape of the fibers. Density, birefringence, and crystallinity data indicated that a higher cellulose concentration and lower coagulation temperature favored development of a fiber with a denser and more oriented structure. Under optimum conditions, a welldefined fibrillar structure was obtained. Fiber tensile property measurements suggested the existence of a linear relationship between the fiber breaking tenacity and the product of the square of the Hermans' orientation factor and the infrared crystallinity index.