Search

We study the turbulent evolution originated from a system subjected to a Rayleigh-Taylor instability with a double density at high resolution in a two-dimensional geometry using a highly optimized thermal lattice-Boltzmann code for GPUs. Our investigation's initial condition, given by the superposition of three layers with three different densities, leads to the development of two Rayleigh-Taylor fronts that expand upward and downward and collide in the middle of the cell. By using high-resolution numerical data we highlight the effects induced by the collision of the two turbulent fronts in the long-time asymptotic regime. We also provide details on the optimized lattice-Boltzmann code that we have run on a cluster of GPUs. [ABSTRACT FROM AUTHOR]

The parametrization of small-scale turbulent fluctuations in convective systems and in the presence of strong stratification is a key issue for many applied problems in oceanography, atmospheric science, and planetology. In the presence of stratification, one needs to cope with bulk turbulent fluctuations and with inversion regions, where temperature, density, or both develop highly nonlinear mean profiles due to the interactions between the turbulent boundary layer and the unmixed--stable--flow above or below it. We present a second-order closure able to cope simultaneously with both bulk and boundary layer regions, and we test it against high-resolution state-of-the-art two-dimensional numerical simulations in a convective and stratified belt for values of the Rayleigh number up to Ra ∼ 1010. Data are taken from a Rayleigh-Taylor system confined by the existence of an adiabatic gradient. [ABSTRACT FROM AUTHOR]

We introduce a variant of the Hybrid Monte Carlo (HMC) algorithm to address large-deviation statistics in stochastic hydrodynamics. Based on the path-integral approach to stochastic (partial) differential equations, our HMC algorithm samples space-time histories of the dynamical degrees of freedom under the influence of random noise. First, we validate and benchmark the HMC algorithm by reproducing multiscale properties of the one-dimensional Burgers equation driven by Gaussian and white-in-time noise. Second, we show how to implement an importance sampling protocol to significantly enhance, by orders of magnitudes, the probability to sample extreme and rare events, making it possible to estimate moments of field variables of extremely high order (up to 30 and more). By employing reweighting techniques, we map the biased configurations back to the original probability measure in order to probe their statistical importance. Finally, we show that by biasing the system towards very intense negative gradients, the HMC algorithm is able to explore the statistical fluctuations around instanton configurations. Our results will also be interesting and relevant in lattice gauge theory since they provide unique insights into reweighting techniques. [ABSTRACT FROM AUTHOR]

We present an analytical derivation of the transport coefficients of a relativistic gas in (2+1) dimensions for both Chapman-Enskog (CE) asymptotics and Grad's expansion methods. We further develop a systematic calibration method, connecting the relaxation time of relativistic kinetic theory to the transport parameters of the associated dissipative hydrodynamic equations. Comparison of our analytical results and numerical simulations shows that the CE method correctly captures dissipative effects, while Grad's method does not, in agreement with previous analyses performed in the (3+1)-dimensional case. These results provide a solid basis for accurately calibrated computational studies of relativistic dissipative flows. [ABSTRACT FROM AUTHOR]

Based on extensive numerical simulations, accounting for electrostatic interactions and dissipative electron-phonon scattering, we propose experimentally realizable geometries capable of sustaining electronic preturbulence in graphene samples. In particular, preturbulence is predicted to occur at experimentally attainable values of the Reynolds number between 10 and 50, over a broad spectrum of frequencies between 10 and 100 GHz. [ABSTRACT FROM AUTHOR]

Despite a long record of intense effort, the basic mechanisms by which dissipation emerges from the microscopic dynamics of a relativistic fluid still elude complete understanding. In particular, several details must still be finalized in the pathway from kinetic theory to hydrodynamics mainly in the derivation of the values of the transport coefficients. In this paper, we approach the problem by matching data from lattice-kinetic simulations with analytical predictions. Our numerical results provide neat evidence in favor of the Chapman-Enskog [The Mathematical Theory of Non-Uniform Gases, 3rd ed. (Cambridge University Press, Cambridge, U.K., 1970)] procedure as suggested by recent theoretical analyses along with qualitative hints at the basic reasons why the Chapman-Enskog expansion might be better suited than Grad's method [Commun. Pure Appl. Math. 2, 331 (1949)] to capture the emergence of dissipative effects in relativistic fluids. [ABSTRACT FROM AUTHOR]

