Modeling harvest rates and numbers from age and sex ratios: A demonstration for elephant populations

Illegal harvest rates of wildlife populations are often unknown or difficult to estimate from field data due to under-reporting or incomplete detection of carcasses. This is especially true for elephants that are killed for ivory or in conflicts with people. We describe a method to infer harvest rates from coarse field data of three population parameters, namely, adult female to male ratio, male old-adult to young-adult ratio, and proportion of adult males in the population using Jensen's (2000) 2-sex, density-dependent Leslie matrix model. The specific combination of male and female harvest rates and numbers can be determined from the history of harvest and estimate of population size. We applied this technique to two populations of elephants for which data on age structure and records of mortality were available-a forest-dwelling population of the Asian elephant (at Nagarahole, India) and an African savannah elephant population (at Samburu, Kenya) that had experienced male-biased harvest regimes over 2-3 decades. For the Nagarahole population, the recorded numbers of male and female elephants killed illegally during 1981-2000 were 64% and 88% of the values predicted by the model, respectively, implying some non-detection or incomplete reporting while for the Samburu population the recorded and modeled numbers of harvest during 1990-1999 closely matched. This technique, applicable to any animal population following logistic growth model, can be especially useful for inferring illegal harvest numbers of forest elephants in Africa and Asia. (C) 2013 Elsevier Ltd. All rights reserved.

A Comparative Study of Early Afterdepolarization-Mediated Fibrillation in Two Mathematical Models for Human Ventricular Cells

Early afterdepolarizations (EADs), which are abnormal oscillations of the membrane potential at the plateau phase of an action potential, are implicated in the development of cardiac arrhythmias like Torsade de Pointes. We carry out extensive numerical simulations of the TP06 and ORd mathematical models for human ventricular cells with EADs. We investigate the different regimes in both these models, namely, the parameter regimes where they exhibit (1) a normal action potential (AP) with no EADs, (2) an AP with EADs, and (3) an AP with EADs that does not go back to the resting potential. We also study the dependence of EADs on the rate of at which we pace a cell, with the specific goal of elucidating EADs that are induced by slow or fast rate pacing. In our simulations in two-and three-dimensional domains, in the presence of EADs, we find the following wave types: (A) waves driven by the fast sodium current and the L-type calcium current (Na-Ca-mediated waves); (B) waves driven only by the L-type calcium current (Ca-mediated waves); (C) phase waves, which are pseudo-travelling waves. Furthermore, we compare the wave patterns of the various wave-types (Na-Ca-mediated, Ca-mediated, and phase waves) in both these models. We find that the two models produce qualitatively similar results in terms of exhibiting Na-Ca-mediated wave patterns that are more chaotic than those for the Ca-mediated and phase waves. However, there are quantitative differences in the wave patterns of each wave type. The Na-Ca-mediated waves in the ORd model show short-lived spirals but the TP06 model does not. The TP06 model supports more Ca-mediated spirals than those in the ORd model, and the TP06 model exhibits more phase-wave patterns than does the ORd model.

Internal noise-driven generalized Langevin equation from a nonlocal continuum model

Starting with a micropolar formulation, known to account for nonlocal microstructural effects at the continuum level, a generalized Langevin equation (GLE) for a particle, describing the predominant motion of a localized region through a single displacement degree of freedom, is derived. The GLE features a memory-dependent multiplicative or internal noise, which appears upon recognizing that the microrotation variables possess randomness owing to an uncertainty principle. Unlike its classical version, the present GLE qualitatively reproduces the experimentally measured fluctuations in the steady-state mean square displacement of scattering centers in a polyvinyl alcohol slab. The origin of the fluctuations is traced to nonlocal spatial interactions within the continuum, a phenomenon that is ubiquitous across a broad class of response regimes in solids and fluids. This renders the proposed GLE a potentially useful model in such cases.

THE EFFECTS OF SURFACTANT ON SIMPLEX NOZZLE SPRAY BEHAVIOR AND ITS COMPARISON TO LIQUID FUELS

