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     

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.     

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.     

Theory of Non-Equilibrium Stationary States as a Theory of Resonances

We study a small quantum system (e.g., a simplified model for an atom or molecule) interacting with two bosonic or fermionic reservoirs (say, photon or phonon fields). We show that the combined system has a family of stationary states parametrized by two numbers, T ₁ and T ₂ (‘reservoir temperatures’). If T ₁ ≠ T ₂, then these states are non-equilibrium stationary states (NESS). In the latter case we show that they have nonvanishing heat fluxes and positive entropy production and are dynamically asymptotically stable. The latter means that the evolution with an initial condition, normal with respect to any state where the reservoirs are in equilibria at temperatures T ₁ and T ₂, converges to the corresponding NESS. Our results are valid for the temperatures satisfying the bound min (T ₁,T ₂) > g ^{2 + α}, where g is the coupling constant and 0 < α < 1 is a power related to the infra-red behaviour of the coupling functions. Submitted: March 20, 2006. Revised: March 19, 2007. Accepted: May 11, 2007. Marco Merkli: Partly supported by an NSERC PDF, the Institute of Theoretical Physics of ETH Zürich, Switzerland, the Departments of Mathematics of McGill University and the University of Toronto, Canada. Matthias Mück: Supported by DAAD under grant HSP III. Israel Michael Sigal: Supported by NSERC under grant NA7901.     

Solubility of carbon dioxide in aqueous solutions of sodium chloride: Experimental results and correlation

The solubility of carbon dioxide in aqueous solutions of sodium chloride was measured in the temperature range from 40 to 160°C, up to 6 mol-kg salt solutions and total pressures up to 10 MPa. Pitzer's⁽¹⁾ equations as well as the Chen and Evans⁽²⁾ model were used to correlate the new data. Results are reported and compared to literature data and correlations.     

Terrestrial digital video broadcasting for mobile reception using OFDM

Performance of a system design for digital video broadcasting is examined with emphasis on mobile reception. Orthogonal frequency division multiplexing (OFDM) is used to achieve good bandwidth efficiency and to mitigate the intersymbol interference resulting from the channel delay spread. The resulting equivalent channel including OFDM can be modeled as a flat Rayleigh fading channel plus an interchannel interference (ICI) term due to the channel Doppler spread. This ICI term is analyzed and shown to result in an error floor. Performance improvements due to antenna diversity and trellis-coded modulation (TCM) are given. Finally, multiresolution modulation is discussed as a means of achieving graceful degradation and giving degrees of freedom for further performance improvement.This research was supported by the Multimedia Systems R & D Laboratory, Hitachi, Ltd.     

The excited states of pyrazine: A basis set study

Summary Multiconfigurational second order perturbation theory (CASSCF/CASPT2) has been used to investigate the dependence of computed valence excitation energies and transition moments on the basis sets. Pyrazine has been selected as the test molecule. Atomic normal orbital (ANO) type basis sets are used throughout. Contractions of the structure (4s3p1d/2s) are found to be an optimal compromise between the quality and the size of the calculations and are capable of yielding results virtually identical to more extended basis sets.     

Microwave cavity perturbation technique: Part I: Principles

This report reviews the analysis used to extract the complex conductivity of a compound from a microwave cavity perturbation measurement. We intend to present a generalized treatment valid for any spheroidally shaped sample of arbitrary conductivity which is placed at either the electric or magnetic field antinode of the cavity. To begin with, we establish the relationship between the measured parameters and the conductivity for a spherical sample. Next, we extend these results to the case of spheroids; and for the first time, we cover all different configurations that one can possibly use to study an arbitrary conducting sample inside a cavity: in particular, all possible orientations of the sample with respect to the applied field are solved.