Universality and scaling behavior of injected power in elastic turbulence in wormlike micellar gel

We study the statistical properties of spatially averaged global injected power fluctuations for Taylor-Couette flow of a wormlike micellar gel formed by surfactant cetyltrimethylammonium tosylate. At sufficiently high Weissenberg numbers the shear rate, and hence the injected power p(t), at a constant applied stress shows large irregular fluctuations in time. The nature of the probability distribution function (PDF) of p(t) and the power-law decay of its power spectrum are very similar to that observed in recent studies of elastic turbulence for polymer solutions. Remarkably, these non-Gaussian PDFs can be well described by a universal, large deviation functional form given by the generalized Gumbel distribution observed in the context of spatially averaged global measures in diverse classes of highly correlated systems. We show by in situ rheology and polarized light scattering experiments that in the elastic turbulent regime the flow is spatially smooth but random in time, in agreement with a recent hypothesis for elastic turbulence.

On the k(1)(-1) scaling in sink-flow turbulent boundary layers

Scaling of the streamwise velocity spectrum phi(11)(k(1)) in the so-called sink-flow turbulent boundary layer is investigated in this work. The present experiments show strong evidence for the k(1)(-1) scaling i.e. phi(11)(k(1)) = Lambda(1)U(tau)(2)k(1)(-1), where k(1)(-1) is the streamwise wavenumber and U-tau is the friction velocity. Interestingly, this k(1)(-1) scaling is observed much farther from the wall and at much lower flow Reynolds number (both differing by almost an order of magnitude) than what the expectations from experiments on a zero-pressure-gradient turbulent boundary layer flow would suggest. Furthermore, the coefficient A(1) in the present sink-flow data is seen to be non-universal, i.e. A(1) varies with height from the wall; the scaling exponent -1 remains universal. Logarithmic variation of the so-called longitudinal structure function, which is the physical-space counterpart of spectral k(1)(-1) scaling, is also seen to be non-universal, consistent with the non-universality of A(1). These observations are to be contrasted with the universal value of A(1) (along with the universal scaling exponent of 1) reported in the literature on zero-pressure-gradient turbulent boundary layers. Theoretical arguments based on dimensional analysis indicate that the presence of a streamwise pressure gradient in sink-flow turbulent boundary layers makes the coefficient A(1) non-universal while leaving the scaling exponent -1 unaffected. This effect of the pressure gradient on the streamwise spectra, as discussed in the present study (experiments as well as theory), is consistent with other recent studies in the literature that are focused on the structural aspects of turbulent boundary layer flows in pressure gradients (Harun etal., J. Flui(d) Mech., vol. 715, 2013, pp. 477-498); the present paper establishes the link between these two. The variability of A(1) accommodated in the present framework serves to clarify the ideas of universality of the k(1)(-1) scaling.

Extended self-similar scaling law of multi-scale eddy structure in wall turbulence

The longitudinal fluctuating velocity of a turbulent boundary layer was measured in a water channel at a moderate Reynolds number. The extended self-similar scaling law of structure function proposed by Benzi was verified. The longitudinal fluctuating velocity, in the turbulent boundary layer was decomposed into many multi-scale eddy structures by wavelet transform. The extended self-similar scaling law of structure function for each scale eddy velocity was investigated. The conclusions are I) The statistical properties of turbulence could be self-similar not only at high Reynolds number, but also at moderate and low Reynolds number, and they could be characterized by the same set of scaling exponents xi (1)(n) = n/3 and xi (2)(n) = n/3 of the fully developed regime. 2) The range of scales where the extended self-similarity valid is much larger than the inertial range and extends far deep into the dissipation range,vith the same set of scaling exponents. 3) The extended selfsimilarity is applicable not only for homogeneous turbulence, but also for shear turbulence such as turbulent boundary layers.

Asymptotic scaling in the two-dimensional O(3) Nonlinear sigma model: a Monte Carlo study on parallel computers