Pressure-swirl nozzles (simplex nozzles) are used in various field applications such as aero-engines, power generation, spray painting and agricultural irrigation. For this particular nozzle, research in the past decade has dealt with the development of numerical models for predicting droplet distribution profiles. Although these results have been valuable, the experimental results have been contradictory, therefore fundamental understanding of the influence of properties in nozzle is important. This paper experimentally investigates the effect of surfactants on breakup and coalescence. Since most of the fuels and biofuels have low surface tension compared to water, a comparative analysis between a surfactant solution and a liquid fuel is imperative. For this experimental study, a simplex nozzle characterized as flow number 0.4 will be utilized. The injection pressures will range from 0.3 - 4Mpa while altering the surface tension from 72 to 28mN/m. By applying Phase Doppler Particle Anemometry (PDPA) which is a non-intrusive laser diagnostic technique, the differences in spray characteristics due to spray surface tension can be highlighted. The average droplet diameter decreases for a low surface tension fluid in the axial direction in comparison to pure water. The average velocity of droplets is surprisingly lower in the same spray zone. Measurements made in the radial direction show no significant changes, but at the locations close to the nozzle, water droplets have larger diameter and velocity. The results indicate the breakup and coalescence regimes have been altered when surface tension is lowered. A decrease in surface tension alters the breakup length while increasing the spray angle. Moreover, higher injection pressure shortens the breakup length and decrease in overall diameter of the droplets. By performing this experimental study the fundamentals of spray dynamics, such as spray formation, liquid breakup length, and droplet breakup regimes can be observed as a function of surface tension and how a surrogate fuel compares with a real fuel for experimental purposes. This knowledge potentially will lead to designing a better atomizer or new biofuels.

Homogeneous isotropic superfluid turbulence in two dimensions: Inverse and forward cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the simplified, Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normalfluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli.

Universal buckling kinetics in drying nanoparticle-laden droplets on a hydrophobic substrate

We provide a comprehensive physical description of the vaporization, self-assembly, agglomeration, and buckling kinetics of sessile nanofluid droplets pinned on a hydrophobic substrate. We have deciphered five distinct regimes of the droplet life cycle. Regimes I-III consists of evaporation-induced preferential agglomeration that leads to the formation of a unique dome-shaped inhomogeneous shell with a stratified varying-density liquid core. Regime IV involves capillary-pressure-initiated shell buckling and stress-induced shell rupture. Regime V marks rupture-induced cavity inception and growth. We demonstrate through scaling arguments that the growth of the cavity (which controls the final morphology or structure) can be described by a universal function.

General framework for acoustic emission during plastic deformation

Despite the long history, so far there is no general theoretical framework for calculating the acoustic emission spectrum accompanying any plastic deformation. We set up a discrete wave equation with plastic strain rate as a source term and include the Rayleigh-dissipation function to represent dissipation accompanying acoustic emission. We devise a method of bridging the widely separated time scales of plastic deformation and elastic degrees of freedom. While this equation is applicable to any type of plastic deformation, it should be supplemented by evolution equations for the dislocation microstructure for calculating the plastic strain rate. The efficacy of the framework is illustrated by considering three distinct cases of plastic deformation. The first one is the acoustic emission during a typical continuous yield exhibiting a smooth stress-strain curve. We first construct an appropriate set of evolution equations for two types of dislocation densities and then show that the shape of the model stress-strain curve and accompanying acoustic emission spectrum match very well with experimental results. The second and the third are the more complex cases of the Portevin-Le Chatelier bands and the Luders band. These two cases are dealt with in the context of the Ananthakrishna model since the model predicts the three types of the Portevin-Le Chatelier bands and also Luders-like bands. Our results show that for the type-C bands where the serration amplitude is large, the acoustic emission spectrum consists of well-separated bursts of acoustic emission. At higher strain rates of hopping type-B bands, the burst-type acoustic emission spectrum tends to overlap, forming a nearly continuous background with some sharp acoustic emission bursts. The latter can be identified with the nucleation of new bands. The acoustic emission spectrum associated with the continuously propagating type-A band is continuous. These predictions are consistent with experimental results. More importantly, our study shows that the low-amplitude continuous acoustic emission spectrum seen in both the type-B and type-A band regimes is directly correlated to small-amplitude serrations induced by propagating bands. The acoustic emission spectrum of the Luders-like band matches with recent experiments as well. In all of these cases, acoustic emission signals are burstlike, reflecting the intermittent character of dislocation-mediated plastic flow.

