Artificial time integration

Many recent algorithmic approaches involve the construction of a differential equation model for computational purposes, typically by introducing an artificial time variable. The actual computational model involves a discretization of the now time-dependent differential system, usually employing forward Euler. The resulting dynamics of such an algorithm is then a discrete dynamics, and it is expected to be “close enough” to the dynamics of the continuous system (which is typically easier to analyze) provided that small – hence many – time steps, or iterations, are taken. Indeed, recent papers in inverse problems and image processing routinely report results requiring thousands of iterations to converge. This makes one wonder if and how the computational modeling process can be improved to better reflect the actual properties sought. In this article we elaborate on several problem instances that illustrate the above observations. Algorithms may often lend themselves to a dual interpretation, in terms of a simply discretized differential equation with artificial time and in terms of a simple optimization algorithm; such a dual interpretation can be advantageous. We show how a broader computational modeling approach may possibly lead to algorithms with improved efficiency. AMS subject classification (2000)  65L05, 65M32, 65N21, 65N22, 65D18     

Estimation of the convergence rate in a transfer theorem for centered random sums

Estimates of the convergence rate in a transfer theorem for random sums centered by constants are presented. The estimates are obtained as a consequence of the estimates of the stability of the distributions of the sums represented in the form of special mixtures to small variations in the mixing and miscible distributions.     

Existence Results for Semilinear Neutral Functional Differential Inclusions via Analytic Semigroups

In this paper we prove existence results for semilinear neutral functional differential inclusions with finite or infinite delay in Banach spaces. Our theory makes use of analytic semigroups and fractional powers of closed operators, integrated semigroups and cosine families.      

Stored energy release behaviour of disordered carbon

The use of graphite as a moderator in a low temperature thermal nuclear reactor is restricted due to accumulation of energy caused by displacement of atoms by neutrons and high energetic particles. Thermal transients may lead to a release of stored energy that may raise the temperature of the fuel clad above the design limit. Disordered carbon is thought to be an alternative choice for this purpose. Two types of disordered carbon composites, namely, CB (made up of 15 wt. % carbon black dispersed in carbonized phenolic resin) and PAN (made up of 20 vol. % chopped polyacrylonitrile carbon fibre dispersed in carbonized phenolic resin matrix) have been irradiated with 145 MeV Ne⁶⁺ ions at three fluence levels of 1.0×10¹³, 5.0×10¹³ and 1.5×10¹⁴ Ne⁶⁺/cm², respectively. The XRD patterns revealed that both the samples remained disordered even after irradiation. The maximum release of stored energy for CB was 212 J/g and that of PAN was 906 J/g. For CB, the release of stored energy was a first order reaction with activation energy of 2.79 eV and a frequency factor of 3.72×10²⁸ per second. 13% of the defects got annealed by heating up to 700 °C. PAN showed a third-order release rate with activation energy of 1.69 eV and a frequency factor of 1.77×10¹⁴ per second. 56% of the total defects got annealed by heating it up to 700 °C. CB seems to be the better choice than PAN as it showed less energy release with a slower rate. PACS 61.80.Jh; 61.80.-x; 61.43Er; 61.43.-j; 68.43.Vx     

Dynamic programming in cricket: optimizing batting order for a sticky wicket

In cricket, a rain-affected pitch can make batting more difficult than normal. Several other conditions such as poor light or an initially lively pitch, may also result in difficulties for the batsmen. In this note, we refer to all of them as ‘sticky wickets’. On sticky wickets, lower order batsmen are often sent into ‘hold the fort’ until conditions improve. In this paper, a stochastic dynamic programming model is used to examine the appropriateness of this policy. The model suggests that the tactic is often optimal when the sticky wicket persists until the end of the day's play, but not often when the sticky wicket is transitory. In some circumstances, it is worthwhile, on a normal wicket near the end of the day, to send in a lower order batsman to hold the fort (a night watchman): when the wicket is sticky, this tactic is even more worthwhile.     

On the excitation of a closely spaced array by a line source

An infinite row of periodically spaced, identical rigid circularcylinders is excited by an acoustic line source which is parallelto the generators of the cylinders. A method for calculatingthe scattered field accurately and efficiently is presented.When the cylinders are sufficiently close together, Rayleigh–Blochsurface waves that propagate energy to infinity along the arrayare excited. An expression is derived which enables the amplitudesof these surface waves to be computed without requiring thesolution to the full scattering problem.     

A few points on omnidirectional band gaps in one-dimensional photonic crystals

The universal condition for the formation of omnidirectional band gaps (OBG) in photonic crystal (PC) was derived with consideration of permeability of the materials. And it was found that there are four kinds of PCs: one of them has no OBG, and one always possesses OBG. For the other two kinds of PCs, there are OBG for only TM or TE waves respectively. Moreover, in all PCs, the OBG can be broadened by decreasing the refractive index of the ambient medium or/and increasing the contrast between the wave impedances of the component materials of the PC. PACS 42.70.Qs; 71.20.Tx     

Sympathetic cooling rate of gas-phase ions in a radio-frequency-quadrupole ion trap

We theoretically investigated the mass dependence of the sympathetic cooling rate of gas-phase ions trapped in a linear radio-frequency-quadrupole ion trap. Using an a priori molecular dynamical calculation, tracing numerically with Newtonian equations of motion, we found that ions with a mass greater than 0.54±0.04 times that of the laser-cooled ions are sympathetically cooled; otherwise, they are heated. To understand the mass dependence obtained using the molecular-dynamical calculation, we made a heat-exchange model of sympathetic cooling, which shows that the factor of 0.54±0.04 is a consequence of absence of micro-motion along the axis of the linear ion trap. Received: 10 December 2001 / Revised version: 28 January 2002 / Published online: 14 March 2002     

Low-lying states of trinuclear mixed-valence cluster: [Fe3S4]0

Low-lying energy states of the [Fe₃S₄]⁰ cluster have been calculated by taking into account the double exchange, superexchange and vibronic interaction. It was found that the adiabatic potential of the excited state withS=0 corresponds to the full delocalization of the “excess” charge. From the analysis of experimental data of Mössbauer spectroscopy and the temperature dependence of the magnetic susceptibility the double exchange parametert≥4000 cm^?1 and the vibronic interaction parameter λ₂/2k cm^?1 have been estimated.