Nonlinear stability analysis of plane Poiseuille flow by normal forms

In the subcritical interval of the Reynolds number 4320\leq R\leq R_c\equiv 5772, the Navier--Stokes equations of the two--dimensional plane Poiseuille flow are approximated by a 22--dimensional Galerkin representation formed from eigenfunctions of the Orr--Sommerfeld equation. The resulting dynamical system is brought into a generalized normal form which is characterized by a disposable parameter controlling the magnitude of denominators of the normal form transformation. As rigorously proved, the generalized normal form decouples into a low--dimensional dominant and a slaved subsystem. {}From the dominant system the critical amplitude is calculated as a function of the Reynolds number. As compared with the Landau method, which works down to R=5300, the phase velocity of the critical mode agrees within 1 per cent; the critical amplitude is reproduced similarly well except close to the critical point, where the maximal error is about 16 per cent. We also examine boundary conditions which partly differ from the usual ones.Comment: latex file; 4 Figures will be sent, on request, by airmail or by fax (e-mail address: rauh at beta.physik.uni-oldenburg.de

A risk profile for information fusion algorithms

E.T. Jaynes, originator of the maximum entropy interpretation of statistical mechanics, emphasized that there is an inevitable trade-off between the conflicting requirements of robustness and accuracy for any inferencing algorithm. This is because robustness requires discarding of information in order to reduce the sensitivity to outliers. The principal of nonlinear statistical coupling, which is an interpretation of the Tsallis entropy generalization, can be used to quantify this trade-off. The coupled-surprisal, -ln_k (p)=-(p^k-1)/k, is a generalization of Shannon surprisal or the logarithmic scoring rule, given a forecast p of a true event by an inferencing algorithm. The coupling parameter k=1-q, where q is the Tsallis entropy index, is the degree of nonlinear coupling between statistical states. Positive (negative) values of nonlinear coupling decrease (increase) the surprisal information metric and thereby biases the risk in favor of decisive (robust) algorithms relative to the Shannon surprisal (k=0). We show that translating the average coupled-surprisal to an effective probability is equivalent to using the generalized mean of the true event probabilities as a scoring rule. The metric is used to assess the robustness, accuracy, and decisiveness of a fusion algorithm. We use a two-parameter fusion algorithm to combine input probabilities from N sources. The generalized mean parameter 'alpha' varies the degree of smoothing and raising to a power N^beta with beta between 0 and 1 provides a model of correlation.Comment: 15 pages, 4 figure

Radiative damping: a case study

We are interested in the motion of a classical charge coupled to the Maxwell self-field and subject to a uniform external magnetic field, B. This is a physically relevant, but difficult dynamical problem, to which contributions range over more than one hundred years. Specifically, we will study the Sommerfeld-Page approximation which assumes an extended charge distribution at small velocities. The memory equation is then linear and many details become available. We discuss how the friction equation arises in the limit of "small" B and contrast this result with the standard Taylor expansion resulting in a second order equation for the velocity of the charge.Comment: 4 figure

Electrodynamic Radiation Reaction and General Relativity

We argue that the well-known problem of the instabilities associated with the self-forces (radiation reaction forces) in classical electrodynamics are possibly stabilized by the introduction of gravitational forces via general relativity

Susceptibility divergence, phase transition and multistability of a highly turbulent closed flow

Using time-series of stereoscopic particle image velocimetry data, we study the response of a turbulent von K\'{a}rm\'{a}n swirling flow to a continuous breaking of its forcing symmetry. Experiments are carried over a wide Reynolds number range, from laminar regime at $Re = 10^{2}$ to highly turbulent regime near $Re = 10^{6}$. We show that the flow symmetry can be quantitatively characterized by two scalars, the global angular momentum $I$ and the mixing layer altitude $z_s$, which are shown to be statistically equivalent. Furthermore, we report that the flow response to small forcing dissymetry is linear, with a slope depending on the Reynolds number: this response coefficient increases non monotonically from small to large Reynolds number and presents a divergence at a critical Reynolds number $Re_c = 40\,000 \pm 5\,000$. This divergence coincides with a change in the statistical properties of the instantaneous flow symmetry $I(t)$: its pdf changes from Gaussian to non-Gaussian with multiple maxima, revealing metastable non-symmetrical states. For symmetric forcing, a peak of fluctuations of $I(t)$ is also observed at $Re_c$: these fluctuations correspond to time-intermittencies between metastable states of the flow which, contrary to the very-long-time-averaged mean flow, spontaneously and dynamically break the system symmetry. We show that these observations can be interpreted in terms of divergence of the susceptibility to symmetry breaking, revealing the existence of a phase transition. An analogy with the ferromagnetic-paramagnetic transition in solid-state physics is presented and discussed.Comment: to appear in Journal of Statistical Mechanic

A Low Mach Number Solver: Enhancing Applicability

In astrophysics and meteorology there exist numerous situations where flows exhibit small velocities compared to the sound speed. To overcome the stringent timestep restrictions posed by the predominantly used explicit methods for integration in time the Euler (or Navier-Stokes) equations are usually replaced by modified versions. In astrophysics this is nearly exclusively the anelastic approximation. Kwatra et al. have proposed a method with favourable time-step properties integrating the original equations (and thus allowing, for example, also the treatment of shocks). We describe the extension of the method to the Navier-Stokes and two-component equations. - However, when applying the extended method to problems in convection and double diffusive convection (semiconvection) we ran into numerical difficulties. We describe our procedure for stabilizing the method. We also investigate the behaviour of Kwatra et al.'s method for very low Mach numbers (down to Ma = 0.001) and point out its very favourable properties in this realm for situations where the explicit counterpart of this method returns absolutely unusable results. Furthermore, we show that the method strongly scales over 3 orders of magnitude of processor cores and is limited only by the specific network structure of the high performance computer we use.Comment: author's accepted version at Elsevier,JCP; 42 pages, 14 figure

Wheeler-DeWitt Equation in 2 + 1 Dimensions

The infrared structure of quantum gravity is explored by solving a lattice version of the Wheeler-DeWitt equations. In the present paper only the case of 2+1 dimensions is considered. The nature of the wavefunction solutions is such that a finite correlation length emerges and naturally cuts off any infrared divergences. Properties of the lattice vacuum are consistent with the existence of an ultraviolet fixed point in $G$ located at the origin, thus precluding the existence of a weak coupling perturbative phase. The correlation length exponent is determined exactly and found to be $\nu=6/11$. The results obtained so far lend support to the claim that the Lorentzian and Euclidean formulations belong to the same field-theoretic universality class.Comment: 56 pages, 7 figures, typos fixed, references adde

Ab-Initio Molecular Dynamics

Computer simulation methods, such as Monte Carlo or Molecular Dynamics, are very powerful computational techniques that provide detailed and essentially exact information on classical many-body problems. With the advent of ab-initio molecular dynamics, where the forces are computed on-the-fly by accurate electronic structure calculations, the scope of either method has been greatly extended. This new approach, which unifies Newton's and Schr\"odinger's equations, allows for complex simulations without relying on any adjustable parameter. This review is intended to outline the basic principles as well as a survey of the field. Beginning with the derivation of Born-Oppenheimer molecular dynamics, the Car-Parrinello method and the recently devised efficient and accurate Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics, which unifies best of both schemes are discussed. The predictive power of this novel second-generation Car-Parrinello approach is demonstrated by a series of applications ranging from liquid metals, to semiconductors and water. This development allows for ab-initio molecular dynamics simulations on much larger length and time scales than previously thought feasible.Comment: 13 pages, 3 figure

Coherent Compton scattering on light nuclei in the delta resonance region

Coherent Compton scattering on light nuclei in the delta resonance region is studied in the impulse approximation and is shown to be a sensitive probe of the in-medium properties of the delta resonance. The elementary amplitude on a single nucleon is calculated from the unitary K-matrix approach developed previously. Modifications of the properties of the delta resonance due to the nuclear medium are accounted for through the self-energy operator of the delta, calculated from the one-pion loop. The dominant medium effects such as the Pauli blocking, mean-field modification of the nucleon and delta masses, and particle-hole excitations in the pion propagator are consistently included in nuclear matter.Comment: 30 pages, 11 figures, accepted for publication in Phys. Rev.