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

A heat transfer model of the human upper limbs

A steady state multi-segmented heat transfer model of the human upper limbs was developed. The main purpose was to evaluate the impact of blood flow through superficial veins and subcutaneous vascular structures in the palm of the hands over the heat transfer between the limbs and the environment. The distinguishing feature of the model is the inclusion of a detailed circulatory network composed of vessels with diameter larger than 1 mm. The model was validated by comparing its results from exposures to a hot, a neutral, and a cold environment to experimental data presented in the literature. (C) 2011 Elsevier Ltd. All rights reserved.

Modelagem dos efeitos da irradiação luminosa no cérebro de camundongos e rastreamento de neurônios durante experimentos de microscopia de fluorescência

The fluorescent proteins are an essential tool in many fields of biology, since they allow us to watch the development of structures and dynamic processes of cells in living tissue, with the aid of fluorescence microscopy. Optogenectics is another technique that is currently widely used in Neuroscience. In general, this technique allows to activate/deactivate neurons with the radiation of certain wavelengths on the cells that have ion channels sensitive to light, at the same time that can be used with fluorescent proteins. This dissertation has two main objectives. Initially, we study the interaction of light radiation and mice brain tissue to be applied in optogenetic experiments. In this step, we model absorption and scattering effects using mice brain tissue characteristics and Kubelka-Munk theory, for specific wavelengths, as a function of light penetration depth (distance) within the tissue. Furthermore, we model temperature variations using the finite element method to solve Pennes’ bioheat equation, with the aid of COMSOL Multiphysics Modeling Software 4.4, where we simulate protocols of light stimulation tipically used in optogenetics. Subsequently, we develop some computational algorithms to reduce the exposure of neuron cells to the light radiation necessary for the visualization of their emitted fluorescence. At this stage, we describe the image processing techniques developed to be used in fluorescence microscopy to reduce the exposure of the brain samples to continuous light, which is responsible for fluorochrome excitation. The developed techniques are able to track, in real time, a region of interest (ROI) and replace the fluorescence emitted by the cells by a virtual mask, as a result of the overlay of the tracked ROI and the fluorescence information previously stored, preserving cell location, independently of the time exposure to fluorescent light. In summary, this dissertation intends to investigate and describe the effects of light radiation in brain tissue, within the context of Optogenetics, in addition to providing a computational tool to be used in fluorescence microscopy experiments to reduce image bleaching and photodamage due to the intense exposure of fluorescent cells to light radiation.

SAR ed effetti termici su Cheratinociti sottoposti a Radiofrequenze

