Well-posed absorbing layer for hyperbolic problems

Summary. The perfectly matched layer (PML) is an efficient tool to simulate propagation phenomena in free space on unbounded domain. In this paper we consider a new type of absorbing layer for Maxwell's equations and the linearized Euler equations which is also valid for several classes of first order hyperbolic systems. The definition of this layer appears as a slight modification of the PML technique. We show that the associated Cauchy problem is well-posed in suitable spaces. This theory is finally illustrated by some numerical results. It must be underlined that the discretization of this layer leads to a new discretization of the classical PML formulation. Received May 5, 2000 / Published online November 15, 2001     

On the Point-Particle (Newtonian) Limit¶of the Non-Linear Hartree Equation

We consider the nonlinear Hartree equation describing the dynamics of weakly interacting non-relativistic Bosons. We show that a nonlinear M?ller wave operator describing the scattering of a soliton and a wave can be defined. We also consider the dynamics of a soliton in a slowly varying background potential W(ɛx). We prove that the soliton decomposes into a soliton plus a scattering wave (radiation) up to times of order ɛ⁻¹. To leading order, the center of the soliton follows the trajectory of a classical particle in the potential W(ɛx). Received: 30 June 2000 / Accepted: 25 June 2001     

A posteriori error estimators for the Stokes equations

Summary We present two a posteriori error estimators for the mini-element discretization of the Stokes equations. One is based on a suitable evaluation of the residual of the finite element solution. The other one is based on the solution of suitable local Stokes problems involving the residual of the finite element solution. Both estimators are globally upper and locally lower bounds for the error of the finite element discretization. Numerical examples show their efficiency both in estimating the error and in controlling an automatic, self-adaptive mesh-refinement process. The methods presented here can easily be generalized to the Navier-Stokes equations and to other discretization schemes.This work was accomplished at the Universität Heidelberg with the support of the Deutsche Forschungsgemeinschaft     

A simple proof of Sanov’s theorem*

A simple self-contained proof of Sanov's theorem in τ-topology is given, well suited for a first course on large deviations. *This work was supported by the Hungarian National Foundation for Scientific Research under Grant T046376     

SPM-induced spectral compression of picosecond pulses in a single-mode Yb-doped fiber amplifier

In a fiber amplifier, spectral compression due to self-phase modulation is demonstrated for ultrashort pulses. We report the generation of near-transform-limited picosecond pulses with peak powers of several kW at a repetition rate of 74 MHz and diffraction-limited beam quality in a Yb-doped fiber amplifier when seeding with a negative chirped pulse. Received: 17 September 2001 / Revised version: 22 November 2001 / Published online: 17 January 2002     

Mixing mechanisms of vortex ring formed by gravity slumping motion

 This paper is concerned with an experimental investigation of the mixing inside the vortex ring formed by the gravity slumping motion of a dense cloud in a less dense atmosphere. The dynamics of the spreading and instantaneous structures of the turbulent flow were examined by visualization, single and multi-point measurements of velocity and concentration for two heavy gases, carbondioxide (CO₂) and dichlorodifluoromethane (CCl₂F₂), in a configuration in that heavy gas, initially trapped in a reservoir, was released with the rise of a shutter into calm air of a sector-shaped dispersion channel. Visualization of the cloud as a whole showed a spreading motion in which an advancing frontal structure was followed by a stratified flow with a layer of dense fluid of higher velocities near the wall and, on top of it, a layer of dilute fluid whose concentration is controlled by the mixing mechanisms within the head. During the course of spreading, there was always a phase in which the head attained to a constant speed of advance, which occurred as 0.13 m/s for CO₂ and 0.48 m/s for CCl₂F₂. It was interesting to observe for CO₂ that the phase of constant speed took place in between two acceleration phases; the former was due to the initial slumping of the cloud at the exit of the reservoir, and the latter was attributed to the collapse of the head on the transition to the passive dispersion phase. Instantaneous two-dimensional velocity field, measured with particle image velocimeter (PIV), showed that the cloud overran the ambient air which caused the approaching dense fluid deflected away from the wall with significant vertical velocities and downstream-moving separation, and the air trapped under the head resulted in the density inversion which introduced further intricacy to the turbulent structure of the head. Instabilities at the upper free shear layer due to density and velocity discontinuity rolled into periodic array of vortices which engulfed a considerable amount of air as they were convected backwards over the head, but the incorporation of heavy and light fluids was completed with the appearance of microscales after the collapse on the stratified layer. Analyses of the cloud head at different downstream locations also revealed that its size remained unchanged when the speed of advance was constant, allowing the rate of change of the cloud volume being modeled with the rate of spreading. Contours of concentration obtained from digitized PIV pictures confirmed the kinematic features of the mixing revealed by the velocity field and that the concentration values within the large structures were higher than those at the upper part of the stratified layer. Motivated by the experimental observations, a semi-empirical analysis was presented to describe the results and based on local values of the Richardson and Reynolds numbers. Received: 4 October 1995 / Accepted: 4 July 1996     

Fibre stability in anisotropic matrix at small strains

As acknowledged in almost all monographs on the fracture of composite materials, one of the major fracture mechanism in unidirectional fibrous composites under uniaxial compression along the reinforcing elements is the stability loss of the material structure (the structural instability). According to this mechanism, theoretical investigations of the fracture along the fibres are reduced to those of the stability loss in the material structure, and the value of external critical forces is accepted as the value of failure forces. At present, numerous theoretical investigations have been carried out in this field with the use of the three-dimensional linearized theory of stability in the framework of the piecewise-homogeneous body model. However, in all the investigations it is assumed that the matrix and the fibre material are isotropic. It is evident that in many cases it is necessary to take into account the anisotropy of the matrix material when investigating the stability loss of fibres. In view of the above, in the framework of the piecewise-homogeneous body model using the three-dimensional linearized theory of stability, the present paper considers the stability loss of the fibre in the anisotropic (transversally isotropic) matrix. The effect of the properties of the matrix material on the critical values of the external loading is examined.Submitted to the 10th International Conference of Mechanics of Composite Materials (Riga, April 20–23, 1998).Published in Mekhanika Kompozitnykh Materialov, Vol. 33, No. 5, pp. 603–611. September–October, 1997.     

Acousto-optic mode coupling of partially coherent light in two-mode optical fibers

An acoustic pulse propagating on a two-mode fiber can act as a beam splitter in a scanning interferometer. When this device is employed in white-light interferometry, the effects of distributed coupling and dispersive interferometer arms need to be considered. A theory suitable for treating acousto-optic interaction of partially coherent light in a moving interaction region was developed. It was found that differential optical dispersion should be negligible and the acoustic pulse length short. Also the coherence time should be short but long compared to the intermodal group delay difference over a pulse length. Experiments with long acoustic pulses were performed, and fairly good agreement with theory was obtained.     

Diffusion-limited annihilation in inhomogeneous environments

We study diffusion-limited (on-site) pair annihilation A + A → 0 and (on-site) fusion A + A → A which we show to be equivalent for arbitrary space-dependent diffusion and reaction rates. For one-dimensional lattices with nearest neighbour hopping we find that in the limit of infinite reaction rate the time-dependent n-point density correlations for many-particle initial states are determined by the correlation functions of a dual diffusion-limited annihilation process with at most 2n particles initially. Furthermore, by reformulating general properties of annihilating random walks in one dimension in terms of fermionic anticommutation relations we derive an exact representation for these correlation functions in terms of conditional probabilities for a single particle performing a random walk with dual hopping rates. This allows for the exact and explicit calculation of a wide range of universal and non-universal types of behaviour for the decay of the density and density correlations.