Averting biodiversity collapse in tropical forest protected areas

The rapid disruption of tropical forests probably imperils global biodiversity more than any other contemporary phenomenon(1-3). With deforestation advancing quickly, protected areas are increasingly becoming final refuges for threatened species and natural ecosystem processes. However, many protected areas in the tropics are themselves vulnerable to human encroachment and other environmental stresses(4-9). As pressures mount, it is vital to know whether existing reserves can sustain their biodiversity. A critical constraint in addressing this question has been that data describing a broad array of biodiversity groups have been unavailable for a sufficiently large and representative sample of reserves. Here we present a uniquely comprehensive data set on changes over the past 20 to 30 years in 31 functional groups of species and 21 potential drivers of environmental change, for 60 protected areas stratified across the world's major tropical regions. Our analysis reveals great variation in reserve `health': about half of all reserves have been effective or performed passably, but the rest are experiencing an erosion of biodiversity that is often alarmingly widespread taxonomically and functionally. Habitat disruption, hunting and forest-product exploitation were the strongest predictors of declining reserve health. Crucially, environmental changes immediately outside reserves seemed nearly as important as those inside in determining their ecological fate, with changes inside reserves strongly mirroring those occurring around them. These findings suggest that tropical protected areas are often intimately linked ecologically to their surrounding habitats, and that a failure to stem broad-scale loss and degradation of such habitats could sharply increase the likelihood of serious biodiversity declines.

Dynamical configurations and bistability of helical nanostructures under external torque

We study the motion of a ferromagnetic helical nanostructure under the action of a rotating magnetic field. A variety of dynamical configurations were observed that depended strongly on the direction of magnetization and the geometrical parameters, which were also confirmed by a theoretical model, based on the dynamics of a rigid body under Stokes flow. Although motion at low Reynolds numbers is typically deterministic, under certain experimental conditions the nanostructures showed a surprising bistable behavior, such that the dynamics switched randomly between two configurations, possibly induced by thermal fluctuations. The experimental observations and the theoretical results presented in this paper are general enough to be applicable to any system of ellipsoidal symmetry under external force or torque.

Cumulus clouds and convective boundary layers: a tropical perspective on two turbulent shear flows

This paper is a review prepared for the second Marseille Colloquium on the mechanics of turbulence, held in 2011, 50 years after the first. The review covers recent developments in our understanding of the large-scale dynamics of cumulus cloud flows and of the atmospheric boundary layer in the low-wind convective regime that is often encountered in the tropics. It has recently been shown that a variety of cumulus cloud forms and life cycles can be experimentally realized in the laboratory, with the transient diabatic plume taken as the flow model for a cumulus cloud. The plume is subjected to diabatic heating scaled to be dynamically similar to heat release from phase changes in clouds. The experiments are complemented by exact numerical solutions of the Navier-Stokes-Boussinesq equations for plumes with scaled off-source heating. The results show that the Taylor entrainment coefficient first increases with heating, reaches a positive maximum and then drops rapidly to zero or even negative values. This reduction in entrainment is a consequence of structural changes in the flow, smoothing out the convoluted boundaries in the non-diabatic plume, including the tongues engulfing the ambient flow. This is accompanied by a greater degree of mixedness in the core flow because of lower dilution by the ambient fluid. The cloud forms generated depend strongly on the history of the diabatic heating profile in the vertical direction. The striking effects of heating on the flow are attributable to the operation of the baroclinic torque due to the temperature field. The mean baroclinic torque is shown to peak around a quasi-cylindrical sheet situated midway between the axis of the flow and the edges. This torque is shear-enhancing and folds down the engulfment tongues. The increase in mixedness can be traced to an explosive growth in the enstrophy, triggered by a strong fluctuating baroclinic torque that acts as a source, especially at the higher wave numbers, thus enhancing the mixedness. In convective boundary layers field measurements show that, under conditions prevailing in the tropics, the eddy fluxes of momentum and energy do not follow the Monin-Obukhov similarity. Instead, the eddy momentum flux is found to be linear in the wind speed at low winds; and the eddy heat flux is, to a first approximation, governed by free convection laws, with wind acting as a small perturbation on a regime of free convection. A new boundary layer code, based on heat flux scaling rather than wall-stress scaling, shows promising improvements in predictive skills of a general circulation model.

