Analysis of an asymptotic preserving scheme for linear kinetic equations in the diffusion limit

We present a mathematical analysis of the asymptotic preserving scheme proposed in [M. Lemou and L. Mieussens, SIAM J. Sci. Comput., 31 (2008), pp. 334-368] for linear transport equations in kinetic and diffusive regimes. We prove that the scheme is uniformly stable and accurate with respect to the mean free path of the particles. This property is satisfied under an explicitly given CFL condition. This condition tends to a parabolic CFL condition for small mean free paths and is close to a convection CFL condition for large mean free paths. Our analysis is based on very simple energy estimates. © 2010 Society for Industrial and Applied Mathematics.

Local kinetic interpretation of entropy production through reversed diffusion.

The time reversal of stochastic diffusion processes is revisited with emphasis on the physical meaning of the time-reversed drift and the noise prescription in the case of multiplicative noise. The local kinematics and mechanics of free diffusion are linked to the hydrodynamic description. These properties also provide an interpretation of the Pope-Ching formula for the steady-state probability density function along with a geometric interpretation of the fluctuation-dissipation relation. Finally, the statistics of the local entropy production rate of diffusion are discussed in the light of local diffusion properties, and a stochastic differential equation for entropy production is obtained using the Girsanov theorem for reversed diffusion. The results are illustrated for the Ornstein-Uhlenbeck process.

Quantitative mapping of trimethyltin injury in the rat brain using magnetic resonance histology.

The growing exposure to chemicals in our environment and the increasing concern over their impact on health have elevated the need for new methods for surveying the detrimental effects of these compounds. Today's gold standard for assessing the effects of toxicants on the brain is based on hematoxylin and eosin (H&E)-stained histology, sometimes accompanied by special stains or immunohistochemistry for neural processes and myelin. This approach is time-consuming and is usually limited to a fraction of the total brain volume. We demonstrate that magnetic resonance histology (MRH) can be used for quantitatively assessing the effects of central nervous system toxicants in rat models. We show that subtle and sparse changes to brain structure can be detected using magnetic resonance histology, and correspond to some of the locations in which lesions are found by traditional pathological examination. We report for the first time diffusion tensor image-based detection of changes in white matter regions, including fimbria and corpus callosum, in the brains of rats exposed to 8 mg/kg and 12 mg/kg trimethyltin. Besides detecting brain-wide changes, magnetic resonance histology provides a quantitative assessment of dose-dependent effects. These effects can be found in different magnetic resonance contrast mechanisms, providing multivariate biomarkers for the same spatial location. In this study, deformation-based morphometry detected areas where previous studies have detected cell loss, while voxel-wise analyses of diffusion tensor parameters revealed microstructural changes due to such things as cellular swelling, apoptosis, and inflammation. Magnetic resonance histology brings a valuable addition to pathology with the ability to generate brain-wide quantitative parametric maps for markers of toxic insults in the rodent brain.

Effects of chronic mild traumatic brain injury on white matter integrity in iraq and afghanistan war veterans

Mild traumatic brain injury (TBI) is a common source of morbidity from the wars in Iraq and Afghanistan. With no overt lesions on structural MRI, diagnosis of chronic mild TBI in military veterans relies on obtaining an accurate history and assessment of behavioral symptoms that are also associated with frequent comorbid disorders, particularly posttraumatic stress disorder (PTSD) and depression. Military veterans from Iraq and Afghanistan with mild TBI (n = 30) with comorbid PTSD and depression and non-TBI participants from primary (n = 42) and confirmatory (n = 28) control groups were assessed with high angular resolution diffusion imaging (HARDI). White matter-specific registration followed by whole-brain voxelwise analysis of crossing fibers provided separate partial volume fractions reflecting the integrity of primary fibers and secondary (crossing) fibers. Loss of white matter integrity in primary fibers (P < 0.05; corrected) was associated with chronic mild TBI in a widely distributed pattern of major fiber bundles and smaller peripheral tracts including the corpus callosum (genu, body, and splenium), forceps minor, forceps major, superior and posterior corona radiata, internal capsule, superior longitudinal fasciculus, and others. Distributed loss of white matter integrity correlated with duration of loss of consciousness and most notably with "feeling dazed or confused," but not diagnosis of PTSD or depressive symptoms. This widespread spatial extent of white matter damage has typically been reported in moderate to severe TBI. The diffuse loss of white matter integrity appears consistent with systemic mechanisms of damage shared by blast- and impact-related mild TBI that involves a cascade of inflammatory and neurochemical events. © 2012 Wiley Periodicals, Inc.

