Improvign mixing processes with the help of numerical tools

In several industrial applications, highly complex behaviour materials are used together with intricate mixing processes, which difficult the achievement of the desired properties for the produced materials. This is the case of the well-known dispersion of nano-sized fillers in a melt polymer matrix, used to improve the nanocomposite mechanical and/or electrical properties. This mixing is usually performed in twin-screw extruders, that promote complex flow patterns, and, since an in loco analysis of the material evolution and mixing is difficult to perform, numerical tools can be very useful to predict the evolution and behaviour of the material. This work presents a numerical based study to improve the understanding of mixing processes. Initial numerical studies were performed with generalized Newtonian fluids, but, due to the null relaxation time that characterize this type of fluids, the assumption of viscoelastic behavior was required. Therefore, the polymer melt was rheologically characterized, and, a six mode Phan-Thien-Tanner and Giesekus models were used to fit the rheological data. These viscoelastic rheological models were used to model the process. The conclusions obtained in this work provide additional and useful data to correlate the type and intensity of the deformation history promoted to the polymer nanocomposite and the quality of the mixing obtained.

Modelação, simulação e implementação de agitação acústica para aplicações em microfluídica

Tese de Doutoramento (Programa Doutoral em Engenharia Biomédica)

A Real-Time intelligent system for tracking patient condition

Hospitals have multiple data sources, such as embedded systems, monitors and sensors. The number of data available is increasing and the information are used not only to care the patient but also to assist the decision processes. The introduction of intelligent environments in health care institutions has been adopted due their ability to provide useful information for health professionals, either in helping to identify prognosis or also to understand patient condition. Behind of this concept arises this Intelligent System to track patient condition (e.g. critic events) in health care. This system has the great advantage of being adaptable to the environment and user needs. The system is focused in identifying critic events from data streaming (e.g. vital signs and ventilation) which is particularly valuable for understanding the patient’s condition. This work aims to demonstrate the process of creating an intelligent system capable of operating in a real environment using streaming data provided by ventilators and vital signs monitors. Its development is important to the physician because becomes possible crossing multiple variables in real-time by analyzing if a value is critic or not and if their variation has or not clinical importance.

Development of nanohydrogels as a vehicle for the release of bioactive compounds in food matrices

PhD in Chemical and Biological Engineering

A multi–physics approach applied to masonry structures with non–hydraulic lime mortars

Tese de Doutoramento em Engenharia Civil (área de especialização em Engenharia de Estruturas).

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.

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.

Modeling and simulation of population balances for particulate processes

This work focuses on the modeling and numerical approximations of population balance equations (PBEs) for the simulation of different phenomena occurring in process engineering. The population balance equation (PBE) is considered to be a statement of continuity. It tracks the change in particle size distribution as particles are born, die, grow or leave a given control volume. In the population balance models the one independent variable represents the time, the other(s) are property coordinate(s), e.g., the particle volume (size) in the present case. They typically describe the temporal evolution of the number density functions and have been used to model various processes such as granulation, crystallization, polymerization, emulsion and cell dynamics. The semi-discrete high resolution schemes are proposed for solving PBEs modeling one and two-dimensional batch crystallization models. The schemes are discrete in property coordinates but continuous in time. The resulting ordinary differential equations can be solved by any standard ODE solver. To improve the numerical accuracy of the schemes a moving mesh technique is introduced in both one and two-dimensional cases ...

Coupling and monotonicity of queueing processes

The main purpose of this work is to give a survey of main monotonicity properties of queueing processes based on the coupling method. The literature on this topic is quite extensive, and we do not consider all aspects of this topic. Our more concrete goal is to select the most interesting basic monotonicity results and give simple and elegant proofs. Also we give a few new (or revised) proofs of a few important monotonicity properties for the queue-size and workload processes both in single-server and multi- server systems. The paper is organized as follows. In Section 1, the basic notions and results on coupling method are given. Section 2 contains known coupling results for renewal processes with focus on construction of synchronized renewal instants for a superposition of independent renewal processes. In Section 3, we present basic monotonicity results for the queue-size and workload processes. We consider both discrete-and continuous-time queueing systems with single and multi servers. Less known results on monotonicity of queueing processes with dependent service times and interarrival times are also presented. Section 4 is devoted to monotonicity of general Jackson-type queueing networks with Markovian routing. This section is based on the notable paper [17]. Finally, Section 5 contains elements of stability analysis of regenerative queues and networks, where coupling and monotonicity results play a crucial role to establish minimal suficient stability conditions. Besides, we present some new monotonicity results for tandem networks.

In situ evaluation of single-cell lysis by cytosol extraction observation through fluorescence decay and dielectrophoretic trapping time

