Movement effects on the flow physics and nutrient delivery in engineered valvular tissues

Mechanical conditioning has been shown to promote tissue formation in a wide variety of tissue engineering efforts. However the underlying mechanisms by which external mechanical stimuli regulate cells and tissues are not known. This is particularly relevant in the area of heart valve tissue engineering (HVTE) owing to the intense hemodynamic environments that surround native valves. Some studies suggest that oscillatory shear stress (OSS) caused by steady flow and scaffold flexure play a critical role in engineered tissue formation derived from bone marrow derived stem cells (BMSCs). In addition, scaffold flexure may enhance nutrient (e.g. oxygen, glucose) transport. In this study, we computationally quantified the i) magnitude of fluid-induced shear stresses; ii) the extent of temporal fluid oscillations in the flow field using the oscillatory shear index (OSI) parameter, and iii) glucose and oxygen mass transport profiles. Noting that sample cyclic flexure induces a high degree of oscillatory shear stress (OSS), we incorporated moving boundary computational fluid dynamic simulations of samples housed within a bioreactor to consider the effects of: 1) no flow, no flexure (control group), 2) steady flow-alone, 3) cyclic flexure-alone and 4) combined steady flow and cyclic flexure environments. We also coupled a diffusion and convention mass transport equation to the simulated system. We found that the coexistence of both OSS and appreciable shear stress magnitudes, described by the newly introduced parameter OSI-:τ: explained the high levels of engineered collagen previously observed from combining cyclic flexure and steady flow states. On the other hand, each of these metrics on its own showed no association. This finding suggests that cyclic flexure and steady flow synergistically promote engineered heart valve tissue production via OSS, so long as the oscillations are accompanied by a critical magnitude of shear stress. In addition, our simulations showed that mass transport of glucose and oxygen is enhanced by sample movement at low sample porosities, but did not play a role in highly porous scaffolds. Preliminary in-house in vitro experiments showed that cell proliferation and phenotype is enhanced in OSI-:τ: environments.^

Mathematical equation correction to spectral and transport interferences in high-resolution continuum source flame atomic absorption spectrometry: determination of lead in phosphoric acid

In this work, a new mathematical equation correction approach for overcoming spectral and transport interferences was proposed. The proposal was applied to eliminate spectral interference caused by PO molecules at the 217.0005 nm Pb line, and the transport interference caused by variations in phosphoric acid concentrations. Correction may be necessary at 217.0005 nm to account for the contribution of PO, since Atotal217.0005 nm = A Pb217.0005 nm + A PO217.0005 nm. This may be easily done by measuring other PO wavelengths (e.g. 217.0458 nm) and calculating the relative contribution of PO absorbance (A PO) to the total absorbance (Atotal) at 217.0005 nm: A Pb217.0005 nm = Atotal217.0005 nm - A PO217.0005 nm = Atotal217.0005 nm - k (A PO217.0458 nm). The correction factor k is calculated from slopes of calibration curves built up for phosphorous (P) standard solutions measured at 217.0005 and 217.0458 nm, i.e. k = (slope217.0005 nm/slope217.0458 nm). For wavelength integrated absorbance of 3 pixels, sample aspiration rate of 5.0 ml min-1, analytical curves in the 0.1 - 1.0 mg L-1 Pb range with linearity better than 0.9990 were consistently obtained. Calibration curves for P at 217.0005 and 217.0458 nm with linearity better than 0.998 were obtained. Relative standard deviations (RSD) of measurements (n = 12) in the range of 1.4 - 4.3% and 2.0 - 6.0% without and with mathematical equation correction approach were obtained respectively. The limit of detection calculated to analytical line at 217.0005 nm was 10 µg L-1 Pb. Recoveries for Pb spikes were in the 97.5 - 100% and 105 - 230% intervals with and without mathematical equation correction approach, respectively.

Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tesselations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique.

Transport properties and equation of state of quark-nuclear matter

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Mathematical equation correction to spectral and transport interferences in high-resolution continuum source flame atomic absorption spectrometry: determination of lead in phosphoric acid

Um método de correção de interferência espectral e de transporte é proposto, e foi aplicado para minimizar interferências por moléculas de PO produzidas em chama ar-acetileno e de transporte causada pela variação da concentração de ácido fosfórico. Átomos de Pb e moléculas de PO absorvem a 217,0005 nm, então Atotal217,0005 nm = A Pb217,0005 nm + A PO217,0005 nm. Monitorando o comprimento de onda alternativo de PO em 217,0458 nm, é possível calcular a contribuição relativa de PO na absorbância total a 217,0005 nm: A Pb217,0005 nm = Atotal217,0005 nm - A PO217,0005 nm = Atotal217,0005 nm - k (A PO217,0458 nm). O fator de correção k é a razão entre os coeficientes angulares de duas curvas analíticas para P obtidas a 217,0005 e 217,0458 nm (k = b217,0005 nm/b217,0458 nm). Fixando-se a taxa de aspiração da amostra em 5,0 ml min-1, e integrando-se a absorbância no comprimento de onda a 3 pixels, curvas analíticas para Pb (0,1 - 1,0 mg L-1) foram obtidas com coeficientes de correlação típicos > 0,9990. As correlações lineares entre absorbância e concentração de P nos comprimentos de onda 217,0005 e 217,0458 foram > 0,998. O limite de detecção de Pb foi 10 µg L-1. O método de correção proposto forneceu desvios padrão relativos (n=12) de 2,0 a 6,0%, ligeiramente menores que os obtidos sem correção (1,4-4,3%). As recuperações de Pb adicionado às amostras de ácido fosfórico variaram de 97,5 a 100% (com correção pelo método proposto) e de 105 a 230% (sem correção).

