Transfer of a Passive Vector Admixture by a Two-Dimensional Turbulent Flow

We consider a model of a passive vector field transfer by a random two-dimensional transverse velocity field that is uncorrelated in time and has Gaussian spatial statistics given by a powerlike correlator. We use the renormalization group and the operator product expansion techniques to show that the asymptotic approximation of the structure functions of a vector field in the inertial range is determined by the energy dissipation fluctuations. The dependence of the asymptotic approximation on the external scale of turbulence is essential and has a powerlike form (the case of an anomalous scaling). The corresponding exponents are calculated in the one-loop approximation for structure functions of an arbitrary order.     

Renormalization group in the theory of turbulence: Three-loop approximation as <Emphasis Type="Italic">d</Emphasis> → ∞

We use the renormalization group method to study the stochastic Navier-Stokes equation with a random force correlator of the form k ^4−d−2ɛ in a d-dimensional space in connection with the problem of constructing a 1/d-expansion and going beyond the framework of the standard ɛ-expansion in the theory of fully developed hydrodynamic turbulence. We find a sharp decrease in the number of diagrams of the perturbation theory for the Green’s function in the large-d limit and develop a technique for calculating the diagrams analytically. We calculate the basic ingredients of the renormalization group approach (renormalization constant, β-function, fixed-point coordinates, and ultraviolet correction index ω) up to the order ɛ ³ (three-loop approximation). We use the obtained results to propose hypothetical exact expressions (i.e., not in the form of ɛ-expansions) for the fixed-point coordinate and the index ω. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 158, No. 3, pp. 460–477, March, 2009.     

Critical behavior of the Gaussian model on fractal lattices in external magnetic field

For inhomogeneous lattices we generalize the classical Gaussian model, i.e. it is proposed that the Gaussian type distribution constant and the external magnetic field of site i in this model depend on the coordination number q_i of site i, and that the relation holds among b_q's, where b_q is the Gaussian type distribution constant of site j. Using the decimation real-space renormalization group following the spin-rescaling method, the critical points and critical exponents of the Gaussian model are calculated on some Koch type curves and a family of the diamond-type hierarchical (or DH) lattices. At the critical points, it is found that the nearest-neighbor interaction and the magnetic field of site j can be expressed in the form respectively. It is also found that most critical exponents depend on the fractal dimensionality of a fractal system. For the family of the DH lattices, the results are identical with the exact results on translation symmetric lattices, and if the fractal dimensionalityd _f=4, the Gaussian model and the mean field theories give the same results.     

Anomalous scaling in statistical models of passively advected vector fields

We use the methods of the renormalization group and the operator product expansion to consider the problem of the stochastic advection of a passive vector field with the most general form of the nonlinear term allowed by the Galilean symmetry. The external velocity field satisfies the Navier-Stokes equation. We show that the correlation functions have anomalous scaling in the inertial range. The corresponding anomalous exponents are determined by the critical dimensions of tensor composite fields (operators) built from only the fields themselves. We calculate the anomalous dimensions in the leading order of the expansion in the exponent in the correlator of the external force in the Navier-Stokes equation (the oneloop approximation of the renormalization group). The anomalous exponents exhibit a hierarchy related to the anisotropy degree: the lower the rank of the tensor operator is, the lower its dimension. The leading asymptotic terms are determined by the scalar operators in both the isotropic and the anisotropic cases, which completely agrees with Kolmogorov’s hypothesis of local isotropy restoration.     

DNS and scaling law analysis of compressible turbulent channel flow

Fully developed compressible turbulent channel flow (Ma = 0.8, Re = 3300) is numerically simulated, and the data base of turbulence is established. The statistics such as density-weighted mean velocity and RMS velocity fluctuations in semi-local coordinates agree well with those from other DNS data. High order statistics (skewness and flatness factors) of velocity fluctuations of compressible turbulence are reported for the first time. Compressibility effects are also discussed. Pressure-dilatation absorbs part of the kinetic energy and makes the streaks of compressible channel flow more smooth. The scaling laws of compressible channel flow are also discussed. The conclusions are: (a) Scaling law is found in the center area of the channel, (b) In this area, ESS is also found, (c) When Mach number is not very high, compressibility has little effect on scaling exponents.     

Characterization of scaling functions in a multiresolution analysis

