A finite volume method for solving the two-sided time-space fractional advection-dispersion equation

The field of fractional differential equations provides a means for modelling transport processes within complex media which are governed by anomalous transport. Indeed, the application to anomalous transport has been a significant driving force behind the rapid growth and expansion of the literature in the field of fractional calculus. In this paper, we present a finite volume method to solve the time-space two-sided fractional advection dispersion equation on a one-dimensional domain. Such an equation allows modelling different flow regime impacts from either side. The finite volume formulation provides a natural way to handle fractional advection-dispersion equations written in conservative form. The novel spatial discretisation employs fractionally-shifted Gr¨unwald formulas to discretise the Riemann-Liouville fractional derivatives at control volume faces in terms of function values at the nodes, while the L1-algorithm is used to discretise the Caputo time fractional derivative. Results of numerical experiments are presented to demonstrate the effectiveness of the approach.

An exploratory study of the potential of resurfacing articular cartilage with synthetic phospholipids

This thesis is aimed at further understanding the uppermost lipid-filled membranous layer (i.e. surface amorphous layer (SAL)) of articular cartilage and to develop a scientific framework for re-introducing lipids onto the surface of lipid-depleted articular cartilage (i.e. "resurfacing"). The outcome will potentially contribute to knowledge that will facilitate the repair of the articular surface of cartilage where degradation is limited to the loss of the lipids of the SAL only. The surface amorphous layer is of utmost importance to the effective load-spreading, lubrication, and semipermeability (which controls its fluid management, nutrient transport and waste removal) of articular cartilage in the mammalian joints. However, because this uppermost layer of cartilage is often in contact during physiological function, it is prone to wear and tear, and thus, is the site for damage initiation that can lead to the early stages of joint condition like osteoarthritis, and related conditions that cause pain and discomfort leading to low quality of life in patients. It is therefore imperative to conduct a study which offers insight into remedying this problem. It is hypothesized that restoration (resurfacing) of the surface amorphous layer can be achieved by re-introducing synthetic surface-active phospholipids (SAPL) into the joint space. This hypothesis was tested in this thesis by exposing cartilage samples whose surface lipids had been depleted to individual and mixtures of synthetic saturated and unsaturated phospholipids. The surfaces of normal, delipidized, and relipidized samples of cartilage were characterized for their structural integrity and functionality using atomic force microscope (AFM), confocal microscope (COFM), Raman spectroscopy, magnetic resonance imaging (MRI) with image processing in the MATLAB® environment and mechanical loading experiments. The results from AFM imaging, confocal microscopy, and Raman spectroscopy revealed a successful deposition of new surface layer on delipidized cartilage when incubated in synthetic phospholipids. The relipidization resulted in a significant improvement in the surface nanostructure of the artificially degraded cartilage, with the complete SAPL mixture providing better outcomes in comparison to those created with the single SAPL components (palmitoyl-oleoyl-phosphatidylcholine, POPC and dipalmitoyl-phosphatidylcholine, DPPC). MRI analysis revealed that the surface created with the complete mixture of synthetic lipids was capable of providing semipermeability to the surface layer of the treated cartilage samples relative to the normal intact surface. Furthermore, deformation energy analysis revealed that the treated samples were capable of delivering the elastic properties required for load bearing and recovery of the tissue relative to the normal intact samples, with this capability closer between the normal and the samples incubated in the complete lipid mixture. In conclusion, this thesis has established that it is possible to deposit/create a potentially viable layer on the surface of cartilage following degradation/lipid loss through incubation in synthetic lipid solutions. However, further studies will be required to advance the ideas developed in this thesis, for the development of synthetic lipid-based injections/drugs for treatment of osteoarthritis and other related joint conditions.

A finite volume method for solving the two-sided time-space fractional advection-dispersion equation

We present a finite volume method to solve the time-space two-sided fractional advection-dispersion equation on a one-dimensional domain. The spatial discretisation employs fractionally-shifted Grünwald formulas to discretise the Riemann-Liouville fractional derivatives at control volume faces in terms of function values at the nodes. We demonstrate how the finite volume formulation provides a natural, convenient and accurate means of discretising this equation in conservative form, compared to using a conventional finite difference approach. Results of numerical experiments are presented to demonstrate the effectiveness of the approach.

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.