Injection barrier induced deviations in space charge limited conduction in doped poly(3-methylthiophene) based devices

The carrier density dependent current-voltage (J V) characteristics of electrochemically prepared poly(3-methylthiophene) (P3MeT) have been investigated in Pt/P3MeT/Al devices, as a function of temperature from 280 to 84 K. In these devices, the charge transport is found to be mainly governed by different transport regimes of space charge limited conduction (SCLC). In a lightly doped device, SCLC controlled by exponentially distributed traps (Vl+1 law, l > 1) is observed in the intermediate voltage range (0.5-2 V) at all temperatures. However, at higher bias (> 2 V), the current deviates from the usual Vl+1 law where the slope is found to be less than 2 of the logJ-logV plot, which is attributed to the presence of the injection barrier. These deviations gradually disappear at higher doping level due to reduction in the injection barrier. Numerical simulations of the Vl+1 law by introducing the injection barrier show good agreement with experimental data. The results show that carrier density can tune the charge transport mechanism in Pt/P3MeT/Al devices to understand the non-Ohmic behavior. The plausible reasons for the origin of injection barrier and the transitions in the transport mechanism with carrier density are discussed. (C) 2015 AIP Publishing LLC.

Onset of plane layer magnetoconvection at low Ekman number

We study the onset of magnetoconvection between two infinite horizontal planes subject to a vertical magnetic field aligned with background rotation. In order to gain insight into the convection taking place in the Earth's tangent cylinder, we target regimes of asymptotically strong rotation. The critical Rayleigh number Ra-c and critical wavenumber a(c) are computed numerically by solving the linear stability problem in a systematic way, with either stress-free or no-slip kinematic boundary conditions. A parametric study is conducted, varying the Ekman number E (ratio of viscous to Coriolis forces) and the Elsasser number. (ratio of the Lorentz force to the Coriolis force). E is varied from 10(-9) to 10(-2) and. from 10(-3) to 1. For a wide range of thermal and magnetic Prandtl numbers, our results verify and confirm previous experimental and theoretical results showing the existence of two distinct unstable modes at low values of E-one being controlled by the magnetic field, the other being controlled by viscosity (often called the viscous mode). It is shown that oscillatory onset does not occur in the range of parameters we are interested in. Asymptotic scalings for the onset of these modes are numerically confirmed and their domain of validity is precisely quantified. We show that with no-slip boundary conditions, the asymptotic behavior is reached for E < 10(-6) and establish a map in the (E, Lambda) plane. We distinguish regions where convection sets in either through the magnetic mode or through the viscous mode. Our analysis gives the regime in which the transition between magnetic and viscous modes may be observed. We also show that within the asymptotic regime, the role played by the kinematic boundary conditions is minimal. (C) 2015 AIP Publishing LLC.

Efficient coupled polynomial interpolation scheme for out-of-plane free vibration analysis of curved beams

The performance of two curved beam finite element models based on coupled polynomial displacement fields is investigated for out-of-plane vibration of arches. These two-noded beam models employ curvilinear strain definitions and have three degrees of freedom per node namely, out-of-plane translation (v), out-of-plane bending rotation (theta(z)) and torsion rotation (theta(s)). The coupled polynomial interpolation fields are derived independently for Timoshenko and Euler-Bernoulli beam elements using the force-moment equilibrium equations. Numerical performance of these elements for constrained and unconstrained arches is compared with the conventional curved beam models which are based on independent polynomial fields. The formulation is shown to be free from any spurious constraints in the limit of `flexureless torsion' and `torsionless flexure' and hence devoid of flexure and torsion locking. The resulting stiffness and consistent mass matrices generated from the coupled displacement models show excellent convergence of natural frequencies in locking regimes. The accuracy of the shear flexibility added to the elements is also demonstrated. The coupled polynomial models are shown to perform consistently over a wide range of flexure-to-shear (EI/GA) and flexure-to-torsion (EI/GJ) stiffness ratios and are inherently devoid of flexure, torsion and shear locking phenomena. (C) 2015 Elsevier B.V. All rights reserved.

