Computing and Visualizing Solution Sets of Interval Linear Systems

The paper has been presented at the 12th International Conference on Applications of Computer Algebra, Varna, Bulgaria, June, 2006

Cardiac-induced localized thoracic motion detected by a fiber optic sensing scheme

The cardiovascular health of the human population is a major concern for medical clinicians, with cardiovascular diseases responsible for 48% of all deaths worldwide, according to the World Health Organization. The development of new diagnostic tools that are practicable and economical to scrutinize the cardiovascular health of humans is a major driver for clinicians. We offer a new technique to obtain seismocardiographic signals up to 54 Hz covering both ballistocardiography (below 20 Hz) and audible heart sounds (20 Hz upward), using a system based on curvature sensors formed from fiber optic long period gratings. This system can visualize the real-time three-dimensional (3-D) mechanical motion of the heart by using the data from the sensing array in conjunction with a bespoke 3-D shape reconstruction algorithm. Visualization is demonstrated by adhering three to four sensors on the outside of the thorax and in close proximity to the apex of the heart; the sensing scheme revealed a complex motion of the heart wall next to the apex region of the heart. The detection scheme is low-cost, portable, easily operated and has the potential for ambulatory applications.

An integrated functional mapping and source localization platform for brain research

The premise of this dissertation is to create a highly integrated platform that combines the most current recording technologies for brain research through the development of new algorithms for three-dimensional (3D) functional mapping and 3D source localization. The recording modalities that were integrated include: Electroencephalography (EEG), Optical Topographic Maps (OTM), Magnetic Resonance Imaging (MRI), and Diffusion Tensor Imaging (DTI). This work can be divided into two parts: The first part involves the integration of OTM with MRI, where the topographic maps are mapped to both the skull and cortical surface of the brain. This integration process is made possible through the development of new algorithms that determine the probes location on the MRI head model and warping the 2D topographic maps onto the 3D MRI head/brain model. Dynamic changes of the brain activation can be visualized on the MRI head model through a graphical user interface. The second part of this research involves augmenting a fiber tracking system, by adding the ability to integrate the source localization results generated by commercial software named Curry. This task involved registering the EEG electrodes and the dipole results to the MRI data. Such Integration will allow the visualization of fiber tracts, along with the source of the EEG, in a 3D transparent brain structure. The research findings of this dissertation were tested and validated through the participation of patients from Miami Children Hospital (MCH). Such an integrated platform presented to the medical professionals in the form of a user-friendly graphical interface is viewed as a major contribution of this dissertation. It should be emphasized that there are two main aspects to this research endeavor: (1) if a dipole could be situated in time at its different positions, its trajectory may reveal additional information on the extent and nature of the brain malfunction; (2) situating such a dipole trajectory with respect to the fiber tracks could ensure the preservation of these fiber tracks (axons) during surgical interventions, preserving as a consequence these parts of the brain that are responsible for information transmission.

Estratégias para obtenção de adequada iluminação natural em escolas: uma análise de sistemas de aberturas para Natal/RN

The purpose of this research is to analyze different daylighting systems in schools in the city of Natal/RN. Although with the abundantly daylight available locally, there are a scarce and diffuse architectural recommendations relating sky conditions, dimensions of daylight systems, shading, fraction of sky visibility, required illuminance, glare, period of occupation and depth of the lit area. This research explores different selected apertures systems to explore the potential of natural light for each system. The method has divided into three phases: The first phase is the modeling which involves the construction of three-dimensional model of a classroom in Sketchup software 2014, which is featured in follow recommendations presented in the literature to obtain a good quality of environmental comfort in school settings. The second phase is the dynamic performance computer simulation of the light through the Daysim software. The input data are the climate file of 2009 the city of Natal / RN, the classroom volumetry in 3ds format with the assignment of optical properties of each surface, the sensor mapping file and the user load file . The results produced in the simulation are organized in a spreadsheet prepared by Carvalho (2014) to determine the occurrence of useful daylight illuminance (UDI) in the range of 300 to 3000lux and build graphics illuminance curves and contours of UDI to identify the uniformity of distribution light, the need of the minimum level of illuminance and the occurrence of glare.

Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.

Axis Mapping: The Estimation of Surface Orientations and its Applications in Vehicle Localization and Structural Geology

