Segmentation and grid generation for numerical simulations of vad connections with patient-specific data

Objectives: We are interested in the numerical simulation of the anastomotic region comprised between outflow canula of LVAD and the aorta. Segmenta¬tion, geometry reconstruction and grid generation from patient-specific data remain an issue because of the variable quality of DICOM images, in particular CT-scan (e.g. metallic noise of the device, non-aortic contrast phase). We pro¬pose a general framework to overcome this problem and create suitable grids for numerical simulations.Methods: Preliminary treatment of images is performed by reducing the level window and enhancing the contrast of the greyscale image using contrast-limited adaptive histogram equalization. A gradient anisotropic diffusion filter is applied to reduce the noise. Then, watershed segmentation algorithms and mathematical morphology filters allow reconstructing the patient geometry. This is done using the InsightToolKit library (www.itk.org). Finally the Vascular Model¬ing ToolKit (www.vmtk.org) and gmsh (www.geuz.org/gmsh) are used to create the meshes for the fluid (blood) and structure (arterial wall, outflow canula) and to a priori identify the boundary layers. The method is tested on five different patients with left ventricular assistance and who underwent a CT-scan exam.Results: This method produced good results in four patients. The anastomosis area is recovered and the generated grids are suitable for numerical simulations. In one patient the method failed to produce a good segmentation because of the small dimension of the aortic arch with respect to the image resolution.Conclusions: The described framework allows the use of data that could not be otherwise segmented by standard automatic segmentation tools. In particular the computational grids that have been generated are suitable for simulations that take into account fluid-structure interactions. Finally the presented method features a good reproducibility and fast application.

Numerical grid generation techniques : proceedings of a workshop held at NASA Langley Research Center, Hampton, Virginia, October 6-7, 1980 /

Includes bibliographical references.

Mesh Generation for Isogeometric Analysis with T-spline

[EN]The application of the Isogeometric Analysis (IA) with T-splines [1] demands a partition of the parametric space, C, in a tiling containing T-junctions denominated T-mesh. The T-splines are used both for the geometric modelization of the physical domain, D, and the basis of the numerical approximation. They have the advantage over the NURBS of allowing local refinement. In this work we propose a procedure to construct T-spline representations of complex domains in order to be applied to the resolution of elliptic PDE with IA. In precedent works [2, 3] we accomplished this task by using a tetrahedral parametrization…

Mortar-Techniken zur Behandlung von Grenzschichtproblemen

Elliptic differential equations, finite element method, mortar element method, streamline diffusion FEM, upwind method, numerical method, error estimate, interpolation operator, grid generation, adaptive refinement

Virtauslaskennankäyttö keskipakopumpun eroosiokulumisen mallintamisessa

Diplomityössä kehitetään malli tiheän kaksifaasivirtauksen aiheuttaman eroosiokulumi-sen mallintamiseksi, ratkaistaan virtauskenttä kahdessa erilaisessa keskipakopumpussa, sovelletaan kehitettyä mallia, sekävertaillaan mallin antamia tuloksia käytännön kokeissa saavutettuihin tuloksiin. Työssä on erityisenä mielenkiinnon kohteena savukaasupesurin pumppu. Työn alkuosa sisältää tarkemman kuvauksen savukaasupesurin toi-minnasta. Numeerinen ratkaisu ja laskentahilan generointi suoritetaan ANSYS CFX- ja Turbo-Grid-ohjelmistoilla. Laskennassa virtauksen Navier-Stokesin yhtälöt on aikakeskiarvo-tettu ja ratkaistu käyttäen kontrollitilavuusmenetelmää. Tiheiden kaksifaasivirtausten eroosiokulumista on mallinnettu tekijän kehittämällä mallilla, jonka käytännön toteutus-ta ei kuitenkaan saateta julkiseksi, koodin kehittämisessä käytetty teoria on kuitenkin esitetty työssä. Työn piirissä tehtiin myös kulumiskokeita Sulzer Pumps Finlandin Karhulan tehtailla, sekä vertailtiin simuloinnin tuloksia aikaisemmissa kokeissa saavutettuihin. Koejärjestelyt kuvataan työssä.

Numerical and ExperimentalModelling of Gas Flow and Heat Transfer in the Air Gap of An Electric Machine

