Climatological features of stratospheric streamers in the FUB-CMAM with increased horizontal resolution

International audienceThe purpose of this study is to investigate horizontal transport processes in the winter stratosphere using data with a resolution relevant for chemistry and climate modeling. For this reason the Freie Universität Berlin Climate Middle Atmosphere Model (FUB-CMAM) with its model top at 83 km altitude, increased horizontal resolution T42 and the semi-Lagrangian transport scheme for advecting passive tracers is used. A new approach of this paper is the classification of specific transport phenomena within the stratosphere into tropical-subtropical streamers (e.g. Offermann et al., 1999) and polar vortex extrusions hereafter called polar vortex streamers. To investigate the role played by these large-scale structures on the inter-annual and seasonal variability of transport processes in northern mid-latitudes, the global occurrence of such streamers was calculated based on a 10-year model climatology, concentrating on the existence of the Arctic polar vortex. For the identification and counting of streamers, the new method of zonal anomaly was chosen. The analysis of the months October-May yielded a maximum occurrence of tropical-subtropical streamers during Arctic winter and spring in the middle and upper stratosphere. Synoptic maps revealed highest intensities in the subtropics over East Asia with a secondary maximum over the Atlantic in the northern hemisphere. Furthermore, tropical-subtropical streamers exhibited a higher occurrence than polar vortex streamers, indicating that the subtropical barrier is more permeable than the polar vortex barrier (edge) in the model, which is in good correspondence with observations (e.g. Plumb, 2002; Neu et al., 2003). Interesting for the total ozone decrease in mid-latitudes is the consideration of the lower stratosphere for tropical-subtropical streamers and the stratosphere above ~20 km altitude for polar vortex streamers, where strongest ozone depletion is observed at polar latitudes (WMO, 2003). In the lower stratosphere the FUB-CMAM simulated a climatological maximum of 10% occurrence of tropical-subtropical streamers over East-Asia/West Pacific and the Atlantic during early- and mid-winter. The results of this paper demonstrate that stratospheric streamers e.g. large-scale, tongue-like structures transporting tropical-subtropical and polar vortex air masses into mid-latitudes occur frequently during Arctic winter. They can therefore play a significant role on the strength and variability of the observed total ozone decrease at mid-latitudes and should not be neglected in future climate change studies

Colloid-colloid and colloid-wall interactions in driven suspensions

We investigate the non-equilibrium fluid structure mediated forces between two colloids driven through a suspension of mutually non-interacting Brownian particles as well as between a colloid and a wall in stationary situations. We solve the Smoluchowski equation in bispherical coordinates as well as with a method of reflections, both in linear approximation for small velocities and numerically for intermediate velocities, and we compare the results to a superposition approximation considered previously. In particular we find an enhancement of the friction (compared to the friction on an isolated particle) for two colloids driven side by side as well as for a colloid traveling along a wall. The friction on tailgating colloids is reduced. Colloids traveling side by side experience a solute induced repulsion while tailgating colloids are attracted to each other.Comment: 8 Pages, 8 figure

Non-universal ordering of spin and charge in stripe phases

We study the interplay of topological excitations in stripe phases: charge dislocations, charge loops, and spin vortices. In two dimensions these defects interact logarithmically on large distances. Using a renormalization-group analysis in the Coulomb gas representation of these defects, we calculate the phase diagram and the critical properties of the transitions. Depending on the interaction parameters, spin and charge order can disappear at a single transition or in a sequence of two transitions (spin-charge separation). These transitions are non-universal with continuously varying critical exponents. We also determine the nature of the points where three phases coexist.Comment: 4 pages, 3 figure

From Equilibrium to Steady-State Dynamics after Switch-On of Shear

A relation between equilibrium, steady-state, and waiting-time dependent dynamical two-time correlation functions in dense glass-forming liquids subject to homogeneous steady shear flow is discussed. The systems under study show pronounced shear thinning, i.e., a significant speedup in their steady-state slow relaxation as compared to equilibrium. An approximate relation that recovers the exact limit for small waiting times is derived following the integration through transients (ITT) approach for the nonequilibrium Smoluchowski dynamics, and is exemplified within a schematic model in the framework of the mode-coupling theory of the glass transition (MCT). Computer simulation results for the tagged-particle density correlation functions corresponding to wave vectors in the shear-gradient directions from both event-driven stochastic dynamics of a two-dimensional hard-disk system and from previously published Newtonian-dynamics simulations of a three-dimensional soft-sphere mixture are analyzed and compared with the predictions of the ITT-based approximation. Good qualitative and semi-quantitative agreement is found. Furthermore, for short waiting times, the theoretical description of the waiting time dependence shows excellent quantitative agreement to the simulations. This confirms the accuracy of the central approximation used earlier to derive fluctuation dissipation ratios (Phys. Rev. Lett. 102, 135701). For intermediate waiting times, the correlation functions decay faster at long times than the stationary ones. This behavior is predicted by our theory and observed in simulations.Comment: 16 pages, 12 figures, submitted to Phys Rev

