Anomalous scaling regimes of a passive scalar advected by the synthetic velocity field

Amplification of the anomalous scaling in the Kazantsev-Kraichnan model with finite-time correlations and spatial parity violation

By using the field theoretic renormalization group technique together with the operator product expansion, simultaneous influence of the spatial parity violation and finite-time correlations of an electrically conductive turbulent environment on the inertial-range scaling behavior of correlation functions of a passively advected weak magnetic field is investigated within the corresponding generalized Kazantsev-Kraichnan model in the second order of the perturbation theory (in the two-loop approximation). The explicit dependence of the anomalous dimensions of the leading composite operators on the fixed point value of the parameter that controls the presence of finite-time correlations of the turbulent field as well as on the parameter that drives the amount of the spatial parity violation (helicity) in the system is found even in the case with the presence of the large-scale anisotropy. In accordance with the Kolmogorov's local isotropy restoration hypothesis, it is shown that, regardless of the amount of the spatial parity violation, the scaling properties of the model are always driven by the anomalous dimensions of the composite operators near the isotropic shell. The asymptotic (inertial-range) scaling form of all single-time two-point correlation functions of arbitrary order of the passively advected magnetic field is found. The explicit dependence of the corresponding scaling exponents on the helicity parameter as well as on the parameter that controls the finite-time velocity correlations is determined. It is shown that, regardless of the amount of the finite-time correlations of the given Gaussian turbulent environment, the presence of the spatial parity violation always leads to more negative values of the scaling exponents, i.e., to the more pronounced anomalous scaling of the magnetic correlation functions. At the same time, it is shown that the stronger the violation of spatial parity, the larger the anomalous behavior of magnetic correlations.

Anomalous scaling in kinematic magnetohydrodynamic turbulence: Two-loop anomalous dimensions of leading composite operators

Using the field theoretic formulation of the kinematic magnetohydrodynamic turbulence, the explicit expressions for the anomalous dimensions of leading composite operators, which govern the inertial-range scaling properties of correlation functions of the weak magnetic field passively advected by the electrically conductive turbulent environment driven by the Navier-Stokes velocity field, are derived and analyzed in the second order of the corresponding perturbation expansion (in the two-loop approximation). Their properties are compared to the properties of the same anomalous dimensions obtained in the framework of the Kazantsev-Kraichnan model of the kinematic magnetohydrodynamics with the Gaussian statistics of the turbulent velocity field as well as to the analogous anomalous dimensions of the leading composite operators in the problem of the passive scalar advection by the Gaussian (the Kraichnan model) and non-Gaussian (driven by the Navier-Stokes equation) turbulent velocity field. It is shown that, regardless of the Gaussian or non-Gaussian statistics of the turbulent velocity field, the two-loop corrections to the leading anomalous dimensions are much more important in the case of the problem of the passive advection of the vector (magnetic) field than in the case of the problem of the passive advection of scalar fields. At the same time, it is also shown that, in phenomenologically the most interesting case with three spatial dimensions, higher velocity correlations of the turbulent environment given by the Navier-Stokes velocity field play a rather limited role in the anomalous scaling of passive scalar as well as passive vector quantities, i.e., that the two-loop corrections to the corresponding leading anomalous dimensions are rather close to those obtained in the framework of the Gaussian models, especially as for the problem of scalar field advection. On the other hand, the role of the non-Gaussian statistics of the turbulent velocity field becomes dominant for higher spatial dimensions in the case of the kinematic magnetohydrodynamic turbulence but remains negligible in the problem of the passive scalar advection.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Strongly Nonlinear Diffusion in Turbulent Environment: A Problem with Infinitely Many Couplings

The field theoretic renormalization group is applied to the strongly nonlinear stochastic advection-diffusion equation. The turbulent advection is modelled by the Kazantsev–Kraichnan “rapid-change” ensemble. As a requirement of the renormalizability, the model necessarily involves infinite number of coupling constants (“charges”). The one-loop counterterm is calculated explicitly. The corresponding renormalization group equation demonstrates existence of a pair of two-dimensional surfaces of fixed points in the infinite-dimensional parameter space. If the surfaces contain infrared attractive regions, the problem allows for the large-scale, long-time scaling behaviour. For the first surface (advection is irrelevant), the critical dimensions of the scalar field Δθ, the response field Δθ′ and the frequency Δω are nonuniversal (through the dependence on the effective couplings) but satisfy certain exact identities. For the second surface (advection is relevant), the dimensions are universal and they are found exactly.

Evidence for enhancement of anisotropy persistence in kinematic magnetohydrodynamic turbulent systems with finite-time correlations

Using the field-theoretic renormalization group approach and the operator product expansion technique in the second order of the corresponding perturbative expansion, the influence of finite-time correlations of the turbulent velocity field on the scaling properties of the magnetic field correlation functions as well as on the anisotropy persistence deep inside the inertial range are investigated in the framework of the generalized Kazantsev-Kraichnan model of kinematic magnetohydrodynamic turbulence. Explicit two-loop expressions for the scaling exponents of the single-time two-point correlation functions of the magnetic field are derived and it is shown that the presence of the finite-time velocity correlations has a nontrivial impact on their inertial-range behavior and can lead, in general, to significantly more pronounced anomalous scaling of the magnetic field correlation functions in comparison to the rapid-change limit of the model, especially for the most interesting three-dimensional case. Moreover, by analyzing the asymptotic behavior of appropriate dimensionless ratios of the magnetic field correlation functions, it is also shown that the presence of finite-time correlations of the turbulent velocity field has a strong impact on the large-scale anisotropy persistence deep inside the inertial interval. Namely, it leads to a significant enhancement of the anisotropy persistence, again, especially in three spatial dimensions.