We characterize the scaling functions of a multiresolution analysis in a general context, where instead of the dyadic dilation one considers the dilation given by a fixed linear map such that and all (complex) eigenvalues of have absolute value greater than In the general case the conditions depend on the map We identify some maps for which the obtained condition is equivalent to the dyadic case, i.e., when is a diagonal matrix with all numbers in the diagonal equal to There are also easy examples of expanding maps for which the obtained condition is not compatible with the dyadic case. The complete characterization of the maps for which the obtained conditions are equivalent is out of the scope of the present note.

     

Hierarchical structure analysis of forced two-dimensional turbulence simulated by lattice BGK model

Introduction\\hether there exist universal statistic behaviors in small-scale turbulence of higll Re}'nolds-llunlber and what are their fundamenta1 forms have been an important issue both fOr under-standing physics of turbulence and fOr modeling turbu1ent flows in engineering. A powerfultool to these is direct numerical simulation. Besides spectra method and other high-precisiondifference schemes. lattice Bo1tzmann methodl1] developed in the past decade promises anothertechnique to simulate t…     

Renormalization group, operator expansion, and anomalous scaling in a simple model of turbulent diffusion

Using the renormalization group method and the operator expansion in the Obukhov-Kraichnan model that describes the intermixing of a passive scalar admixture by a random Gaussian field of velocities with the correlator 〈v(t,x)v(t′,x)〉−〈v(t,x)v(t′,x′)〉∝δ(t−t′)|x−x′|^ε, we prove that the anomalous scaling in the inertial interval is caused by the presence of “dangerous” composite operators (powers of the local dissipation rate) whose negative critical dimensions determine the anomalous exponents. These exponents are calculated up to the second order of the ε expansion. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 120, No. 2, pp. 309–314, August, 1999.     

Simulation of turbulent dispersion using a simple random model of the flow field

A method is devised to simulate the movement and spreading of a patch of contaminant in two-dimensional turbulent flow. The turbulent motion is exponentially divided into components of differing wave number, adjacent components being made to have correlation times differing by a factor of two. The turbulent motion is then reconstructed by replacing each component with a sinusoidal advection field having a randomly directed wave number. Contaminant particles are advected by each of the reconstructed components, the smallest scale components being applied first. A computer simulation was performed, using a Kolmogorov k^{-$\text{5}\text{3}$} turbulent energy spectrum. Batchelor's σ ∝ t^{$\text{3}\text{2}$} law for the spreading of a contaminant patch was reproduced, approximately, as was Richardson's non-Gaussian asymptotic form of the distance-neighbour function.     

A scaling limit for the length of the longest cycle in a sparse random digraph

We discuss the length of the longest directed cycle in the sparse random digraph , constant. We show that for large there exists a function such that a.s. The function where is a polynomial in . We are only able to explicitly give the values , although we could in principle compute any .     

Structure of Lorentzian tori with a killing vector field

All Lorentzian tori with a non-discrete group of isometries are characterized and explicitly obtained. They can lie into three cases: (a) flat, (b) conformally flat but non-flat, and (c) geodesically incomplete. A detailed study of many of their properties (including results on the logical dependence of the three kinds of causal completeness, on geodesic connectedness and on prescribed curvature) is carried out. The incomplete case is specially analyzed, and several known examples and results in the literature are generalized from a unified point of view.

     

Simulation of a turbulent premixed open V-shaped flame using contour advection with surgery

Despite its capability of high spatial resolution, simulation of turbulent flows with traditional Lagrangian (front tracking) scheme is often discouraged by numerical instability caused by clustering of marker nodes and topological changes of fronts. Contour advection surgery (CAS), being a robust front tracking scheme, can limit the growth of front complexity during simulation without jeopardizing accuracy or efficiency. This allows it to open up an advantage over traditional front-tracking schemes. It has already been demonstrated that CAS, with incorporation of the reaction sheet model, can accurately simulate the propagation and advection of a turbulent premixed V-shaped flame. In this study, it is further tested with 10 values of vortex circulation. A range of upstream turbulence levels of 1.8–19.8% was obtained. Results indicate that upstream turbulence increase the average flame length, flame zone area and the overall burning rate. Flame surface density Σ was also estimated. Maximum values of Σ obtained lie in the range 0.1–1.4 mm⁻¹. Skewness towards the burnt region was observed in all profiles of Σ. Similar to results from laboratory experiments, it was found that Σ values decreases with upstream turbulence. From this study, the ability of CAS to cope with intense turbulence is demonstrated and a better quantitative understanding on the scheme has also been acquired.     

A simple proof of a primal affine scaling method

