Stochastic Modeling of Reversible Biochemical Reaction-Diffusion Systems and High-Resolution Shock-Capturing Methods for Fluid Interfaces

Thesis (Ph.D.)--University of Washington, 2016-08

Diffusion de la lumière dans les milieux complexes à base de cristaux liquides

Cette thèse présente les étapes de fabrication et les résultats de la caractérisation des modèles de tissu biologique fabriqués en utilisant des cristaux liquides et des polymères qui permettent de contrôler la diffusion et la transmission de la lumière par un champ électrique. Ce champ électrique entraîne un changement d’orientation des molécules de cristaux liquides, ce qui induit des variations locales de l’indice de réfraction dans ces matériaux. Nous avons utilisé ce contrôle de l’indice de réfraction pour contrôler le niveau de diffusion et de transmission dans différents types d’échantillons. Les échantillons utilisés sont faits a) de cristaux liquides purs, b) de sphères de polymère dans un environnement de cristaux liquides, c) de cristaux liquides ou de gouttelettes de cristaux liquides dans un environnement de polymère. Les travaux réalisés nous ont permis de proposer une méthode électro-optique pour simuler, à l’aide de cristaux liquides, des milieux diffusants tels que les tissus biologiques. Dans la recherche de modèles, nous avons montré qu’il est possible de contrôler, par la fréquence du champ électrique, des gouttelettes de cristaux liquides dispersées dans une matrice de polymère contenant des ions lithium. Une sensibilité sélective à la fréquence de la tension électrique appliquée en fonction de la taille des gouttelettes est observée. Par la suite, nous avons démontré l’effet de « quenching » interférentiel de la diffusion pour une sonde optique cohérente en étudiant des modèles non alignés (clustérisés) de cristaux liquides. Nous avons aussi démontré des applications potentielles de ce phénomène tel le filtrage en longueur d’onde par l’ajustement d’un champ électrique. Nous avons ensuite démontré que l’utilisation des cristaux liquides - dits « dual-frequency » - permet une modulation rapide des états ou des structures du matériau ainsi que de la transmission et de la diffusion de lumière de ce dernier. Enfin, nous avons aussi découvert que ces modèles peuvent être contrôlés pour retrouver le degré de polarisation de la lumière diffusée.

Formulation and analysis of a diffusion-velocity particle model for transport-dispersion equations

The modelling of diffusive terms in particle methods is a delicate matter and several models were proposed in the literature to take such terms into account. The diffusion velocity method (DVM), originally designed for the diffusion of passive scalars, turns diffusive terms into convective ones by expressing them as a divergence involving a so-called diffusion velocity. In this paper, DVM is extended to the diffusion of vectorial quantities in the three-dimensional Navier–Stokes equations, in their incompressible, velocity–vorticity formulation. The integration of a large eddy simulation (LES) turbulence model is investigated and a DVM general formulation is proposed. Either with or without LES, a novel expression of the diffusion velocity is derived, which makes it easier to approximate and which highlights the analogy with the original formulation for scalar transport. From this statement, DVM is then analysed in one dimension, both analytically and numerically on test cases to point out its good behaviour.

Development of computer application for determination of chemical diffusion properties in biological tissues

The main objective of this work was to develop an application capable of determining the diffusion times and diffusion coefficients of optical clearing agents and water inside a known type of muscle. Different types of chemical agents can also be used with the method implemented, such as medications or metabolic products. Since the diffusion times can be calculated, it is possible to describe the dehydration mechanism that occurs in the muscle. The calculation of the diffusion time of an optical clearing agent allows to characterize the refractive index matching mechanism of optical clearing. By using both the diffusion times and diffusion of water and clearing agents not only the optical clearing mechanisms are characterized, but also information about optical clearing effect duration and magnitude is obtained. Such information is crucial to plan a clinical intervention in cooperation with optical clearing. The experimental method and equations implemented in the developed application are described in throughout this document, demonstrating its effectiveness. The application was developed in MATLAB code, but the method was personalized so it better fits the application needs. This process significantly improved the processing efficiency, reduced the time to obtain he results, multiple validations prevents common errors and some extra functionalities were added such as saving application progress or export information in different formats. Tests were made using glucose measurements in muscle. Some of the data, for testing purposes, was also intentionally changed in order to obtain different simulations and results from the application. The entire project was validated by comparing the calculated results with the ones found in literature, which are also described in this document.

A survey on prevalence rate of aflatoxigenic species of Aspergillus and residues of aflatoxins by ELISA method in cold-water cultural fish (rainbow trout) feeds in Tehran and West Azarbayjan provinces-Iran