Effects of turbulent environment and random noise on self-organized critical behavior: Universality versus nonuniversality

Self-organized criticality in the Hwa-Kardar model of a "running sandpile" [Phys. Rev. Lett. 62, 1813 (1989)10.1103/PhysRevLett.62.1813; Phys. Rev. A 45, 7002 (1992)10.1103/PhysRevA.45.7002] with a turbulent motion of the environment taken into account is studied with the field theoretic renormalization group (RG). The turbulent flow is modeled by the synthetic d-dimensional generalization of the anisotropic Gaussian velocity ensemble with finite correlation time, introduced by Avellaneda and Majda [Commun. Math. Phys. 131, 381 (1990)10.1007/BF02161420; Commun. Math. Phys. 146, 139 (1992)10.1007/BF02099212]. The Hwa-Kardar model with time-independent (spatially quenched) random noise is considered alongside the original model with white noise. The aim of the present paper is to explore fixed points of the RG equations which determine the possible types of universality classes (regimes of critical behavior of the system) and critical dimensions of the measurable quantities. Our calculations demonstrate that influence of the type of random noise is extremely large: in contrast to the case of white noise where the system possesses three fixed points, the case of spatially quenched noise involves four fixed points with overlapping stability regions. This means that in the latter case the critical behavior of the system depends not only on the global parameters of the system, which is the usual case, but also on the initial values of the charges (coupling constants) of the system. These initial conditions determine the specific fixed point which will be reached by the RG flow. Since now the critical properties of the system are not defined strictly by its parameters, the situation may be interpreted as a universality violation. Such systems are not forbidden, but they are rather rare. It is especially interesting that the same model without turbulent motion of the environment does not predict this nonuniversal behavior and demonstrates the usual one with prescribed universality classes instead [J. Stat. Phys. 178, 392 (2020)10.1007/s10955-019-02436-8].

Renormalization group study of superfluid phase transition: Effect of compressibility

Dynamic critical behavior in superfluid systems is considered in the presence of external stirring and advecting processes. The latter are generated by means of the Gaussian random velocity ensemble with white-noise character in time variable and self-similar spatial dependence. The main focus of this work is to analyze an effect of compressible modes on the critical behavior. The model is formulated through stochastic Langevin equations, which are then recast into the Janssen-De Dominicis response formalism. Employing the field-theoretic perturbative renormalization group method we analyze large-scale properties of the model. Explicit calculations are performed to the leading one-loop approximation in the double (ɛ,y) expansion scheme, where ɛ is a deviation from the upper critical dimension d_{c}=4 and y describes a scaling property of the velocity ensemble. Altogether five distinct universality classes are expected to be macroscopically observable. In contrast to the incompressible case, we find that compressibility leads to an enhancement and stabilization of nontrivial asymptotic regimes.

Effects of Turbulent Environment on Self-Organized Critical Behavior: Isotropy vs. Anisotropy

We study a self-organized critical system under the influence of turbulent motion of the environment. The system is described by the anisotropic continuous stochastic equation proposed by Hwa and Kardar [Phys. Rev. Lett.62: 1813 (1989)]. The motion of the environment is modelled by the isotropic Kazantsev–Kraichnan “rapid-change” ensemble for an incompressible fluid: it is Gaussian with vanishing correlation time and the pair correlation function of the form ∝δ(t−t′)/kd+ξ, where k is the wave number and ξ is an arbitrary exponent with the most realistic values ξ=4/3 (Kolmogorov turbulence) and ξ→2 (Batchelor’s limit). Using the field-theoretic renormalization group, we find infrared attractive fixed points of the renormalization group equation associated with universality classes, i.e., with regimes of critical behavior. The most realistic values of the spatial dimension d=2 and the exponent ξ=4/3 correspond to the universality class of pure turbulent advection where the nonlinearity of the Hwa–Kardar (HK) equation is irrelevant. Nevertheless, the universality class where both the (anisotropic) nonlinearity of the HK equation and the (isotropic) advecting velocity field are relevant also exists for some values of the parameters ε=4−d and ξ. Depending on what terms (anisotropic, isotropic, or both) are relevant in specific universality class, different types of scaling behavior (ordinary one or generalized) are established.

Universality in incompressible active fluid: Effect of nonlocal shear stress

Phase transitions in active fluids attracted significant attention within the last decades. Recent results show [L. Chen et al., New J. Phys. 17, 042002 (2015)10.1088/1367-2630/17/4/042002] that an order-disorder phase transition in incompressible active fluids belongs to a new universality class. In this work, we further investigate this type of phase transition and focus on the effect of long-range interactions. This is achieved by introducing a nonlocal shear stress into the hydrodynamic description, which leads to superdiffusion of the velocity field, and can be viewed as a result of the active particles performing Lévy walks. The universal properties in the critical region are derived by performing a perturbative renormalization group analysis of the corresponding response functional within the one-loop approximation. We show that the effect of nonlocal shear stress decreases the upper critical dimension of the model, and can lead to the irrelevance of the active fluid self-advection with the resulting model belonging to an unusual long-range Model A universality class not reported before, to our knowledge. Moreover, when the degree of nonlocality is sufficiently high all nonlinearities become irrelevant and the mean-field description is valid in any spatial dimension.