This work deals with the cooling of high-speed electric machines, such as motors and generators, through an air gap. It consists of numerical and experimental modelling of gas flow and heat transfer in an annular channel. Velocity and temperature profiles are modelled in the air gap of a high-speed testmachine. Local and mean heat transfer coefficients and total friction coefficients are attained for a smooth rotor-stator combination at a large velocity range. The aim is to solve the heat transfer numerically and experimentally. The FINFLO software, developed at Helsinki University of Technology, has been used in the flow solution, and the commercial IGG and Field view programs for the grid generation and post processing. The annular channel is discretized as a sector mesh. Calculation is performed with constant mass flow rate on six rotational speeds. The effect of turbulence is calculated using three turbulence models. The friction coefficient and velocity factor are attained via total friction power. The first part of experimental section consists of finding the proper sensors and calibrating them in a straight pipe. After preliminary tests, a RdF-sensor is glued on the walls of stator and rotor surfaces. Telemetry is needed to be able to measure the heat transfer coefficients at the rotor. The mean heat transfer coefficients are measured in a test machine on four cooling air mass flow rates at a wide Couette Reynolds number range. The calculated values concerning the friction and heat transfer coefficients are compared with measured and semi-empirical data. Heat is transferred from the hotter stator and rotor surfaces to the coolerair flow in the air gap, not from the rotor to the stator via the air gap, althought the stator temperature is lower than the rotor temperature. The calculatedfriction coefficients fits well with the semi-empirical equations and precedingmeasurements. On constant mass flow rate the rotor heat transfer coefficient attains a saturation point at a higher rotational speed, while the heat transfer coefficient of the stator grows uniformly. The magnitudes of the heat transfer coefficients are almost constant with different turbulence models. The calibrationof sensors in a straight pipe is only an advisory step in the selection process. Telemetry is tested in the pipe conditions and compared to the same measurements with a plain sensor. The magnitudes of the measured data and the data from the semi-empirical equation are higher for the heat transfer coefficients than thenumerical data considered on the velocity range. Friction and heat transfer coefficients are presented in a large velocity range in the report. The goals are reached acceptably using numerical and experimental research. The next challenge is to achieve results for grooved stator-rotor combinations. The work contains also results for an air gap with a grooved stator with 36 slots. The velocity field by the numerical method does not match in every respect the estimated flow mode. The absence of secondary Taylor vortices is evident when using time averagednumerical simulation.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.

Modelagem de distorções entre sad 69 e sirgas 2000 pelo método de shepard e grades regulares

Brazil follows the tendency of some countries to update and/or review their fundamental geodetic network. The adoption of geocentric referentials like SIRGAS 2000, the new official reference system of the Geodetic Brazilian System has been an advance. Changes in referential implies in coordinates changes on the network stations as well as the network geometry. To make use of analogical and digital products which are already known in the old referentials are necessary approaches to the coordinate conversion, which minimize the distortions between the used reference frames. This paper presents a distortion modeling approach between reference frames, based on distortion grid generation by using the Shepard's method. To analyze the approach some experiments were performed with the generation of a 1 degrees x1 degrees distortion grid to model the distortions between SAD 69 (1996) and SIRGAS (2000) frames. The results in the test stations were promising, with an average reduction of 50% in the RMS coordinates after the distortions modeling.

Modelagem de distorções entre realizações de referenciais geodésicos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

An Adaptive Quadrilateral Mesh in Curved Domains

An nonlinear elliptic system for generating adaptive quadrilateral meshes in curved domains is presented. The presented technique has been implemented in the C++ language with the help of the standard template library. The software package writes the converged meshes in the GMV and the Matlab formats. Grid generation is the first very important step for numerically solving partial differential equations. Thus, the presented C++ grid generator is extremely important to the computational science community.

Load-balancing for parallel adaptive unstructured grids

A parallel method for the dynamic partitioning of unstructured meshes is outlined. The method includes diffusive load-balancing techniques and an iterative optimisation technique known as relative gain optimisationwhich both balances theworkload and attempts to minimise the interprocessor communications overhead. It can also optionally include amultilevel strategy. Experiments on a series of adaptively refined meshes indicate that the algorithmprovides partitions of an equivalent or higher quality to static partitioners (which do not reuse the existing partition) and much more rapidly. Perhaps more importantly, the algorithm results in only a small fraction of the amount of data migration compared to the static partitioners.