Correction for Eddy Current-Induced Echo-Shifting Effect in Partial-Fourier Diffusion Tensor Imaging.

In most diffusion tensor imaging (DTI) studies, images are acquired with either a partial-Fourier or a parallel partial-Fourier echo-planar imaging (EPI) sequence, in order to shorten the echo time and increase the signal-to-noise ratio (SNR). However, eddy currents induced by the diffusion-sensitizing gradients can often lead to a shift of the echo in k-space, resulting in three distinct types of artifacts in partial-Fourier DTI. Here, we present an improved DTI acquisition and reconstruction scheme, capable of generating high-quality and high-SNR DTI data without eddy current-induced artifacts. This new scheme consists of three components, respectively, addressing the three distinct types of artifacts. First, a k-space energy-anchored DTI sequence is designed to recover eddy current-induced signal loss (i.e., Type 1 artifact). Second, a multischeme partial-Fourier reconstruction is used to eliminate artificial signal elevation (i.e., Type 2 artifact) associated with the conventional partial-Fourier reconstruction. Third, a signal intensity correction is applied to remove artificial signal modulations due to eddy current-induced erroneous T2(∗) -weighting (i.e., Type 3 artifact). These systematic improvements will greatly increase the consistency and accuracy of DTI measurements, expanding the utility of DTI in translational applications where quantitative robustness is much needed.

Modelling meso-scale diffusion processes in stochastic fluid bio-membranes

The space–time dynamics of rigid inhomogeneities (inclusions) free to move in a randomly fluctuating fluid bio-membrane is derived and numerically simulated as a function of the membrane shape changes. Both vertically placed (embedded) inclusions and horizontally placed (surface) inclusions are considered. The energetics of the membrane, as a two-dimensional (2D) meso-scale continuum sheet, is described by the Canham–Helfrich Hamiltonian, with the membrane height function treated as a stochastic process. The diffusion parameter of this process acts as the link coupling the membrane shape fluctuations to the kinematics of the inclusions. The latter is described via Ito stochastic differential equation. In addition to stochastic forces, the inclusions also experience membrane-induced deterministic forces. Our aim is to simulate the diffusion-driven aggregation of inclusions and show how the external inclusions arrive at the sites of the embedded inclusions. The model has potential use in such emerging fields as designing a targeted drug delivery system.

A hybrid Laplace transform/finite difference boundary element method for diffusion problems

The solution process for diffusion problems usually involves the time development separately from the space solution. A finite difference algorithm in time requires a sequential time development in which all previous values must be determined prior to the current value. The Stehfest Laplace transform algorithm, however, allows time solutions without the knowledge of prior values. It is of interest to be able to develop a time-domain decomposition suitable for implementation in a parallel environment. One such possibility is to use the Laplace transform to develop coarse-grained solutions which act as the initial values for a set of fine-grained solutions. The independence of the Laplace transform solutions means that we do indeed have a time-domain decomposition process. Any suitable time solver can be used for the fine-grained solution. To illustrate the technique we shall use an Euler solver in time together with the dual reciprocity boundary element method for the space solution

Modeling the diffusion of solid copper into liquid solder alloys

During the soldering process, the copper atoms diffuse into liquid solders. The diffusion process determines integrity and the reworking possibility of a solder joint. In order to capture the diffusion scenarios of solid copper into liquid Sn–Pb and Sn–Cu solders, a computer modeling has been performed for 10 s. An analytical model has also been proposed for calculating the diffusion coefficient of copper into liquid solders. It is found that the diffusion coefficient for Sn–Pb solder is 2.74 × 10− 10 m2/s and for Sn–Cu solder is 6.44 × 10−9 m2/s. The modeling results reveal that the diffusion coefficient is one of the major factors that govern the rate at which solid Cu dissolve in the molten solder. The predicted dissolved amounts of copper into solders have been validated with the help of scanning electron microscopic analysis.