We investigate the 2d O(3) model with the standard action by Monte Carlo simulation at couplings β up to 2.05. We measure the energy density, mass gap and susceptibility of the model, and gather high statistics on lattices of size L ≤ 1024 using the Floating Point Systems T-series vector hypercube and the Thinking Machines Corp.'s Connection Machine 2. Asymptotic scaling does not appear to set in for this action, even at β = 2.10, where the correlation length is 420. We observe a 20% difference between our estimate m/Λ^─_(Ms) = 3.52(6) at this β and the recent exact analytical result . We use the overrelaxation algorithm interleaved with Metropolis updates and show that decorrelation time scales with the correlation length and the number of overrelaxation steps per sweep. We determine its effective dynamical critical exponent to be z' = 1.079(10); thus critical slowing down is reduced significantly for this local algorithm that is vectorizable and parallelizable.

We also use the cluster Monte Carlo algorithms, which are non-local Monte Carlo update schemes which can greatly increase the efficiency of computer simulations of spin models. The major computational task in these algorithms is connected component labeling, to identify clusters of connected sites on a lattice. We have devised some new SIMD component labeling algorithms, and implemented them on the Connection Machine. We investigate their performance when applied to the cluster update of the two dimensional Ising spin model.

Finally we use a Monte Carlo Renormalization Group method to directly measure the couplings of block Hamiltonians at different blocking levels. For the usual averaging block transformation we confirm the renormalized trajectory (RT) observed by Okawa. For another improved probabilistic block transformation we find the RT, showing that it is much closer to the Standard Action. We then use this block transformation to obtain the discrete β-function of the model which we compare to the perturbative result. We do not see convergence, except when using a rescaled coupling β_E to effectively resum the series. For the latter case we see agreement for m/ Λ^─_(Ms) at , β = 2.14, 2.26, 2.38 and 2.50. To three loops m/Λ^─_(Ms) = 3.047(35) at β = 2.50, which is very close to the exact value m/ Λ^─_(Ms) = 2.943. Our last point at β = 2.62 disagrees with this estimate however.

Optical interconnection for neural networks by use of a self-imaging function

A grating-lens combination unit is developed to form a scaling self-transform function that can self-image on scale. Then an array of many such grating-lens units is used for the optical interconnection of a two-dimensional neural network, and experiments are carried out. We find that our idea is feasible, the optical interconnection system is simple, and optical adjustment is easy. (C) 1998 Optical Society of America.

Scaling the Indian Buffet Process via Submodular Maximization

Inference for latent feature models is inherently difficult as the inference space grows exponentially with the size of the input data and number of latent features. In this work, we use Kurihara & Welling (2008)'s maximization-expectation framework to perform approximate MAP inference for linear-Gaussian latent feature models with an Indian Buffet Process (IBP) prior. This formulation yields a submodular function of the features that corresponds to a lower bound on the model evidence. By adding a constant to this function, we obtain a nonnegative submodular function that can be maximized via a greedy algorithm that obtains at least a one-third approximation to the optimal solution. Our inference method scales linearly with the size of the input data, and we show the efficacy of our method on the largest datasets currently analyzed using an IBP model.

Scaling arguments for the fluxes in turbulent miscible fountains

For established axisymmetric turbulent miscible Boussinesq fountains in quiescent uniform environments, expressions are developed for the fluxes of volume, momentum and buoyancy at the outflow from the fountain: the outflow referring to the counterflow at the horizontal plane of the source. The fluxes are expressed in terms of the fountain source conditions and two dimensionless functions of the source Froude number, Fr0: a radial function (relating a horizontal scale of the outflow to the source radius) and a volume flux function (relating the outflow and source volume fluxes). The forms taken by these two functions at low Fr0 and high Fr0 are deduced, thereby providing the outflow fluxes and outflow Froude number solely in terms of the source conditions. For high Fr0, the outflow Froude number, Frout, is shown to be invariant, indicating (by analogy with plumes for which the 'far-field' Froude number is invariant with source Froude number) that the outflow may be regarded as 'far-field' since the fluxes within the fountain have adjusted to attain a balance which is independent of the source conditions. Based on Frout, the fluxes in the plume that forms beyond the fountain outflow are deduced. Finally, from the results of previously published studies, we show that the scalings deduced for fountains are valid for 0.0025 ≲ Fr0 ≲ 1.0 for low Fr0 and Fr0≳ 3.0 for high Fr0. © 2014 Cambridge University Press.

