Some recent advances in the theory of homogeneous isotropic turbulence

We review some advances in the theory of homogeneous, isotropic turbulence. Our emphasis is on the new insights that have been gained from recent numerical studies of the three-dimensional Navier Stokes equation and simpler shell models for turbulence. In particular, we examine the status of multiscaling corrections to Kolmogorov scaling, extended self similarity, generalized extended self similarity, and non-Gaussian probability distributions for velocity differences and related quantities. We recount our recent proposal of a wave-vector-space version of generalized extended self similarity and show how it allows us to explore an intriguing and apparently universal crossover from inertial- to dissipation-range asymptotics.

Thermal sensitivity analysis of spacecraft battery

Transient thermal sensitivity is studied for systems that are subjected to conductive heat transfer within themselves and radiative heat transfer with the surrounding environment, including solar heat radiation, The battery in the Indian national communication satellite is one such system for which the studies are conducted with respect to panel conduction, conductance of insulating blanket, power dissipation within the battery, and absorptance and emittance of various elements, Comparison of sensitivities revealed that battery temperature is sensitive to its power dissipation during the beginning of life of the spacecraft, whereas toward the end of life of the spacecraft mission, the effect of absorptance of optical solar reflector is dominating, The influence of optical property values of the multilayer insulation blanket is almost negligible. Among the parameters studied in this analysis, the battery temperature is found to be mast sensitive to emittance of the optical solar reflector.

Antiferroelectric lead zirconate thin films by pulsed laser ablation

Lead Zirconate (PbZrO3) thin films were deposited by pulsed laser ablation method. Pseudocubic (110) oriented in-situ films were grown at low pressure. The field enforced anti-ferroelectric (AFE) to ferroelectric (FE) phase transformation behaviour was investigated by means of a modified Sawyer Tower circuit as well as capacitance versus applied voltage measurements. The maximum polarisation obtained was 36 mu C cm(-2) and the critical field to induce ferroelectric state and to reverse the antiferroelectric slates were 65 and 90 kV cm(-1) respectively. The dielectric properties were investigated as a function of frequency and temperature. The dielectric constant of the AFE lead zirconate thin him was 190 at 100 kHz which is more than the bulk ceramic value (120) with a dissipation factor of less than 0.07. The polarisation switching kinetics of the antiferroelectric PbZrO3 thin films showed that the switching time to be around 275 ns between antipolar state to polar states. (C) 1999 Elsevier Science S.A. All rights reserved.

Low-power pipelined LMS adaptive filter architectures with minimal adaptation delay

The use of delayed coefficient adaptation in the least mean square (LMS) algorithm has enabled the design of pipelined architectures for real-time transversal adaptive filtering. However, the convergence speed of this delayed LMS (DLMS) algorithm, when compared with that of the standard LMS algorithm, is degraded and worsens with increase in the adaptation delay. Existing pipelined DLMS architectures have large adaptation delay and hence degraded convergence speed. We in this paper, first present a pipelined DLMS architecture with minimal adaptation delay for any given sampling rate. The architecture is synthesized by using a number of function preserving transformations on the signal flow graph representation of the DLMS algorithm. With the use of carry-save arithmetic, the pipelined architecture can support high sampling rates, limited only by the delay of a full adder and a 2-to-1 multiplexer. In the second part of this paper, we extend the synthesis methodology described in the first part, to synthesize pipelined DLMS architectures whose power dissipation meets a specified budget. This low-power architecture exploits the parallelism in the DLMS algorithm to meet the required computational throughput. The architecture exhibits a novel tradeoff between algorithmic performance (convergence speed) and power dissipation. (C) 1999 Elsevier Science B.V. All rights resented.

Irreversibility-to-reversibility crossover in transient response of an optically trapped particle

We study the transient response of a colloidal bead which is released from different heights and allowed to relax in the potential well of an optical trap. Depending on the initial potential energy, the system's time evolution shows dramatically different behaviors. Starting from the short-time reversible to long-time irreversible transition, a stationary reversible state with zero net dissipation can be achieved as the release point energy is decreased. If the system starts with even lower energy, it progressively extracts useful work from thermal noise and exhibits an anomalous irreversibility. In addition, we have verified the Transient Fluctuation Theorem and the Integrated Transient Fluctuation Theorem even for the non-ergodic descriptions of our system. Copyright (C) EPLA, 2011

Nonequilibrium stationary states of a particle in a gravitational field driven by a vibrating surface

Stationary velocity distribution functions are determined for a particle in a gravitational field driven by a vibrating surface in the limit of small dissipation. It is found that the form of the distribution function is sensitive to the mechanism of energy dissipation, inelastic collisions or viscous drag, and also to the form of the amplitude function of the vibrating surface. The velocity distributions obtained analytically are found to be in excellent agreement with the results of computer simulations in the limit of low dissipation. [S0031-9007(99)08898-5].