Stability of developing film flow down an inclined surface

Film flows on inclined surfaces are often assumed to be of constant thickness, which ensures that the velocity profile is half-Poiseuille. It is shown here that by shallow water theory, only flows in a portion of Reynolds number-Froude number (Re-Fr) plane can asymptotically attain constant film thickness. In another portion on the plane, the constant thickness solution appears as an unstable fixed point, while in other regions the film thickness seems to asymptote to a positive slope. Our simulations of the Navier-Stokes equations confirm the predictions of shallow water theory at higher Froude numbers, but disagree with them at lower Froude numbers. We show that different regimes of film flow show completely different stability behaviour from that predicted earlier. Supercritical decelerating flows are shown to be always unstable, whereas accelerating flows become unstable below a certain Reynolds number for a given Froude number. Subcritical flows on the other hand are shown to be unstable above a certain Reynolds number. In some range of parameters, two solutions for the base flowexist, and the attached profile is found to be more stable. All flows except those with separation become more stable as they proceed downstream. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4758299]

Chapman-Enskog analysis for finite-volume formulation of lattice Boltzmann equation

The classical Chapman-Enskog expansion is performed for the recently proposed finite-volume formulation of lattice Boltzmann equation (LBE) method D.V. Patil, K.N. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, J. Comput. Phys. 228 (2009) 5262-5279]. First, a modified partial differential equation is derived from a numerical approximation of the discrete Boltzmann equation. Then, the multi-scale, small parameter expansion is followed to recover the continuity and the Navier-Stokes (NS) equations with additional error terms. The expression for apparent value of the kinematic viscosity is derived for finite-volume formulation under certain assumptions. The attenuation of a shear wave, Taylor-Green vortex flow and driven channel flow are studied to analyze the apparent viscosity relation.

Analytical theory and stability analysis of an elongated nanoscale object under external torque

We consider the rotational motion of an elongated nanoscale object in a fluid under an external torque. The experimentally observed dynamics could be understood from analytical solutions of the Stokes equation, with explicit formulae derived for the dynamical states as a function of the object dimensions and the parameters defining the external torque. Under certain conditions, multiple analytical solutions to the Stokes equations exist, which have been investigated through numerical analysis of their stability against small perturbations and their sensitivity towards initial conditions. These experimental results and analytical formulae are general enough to be applicable to the rotational motion of any isolated elongated object at low Reynolds numbers, and could be useful in the design of non-spherical nanostructures for diverse applications pertaining to microfluidics and nanoscale propulsion technologies.

Loss and gain signals in broadband stimulated-Raman spectra: Theoretical analysis

Stimulated optical signals obtained by subjecting the system to a narrow band and a broadband pulse show both gain and loss Raman features at the red and blue side of the narrow beam, respectively. Recently observed temperature-dependent asymmetry in these features Mallick et al., J. Raman Spectrosc. 42, 1883 (2011); Dang et al., Phys. Rev. Lett. 107, 043001 (2011)] has been attributed to the Stokes and anti-Stokes components of the third-order susceptibility, chi((3)). By treating the setup as a steady state of an open system coupled to four quantum radiation field modes, we show that Stokes and anti-Stokes processes contribute to both the loss and gain resonances. chi((3)) predicts loss and gain signals with equal intensity for electronically off-resonant excitation. Some asymmetry may exist for resonant excitation. However, this is unrelated to the Stokes vs anti-Stokes processes. Any observed temperature-dependent asymmetry must thus originate from effects lying outside the chi((3)) regime.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed.

Penetrative phototactic bioconvection in an isotropic scattering suspension

