Real-Time High-Quality Image to Mesh Conversion for Finite Element Simulations

Technological Advances in Medical Imaging have enabled the acquisition of images accurately describing biological tissues. Finite Element (FE) methods on these images provide the means to simulate biological phenomena such as brain shift registration, respiratory organ motion, blood flow pressure in vessels, etc. FE methods require the domain of tissues be discretized by simpler geometric elements, such as triangles in two dimensions, tetrahedra in three, and pentatopes in four. This exact discretization is called a mesh . The accuracy and speed of FE methods depend on the quality and fidelity of the mesh used to describe the biological object. Elements with bad quality introduce numerical errors and slower solver convergence. Also, analysis based on poor fidelity meshes do not yield accurate results specially near the surface. In this dissertation, we present the theory and the implementation of both a sequential and a parallel Delaunay meshing technique for 3D and ---for the first time--- 4D space-time domains. Our method provably guarantees that the mesh is a faithful representation of the multi-tissue domain in topological and geometric sense. Moreover, we show that our method generates graded elements of bounded radius-edge and aspect ratio, which renders our technique suitable for Finite Element analysis. A notable feature of our implementation is speed and scalability. The single-threaded performance of our 3D code is faster than the state of the art open source meshing tools. Experimental evaluation shows a more than 82% weak scaling efficiency for up to 144 cores, reaching a rate of more than 14.3 million elements per second. This is the first 3D parallel Delaunay refinement method to achieve such a performance, on either distributed or shared-memory architectures. Lastly, this dissertation is the first to develop and examine the sequential and parallel high-quality and fidelity meshing of general space-time 4D multi-tissue domains

Ultra-high speed particle image velocimetry on drop-on-demand jetting

Se estudian los conceptos claves de fortuna, virtù y gloria, con elobjeto de buscar algunas aclaraciones sobre el problema de larelación entre moral y política en Maquiavelo. La virtù maquiavelianaincluye, sin duda, múltiples componentes de energía, talento,pero no está completamente exenta de elementos morales. Delconcepto de fortuna podemos extraer una serie de criterios parala acción política que forman parte de la virtù maquiaveliana, lacual, si bien se aleja de la moral cristiana vigente en su tiempo, lesitúa dentro de una moral republicana. Del examen del conceptode gloria descubrimos cómo la gloria es recompensa de la virtù,pero no se reconoce a cualquier político que tiene éxito en susempresas, sino solo a quien salva a la patria y lo hace con mediosbenignos.</jats:p

Simulation of ultrasonic lamb wave generation, propagation and detection for an air coupled robotic scanner

A computer simulator, to facilitate the design and assessment of a reconfigurable, air-coupled ultrasonic scanner is described and evaluated. The specific scanning system comprises a team of remote sensing agents, in the form of miniature robotic platforms that can reposition non-contact Lamb wave transducers over a plate type of structure, for the purpose of non-destructive evaluation (NDE). The overall objective is to implement reconfigurable array scanning, where transmission and reception are facilitated by different sensing agents which can be organised in a variety of pulse-echo and pitch-catch configurations, with guided waves used to generate data in the form of 2-D and 3-D images. The ability to reconfigure the scanner adaptively requires an understanding of the ultrasonic wave generation, its propagation and interaction with potential defects and boundaries. Transducer behaviour has been simulated using a linear systems approximation, with wave propagation in the structure modelled using the local interaction simulation approach (LISA). Integration of the linear systems and LISA approaches are validated for use in Lamb wave scanning by comparison with both analytic techniques and more computationally intensive commercial finite element/difference codes. Starting with fundamental dispersion data, the paper goes on to describe the simulation of wave propagation and the subsequent interaction with artificial defects and plate boundaries, before presenting a theoretical image obtained from a team of sensing agents based on the current generation of sensors and instrumentation

Phenomenological model of diffuse global and regional atrophy using finite-element methods

