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.

Low-temperature plasmas in carbon nanostructure synthesis

Plasma-based techniques offer many unique possibilities for the synthesis of various nanostructures both on the surface and in the plasma bulk. In contrast to the conventional chemical vapor deposition and some other techniques, plasma-based processes ensure high level of controllability, good quality of the produced nanomaterials, and reduced environmental risk. In this work, the authors briefly review the unique features of the plasma-enhanced chemical vapor deposition approaches, namely, the techniques based on inductively coupled, microwave, and arc discharges. Specifically, the authors consider the plasmas with the ion/electron density ranging from 10^10 to 10^14 cm−3, electron energy in the discharge up to ∼10 eV, and the operating pressure ranging from 1 to 10^4 Pa (up to 105 Pa for the atmospheric-pressure arc discharges). The operating frequencies of the discharges considered range from 460 kHz for the inductively coupled plasmas, and up to 2.45 GHz for the microwave plasmas. The features of the direct-current arc discharges are also examined. The authors also discuss the principles of operation of these systems, as well as the effects of the key plasma parameters on the conditions of nucleation and growth of the carbon nanostructures, mainly carbon nanotubes and graphene. Advantages and disadvantages of these plasma systems are considered. Future trends in the development of these plasma-based systems are also discussed.

Ion-assisted functional monolayer coating of nanorod arrays in hydrogen plasmas

Uniformity of postprocessing of large-area, dense nanostructure arrays is currently one of the greatest challenges in nanoscience and nanofabrication. One of the major issues is to achieve a high level of control in specie fluxes to specific surface areas of the nanostructures. As suggested by the numerical experiments in this work, this goal can be achieved by manipulating microscopic ion fluxes by varying the plasma sheath and nanorod array parameters. The dynamics of ion-assisted deposition of functional monolayer coatings onto two-dimensional carbon nanorod arrays in a hydrogen plasma is simulated by using a multiscale hybrid numerical simulation. The numerical results show evidence of a strong correlation between the aspect ratios and nanopattern positioning of the nanorods, plasma sheath width, and densities and distributions of microscopic ion fluxes. When the spacing between the nanorods and/or their aspect ratios are larger, and/or the plasma sheath is wider, the density of microscopic ion current flowing to each of the individual nanorods increases, thus reducing the time required to apply a functional monolayer coating down to 11 s for a 7-μm-wide sheath, and to 5 s for a 50-μm-wide sheath. The computed monolayer coating development time is consistent with previous experimental reports on plasma-assisted functionalization of related carbon nanostructures [B. N. Khare et al., Appl. Phys. Lett. 81, 5237 (2002)]. The results are generic in that they can be applied to a broader range of plasma-based processes and nanostructures, and contribute to the development of deterministic strategies of postprocessing and functionalization of various nanoarrays for nanoelectronic, biomedical, and other emerging applications.

Genotoxicity of inorganic lead salts and disturbance of microtubule function

Lead compounds are known genotoxicants, principally affecting the integrity of chromosomes. Lead chloride and lead acetate induced concentration-dependent increases in micronucleus frequency in V79 cells, starting at 1.1 μM lead chloride and 0.05 μM lead acetate. The difference between the lead salts, which was expected based on their relative abilities to form complex acetato-cations, was confirmed in an independent experiment. CREST analyses of the micronuclei verified that lead chloride and acetate were predominantly aneugenic (CREST-positive response), which was consistent with the morphology of the micronuclei (larger micronuclei, compared with micronuclei induced by a clastogenic mechanism). The effects of high concentrations of lead salts on the microtubule network of V79 cells were also examined using immunofluorescence staining. The dose effects of these responses were consistent with the cytotoxicity of lead(II), as visualized in the neutral-red uptake assay. In a cell-free system, 20-60 μM lead salts inhibited tubulin assembly dose-dependently. The no-observed-effect concentration of lead(II) in this assay was 10 μM. This inhibitory effect was interpreted as a shift of the assembly/disassembly steady-state toward disassembly, e.g., by reducing the concentration of assembly-competent tubulin dimers. The effects of lead salts on microtubule-associated motor-protein functions were studied using a kinesin-gliding assay that mimics intracellular transport processes in vitro by quantifying the movement of paclitaxel-stabilized microtubules across a kinesin-coated glass surface. There was a dose-dependent effect of lead nitrate on microtubule motility. Lead nitrate affected the gliding velocities of microtubules starting at concentrations above 10 μM and reached half-maximal inhibition of motility at about 50 μM. The processes reported here point to relevant interactions of lead with tubulin and kinesin at low dose levels.