We present a new method for lysis of single cells in continuous flow, where cells are sequentially trapped, lysed and released in an automatic process. Using optimized frequencies, dielectrophoretic trapping allows exposing cells in a reproducible way to high electrical fields for long durations, thereby giving good control on the lysis parameters. In situ evaluation of cytosol extraction on single cells has been studied for Chinese hamster ovary (CHO) cells through out-diffusion of fluorescent molecules for different voltage amplitudes. A diffusion model is proposed to correlate this out-diffusion to the total area of the created pores, which is dependent on the potential drop across the cell membrane and enables evaluation of the total pore area in the membrane. The dielectrophoretic trapping is no longer effective after lysis because of the reduced conductivity inside the cells, leading to cell release. The trapping time is linked to the time required for cytosol extraction and can thus provide additional validation of the effective cytosol extraction for non-fluorescent cells. Furthermore, the application of one single voltage for both trapping and lysis provides a fully automatic process including cell trapping, lysis, and release, allowing operating the device in continuous flow without human intervention.

Adaptive Continuous time Markov Chain Approximation Model to General Jump-Diusions

We propose a non-equidistant Q rate matrix formula and an adaptive numerical algorithm for a continuous time Markov chain to approximate jump-diffusions with affine or non-affine functional specifications. Our approach also accommodates state-dependent jump intensity and jump distribution, a flexibility that is very hard to achieve with other numerical methods. The Kolmogorov-Smirnov test shows that the proposed Markov chain transition density converges to the one given by the likelihood expansion formula as in Ait-Sahalia (2008). We provide numerical examples for European stock option pricing in Black and Scholes (1973), Merton (1976) and Kou (2002).

Quantitative comparison of reconstruction methods for intra-voxel fiber recovery from diffusion MRI.

Validation is arguably the bottleneck in the diffusion magnetic resonance imaging (MRI) community. This paper evaluates and compares 20 algorithms for recovering the local intra-voxel fiber structure from diffusion MRI data and is based on the results of the "HARDI reconstruction challenge" organized in the context of the "ISBI 2012" conference. Evaluated methods encompass a mixture of classical techniques well known in the literature such as diffusion tensor, Q-Ball and diffusion spectrum imaging, algorithms inspired by the recent theory of compressed sensing and also brand new approaches proposed for the first time at this contest. To quantitatively compare the methods under controlled conditions, two datasets with known ground-truth were synthetically generated and two main criteria were used to evaluate the quality of the reconstructions in every voxel: correct assessment of the number of fiber populations and angular accuracy in their orientation. This comparative study investigates the behavior of every algorithm with varying experimental conditions and highlights strengths and weaknesses of each approach. This information can be useful not only for enhancing current algorithms and develop the next generation of reconstruction methods, but also to assist physicians in the choice of the most adequate technique for their studies.

Effects of anxiety on visuospatial pattern discrimination and decision-making processes: what can we learn from brain dynamics ?

Detection and discrimination of visuospatial input involve at least extracting, selecting and encoding relevant information and decision-making processes allowing selecting a response. These two operations are altered, respectively, by attentional mechanisms that change discrimination capacities, and by beliefs concerning the likelihood of uncertain events. Information processing is tuned by the attentional level that acts like a filter on perception, while decision-making processes are weighed by subjective probability of risk. In addition, it has been shown that anxiety could affect the detection of unexpected events through the modification of the level of arousal. Consequently, purpose of this study concerns whether and how decision-making and brain dynamics are affected by anxiety. To investigate these questions, the performance of women with either a high (12) or a low (12) STAI-T (State-Trait Anxiety Inventory, Spielberger, 1983) was examined in a decision-making visuospatial task where subjects have to recognize a target visual pattern from non-target patterns. The target pattern was a schematic image of furniture arranged in such a way as to give the impression of a living room. Non-target patterns were created by either the compression or the dilatation of the distances between objects. Target and non-target patterns were always presented in the same configuration. Preliminary behavioral results show no group difference in reaction time. In addition, visuo-spatial abilities were analyzed trough the signal detection theory for quantifying perceptual decisions in the presence of uncertainty (Green and Swets, 1966). This theory treats detection of a stimulus as a decision-making process determined by the nature of the stimulus and cognitive factors. Astonishingly, no difference in d' (corresponding to the distance between means of the distributions) and c (corresponds to the likelihood ratio) indexes was observed. Comparison of Event-related potentials (ERP) reveals that brain dynamics differ according to anxiety. It shows differences in component latencies, particularly a delay in anxious subjects over posterior electrode sites. However, these differences are compensated during later components by shorter latencies in anxious subjects compared to non-anxious one. These inverted effects seem indicate that the absence of difference in reaction time rely on a compensation of attentional level that tunes cortical activation in anxious subjects, but they have to hammer away to maintain performance.

