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.

Projecto de fundações e estrutura de um edifício destinado a pavilhão gimnodesportivo

Este trabalho tem como objectivo a elaboração do projecto de estruturas de um edifício destinado a pavilhão gimnodesportivo, caracterizando as suas diferentes fases de execução, desde a etapa inicial de concepção até à fase final de dimensionamento. Trata-se de um projecto complexo de uma estrutura com elementos estruturais em betão armado e pré-esforçado, e com muros de contenção. Na concepção do edifício foram utilizados os critérios gerais de dimensionamento presentes na regulamentação Europeia (Eurocódigos), uma vez que estes elementos representam o futuro da regulamentação de estruturas em termos Europeus, vindo substituir a nível nacional o “Regulamento de Segurança e Acções para Estruturas de Betão Armado (RSA)” e o “Regulamento para Estruturas de Betão Armado e Pré- Esforçado (REBAP)”. A adopção das normas europeias representam assim um elevado desafio devido ao aumento da complexidade na concepção e dimensionamento de estruturas que estes regulamentos traduzem, principalmente o Eurocódigo 8, que define de um modo mais detalhado e complexo a análise sísmica, relativamente à regulamentação actual em vigor. Devido à elevada complexidade que os projectos de estruturas apresentam, utilizam-se actualmente ferramentas de cálculo automático. No dimensionamento deste edifício foi utilizado um programa tridimensional de elementos finitos para a modelação da estrutura. Pretende-se com a escolha deste projecto e dos métodos de dimensionamento presentes nos Eurocódigos, o desenvolvimento de um trabalho detalhado e correcto, permitindo assim adquirir conhecimentos importantes relativamente às futuras normas, e pôr em prática as competências e os conhecimentos obtidos ao longo curso.

The shape of soap films and Plateau borders

We have calculated the shapes of flat liquid films, and of the transition region to the associated Plateau borders (PBs), by integrating the Laplace equation with a position-dependent surface tension γ(x), where 2x is the local film thickness. We discuss films in either zero or non-zero gravity, using standard γ(x) potentials for the interaction between the two bounding surfaces. We have investigated the effects of the film flatness, liquid underpressure, and gravity on the shape of films and their PBs. Films may exhibit 'humps' and/or 'dips' associated with inflection points and minima of the film thickness. Finally, we propose an asymptotic analytical solution for the film width profile.

Contact angle of a hemispherical bubble: An analytical approach

We have calculated the equilibrium shape of the axially symmetric Plateau border along which a spherical bubble contacts a flat wall, by analytically integrating Laplace's equation in the presence of gravity, in the limit of small Plateau border sizes. This method has the advantage that it provides closed-form expressions for the positions and orientations of the Plateau border surfaces. Results are in very good overall agreement with those obtained from a numerical solution procedure, and are consistent with experimental data. In particular we find that the effect of gravity on Plateau border shape is relatively small for typical bubble sizes, leading to a widening of the Plateau border for sessile bubbles and to a narrowing for pendant bubbles. The contact angle of the bubble is found to depend even more weakly on gravity. (C) 2009 Elsevier Inc. All rights reserved.

Positive Solutions of the Dirichlet Problem for the One-dimensional Minkowski-Curvature Equation

We discuss existence and multiplicity of positive solutions of the Dirichlet problem for the quasilinear ordinary differential equation-(u' / root 1 - u'(2))' = f(t, u). Depending on the behaviour of f = f(t, s) near s = 0, we prove the existence of either one, or two, or three, or infinitely many positive solutions. In general, the positivity of f is not required. All results are obtained by reduction to an equivalent non-singular problem to which variational or topological methods apply in a classical fashion.

A one-dimensional prescribed curvature equation modeling the corneal shape

We prove existence, uniqueness, and stability of solutions of the prescribed curvature problem (u'/root 1 + u'(2))' = au - b/root 1 + u'(2) in [0, 1], u'(0) = u(1) = 0, for any given a > 0 and b > 0. We also develop a linear monotone iterative scheme for approximating the solution. This equation has been proposed as a model of the corneal shape in the recent paper (Okrasinski and Plociniczak in Nonlinear Anal., Real World Appl. 13:1498-1505, 2012), where a simplified version obtained by partial linearization has been investigated.

A one-dimensional prescribed curvature equation modeling the corneal shape.

We prove existence, uniqueness, and stability of solutions of the prescribed curvature problem (u'/root 1 + u'(2))' = au - b/root 1 + u'(2) in [0, 1], u'(0) = u(1) = 0, for any given a > 0 and b > 0. We also develop a linear monotone iterative scheme for approximating the solution. This equation has been proposed as a model of the corneal shape in the recent paper (Okrasinski and Plociniczak in Nonlinear Anal., Real World Appl. 13:1498-1505, 2012), where a simplified version obtained by partial linearization has been investigated.

Scaling law in Saddle-node bifurcations for one-dimensional maps: A complex variable approach

The study of transient dynamical phenomena near bifurcation thresholds has attracted the interest of many researchers due to the relevance of bifurcations in different physical or biological systems. In the context of saddle-node bifurcations, where two or more fixed points collide annihilating each other, it is known that the dynamics can suffer the so-called delayed transition. This phenomenon emerges when the system spends a lot of time before reaching the remaining stable equilibrium, found after the bifurcation, because of the presence of a saddle-remnant in phase space. Some works have analytically tackled this phenomenon, especially in time-continuous dynamical systems, showing that the time delay, tau, scales according to an inverse square-root power law, tau similar to (mu-mu (c) )(-1/2), as the bifurcation parameter mu, is driven further away from its critical value, mu (c) . In this work, we first characterize analytically this scaling law using complex variable techniques for a family of one-dimensional maps, called the normal form for the saddle-node bifurcation. We then apply our general analytic results to a single-species ecological model with harvesting given by a unimodal map, characterizing the delayed transition and the scaling law arising due to the constant of harvesting. For both analyzed systems, we show that the numerical results are in perfect agreement with the analytical solutions we are providing. The procedure presented in this work can be used to characterize the scaling laws of one-dimensional discrete dynamical systems with saddle-node bifurcations.

Explicit series solution for a glucose-induced electrical activity model of pancreatic beta-cells

In this work, we present the explicit series solution of a specific mathematical model from the literature, the Deng bursting model, that mimics the glucose-induced electrical activity of pancreatic beta-cells (Deng, 1993). To serve to this purpose, we use a technique developed to find analytic approximate solutions for strongly nonlinear problems. This analytical algorithm involves an auxiliary parameter which provides us with an efficient way to ensure the rapid and accurate convergence to the exact solution of the bursting model. By using the homotopy solution, we investigate the dynamical effect of a biologically meaningful bifurcation parameter rho, which increases with the glucose concentration. Our analytical results are found to be in excellent agreement with the numerical ones. This work provides an illustration of how our understanding of biophysically motivated models can be directly enhanced by the application of a newly analytic method.

Guiding Monte Carlo tree searches with neural networks in the game of go

A dissertation submitted in fulfillment of the requirements to the degree of Master in Computer Science and Computer Engineering