Spectral analysis of structure functions and their scaling exponents in forced isotropic turbulence

The pseudospectral method, in conjunction with a new technique for obtaining scaling exponents $\zeta_n$ from the structure functions $S_n(r)$, is presented as an alternative to the extended self-similarity (ESS) method and the use of generalized structure functions. We propose plotting the ratio $|S_n(r)/S_3(r)|$ against the separation $r$ in accordance with a standard technique for analysing experimental data. This method differs from the ESS technique, which plots $S_n(r)$ against $S_3(r)$, with the assumption $S_3(r) \sim r$. Using our method for the particular case of $S_2(r)$ we obtain the new result that the exponent $\zeta_2$ decreases as the Taylor-Reynolds number increases, with $\zeta_2 \to 0.679 \pm 0.013$ as $R_{\lambda} \to \infty$. This supports the idea of finite-viscosity corrections to the K41 prediction for $S_2$, and is the opposite of the result obtained by ESS. The pseudospectral method also permits the forcing to be taken into account exactly through the calculation of the energy input in real space from the work spectrum of the stirring forces.Comment: 31 pages including appendices, 10 figure

Energy transfer and dissipation in forced isotropic turbulence

A model for the Reynolds number dependence of the dimensionless dissipation rate $C_{\varepsilon}$ was derived from the dimensionless K\'{a}rm\'{a}n-Howarth equation, resulting in $C_{\varepsilon}=C_{\varepsilon, \infty} + C/R_L + O(1/R_L^2)$, where $R_L$ is the integral scale Reynolds number. The coefficients $C$ and $C_{\varepsilon,\infty}$ arise from asymptotic expansions of the dimensionless second- and third-order structure functions. This theoretical work was supplemented by direct numerical simulations (DNSs) of forced isotropic turbulence for integral scale Reynolds numbers up to $R_L=5875$ ($R_\lambda=435$), which were used to establish that the decay of dimensionless dissipation with increasing Reynolds number took the form of a power law $R_L^n$ with exponent value $n = -1.000\pm 0.009$, and that this decay of $C_{\varepsilon}$ was actually due to the increase in the Taylor surrogate $U^3/L$. The model equation was fitted to data from the DNS which resulted in the value $C=18.9\pm 1.3$ and in an asymptotic value for $C_\varepsilon$ in the infinite Reynolds number limit of $C_{\varepsilon,\infty} = 0.468 \pm 0.006$.Comment: 26 pages including references and 6 figures. arXiv admin note: text overlap with arXiv:1307.457

Re-examination of the infra-red properties of randomly stirred hydrodynamics

Dynamic renormalization group (RG) methods were originally used by Forster, Nelson and Stephen (FNS) to study the large-scale behaviour of randomly-stirred, incompressible fluids governed by the Navier-Stokes equations. Similar calculations using a variety of methods have been performed since, but have led to a discrepancy in results. In this paper, we carefully re-examine in $d$-dimensions the approaches used to calculate the renormalized viscosity increment and, by including an additional constraint which is neglected in many procedures, conclude that the original result of FNS is correct. By explicitly using step functions to control the domain of integration, we calculate a non-zero correction caused by boundary terms which cannot be ignored. We then go on to analyze how the noise renormalization, absent in many approaches, contributes an ${\mathcal O}(k^2)$ correction to the force autocorrelation and show conditions for this to be taken as a renormalization of the noise coefficient. Following this, we discuss the applicability of this RG procedure to the calculation of the inertial range properties of fluid turbulence.Comment: 16 pages, 6 figure

A formal derivation of the local energy transfer (LET) theory of homogeneous turbulence

A statistical closure of the Navier-Stokes hierarchy which leads to equations for the two-point, two-time covariance of the velocity field for stationary, homogeneous isotropic turbulence is presented. It is a generalisation of the self-consistent field method due to Edwards (1964) for the stationary, single-time velocity covariance. The probability distribution functional P [u, t] is obtained, in the form of a series, from the Liouville equation by means of a perturbation expansion about a Gaussian distribution, which is chosen to give the exact two-point, two-time covariance. The triple moment is calculated in terms of an ensemble-averaged infinitesimal velocity-field propagator, and shown to yield the Edwards result as a special case. The use of a Gaussian zero-order distribution has been found to justify the introduction of a fluctuation-response relation, which is in accord with modern dynamical theories. In a sense this work completes the analogy drawn by Edwards between turbulence and Brownian motion. Originally Edwards had shown that the noise input was determined by the correlation of the velocity field with the externally applied stirring forces but was unable to determine the system response. Now we find that the system response is determined by the correlation of the velocity field with internal quasi-entropic forces. This analysis is valid to all orders of perturbation theory, and allows the recovery of the Local Energy Transfer (LET) theory, which had previously been derived by more heuristical methods. The LET theory is known to be in good agreement with experimental results. It is also unique among two-point statistical closures in displaying an acceptable (i.e. non-Markovian) relationship between the transfer spectrum and the system response, in accordance with experimental results. As a result of the latter property, it is compatible with the Kolmogorov (K41) spectral phenomenology