Scaling and the thermal conductivity of the Frenkel-Kontorova model

We have studied the dependence of the thermal conductivity kappa on the strength of the interparticle potential lambda and the strength of the external potential beta in the Frenkel-Kontorova model. We found that the functional relation can be expressed in a scaling form, kappa(proportional to) lambda 3/2/beta(2 center dot). This result is first obtained by nonequilibrium molecular dynamics. It is then confirmed by two analytical methods, the self-consistent phonon theory and the self-consistent stochastic reservoirs method. The thermal conductivity kappa is therefore a decreasing functon of beta and an increasing function of lambda.

DYNAMIC SCALING AND FRACTAL BEHAVIOR OF SPINODAL PHASE-SEPARATION IN POLYMER MIXTURES OF POLY(METHYL METHACRYLATE) WITH POLY(STYRENE-CO-ACRYLONITRILE)

Dynamic scaling and fractal behaviour of spinodal phase separation is studied in a binary polymer mixture of poly(methyl methacrylate) (PMMA) and poly(styrene-co-acrylonitrile) (SAN). In the later stages of spinodal phase separation, a simple dynamic scaling law was found for the scattering function S(q,t):S(q,t) approximately q(m)-3S approximately (q/q(m)). The possibility of using fractal theory to describe the complex morphology of spinodal phase separation is discussed. In phase separation, morphology exhibits strong self-similarity. The two-dimensional image obtained by optical microscopy can be analysed within the framework of fractal concepts. The results give a fractal dimension of 1.64. This implies that the fractal structure may be the reason for the dynamic scaling behaviour of the structure function.

Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.

Single-ionization of helium at Ti:sapphire wavelengths: rates and scaling laws

We present a numerical and theoretical study of intense-field single-electron ionization of helium at 390 nm and 780 nm. Accurate ionization rates (over an intensity range of (0.175-34) X10^14 W/ cm^2 at 390 nm, and (0.275 - 14.4) X 10^14 W /cm^2 at 780 nm) are obtained from full-dimensionality integrations of the time-dependent helium-laser Schroedinger equation. We show that the power law of lowest order perturbation theory, modified with a ponderomotive-shifted ionization potential, is capable of modelling the ionization rates over an intensity range that extends up to two orders of magnitude higher than that applicable to perturbation theory alone. Writing the modified perturbation theory in terms of scaled wavelength and intensity variables, we obtain to first approximation a single ionization law for both the 390 nm and 780 nm cases. To model the data in the high intensity limit as well as in the low, a new function is introduced for the rate. This function has, in part, a resemblance to that derived from tunnelling theory but, importantly, retains the correct frequency-dependence and scaling behaviour derived from the perturbative-like models at lower intensities. Comparison with the predictions of classical ADK tunnelling theory confirms that ADK performs poorly in the frequency and intensity domain treated here.

Wall thickness dependence of the scaling law for ferroic stripe domains

The periodicity of 180 degrees. stripe domains as a function of crystal thickness scales with the width of the domain walls, both for ferroelectric and for ferromagnetic materials. Here we derive an analytical expression for the generalized ferroic scaling factor and use this to calculate the domain wall thickness and gradient coefficients ( exchange constants) in some ferroelectric and ferromagnetic materials. We then use these to discuss some of the wider implications for the physics of ferroelectric nanodevices and periodically poled photonic crystals.

The top-down mechanism for body-mass-abundance scaling