Subduction and obduction processes in the Swiss Alps

The significance of the Brianconnais domain in the Alpine orogen is reviewed in the light of data concerning its collision with the active Adriatic margin and the passive Helvetic margin. The Brianconnais which formerly belonged to the Iberian plate, was located on the northern margin of the Alpine Tethys (Liguro-Piemont ocean) since its opening in the early-Middle Jurassic. Together with the Iberian plate the Brianconnais terrane was separated from the European plate in the Late Jurassic-Early Cretaceous, following the northern Atlantic, Bay of Biscay, Valais ocean opening. This was accompanied by the onset of subduction along the northern margin of Adria and the closure of the Alpine Tethys. Stratigraphic and metamorphic data regarding this subduction and the geohistory of the Brianconnais allows the scenario of subduction-obduction processes during the Late Cretaceous-early Tertiary in the eastern and western Alps to be specified. HP-LT metamorphism record a long-lasting history of oceanic subduction-accretion, followed in the Middle Eocene by the incorporation of the Brianconnais as an exotic terrane into the accretionary prism. Middle to Late Eocene cooling ages of the Brianconnais basement and the presence of pelagic, anorogenic sedimentation lasting until the Middle Eocene on the Brianconnais preclude any sort of collision before that time between this domain and the active Adria margin or the Helvetic margin. This is confirmed by plate reconstructions constrained by magnetic anomalies in the Atlantic domain. Only a small percentage of the former Brianconnais domain was obducted, most of the crust and lithospheric roots were subducted. This applies also to domains formerly belonging to the southern Alpine Tethys margin (Austroalpine-inner Carpathian domain). It is proposed that there was a single Palaeogene subduction zone responsible for the Alpine orogen formation (from northern Spain to the East Carpathians), with the exception of a short-lived Late Cretaceous partial closure of the Valais ocean. Subduction in the western Tethyan domain originated during the closure of the Meliata ocean during the Jurassic incorporating the Austroalpine-Carpathian domain as terranes during the Cretaceous. The subduction zone propagated into the northern margin of Adria and then to the northern margin of the Iberian plate, where it gave birth to the Pyrenean-Provencal orogenic belt. This implies the absence of a separated Cretaceous subduction zone within the Austro-Carpathian Penninic ocean. Collision of Iberia with Europe forced the subduction to jump to the SE margin of Iberia in the Eocene, creating the Apenninic orogenic wedge and inverting the vergence of subduction from south- to north-directed. (C) 1998 Elsevier Science B.V. All rights reserved.

Prognostic accuracy of cerebral blood flow measurement by perfusion computed tomography, at the time of emergency room admission, in acute stroke patients.

The purpose of this study was to determine the prognostic accuracy of perfusion computed tomography (CT), performed at the time of emergency room admission, in acute stroke patients. Accuracy was determined by comparison of perfusion CT with delayed magnetic resonance (MR) and by monitoring the evolution of each patient's clinical condition. Twenty-two acute stroke patients underwent perfusion CT covering four contiguous 10mm slices on admission, as well as delayed MR, performed after a median interval of 3 days after emergency room admission. Eight were treated with thrombolytic agents. Infarct size on the admission perfusion CT was compared with that on the delayed diffusion-weighted (DWI)-MR, chosen as the gold standard. Delayed magnetic resonance angiography and perfusion-weighted MR were used to detect recanalization. A potential recuperation ratio, defined as PRR = penumbra size/(penumbra size + infarct size) on the admission perfusion CT, was compared with the evolution in each patient's clinical condition, defined by the National Institutes of Health Stroke Scale (NIHSS). In the 8 cases with arterial recanalization, the size of the cerebral infarct on the delayed DWI-MR was larger than or equal to that of the infarct on the admission perfusion CT, but smaller than or equal to that of the ischemic lesion on the admission perfusion CT; and the observed improvement in the NIHSS correlated with the PRR (correlation coefficient = 0.833). In the 14 cases with persistent arterial occlusion, infarct size on the delayed DWI-MR correlated with ischemic lesion size on the admission perfusion CT (r = 0.958). In all 22 patients, the admission NIHSS correlated with the size of the ischemic area on the admission perfusion CT (r = 0.627). Based on these findings, we conclude that perfusion CT allows the accurate prediction of the final infarct size and the evaluation of clinical prognosis for acute stroke patients at the time of emergency evaluation. It may also provide information about the extent of the penumbra. Perfusion CT could therefore be a valuable tool in the early management of acute stroke patients.