The map representation of an environment should be selected based on its intended application. For example, a geometrically accurate map describing the Euclidean space of an environment is not necessarily the best choice if only a small subset its features are required. One possible subset is the orientations of the flat surfaces in the environment, represented by a special parameterization of normal vectors called axes. Devoid of positional information, the entries of an axis map form a non-injective relationship with the flat surfaces in the environment, which results in physically distinct flat surfaces being represented by a single axis. This drastically reduces the complexity of the map, but retains important information about the environment that can be used in meaningful applications in both two and three dimensions. This thesis presents axis mapping, which is an algorithm that accurately and automatically estimates an axis map of an environment based on sensor measurements collected by a mobile platform. Furthermore, two major applications of axis maps are developed and implemented. First, the LiDAR compass is a heading estimation algorithm that compares measurements of axes with an axis map of the environment. Pairing the LiDAR compass with simple translation measurements forms the basis for an accurate two-dimensional localization algorithm. It is shown that this algorithm eliminates the growth of heading error in both indoor and outdoor environments, resulting in accurate localization over long distances. Second, in the context of geotechnical engineering, a three-dimensional axis map is called a stereonet, which is used as a tool to examine the strength and stability of a rock face. Axis mapping provides a novel approach to create accurate stereonets safely, rapidly, and inexpensively compared to established methods. The non-injective property of axis maps is leveraged to probabilistically describe the relationships between non-sequential measurements of the rock face. The automatic estimation of stereonets was tested in three separate outdoor environments. It is shown that axis mapping can accurately estimate stereonets while improving safety, requiring significantly less time and effort, and lowering costs compared to traditional and current state-of-the-art approaches.

A framework for experimental determination of localised vertical pedestrian forces on full-scale structures using wireless attitude and heading reference systems

A major weakness among loading models for pedestrians walking on flexible structures proposed in recent years is the various uncorroborated assumptions made in their development. This applies to spatio-temporal characteristics of pedestrian loading and the nature of multi-object interactions. To alleviate this problem, a framework for the determination of localised pedestrian forces on full-scale structures is presented using a wireless attitude and heading reference systems (AHRS). An AHRS comprises a triad of tri-axial accelerometers, gyroscopes and magnetometers managed by a dedicated data processing unit, allowing motion in three-dimensional space to be reconstructed. A pedestrian loading model based on a single point inertial measurement from an AHRS is derived and shown to perform well against benchmark data collected on an instrumented treadmill. Unlike other models, the current model does not take any predefined form nor does it require any extrapolations as to the timing and amplitude of pedestrian loading. In order to assess correctly the influence of the moving pedestrian on behaviour of a structure, an algorithm for tracking the point of application of pedestrian force is developed based on data from a single AHRS attached to a foot. A set of controlled walking tests with a single pedestrian is conducted on a real footbridge for validation purposes. A remarkably good match between the measured and simulated bridge response is found, indeed confirming applicability of the proposed framework.

Stochastic Modeling of Reversible Biochemical Reaction-Diffusion Systems and High-Resolution Shock-Capturing Methods for Fluid Interfaces

Thesis (Ph.D.)--University of Washington, 2016-08

LARGE MOTION VISUALIZATION AND ESTIMATION FOR FLAPPING WING SYSTEMS

Studies of fluid-structure interactions associated with flexible structures such as flapping wings require the capture and quantification of large motions of bodies that may be opaque. Motion capture of a free flying insect is considered by using three synchronized high-speed cameras. A solid finite element representation is used as a reference body and successive snapshots in time of the displacement fields are reconstructed via an optimization procedure. An objective function is formulated, and various shape difference definitions are considered. The proposed methodology is first studied for a synthetic case of a flexible cantilever structure undergoing large deformations, and then applied to a Manduca Sexta (hawkmoth) in free flight. The three-dimensional motions of this flapping system are reconstructed from image date collected by using three cameras. The complete deformation geometry of this system is analyzed. Finally, a computational investigation is carried out to understand the flow physics and aerodynamic performance by prescribing the body and wing motions in a fluid-body code. This thesis work contains one of the first set of such motion visualization and deformation analyses carried out for a hawkmoth in free flight. The tools and procedures used in this work are widely applicable to the studies of other flying animals with flexible wings as well as synthetic systems with flexible body elements.

Shaping surface acoustic waves for cardiac tissue engineering