Genotoxicity of inorganic mercury salts based on disturbed microtubule function

This study investigated the hypothesis that the chromosomal genotoxicity of inorganic mercury results from interaction(s) with cytoskeletal proteins. Effects of Hg2+ salts on functional activities of tubulin and kinesin were investigated by determining tubulin assembly and kinesin-driven motility in cell-free systems. Hg2+ inhibits microtubule assembly at concentrations above 1 μM, and inhibition is complete at about 10 μM. In this range, the tubulin assembly is fully (up to 6 μM) or partially (∼6-10 μM) reversible. The inhibition of tubulin assembly by mercury is independent of the anion, chloride or nitrate. The no-observed-effect- concentration for inhibition of microtubule assembly in vitro was 1 μM Hg2+, the IC50 5.8 μM. Mercury(II) salts at the IC 50 concentrations partly inhibiting tubulin assembly did not cause the formation of aberrant microtubule structures. Effects of mercury salts on the functionality of the microtubule motility apparatus were studied with the motor protein kinesin. By using a "gliding assay" mimicking intracellular movement and transport processes in vitro, HgCl2 affected the gliding velocity of paclitaxel-stabilised microtubules in a clear dose-dependent manner. An apparent effect is detected at a concentration of 0.1 μM and a complete inhibition is reached at 1 μM. Cytotoxicity of mercury chloride was studied in V79 cells using neutral red uptake, showing an influence above 17 μM HgCl2. Between 15 and 20 μM HgCl2 there was a steep increase in cell toxicity. Both mercury chloride and mercury nitrate induced micronuclei concentration-dependently, starting at concentrations above 0.01 μM. CREST analyses on micronuclei formation in V79 cells demonstrated both clastogenic (CREST-negative) and aneugenic effects of Hg2+, with some preponderance of aneugenicity. A morphological effect of high Hg2+ concentrations (100 μM HgCl2) on the microtubule cytoskeleton was verified in V79 cells by immuno-fluorescence staining. The overall data are consistent with the concept that the chromosomal genotoxicity could be due to interaction of Hg2+ with the motor protein kinesin mediating cellular transport processes. Interactions of Hg 2+ with the tubulin shown by in vitro investigations could also partly influence intracellular microtubule functions leading, together with the effects on the kinesin, to an impaired chromosome distribution as shown by the micronucleus test.

Levels and trends of PBDEs and HBCDs in the global environment : status at the end of 2012

In this paper, we have compiled and reviewed the most recent literature, published from January2010 to December 2012, relating to the human exposure, environmental distribution, behaviour, fate and concentration time trends of polybrominated diphenyl ether (PBDE) and hexabromocyclododecane (HBCD) flame retardants, in order to establish their current trends and priorities for future study. Due to the large volume of literature included, we have provided full detail of the reviewed studies as Electronic Supplementary Information and here summarise the most relevant findings. Decreasing time trends for penta-mix PBDE congeners were seen for soils in northern Europe, sewage sludge in Sweden and the USA, carp from a US river, trout from three of the Great Lakes and in Arctic and UK marine mammals and many birds, but increasing time trends continue in Arctic polar bears and some birds at high trophic levels in northern Europe. This is a result of the time delay inherent in long-range atmospheric transport processes. In general, concentrations of BDE209 (the major component of the deca-mix PBDE product) are continuing to increase. Of major concern is the possible/likely debromination of the large reservoir of BDE209 in soils and sediments worldwide, to yield lower brominated congeners which are both more mobile and more toxic, and we have compiled the most recent evidence for the occurrence of this degradation process. Numerous studies reported here reinforce the importance o f this future concern. Time trends for HBCDs are mixed, with both increases and decreases evident in different matrices and locations and, notably, with increasing occurrence in birds of prey.