In this paper, we present a simpler proof of the result of Tsuchiya and Muramatsu on the convergence of the primal affine scaling method. We show that the primal sequence generated by the method converges to the interior of the optimum face and the dual sequence to the analytic center of the optimal dual face, when the step size implemented in the procedure is bounded by 2/3. We also prove the optimality of the limit of the primal sequence for a slightly larger step size of 2q/(3q–1), whereq is the number of zero variables in the limit. We show this by proving the dual feasibility of a cluster point of the dual sequence.Partially supported by the grant CCR-9321550 from NSF.     

Simulation of turbulent magnetic reconnection in the small-scale solar wind

Some observational examples for the possible occurrence of the turbulent magnetic reconnection in the solar wind are found by analysing Helios spacecraft's high resolution data. The phenomena of turbulent magnetic reconnections in small scale solar wind are simulated by introducing a third order accuracy upwind compact difference scheme to the compressible two_dimensional MHD flow. Numerical results verify that the turbulent magnetic reconnection process could occur in small scale interplanetary solar wind, which is a basic feature characterizing the magnetic reconnection in high_magnetic Reynolds number (RM=2 000-10 000) solar wind. The configurations of the magnetic reconnection could evolve from a single X_line to a multiple X-line reconnection, exhibiting a complex picture of the formation, merging and evolution of magnetic islands, and finally the magnetic reconnection would evolve into a low_energy state. Its life_span of evolution is about one hour order of magnitude. Various magnetic and flow signatures are recorded in the numerical test for different evolution stages and along different crossing paths, which could in principle explain and confirm the observational samples from the Helios spacecraft. These results are helpful for revealing the basic physical processes in the solar wind turbulence.     

The heat transfer problem in the wall region of a turbulent boundary layer

The problem of heat transfer in the wall region of a turbulent boundary layer has been investigated. The resonant triad in the theory of hydrodynamic stability was used to obtain the velocity field induced by the coherent structures in the wall region of the turbulent boundary layer, while the small scale turbulence was represented by a simple model. By such a new approach of modeling, the 3-D temperature field is calculated, the mean temperature profile in the wall region and the Nusselt number characterizing the heat flux, which was found to be in good agreement with the experimental observations are obtained. The instantaneous temperature field had streaky structures, thus offering a mechanism for their generation found in numerical simulations. Project supported by the National Natural Science Foundation of China (Grant No. 19132011).     

Convergence analysis of the projective scaling algorithm based on a long-step homogeneous affine scaling algorithm

In this paper we give a new convergence analysis of a projective scaling algorithm. We consider a long-step affine scaling algorithm applied to a homogeneous linear programming problem obtained from the original linear programming problem. This algorithm takes a fixed fraction λ≤2/3 of the way towards the boundary of the nonnegative orthant at each iteration. The iteration sequence for the original problem is obtained by pulling back the homogeneous iterates onto the original feasible region with a conical projection, which generates the same search direction as the original projective scaling algorithm at each iterate. The recent convergence results for the long-step affine scaling algorithm by the authors are applied to this algorithm to obtain some convergence results on the projective scaling algorithm. Specifically, we will show (i) polynomiality of the algorithm with complexities of O(nL) and O(n ² L) iterations for λ<2/3 and λ=2/3, respectively; (ii) global covnergence of the algorithm when the optimal face is unbounded; (iii) convergence of the primal iterates to a relative interior point of the optimal face; (iv) convergence of the dual estimates to the analytic center of the dual optimal face; and (v) convergence of the reduction rate of the objective function value to 1−λ.     

Critical behavior of Gaussian model on diamond-type hierarchical lattices

It is proposed that the Gaussian type distribution constantb _qi in the Gaussian model depends on the coordination numberq _i of sitei, and that the relation holds amongb _qi ’s. The Gaussian model is then studied on a family of the diamond-type hierarchical (or DH) lattices, by the decimation real-space renormalization group following spin-rescaling method. It is found that the magnetic property of the Gaussian model belongs to the same universal class, and that the critical pointK* and the critical exponentv are given by and , respectively. Project supported by the National Natural Science Foundation of China (Grant No. 19775008), the National Basic Research Project supported by the National Natural Science Foundation of China (Grant No. 19775008), the National Basic Research Project supported by the National Natural Science Foundation of China (Grant No. 19775008), the National Basic Research     

Decoupling of the longitudinal polarization of the vector field in the massless Higgs-Kibble model

We show that the three-dimensionally longitudinal component of the vector field decouples in the massless limit of non-Abelian Higgs model.