The heart is a non-regenerating organ that gradually suffers a loss of cardiac cells and functionality. Given the scarcity of organ donors and complications in existing medical implantation solutions, it is desired to engineer a three-dimensional architecture to successfully control the cardiac cells in vitro and yield true myocardial structures similar to native heart. This thesis investigates the synthesis of a biocompatible gelatin methacrylate hydrogel to promote growth of cardiac cells using biotechnology methodology: surface acoustic waves, to create cell sheets. Firstly, the synthesis of a photo-crosslinkable gelatin methacrylate (GelMA) hydrogel was investigated with different degree of methacrylation concentration. The porous matrix of the hydrogel should be biocompatible, allow cell-cell interaction and promote cell adhesion for growth through the porous network of matrix. The rheological properties, such as polymer concentration, ultraviolet exposure time, viscosity, elasticity and swelling characteristics of the hydrogel were investigated. In tissue engineering hydrogels have been used for embedding cells to mimic native microenvironments while controlling the mechanical properties. Gelatin methacrylate hydrogels have the advantage of allowing such control of mechanical properties in addition to easy compatibility with Lab-on-a-chip methodologies. Secondly in this thesis, standing surface acoustic waves were used to control the degree of movement of cells in the hydrogel and produce three-dimensional engineered scaffolds to investigate in-vitro studies of cardiac muscle electrophysiology and cardiac tissue engineering therapies for myocardial infarction. The acoustic waves were characterized on a piezoelectric substrate, lithium niobate that was micro-fabricated with slanted-finger interdigitated transducers for to generate waves at multiple wavelengths. This characterization successfully created three-dimensional micro-patterning of cells in the constructs through means of one- and two-dimensional non-invasive forces. The micro-patterning was controlled by tuning different input frequencies that allowed manipulation of the cells spatially without any pre- treatment of cells, hydrogel or substrate. This resulted in a synchronous heartbeat being produced in the hydrogel construct. To complement these mechanical forces, work in dielectrophoresis was conducted centred on a method to pattern micro-particles. Although manipulation of particles were shown, difficulties were encountered concerning the close proximity of particles and hydrogel to the microfabricated electrode arrays, dependence on conductivity of hydrogel and difficult manoeuvrability of scaffold from the surface of electrodes precluded measurements on cardiac cells. In addition, COMSOL Multiphysics software was used to investigate the mechanical and electrical forces theoretically acting on the cells. Thirdly, in this thesis the cardiac electrophysiology was investigated using immunostaining techniques to visualize the growth of sarcomeres and gap junctions that promote cell-cell interaction and excitation-contraction of heart muscles. The physiological response of beating of co-cultured cardiomyocytes and cardiac fibroblasts was observed in a synchronous and simultaneous manner closely mimicking the native cardiac impulses. Further investigations were carried out by mechanically stimulating the cells in the three-dimensional hydrogel using standing surface acoustic waves and comparing with traditional two-dimensional flat surface coated with fibronectin. The electrophysiological responses of the cells under the effect of the mechanical stimulations yielded a higher magnitude of contractility, action potential and calcium transient.

Stenosis triggers spread of helical Pseudomonas biofilms in cylindrical flow systems

Biofilms are multicellular bacterial structures that adhere to surfaces and often endow the bacterial population with tolerance to antibiotics and other environmental insults. Biofilms frequently colonize the tubing of medical devices through mechanisms that are poorly understood. Here we studied the helicoidal spread of Pseudomonas putida biofilms through cylindrical conduits of varied diameters in slow laminar flow regimes. Numerical simulations of such flows reveal vortical motion at stenoses and junctions, which enhances bacterial adhesion and fosters formation of filamentous structures. Formation of long, downstream-flowing bacterial threads that stem from narrowings and connections was detected experimentally, as predicted by our model. Accumulation of bacterial biomass makes the resulting filaments undergo a helical instability. These incipient helices then coarsened until constrained by the tubing walls, and spread along the whole tube length without obstructing the flow. A three-dimensional discrete filament model supports this coarsening mechanism and yields simulations of helix dynamics in accordance with our experimental observations. These findings describe an unanticipated mechanism for bacterial spreading in tubing networks which might be involved in some hospital-acquired infections and bacterial contamination of catheters.

Collaborative Systems in Crisis Management: A Proposal for a Conceptual Framework

Part 11: Reference and Conceptual Models

Energy performance comparison of three innovative HVAC systems for renovation through dynamic simulation

In this paper, dynamic simulation was used to compare the energy performance of three innovativeHVAC systems: (A) mechanical ventilation with heat recovery (MVHR) and micro heat pump, (B) exhaustventilation with exhaust air-to-water heat pump and ventilation radiators, and (C) exhaust ventilationwith air-to-water heat pump and ventilation radiators, to a reference system: (D) exhaust ventilation withair-to-water heat pump and panel radiators. System A was modelled in MATLAB Simulink and systems Band C in TRNSYS 17. The reference system was modelled in both tools, for comparison between the two.All systems were tested with a model of a renovated single family house for varying U-values, climates,infiltration and ventilation rates.It was found that A was the best system for lower heating demand, while for higher heating demandsystem B would be preferable. System C was better than the reference system, but not as good as A or B.The difference in energy consumption of the reference system was less than 2 kWh/(m2a) betweenSimulink and TRNSYS. This could be explained by the different ways of handling solar gains, but also bythe fact that the TRNSYS systems supplied slightly more than the ideal heating demand.

A finite element framework for the interplay between delamination and buckling of rubber-like bi-material systems and stretchable electronics

In this study, a finite element (FE) framework for the analysis of the interplay between buckling and delamination of thin layers bonded to soft substrates is proposed. The current framework incorporates the following modeling features: (i) geometrically nonlinear solid shell elements, (ii) geometrically nonlinear cohesive interface elements, and (iii) hyperelastic material constitutive response for the bodies that compose the system. A fully implicit Newton–Raphson solution strategy is adopted to deal with the complex simultaneous presence of geometrical and material nonlinearities through the derivation of the consistent FE formulation. Applications to a rubber-like bi-material system under finite bending and to patterned stiff islands resting on soft substrate for stretchable solar cells subjected to tensile loading are proposed. The results obtained are in good agreement with benchmark results available in the literature, confirming the accuracy and the capabilities of the proposed numerical method for the analysis of complex three-dimensional fracture mechanics problems under finite deformations.