The role of nuclear physics in supernovae and the evolution of neutron stars, Neutrino Opacities, Equation of State, Transport Coefficients, and Dark Matter Production

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

A study of electrochemical transport and diffuse charge dynamics using a Langevin equation

This thesis aims to develop new numerical and computational tools to study electrochemical transport and diffuse charge dynamics at small scales. Previous efforts at modeling electrokinetic phenomena at scales where the noncontinuum effects become significant have included continuum models based on the Poisson-Nernst-Planck equations and atomic simulations using molecular dynamics algorithms. Neither of them is easy to use or conducive to electrokinetic transport modeling in strong confinement or over long time scales. This work introduces a new approach based on a Langevin equation for diffuse charge dynamics in nanofluidic devices, which incorporates features from both continuum and atomistic methods. The model is then extended to include steric effects resulting from finite ion size, and applied to the phenomenon of double layer charging in a symmetric binary electrolyte between parallel-plate blocking electrodes, between which a voltage is applied. Finally, the results of this approach are compared to those of the continuum model based on the Poisson-Nernst-Planck equations.

Nonlinear transport and oscillating magnetoresistance in double quantum wells

The nonlinear regime of low-temperature magnetoresistance of double quantum wells in the region of magnetic fields below 1 T is studied both experimentally and theoretically. The observed inversion of the magnetointersubband oscillation peaks with increasing electric current and splitting of these peaks are described by the theory based on the kinetic equation for the isotropic nonequilibrium part of electron distribution function. The inelastic-scattering time of electrons is determined from the current dependence of the inversion field.

Solute transport by surface runoff from low-angle slopes: theory and application

The removal of chemicals in solution by overland how from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface-mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff: to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright (C) 2000 John Wiley & Sons, Ltd.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Periodic electric field enhanced transport through membranes

We examine the mean flux across a homogeneous membrane of a charged tracer subject to an alternating, symmetric voltage waveform. The analysis is based on the Nernst-Planck flux equation, with electric field subject to time dependence only. For low frequency electric fields the quasi steady-state flux can be approximated using the Goldman model, which has exact analytical solutions for tracer concentration and flux. No such closed form solutions can be found for arbitrary frequencies, however we find approximations for high frequency. An approximation formula for the average flux at all frequencies is also obtained from the two limiting approximations. Numerical integration of the governing equation is accomplished by use of the numerical method of lines and is performed for four different voltage waveforms. For the different voltage profiles, comparisons are made with the approximate analytical solutions which demonstrates their applicability. (c) 2005 Elsevier B.V. All rights reserved.

Proton transport model in the ionosphere .1. Multistream approach of the transport equations

The suprathermal particles, electrons and protons, coming from the magnetosphere and precipitating into the high-latitude atmosphere are an energy source of the Earth's ionosphere. They interact with ambient thermal gas through inelastic and elastic collisions. The physical quantities perturbed by these precipitations, such as the heating rate, the electron production rate, or the emission intensities, can be provided in solving the kinetic stationary Boltzmann equation. This equation yields particle fluxes as a function of altitude, energy, and pitch angle. While this equation has been solved through different ways for the electron transport and fully tested, the proton transport is more complicated. Because of charge-changing reactions, the latter is a set of two-coupled transport equations that must be solved: one for protons and the other for H atoms. We present here a new approach that solves the multistream proton/hydrogen transport equations encompassing the collision angular redistributions and the magnetic mirroring effect. In order to validate our model we discuss the energy conservation and we compare with another model under the same inputs and with rocket observations. The influence of the angular redistributions is discussed in a forthcoming paper.

Contractive metrics for a Boltzmann equation for granular gases: diffusive equilibriaThe Patterson-Sullivan embedding and minimal volume entropy for outer space

We quantify the long-time behavior of a system of (partially) inelastic particles in a stochastic thermostat by means of the contractivity of a suitable metric in the set of probability measures. Existence, uniqueness, boundedness of moments and regularity of a steady state are derived from this basic property. The solutions of the kinetic model are proved to converge exponentially as t→ ∞ to this diffusive equilibrium in this distance metrizing the weak convergence of measures. Then, we prove a uniform bound in time on Sobolev norms of the solution, provided the initial data has a finite norm in the corresponding Sobolev space. These results are then combined, using interpolation inequalities, to obtain exponential convergence to the diffusive equilibrium in the strong L¹-norm, as well as various Sobolev norms.

Convergence of the mass-transport steepest descent scheme for the sub-critical Patlak-Keller-Segel model

Variational steepest descent approximation schemes for the modified Patlak-Keller-Segel equation with a logarithmic interaction kernel in any dimension are considered. We prove the convergence of the suitably interpolated in time implicit Euler scheme, defined in terms of the Euclidean Wasserstein distance, associated to this equation for sub-critical masses. As a consequence, we recover the recent result about the global in time existence of weak-solutions to the modified Patlak-Keller-Segel equation for the logarithmic interaction kernel in any dimension in the sub-critical case. Moreover, we show how this method performs numerically in one dimension. In this particular case, this numerical scheme corresponds to a standard implicit Euler method for the pseudo-inverse of the cumulative distribution function. We demonstrate its capabilities to reproduce easily without the need of mesh-refinement the blow-up of solutions for super-critical masses.