Aflatoxins are one kind of fungal toxins produced by species of toxigenic Aspergillus (A. flavus and A. parasiticus) and in other words they are secondary metabolites which are considered as one of the threatening factors of food consumer's health. In this research 96 samples of cold-water cultural fish feed, rainbow trout, during the seasons of spring and summer of 2007 (every fifteenth of the month) were randomized (by simple and stratified random) to determine: 1. The prevalence rate of aflatoxigenic species of Aspergillus in stored feed of cold-water cultural fish in West Azarbayjan cultural fish farms in both seasons (spring and summer); 2. The residues of total aflatoxin in stored feed of fish in cultural fish farms of West Azarbayjan in both seasons by ELISA method; and 3. The residues of that toxin in feed produced in aquatic feed factories in Tehran and West Azarbyjan provinces with the same method. In order to study prevalence rate of toxigenic species of Aspergillus, pour-plate culture method by general medium such as Malt Extract Agar (M.E.A.) and Sabouraud-Dextrose Agar (S.D.A.) and by standard No.997 of Iranian Standard Institute were used. The produced colonies were examined microscopically. To determine the aflatoxins residues, ELISA method using Agra-Quant kit of Romer Lab company, were applied. The results of this survey indicated that only 8.3% of the samples were infected by A. flavus. A. parasiticus was not observed. There were no significant differences between the prevalence rate of AFT and seasons/months, either (P<0.05). Evaluating mean of aflatoxin rate showed that the rates of this variable are lower than the tolerance levels designated by the joint FAO/WHO expert committee (The mean of AFT in all data was lower than 11 ppb). Furthermore, mean of total AFT residues rates of stored feed of various cultural center of West Azarbayjan and Tehran factories were comparable in spring and summer, and no significant differences were observed (P<0.05). But there were significant differences between the total aflatoxin rates in the feed of West Azarbayjan factory and spring and summer (P<0.05), and AFT residues in spring (8.6 ppb) were higher than summer (6.1 ppb). Prevalence rates of AFT in Tehran feed factories (9.2 ppb) are higher than W. Azarbayjan (7.4 ppb). In other words, location was considered as a decisive factor in total AFT rates of samples. Moreover, the results indicated that there was significant difference between total aflatoxin rates of feed and cultural centers (P<0.05). The mean of AFT rates in embankment dam cultural fish farms (6.75 ppb) and multi-functions cultural fish farms (6.25 ppb) was higher than individual cultural pond (4.67 ppb). In conclusion, the finally results of this survey indicated that the lower rates of Aspergillus is not effective on the presence of total aflatoxin rates in trout feed. Due to low levels of aflatoxin rates (lower than 20 ppb), the produced feed of cold-cultural fishes, Rainbow Trout, in Tehran and West Azarbayjan provinces, in spring and summer of 2007, were safe and healthy both for fish and their consumers.

Fluorinated surfactants in solution: Diffusion coefficients of fluorinated alcohols in water

Intra-diffusion coefficients of three fluorinated alcohols, 2,2,3,3,3-pentafluoropropan-1-ol (PFP), 2,2,3,3,4,4,4-heptafluorobutan-1-ol (HFB) and 2,2,3,3,4,4,5,5,5-nonafluoropentan-1-ol (NFP) in water have been measured by the PFG–NMR spin-echo technique as a function of temperature and composition, focusing on the alcohol dilute region. For comparison, intra-diffusion coefficients of 2,2,2- trifluoroethanol (TFE) and HFB have also been measured in heavy water using the same method and conditions. As far as we know, these are the first experimental measurements of this property for these binary systems. Intra-diffusion coefficients for NFP in water and for TFE and HFB in heavy water have also been obtained by molecular dynamics simulation, complementing those for TFE, PFP and HFB reported in a previous work. The agreement between experimental and simulated results for PFP, HFB and NFP in water is reasonable, although presenting higher deviations than for the TFE/water system. From the dependence of the intra-diffusion coefficients on temperature, diffusion activation energies were estimated for all the solutes in water and heavy water.

Asymptotic stability of a coupled Advection-Diffusion-Reaction system arising in bioreactor processes.

In this work, we perform an asymptotic analysis of a coupled system of two Advection-Diffusion-Reaction equations with Danckwerts boundary conditions, which models the interaction between a microbial population (e.g., bacterias), called biomass, and a diluted organic contaminant (e.g., nitrates), called substrate, in a continuous flow bioreactor. This system exhibits, under suitable conditions, two stable equilibrium states: one steady state in which the biomass becomes extinct and no reaction is produced, called washout, and another steady state, which corresponds to the partial elimination of the substrate. We use the method of linearization to give sufficient conditions for the asymptotic stability of the two stable equilibrium configurations. Finally, we compare our asymptotic analysis with the usual asymptotic analysis associated to the continuous bioreactor when it is modeled with ordinary differential equations.

Adaptive coatings based on polyaniline for direct 2D observation of diffusion processes in microfluidic systems

[EN] Therefore the understanding and proper evaluation of the flow and mixing behaviour at microscale becomes a very important issue. In this study, the diffusion behaviour of two reacting solutions of HCI and NaOH were directly observed in a glass/polydimethylsiloxane microfluidic device using adaptive coatings based on the conductive polymer polyaniline that are covalently attached to the microchannel walls. The two liquid streams were combined at the junction of a Y-shaped microchannel, and allowed to diffuse into each other and react. The results showed excellent correlation between optical observation of the diffusion process and the numerical results. A numerical model which is based on finite volume method (FVM) discretisation of steady Navier-Stokes (fluid flow) equations and mass transport equations without reactions was used to calculate the flow variables at discrete points in the finite volume mesh element. The high correlation between theory and practical data indicates the potential of such coatings to monitor diffusion processes and mixing behaviour inside microfluidic channels in a dye free environment.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.

An approximate method for the solution of an influence function foil rolling model

Fleck and Johnson (Int. J. Mech. Sci. 29 (1987) 507) and Fleck et al. (Proc. Inst. Mech. Eng. 206 (1992) 119) have developed foil rolling models which allow for large deformations in the roll profile, including the possibility that the rolls flatten completely. However, these models require computationally expensive iterative solution techniques. A new approach to the approximate solution of the Fleck et al. (1992) Influence Function Model has been developed using both analytic and approximation techniques. The numerical difficulties arising from solving an integral equation in the flattened region have been reduced by applying an Inverse Hilbert Transform to get an analytic expression for the pressure. The method described in this paper is applicable to cases where there is or there is not a flat region.