Modelling intermittent microwave convective drying (IMCD) of food materials

This thesis develops comprehensive mathematical models for an advanced drying technology Intermittent Microwave Convective Drying (IMCD). The models provide an improved physical understanding of the heat and mass transport during the drying process, which will help to improve the quality of dried food and energy efficiency of the process, as well as will increase the ability of automation and optimization. The final model in this thesis represents the most comprehensive fundamental multiphase model for IMCD that considers 3D electromagnetics coupled with multiphase porous media transport processes. The 3D electromagnetics considers Maxwell's equation and multiphase transport model considers three different phases: solid matrix, liquid water and gas consisting water vapour and air. The multiphase transport includes pressure-driven flow, capillary diffusion, binary diffusion, and evaporation. The models developed in this thesis were validated with extensive experimental investigations.

On the order of the fractional Laplacian in determining the spatio-temporal evolution of a space-fractional model of cardiac electrophysiology

Space-fractional operators have been used with success in a variety of practical applications to describe transport processes in media characterised by spatial connectivity properties and high structural heterogeneity altering the classical laws of diffusion. This study provides a systematic investigation of the spatio-temporal effects of a space-fractional model in cardiac electrophysiology. We consider a simplified model of electrical pulse propagation through cardiac tissue, namely the monodomain formulation of the Beeler-Reuter cell model on insulated tissue fibres, and obtain a space-fractional modification of the model by using the spectral definition of the one-dimensional continuous fractional Laplacian. The spectral decomposition of the fractional operator allows us to develop an efficient numerical method for the space-fractional problem. Particular attention is paid to the role played by the fractional operator in determining the solution behaviour and to the identification of crucial differences between the non-fractional and the fractional cases. We find a positive linear dependence of the depolarization peak height and a power law decay of notch and dome peak amplitudes for decreasing orders of the fractional operator. Furthermore, we establish a quadratic relationship in conduction velocity, and quantify the increasingly wider action potential foot and more pronounced dispersion of action potential duration, as the fractional order is decreased. A discussion of the physiological interpretation of the presented findings is made.

Calculating the Fickian diffusivity for a lattice-based random walk with agents and obstacles of different shapes and sizes

Random walk models are often used to interpret experimental observations of the motion of biological cells and molecules. A key aim in applying a random walk model to mimic an in vitro experiment is to estimate the Fickian diffusivity (or Fickian diffusion coefficient),D. However, many in vivo experiments are complicated by the fact that the motion of cells and molecules is hindered by the presence of obstacles. Crowded transport processes have been modeled using repeated stochastic simulations in which a motile agent undergoes a random walk on a lattice that is populated by immobile obstacles. Early studies considered the most straightforward case in which the motile agent and the obstacles are the same size. More recent studies considered stochastic random walk simulations describing the motion of an agent through an environment populated by obstacles of different shapes and sizes. Here, we build on previous simulation studies by analyzing a general class of lattice-based random walk models with agents and obstacles of various shapes and sizes. Our analysis provides exact calculations of the Fickian diffusivity, allowing us to draw conclusions about the role of the size, shape and density of the obstacles, as well as examining the role of the size and shape of the motile agent. Since our analysis is exact, we calculateDdirectly without the need for random walk simulations. In summary, we find that the shape, size and density of obstacles has a major influence on the exact Fickian diffusivity. Furthermore, our results indicate that the difference in diffusivity for symmetric and asymmetric obstacles is significant.

Modeling of heat and mass transfer for solid state fermentation process in tray bioreactor

A mathematical model is developed to simulate oxygen consumption, heat generation and cell growth in solid state fermentation (SSF). The fungal growth on the solid substrate particles results in the increase of the cell film thickness around the particles. The model incorporates this increase in the biofilm size which leads to decrease in the porosity of the substrate bed and diffusivity of oxygen in the bed. The model also takes into account the effect of steric hindrance limitations in SSF. The growth of cells around single particle and resulting expansion of biofilm around the particle is analyzed for simplified zero and first order oxygen consumption kinetics. Under conditions of zero order kinetics, the model predicts upper limit on cell density. The model simulations for packed bed of solid particles in tray bioreactor show distinct limitations on growth due to simultaneous heat and mass transport phenomena accompanying solid state fermentation process. The extent of limitation due to heat and/or mass transport phenomena is analyzed during different stages of fermentation. It is expected that the model will lead to better understanding of the transport processes in SSF, and therefore, will assist in optimal design of bioreactors for SSF.