Triangular mesh parameterization with trimmed surfaces

Given a 2manifold triangular mesh \(M \subset {\mathbb {R}}^3\), with border, a parameterization of \(M\) is a FACE or trimmed surface \(F=\{S,L_0,\ldots, L_m\}\) -- \(F\) is a connected subset or region of a parametric surface \(S\), bounded by a set of LOOPs \(L_0,\ldots ,L_m\) such that each \(L_i \subset S\) is a closed 1manifold having no intersection with the other \(L_j\) LOOPs -- The parametric surface \(S\) is a statistical fit of the mesh \(M\) -- \(L_0\) is the outermost LOOP bounding \(F\) and \(L_i\) is the LOOP of the ith hole in \(F\) (if any) -- The problem of parameterizing triangular meshes is relevant for reverse engineering, tool path planning, feature detection, redesign, etc -- Stateofart mesh procedures parameterize a rectangular mesh \(M\) -- To improve such procedures, we report here the implementation of an algorithm which parameterizes meshes \(M\) presenting holes and concavities -- We synthesize a parametric surface \(S \subset {\mathbb {R}}^3\) which approximates a superset of the mesh \(M\) -- Then, we compute a set of LOOPs trimming \(S\), and therefore completing the FACE \(F=\ {S,L_0,\ldots ,L_m\}\) -- Our algorithm gives satisfactory results for \(M\) having low Gaussian curvature (i.e., \(M\) being quasi-developable or developable) -- This assumption is a reasonable one, since \(M\) is the product of manifold segmentation preprocessing -- Our algorithm computes: (1) a manifold learning mapping \(\phi : M \rightarrow U \subset {\mathbb {R}}^2\), (2) an inverse mapping \(S: W \subset {\mathbb {R}}^2 \rightarrow {\mathbb {R}}^3\), with \ (W\) being a rectangular grid containing and surpassing \(U\) -- To compute \(\phi\) we test IsoMap, Laplacian Eigenmaps and Hessian local linear embedding (best results with HLLE) -- For the back mapping (NURBS) \(S\) the crucial step is to find a control polyhedron \(P\), which is an extrapolation of \(M\) -- We calculate \(P\) by extrapolating radial basis functions that interpolate points inside \(\phi (M)\) -- We successfully test our implementation with several datasets presenting concavities, holes, and are extremely nondevelopable -- Ongoing work is being devoted to manifold segmentation which facilitates mesh parameterization

Using demand response to deal with unexpected low wind power generation in the context of smart grid

Demand response is assumed an essential resource to fully achieve the smart grids operating benefits, namely in the context of competitive markets. Some advantages of Demand Response (DR) programs and of smart grids can only be achieved through the implementation of Real Time Pricing (RTP). The integration of the expected increasing amounts of distributed energy resources, as well as new players, requires new approaches for the changing operation of power systems. The methodology proposed aims the minimization of the operation costs in a smart grid operated by a virtual power player. It is especially useful when actual and day ahead wind forecast differ significantly. When facing lower wind power generation than expected, RTP is used in order to minimize the impacts of such wind availability change. The proposed model application is here illustrated using the scenario of a special wind availability reduction day in the Portuguese power system (8th February 2012).

Improving the grid operation and reliability cost of distribution systems with dispersed generation

In this work, a mathematical model to analyze the impact of the installation and operation of dispersed generation units in power distribution systems is proposed. The main focus is to determine the trade-off between the reliability and operational costs of distribution networks when the operation of isolated areas is allowed. In order to increase the system operator revenue, an optimal power flow makes use of the different energy prices offered by the dispersed generation connected to the grid. Simultaneously, the type and location of the protective devices initially installed on the protection system are reconfigured in order to minimize the interruption and expenditure of adjusting the protection system to conditions imposed by the operation of dispersed units. The interruption cost regards the unsupplied energy to customers in secure systems but affected by the normal tripping of protective devices. Therefore, the tripping of fuses, reclosers, and overcurrent relays aims to protect the system against both temporary and permanent fault types. Additionally, in order to reduce the average duration of the system interruption experienced by customers, the isolated operation of dispersed generation is allowed by installing directional overcurrent relays with synchronized reclose capabilities. A 135-bus real distribution system is used in order to show the advantages of using the mathematical model proposed. © 1969-2012 IEEE.