Raman backscattering saturation due to coupling between ωp and 2ωp modes in plasma

Raman backscattering (RBS) in plasma is the basis of plasma-based amplifiers and is important in laser-driven fusion experiments. We show that saturation can arise from nonlinearities due to coupling between the fundamental and harmonic plasma wave modes for sufficiently intense pump and seed pulses. We present a time-dependent analysis that shows that plasma wave phase shifts reach a maximum close to wavebreaking. The study contributes to a new understanding of RBS saturation for counter-propagating laser pulses

Distribution of graph-distances in Boltzmann ensembles of RNA secondary structures

Large RNA molecules often carry multiple functional domains whose spatial arrangement is an important determinant of their function. Pre-mRNA splicing, furthermore, relies on the spatial proximity of the splice junctions that can be separated by very long introns. Similar effects appear in the processing of RNA virus genomes. Albeit a crude measure, the distribution of spatial distances in thermodynamic equilibrium therefore provides useful information on the overall shape of the molecule can provide insights into the interplay of its functional domains. Spatial distance can be approximated by the graph-distance in RNA secondary structure. We show here that the equilibrium distribution of graph-distances between arbitrary nucleotides can be computed in polynomial time by means of dynamic programming. A naive implementation would yield recursions with a very high time complexity of O(n^11). Although we were able to reduce this to O(n^6) for many practical applications a further reduction seems difficult. We conclude, therefore, that sampling approaches, which are much easier to implement, are also theoretically favorable for most real-life applications, in particular since these primarily concern long-range interactions in very large RNA molecules.Comment: Peer-reviewed and presented as part of the 13th Workshop on Algorithms in Bioinformatics (WABI2013

Nonequilibrium brittle fracture propagation: Steady state, oscillations and intermittency

A minimal model is constructed for two-dimensional fracture propagation. The heterogeneous process zone is presumed to suppress stress relaxation rate, leading to non-quasistatic behavior. Using the Yoffe solution, I construct and solve a dynamical equation for the tip stress. I discuss a generic tip velocity response to local stress and find that noise-free propagation is either at steady state or oscillatory, depending only on one material parameter. Noise gives rise to intermittency and quasi-periodicity. The theory explains the velocity oscillations and the complicated behavior seen in polymeric and amorphous brittle materials. I suggest experimental verifications and new connections between velocity measurements and material properties.Comment: To appear in Phys. Rev. Lett., 6 pages, self-contained TeX file, 3 postscript figures upon request from author at [email protected] or [email protected], http://cnls-www.lanl.gov/homepages/rafi/rafindex.htm

An accurate description of quantum size effects in InP nanocrystallites over a wide range of sizes

We obtain an effective parametrization of the bulk electronic structure of InP within the Tight Binding scheme. Using these parameters, we calculate the electronic structure of InP clusters with the size ranging upto 7.5 nm. The calculated variations in the electronic structure as a function of the cluster size is found to be in excellent agreement with experimental results over the entire range of sizes, establishing the effectiveness and transferability of the obtained parameter strengths.Comment: 9 pages, 3 figures, pdf file available at http://sscu.iisc.ernet.in/~sampan/publications.htm

Oscillating Fracture in Rubber

We have found an oscillating instability of fast-running cracks in thin rubber sheets. A well-defined transition from straight to oscillating cracks occurs as the amount of biaxial strain increases. Measurements of the amplitude and wavelength of the oscillation near the onset of this instability indicate that the instability is a Hopf bifurcation

Towards attosecond high-energy electron bunches : controlling self-injection in laser wakefield accelerators through plasma density modulation

Self-injection in a laser-plasma wakefield accelerator (LWFA) is usually achieved by increasing the laser intensity until the threshold for injection is exceeded. Alternatively, the velocity of the bubble accelerating structure can be controlled using plasma density ramps, reducing the electron velocity required for injection. We present a model describing self-injection in the short bunch regime for arbitrary changes in the plasma density. We derive the threshold condition for injection due to a plasma density gradient, which is confirmed using particle-in-cell (PIC) simulations that demonstrate injection of sub-femtosecond bunches. It is shown that the bunch charge, bunch length and separation of bunches in a bunch train can be controlled by tailoring the plasma density profile