Study of the electrical properties of pulsed laser ablated (Ba0.5Sr0.5)TiO3 thin films

The laser ablated barium strontium titanate (BST) thin films were characterized in terms of composition, structure, microstructure and electrical properties. Films deposited at 300 degrees C under 50 mTorr oxygen pressure and 3 J cm(-2) laser fluence and further annealed at 600 degrees C in flowing oxygen showed a dielectric constant of 467 and a dissipation factor of 0.02. The room-temperature current-voltage characteristics revealed a space charge limited conduction (SCLC) mechanism, though at low fields the effect of the electrodes was predominant. The conduction mechanism was thoroughly-investigated in terms of Schottky emission at low fields, and bulk-limited SCLC at high fields. The change over to the bulk-limited conduction process from the electrode-limited Schottky emission was, attributed to the process of tunneling through the electrode interface at high fields resulting into the lowering of the electrode contact resistance and consequently giving rise to a bulk limited conduction process. The predominance of SCLC mechanism in the films suggests that the bulk properties are only revealed if the depletion width at the electrode interface is thin enough to allow the tunneling process to take place. This condition is only favorable if the him thickness is high or if the doping concentration is high enough. In the present case the film thickness ranged from 0.3 to 0.7 mu m which was suitable to show the transition mentioned above. (C) 1999 Elsevier Science S.A. All rights reserved.

Temperature scaling in a dense vibrofluidized granular material

The leading order "temperature" of a dense two-dimensional granular material fluidized by external vibrations is determined. The grain interactions are characterized by inelastic collisions, but the coefficient of restitution is considered to be close to 1, so that the dissipation of energy during a collision is small compared to the average energy of a particle. An asymptotic solution is obtained where the particles are considered to be elastic in the leading approximation. The velocity distribution is a Maxwell-Boltzmann distribution in the leading approximation,. The density profile is determined by solving the momentum balance equation in the vertical direction, where the relation between the pressure and density is provided by the virial equation of state. The temperature is determined by relating the source of energy due to the vibrating surface and the energy dissipation due to inelastic collisions. The predictions of the present analysis show good agreement with simulation results at higher densities where theories for a dilute vibrated granular material, with the pressure-density relation provided by the ideal gas law, sire in error. [:S1063-651X(99)04408-6].

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

Assessment of conditional moment closure models of turbulent autoignition using DNS data

A systematic assessment of the submodels of conditional moment closure (CMC) formalism for the autoignition problem is carried out using direct numerical simulation (DNS) data. An initially non-premixed, n-heptane/air system, subjected to a three-dimensional, homogeneous, isotropic, and decaying turbulence, is considered. Two kinetic schemes, (1) a one-step and (2) a reduced four-step reaction mechanism, are considered for chemistry An alternative formulation is developed for closure of the mean chemical source term , based on the condition that the instantaneous fluctuation of excess temperature is small. With this model, it is shown that the CMC equations describe the autoignition process all the way up to near the equilibrium limit. The effect of second-order terms (namely, conditional variance of temperature excess sigma(2) and conditional correlations of species q(ij)) in modeling is examined. Comparison with DNS data shows that sigma(2) has little effect on the predicted conditional mean temperature evolution, if the average conditional scalar dissipation rate is properly modeled. Using DNS data, a correction factor is introduced in the modeling of nonlinear terms to include the effect of species fluctuations. Computations including such a correction factor show that the species conditional correlations q(ij) have little effect on model predictions with a one-step reaction, but those q(ij) involving intermediate species are found to be crucial when four-step reduced kinetics is considered. The "most reactive mixture fraction" is found to vary with time when a four-step kinetics is considered. First-order CMC results are found to be qualitatively wrong if the conditional mean scalar dissipation rate is not modeled properly. The autoignition delay time predicted by the CMC model compares excellently with DNS results and shows a trend similar to experimental data over a range of initial temperatures.

Autoignition in a non-premixed medium: DNS studies on the effects of three-dimensional turbulence