The main goal of this work is the generation of ground-truth data for the validation of atrophy measurement techniques, commonly used in the study of neurodegenerative diseases such as dementia. Several techniques have been used to measure atrophy in cross-sectional and longitudinal studies, but it is extremely difficult to compare their performance since they have been applied to different patient populations. Furthermore, assessment of performance based on phantom measurements or simple scaled images overestimates these techniques' ability to capture the complexity of neurodegeneration of the human brain. We propose a method for atrophy simulation in structural magnetic resonance (MR) images based on finite-element methods. The method produces cohorts of brain images with known change that is physically and clinically plausible, providing data for objective evaluation of atrophy measurement techniques. Atrophy is simulated in different tissue compartments or in different neuroanatomical structures with a phenomenological model. This model of diffuse global and regional atrophy is based on volumetric measurements such as the brain or the hippocampus, from patients with known disease and guided by clinical knowledge of the relative pathological involvement of regions and tissues. The consequent biomechanical readjustment of structures is modelled using conventional physics-based techniques based on biomechanical tissue properties and simulating plausible tissue deformations with finite-element methods. A thermoelastic model of tissue deformation is employed, controlling the rate of progression of atrophy by means of a set of thermal coefficients, each one corresponding to a different type of tissue. Tissue characterization is performed by means of the meshing of a labelled brain atlas, creating a reference volumetric mesh that will be introduced to a finite-element solver to create the simulated deformations. Preliminary work on the simulation of acquisition artefa- - cts is also presented. Cross-sectional and

Developing Guidelines for Two-Dimensional Model Review and Acceptance

Two independent modelers ran two hydraulic models, SRH-2D and HEC-RAS 2D. The models were applied to the Lakina River (MP 44 McCarthy Road) and to Quartz Creek (MP 0.7 Quartz Creek Road), which approximately represent straight and bend flow conditions, respectively. We compared the results, including water depth, depth averaged velocity, and bed shear stress, from the two models for both modelers. We found that the extent and density of survey data were insufficient for Quartz Creek. Neither model was calibrated due to the lack of basic field data (i.e., discharge, water surface elevation, and sediment characteristics). Consequently, we were unable to draw any conclusion about the accuracy of the models. Concerning the time step and the equations used (simplified or full) to solve the momentum equation in the HEC-RAS 2D model, we found that the minimum time step allowed by the model must be used if the diffusion wave equation is used in the simulations. A greater time step can be used if the full momentum equation is used in the simulations. We developed a set of guidelines for reviewing model results, and developed and provided a two-day training workshop on the two models for ADOT&PF hydraulic engineers

Prediction of risk of fracture in the tibia due to altered bone mineral density distribution resulting from disuse : a finite element study

The disuse-related bone loss that results from immobilisation following injury shares characteristics with osteoporosis in postmenopausal women and the aged, with decreases in bone mineral density (BMD) leading to weakening of the bone and increased risk of fracture. The aim of the study was to use the finite element method to: (i) calculate the mechanical response of the tibia under mechanical load and (ii) estimate the risk of fracture; comparing between two groups, an able bodied (AB) group and spinal cord injury (SCI) patients group suffering from varying degree of bone loss. The tibiae of eight male subjects with chronic SCI and those of four able-bodied (AB) age-matched controls were scanned using multi-slice peripheral Quantitative Computed Tomography. Images were used to develop full three-dimensional models of the tibiae in Mimics (Materialise) and exported into Abaqus (Simulia) for calculation of stress distribution and fracture risk in response to specified loading conditions – compression, bending and torsion. The percentage of elements that exceeded a calculated value of the ultimate stress provided an estimate of the risk of fracture for each subject, which differed between SCI subjects and their controls. The differences in BMD distribution along the tibia in different subjects resulted in different regions of the bone being at high risk of fracture under set loading conditions, illustrating the benefit of creating individual material distribution models. A predictive tool can be developed based on these models, to enable clinicians to estimate the amount of loading that can be safely allowed onto the skeletal frame of individual patients who suffer from extensive musculoskeletal degeneration (including SCI, multiple sclerosis and the ageing population). The ultimate aim would be to reduce fracture occurrence in these vulnerable groups

