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.

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

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

Self-monitoring of freeze–thaw damage using triphasic electric conductive concrete

The effect of freeze–thaw cycles on concrete is of great importance for durability evaluation of concrete structures in cold regions. In this paper, damage accumulation was studied by following the fractional change of impedance (FCI) with number of freeze–thaw cycles (N). The nano-carbon black (NCB), carbon fiber (CF) and steel fiber (SF) were added to plain concrete to produce the triphasic electrical conductive (TEC) and ductile concrete. The effects of NCB, CF and SF on the compressive strength, flexural properties, electrical impedance were investigated. The concrete beams with different dosages of conductive materials were studied for FCI, N and mass loss (ML), the relationship between FCI and N of conductive concrete can be well defined by a first order exponential decay curve. It is noted that this nondestructive and sensitive real-time testing method is meaningful for evaluating of freeze–thaw damage in concrete.

Annular flow of viscoelastic fluids: Analytical and numerical solutions

This work provides analytical and numerical solutions for the linear, quadratic and exponential Phan–Thien–Tanner (PTT) viscoelastic models, for axial and helical annular fully-developed flows under no slip and slip boundary conditions, the latter given by the linear and nonlinear Navier slip laws. The rheology of the three PTT model functions is discussed together with the influence of the slip velocity upon the flow velocity and stress fields. For the linear PTT model, full analytical solutions for the inverse problem (unknown velocity) are devised for the linear Navier slip law and two different slip exponents. For the linear PTT model with other values of the slip exponent and for the quadratic PTT model, the polynomial equation for the radial location (β) of the null shear stress must be solved numerically. For both models, the solution of the direct problem is given by an iterative procedure involving three nonlinear equations, one for β, other for the pressure gradient and another for the torque per unit length. For the exponential PTT model we devise a numerical procedure that can easily compute the numerical solution of the pure axial flow problem

A sixth-order finite volume method for the 1D biharmonic operator: application to intramedullary nail simulation

A new very high-order finite volume method to solve problems with harmonic and biharmonic operators for one- dimensional geometries is proposed. The main ingredient is polynomial reconstruction based on local interpolations of mean values providing accurate approximations of the solution up to the sixth-order accuracy. First developed with the harmonic operator, an extension for the biharmonic operator is obtained, which allows designing a very high-order finite volume scheme where the solution is obtained by solving a matrix-free problem. An application in elasticity coupling the two operators is presented. We consider a beam subject to a combination of tensile and bending loads, where the main goal is the stress critical point determination for an intramedullary nail.

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.

Expressão e caracterização de uma lectina recombinante de interesse biomédico

Dissertação de mestrado integrado em Engenharia Biomédica (área de especialização em Engenharia Clínica)

Concrete with triphasic conductive materials for self-monitoring of cracking development subjected to flexure

In this study, the macro steel fiber (SF), carbon fiber (CF) and nano carbon black (NCB) as triphasic conductive materials were added into concrete, in order to improve the conductivity and ductility of concrete. The influence of NCB, SF and CF on the post crack behavior and conductivity of concrete was explored. The effect of the triphasic conductive materials on the self-diagnosing ability to the load–deflection property and crack widening of conductive concrete member subjected to bending were investigated. The relationship between the fractional change in surface impedance (FCR) and the crack opening displacement (COD) of concrete beams with conductive materials has been established. The results illustrated that there is a linear relationship between COD and FCR.

Finite integration methods for isospectral flows

In this paper we consider the approximate computation of isospectral flows based on finite integration methods( FIM) with radial basis functions( RBF) interpolation,a new algorithm is developed. Our method ensures the symmetry of the solutions. Numerical experiments demonstrate that the solutions have higher accuracy by our algorithm than by the second order Runge- Kutta( RK2) method.

Análise técnica - económica da produção de energia em sistemas de abastecimento de água: caso do aproveitamento de fins múltiplos dos socorridos

Dissertação de mestrado integrado em Engenharia Civil