Direct numerical simulation (DNS) results of autoignition in anon-premixed medium under an isotropic, homogeneous, and decaying turbulence are presented. The initial mixture consists of segregated fuel parcels randomly distributed within warm air, and the entire medium is subjected to a three-dimensional turbulence. Chemical kinetics is modeled by a four-step reduced reaction mechanism for autoignition of n-heptane/air mixture. Thus, this work overcomes the principal limitations of a previous contribution of the authors on two-dimensional DNS of autoignition with a one-step reaction model. Specific attention is focused on the differences in the effects of two- and three-dimensional turbulence on autoignition characteristics. The three-dimensional results show that ignition spots are most likely to originate at locations jointly corresponding to the most reactive mixture fraction and low scalar dissipation rate. Further, these ignition spots are found to originate at locations corresponding to the core of local vortical structures, and after ignition, the burning gases move toward the vortex periphery Such a movement is explained as caused by the cyclostrophic imbalance developed when the local gas density is variable. These results lead to the conclusion that the local ignition-zone structure does not conform to the classical stretched flamelet description. Parametric studies show that the ignition delay time decreases with an increase in turbulence intensity. Hence, these three-dimensional simulation results resolve the discrepancy between trends in experimental data and predictions from DNSs of two-dimensional turbulence. This qualitative difference between DNS results from three- and two-dimensional simulations is discussed and attributed to the effect of vortex stretching that is present in the former, but not in the latter.

Phase diagram of a two-species lattice model with a linear instability

We discuss the properties of a one-dimensional lattice model of a driven system with two species of particles in which the mobility of one species depends on the density of the other. This model was introduced by Lahiri and Ramaswamy (Phys. Rev. Lett., 79, 1150 (1997)) in the context of sedimenting colloidal crystals, and its continuum version was shown to exhibit an instability arising from linear gradient couplings. In this paper we review recent progress in understanding the full phase diagram of the model. There are three phases. In the first, the steady state can be determined exactly along a representative locus using the condition of detailed balance. The system shows phase separation of an exceptionally robust sort, termed strong phase separation, which survives at all temperatures. The second phase arises in the threshold case where the first species evolves independently of the second, but the fluctuations of the first influence the evolution of the second, as in the passive scalar problem. The second species then shows phase separation of a delicate sort, in which long-range order coexists with fluctuations which do not damp down in the large-size limit. This fluctuation-dominated phase ordering is associated with power law decays in cluster size distributions and a breakdown of the Porod law. The third phase is one with a uniform overall density, and along a representative locus the steady state is shown to have product measure form. Density fluctuations are transported by two kinematic waves, each involving both species and coupled at the nonlinear level. Their dissipation properties are governed by the symmetries of these couplings, which depend on the overall densities. In the most interesting case,, the dissipation of the two modes is characterized by different critical exponents, despite the nonlinear coupling.

Direct numerical simulation of turbulent spots

A group of high-order finite-difference schemes for incompressible flow was implemented to simulate the evolution of turbulent spots in channel flows. The long-time accuracy of these schemes was tested by comparing the evolution of small disturbances to a plane channel flow against the growth rate predicted by linear theory. When the perturbation is the unstable eigenfunction at a Reynolds number of 7500, the solution grows only if there are a comparatively large number of (equispaced) grid points across the channel. Fifth-order upwind biasing of convection terms is found to be worse than second-order central differencing. But, for a decaying mode at a Reynolds number of 1000, about a fourth of the points suffice to obtain the correct decay rate. We show that this is due to the comparatively high gradients in the unstable eigenfunction near the walls. So, high-wave-number dissipation of the high-order upwind biasing degrades the solution especially. But for a well-resolved calculation, the weak dissipation does not degrade solutions even over the very long times (O(100)) computed in these tests. Some new solutions of spot evolution in Couette flows with pressure gradients are presented. The approach to self-similarity at long times can be seen readily in contour plots.

Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil. (C) 2011 Elsevier Ltd. All rights reserved.

Turbulence measurements in an ‘equilibrium’ axisymmetric wall jet Journal

This paper reports measurements of turbulent quantities in an axisymmetric wall jet subjected to an adverse pressure gradient in a conical diffuser, in such a way that a suitably defined pressure-gradient parameter is everywhere small. Self-similarity is observed in the mean velocity profile, as well as the profiles of many turbulent quantities at sufficiently large distances from the injection slot. Autocorrelation measurements indicate that, in the region of turbulent production, the time scale of ν fluctuations is very much smaller than the time scale of u fluctuations. Based on the data on these time scales, a possible model is proposed for the Reynolds stress. One-dimensional energy spectra are obtained for the u, v and w components at several points in the wall jet. It is found that self-similarity is exhibited by the one-dimensional wavenumber spectrum of $\overline{q^2}(=\overline{u^2}+\overline{v^2}+\overline{w^2})$, if the half-width of the wall jet and the local mean velocity are used for forming the non-dimensional wavenumber. Both the autocorrelation curves and the spectra indicate the existence of periodicity in the flow. The rate of dissipation of turbulent energy is estimated from the $\overline{q^2}$ spectra, using a slightly modified version of a previously suggested method.