A fast semi-implicit difference method for a nonlinear two-sided space-fractional diffusion equation with variable diffusivity coefficients

In this paper, we derive a new nonlinear two-sided space-fractional diffusion equation with variable coefficients from the fractional Fick’s law. A semi-implicit difference method (SIDM) for this equation is proposed. The stability and convergence of the SIDM are discussed. For the implementation, we develop a fast accurate iterative method for the SIDM by decomposing the dense coefficient matrix into a combination of Toeplitz-like matrices. This fast iterative method significantly reduces the storage requirement of O(n2)O(n2) and computational cost of O(n3)O(n3) down to n and O(nlogn)O(nlogn), where n is the number of grid points. The method retains the same accuracy as the underlying SIDM solved with Gaussian elimination. Finally, some numerical results are shown to verify the accuracy and efficiency of the new method.

Development of chitosan packaging films with carvacrol and grape seed extract as antimicrobial and antioxidant agents

In order to produce packaging films with a broad spectrum of action on microorganisms, the effect of two antimicrobial (AM) to be included in the films, carvacrol and GSE were studied separately on different microorganisms. Carvacrol was more effective against the grampositive bacteria than against the gram-negative bacterium. GSE was not effective against yeast. Subsequently, a search for optimal combinations of carvacrol, GSE and the addition of chitosan (as a third component with film forming properties) was carried out. Response surface analysis showed several synergetic effects and three optimal AM combinations (OAMC) were obtained for each microorganism. The experimental validation confirmed that the optimal solutions found can successfully predict the response for each microorganism. The optimization of mixtures of the three components, but this time, using the same concentration for all microorganisms, was also studied to obtain an OAMC with wide spectrum of activity. The results of the response surface analysis showed several synergistic effects for all microorganisms. Three OAMC, OAMC-1, OAMC-2, OAMC-3, were found to be the optimal mixtures for all microorganisms. The radical scavenging activity (RSA) of the different agents was then compared with a standard antioxidant (AOX) BHT, at different concentrations; as also at the OAMC. The RSA increased in the following order: chitosanlaw successfully fitted the experimental data and GA diffusion coefficients, assumed to be Arrhenius temperature dependent were obtained for each film. This work will contribute to the development of environmentally friendly packaging films, resulting in improved food preservation and shelf-life extension.

Estudo da secagem do pedúnculo do caju em sistemas convencional e solar: modelagem e simulação do processo

This work aims to study the drying of cashew-nut pulp with different lay-out of dryers using conventional and solar energy. It concerns with the use of exceeding of the regional raw material and the suitable knowledge for the applicability of the drying systems as pathway for food conservation. Besides, it used renewable sources as solar energy to dry these agroindustrial products. Runs were carried out using a conventional tray-dryer with temperature, air velocity control and cashew slice thickness of 55°C, 65°C, 75°C; 3.0; 4.5, 6.0 m s-1; 1.0; 1.5 and 2.0 cm, respectively, in order to compare the studied systems. To evaluate the conventional tray-dryer, it was used a diffusional model of 2nd Fick´s law, where the drying curves were quite well fitted to an infinite flat plate design. For the drying runs where the room temperature had no control, it was developed a phenomenological-mathematical model for the solar dryer with indirect radiation under natural and forced convection based on material and energy balances of the system. Besides, it was carried out assays in the in natura as well as dehydrated, statistic analysis of the experimental drying data, sensorial analysis of the final dry product and a simplified economical analysis of the systems studied

EVALUATION OF DIFFUSION-COEFFICIENTS IN OSMOTIC CONCENTRATION OF BANANAS (MUSA-CAVENDISH LAMBERT)

The effect of time of exposure, solution concentration and temperature on the osmotic concentration of banana (slices of 11 mm thickness) was studied in aqueous sucrose solutions. The selectivity of the cellular tissues was reduced by steam blanching the banana slices before osmotic treatment. Effective diffusion coefficients for the loss of water and the increase in sucrose content were determined according to Fick's Law applied to a two-dimensional body; calculated on the basis of the concentration of various components in the liquid phase impenetrating the fruit. These coefficients revealed values similar to binary diffusion coefficients for pure sucrose solutions.

