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 $\propto \delta(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.

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.

A system of stochastic differential equations for the velocity and density of a classical self-gravitating matter is investigated by means of the field theoretic renormalization group. The existence of two types of large-scale scaling behavior, associated to physically admissible fixed points of the renormalization-group equations, is established. Their regions of stability are identified and the corresponding scaling dimensions are calculated in the one-loop approximation (first order of the epsilon expansion). The velocity and density fields have independent scaling dimensions. Our analysis supports the importance of the rotational (nonpotential) components of the velocity field in the formation of those scaling laws.