Insight into the Evaporation Dynamics of a Pair of Sessile Droplets on a Hydrophobic Substrate

In this work, we have demonstrated three unique regimes in the evaporation lifecycle of a pair of sessile droplets placed in variable proximity on a hydrophobic substrate. For small separation distance, the droplets undergo asymmetric spatiotemporal,evaporation leading to contact angle hysteresis and suppressed vaporization. The reduced evaporation has been attributed quantitatively to the existence of a constrained vapor-rich dome between the two droplets. However, a dynamic decrease in the droplet radius due to solvent removal marks a return to symmetry in terms of evaporation and contact angle. We have described the variation in evaporation flux using a universal correction factor. We have also demonstrated the existence of a critical separation distance beyond which the droplets in the, droplet pair do not affect each other. The results are crucial to a plethora of applications ranging from surface patterning to lab-on-a-chip devices.

HOMOGENIZATION OF A HYPERBOLIC EQUATION WITH HIGHLY CONTRASTING DIFFUSIVITY COEFFICIENTS

We study a hyperbolic problem in the framework of periodic homogenization assuming a high contrast between the diffusivity coefficients of the two components M-epsilon and B-epsilon of the heterogeneous medium. There are three regimes depending on the ratio between the size of the period and the amplitude a, of the diffusivity in B-epsilon. For the critical regime alpha(epsilon) similar or equal to epsilon, the limit problem is a strongly coupled system involving both the macroscopic and the microscopic variables. We also include the results in the non critical case.

Insight into the Evaporation Dynamics of a Pair of Sessile Droplets on a Hydrophobic Substrate

In this work, we have demonstrated three unique regimes in the evaporation lifecycle of a pair of sessile droplets placed in variable proximity on a hydrophobic substrate. For small separation distance, the droplets undergo asymmetric spatiotemporal,evaporation leading to contact angle hysteresis and suppressed vaporization. The reduced evaporation has been attributed quantitatively to the existence of a constrained vapor-rich dome between the two droplets. However, a dynamic decrease in the droplet radius due to solvent removal marks a return to symmetry in terms of evaporation and contact angle. We have described the variation in evaporation flux using a universal correction factor. We have also demonstrated the existence of a critical separation distance beyond which the droplets in the, droplet pair do not affect each other. The results are crucial to a plethora of applications ranging from surface patterning to lab-on-a-chip devices.

HOMOGENIZATION OF A HYPERBOLIC EQUATION WITH HIGHLY CONTRASTING DIFFUSIVITY COEFFICIENTS

We study a hyperbolic problem in the framework of periodic homogenization assuming a high contrast between the diffusivity coefficients of the two components M-epsilon and B-epsilon of the heterogeneous medium. There are three regimes depending on the ratio between the size of the period and the amplitude a, of the diffusivity in B-epsilon. For the critical regime alpha(epsilon) similar or equal to epsilon, the limit problem is a strongly coupled system involving both the macroscopic and the microscopic variables. We also include the results in the non critical case.

Deformation pathways and breakup modes in acoustically levitated bicomponent droplets under external heating

Controlled breakup of droplets using heat or acoustics is pivotal in applications such as pharmaceutics, nanoparticle production, and combustion. In the current work we have identified distinct thermal acoustics-induced deformation regimes (ligaments and bubbles) and breakup dynamics in externally heated acoustically levitated bicomponent (benzene-dodecane) droplets with a wide variation in volatility of the two components (benzene is significantly more volatile than dodecane). We showcase the physical mechanism and universal behavior of droplet surface caving in leading to the inception and growth of ligaments. The caving of the top surface is governed by a balance between the acoustic pressure field and the restrictive surface tension of the droplet. The universal collapse of caving profiles for different benzene concentration (<70% by volume) is shown by using an appropriate time scale obtained from force balance. Continuous caving leads to the formation of a liquid membrane-type structure which undergoes radial extension due to inertia gained during the precursor phase. The membrane subsequently closes at the rim and the kinetic energy leads to ligament formation and growth. Subsequent ligament breakup is primarily Rayleigh-Plateau type. The breakup mode shifts to diffusional entrapment-induced boiling with an increase in concentration of the volatile component (benzene >70% by volume). The findings are portable to any similar bicomponent systems with differential volatility.