Modelling boreal forest CO2 exchange and seasonality

Man-induced climate change has raised the need to predict the future climate and its feedback to vegetation. These are studied with global climate models; to ensure the reliability of these predictions, it is important to have a biosphere description that is based upon the latest scientific knowledge. This work concentrates on the modelling of the CO2 exchange of the boreal coniferous forest, studying also the factors controlling its growing season and how these can be used in modelling. In addition, the modelling of CO2 gas exchange at several scales was studied. A canopy-level CO2 gas exchange model was developed based on the biochemical photosynthesis model. This model was first parameterized using CO2 exchange data obtained by eddy covariance (EC) measurements from a Scots pine forest at Sodankylä. The results were compared with a semi-empirical model that was also parameterized using EC measurements. Both of the models gave satisfactory results. The biochemical canopy-level model was further parameterized at three other coniferous forest sites located in Finland and Sweden. At all the sites, the two most important biochemical model parameters showed seasonal behaviour, i.e., their temperature responses changed according to the season. Modelling results were improved when these changeover dates were related to temperature indices. During summer-time the values of the biochemical model parameters were similar at all the four sites. Different control factors for CO2 gas exchange were studied at the four coniferous forests, including how well these factors can be used to predict the initiation and cessation of the CO2 uptake. Temperature indices, atmospheric CO2 concentration, surface albedo and chlorophyll fluorescence (CF) were all found to be useful and have predictive power. In addition, a detailed simulation study of leaf stomata in order to separate physical and biochemical processes was performed. The simulation study brought to light the relative contribution and importance of the physical transport processes. The results of this work can be used in improving CO2 gas exchange models in boreal coniferous forests. The meteorological and biological variables that represent the seasonal cycle were studied, and a method for incorporating this cycle into a biochemical canopy-level model was introduced.

Segmental relaxations and crystallization-induced phase separation in PVDF/PMMA blends in the presence of surface-functionalized multiwall carbon nanotubes

