A comparison of finite difference and finite volume methods for solving the space-fractional advection-dispersion equation with variable coefficients

Transport processes within heterogeneous media may exhibit non- classical diffusion or dispersion which is not adequately described by the classical theory of Brownian motion and Fick’s law. We consider a space-fractional advection-dispersion equation based on a fractional Fick’s law. Zhang et al. [Water Resources Research, 43(5)(2007)] considered such an equation with variable coefficients, which they dis- cretised using the finite difference method proposed by Meerschaert and Tadjeran [Journal of Computational and Applied Mathematics, 172(1):65-77 (2004)]. For this method the presence of variable coef- ficients necessitates applying the product rule before discretising the Riemann–Liouville fractional derivatives using standard and shifted Gru ̈nwald formulas, depending on the fractional order. As an alternative, we propose using a finite volume method that deals directly with the equation in conservative form. Fractionally-shifted Gru ̈nwald formulas are used to discretise the Riemann–Liouville fractional derivatives at control volume faces, eliminating the need for product rule expansions. We compare the two methods for several case studies, highlighting the convenience of the finite volume approach.

Ion-acoustic waves in a complex plasma with negative ions

A wave propagation in a complex dusty plasma with negative ions was considered. The relevant processes such as ionization, electron attachment, diffusion, positive-negative ion recombination, plasma particle collisions, as well as elastic Coulomb and inelastic dust-charging collisions were taken self-consistently. It was found that the equilibrium of the plasma as well as the propagation of ion waves were modified to various degrees by these effects.

Charging and trapping of macroparticles in near-electrode regions of fluorocarbon plasmas with negative ions

Charging and trapping of macroparticles in the near-electrode region of fluorocarbon etching plasmas with negative ions is considered. The equilibrium charge and forces on particles are computed as a function of the local position in the plasma presheath and sheath. The ionic composition of the plasma corresponds to the etching experiments in 2.45 GHz surface-wave sustained and 13.56 MHz inductively coupled C4F8+Ar plasmas. It is shown that despite negligible negative ion currents collected by the particles, the negative fluorine ions affect the charging and trapping of particulates through modification of the sheath/presheath structure.

Ion-acoustic waves in a complex plasma with negative ions

Negative ions and negatively charged micro- to nano-meter sized dust grains are ubiquitous in astrophysical as well as industrial processing plasmas. The negative ions can appear in electro-negative plasmas as a result of elementary processes such as dissociative or non-dissociative electron attachment to neutrals. They are usually rather small in number, and in general do not affect the overall plasma behavior. On the other hand, since the dust grains are almost always highly negative, even in small numbers they can take up a considerable proportion of the total negative charge in the system. The presence of dusts can affect the characteristics of most collective processes of the plasma since the charge balance in both the steady and dynamic states can be significantly altered. Another situation that often occurs is that the electron number density becomes small because of their absorption by the dust grains or the discharge walls. In this case the negative ions in the plasma can play a very important role. Here, a self-consistent theory of linear waves in complex laboratory plasmas containing dust grains and negative ions is presented. A comprehensive model for such plasmas including source and sink effects associated with the presence of dust grains and negative ions is introduced. The stationary state of the plasma as well as the dispersion and damping characteristics of the waves are investigated. All relevant processes, such as ionization, diffusion, electron attachment, negative-positive ion recombination, dust charge relaxation, and dissipation due to electron and ion elastic collisions with neutrals and dust particles, as well as charging collisions with the dusts, are taken into consideration.

Shared generation of pseudo-random functions with cumulative maps

In Crypto’95, Micali and Sidney proposed a method for shared generation of a pseudo-random function f(·) among n players in such a way that for all the inputs x, any u players can compute f(x) while t or fewer players fail to do so, where 0 ≤ t < u ≤ n. The idea behind the Micali-Sidney scheme is to generate and distribute secret seeds S = s1, . . . , sd of a poly-random collection of functions, among the n players, each player gets a subset of S, in such a way that any u players together hold all the secret seeds in S while any t or fewer players will lack at least one element from S. The pseudo-random function is then computed as where f s i (·)’s are poly-random functions. One question raised by Micali and Sidney is how to distribute the secret seeds satisfying the above condition such that the number of seeds, d, is as small as possible. In this paper, we continue the work of Micali and Sidney. We first provide a general framework for shared generation of pseudo-random function using cumulative maps. We demonstrate that the Micali-Sidney scheme is a special case of this general construction.We then derive an upper and a lower bound for d. Finally we give a simple, yet efficient, approximation greedy algorithm for generating the secret seeds S in which d is close to the optimum by a factor of at most u ln 2.