Finite Time Correlations and Compressibility Effects in the Three-Dimensional Kraichnan Model

Using the field theoretic renormalization group technique the simultaneous influence of the compressibility and finite time correlations of the non-solenoidal Gaussian velocity field on the advection of a passive scalar field is studied within the generalized Kraichnan model in three spatial dimensions up to the second-order approximation in the corresponding perturbative expansion. All possible infrared stable fixed points of the model, which drive the corresponding scaling regimes of the model, are identified and their regions of the infrared stability in the model parametric space are discussed. It is shown that, depending on the value of the parameter that drives the compressibility of the system, there exists a gap in the parametric space between the regions where the model with the frozen velocity field and the model with finite-time correlations of the velocity field are stable or there exists an overlap between them.

Symmetry Breaking in Stochastic Dynamics and Turbulence

Symmetries play paramount roles in dynamics of physical systems. All theories of quantum physics and microworld including the fundamental Standard Model are constructed on the basis of symmetry principles. In classical physics, the importance and weight of these principles are the same as in quantum physics: dynamics of complex nonlinear statistical systems is straightforwardly dictated by their symmetry or its breaking, as we demonstrate on the example of developed (magneto)hydrodynamic turbulence and the related theoretical models. To simplify the problem, unbounded models are commonly used. However, turbulence is a mesoscopic phenomenon and the size of the system must be taken into account. It turns out that influence of outer length of turbulence is significant and can lead to intermittency. More precisely, we analyze the connection of phenomena such as behavior of statistical correlations of observable quantities, anomalous scaling, and generation of magnetic field by hydrodynamic fluctuations with symmetries such as Galilean symmetry, isotropy, spatial parity and their violation and finite size of the system.

Perspectives on scaling and multiscaling in passive scalar turbulence

We revisit the well-known problem of multiscaling in substances passively advected by homogeneous and isotropic turbulent flows or passive scalar turbulence. To that end we propose a two-parameter continuum hydrodynamic model for an advected substance concentration θ, parametrized jointly by y and y[over ¯], that characterize the spatial scaling behavior of the variances of the advecting stochastic velocity and the stochastic additive driving force, respectively. We analyze it within a one-loop dynamic renormalization group method to calculate the multiscaling exponents of the equal-time structure functions of θ. We show how the interplay between the advective velocity and the additive force may lead to simple scaling or multiscaling. In one limit, our results reduce to the well-known results from the Kraichnan model for passive scalar. Our framework of analysis should be of help for analytical approaches for the still intractable problem of fluid turbulence itself.

Statistical symmetry restoration in fully developed turbulence: Renormalization group analysis of two models

In this paper we consider the model of incompressible fluid described by the stochastic Navier-Stokes equation with finite correlation time of a random force. Inertial-range asymptotic behavior of fully developed turbulence is studied by means of the field theoretic renormalization group within the one-loop approximation. It is corroborated that regardless of the values of model parameters and initial data the inertial-range behavior of the model is described by the limiting case of vanishing correlation time. This indicates that the Galilean symmetry of the model violated by the "colored" random force is restored in the inertial range. This regime corresponds to the only nontrivial fixed point of the renormalization group equation. The stability of this point depends on the relation between the exponents in the energy spectrum E∝k^{1-y} and the dispersion law ω∝k^{2-η}. The second analyzed problem is the passive advection of a scalar field by this velocity ensemble. Correlation functions of the scalar field exhibit anomalous scaling behavior in the inertial-convective range. We demonstrate that in accordance with Kolmogorov's hypothesis of the local symmetry restoration the main contribution to the operator product expansion is given by the isotropic operator, while anisotropic terms should be considered only as corrections.

Influence of turbulent mixing on critical behavior of directed percolation process: Effect of compressibility

Universal behavior is a typical emergent feature of critical systems. A paramount model of the nonequilibrium critical behavior is the directed bond percolation process that exhibits an active-to-absorbing state phase transition in the vicinity of a percolation threshold. Fluctuations of the ambient environment might affect or destroy the universality properties completely. In this work, we assume that the random environment can be described by means of compressible velocity fluctuations. Using field-theoretic models and renormalization group methods, we investigate large-scale and long-time behavior. Altogether, 11 universality classes are found, out of which 4 are stable in the infrared limit and thus macroscopically accessible. In contrast to the model without velocity fluctuations, a possible candidate for a realistic three-dimensional case, a regime with relevant short-range noise, is identified. Depending on the dimensionality of space and the structure of the turbulent flow, we calculate critical exponents of the directed percolation process. In the limit of the purely transversal random force, critical exponents comply with the incompressible results obtained by previous authors. We have found intriguing nonuniversal behavior related to the mutual effect of compressibility and advection.

Passive scalar transport by a non-Gaussian turbulent flow in the Batchelor regime

We analyze passive scalar advection by a turbulent flow in the Batchelor regime. No restrictions on the velocity statistics of the flow are assumed. The properties of the scalar are derived from the statistical properties of velocity; analytic expressions for the moments of scalar density are obtained. We show that the scalar statistics can differ significantly from that obtained in the frames of the Kraichnan model.