Crystallization-induced phase separation and segmental relaxations in poly(vinylidene fluoride)/poly(methyl methacrylate) (PVDF/PMMA) blends was systematically investigated by melt-rheology and broadband dielectric spectroscopy in the presence of multiwall carbon nanotubes (MWNTs). Different functionalized MWNTs (amine, -NH2; acid, -COOH) were incorporated in the blends by melt-mixing above the melting temperature of PVDF, where the blends are miscible, and the crystallization induced phase separation was probed in situ by shear rheology. Interestingly, only -NH2 functionalized MWNTs (a-MWNTs) aided in the formation of beta-phase (trans-trans) crystals in PVDF, whereas both the neat blends and the blends with -COOH functionalized MWNTs (c-MWNTs) showed only alpha-phase (trans-gauche-trans-gauche') crystals as inferred from wide-angle X-ray diffraction (WXRD) and Fourier transform infrared (FTIR). Furthermore, blends with only a-MWNTs facilitated in heterogeneous nucleation in the blends manifesting in an increase in the calorimetric crystallization temperature and hence, augmented the theologically determined crystallintion induced phase separation temperature. The dielectric relaxations associated with the crystalline phase of PVDF (alpha(c)) was completely absent in the blends with a-MWNTs in contrast to neat blends and the blends with c-MWNTs in the dielectric loss spectra. The relaxations in the blends investigated here appeared to follow Havriliak-Negami (HN) empirical equations, and, more interestingly, the dynamic heterogeneity in the system could be mapped by an extra relaxation at higher frequency at the crystallization-induced phase separation temperature. The mean relaxation time (tau(HN)) was evaluated and observed to be delayed in the presence of MWNTs in the blends, more prominently in the case of blends with a-MWNTs. The latter also showed a significant increase in the dielectric relaxation strength (Delta epsilon). Electron microscopy and selective etching was used to confirm the localization of MWNTs in the amorphous phases of the interspherulitic regions as observed from scanning electron microscopy (SEM). The evolved crystalline morphology, during crystallization-induced phase separation, was observed to have a strong influence on the charge transport processes in the blends. These observations were further supported by the specific interactions (like dipole induced dipole interaction) between a-MWNTs and PVDF, as inferred from FTIR, and the differences in the crystalline morphology as observed from WXRD and polarized optical microscopy (POM).

Analytical solutions for transient temperature distribution in a geothermal reservoir due to cold water injection

An analytical solution to describe the transient temperature distribution in a geothermal reservoir in response to injection of cold water is presented. The reservoir is composed of a confined aquifer, sandwiched between rocks of different thermo-geological properties. The heat transport processes considered are advection, longitudinal conduction in the geothermal aquifer, and the conductive heat transfer to the underlying and overlying rocks of different geological properties. The one-dimensional heat transfer equation has been solved using the Laplace transform with the assumption of constant density and thermal properties of both rock and fluid. Two simple solutions are derived afterwards, first neglecting the longitudinal conductive heat transport and then heat transport to confining rocks. Results show that heat loss to the confining rock layers plays a vital role in slowing down the cooling of the reservoir. The influence of some parameters, e.g. the volumetric injection rate, the longitudinal thermal conductivity and the porosity of the porous media, on the transient heat transport phenomenon is judged by observing the variation of the transient temperature distribution with different values of the parameters. The effects of injection rate and thermal conductivity have been found to be profound on the results.

Temperature dependent current-voltage characteristics of zinc oxide nanowire/polypyrrole nanocomposite

Temperature dependent current-voltage (I-V) measurements of electrochemically prepared zinc oxide nanowire/polypyrrole (ZnONW/PPy) nanocomposite yielded non-linear I-V characteristics at temperatures between 300 and 4.5 K. The low-field conductance (G) of the ZnONW/PPy film exhibits pronounced temperature dependence with room temperature conductance (G(300K)) similar to 10(-3) S and a conductance ratio (G(300)K/G(4.5K)) of similar to 10(4), indicating dominance of significant temperature dependent charge transport processes. The conduction mechanism of the film is satisfactorily understood by extended fluctuation induced tunneling (FIT) model as the non-linear I-V characteristics fit fairly well to the extended FIT model. Further, the temperature dependence of G(o) obtained from fitting followed Sheng's model also. (C) 2014 AIP Publishing LLC.

Thermally developing combined electroosmotic and pressure-driven flow of nanofluids in a microchannel under the effect of magnetic field

In the present study, the heat transfer characteristics of thermally developing magnetohydroclynamic flow of nanofluid through microchannel are delineated by following a semi analytical approach. The combined influences of pressure driven flow, electroosmotic transport and magnetic field is taken into account for the analysis of the complex microscale thermal transport processes. Solutions for the normalized temperature distributions and the Nusselt number variations, considering the simultaneous interplay of electrokinetic effects (electroosmosis), magnetic effects, Joule heating and viscous dissipation are obtained, for constant wall temperature condition. Particular attention is paid to assess the role of nanolluids in altering the transport phenomena, through variations in the effective nanoparticle volume fractions, as well as the aggregate structure of the particulate phases. It is observed that magnetohydrodynamic effect reduces advective transport of the liquid resulting in gradual reduction of heat transfer. Increase in nanoparticle volume fraction shows decrease in heat transfer. Similar effects are observed with increase in aggregate sizes of the nanoparticles. The effect of the nanofluids on system irreversibility is also studied through entropy generation analysis due to flow and heat transfer in the microchannel. Total entropy generation is found to be dominant at the thermally developing region of the microchannel, whereas it drops sharply at the thermally developed region. Presence of nanoparticles in the base fluid reduces the total entropy generation in the microchannel, thereby indicating decrease in thermodynamic irreversibility with increasing nanoparticle volume fraction. (C) 2014 Elsevier Ltd. All rights reserved.