Protein-energy malnutrition exists and is associated with negative outcomes in morbidly obese hospital patients

Prevalence of protein-energy malnutrition (PEM), food intake inadequacy and associated health-related outcomes in morbidly obese (Body Mass Index ≥ 40 kg/m2) acute care patients are unknown. This study reports findings in morbidly obese participants from the Australasian Nutrition Care Day Survey (ANCDS) conducted in 2010. The ANCDS was a cross-sectional survey involving acute care patients from 56 Australian and New Zealand hospitals. Hospital-based dietitians evaluated participants’ nutritional status (defined by Subjective Global Assessment, SGA) and 24-hour food intake (as 0%, 25%, 50%, 75%, and 100% of the offered food). Three months later, outcome data, including length of stay (LOS) and 90-day in-hospital mortality, were collected. Of the 3122 participants, 4% (n = 136) were morbidly obese (67% females, 55 ± 14 years, BMI: 48 ± 8 kg/m2). Eleven percent (n = 15) of the morbidly obese patients were malnourished, and most (n = 11/15, 73%)received standard hospital diets without additional nutritional support. Malnourished morbidly obese patients had significantly longer LOS and greater 90-day in-hospital mortality than well-nourished counterparts (23 days vs. 9 days, p = 0.036; 14% vs. 0% mortality, p = 0.011 respectively). Thirteen morbidly obese patients (10%) consumed only 25% of the offered meals with a significantly greater proportion of malnourished (n = 4, 27%) versus well-nourished (n = 9, 7%) (p = 0.018). These results provide new knowledge on the prevalence of PEM and poor food intake in morbidly obese patients in Australian and New Zealand hospitals. For the first time internationally, the study establishes that PEM is significantly associated with negative outcomes in morbidly obese patients and warrants timely nutritional support during hospitalisation.

Plasmon nanofocusing with negative refraction in a high-index dielectric wedge

Here, we describe a metal-insulator-insulator nanofocusing structure formed by a high-permittivity dielectric wedge on a metal substrate. The structure is shown to produce nanofocusing of surface plasmon polaritons (SPPs) in the direction opposite to the taper of the wedge, including a range of nanoplasmonic effects such as nanofocusing of SPPs with negative refraction, formation of plasmonic caustics within a nanoscale distance from the wedge tip, mutual transformation of SPP modes, and significant local field enhancements in the adiabatic and strongly nonadiabatic regimes. A combination of approximate analytical and rigorous numerical approaches is used to analyze the strength and position of caustics in the structure. In particular, it is demonstrated that strong SPP localization within spatial regions as small as a few tens of nanometers near the caustic is achievable in the considered structures. Contrary to other nanofocusing configurations, efficient nanofocusing is shown to occur in the strongly nonadiabatic regime with taper angles of the dielectric wedge as large as ∼40° and within uniquely short distances (as small as a few dozens of nanometers) from the tip of the wedge. Physical interpretations of the obtained results are also presented and discussed.

Korteweg-De Vries Equation for Non-Linear Ion-Acoustic Waves in a Plasma with Negative Ions

The theoretical analysis, based on the perturbation technique, of ion-acoustic waves in the vicinity of a Korteweg-de Vries (K-dV) equation derived in a plasma with some negative ions has been made. The investigation shows that the negative ions in plasma with isothermal electrons introduced a critical concentration at which the ion-acoustic wave plays an important role of wave-breaking and forming a precursor while the plasma with non-isothermal electrons has no such singular behaviour of the wave. These two distinct features of ion waves lead to an overall different approach of present study of ion-waves. A distinct feature of non-uniform transition from the nonisothermal case to isothermal case has been shown. Few particular plasma models have been chosen to show the characteristics behaviour of the ion-waves existing in different cases

