Identification of nonstationary axisymmetric load distributed along a cylindrical shell

The paper outlines a procedure to identify the space-and time-dependent external nonstationary load acting on a closed circular cylindrical shell of medium thickness. Time-dependent deflections at several points of the shell are used as input data to solve the inverse problem. Examples of numerical identification of various nonstationary loads, including moving ones are presented. The relationship between the external load and the stress-strain state of the shell is described by the Volterra equation of the first kind. The identification problem is solved using Tikhonov's regularization method and Apartsin's h-regularization method __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 7, pp. 91–100, July 2008.     

A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

A posteriori H1 error estimation for a singularly perturbed reaction diffusion problem on anisotropic meshes

The paper deals with a singularly perturbed reaction diffusionmodel problem. The focus is on reliable a posteriori error estimatorsfor the H¹ seminorm that can be applied to anisotropic finiteelement meshes. A residual error estimator and a local problemerror estimator are proposed and rigorously analysed. They arelocally equivalent, and both bound the error reliably. Threemodifications of these estimators are introduced and discussed. Much attention is given to the performance of the error estimatorin numerical experiments. This helps to identify those estimatorsthat are suitable for practical applications.     

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     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

On a problem of Zygmund

It is proved that there exists an integrable function on [0, 1]² whose integral is nondifferentiable in each direction belonging to a set everywhere dense in [0, 2π] but is strongly differentiable. Translated fromMatematicheskie Zametki, Vol. 64, No. 5, pp. 749–762, November, 1998.     

Minimum energy configurations on a toric lattice as a quadratic assignment problem

We consider three known bounds for the quadratic assignment problem (QAP): an eigenvalue, a convex quadratic programming (CQP), and a semidefinite programming (SDP) bound. Since the last two bounds were not compared directly before, we prove that the SDP bound is stronger than the CQP bound. We then apply these to improve known bounds on a discrete energy minimization problem, reformulated as a QAP, which aims to minimize the potential energy between repulsive particles on a toric grid. Thus we are able to prove optimality for several configurations of particles and grid sizes, complementing earlier results by Bouman et al. (2013). The semidefinite programs in question are too large to solve without pre-processing, and we use a symmetry reduction method by Permenter and Parrilo (2020) to make computation of the SDP bounds possible.     

Accurate electrical and spectroscopic properties ofX 1Σ+ BeO from coupled-cluster methods

Summary This paper reports a series of coupled-cluster (CC) calculations through CCSDT on the theoretically challenging ground state of the BeO molecule. Along with CC methods, quadratic configuration interaction (QCI) approximations to CC theory have been used (QCISD and QCISD(T)), which show several dramatic failings. Equilibrium electrical properties (, _xx , and _zz ) and basic spectroscopic properties (r _e, _e,D _e, and infrared intensity (I)) have been computed. Basis set and electron correlation effects are analyzed in order to arrive at accurate values of the dipole moment and polarizability, which are not known experimentally. For the dipole moment, we obtain a value of 6.25 D, with an uncertainty of about 0.1 D. For _xx and _zz , we suggest respective values of 32 and 36 atomic units (a.u.) and error bars of about 1 and 2 a.u. With extended basis sets, the spectroscopic propertiesr _e, _e, andD _e are reproduced to high accuracy, which is the first time this has been achieved for this species byab initio methods. At the highest calculation levels,I is predicted to be very small. AlthoughI has not been measured, some support for this prediction comes from a recent infrared study of BeO-rare gas complexes. The QCI methods are shown to be much more sensitive to basis set, and even with large basis sets yield values of _zz andI which differ from CC results by an order of magnitude and three orders of magnitude, respectively. These differences doubtless arise from the importance of single excitations (T ₁) for this molecule, as several terms involvingT ₁ are neglected in the QCISD approximation compared with CCSD. We also report CC calculations with Brueckner orbitals, which yield results similar to those obtained with restricted Hartree-Fock orbitals.     

Necessary and sufficient conditions for the existence of positive solutions to algebraic Riccati equations with indefinite quadratic term

Necessary and sufficient conditions are established in this paper for the existence of positive- and/or negative-definite solutions to the algebraic Riccati equation with indefinite coefficient. An iterative procedure is also given for computing such a solution.Project supported by the National Science Foundation of China and by the special program of the State Education Commission of China under grant 9033507.     

On a class of least-element complementarity problems

The present paper studies the linear complementarity problem of finding vectorsx andy inR ₊ ⁿ such thatc + Dx + y 0,b – x 0 andx ^T (c + Dx + y) = y ^T (b – x) = 0 whereD is aZ-matrix andb > 0. Complementarity problems of this nature arise, for example, from the minimization of certain quadratic functions subject to upper and lower bounds on the variables. Two least-element characterizations of solutions to the above linear complementarity problem are established first. Next, a new and direct method to solve this class of problems, which depends on the idea of least-element solution is presented. Finally, applications and computational experience with its implementation are discussed.Research partially supported by the National Science Foundation Grant MCS 71-03341 A04 and the Air Force Office of Scientific Research Contract F 44620 14 C 0079.