Tetrahedral Image-to-Mesh Conversion Software for Anatomic Modeling of Arteriovenous Malformations

We describe a new implementation of an adaptive multi-tissue tetrahedral mesh generator targeting anatomic modeling of Arteriovenous Malformation (AVM) for surgical simulations. Our method, initially constructs an adaptive Body-Centered Cubic (BCC) mesh of high quality elements. Then, it deforms the mesh surfaces to their corresponding physical image boundaries, hence, improving the mesh fidelity and smoothness. Our deformation scheme, which builds upon the ITK toolkit, is based on the concept of energy minimization, and relies on a multi-material point-based registration. It uses non-connectivity patterns to implicitly control the number of the extracted feature points needed for the registration, and thus, adjusts the trade-off between the achieved mesh fidelity and the deformation speed. While many medical imaging applications require robust mesh generation, there are few codes available to the public. We compare our implementation with two similar open-source image-to-mesh conversion codes: (1) Cleaver from US, and (2) CGAL from EU. Our evaluation is based on five isotropic/anisotropic segmented images, and relies on metrics like geometric & topologic fidelity, mesh quality, gradation and smoothness. The implementation we describe is open- source and it will be available within: (i) the 3D Slicer package for visualization and image analysis from Harvard Medical School, and (ii) an interactive simulator for neurosurgical procedures involving vasculature using SOFA, a framework for real-time medical simulation developed by INRIA

Parallel Anisotropic Unstructured Grid Adaptation

Computational Fluid Dynamics (CFD) has become critical to the design and analysis of aerospace vehicles. Parallel grid adaptation that resolves multiple scales with anisotropy is identified as one of the challenges in the CFD Vision 2030 Study to increase the capacity and capability of CFD simulation. The Study also cautions that computer architectures are undergoing a radical change and dramatic increases in algorithm concurrency will be required to exploit full performance. This paper reviews four different methods to parallel anisotropic grid generation. They cover both ends of the spectrum: (i) using existing state-of-the-art software optimized for a single core and modifying it for parallel platforms and (ii) designing and implementing scalable software with incomplete, but rapidly maturating functionality. A brief overview for each grid adaptation system is presented in the context of a telescopic approach for multilevel concurrency. These methods employ different approaches to enable parallel execution, which provides a unique opportunity to illustrate the relative behavior of each approach. Qualitative and quantitative metric evaluations are used to draw lessons for future developments in this critical area for parallel CFD simulation

Multiscale Finite-Difference-Diffusion-Monte-Carlo Method for Simulating Dendritic Solidification

We present a novel hybrid computational method to simulate accurately dendritic solidification in the low undercooling limit where the dendrite tip radius is one or more orders of magnitude smaller than the characteristic spatial scale of variation of the surrounding thermal or solutal diffusion field. The first key feature of this method is an efficient multiscale diffusion Monte-Carlo (DMC) algorithm which allows off-lattice random walkers to take longer and concomitantly rarer steps with increasing distance away from the solid-liquid interface. As a result, the computational cost of evolving the large scale diffusion field becomes insignificant when compared to that of calculating the interface evolution. The second key feature is that random walks are only permitted outside of a thin liquid layer surrounding the interface. Inside this layer and in the solid, the diffusion equation is solved using a standard finite-difference algorithm that is interfaced with the DMC algorithm using the local conservation law for the diffusing quantity. Here we combine this algorithm with a previously developed phase-field formulation of the interface dynamics and demonstrate that it can accurately simulate three-dimensional dendritic growth in a previously unreachable range of low undercoolings that is of direct experimental relevance.Comment: RevTeX, 16 pages, 10 eps figures, submitted to J. Comp. Phy