On the stability of a nonlinear system with periodic coefficients

A quasi-geometric stability criterion for feedback systems with a linear time invariant forward block and a periodically time varying nonlinear gain in the feedback loop is developed.

Exact-Solutions Of Linear Partial-Differential Equations With Variable-Coefficients

Backlund transformations relating the solutions of linear PDE with variable coefficients to those of PDE with constant coefficients are found, generalizing the study of Varley and Seymour [2]. Auto-Backlund transformations are also determined. To facilitate the generation of new solutions via Backlund transformation, explicit solutions of both classes of the PDE just mentioned are found using invariance properties of these equations and other methods. Some of these solutions are new.

Fracture analysis of concrete using singular fractal functions with lattice beam network and confirmation with acoustic emission study

In the present study singular fractal functions (SFF) were used to generate stress-strain plots for quasibrittle material like concrete and cement mortar and subsequently stress-strain plot of cement mortar obtained using SFF was used for modeling fracture process in concrete. The fracture surface of concrete is rough and irregular. The fracture surface of concrete is affected by the concrete's microstructure that is influenced by water cement ratio, grade of cement and type of aggregate 11-41. Also the macrostructural properties such as the size and shape of the specimen, the initial notch length and the rate of loading contribute to the shape of the fracture surface of concrete. It is known that concrete is a heterogeneous and quasi-brittle material containing micro-defects and its mechanical properties strongly relate to the presence of micro-pores and micro-cracks in concrete 11-41. The damage in concrete is believed to be mainly due to initiation and development of micro-defects with irregularity and fractal characteristics. However, repeated observations at various magnifications also reveal a variety of additional structures that fall between the `micro' and the `macro' and have not yet been described satisfactorily in a systematic manner [1-11,15-17]. The concept of singular fractal functions by Mosolov was used to generate stress-strain plot of cement concrete, cement mortar and subsequently the stress-strain plot of cement mortar was used in two-dimensional lattice model [28]. A two-dimensional lattice model was used to study concrete fracture by considering softening of matrix (cement mortar). The results obtained from simulations with lattice model show softening behavior of concrete and fairly agrees with the experimental results. The number of fractured elements are compared with the acoustic emission (AE) hits. The trend in the cumulative fractured beam elements in the lattice fracture simulation reasonably reflected the trend in the recorded AE measurements. In other words, the pattern in which AE hits were distributed around the notch has the same trend as that of the fractured elements around the notch which is in support of lattice model. (C) 2011 Elsevier Ltd. All rights reserved.

On Convergence of Differential Evolution Over a Class of Continuous Functions With Unique Global Optimum

Differential evolution (DE) is arguably one of the most powerful stochastic real-parameter optimization algorithms of current interest. Since its inception in the mid 1990s, DE has been finding many successful applications in real-world optimization problems from diverse domains of science and engineering. This paper takes a first significant step toward the convergence analysis of a canonical DE (DE/rand/1/bin) algorithm. It first deduces a time-recursive relationship for the probability density function (PDF) of the trial solutions, taking into consideration the DE-type mutation, crossover, and selection mechanisms. Then, by applying the concepts of Lyapunov stability theorems, it shows that as time approaches infinity, the PDF of the trial solutions concentrates narrowly around the global optimum of the objective function, assuming the shape of a Dirac delta distribution. Asymptotic convergence behavior of the population PDF is established by constructing a Lyapunov functional based on the PDF and showing that it monotonically decreases with time. The analysis is applicable to a class of continuous and real-valued objective functions that possesses a unique global optimum (but may have multiple local optima). Theoretical results have been substantiated with relevant computer simulations.

Ion-Acoustic Solitary Waves in Multicomponent Plasmas with Negative Ions

By using the perturbation technique, a Kortewege-de-Vries (K-dV) equation for a multicomponent plasma with negative ions and isothermal electrons has been derived. We have discussed the stationary solutions of K-dV equation and it has shown that in the presece of multiple ions, the amplitude of solitons exhibits interesting behaviour, especiallY when the negative ions are present.