Scaling relationships between mean body masses and abundances of species in multitrophic communities continue to be a subject of intense research and debate. The top-down mechanism explored in this paper explains the frequently observed inverse linear relationship between body mass and abundance (i.e., constant biomass) in terms of a balancing of resource biomasses by behaviorally and evolutionarily adapting foragers, and the evolutionary response of resources to this foraging pressure. The mechanism is tested using an allometric, multitrophic community model with a complex food web structure. It is a statistical model describing the evolutionary and population dynamics of tens to hundreds of species in a uniform way. Particularities of the model are the detailed representation of the evolution and interaction of trophic traits to reproduce topological food web patterns, prey switching behavior modeled after experimental observations, and the evolutionary adaptation of attack rates. Model structure and design are discussed. For model states comparable to natural communities, we find that (1) the body-mass-abundance scaling does not depend on the allometric scaling exponent of physiological rates in the form expected from the energetic equivalence rule or other bottom-up theories; (2) the scaling exponent of abundance as a function of body mass is approximately -1, independent of the allometric exponent for physiological rates assumed; (3) removal of top-down control destroys this pattern, and energetic equivalence is recovered. We conclude that the top-down mechanism is active in the model, and that it is a viable alternative to bottom-up mechanisms for controlling body-mass-abundance relations in natural communities.

Scaling of Food-Web Properties with Diversity and Complexity Across Ecosystems

Trophic scaling models describe how topological food-web properties such as the number of predator prey links scale with species richness of the community. Early models predicted that either the link density (i.e. the number of links per species) or the connectance (i.e. the linkage probability between any pair of species) is constant across communities. More recent analyses, however, suggest that both these scaling models have to be rejected, and we discuss several hypotheses that aim to explain the scale dependence of these complexity parameters. Based on a recent, highly resolved food-web compilation, we analysed the scaling behaviour of 16 topological parameters and found significant power law scaling relationships with diversity (i.e. species richness) and complexity (i.e. connectance) for most of them. These results illustrate the lack of universal constants in food-web ecology as a function of diversity or complexity. Nonetheless, our power law scaling relationships suggest that fundamental processes determine food-web topology, and subsequent analyses demonstrated that ecosystem-specific differences in these relationships were of minor importance. As such, these newly described scaling relationships provide robust and testable cornerstones for future structural food-web models.

Testing metabolic scaling theory using intraspecific allometries in Antarctic microarthropods

Quantitative scaling relationships among body mass, temperature and metabolic rate of organisms are still controversial, while resolution may be further complicated through the use of different and possibly inappropriate approaches to statistical analysis. We propose the application of a modelling strategy based on the theoretical approach of Akaike's information criteria and non-linear model fitting (nlm). Accordingly, we collated and modelled available data at intraspecific level on the individual standard metabolic rate of Antarctic microarthropods as a function of body mass (M), temperature (T), species identity (S) and high rank taxa to which species belong (G) and tested predictions from metabolic scaling theory (mass-metabolism allometric exponent b = 0.75, activation energy range 0.2-1.2 eV). We also performed allometric analysis based on logarithmic transformations (lm). Conclusions from lm and nlm approaches were different. Best-supported models from lm incorporated T, M and S. The estimates of the allometric scaling exponent linking body mass and metabolic rate resulted in a value of 0.696 +/- 0.105 (mean +/- 95% CI). In contrast, the four best-supported nlm models suggested that both the scaling exponent and activation energy significantly vary across the high rank taxa (Collembola, Cryptostigmata, Mesostigmata and Prostigmata) to which species belong, with mean values of b ranging from about 0.6 to 0.8. We therefore reached two conclusions: 1, published analyses of arthropod metabolism based on logarithmic data may be biased by data transformation; 2, non-linear models applied to Antarctic microarthropod metabolic rate suggest that intraspecific scaling of standard metabolic rate in Antarctic microarthropods is highly variable and can be characterised by scaling exponents that greatly vary within taxa, which may have biased previous interspecific comparisons that neglected intraspecific variability.