Simultaneous influence of helicity and compressibility on anomalous scaling of the magnetic field in the Kazantsev-Kraichnan model

Using the field theoretic renormalization group technique and the operator product expansion, the systematic investigation of the influence of the spatial parity violation on the anomalous scaling behavior of correlation functions of the weak passive magnetic field in the framework of the compressible Kazantsev-Kraichnan model with the presence of a large-scale anisotropy is performed up to the second order of the perturbation theory (two-loop approximation). The renormalization group analysis of the model is done and the two-loop explicit expressions for the anomalous and critical dimensions of the leading composite operators are found as functions of the helicity and compressibility parameters and their anisotropic hierarchies are discussed. It is shown that for arbitrary values of the helicity parameter and for physically acceptable (small enough) values of the compressibility parameter, the main role is played by the composite operators near the isotropic shell in accordance with the Kolmogorov's local isotropy restoration hypothesis. The anomalous dimensions of the relevant composite operators are then compared with the anomalous dimensions of the corresponding leading composite operators in the Kraichnan model of passively advected scalar field. The significant difference between these two sets of anomalous dimensions is discussed. The two-loop inertial-range scaling exponents of the single-time two-point correlation functions of the magnetic field are found and their dependence on the helicity and compressibility parameters is studied in detail. It is shown that while the presence of the helicity leads to more pronounced anomalous scaling for correlation functions of arbitrary order, the compressibility, in general, makes the anomalous scaling more pronounced in comparison to the incompressible case only for low-order correlation functions. The persistence of the anisotropy deep inside the inertial interval is investigated using the appropriate odd ratios of the correlation functions. It is shown that, in general, the persistence of the anisotropy is much more pronounced in the helical systems, while in the compressible turbulent environments this is true only for low-order odd ratios of the correlation functions.

Turbulent compressible fluid: Renormalization group analysis, scaling regimes, and anomalous scaling of advected scalar fields

We study a model of fully developed turbulence of a compressible fluid, based on the stochastic Navier-Stokes equation, by means of the field-theoretic renormalization group. In this approach, scaling properties are related to the fixed points of the renormalization group equations. Previous analysis of this model near the real-world space dimension 3 identified a scaling regime [N. V. Antonov et al., Theor. Math. Phys. 110, 305 (1997)TMPHAH0040-577910.1007/BF02630456]. The aim of the present paper is to explore the existence of additional regimes, which could not be found using the direct perturbative approach of the previous work, and to analyze the crossover between different regimes. It seems possible to determine them near the special value of space dimension 4 in the framework of double y and ɛ expansion, where y is the exponent associated with the random force and ɛ=4-d is the deviation from the space dimension 4. Our calculations show that there exists an additional fixed point that governs scaling behavior. Turbulent advection of a passive scalar (density) field by this velocity ensemble is considered as well. We demonstrate that various correlation functions of the scalar field exhibit anomalous scaling behavior in the inertial-convective range. The corresponding anomalous exponents, identified as scaling dimensions of certain composite fields, can be systematically calculated as a series in y and ɛ. All calculations are performed in the leading one-loop approximation.

Directed percolation process in the presence of velocity fluctuations: Effect of compressibility and finite correlation time

The direct bond percolation process (Gribov process) is studied in the presence of random velocity fluctuations generated by the Gaussian self-similar ensemble with finite correlation time. We employ the renormalization group in order to analyze a combined effect of the compressibility and finite correlation time on the long-time behavior of the phase transition between an active and an absorbing state. The renormalization procedure is performed to the one-loop order. Stable fixed points of the renormalization group and their regions of stability are calculated in the one-loop approximation within the three-parameter (ɛ,y,η) expansion. Different regimes corresponding to the rapid-change limit and frozen velocity field are discussed, and their fixed points' structure is determined in numerical fashion.

Anomalous scaling in magnetohydrodynamic turbulence: Effects of anisotropy and compressibility in the kinematic approximation

The field-theoretic renormalization group and the operator product expansion are applied to the model of passive vector (magnetic) field advected by a random turbulent velocity field. The latter is governed by the Navier-Stokes equation for compressible fluid, subject to external random force with the covariance ∝ δ(t-t')k(4-d-y), where d is the dimension of space and y is an arbitrary exponent. From physics viewpoints, the model describes magnetohydrodynamic turbulence in the so-called kinematic approximation, where the effects of the magnetic field on the dynamics of the fluid are neglected. The original stochastic problem is reformulated as a multiplicatively renormalizable field-theoretic model; the corresponding renormalization group equations possess an infrared attractive fixed point. It is shown that various correlation functions of the magnetic field and its powers demonstrate anomalous scaling behavior in the inertial-convective range already for small values of y. The corresponding anomalous exponents, identified with scaling (critical) dimensions of certain composite fields ("operators" in the quantum-field terminology), can be systematically calculated as series in y. The practical calculation is performed in the leading one-loop approximation, including exponents in anisotropic contributions. It should be emphasized that, in contrast to Gaussian ensembles with finite correlation time, the model and the perturbation theory presented here are manifestly Galilean covariant.

Passive advection of a vector field: Anisotropy, finite correlation time, exact solution, and logarithmic corrections to ordinary scaling

In this work we study the generalization of the problem considered in [Phys. Rev. E 91, 013002 (2015)] to the case of finite correlation time of the environment (velocity) field. The model describes a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow. Inertial-range asymptotic behavior is studied by means of the field theoretic renormalization group and the operator product expansion. The advecting velocity field is Gaussian, with finite correlation time and preassigned pair correlation function. Due to the presence of distinguished direction n, all the multiloop diagrams in this model vanish, so that the results obtained are exact. The inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite correlation time) that correspond to the two nontrivial fixed points of the RG equations. Their stability depends on the relation between the exponents in the energy spectrum E∝k(⊥)(1-ξ) and the dispersion law ω∝k(⊥)(2-η). In contrast to the well-known isotropic Kraichnan's model, where various correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, here the corrections to ordinary scaling are polynomials of logarithms of the integral turbulence scale L.