Local examination of skin diffusion using FTIR spectroscopic imaging and multivariate target factor analysis

In the context of trans-dermal drug delivery it is very important to have mechanistic insight into the barrier function of the skin's stratum corneum and the diffusion mechanisms of topically applied drugs. Currently spectroscopic imaging techniques are evolving which enable a spatial examination of various types of samples in a dynamic way. ATR-FTIR imaging opens up the possibility to monitor spatial diffusion profiles across the stratum corneum of a skin sample. Multivariate data analyses methods based on factor analysis are able to provide insight into the large amount of spectroscopically complex and highly overlapping signals generated. Multivariate target factor analysis was used for spectral resolution and local diffusion profiles with time through stratum corneum. A model drug, 4-cyanophenol in polyethylene glycol 600 and water was studied. Results indicate that the average diffusion profiles between spatially different locations show similar profiles despite the heterogeneous nature of the biological sample and the challenging experimental set-up.

Simultaneous monitoring of drug and solvent diffusion across a model membrane using ATR-FTIR spectroscopy

Attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy has been used to simultaneously follow the diffusion of model drugs and solvent across polydimethylsiloxane (silicone) membrane. Three model drugs, cyanophenol (CNP), methyl nicotinate (MN) and butyl paraben (BP) were selected to cover a range of lipophilicities. Isostearyl isostearate (ISIS) was chosen as the solvent because its large molecular weight should facilitate observation of whether the drug molecules are able to diffuse through the membrane independently of the solvent. The diffusion of the three drugs and the solvent was successfully described by a Fickian model. The effects of parameters such as the absorption wavelength used to follow diffusion on the calculated diffusion coefficient were investigated. Absorption wavelength which affects the depth of penetration of the infrared radiation into the membrane did not significantly affect the calculated diffusion coefficient over the wavelength range tested. Each of the model drugs was observed to diffuse independently of the solvent across the membrane. The diffusion of a CNP-ISIS hydrogen bonded complex across the membrane was also monitored. The relative diffusion rates of the solute and solvent across the membrane can largely be accounted for by the molecular size of the permeant.

ATR-FTIR spectroscopy and spectroscopic imaging of solvent and permeant diffusion across model membranes

The uptake and diffusion of solvents across polymer membranes is important in controlled drug delivery, effects on drug uptake into, for example, infusion bags and containers, as well as transport across protective clothing. Attenuated Total Reflectance Fourier Transform Infrared (ATR-FTIR) spectroscopy has been used to monitor the effects of different solvents on the diffusion of a model compound, 4-cyanophenol (CNP) across silicone membrane and on the equilibrium concentration of CNP obtained in the membrane following diffusion. ATR-FTIR spectroscopic imaging of membrane diffusion was used to gain an understanding of when the boundary conditions applied to Fick's second law, used to model the diffusion of permeants across the silicone membrane do not hold. The imaging experiments indicated that when the solvent was not taken up appreciably into the membrane, the presence of discrete solvent pools between the ATR crystal and the silicone membrane can affect the diffusion profile of the permeant. This effect is more significant if the permeant has a high solubility in the solvent. In contrast, solvents that are taken up into the membrane to a greater extent, or those where the solubility of the permeant in the vehicle is relatively low, were found to show a good fit to the diffusion model. As such these systems allow the ATR-FTIR spectroscopic approach to give mechanistic insight into how the particular solvents enhance permeation. The solubility of CNP in the solvent and the uptake of the solvent into the membrane were found to be important influences on the equilibrium concentration of the permeant obtained in the membrane following diffusion. In general, solvents which were taken up to a significant extent into the membrane and which caused the membrane to swell increased the diffusion coefficient of the permeant in the membrane though other factors such as solvent viscosity may also be important.