We present a systematic derivation of relativistic lattice kinetic equations for finite-mass particles, reaching close to the zero-mass ultrarelativistic regime treated in the previous literature. Starting from an expansion of the Maxwell-Jüttner distribution on orthogonal polynomials, we perform a Gauss-type quadrature procedure and discretize the relativistic Boltzmann equation on space-filling Cartesian lattices. The model is validated through numerical comparison with standard tests and solvers in relativistic fluid dynamics such as Boltzmann approach multiparton scattering and previous relativistic lattice Boltzmann models. This work provides a significant step towards the formulation of a unified relativistic lattice kinetic scheme, covering both massive and near-massless particles regimes. [ABSTRACT FROM AUTHOR]

The correlation length ξ, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length. The scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer.

Experiments on spin glasses can now make precise measurements of the exponent z(T) governing the growth of glassy domains, while our computational capabilities allow us to make quantitative predictions for experimental scales. However, experimental and numerical values for z(T) have differed. We use new simulations on the Janus II computer to resolve this discrepancy, finding a time-dependent z(T,t_{w}), which leads to the experimental value through mild extrapolations. Furthermore, theoretical insight is gained by studying a crossover between the T=T_{c} and T=0 fixed points.

We first reproduce on the Janus and Janus II computers a milestone experiment that measures the spin-glass coherence length through the lowering of free-energy barriers induced by the Zeeman effect. Secondly, we determine the scaling behavior that allows a quantitative analysis of a new experiment reported in the companion Letter [S. Guchhait and R. Orbach, Phys. Rev. Lett. 118, 157203 (2017)].PRLTAO0031-900710.1103/PhysRevLett.118.157203 The value of the coherence length estimated through the analysis of microscopic correlation functions turns out to be quantitatively consistent with its measurement through macroscopic response functions. Further, nonlinear susceptibilities, recently measured in glass-forming liquids, scale as powers of the same microscopic length.

We study the off-equilibrium dynamics of the three-dimensional Ising spin glass in the presence of an external magnetic field. We have performed simulations both at fixed temperature and with an annealing protocol. Thanks to the Janus special-purpose computer, based on field-programmable gate array (FPGAs), we have been able to reach times equivalent to 0.01 s in experiments. We have studied the system relaxation both for high and for low temperatures, clearly identifying a dynamical transition point. This dynamical temperature is strictly positive and depends on the external applied magnetic field. We discuss different possibilities for the underlying physics, which include a thermodynamical spin-glass transition, a mode-coupling crossover, or an interpretation reminiscent of the random first-order picture of structural glasses.

We numerically study the aging properties of the dynamical heterogeneities in the Ising spin glass. We find that a phase transition takes place during the aging process. Statics-dynamics correspondence implies that systems of finite size in equilibrium have static heterogeneities that obey finite-size scaling, thus signaling an analogous phase transition in the thermodynamical limit. We compute the critical exponents and the transition point in the equilibrium setting, and use them to show that aging in dynamic heterogeneities can be described by a finite-time scaling ansatz, with potential implications for experimental work.

We study numerically the nonequilibrium dynamics of the Ising spin glass, for a time spanning 11 orders of magnitude, thus approaching the experimentally relevant scale (i.e., seconds). We introduce novel analysis techniques to compute the coherence length in a model-independent way. We present strong evidence for a replicon correlator and for overlap equivalence. The emerging picture is compatible with noncoarsening behavior.

We present different results from high-resolution high-statistics direct numerical simulations of a three-dimensional convective cell. We test the fundamental physical picture of the presence of both a Bolgiano-Obhukhov-like and a Kolmogorov-like regime. We find that the dimensional predictions for these two distinct regimes (characterized, respectively, by an active and passive role of the temperature field) are consistent with our analysis.