Fractional pennes' bioheat equation: theoretical and numerical studies

In this work we provide a new mathematical model for the Pennes’ bioheat equation, assuming a fractional time derivative of single order. Alternative versions of the bioheat equation are studied and discussed, to take into account the temperature-dependent variability in the tissue perfusion, and both finite and infinite speed of heat propagation. The proposed bioheat model is solved numerically using an implicit finite difference scheme that we prove to be convergent and stable. The numerical method proposed can be applied to general reaction diffusion equations, with a variable diffusion coefficient. The results obtained with the single order fractional model, are compared with the original models that use classical derivatives.

Comparison of different numerical methods for the solution of the time-fractional reaction-diffusion equation with variable diffusion coefficient

In this work we perform a comparison of two different numerical schemes for the solution of the time-fractional diffusion equation with variable diffusion coefficient and a nonlinear source term. The two methods are the implicit numerical scheme presented in [M.L. Morgado, M. Rebelo, Numerical approximation of distributed order reaction- diffusion equations, Journal of Computational and Applied Mathematics 275 (2015) 216-227] that is adapted to our type of equation, and a colocation method where Chebyshev polynomials are used to reduce the fractional differential equation to a system of ordinary differential equations

Fractional bioheat equation

In this work we develop a new mathematical model for the Pennes’ bioheat equation assuming a fractional time derivative of single order. A numerical method for the solu- tion of such equations is proposed, and, the suitability of the new model for modelling real physical problems is studied and discussed

Nonsymmetrical separations in thermal diffusion columns

The influence of the feed composition upon the actual degrees of separation attained at the top and bottom sections of a thermogravitational column is discussed using the classical phenomenological theory of Furry, Jones, and Onsager. It is shown that, except for a feed composition of C 0 = 0.5 (mass fraction), the separation profile is nonsymmetric, i.e., the separations in the top and bottom sections of the column are nonsymmetric with respect to the feed composition, the asymmetry increasing as the feed composition moves away from C 0 = 0.5. An equation for the determination of the optimum feed location as a function of the feed composition is derived.

The application of the equivalent annulus width concept in thermal diffusion separations

The equivalent annulus width concept is used to characterize a small commercial thermogravitational hermal diffusion column and its validity checked experimentally by separating batchwise in the column mixtures of n-heptanebenzene with different initial concentrations. The equation of Ruppell and Coull was used to analyse the data in the short separation times range and determine the equivalent annulus width. Good agreement was obtained between the experimental and predicted time-separation curves when using the equivalent annulus width value and on averaged value of the thermal diffusion constant. A new method is presented for the simultaneous determination of the equivalent annulus width and the thermal diffusion constant of a binary mixture from a single set of experimental data.

Masonry infill walls under blast loading using confined underwater blast wave generators (WBWG)

The vulnerability of the masonry envelop under blast loading is considered critical due to the risk of loss of lives. The behaviour of masonry infill walls subjected to dynamic out-of-plane loading was experimentally investigated in this work. Using confined underwater blast wave generators (WBWG), applying the extremely high rate conversion of the explosive detonation energy into the kinetic energy of a thick water confinement, allowed a surface area distribution avoiding also the generation of high velocity fragments and reducing atmospheric sound wave. In the present study, water plastic containers, having in its centre a detonator inside a cylindrical explosive charge, were used in unreinforced masonry infills panels with 1.7m by 3.5m. Besides the usage of pressure and displacement transducers, pictures with high-speed video cameras were recorded to enable processing of the deflections and identification of failure modes. Additional numerical studies were performed in both unreinforced and reinforced walls. Bed joint reinforcement and grid reinforcement were used to strengthen the infill walls, and the results are presented and compared, allowing to obtain pressure-impulse diagrams for design of masonry infill walls.

Simulation of humidity fields in concrete: experimental validation and parameter estimation

The moisture content in concrete structures has an important influence in their behavior and performance. Several vali-dated numerical approaches adopt the governing equation for relative humidity fields proposed in Model Code 1990/2010. Nevertheless there is no integrative study which addresses the choice of parameters for the simulation of the humidity diffusion phenomenon, particularly in concern to the range of parameters forwarded by Model Code 1990/2010. A software based on a Finite Difference Method Algorithm (1D and axisymmetric cases) is used to perform sensitivity analyses on the main parameters in a normal strength concrete. Then, based on the conclusions of the sensi-tivity analyses, experimental results from nine different concrete compositions are analyzed. The software is used to identify the main material parameters that better fit the experimental data. In general, the model was able to satisfactory fit the experimental results and new correlations were proposed, particularly focusing on the boundary transfer coeffi-cient.

Assessing the application of polypropylene multifilament yarns in brain phantoms