Negli ultimi anni, la popolazione è stata esposta a vari tipi di campi elettromagnetici generati da strumentazioni elettroniche e dispositivi di telecomunicazione. In questa tesi, si valutano SAR ed effetti termici prodotti da tre antenne patch a radiofrequenza sia su Cheratinociti (cellule dell'epidermide) in vitro che sull'epidermide umana in vivo caratterizzata in un modello multistrato contenente tessuti biologici. Le antenne progettate hanno frequenze di risonanza di 1.8 e 2.4 GHz, tipiche delle bande utilizzate rispettivamente da LTE (Long Term Evolution, la più recente evoluzione degli standard di telefonia mobile cellulare) e dalle moderne tecnologie Wi-Fi, e di 60 GHz, propria delle cosiddette onde millimetriche. Vengono valutati quindi il SAR (Specific Absorption Rate, grandezza che fornisce una misura dell'assorbimento delle onde da parte dei tessuti biologici) e le variazioni di temperatura prodotte dall'applicazione del campo elettromagnetico: ciò viene realizzato attraverso l'equazione del calore stazionaria e, nel caso dell'epidermide in vivo, con la Bioheat Equation, che contempla anche la circolazione sanguigna ed il calore generato nei processi metabolici che avvengono nell'organismo.

Study of the one dimensional and transient bioheat transfer equation: Multi-layer solution development and applications

: In this work we derive an analytical solution given by Bessel series to the transient and one-dimensional (1D) bioheat transfer equation in a multi-layer region with spatially dependent heat sources. Each region represents an independent biological tissue characterized by temperature-invariant physiological parameters and a linearly temperature dependent metabolic heat generation. Moreover, 1D Cartesian, cylindrical or spherical coordinates are used to define the geometry and temperature boundary conditions of first, second and third kinds are assumed at the inner and outer surfaces. We present two examples of clinical applications for the developed solution. In the first one, we investigate two different heat source terms to simulate the heating in a tumor and its surrounding tissue, induced during a magnetic fluid hyperthermia technique used for cancer treatment. To obtain an accurate analytical solution, we determine the error associated with the truncated Bessel series that defines the transient solution. In the second application, we explore the potential of this model to study the effect of different environmental conditions in a multi-layered human head model (brain, bone and scalp). The convective heat transfer effect of a large blood vessel located inside the brain is also investigated. The results are further compared with a numerical solution obtained by the Finite Element Method and computed with COMSOL Multi-physics v4.1 (c). (c) 2013 Elsevier Ltd. All rights reserved.

A method of fundamental solutions for theaxisymmetric backward heat equation

We propose and investigate an application of the method of fundamental solutions (MFS) to the radially symmetric and axisymmetric backward heat conduction problem (BHCP) in a solid or hollow cylinder. In the BHCP, the initial temperature is to be determined from the temperature measurements at a later time. This is an inverse and ill-posed problem, and we employ and generalize the MFS regularization approach [B.T. Johansson and D. Lesnic, A method of fundamental solutions for transient heat conduction, Eng. Anal. Boundary Elements 32 (2008), pp. 697–703] for the time-dependent heat equation to obtain a stable and accurate numerical approximation with small computational cost.

Properties of a method of fundamental solutions for the parabolic heat equation

We show that a set of fundamental solutions to the parabolic heat equation, with each element in the set corresponding to a point source located on a given surface with the number of source points being dense on this surface, constitute a linearly independent and dense set with respect to the standard inner product of square integrable functions, both on lateral- and time-boundaries. This result leads naturally to a method of numerically approximating solutions to the parabolic heat equation denoted a method of fundamental solutions (MFS). A discussion around convergence of such an approximation is included.

Avaliação do crescimento físico e maturação somática de escolares da rede de ensino fundamental de Campinas, SP

Universidade Estadual de Campinas . Faculdade de Educação Física

Fundamental studies on natural deep eutectic solvents: physico-chemical, thermal and rheological properties

When combined at particular molar fractions, sugars, aminoacids or organic acids a present a high melting point depression, becoming liquids at room temperature. These are called Natural Deep Eutectic Solvents – NADES and are envisaged to play a major role on the chemical engineering processes of the future. Nonetheless, there is a significant lack of knowledge of its fundamental and basic properties, which is hindering their industrial applications. For this reason it is important to extend the knowledge on these systems, boosting their application development [1]. In this work, we have developed and characterized NADES based on choline chloride, organic acids, amino acids and sugars. Their density, thermal behavior, conductivity and polarity were assessed for different compositions. The conductivity was measured from 0 to 40 °C and the temperature effect was well described by the Vogel-Fulcher-Tammann equation. The morphological characterization of the crystallizable materials was done by polarized optical microscopy that provided also evidence of homogeneity/phase separation. Additionally, the rheological and thermodynamic properties of the NADES and the effect of water content were also studied. The results show these systems have Newtonian behavior and present significant viscosity decrease with temperature and water content, due to increase on the molecular mobility. The anhydrous systems present viscosities that range from higher than 1000Pa.s at 20°C to less than 1Pa.s at 70°C. DSC characterization confirms that for water content as high as 1:1:1 molar ratio, the mixture retains its single phase behavior. The results obtained demonstrate that the NADES properties can be finely tunned by careful selection of its constituents. NADES present the necessary properties for use as extraction solvents. They can be prepared from inexpensive raw materials and tailored for the selective extraction of target molecules. The data produced in this work is hereafter importance for the selection of the most promising candidates avoiding a time consuming and expensive trial and error phase providing also data for the development of models able to predict their properties and the mechanisms that allow the formation of the deep eutectic mixtures.

A well-posed integral equation formulation for three-dimensional rough surface scattering

We consider the problem of scattering of time-harmonic acoustic waves by an unbounded sound-soft rough surface. Recently, a Brakhage Werner type integral equation formulation of this problem has been proposed, based on an ansatz as a combined single- and double-layer potential, but replacing the usual fundamental solution of the Helmholtz equation with an appropriate half-space Green's function. Moreover, it has been shown in the three-dimensional case that this integral equation is uniquely solvable in the space L-2 (Gamma) when the scattering surface G does not differ too much from a plane. In this paper, we show that this integral equation is uniquely solvable with no restriction on the surface elevation or slope. Moreover, we construct explicit bounds on the inverse of the associated boundary integral operator, as a function of the wave number, the parameter coupling the single- and double-layer potentials, and the maximum surface slope. These bounds show that the norm of the inverse operator is bounded uniformly in the wave number, kappa, for kappa > 0, if the coupling parameter h is chosen proportional to the wave number. In the case when G is a plane, we show that the choice eta = kappa/2 is nearly optimal in terms of minimizing the condition number.

The omega equation and potential vorticity

Two fundamental perspectives on the dynamics of midlatitude weather systems are provided by potential vorticity (PV) and the omega equation. The aim of this paper is to investigate the link between the two perspectives, which has so far received very little attention in the meteorological literature. It also aims to give a quantitative basis for discussion of quasi-geostrophic vertical motion in terms of components associated with system movement, maintaining a constant thermal structure, and with the development of that structure. The former links with the isentropic relative-flow analysis technique. Viewed in a moving frame of reference, the measured development of a system depends on the velocity of that frame of reference. The requirement that the development should be a minimum provides a quantitative method for determining the optimum system velocity. The component of vertical velocity associated with development is shown to satisfy an omega equation with forcing determined from the relative advection of interior PV and boundary temperature. The analysis carries through in the presence of diabatic heating provided the omega equation forcing is based on the interior PV and boundary thermal tendencies, including the heating effect. The analysis is shown to be possible also at the level of the semi-geostrophic approximation. The analysis technique is applied to a number of idealized problems that can be considered to be building blocks for midlatitude synoptic-scale dynamics. They focus on the influences of interior PV, boundary temperature, an interior boundary, baroclinic instability associated with two boundaries, and also diabatic heating. In each case, insights yielded by the new perspective are sought into the dynamical behaviour, especially that related to vertical motion. Copyright © 2003 Royal Meteorological Society

Scattering by rough surfaces: the Dirichlet problem for the Helmholtz equation in a non-locally perturbed half-plane

We consider the Dirichlet boundary value problem for the Helmholtz equation in a non-locally perturbed half-plane, this problem arising in electromagnetic scattering by one-dimensional rough, perfectly conducting surfaces. We propose a new boundary integral equation formulation for this problem, utilizing the Green's function for an impedance half-plane in place of the standard fundamental solution. We show, at least for surfaces not differing too much from the flat boundary, that the integral equation is uniquely solvable in the space of bounded and continuous functions, and hence that, for a variety of incident fields including an incident plane wave, the boundary value problem for the scattered field has a unique solution satisfying the limiting absorption principle. Finally, a result of continuous dependence of the solution on the boundary shape is obtained.

An historical perspective on the development of the thermodynamic equation of seawater - 2010

Oceanography is concerned with understanding the mechanisms controlling the movement of seawater and its contents. A fundamental tool in this process is the characterization of the thermophysical properties of seawater as functions of measured temperature and electrical conductivity, the latter used as a proxy for the concentration of dissolved matter in seawater. For many years a collection of algorithms denoted the Equation of State 1980 (EOS-80) has been the internationally accepted standard for calculating such properties. However, modern measurement technology now allows routine observations of temperature and electrical conductivity to be made to at least one order of magnitude more accurately than the uncertainty in this standard. Recently, a new standard has been developed, the Thermodynamical Equation of Seawater 2010 (TEOS-10). This new standard is thermodynamically consistent, valid over a wider range of temperature and salinity, and includes a mechanism to account for composition variations in seawater. Here we review the scientific development of this standard, and describe the literature involved in its development, which includes many of the articles in this special issue.

MODELLING OF HIGH-PRESSURE PHASE EQUILIBRIUM IN SYSTEMS OF INTEREST IN THE FOOD ENGINEERING FIELD USING THE PENG-ROBINSON EQUATION OF STATE WITH TWO DIFFERENT MIXING RULES

In this work, thermodynamic models for fitting the phase equilibrium of binary systems were applied, aiming to predict the high pressure phase equilibrium of multicomponent systems of interest in the food engineering field, comparing the results generated by the models with new experimental data and with those from the literature. Two mixing rules were used with the Peng-Robinson equation of state, one with the mixing rule of van der Waals and the other with the composition-dependent mixing rule of Mathias et al. The systems chosen are of fundamental importance in food industries, such as the binary systems CO(2)-limonene, CO(2)-citral and CO(2)-linalool, and the ternary systems CO(2)-Limonene-Citral and CO(2)-Limonene-Linalool, where high pressure phase equilibrium knowledge is important to extract and fractionate citrus fruit essential oils. For the CO(2)-limonene system, some experimental data were also measured in this work. The results showed the high capability of the model using the composition-dependent mixing rule to model the phase equilibrium behavior of these systems.