Anomalous scaling of passive scalar fields advected by the Navier-Stokes velocity ensemble: effects of strong compressibility and large-scale anisotropy

The field theoretic renormalization group and the operator product expansion are applied to two models of passive scalar quantities (the density and the tracer fields) advected by a random turbulent velocity field. The latter is governed by the Navier-Stokes equation for compressible fluid, subject to external random force with the covariance ∝δ(t-t')k(4-d-y), where d is the dimension of space and y is an arbitrary exponent. The original stochastic problems are reformulated as multiplicatively renormalizable field theoretic models; the corresponding renormalization group equations possess infrared attractive fixed points. It is shown that various correlation functions of the scalar field, its powers and gradients, demonstrate anomalous scaling behavior in the inertial-convective range already for small values of y. The corresponding anomalous exponents, identified with scaling (critical) dimensions of certain composite fields ("operators" in the quantum-field terminology), can be systematically calculated as series in y. The practical calculation is performed in the leading one-loop approximation, including exponents in anisotropic contributions. It should be emphasized that, in contrast to Gaussian ensembles with finite correlation time, the model and the perturbation theory presented here are manifestly Galilean covariant. The validity of the one-loop approximation and comparison with Gaussian models are briefly discussed.

Anomalous scaling of a passive scalar field near two dimensions

The anomalous scaling of the single-time structure functions of a passive scalar field advected by the velocity field governed by the stochastic Navier-Stokes equation is investigated by the field theoretic renormalization group and the operator-product expansion with inclusion of additional ultraviolet divergences related to the spatial dimension d=2. Some modification of the standard minimal subtraction scheme is used to calculate the turbulent Prandtl number and the anomalous exponents within the one-loop approximation of a perturbation theory. It is shown that the taking into account of these additional divergences is in full agreement with results obtained in the framework of the usual analytical expansion technique, which are valid for d>2.

Anomalous scaling and large-scale anisotropy in magnetohydrodynamic turbulence: two-loop renormalization-group analysis of the Kazantsev-Kraichnan kinematic model

The field theoretic renormalization group and operator product expansion are applied to the Kazantsev-Kraichnan kinematic model for the magnetohydrodynamic turbulence. The anomalous scaling emerges as a consequence of the existence of certain composite fields ("operators") with negative dimensions. The anomalous exponents for the correlation functions of arbitrary order are calculated in the two-loop approximation (second order of the renormalization-group expansion), including the anisotropic sectors. The anomalous scaling and the hierarchy of anisotropic contributions become stronger due to those second-order contributions.

Influence of anisotropy on anomalous scaling of a passive scalar advected by the Navier-Stokes velocity field

The influence of weak uniaxial small-scale anisotropy on the stability of the scaling regime and on the anomalous scaling of the single-time structure functions of a passive scalar advected by the velocity field governed by the stochastic Navier-Stokes equation is investigated by the field theoretic renormalization group and operator-product expansion within one-loop approximation of a perturbation theory. The explicit analytical expressions for coordinates of the corresponding fixed point of the renormalization-group equations as functions of anisotropy parameters are found, the stability of the three-dimensional Kolmogorov-like scaling regime is demonstrated, and the dependence of the borderline dimension d(c) is an element of (2,3] between stable and unstable scaling regimes is found as a function of the anisotropy parameters. The dependence of the turbulent Prandtl number on the anisotropy parameters is also briefly discussed. The influence of weak small-scale anisotropy on the anomalous scaling of the structure functions of a passive scalar field is studied by the operator-product expansion and their explicit dependence on the anisotropy parameters is present. It is shown that the anomalous dimensions of the structure functions, which are the same (universal) for the Kraichnan model, for the model with finite time correlations of the velocity field, and for the model with the advection by the velocity field driven by the stochastic Navier-Stokes equation in the isotropic case, can be distinguished by the assumption of the presence of the small-scale anisotropy in the systems even within one-loop approximation. The corresponding comparison of the anisotropic anomalous dimensions for the present model with that obtained within the Kraichnan rapid-change model is done.

Self-consistent four-point closure for transport in steady random flows

Ensemble averaging of advection-dispersion equations describing transport of a passive scalar in incompressible random velocity fields requires a closure approximation. Commonly used approaches, such as the direct interaction approximation and large-eddy simulations as well as equivalent renormalization schemes, employ so-called two-point (or one-loop) closures. These approaches have proven to be adequate for transport in zero-mean (unbiased) time-dependent random velocity fields with increasing accuracy for decreasing temporal coherence. In the opposite limit of steady velocity fields with finite bias, however, these schemes fail to predict effective transport properties both quantitatively and qualitatively, leading to an obvious inconsistency for transverse dispersion in two spatial dimensions. For this case, two-point closures predict that macroscopic transverse dispersion increases as the square root of the disorder variance while it has been proven rigorously that there is no disorder-induced contribution to macroscopic transverse dispersion for purely advective transport. Furthermore, two-point closures significantly underestimate the disorder-induced contribution to longitudinal dispersion. We derive a four-point closure for stochastically averaged transport equations that goes beyond classical one-loop schemes and demonstrate that it is exact for transverse dispersion and correctly predicts an increase of the longitudinal disorder-induced dispersion coefficient with the square of the variance of the strong disorder. The predicted values of asymptotic longitudinal dispersion coefficients are consistent with those obtained via Monte Carlo random walk simulations.

Anomalous scaling of a passive scalar advected by a turbulent velocity field with finite correlation time and uniaxial small-scale anisotropy

The influence of uniaxial small-scale anisotropy on the stability of the scaling regimes and on the anomalous scaling of the structure functions of a passive scalar advected by a Gaussian solenoidal velocity field with finite correlation time is investigated by the field theoretic renormalization group and operator product expansion within one-loop approximation. Possible scaling regimes are found and classified in the plane of exponents epsilon-eta , where epsilon characterizes the energy spectrum of the velocity field in the inertial range E proportional, variantk;{1-2epsilon} , and eta is related to the correlation time of the velocity field at the wave number k which is scaled as k;{-2+eta} . It is shown that the presence of anisotropy does not disturb the stability of the infrared fixed points of the renormalization group equations, which are directly related to the corresponding scaling regimes. The influence of anisotropy on the anomalous scaling of the structure functions of the passive scalar field is studied as a function of the fixed point value of the parameter u , which represents the ratio of turnover time of scalar field and velocity correlation time. It is shown that the corresponding one-loop anomalous dimensions, which are the same (universal) for all particular models with a concrete value of u in the isotropic case, are different (nonuniversal) in the case with the presence of small-scale anisotropy and they are continuous functions of the anisotropy parameters, as well as the parameter u . The dependence of the anomalous dimensions on the anisotropy parameters of two special limits of the general model, namely, the rapid-change model and the frozen velocity field model, are found when u-->infinity and u-->0 , respectively.

Influence of helicity on anomalous scaling of a passive scalar advected by the turbulent velocity field with finite correlation time: two-loop approximation

The influence of helicity on the stability of scaling regimes, on the effective diffusivity, and on the anomalous scaling of structure functions of a passive scalar advected by a Gaussian solenoidal velocity field with finite correlation time is investigated by the field theoretic renormalization group and operator-product expansion within the two-loop approximation. The influence of helicity on the scaling regimes is discussed and shown in the plane of exponents epsilon-eta, where epsilon characterizes the energy spectrum of the velocity field in the inertial range E proportional to k(1-2epsilon), and eta is related to the correlation time at the wave number k, which is scaled as k(-2+eta). The restrictions given by nonzero helicity on the regions with stable fixed points that correspond to the scaling regimes are analyzed in detail. The dependence of the effective diffusivity on the helicity parameter is discussed. The anomalous exponents of the structure functions of the passive scalar field which define their anomalous scaling are calculated and it is shown that, although the separate composite operators which define them strongly depend on the helicity parameter, the resulting two-loop contributions to the critical dimensions of the structure functions are independent of helicity. Details of calculations are shown.

Anomalous scaling of passively advected magnetic field in the presence of strong anisotropy

Inertial-range scaling behavior of high-order (up to order N=51 ) two-point correlation functions of a passively advected vector field has been analyzed in the framework of the rapid-change model with strong small-scale anisotropy with the aid of the renormalization group and the operator-product expansion. Exponents of the power-like asymptotic behavior of the correlation functions have been calculated in the one-loop approximation. These exponents are shown to depend on anisotropy parameters in such a way that a specific hierarchy related to the degree of anisotropy is observed. Deviations from power-law behavior like oscillations or logarithmic behavior in the corrections to correlation functions have not been found.

Anomalous scaling of a passive scalar advected by the Navier-Stokes velocity field: two-loop approximation

The field theoretic renormalization group and operator-product expansion are applied to the model of a passive scalar quantity advected by a non-Gaussian velocity field with finite correlation time. The velocity is governed by the Navier-Stokes equation, subject to an external random stirring force with the correlation function proportional to delta(t- t')k(4-d-2epsilon). It is shown that the scalar field is intermittent already for small epsilon, its structure functions display anomalous scaling behavior, and the corresponding exponents can be systematically calculated as series in epsilon. The practical calculation is accomplished to order epsilon2 (two-loop approximation), including anisotropic sectors. As for the well-known Kraichnan rapid-change model, the anomalous scaling results from the existence in the model of composite fields (operators) with negative scaling dimensions, identified with the anomalous exponents. Thus the mechanism of the origin of anomalous scaling appears similar for the Gaussian model with zero correlation time and the non-Gaussian model with finite correlation time. It should be emphasized that, in contrast to Gaussian velocity ensembles with finite correlation time, the model and the perturbation theory discussed here are manifestly Galilean covariant. The relevance of these results for real passive advection and comparison with the Gaussian models and experiments are briefly discussed.

Turbulence with pressure: anomalous scaling of a passive vector field

The field theoretic renormalization group (RG) and the operator-product expansion are applied to the model of a transverse (divergence-free) vector quantity, passively advected by the "synthetic" turbulent flow with a finite (and not small) correlation time. The vector field is described by the stochastic advection-diffusion equation with the most general form of the inertial nonlinearity; it contains as special cases the kinematic dynamo model, linearized Navier-Stokes (NS) equation, the special model without the stretching term that possesses additional symmetries and has a close formal resemblance with the stochastic NS equation. The statistics of the advecting velocity field is Gaussian, with the energy spectrum E(k) proportional to k(1-epsilon) and the dispersion law omega proportional to k(-2+eta), k being the momentum (wave number). The inertial-range behavior of the model is described by seven regimes (or universality classes) that correspond to nontrivial fixed points of the RG equations and exhibit anomalous scaling. The corresponding anomalous exponents are associated with the critical dimensions of tensor composite operators built solely of the passive vector field, which allows one to construct a regular perturbation expansion in epsilon and eta; the actual calculation is performed to the first order (one-loop approximation), including the anisotropic sectors. Universality of the exponents, their (in)dependence on the forcing, effects of the large-scale anisotropy, compressibility, and pressure are discussed. In particular, for all the scaling regimes the exponents obey a hierarchy related to the degree of anisotropy: the more anisotropic is the contribution of a composite operator to a correlation function, the faster it decays in the inertial range. The relevance of these results for the real developed turbulence described by the stochastic NS equation is discussed.

Anomalous scaling of a passive scalar advected by the turbulent velocity field with finite correlation time: two-loop approximation

The renormalization group and operator product expansion are applied to the model of a passive scalar quantity advected by the Gaussian self-similar velocity field with finite, and not small, correlation time. The inertial-range energy spectrum of the velocity is chosen in the form E(k) proportional, variant k(1-2 epsilon ), and the correlation time at the wave number k scales as k(-2+eta). Inertial-range anomalous scaling for the structure functions and other correlation functions emerges as a consequence of the existence in the model of composite operators with negative scaling dimensions, identified with anomalous exponents. For eta> epsilon, these exponents are the same as in the rapid-change limit of the model; for eta< epsilon, they are the same as in the limit of a time-independent (quenched) velocity field. For epsilon =eta (local turnover exponent), the anomalous exponents are nonuniversal through the dependence on a dimensionless parameter, the ratio of the velocity correlation time, and the scalar turnover time. The nonuniversality reveals itself, however, only in the second order of the epsilon expansion and the exponents are derived to order epsilon (2), including anisotropic contributions. It is shown that, for moderate order of the structure function n, and the space dimensionality d, finite correlation time enhances the intermittency in comparison with both the limits: the rapid-change and quenched ones. The situation changes when n and/or d become large enough: the correction to the rapid-change limit due to the finite correlation time is positive (that is, the anomalous scaling is suppressed), it is maximal for the quenched limit and monotonically decreases as the correlation time tends to zero.

Calculation of the anomalous exponents in the rapid-change model of passive scalar advection to order epsilon(3)

The field theoretic renormalization group and operator product expansion are applied to the model of a passive scalar advected by the Gaussian velocity field with zero mean and correlation function approximately equal to delta(t-t('))/k(d + epsilon). Inertial-range anomalous exponents, identified with the critical dimensions of various scalar and tensor composite operators constructed of the scalar gradients, are calculated within the epsilon expansion to order epsilon(3) (three-loop approximation), including the exponents in anisotropic sectors. The main goal of the paper is to give the complete derivation of this third-order result, and to present and explain in detail the corresponding calculational techniques. The character and convergence properties of the epsilon expansion are discussed, the improved "inverse" epsilon expansion is proposed, and the comparison with the existing nonperturbative results is given.

Anomalous scaling, nonlocality, and anisotropy in a model of the passively advected vector field

A model of the passive vector quantity advected by the Gaussian velocity field with the covariance approximately delta(t-t('))|x-x(')|(epsilon) is studied; the effects of pressure and large-scale anisotropy are discussed. The inertial-range behavior of the pair correlation function is described by an infinite family of scaling exponents, which satisfy exact transcendental equations derived explicitly in d dimensions by means of the functional techniques. The exponents are organized in a hierarchical order according to their degree of anisotropy, with the spectrum unbounded from above and the leading (minimal) exponent coming from the isotropic sector. This picture extends to higher-order correlation functions. Like in the scalar model, the second-order structure function appears nonanomalous and is described by the simple dimensional exponent: S2 approximately r(2-epsilon). For the higher-order structure functions, S(2n)approximately r(n(2-epsilon)+ delta(n)), the anomalous scaling behavior is established as a consequence of the existence in the corresponding operator product expansions of "dangerous" composite operators, whose negative critical dimensions determine the anomalous exponents delta(n)<0. A close formal resemblance of the model with the stirred Navier-Stokes equation reveals itself in the mixing of relevant operators and is the main motivation of the paper. Using the renormalization group, the anomalous exponents are calculated in the O(epsilon) approximation, in large d dimensions, for the even structure functions up to the twelfth order.

Anomalous scaling in two models of passive scalar advection: effects of anisotropy and compressibility

The problem of the effects of compressibility and large-scale anisotropy on anomalous scaling behavior is considered for two models describing passive advection of scalar density and tracer fields. The advecting velocity field is Gaussian, delta correlated in time, and scales with a positive exponent epsilon. Explicit inertial-range expressions for the scalar correlation functions are obtained; they are represented by superpositions of power laws with nonuniversal amplitudes and universal anomalous exponents (dependent only on epsilon and alpha, the compressibility parameter). The complete set of anomalous exponents for the pair correlation functions is found nonperturbatively, in any space dimension d, using the zero-mode technique. For higher-order correlation functions, the anomalous exponents are calculated to O(epsilon(2)) using the renormalization group. As in the incompressible case, the exponents exhibit a hierarchy related to the degree of anisotropy: the leading contributions to the even correlation functions are given by the exponents from the isotropic shell, in agreement with the idea of restored small-scale isotropy. As the degree of compressibility increases, the corrections become closer to the leading terms. The small-scale anisotropy reveals itself in the odd ratios of correlation functions: the skewness factor slowly decreases going down to small scales for the incompressible case, but starts to increase if alpha is large enough. The higher odd dimensionless ratios (hyperskewness, etc.) increase, thus signaling persistent small-scale anisotropy; this effect becomes more pronounced for larger values of alpha.

Anomalous exponents to order epsilon 3 in the rapid-change model of passive scalar advection

Field-theoretic renormalization group is applied to the Kraichnan model of a passive scalar advected by the Gaussian velocity field with the covariance <v(t,x)v(t('),x)>-<v(t,x)v(t('),x('))>~delta(t-t('))|x-x(')|(epsilon). Inertial-range anomalous exponents, related to the scaling dimensions of tensor composite operators built of the scalar gradients, are calculated to the order epsilon(3) of the epsilon expansion. The nature and the convergence of the epsilon expansion in the models of turbulence are briefly discussed.

Anomalous scaling of a passive scalar in the presence of strong anisotropy

Field theoretic renormalization group and the operator product expansion are applied to a model of a passive scalar quantity straight theta(t,x), advected by the Gaussian strongly anisotropic velocity field with the covariance infinity delta(t-t('))/x-x(')/(epsilon). Inertial-range anomalous scaling behavior is established, and explicit asymptotic expressions for the structure functions S(n)(r) identical with<[straight theta(t,x+r)-straight theta(t,x)](n)> are obtained. They are represented by superpositions of power laws; the corresponding anomalous exponents, which depend explicitly on the anisotropy parameters, are calculated to the first order in epsilon in any space dimension d. In the limit of vanishing anisotropy, the exponents are associated with tensor composite operators built of the scalar gradients, and exhibit a kind of hierarchy related to the degree of anisotropy: the less is the rank, the less is the dimension and, consequently, the more important is the contribution to the inertial-range behavior. The leading terms of the even (odd) structure functions are given by the scalar (vector) operators. For the finite anisotropy, the exponents cannot be associated with individual operators (which are essentially "mixed" in renormalization), but the aforementioned hierarchy survives for all the cases studied. The second-order structure function S2 is studied in more detail using the renormalization group and zero-mode techniques; the corresponding exponents and amplitudes are calculated within the perturbation theories in epsilon, 1/d, and in the anisotropy parameters. If the anisotropy of the velocity is strong enough, the skewness factor S(3)/S(3/2)(2) increases going down towards the depth of the inertial range; the higher-order odd ratios increase even if the anisotropy is weak.

Manifestation of anisotropy persistence in the hierarchies of magnetohydrodynamical scaling exponents

An example of a turbulent system where the failure of the hypothesis of small-scale isotropy restoration is detectable both in the "flattening" of the inertial-range scaling exponent hierarchy and in the behavior of odd-order dimensionless ratios, e.g., skewness and hyperskewness, is presented. Specifically, within the kinematic approximation in magnetohydrodynamical turbulence, we show that for compressible flows, the isotropic contribution to the scaling of magnetic correlation functions and the first anisotropic ones may become practically indistinguishable. Moreover, the skewness factor now diverges as the Peclet number goes to infinity, a further indication of small-scale anisotropy.

Persistence of small-scale anisotropies and anomalous scaling in a model of magnetohydrodynamics turbulence

The problem of anomalous scaling in magnetohydrodynamics turbulence is considered within the framework of the kinematic approximation, in the presence of a large-scale background magnetic field. The velocity field is Gaussian, delta-correlated in time, and scales with a positive exponent xi. Explicit inertial-range expressions for the magnetic correlation functions are obtained; they are represented by superpositions of power laws with nonuniversal amplitudes and universal (independent of the anisotropy and forcing) anomalous exponents. The complete set of anomalous exponents for the pair correlation function is found nonperturbatively, in any space dimension d, using the zero-mode technique. For higher-order correlation functions, the anomalous exponents are calculated to O(xi) using the renormalization group. The exponents exhibit a hierarchy related to the degree of anisotropy; the leading contributions to the even correlation functions are given by the exponents from the isotropic shell, in agreement with the idea of restored small-scale isotropy. Conversely, the small-scale anisotropy reveals itself in the odd correlation functions: the skewness factor is slowly decreasing going down to small scales and higher odd dimensionless ratios (hyperskewness, etc.) dramatically increase, thus diverging in the r-->0 limit.