The currently available clinical imaging methods do not provide highly detailed information about location and severity of axonal injury or the expected recovery time of patients with traumatic brain injury [1]. High-Definition Fiber Tractography (HDFT) is a novel imaging modality that allows visualizing and quantifying, directly, the degree of axons damage, predicting functional deficits due to traumatic axonal injury and loss of cortical projections. This imaging modality is based on diffusion technology [2]. The inexistence of a phantom able to mimic properly the human brain hinders the possibility of testing, calibrating and validating these medical imaging techniques. Most research done in this area fails in key points, such as the size limit reproduced of the brain fibers and the quick and easy reproducibility of phantoms [3]. For that reason, it is necessary to develop similar structures matching the micron scale of axon tubes. Flexible textiles can play an important role since they allow producing controlled packing densities and crossing structures that match closely the human crossing patterns of the brain. To build a brain phantom, several parameters must be taken into account in what concerns to the materials selection, like hydrophobicity, density and fiber diameter, since these factors influence directly the values of fractional anisotropy. Fiber cross-section shape is other important parameter. Earlier studies showed that synthetic fibrous materials are a good choice for building a brain phantom [4]. The present work is integrated in a broader project that aims to develop a brain phantom made by fibrous materials to validate and calibrate HDFT. Due to the similarity between thousands of hollow multifilaments in a fibrous arrangement, like a yarn, and the axons, low twist polypropylene multifilament yarns were selected for this development. In this sense, extruded hollow filaments were analysed in scanning electron microscope to characterize their main dimensions and shape. In order to approximate the dimensional scale to human axons, five types of polypropylene yarns with different linear density (denier) were used, aiming to understand the effect of linear density on the filament inner and outer areas. Moreover, in order to achieve the required dimensions, the polypropylene filaments cross-section was diminished in a drawing stage of a filament extrusion line. Subsequently, tensile tests were performed to characterize the mechanical behaviour of hollow filaments and to evaluate the differences between stretched and non-stretched filaments. In general, an increase of the linear density causes the increase in the size of the filament cross section. With the increase of structure orientation of filaments, induced by stretching, breaking tenacity increases and elongation at break decreases. The production of hollow fibers, with the required characteristics, is one of the key steps to create a brain phantom that properly mimics the human brain that may be used for the validation and calibration of HDFT, an imaging approach that is expected to contribute significantly to the areas of brain related research.

Inferior frontal gyrus white matter abnormalities in obsessive-compulsive disorder

The aim of the present study is to explore obsessive-compulsive disorder (OCD)-related abnormalities in white matter connectivity in OCD for a core region associated with inhibitory control [i.e. inferior frontal gyrus (IFG)]. Fifteen patients with OCD (11 men) and 15 healthy controls (nine men) underwent diffusion tensor imaging scanning to study four diffusivity indexes of white matter integrity [fractional anisotropy, mean diffusivity (MD), axial diffusivity and radial diffusivity (RD)]. The results showed that persons with OCD manifested significantly lower fractional anisotropy levels in the bilateral IFG as well as its parcellations in the pars opercularis, pars triangularis, and pars orbitalis. Significantly higher levels of MD, RD were evident for the OCD group in the IFG as a whole as well as in the bilateral subregions of the pars triangularis and pars opercularis (for MD and RD), the right side of the pars orbitalis (for RD), and the left side of the pars triangularis and right side pars opercularis (for axial diffusivity). Overall, the results suggest significant alterations in structural connectivity, probably associated with myelination and axonal abnormalities in the IFG of OCD patients.

Pulse wave velocity distribution in a cohort study: from arterial stiffness to early vascular aging

BACKGROUND: By contrast with other southern European people, north Portuguese population registers an especially high prevalence of hypertension and stroke incidence. We designed a cohort study to identify individuals presenting accelerated and premature arterial aging in the Portuguese population. METHOD: Pulse wave velocity (PWV) was measured in randomly sampled population dwellers aged 18-96 years from northern Portugal, and used as a marker of early vascular aging (EVA). Of the 3038 individuals enrolled, 2542 completed the evaluation. RESULTS: Mean PWV value for the entire population was 8.4?m/s (men: 8.6?m/s; women: 8.2?m/s; P??10?m/s). Logistic regression models indicated gender differences concerning the risk of developing large artery damage, with women having the same odds of PWV above 10?m/s 10 years later than men. CONCLUSION: The population PWV values were higher than expected in a low cardiovascular risk area (Portugal). High prevalence rates of EVA and noteworthy large artery damage in young ages were found.

Separation by thermal diffusion in a rotary column

A rotary thermal diffusion column with the inner cylinder rotating and the outer cylinder static was used to separate n-heptane-benzene mixtures at different speeds of rotation. The results show that the column efficiency depends on the speed of rotation. For the optimum speed the increase in efficiency relative to the static column was of the order of 8%. The role of the geometric irregularities in the annulus width on performance of the rotary column is also discussed.

Some observations on asymmetrical separation in thermal diffusion columns

Data have been obtained in steady-state batch operated thermogravitational separation columns using different binary mixtures to test the theory recently published by Morgado et al. The experimental results confirm that separations by thermal diffusion are asymmetrical except when the initial concentration is 0.5 and that the asymmetry is larger as the initial concentration deviates from 0.5 and as the separation potential increases.

Single molecule simulation of diffusion and enzyme kinetics

This work presents a molecular-scale agent-based model for the simulation of enzymatic reactions at experimentally measured concentrations. The model incorporates stochasticity and spatial dependence, using diffusing and reacting particles with physical dimensions. We developed strategies to adjust and validate the enzymatic rates and diffusion coefficients to the information required by the computational agents, i.e., collision efficiency, interaction logic between agents, the time scale associated with interactions (e.g., kinetics), and agent velocity. Also, we tested the impact of molecular location (a source of biological noise) in the speed at which the reactions take place. Simulations were conducted for experimental data on the 2-hydroxymuconate tautomerase (EC 5.3.2.6, UniProt ID Q01468) and the Steroid Delta-isomerase (EC 5.3.3.1, UniProt ID P07445). Obtained results demonstrate that our approach is in accordance to existing experimental data and long-term biophysical and biochemical assumptions.