Phototaxis is a directed swimming response dependent upon the light intensity sensed by micro-organisms. Positive (negative) phototaxis denotes the motion directed towards (away from) the source of light. Using the phototaxis model of Ghorai, Panda, and Hill ''Bioconvection in a suspension of isotropically scattering phototactic algae,'' Phys. Fluids 22, 071901 (2010)], we investigate two-dimensional phototactic bioconvection in an absorbing and isotropic scattering suspension in the nonlinear regime. The suspension is confined by a rigid bottom boundary, and stress-free top and lateral boundaries. The governing equations for phototactic bioconvection consist of Navier-Stokes equations for an incompressible fluid coupled with a conservation equation for micro-organisms and the radiative transfer equation for light transport. The governing system is solved efficiently using a semi-implicit second-order accurate conservative finite-difference method. The radiative transfer equation is solved by the finite volume method using a suitable step scheme. The resulting bioconvective patterns differ qualitatively from those found by Ghorai and Hill ''Penetrative phototactic bioconvection,'' Phys. Fluids 17, 074101 (2005)] at a higher critical wavelength due to the effects of scattering. The solutions show transition from steady state to periodic oscillations as the governing parameters are varied. Also, we notice the accumulation of micro-organisms in two horizontal layers at two different depths via their mean swimming orientation profile for some governing parameters at a higher scattering albedo. (C) 2013 AIP Publishing LLC.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed. (C) 2013 AIP Publishing LLC.

MHD arcjet onset boundary of high specific impulse Plasma Thrusters

Magnetoplasmadynamic thrusters are known to enter a strongly unstable regime, calledas onset in the literature, under high specific impulse operation. This paper probes the early signs of onset in relatively moderate specific impulse operation by a single fluid plasma thruster simulation. The procedure involves solving the combined Maxwell’s-Navier-Stokes equation, with an onset criterion of radial current reaching close to zero values near the electrodes. Thruster parameters are varied starting from voltage potential, plasma temperature and cathodic radius. Onset curves are plotted which can provide important engine-specific information in order to understand the onset performance of the plasma thruster.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Analysis of Annular Plug Nozzle Flowfield

The flowfields associated with truncated annular plug nozzles of varying lengths are studied both experimentally and using computational tools. The nozzles are designed to observe wake structure transition for the range of pressure ratios considered. A classification of the open wake regime is proposed for comparing and analyzing the plug flowfields. The three-dimensional relief experienced by the annular plug flow leads to greater wave interactions on the plug surface as compared with linear plug flow, resulting in a delayed transition of the base wake. The Reynolds averaged Navier-Stokes based solvers employed in the studies could predict the plug surface flow accurately, whereas they exhibited limitations with regard to plug base flow predictions. Based on the experimental data generated, an empirical model for predicting closed wake base pressure is proposed and compared with other models available in literature.

Photophysical, electrochemical and solid state properties of diketopyrrolopyrrole based molecular materials: importance of the donor group

Diketopyrrolopyrrole (DPP) based molecular semiconductors have emerged as promising materials for high performance active layers in organic solar cells. It is imperative to comprehend the origin of such a property by investigating the fundamental structure property correlation. In this report we have investigated the role of the donor group in DPP based donor-acceptor- donor (D-A-D) structure to govern the solid state, photophysical and electrochemical properties. We have prepared three derivatives of DPP with varying strengths of the donor groups, such as phenyl (PDPP-Hex), thiophene (TDPP-Hex) and selenophene (SeDPP-Hex). The influence of the donor units on the solid state packing was studied by single crystal X-ray diffraction. The photophysical, electrochemical and density functional theory ( DFT) results were combined to elucidate the structural and electronic properties of three DPP derivatives. We found that these DPP derivatives crystallized in the monoclinic space group P21/c and show herringbone packing in the crystal lattice. The derivatives exhibit weak p-p stacking interactions as two neighboring molecules slip away from each other with varied torsional angles at the donor units. The high torsional angle of 32 degrees ( PDPP-Hex) between the phenyl and lactam ring results in weak intramolecular interactions between the donor and acceptor, while TDPP-Hex and SeDPP-Hex show lower torsional angles of 9 degrees and 12 degrees with a strong overlap between the donor and acceptor units. The photophysical properties reveal that PDPP-Hex exhibits a high Stokes shift of 0.32 eV and SeDPP- Hex shows a high molar absorption co-efficient of 33 600 L mol -1 1 cm -1 1 with a low band gap of similar to 2.2 eV. The electrochemical studies of SeDPP- Hex indicate the pronounced effect of selenium in stabilizing the LUMO energy levels and this further emphasizes the importance of chalcogens in developing new n-type organic semiconductors for optoelectronic devices.

Anomalous power law decay in solvation dynamics of DNA: a mode coupling theory analysis of ion contribution

Several time dependent fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long-range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect, yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Finally, we propose several experiments to check the prediction of the present theoretical work.