Influência de pré-tratamentos com emulsões a base de lecitina e óleo de soja na cinética de secagem da uva rubi

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Untersuchung der Filmbildung aus Polymerdispersionen mit Hilfe der forcierten Rayleighstreuung

Der Filmbildungsprozess wässriger Polymerdispersionen wurde mit forcierter Rayleighstreuung untersucht.Diffusionskoeffizienten D von Farbstoffsonden sind damit zwischen 10^-21 und 10^-9 m^2s^-1 zugänglich. Das Diffusionsverhalten der hydrophoben FarbstoffsondeAberchrome in feuchten und trockenen Filmen aus Poly(n-butylmethacrylat-co-acrylsäure)-Dispersionen sowie bei der Wiederbefeuchtung trockener Filme wurde untersucht.Die Dynamik von Aberchrome weicht in feuchten Filmen von Fickscher Diffusion ab. Dies äußert sich in Abweichungen vomcharakteristischen q^2-Verhalten der Relaxationsrate tau^-1 (tau^-1 = Dq^2; q:Streuvektor) und im Auftreten gestrecktexponentieller Intensitätsrelaxationskurven. Diese Anomalie wurde mit einem apparenten, längenskalenabhängigen Diffusionskoeffizienten Dapp(Lambda) (Lambda= 2Pi/q) beschrieben, der für Lambda -> 0 Werte annimmt, die einem homogen hydroplastifizierten Polymermaterial entsprechen, während Dapp(Lambda) für Lambda -> Unendlich stark anwächst. Diese Anomalien verschwinden bei Wassergehalten entsprechend der Polymerwasserlöslichkeit. Weiteres Trocknen führt zum Absinken des Fickschen-Diffusionskoeffizienten auf einen Grenzwert für trockene Filme. Die Ergebnisse konnten mit einem Zwei-Zustands-Modell beschrieben werden: Die Sonde diffundiert Ficksch in einer hydrophilen Grenzflächenphase und einer langsamen Polymerphase. Austausch zwischendiesen Phasen ist ohne Einschränkung möglich. Das Modell erlaubt die Quantifizierung des Einflusses des Trocknungsprozesses auf Polymer- und Grenzflächenphaseneigenschaften.Dies wurde durch systematische Veränderungen der Grenzflächeneigenschaften demonstriert. Dies geschah durch Acrylsäuregehaltvariationen in den Poly(n-butylmethacrylat-co-acrylsäure)-Dispersionen und Emulgatorbedeckungsgradvariationen. In beiden Fällen reflektieren sich Dispersionsveränderungen in einer Veränderung der Modellparameter in der Grenzflächenphase.

Diffusion-driven transport in clayrock formations

Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), Callovo-Oxfordian argillite (FR)). The ultimate safety of the corresponding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU). Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clay-solution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological time-space scales. Much effort was placed on improving understanding of coupled sorption-diffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materials for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a D-e value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics. (C) 2011 Elsevier Ltd. All rights reserved.

Ar diffusion in K-feldspar: present and absent

The mechanisms of Ar release from K-feldspar samples in laboratory experiments and during their geological history are assessed here. Modern petrology clearly established that the chemical and isotopic record of minerals is normally dominated by aqueous recrystallization. The laboratory critique is trickier, which explains why so many conflicting approaches have been able to survive long past their expiration date. Current models are evaluated for self-consistency; especially Arrhenian non-linearity leads to paradoxes. The models’ testable geological predictions suggest that temperature-based downslope extrapolations often overestimate observed geological Ar mobility substantially. An updated interpretation is based on the unrelatedness of geological behaviour to laboratory experiments. The isotopic record of K-feldspar in geological samples is not a unique function of temperature, as recrystallisation promoted by aqueous fluids is the predominant mechanism controlling isotope transport. K-feldspar should therefore be viewed as a hygrochronometer. Laboratory degassing proceeds from structural rearrangements and phase transitions such as are observed in situ at high temperature in Na and Pb feldspars. These effects violate the mathematics of an inert Fick’s Law matrix and preclude downslope extrapolation. The similar upward-concave, non-linear shapes of Arrhenius trajectories of many silicates, hydrous and anhydrous, are likely common manifestations of structural rearrangements in silicate structures.