A dynamic density functional theory for particles in a flowing solvent

We present a dynamic density functional theory (dDFT) which takes into accou nt the advection of the particles by a flowing solvent. For potential flows we can use the same closure as in the absence of solvent flow. The structure of the resulting advected dDFT suggests that it could be used for non-potential flows as well. We apply this dDFT to Brownian particles (e.g., polymer coils) in a solvent flowing around a spherical obstacle (e.g., a colloid) and compare the results with direct simulations of the underlying Brownian dynamics. Although numerical limitations do not allow for an accurate quantitative check of the advected dDFT both show the same qualitative features. In contrast to previous works which neglected the deformation of the flow by the obstacle, we find that the bow-wave in the density distribution of particles in front of the obstacle as well as the wake behind it are reduced dramatically. As a consequence the friction force exerted by the (polymer) particles on the colloid can be reduced drastically.Comment: 7 pages, 5 figures, 2 tables, submitte

Colloids dragged through a polymer solution: experiment, theory and simulation

We present micro-rheological measurments of the drag force on colloids pulled through a solution of lambda-DNA (used here as a monodisperse model polymer) with an optical tweezer. The experiments show a violation of the Stokes-Einstein relation based on the independently measured viscosity of the DNA solution: the drag force is larger than expected. We attribute this to the accumulation of DNA infront of the colloid and the reduced DNA density behind the colloid. This hypothesis is corroborated by a simple drift-diffusion model for the DNA molecules, which reproduces the experimental data surprisingly well, as well as by corresponding Brownian dynamics simulations.Comment: 9 pages, 13 figures, 3 table

Three-dimensional central-moments-based lattice Boltzmann method with external forcing: A consistent, concise and universal formulation

The cascaded or central-moments-based lattice Boltzmann method (CM-LBM) is a robust alternative to the more conventional BGK-LBM for the simulation of high-Reynolds number flows. Unfortunately, its original formulation makes its extension to a broader range of physics quite difficult. To tackle this issue, a recent work [A. De Rosis, Phys. Rev. E 95, 013310 (2017)] proposed a more generic way to derive concise and efficient three-dimensional CM-LBMs. Knowing the original model also relies on central moments that are derived in an adhoc manner, i.e., by mimicking those of the Maxwell-Boltzmann distribution to ensure their Galilean invariance a posteriori, a very recent effort [A. De Rosis and K. H. Luo, Phys. Rev. E 99, 013301 (2019)] was proposed to further generalize their derivation. The latter has shown that one could derive Galilean invariant CMs in a systematic and a priori manner by taking into account high-order Hermite polynomials in the derivation of the discrete equilibrium state. Combining these two approaches, a compact and mathematically sound formulation of the CM-LBM with external forcing is proposed. More specifically, the proposed formalism fully takes advantage of the D3Q27 discretization by relying on the corresponding set of 27 Hermite polynomials (up to the sixth order) for the derivation of both the discrete equilibrium state and the forcing term. The present methodology is more consistent than previous approaches, as it properly explains how to derive Galilean invariant CMs of the forcing term in an a priori manner. Furthermore, while keeping the numerical properties of the original CM-LBM, the present work leads to a compact and simple algorithm, representing a universal methodology based on CMs and external forcing within the lattice Boltzmann framework.Comment: Published in Phys. Fluids as Editor's Pic

Residual Stresses in Glasses

The history dependence of the glasses formed from flow-melted steady states by a sudden cessation of the shear rate $\dot\gamma$ is studied in colloidal suspensions, by molecular dynamics simulations, and mode-coupling theory. In an ideal glass, stresses relax only partially, leaving behind a finite persistent residual stress. For intermediate times, relaxation curves scale as a function of $\dot\gamma t$, even though no flow is present. The macroscopic stress evolution is connected to a length scale of residual liquefaction displayed by microscopic mean-squared displacements. The theory describes this history dependence of glasses sharing the same thermodynamic state variables, but differing static properties.Comment: submitted to Physical Revie