Smart Geodesic Tracing in GRworkbench

We have developed a new tool for numerical work in General Relativity: GRworkbench. We discuss how GRworkbench's implementation of a numerically-amenable analogue to Differential Geometry facilitates the development of robust and chart-independent numerical algorithms. We consider, as an example, geodesic tracing on two charts covering the exterior Schwarzschild space-time.     

A primer on perfect simulation for spatial point processes

This primer provides a self-contained exposition of the case where spatial birth-and-death processes are used for perfect simulation of locally stable point processes. Particularly, a simple dominating coupling from the past (CFTP) algorithm and the CFTP algorithms introduced in [13], [14], and [5] are studied. Some empirical results for the algorithms are discussed. Received: 30 June 2002     

Stability of Mappings with Bounded Distortion on the Heisenberg Group

A stability theorem is proved for mappings with bounded distortion over John domains in the Heisenberg group furnished with the Carnot–Caratheodory metric.     

Temperature Stresses in a Piecewise-Homogeneous Cylindrical Shell Caused by the Presence of a Longitudinal Crack

By using the method of generalized conjugation problems, we propose a numerical scheme for investigation of the redistribution of temperature stresses in a piecewise-homogeneous cylindrical shell caused by a longitudinal crack. This scheme is based on systems of integral equations (some of them are singular) to determine the unknown jumps of integral characteristics of the disturbed temperature field and displacements and their derivatives on the line of a crack and on the interface as well as the derivatives of these integral characteristics with respect to the longitudinal coordinate at the interface.     

Barotropic Instability of the Bickley Jet at High Reynolds Numbers

We investigate the linear stability of the Bickley jet in the framework of the beta-plane approximation. Because singular inviscid neutral modes exist in the retrograde case     , it is necessary to add viscosity to interpret them. One of these modes was found in closed form by Howard and Drazin [1] . However, its critical point is at the center of the jet and it was therefore not possible for these authors to ascertain the relationship of this mode to the stability problem or to discuss how to continue the eigenfunction across the singularity.
The viscous critical layer problem associated with this singularity is considerably more difficult than the usual one (which leads to integrals of the Airy function) because     and, consequently, a second-order turning point is involved. Our analysis shows that the Howard–Drazin mode is degenerate in the domain where it is valid as a limit of the viscous problem (wavenumber  α²≤ 9/2  ), that is, it corresponds to both an odd and an even mode. This conclusion is confirmed by direct numerical solution of the Orr–Sommerfeld equation which shows, in addition, that viscosity is destabilizing along portions of the stability boundary. For a retrograde jet, instability is found to occur beyond the inviscid critical value of β, that is, in the region where the flow would be stable according to the Rayleigh–Kuo condition.     

Interference of spin splittings in magneto-oscillation phenomena in two-dimensional systems

The spin splitting caused by the terms linear in wavevector in the effective Hamiltonian containing can give rise to the new magneto-oscillation phenomena in two-dimensional systems. It is shown that the joint action of the spin-dependent contributions due to the heterostructure asymmetry and to the lack of inversion center in the bulk material suppresses beats that arise in the magneto-oscillation phenomena in the presence of the terms of only one of these types.     

Field Theory Methods in Classical Dynamics

A Dirac picture perturbation theory is developed for the time evolution operator in classical dynamics in the spirit of the Schwinger–Feynman–Dyson perturbation expansion and detailed rules are derived for computations. Complexification formalisms are given for the time evolution operator suitable for phase space analyses, and then extended to a two-dimensional setting for a study of the geometrical Berry phase as an example. Finally a direct integration of Hamilton's equations is shown to lead naturally to a path integral expression, as a resolution of the identity, as applied to arbitrary functions of generalized coordinates and momenta.     

Computer simulation of selective absorption of radiation by the components of a light-scattering layer

A method of calculating the partial characteristics of radiation absorption by the components of light-scattering disperse layers is proposed. This method is based on statistical modeling (the Monte Carlo method). The absorptivities of photographic gelatin and silver bromide microcrystals and the corresponding distributions of the absorbed energy over the layer thickness are calculated using the example of an interaction between actinic radiation and silver halide photographic layers in the wavelength range λ=200–440 nm. The following structural parameters of the photographic layer are used in the calculation: the mean size of emulsion crystals d=0.5 μm; the polydispersity C _V =25%; the volume concentrations C _V =10, 20, and 30%; and the thickness of the emulsion layer H=10 μm.