A re-examination of petrogenesis and 40Ar/39Ar systematics in the the Chain of Ponds K-feldspar: "diffusion domain" archetype versus polyphase hygrochronology

K-feldspar (Kfs) from the Chain of Ponds Pluton (CPP) is the archetypal reference material, on which thermochronological modeling of Ar diffusion in discrete “domains” was founded. We re-examine the CPP Kfs using cathodoluminescence and back-scattered electron imaging, transmission electron microscopy, and electron probe microanalysis. 40Ar/39Ar stepwise heating experiments on different sieve fractions, and on handpicked and unpicked aliquots, are compared. Our results reproduce the staircase-shaped age spectrum and the Arrhenius trajectory of the literature sample, confirming that samples collected from the same locality have an identical Ar isotope record. Even the most pristine-looking Kfs from the CPP contains successive generations of secondary, metasomatic/retrograde mineral replacements that post-date magmatic crystallization. These chemically and chronologically distinct phases are responsible for its staircase-shaped age spectra, which are modified by handpicking. While genuine within-grain diffusion gradients are not ruled out by these data, this study demonstrates that the most important control on staircase-shaped age spectra is the simultaneous presence of heterochemical, diachronous post-magmatic mineral growth. At least five distinct mineral species were identified in the Kfs separate, three of which can be traced to external fluids interacting with the CPP in a chemically open system. Sieve fractions have size-shifted Arrhenius trajectories, negating the existence of the smallest “diffusion domains”. Heterochemical phases also play an important role in producing non-linear trajectories. In vacuo degassing rates recovered from Arrhenius plots are neither related to true Fick’s Law diffusion nor to the staircase shape of the age spectra. The CPP Kfs used to define the "diffusion domain" model demonstrates the predominance of metasomatic alteration by hydrothermal fluids and recrystallization in establishing the natural Ar distribution amongst different coexisting phases that gives rise to the staircase-shaped age spectrum. Microbeam imaging of textures is as essential for 40Ar-39Ar hygrochronology as it is for U-Pb geochronology.

Surface, Thermal and Antimicrobial Release Properties of Plasma-Treated Zein Films

The effects of dielectric barrier discharge plasma treatment on zein film containing thymol as an active ingredient were evaluated. The plasma discharge was optically characterized to identify the reactive species. A significant increase in the film roughness (p < 0.05) was observed due to the etching effect of DBD plasma, which was correlated with the increase in the diffusion rate of thymol in the food simulant. The diffusion of thymol from the zein film was measured in aqueous solution. The kinetics of thymol release followed the Fick’s law of diffusion as shown by the high correlation coefficients between experimental and theoretical data. No significant change (p > 0.05) was observed for the thermal properties of the antimicrobial films after DBD plasma treatment.

Procedure for calculating the chloride diffusion coefficient and surface concentration from a profile having a maximum beyond the concrete surface

Chlorides induce local corrosion in the steel reinforcements when reaching the bar surface. The measurement of the rate of ingress of these ions, is made by mathematically fitting the so called “error function equation” into the chloride concentration profile, obtaining so the diffusion coefficient and the chloride concentration at the concrete surface. However, the chloride profiles do not always follow Fick’s law by having the maximum concentration at the concrete surface, but often the profile shows a maximum concentration more in the interior, which indicates a different composition and performance of the most external concrete layer with respect to the internal zones. The paper presents a procedure prepared during the time of the RILEM TC 178-TMC: “Testing and modeling chloride penetration in concrete”, which suggests neglecting the external layer where the chloride concentration increases and using the maximum as an “apparent” surface concentration, called C max and to fit the error function equation into the decreasing concentration profile towards the interior. The prediction of evolution should be made also from the maximum.