The numerical approximation of center manifolds in Hamiltonian systems

In this paper we develop a numerical method for computing higher order local approximations of center manifolds near steady states in Hamiltonian systems. The underlying system is assumed to be large in the sense that a large sparse Jacobian at the equilibrium occurs, for which only a linear solver and a low-dimensional invariant subspace is available. Our method combines this restriction from linear algebra with the requirement that the center manifold is parametrized by a symplectic mapping and that the reduced equation preserves the Hamiltonian form. Our approach can be considered as a special adaptation of a general method from Numer. Math. 80 (1998) 1-38 to the Hamiltonian case such that approximations of the reduced Hamiltonian are obtained simultaneously. As an application we treat a finite difference system for an elliptic problem on an infinite strip.     

Analysis and modelling of overtone singing in the sygyt style

Overtone singing, also called biphonic singing, xöömij or chant diphonique in french; is a special singing style that exhibits two or more separate sounds - one “drone” sound of relatively low pitch and one or more high pitch melody sounds. The perceived pitches of the upper tones are multiples of the drone sound, i.e. taken from its overtone scale. Compared to voiced sounds of western style singers, the relative amplitude of the melody pitches is quite high, and the formant bandwidth of overtone sounds is small. This paper tries to answer the question of how these formant properties are achieved. Experimental investigations and numerical calculations prove the existence of two closely neighboured formants for the production of the melody sound in the sygyt style.     

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°.     

半导体与超晶格的数学模型及其分析Ⅰ:模型

本文论及半导体与超晶格的数学模型方法的一些新近发展动态。介绍了半导体与超晶格的一些背景材料并给出了半导体宏观模型间的层次框架。超晶格中载流子传输的SHE扩散模型也被观察。     

Plastic multipliers as driving variables of numerical simulation in elastoplasticity

Statics of elasto-plastic media is stated in terms of eigenstrains acting upon a background elastic problem with fixed (initial) stiffness. In order to minimize the number of unknowns and to provide computationally cheap algorithms plastic multipliers are used as main driving variables. Fixed-point type iterations are suggested for computing plastic multipliers within the load increment. Numerical experiments results are shown to be in good correlation with other software results. A range of problems is outlined where the usage of the proposed algorithms can be advantageous.     

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.     

Development of a finite-element analysis methodology for the propagation of delaminations in composite structures

Analysing the collapse of skin-stiffened structures requires capturing the critical phenomenon of skin-stiffener separation, which can be considered analogous to interlaminar cracking. This paper presents the development of a numerical approach for simulating the propagation of interlaminar cracks in composite structures. A degradation methodology was introduced in MSC.Marc, which involved the modelling of a structure with shell layers connected by user-defined multiple-point constraints (MPCs). User subroutines were written that employ the virtual crack closure technique (VCCT) to determine the onset of crack growth and modify the properties of the user-defined MPCs to simulate crack propagation. Methodologies for the release of failing MPCs are presented and are discussed with reference to the VCCT assumption of self-similar crack growth. The numerical results obtained by using the release methodologies are then compared with experimental data for a double-cantilever beam specimen. Based on this comparison, recommendations for the future development of the degradation model are made, especially with reference to developing an approach for the collapse analysis of fuselage-representative structures. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 1, pp. 15–42, January–February, 2007.     

Standing waves, clustering, and phase waves in 1D simulations of kinetic relaxation oscillations in NO+NH3 on Pt(1 0 0) coupled by diffusion

The Lombardo–Imbihl–Fink (LFI) ODE model of the NO+NH₃ reaction on a Pt(1 0 0) surface shows stable relaxation oscillations with very sharp transitions for temperatures T between 404 and 433 K. Here we study numerically the effect of linear diffusive coupling of these oscillators in one spatial dimension. Depending on the parameters and initial conditions we find a rich variety of spatio-temporal patterns which we group into four main regimes: bulk oscillations (BOs), standing waves (SW), phase clusters (PC), and phase waves (PW). Two key ingredients for SW and PC are identified, namely the relaxation type of the ODE oscillations and a nonlocal (and nonglobal) coupling due to relatively fast diffusion of the kinetically slaved variables NH₃ and H. In particular, the latter replaces the global coupling through the gas phase used to obtain SW and PC in models of related surface reactions. The PW exist only under the assumption of (relatively) slow diffusion of NH₃ and H.     

Numerical analysis of roughness effect on microtube heat transfer

Roughness effect on the heat transfer and pressure loss performances of microscale tubes and channels are investigated through a finite element CFD code. Surface roughness is explicitly modelled through a set of random generated peaks along the ideal smooth surface. Different peak shapes and distributions are considered; geometrical parameters are representative of tubes in the diameter range from 50 to 150 μm. The use of a fine enough mesh allows the direct computation of tube performances under the assumption of incompressible fully developed flow. As a result, a significant increase in Poiseuille number is detected for all of the configurations considered, while the effect of roughness on heat transfer rate is smaller and highly dependent on the tube shape.