Kainmueller, Dagmar
2014-01-01
? Segmentation of anatomical structures in medical image data is an essential task in clinical practice. Dagmar Kainmueller introduces methods for accurate fully automatic segmentation of anatomical structures in 3D medical image data. The author's core methodological contribution is a novel deformation model that overcomes limitations of state-of-the-art Deformable Surface approaches, hence allowing for accurate segmentation of tip- and ridge-shaped features of anatomical structures. As for practical contributions, she proposes application-specific segmentation pipelines for a range of anatom
Automated hexahedral meshing of anatomic structures using deformable registration.
Grosland, Nicole M; Bafna, Ritesh; Magnotta, Vincent A
2009-02-01
This work introduces a novel method of automating the process of patient-specific finite element (FE) model development using a mapped mesh technique. The objective is to map a predefined mesh (template) of high quality directly onto a new bony surface (target) definition, thereby yielding a similar mesh with minimal user interaction. To bring the template mesh into correspondence with the target surface, a deformable registration technique based on the FE method has been adopted. The procedure has been made hierarchical allowing several levels of mesh refinement to be used, thus reducing the time required to achieve a solution. Our initial efforts have focused on the phalanx bones of the human hand. Mesh quality metrics, such as element volume and distortion were evaluated. Furthermore, the distance between the target surface and the final mapped mesh were measured. The results have satisfactorily proven the applicability of the proposed method. PMID:18688764
Lenz, F; Doll, S; Sohn, C; Brocker, K A
2013-10-01
Purpose: Polypropylene mesh implants are frequently used for pelvic floor reconstruction in women. Yet they vary in size and fixation. The purpose of this study is to compare four mesh products with regard to their anatomical positioning and functionality within the pelvic floor, to determine whether each mesh fits equally well in a female cadaver. Methods: One female pelvis was dissected, opening the retropubic space exposing the endopelvic fascia and demonstrating the arcus tendineus fasciae pelvis (ATFP). Anatomical parameters were measured before and after implanting four meshes via the transobturator approach. Results: The anterior fixation of the ATFP was found to be 5 mm lateral to the symphysis in this cadaver. The endopelvic fascia covered 54.6 cm(2). The obturator nerve was located 35 mm from the white line. The distance of the proximal and lateral points of mesh fixation from the ischial spine or ATFP varied from 0 to 25 mm. The meshes varied in size and anatomical positioning. Conclusion: These observations demonstrate the necessity of developing optimally sized meshes and appropriate introducer techniques that can provide sufficient vaginal support. Surgeons, furthermore, need profound knowledge of anatomy, the patient's pelvic floor defect and the meshes available on the market. PMID:24771893
Accurate solution of the Dirac equation on Lagrange meshes
Baye, Daniel; Godefroid, Michel
2014-01-01
The Lagrange-mesh method is an approximate variational method taking the form of equations on a grid because of the use of a Gauss quadrature approximation. With a basis of Lagrange functions involving associated Laguerre polynomials related to the Gauss quadrature, the method is applied to the Dirac equation. The potential may possess a $1/r$ singularity. For hydrogenic atoms, numerically exact energies and wave functions are obtained with small numbers $n+1$ of mesh points, where $n$ is the principal quantum number. Numerically exact mean values of powers $-2$ to 3 of the radial coordinate $r$ can also be obtained with $n+2$ mesh points. For the Yukawa potential, a 15-digit agreement with benchmark energies of the literature is obtained with 50 mesh points or less.
Accurately simulating anisotropic thermal conduction on a moving mesh
Kannan, Rahul; Pakmor, Rüdiger; Marinacci, Federico; Vogelsberger, Mark
2015-01-01
We present a novel implementation of an extremum preserving anisotropic diffusion solver for thermal conduction on the unstructured moving Voronoi mesh of the AREPO code. The method relies on splitting the one-sided facet fluxes into normal and oblique components, with the oblique fluxes being limited such that the total flux is both locally conservative and extremum preserving. The approach makes use of harmonic averaging points and a simple, robust interpolation scheme that works well for strong heterogeneous and anisotropic diffusion problems. Moreover, the required discretisation stencil is small. Efficient fully implicit and semi-implicit time integration schemes are also implemented. We perform several numerical tests that evaluate the stability and accuracy of the scheme, including applications such as point explosions with heat conduction and calculations of convective instabilities in conducting plasmas. The new implementation is suitable for studying important astrophysical phenomena, such as the co...
Accurately simulating anisotropic thermal conduction on a moving mesh
Kannan, Rahul; Springel, Volker; Pakmor, Rüdiger; Marinacci, Federico; Vogelsberger, Mark
2016-05-01
We present a novel implementation of an extremum preserving anisotropic diffusion solver for thermal conduction on the unstructured moving Voronoi mesh of the AREPO code. The method relies on splitting the one-sided facet fluxes into normal and oblique components, with the oblique fluxes being limited such that the total flux is both locally conservative and extremum preserving. The approach makes use of harmonic averaging points and a simple, robust interpolation scheme that works well for strong heterogeneous and anisotropic diffusion problems. Moreover, the required discretization stencil is small. Efficient fully implicit and semi-implicit time integration schemes are also implemented. We perform several numerical tests that evaluate the stability and accuracy of the scheme, including applications such as point explosions with heat conduction and calculations of convective instabilities in conducting plasmas. The new implementation is suitable for studying important astrophysical phenomena, such as the conductive heat transport in galaxy clusters, the evolution of supernova remnants, or the distribution of heat from black hole-driven jets into the intracluster medium.
Rumple, Christopher; Krane, Michael; Richter, Joseph; Craven, Brent
2013-11-01
The mammalian nose is a multi-purpose organ that houses a convoluted airway labyrinth responsible for respiratory air conditioning, filtering of environmental contaminants, and chemical sensing. Because of the complexity of the nasal cavity, the anatomy and function of these upper airways remain poorly understood in most mammals. However, recent advances in high-resolution medical imaging, computational modeling, and experimental flow measurement techniques are now permitting the study of respiratory airflow and olfactory transport phenomena in anatomically-accurate reconstructions of the nasal cavity. Here, we focus on efforts to manufacture an anatomically-accurate transparent model for stereoscopic particle image velocimetry (SPIV) measurements. Challenges in the design and manufacture of an index-matched anatomical model are addressed. PIV measurements are presented, which are used to validate concurrent computational fluid dynamics (CFD) simulations of mammalian nasal airflow. Supported by the National Science Foundation.
Augustin, Christoph M.; Neic, Aurel; Liebmann, Manfred; Prassl, Anton J.; Niederer, Steven A.; Haase, Gundolf; Plank, Gernot
2016-01-01
Electromechanical (EM) models of the heart have been used successfully to study fundamental mechanisms underlying a heart beat in health and disease. However, in all modeling studies reported so far numerous simplifications were made in terms of representing biophysical details of cellular function and its heterogeneity, gross anatomy and tissue microstructure, as well as the bidirectional coupling between electrophysiology (EP) and tissue distension. One limiting factor is the employed spatial discretization methods which are not sufficiently flexible to accommodate complex geometries or resolve heterogeneities, but, even more importantly, the limited efficiency of the prevailing solver techniques which is not sufficiently scalable to deal with the incurring increase in degrees of freedom (DOF) when modeling cardiac electromechanics at high spatio-temporal resolution. This study reports on the development of a novel methodology for solving the nonlinear equation of finite elasticity using human whole organ models of cardiac electromechanics, discretized at a high para-cellular resolution. Three patient-specific, anatomically accurate, whole heart EM models were reconstructed from magnetic resonance (MR) scans at resolutions of 220 μm, 440 μm and 880 μm, yielding meshes of approximately 184.6, 24.4 and 3.7 million tetrahedral elements and 95.9, 13.2 and 2.1 million displacement DOF, respectively. The same mesh was used for discretizing the governing equations of both electrophysiology (EP) and nonlinear elasticity. A novel algebraic multigrid (AMG) preconditioner for an iterative Krylov solver was developed to deal with the resulting computational load. The AMG preconditioner was designed under the primary objective of achieving favorable strong scaling characteristics for both setup and solution runtimes, as this is key for exploiting current high performance computing hardware. Benchmark results using the 220 μm, 440 μm and 880 μm meshes demonstrate
Creation of Anatomically Accurate Computer-Aided Design (CAD) Solid Models from Medical Images
Stewart, John E.; Graham, R. Scott; Samareh, Jamshid A.; Oberlander, Eric J.; Broaddus, William C.
1999-01-01
Most surgical instrumentation and implants used in the world today are designed with sophisticated Computer-Aided Design (CAD)/Computer-Aided Manufacturing (CAM) software. This software automates the mechanical development of a product from its conceptual design through manufacturing. CAD software also provides a means of manipulating solid models prior to Finite Element Modeling (FEM). Few surgical products are designed in conjunction with accurate CAD models of human anatomy because of the difficulty with which these models are created. We have developed a novel technique that creates anatomically accurate, patient specific CAD solids from medical images in a matter of minutes.
Rumple, C.; Richter, J.; Craven, B. A.; Krane, M.
2012-11-01
A summary of the research being carried out by our multidisciplinary team to better understand the form and function of the nose in different mammalian species that include humans, carnivores, ungulates, rodents, and marine animals will be presented. The mammalian nose houses a convoluted airway labyrinth, where two hallmark features of mammals occur, endothermy and olfaction. Because of the complexity of the nasal cavity, the anatomy and function of these upper airways remain poorly understood in most mammals. However, recent advances in high-resolution medical imaging, computational modeling, and experimental flow measurement techniques are now permitting the study of airflow and respiratory and olfactory transport phenomena in anatomically-accurate reconstructions of the nasal cavity. Here, we focus on efforts to manufacture transparent, anatomically-accurate models for stereo particle image velocimetry (SPIV) measurements of nasal airflow. Challenges in the design and manufacture of index-matched anatomical models are addressed and preliminary SPIV measurements are presented. Such measurements will constitute a validation database for concurrent computational fluid dynamics (CFD) simulations of mammalian respiration and olfaction. Supported by the National Science Foundation.
Parallelization of a high-order accurate unstructured mesh finite-volume solver
A generic solver eliminates the need to write new finite-volume codes for each type of physics. By separating the physics from the numerics of the solver, a modular design is achieved. New physics modules can easily be written with a minimal knowledge of the finite-volume method. A parallel solver allows simulation of complex physics on intricate domains in a timely manner by using numerous processors simultaneously. In this paper we describe the steps needed to adapt a high-order accurate unstructured mesh generic finite-volume solver to a parallel architecture. A message-passing approach is used which allows the solver to operate on a distributed memory system, such as a cluster of workstations. The reconstruction stencil is determined at the preprocessing stage and an appropriate parallel data structure for the solution is formed. Fluxes for faces on the partition boundary are evaluated by communicating the reconstruction coefficients to the adjacent processor. Good performance scalability is achieved for second and fourth-order accurate solutions on cell and vertex centered meshes. (authors)
Lee, Dongwook
2013-01-01
In this paper, we extend the unsplit staggered mesh scheme (USM) for 2D magnetohydrodynamics (MHD) (Lee and Deane, 2009) to a full 3D MHD scheme. The scheme is a finite-volume Godunov method consisting of a constrained transport (CT) method and an efficient and accurate single-step, directionally unsplit multidimensional data reconstruction-evolution algorithm, which extends the original 2D corner transport upwind (CTU) method (Colella, 1990). We present two types of data reconstruction-evolution algorithms for 3D: (1) a reduced CTU scheme and (2) a full CTU scheme. The reduced 3D CTU scheme is a variant of a simple 3D extension of the 2D CTU method by Colella (1990) and is considered as a direct extension from the 2D USM scheme. The full 3D CTU scheme is our primary 3D solver which includes all multidimensional cross-derivative terms for stability. The latter method is logically analogous to the 3D unsplit CTU method by Saltzman. The major novelties in our algorithms are twofold. First, we extend the reduced...
O'Halloran, M.; Lohfeld, S.; Ruvio, G.; Browne, J.; Krewer, F.; Ribeiro, C. O.; Inacio Pita, V. C.; Conceicao, R. C.; Jones, E.; Glavin, M.
2014-05-01
Breast cancer is one of the most common cancers in women. In the United States alone, it accounts for 31% of new cancer cases, and is second only to lung cancer as the leading cause of deaths in American women. More than 184,000 new cases of breast cancer are diagnosed each year resulting in approximately 41,000 deaths. Early detection and intervention is one of the most significant factors in improving the survival rates and quality of life experienced by breast cancer sufferers, since this is the time when treatment is most effective. One of the most promising breast imaging modalities is microwave imaging. The physical basis of active microwave imaging is the dielectric contrast between normal and malignant breast tissue that exists at microwave frequencies. The dielectric contrast is mainly due to the increased water content present in the cancerous tissue. Microwave imaging is non-ionizing, does not require breast compression, is less invasive than X-ray mammography, and is potentially low cost. While several prototype microwave breast imaging systems are currently in various stages of development, the design and fabrication of anatomically and dielectrically representative breast phantoms to evaluate these systems is often problematic. While some existing phantoms are composed of dielectrically representative materials, they rarely accurately represent the shape and size of a typical breast. Conversely, several phantoms have been developed to accurately model the shape of the human breast, but have inappropriate dielectric properties. This study will brie y review existing phantoms before describing the development of a more accurate and practical breast phantom for the evaluation of microwave breast imaging systems.
The Quick Measure of a Nurbs Surface Curvature for Accurate Triangular Meshing
Kniat Aleksander
2014-04-01
Full Text Available NURBS surfaces are the most widely used surfaces for three-dimensional models in CAD/ CAE programs. When a model for FEM calculation is prepared with a CAD program it is inevitable to mesh it finally. There are many algorithms for meshing planar regions. Some of them may be used for meshing surfaces but it is necessary to take the curvature of the surface under consideration to avoid poor quality mesh. The mesh must be denser in the curved regions of the surface. In this paper, instead of analysing a surface curvature, the method to assess how close is a mesh triangle to the surface to which its vertices belong, is presented. The distance between a mesh triangle and a parallel tangent plane through a point on a surface is the measure of the triangle quality. Finding the surface point whose projection is located inside the mesh triangle and which is the tangency point to the plane parallel to this triangle is an optimization problem. Mathematical description of the problem and the algorithm to find its solution are also presented in the paper.
Streett, Craig L.; Venkatachari, Balaji Shankar; Chang, Chau-Lyan; Friedlander, David J.; Wang, Xiao-Yen; Chang, Sin-Chung
2016-01-01
Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.
Background: Numbering of the thoracic spine on MRI can be tedious if C2 and L5-S1 are not included and may lead to errors in lesion level. Purpose: To determine whether anatomic landmarks or external markers are reliable as an aid for accurate numbering of thoracic vertebrae on MRI. Material and Methods: Sixty-seven thoracic spine MR studies of 67 patients (30 males, 37 females, age range 18-83 years) were studied, composed of 52 consecutive MR studies and an additional 15 MRI in which vitamin E markers were placed over the skin. In the 52 thoracic MR examinations potential numbering aids such as the level of the sternal apex, pulmonary artery, aortic arch, and osseous or disc abnormalities were numbered on both cervical localizer (standard of reference) and thoracic sagittal images. The additional 15 examinations in which vitamin E markers were placed over the skin were evaluated for consistency in the level of the markers on different sequences in the same exam. Results: The sternal apex level ranged from T2 to T5 [T3 in 28/51 patients (55%), T2 in 10/51 (20%)]. The aortic arch level ranged from T2 to T4 [T4 in 18/48 (38%) and T3 in 17 (35%)]. Pulmonary artery level ranged from T4 to T6-7 disc [T5 in 20/52 patients (38%) and T6 in 14/52 (27%)]. In 3 of 12 patients who had abnormalities in a vertebral body or disc as definite point reference, the non-localizer image mislabelled the level. In 11/15 (73%) patients with vitamin E markers that were placed over the upper thoracic spine, the results showed consistency in the level of the markers in relation to the reference points or consistent inter-marker gap between the sequences. Conclusion: There are only two reliable ways to accurately define the levels if no landmarking feature is available on the magnet. The first is by including C2 in the thoracic sequence of a diagnostic quality, and the second is by using an abnormality in the discs or vertebral bodies as a point of reference
New approach based on tetrahedral-mesh geometry for accurate 4D Monte Carlo patient-dose calculation
In the present study, to achieve accurate 4D Monte Carlo dose calculation in radiation therapy, we devised a new approach that combines (1) modeling of the patient body using tetrahedral-mesh geometry based on the patient’s 4D CT data, (2) continuous movement/deformation of the tetrahedral patient model by interpolation of deformation vector fields acquired through deformable image registration, and (3) direct transportation of radiation particles during the movement and deformation of the tetrahedral patient model. The results of our feasibility study show that it is certainly possible to construct 4D patient models (= phantoms) with sufficient accuracy using the tetrahedral-mesh geometry and to directly transport radiation particles during continuous movement and deformation of the tetrahedral patient model. This new approach not only produces more accurate dose distribution in the patient but also replaces the current practice of using multiple 3D voxel phantoms and combining multiple dose distributions after Monte Carlo simulations. For routine clinical application of our new approach, the use of fast automatic segmentation algorithms is a must. In order to achieve, simultaneously, both dose accuracy and computation speed, the number of tetrahedrons for the lungs should be optimized. Although the current computation speed of our new 4D Monte Carlo simulation approach is slow (i.e. ∼40 times slower than that of the conventional dose accumulation approach), this problem is resolvable by developing, in Geant4, a dedicated navigation class optimized for particle transportation in tetrahedral-mesh geometry. (paper)
Background and purpose: Low contrast in the cone-beam computed tomography (CBCT) scans hampers fast online evaluation of interfractional changes in the lymph node position on a daily basis. In this study we have investigated whether high-contrast anatomical landmarks in the vicinity of the nodes may be used as surrogates for the lymph node positions. Materials and methods: Forty lung cancer patients were treated with an online CBCT-based setup strategy involving soft-tissue match on the primary tumor. One hundred and sixteen lymph nodes were delineated separately on the planning-CT scans and categorized according to the lymph node stations. Five anatomical landmarks were selected as surrogate structures and assigned to the individual nodes. In addition, the carina was delineated. Registrations between the planning-CT and the daily CBCTs were performed retrospectively and positional deviations between the nodes and the surrogate structures or the carina were registered. Results: The mean displacement between lymph nodes and surrogate structures was 1.6 mm with systematic/random errors of 0.7/0.7 mm, significantly smaller than the mean displacement between nodes and the carina. Conclusions: The position of the lymph nodes can be evaluated using selected anatomical landmarks on a daily basis using CBCT
GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
Hopkins, Philip F
2014-01-01
We present and study two new Lagrangian numerical methods for solving the equations of hydrodynamics, in a systematic comparison with moving-mesh, SPH, and non-moving grid methods. The new methods are designed to capture many advantages of both smoothed-particle hydrodynamics (SPH) and grid-based or adaptive mesh refinement (AMR) schemes. They are based on a kernel discretization of the volume coupled to a high-order matrix gradient estimator and a Riemann solver acting over the volume 'overlap.' We implement and test a parallel, second-order version of the method with coupled self-gravity & cosmological integration, in the code GIZMO: this maintains exact mass, energy and momentum conservation; exhibits superior angular momentum conservation compared to all other methods we study; does not require 'artificial diffusion' terms; and allows fluid elements to move with the flow so resolution is automatically adaptive. We consider a large suite of test problems, and find that on all problems the new methods a...
Hepburn, I; Chen, W; De Schutter, E
2016-08-01
Spatial stochastic molecular simulations in biology are limited by the intense computation required to track molecules in space either in a discrete time or discrete space framework, which has led to the development of parallel methods that can take advantage of the power of modern supercomputers in recent years. We systematically test suggested components of stochastic reaction-diffusion operator splitting in the literature and discuss their effects on accuracy. We introduce an operator splitting implementation for irregular meshes that enhances accuracy with minimal performance cost. We test a range of models in small-scale MPI simulations from simple diffusion models to realistic biological models and find that multi-dimensional geometry partitioning is an important consideration for optimum performance. We demonstrate performance gains of 1-3 orders of magnitude in the parallel implementation, with peak performance strongly dependent on model specification. PMID:27497550
Hepburn, I.; Chen, W.; De Schutter, E.
2016-08-01
Spatial stochastic molecular simulations in biology are limited by the intense computation required to track molecules in space either in a discrete time or discrete space framework, which has led to the development of parallel methods that can take advantage of the power of modern supercomputers in recent years. We systematically test suggested components of stochastic reaction-diffusion operator splitting in the literature and discuss their effects on accuracy. We introduce an operator splitting implementation for irregular meshes that enhances accuracy with minimal performance cost. We test a range of models in small-scale MPI simulations from simple diffusion models to realistic biological models and find that multi-dimensional geometry partitioning is an important consideration for optimum performance. We demonstrate performance gains of 1-3 orders of magnitude in the parallel implementation, with peak performance strongly dependent on model specification.
This paper proposes an alternative approach to enhance localization accuracy of MEG and EEG focal sources. The proposed approach assumes anatomically constrained spatio-temporal dipoles, initial positions of which are estimated from local peak positions of distributed sources obtained from a pre-execution of distributed source reconstruction. The positions of the dipoles are then adjusted on the cortical surface using a novel updating scheme named cortical surface scanning. The proposed approach has many advantages over the conventional ones: (1) as the cortical surface scanning algorithm uses spatio-temporal dipoles, it is robust with respect to noise; (2) it requires no a priori information on the numbers and initial locations of the activations; (3) as the locations of dipoles are restricted only on a tessellated cortical surface, it is physiologically more plausible than the conventional ECD model. To verify the proposed approach, it was applied to several realistic MEG/EEG simulations and practical experiments. From the several case studies, it is concluded that the anatomically constrained dipole adjustment (ANACONDA) approach will be a very promising technique to enhance accuracy of focal source localization which is essential in many clinical and neurological applications of MEG and EEG
Lei Zhang
2015-01-01
Full Text Available Background: Food and Drug Administration announcements have highlighted the standard rate of mesh-related complications. We aimed to report the short-term results and complications of tension-free polypropylene mesh (PROSIMA™ surgical repair of pelvic organ prolapse (POP using the standard category (C, timing (T, and site (S classification system. Methods: A prospective cohort study of 48 patients who underwent PROSIMA™ mesh kit-related surgical repairs were followed for two years at Peking Union Medical College Hospital. Recurrence was defined as symptomatic POP quantification (POP-Q Stage II or higher (leading edge ≥ −1 cm. The Patient Global Impression of Change Questionnaire, the Chinese version of the Pelvic Floor Impact Questionnaire short-form-7 and POP/Urinary Incontinence Sexual Questionnaire short-form-12 were used to evaluate the self-perception and sexual function of each patient. Mesh-related complications conformed to the International Urogynecological Association/International Continence Society joint terminology. The paired-sample t-test, one-way analysis of variance, Fisher′s exact test, Kaplan-Meier survival analysis and log-rank test were used to analyze data. Results: All patients were followed up for ≥12 months; 30 (62.5% patients completed the 24 months study. We observed a 93.8% (45/48 positive anatomical outcome rate at 12 months and 90.0% (27/30 at 24 months. Recurrence most frequently involved the anterior compartment (P < 0.05. Pelvic symptoms improved significantly from baseline (P < 0.05, although the patients′ impressions of change and sexual function were not satisfying. Vaginal complication was the main complication observed (35.4%, 17/48. The survival analysis did not identify any relationship between vaginal complication and anatomical recurrent prolapse (POP-Q ≥ Stage II (P = 0.653. Conclusions: Tension-free polypropylene mesh (PROSIMA™ -related surgical repair of POP has better short
Lei Zhang; Lan Zhu; Juan Chen; Tao Xu; Jing-He Lang
2015-01-01
Background:Food and Drug Administration announcements have highlighted the standard rate of mesh-related complications.We aimed to report the short-term results and complications of tension-free polypropylene mesh (PROSIMATM) surgical repair of pelvic organ prolapse (POP) using the standard category (C),timing (T),and site (S) classification system.Methods:A prospective cohort study of 48 patients who underwent PROSIMATM mesh kit-related surgical repairs were followed for two years at Peking Union Medical College Hospital.Recurrence was defined as symptomatic POP quantification (POP-Q) Stage Ⅱ or higher (leading edge ≥-1 cm).The Patient Global Impression of Change Questionnaire,the Chinese version of the Pelvic Floor Impact Questionnaire short-form-7 and POP/Urinary Incontinence Sexual Questionnaire short-form-12 were used to evaluate the self-perception and sexual function of each patient.Mesh-related complications conformed to the International Urogynecological Association/International Continence Society joint terminology.The paired-sample t-test,one-way analysis of variance,Fisher's exact test,Kaplan-Meier survival analysis and log-rank test were used to analyze data.Results:All patients were followed up for≥1 2 months; 30 (62.5％) patients completed the 24 months study.We observed a 93.8％ (45/48) positive anatomical outcome rate at 12 months and 90.0％ (27/30) at 24 months.Recurrence most frequently involved the anterior compartment (P ＜ 0.05).Pelvic symptoms improved significantly from baseline (P ＜ 0.05),although the patients' impressions of change and sexual function were not satisfying.Vaginal complication was the main complication observed (35.4％,17/48).The survival analysis did not identify any relationship between vaginal complication and anatomical recurrent prolapse (POP-Q ≥ Stage Ⅱ) (P =0.653).Conclusions:Tension-free polypropylene mesh (PROSIMATM)-related surgical repair of POP has better short-term anatomical outcomes at the
The aortic valve exhibits complex three-dimensional (3D) anatomy and heterogeneity essential for the long-term efficient biomechanical function. These are, however, challenging to mimic in de novo engineered living tissue valve strategies. We present a novel simultaneous 3D printing/photocrosslinking technique for rapidly engineering complex, heterogeneous aortic valve scaffolds. Native anatomic and axisymmetric aortic valve geometries (root wall and tri-leaflets) with 12–22 mm inner diameters (ID) were 3D printed with poly-ethylene glycol-diacrylate (PEG-DA) hydrogels (700 or 8000 MW) supplemented with alginate. 3D printing geometric accuracy was quantified and compared using Micro-CT. Porcine aortic valve interstitial cells (PAVIC) seeded scaffolds were cultured for up to 21 days. Results showed that blended PEG-DA scaffolds could achieve over tenfold range in elastic modulus (5.3±0.9 to 74.6±1.5 kPa). 3D printing times for valve conduits with mechanically contrasting hydrogels were optimized to 14 to 45 min, increasing linearly with conduit diameter. Larger printed valves had greater shape fidelity (93.3±2.6, 85.1±2.0 and 73.3±5.2% for 22, 17 and 12 mm ID porcine valves; 89.1±4.0, 84.1±5.6 and 66.6±5.2% for simplified valves). PAVIC seeded scaffolds maintained near 100% viability over 21 days. These results demonstrate that 3D hydrogel printing with controlled photocrosslinking can rapidly fabricate anatomical heterogeneous valve conduits that support cell engraftment. (paper)
ACCURATE BENDING STRENGTH ANALYSIS OF ASYMMETRIC GEARS USING THE NOVEL ES-PIM WITH TRIANGULAR MESH
G. Y. Zhang
2011-12-01
Full Text Available This paper extends the edge-based smoothed point interpolation method (ES-PIM to bending strength analysis of asymmetric gears with complex outlines. Five sets of asymmetric gears with pressure angles of 20°/20°, 25°/20°, 30°/20°, 35°/20° and 40°/20° were generated by a specially designed rack cutter. Four key factors, e.g. accuracy, convergence, the convergence rate and the computational efficiency of the present ES-PIM were checked in great detail on these five models, and the distributions of bending stresses at the fillet of the drive side were carefully investigated. The finite element method (FEM was also used to calculate the abovementioned factors to stress the advantages of ES-PIM. The numerical results indicate that ES-PIM can provide more efficient and accurate solutions in the stress field than the FEM, and is very suitable for stress analysis of complicated asymmetric gears.
There is a need for accurate surface area measurement of internal anatomical structures in order to define light dosimetry in adjunctive intraoperative photodynamic therapy (AIOPDT). The authors investigated whether computer-assisted triangulation of serial sections generated by computed tomography (CT) scanning can give an accurate assessment of the surface area of the walls of the true pelvis after anterior resection and before colorectal anastomosis. They show that the technique of paper density tessellation is an acceptable method of measuring the surface areas of phantom objects, with a maximum error of 0.5%, and is used as the gold standard. Computer-assisted triangulation of CT images of standard geometric objects and accurately-constructed pelvic phantoms gives a surface area assessment with a maximum error of 2.5% compared with the gold standard. The CT images of 20 patients' pelves have been analysed by computer-assisted triangulation and this shows the surface area of the walls varies from 143 cm2 to 392 cm2. (Author)
Bache, Steven T.; Juang, Titania; Belley, Matthew D. [Duke University Medical Physics Graduate Program, Durham, North Carolina 27705 (United States); Koontz, Bridget F.; Yoshizumi, Terry T.; Kirsch, David G.; Oldham, Mark, E-mail: mark.oldham@duke.edu [Duke University Medical Center, Durham, North Carolina 27710 (United States); Adamovics, John [Rider University, Lawrenceville, New Jersey 08648 (United States)
2015-02-15
Purpose: Sophisticated small animal irradiators, incorporating cone-beam-CT image-guidance, have recently been developed which enable exploration of the efficacy of advanced radiation treatments in the preclinical setting. Microstereotactic-body-radiation-therapy (microSBRT) is one technique of interest, utilizing field sizes in the range of 1–15 mm. Verification of the accuracy of microSBRT treatment delivery is challenging due to the lack of available methods to comprehensively measure dose distributions in representative phantoms with sufficiently high spatial resolution and in 3 dimensions (3D). This work introduces a potential solution in the form of anatomically accurate rodent-morphic 3D dosimeters compatible with ultrahigh resolution (0.3 mm{sup 3}) optical computed tomography (optical-CT) dose read-out. Methods: Rodent-morphic dosimeters were produced by 3D-printing molds of rodent anatomy directly from contours defined on x-ray CT data sets of rats and mice, and using these molds to create tissue-equivalent radiochromic 3D dosimeters from Presage. Anatomically accurate spines were incorporated into some dosimeters, by first 3D printing the spine mold, then forming a high-Z bone equivalent spine insert. This spine insert was then set inside the tissue equivalent body mold. The high-Z spinal insert enabled representative cone-beam CT IGRT targeting. On irradiation, a linear radiochromic change in optical-density occurs in the dosimeter, which is proportional to absorbed dose, and was read out using optical-CT in high-resolution (0.5 mm isotropic voxels). Optical-CT data were converted to absolute dose in two ways: (i) using a calibration curve derived from other Presage dosimeters from the same batch, and (ii) by independent measurement of calibrated dose at a point using a novel detector comprised of a yttrium oxide based nanocrystalline scintillator, with a submillimeter active length. A microSBRT spinal treatment was delivered consisting of a 180
Purpose: Sophisticated small animal irradiators, incorporating cone-beam-CT image-guidance, have recently been developed which enable exploration of the efficacy of advanced radiation treatments in the preclinical setting. Microstereotactic-body-radiation-therapy (microSBRT) is one technique of interest, utilizing field sizes in the range of 1–15 mm. Verification of the accuracy of microSBRT treatment delivery is challenging due to the lack of available methods to comprehensively measure dose distributions in representative phantoms with sufficiently high spatial resolution and in 3 dimensions (3D). This work introduces a potential solution in the form of anatomically accurate rodent-morphic 3D dosimeters compatible with ultrahigh resolution (0.3 mm3) optical computed tomography (optical-CT) dose read-out. Methods: Rodent-morphic dosimeters were produced by 3D-printing molds of rodent anatomy directly from contours defined on x-ray CT data sets of rats and mice, and using these molds to create tissue-equivalent radiochromic 3D dosimeters from Presage. Anatomically accurate spines were incorporated into some dosimeters, by first 3D printing the spine mold, then forming a high-Z bone equivalent spine insert. This spine insert was then set inside the tissue equivalent body mold. The high-Z spinal insert enabled representative cone-beam CT IGRT targeting. On irradiation, a linear radiochromic change in optical-density occurs in the dosimeter, which is proportional to absorbed dose, and was read out using optical-CT in high-resolution (0.5 mm isotropic voxels). Optical-CT data were converted to absolute dose in two ways: (i) using a calibration curve derived from other Presage dosimeters from the same batch, and (ii) by independent measurement of calibrated dose at a point using a novel detector comprised of a yttrium oxide based nanocrystalline scintillator, with a submillimeter active length. A microSBRT spinal treatment was delivered consisting of a 180
Zhang, Bin; Liang, Chunlei
2015-08-01
This paper presents a simple, efficient, and high-order accurate sliding-mesh interface approach to the spectral difference (SD) method. We demonstrate the approach by solving the two-dimensional compressible Navier-Stokes equations on quadrilateral grids. This approach is an extension of the straight mortar method originally designed for stationary domains [7,8]. Our sliding method creates curved dynamic mortars on sliding-mesh interfaces to couple rotating and stationary domains. On the nonconforming sliding-mesh interfaces, the related variables are first projected from cell faces to mortars to compute common fluxes, and then the common fluxes are projected back from the mortars to the cell faces to ensure conservation. To verify the spatial order of accuracy of the sliding-mesh spectral difference (SSD) method, both inviscid and viscous flow cases are tested. It is shown that the SSD method preserves the high-order accuracy of the SD method. Meanwhile, the SSD method is found to be very efficient in terms of computational cost. This novel sliding-mesh interface method is very suitable for parallel processing with domain decomposition. It can be applied to a wide range of problems, such as the hydrodynamics of marine propellers, the aerodynamics of rotorcraft, wind turbines, and oscillating wing power generators, etc.
Chang, Chau-Lyan; Venkatachari, Balaji Shankar; Cheng, Gary
2013-01-01
With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.
The clinical differentiation between typical idiopathic Parkinson's disease (IPD) and atypical parkinsonian disorders such as multiple system atrophy (MSA) is complicated by the presence of signs and symptoms common to both forms. The goal of this study was to re-evaluate the contribution of brain perfusion single-photon emission tomography (SPET) with anatomical standardisation and automated analysis in the differentiation of IPD and MSA. This was achieved by discriminant analysis in comparison with a large set of age- and gender-matched healthy volunteers. Technetium-99m ethyl cysteinate dimer SPET was performed on 140 subjects: 81 IPD patients (age 62.6±10.2 years; disease duration 11.0±6.4 years; 50 males/31 females), 15 MSA patients (61.5±9.2 years; disease duration 3.0±2.2 years; 9 males/6 females) and 44 age- and gender-matched healthy volunteers (age 59.2±11.9 years; 27 males/17 females). Patients were matched for severity (Hoehn and Yahr stage). Automated predefined volume of interest (VOI) analysis was carried out after anatomical standardisation. Stepwise discriminant analysis with cross-validation using the leave-one-out method was used to determine the subgroup of variables giving the highest accuracy for this differential diagnosis. Between MSA and IPD, the only regions with highly significant differences in uptake after Bonferroni correction were the putamen VOIs. Comparing MSA versus normals and IPD, with putamen VOI values as discriminating variables, cross-validated performance showed correct classification of MSA patients with a sensitivity of 73.3%, a specificity of 84% and an accuracy of 83.6%. Additional input from the right caudate head and the left prefrontal and left mesial temporal cortex allowed 100% discrimination even after cross-validation. Discrimination between the IPD group alone and healthy volunteers was accurate in 94% of the cases after cross-validation, with a sensitivity of 91.4% and a specificity of 100%. The three
Dumbser, Michael; Loubère, Raphaël
2016-08-01
In this paper we propose a simple, robust and accurate nonlinear a posteriori stabilization of the Discontinuous Galerkin (DG) finite element method for the solution of nonlinear hyperbolic PDE systems on unstructured triangular and tetrahedral meshes in two and three space dimensions. This novel a posteriori limiter, which has been recently proposed for the simple Cartesian grid case in [62], is able to resolve discontinuities at a sub-grid scale and is substantially extended here to general unstructured simplex meshes in 2D and 3D. It can be summarized as follows: At the beginning of each time step, an approximation of the local minimum and maximum of the discrete solution is computed for each cell, taking into account also the vertex neighbors of an element. Then, an unlimited discontinuous Galerkin scheme of approximation degree N is run for one time step to produce a so-called candidate solution. Subsequently, an a posteriori detection step checks the unlimited candidate solution at time t n + 1 for positivity, absence of floating point errors and whether the discrete solution has remained within or at least very close to the bounds given by the local minimum and maximum computed in the first step. Elements that do not satisfy all the previously mentioned detection criteria are flagged as troubled cells. For these troubled cells, the candidate solution is discarded as inappropriate and consequently needs to be recomputed. Within these troubled cells the old discrete solution at the previous time tn is scattered onto small sub-cells (Ns = 2 N + 1 sub-cells per element edge), in order to obtain a set of sub-cell averages at time tn. Then, a more robust second order TVD finite volume scheme is applied to update the sub-cell averages within the troubled DG cells from time tn to time t n + 1. The new sub-grid data at time t n + 1 are finally gathered back into a valid cell-centered DG polynomial of degree N by using a classical conservative and higher order
Augustin, C. M.; Swart, P. K.; Broad, K.
2014-12-01
Geologic carbon capture and storage (CCS) is a feasible solution to the international greenhouse gas (GHG) emissions problem and it has recently been called a "vital" mitigation tool by the International Energy Agency. However, there exists uncertainty concerning the terminal fate of stored carbon dioxide (CO2.) In this regard, reliable monitoring, verification and accounting (MVA) technologies are essential for making CCS publicly acceptable. Chiefly, MVA addresses safety and environmental concerns by providing a warning system to prevent or alleviate CO2 leakages. A secondary purpose of MVA technologies is to prove compliance with CO2 reduction standards through inventory verification. A key MVA tool for tracking CO2 leakages is surface (atmospheric) monitoring. Demonstrating its value, industry actors feel an impetus to invest in surface monitoring as a low-risk, high-value technology to mitigate liability in cases of potential leakages. Despite how necessary this tool is, to date, all surface monitoring mesh designs and best practices have been proposed locally, without discussion of standardization or optimization on a regional, national or international level. We identify the fundamental problem of surface monitoring mesh design as locating the monitoring sites to record CO2 levels over the designated geographic area at lowest cost with maximum impact. We approach this problem from both an operations research (OR) perspective and atmospheric dispersion perspective. From an OR perspective, we approach mesh design using multiobjective optimization models - we specify the relative placement of candidate sites, observation time interval, and optimality criteria. In the second approach, we model CO2 leakage scenarios to test the effectiveness of proposed mesh design from the first approach. We use atmospheric dispersion modeling softwares AERMOD and SCREEN3 - both tools developed by the United States Environmental Protection Agency and codified into law - for
Scandalously Parallelizable Mesh Generation
Bortz, David
2011-01-01
We propose a novel approach which employs random sampling to generate an accurate non-uniform mesh for numerically solving Partial Differential Equation Boundary Value Problems (PDE-BVP's). From a uniform probability distribution U over a 1D domain, we sample M discretizations of size N where M>>N. The statistical moments of the solutions to a given BVP on each of the M ultra-sparse meshes provide insight into identifying highly accurate non-uniform meshes. Essentially, we use the pointwise mean and variance of the coarse-grid solutions to construct a mapping Q(x) from uniformly to non-uniformly spaced mesh-points. The error convergence properties of the approximate solution to the PDE-BVP on the non-uniform mesh are superior to a uniform mesh for a certain class of BVP's. In particular, the method works well for BVP's with locally non-smooth solutions. We present a framework for studying the sampled sparse-mesh solutions and provide numerical evidence for the utility of this approach as applied to a set of e...
On Optimal Bilinear Quadrilateral Meshes
D' Azevedo, E.
1998-10-26
The novelty of this work is in presenting interesting error properties of two types of asymptotically optimal quadrilateral meshes for bilinear approximation. The first type of mesh has an error equidistributing property where the maximum interpolation error is asymptotically the same over all elements. The second type has faster than expected super-convergence property for certain saddle-shaped data functions. The super-convergent mesh may be an order of magnitude more accurate than the error equidistributing mesh. Both types of mesh are generated by a coordinate transformation of a regular mesh of squares. The coordinate transformation is derived by interpreting the Hessian matrix of a data function as a metric tensor. The insights in this work may have application in mesh design near known corner or point singularities.
On Optimal Bilinear Quadrilateral Meshes
D' Azevedo, E
2000-03-17
The novelty of this work is in presenting interesting error properties of two types of asymptotically ''optimal'' quadrilateral meshes for bilinear approximation. The first type of mesh has an error equidistributing property where the maximum interpolation error is asymptotically the same over all elements. The second type has faster than expected ''super-convergence'' property for certain saddle-shaped data functions. The ''superconvergent'' mesh may be an order of magnitude more accurate than the error equidistributing mesh. Both types of mesh are generated by a coordinate transformation of a regular mesh of squares. The coordinate transformation is derived by interpreting the Hessian matrix of a data function as a metric tensor. The insights in this work may have application in mesh design near corner or point singularities.
Finite element (FE) methods are widely used in electrical impedance tomography (EIT) to enable rapid image reconstruction of different tissues based on their electrical conductivity. For EIT of brain function, anatomically-accurate (head-shaped) FE meshes have been shown to improve the quality of the reconstructed images. Unfortunately, given the lack of a computational protocol to generate patient-specific meshes suitable for EIT, production of such meshes is currently ad hoc and therefore very time consuming. Here we describe a robust protocol for rapid generation of patient-specific FE meshes from MRI or CT scan data. Most of the mesh generation process is automated and uses freely available user-friendly software. Other necessary custom scripts are provided as supplementary online data and are fully documented. The patient scan data is segmented into four surfaces: brain, cerebrospinal fluid, skull and scalp. The segmented surfaces are then triangulated and used to generate a global mesh of tetrahedral elements. The resulting meshes exhibit high quality when tested with different criteria and were validated in computational simulations. The proposed protocol provides a rapid and practicable method for generation of patient-specific FE meshes of the human head that are suitable for EIT. This method could eventually be extended to other body regions and might confer benefits with other imaging techniques such as optical tomography or EEG inverse source imaging. (paper)
Parallel Adaptive Mesh Refinement
Diachin, L; Hornung, R; Plassmann, P; WIssink, A
2005-03-04
As large-scale, parallel computers have become more widely available and numerical models and algorithms have advanced, the range of physical phenomena that can be simulated has expanded dramatically. Many important science and engineering problems exhibit solutions with localized behavior where highly-detailed salient features or large gradients appear in certain regions which are separated by much larger regions where the solution is smooth. Examples include chemically-reacting flows with radiative heat transfer, high Reynolds number flows interacting with solid objects, and combustion problems where the flame front is essentially a two-dimensional sheet occupying a small part of a three-dimensional domain. Modeling such problems numerically requires approximating the governing partial differential equations on a discrete domain, or grid. Grid spacing is an important factor in determining the accuracy and cost of a computation. A fine grid may be needed to resolve key local features while a much coarser grid may suffice elsewhere. Employing a fine grid everywhere may be inefficient at best and, at worst, may make an adequately resolved simulation impractical. Moreover, the location and resolution of fine grid required for an accurate solution is a dynamic property of a problem's transient features and may not be known a priori. Adaptive mesh refinement (AMR) is a technique that can be used with both structured and unstructured meshes to adjust local grid spacing dynamically to capture solution features with an appropriate degree of resolution. Thus, computational resources can be focused where and when they are needed most to efficiently achieve an accurate solution without incurring the cost of a globally-fine grid. Figure 1.1 shows two example computations using AMR; on the left is a structured mesh calculation of a impulsively-sheared contact surface and on the right is the fuselage and volume discretization of an RAH-66 Comanche helicopter [35]. Note the
Spherical geodesic mesh generation
Fung, Jimmy [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Kenamond, Mark Andrew [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Burton, Donald E. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Shashkov, Mikhail Jurievich [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2015-02-27
In ALE simulations with moving meshes, mesh topology has a direct influence on feature representation and code robustness. In three-dimensional simulations, modeling spherical volumes and features is particularly challenging for a hydrodynamics code. Calculations on traditional spherical meshes (such as spin meshes) often lead to errors and symmetry breaking. Although the underlying differencing scheme may be modified to rectify this, the differencing scheme may not be accessible. This work documents the use of spherical geodesic meshes to mitigate solution-mesh coupling. These meshes are generated notionally by connecting geodesic surface meshes to produce triangular-prismatic volume meshes. This mesh topology is fundamentally different from traditional mesh topologies and displays superior qualities such as topological symmetry. This work describes the geodesic mesh topology as well as motivating demonstrations with the FLAG hydrocode.
A Lightweight Model for Gear Mesh Dynamics Incorporating Variable Mesh Stiffness
Fokas Nikolaos
2015-01-01
Full Text Available Variable stiffness of the gear tooth mesh for a pair of spur gears is computed using an accurate lightweight mathematical formulation. This is used to simulate gear dynamic behavior. Gear eigenfrequencies are calculated for the SDOF system and correlated with gear physical properties and the effect of stiffness variation during a mesh cycle is studied.
Applications of automatic mesh generation and adaptive methods in computational medicine
Schmidt, J.A.; Macleod, R.S. [Univ. of Utah, Salt Lake City, UT (United States); Johnson, C.R.; Eason, J.C. [Duke Univ., Durham, NC (United States)
1995-12-31
Important problems in Computational Medicine exist that can benefit from the implementation of adaptive mesh refinement techniques. Biological systems are so inherently complex that only efficient models running on state of the art hardware can begin to simulate reality. To tackle the complex geometries associated with medical applications we present a general purpose mesh generation scheme based upon the Delaunay tessellation algorithm and an iterative point generator. In addition, automatic, two- and three-dimensional adaptive mesh refinement methods are presented that are derived from local and global estimates of the finite element error. Mesh generation and adaptive refinement techniques are utilized to obtain accurate approximations of bioelectric fields within anatomically correct models of the heart and human thorax. Specifically, we explore the simulation of cardiac defibrillation and the general forward and inverse problems in electrocardiography (ECG). Comparisons between uniform and adaptive refinement techniques are made to highlight the computational efficiency and accuracy of adaptive methods in the solution of field problems in computational medicine.
Multiple Staggered Mesh Ewald: Boosting the Accuracy of the Smooth Particle Mesh Ewald Method
Wang, Han; Fang, Jun
2016-01-01
The smooth particle mesh Ewald (SPME) method is the standard method for computing the electrostatic interactions in the molecular simulations. In this work, the multiple staggered mesh Ewald (MSME) method is proposed to boost the accuracy of the SPME method. Unlike the SPME that achieves higher accuracy by refining the mesh, the MSME improves the accuracy by averaging the standard SPME forces computed on, e.g. $M$, staggered meshes. We prove, from theoretical perspective, that the MSME is as accurate as the SPME, but uses $M^2$ times less mesh points in a certain parameter range. In the complementary parameter range, the MSME is as accurate as the SPME with twice of the interpolation order. The theoretical conclusions are numerically validated both by a uniform and uncorrelated charge system, and by a three-point-charge water system that is widely used as solvent for the bio-macromolecules.
O'Hara, Ryan P.; Chand, Arpita; Vidiyala, Sowmya; Arechavala, Stacie M.; Mitsouras, Dimitrios; Rudin, Stephen; Ionita, Ciprian N.
2016-03-01
Complex vascular anatomies can cause the failure of image-guided endovascular procedures. 3D printed patient-specific vascular phantoms provide clinicians and medical device companies the ability to preemptively plan surgical treatments, test the likelihood of device success, and determine potential operative setbacks. This research aims to present advanced mesh manipulation techniques of stereolithographic (STL) files segmented from medical imaging and post-print surface optimization to match physiological vascular flow resistance. For phantom design, we developed three mesh manipulation techniques. The first method allows outlet 3D mesh manipulations to merge superfluous vessels into a single junction, decreasing the number of flow outlets and making it feasible to include smaller vessels. Next we introduced Boolean operations to eliminate the need to manually merge mesh layers and eliminate errors of mesh self-intersections that previously occurred. Finally we optimize support addition to preserve the patient anatomical geometry. For post-print surface optimization, we investigated various solutions and methods to remove support material and smooth the inner vessel surface. Solutions of chloroform, alcohol and sodium hydroxide were used to process various phantoms and hydraulic resistance was measured and compared with values reported in literature. The newly mesh manipulation methods decrease the phantom design time by 30 - 80% and allow for rapid development of accurate vascular models. We have created 3D printed vascular models with vessel diameters less than 0.5 mm. The methods presented in this work could lead to shorter design time for patient specific phantoms and better physiological simulations.
Liu, Rong
2009-01-01
Polygonal meshes are ubiquitous in geometric modeling. They are widely used in many applications, such as computer games, computer-aided design, animation, and visualization. One of the important problems in mesh processing and analysis is segmentation, where the goal is to partition a mesh into segments to suit the particular application at hand. In this thesis we study structural-level mesh segmentation, which seeks to decompose a given 3D shape into parts according to human intuition. We t...
SUPERFISH accuracy dependence on mesh size
Merson, J. L.; Boicourt, G. P.
1989-02-01
The RF cavity code SUPERFISH is extensively used for the design of drift-tube linac (DTL), radio-frequency quadrupole (RFQ), and coupled-cavity linac (CCL) structures. It has been known for some time that considerably finer meshes are required near the nose of a drift tube to ensure accurate calculation of the resonant frequency and related secondary quantities. This paper discusses the results of numerical experiments designed to provide rules to set proper mesh sizes for DTL, RFQ, and CCL problems. During this work, SUPERFISH problems involving more than 100,000 mesh points were solved.
Anatomic Total Shoulder System
Full Text Available GLOBAL AP ANATOMIC TOTAL SHOULDER SYSTEM METHODIST HOSPITAL PHILADELPHIA, PA April 17, 2008 00:00:10 ANNOUNCER: ... you'll be able to watch a live global AP anatomic total shoulder surgery from Methodist Hospital ...
Namiot, Dmitry
2015-01-01
With the advances in mobile computing technologies and the growth of the Net, mobile mesh networks are going through a set of important evolutionary steps. In this paper, we survey architectural aspects of mobile mesh networks and their use cases and deployment models. Also, we survey challenging areas of mobile mesh networks and describe our vision of promising mobile services. This paper presents a basic introductory material for Masters of Open Information Technologies Lab, interested in m...
Grid adaptation using chimera composite overlapping meshes
Kao, Kai-Hsiung; Liou, Meng-Sing; Chow, Chuen-Yen
1994-01-01
The objective of this paper is to perform grid adaptation using composite overlapping meshes in regions of large gradient to accurately capture the salient features during computation. The chimera grid scheme, a multiple overset mesh technique, is used in combination with a Navier-Stokes solver. The numerical solution is first converged to a steady state based on an initial coarse mesh. Solution-adaptive enhancement is then performed by using a secondary fine grid system which oversets on top of the base grid in the high-gradient region, but without requiring the mesh boundaries to join in any special way. Communications through boundary interfaces between those separated grids are carried out using trilinear interpolation. Application to the Euler equations for shock reflections and to shock wave/boundary layer interaction problem are tested. With the present method, the salient features are well-resolved.
Grid adaption using Chimera composite overlapping meshes
Kao, Kai-Hsiung; Liou, Meng-Sing; Chow, Chuen-Yen
1993-01-01
The objective of this paper is to perform grid adaptation using composite over-lapping meshes in regions of large gradient to capture the salient features accurately during computation. The Chimera grid scheme, a multiple overset mesh technique, is used in combination with a Navier-Stokes solver. The numerical solution is first converged to a steady state based on an initial coarse mesh. Solution-adaptive enhancement is then performed by using a secondary fine grid system which oversets on top of the base grid in the high-gradient region, but without requiring the mesh boundaries to join in any special way. Communications through boundary interfaces between those separated grids are carried out using tri-linear interpolation. Applications to the Euler equations for shock reflections and to a shock wave/boundary layer interaction problem are tested. With the present method, the salient features are well resolved.
Spacetime Meshing for Discontinuous Galerkin Methods
Thite, Shripad Vidyadhar
2008-01-01
Spacetime discontinuous Galerkin (SDG) finite element methods are used to solve such PDEs involving space and time variables arising from wave propagation phenomena in important applications in science and engineering. To support an accurate and efficient solution procedure using SDG methods and to exploit the flexibility of these methods, we give a meshing algorithm to construct an unstructured simplicial spacetime mesh over an arbitrary simplicial space domain. Our algorithm is the first spacetime meshing algorithm suitable for efficient solution of nonlinear phenomena in anisotropic media using novel discontinuous Galerkin finite element methods for implicit solutions directly in spacetime. Given a triangulated d-dimensional Euclidean space domain M (a simplicial complex) and initial conditions of the underlying hyperbolic spacetime PDE, we construct an unstructured simplicial mesh of the (d+1)-dimensional spacetime domain M x [0,infinity). Our algorithm uses a near-optimal number of spacetime elements, ea...
Pertel, Michael J.
1992-01-01
A table of useful summation formulae are derived, together with a Mathematica package for producing them. The distance distribution in mesh routing networks is derived. The mean and variance of the distance distribution are computed. A program for computing the distance distribution of any mesh is presented.
An Adaptive Mesh Algorithm: Mapping the Mesh Variables
Scannapieco, Anthony J. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2016-07-25
Both thermodynamic and kinematic variables must be mapped. The kinematic variables are defined on a separate kinematic mesh; it is the duel mesh to the thermodynamic mesh. The map of the kinematic variables is done by calculating the contributions of kinematic variables on the old thermodynamic mesh, mapping the kinematic variable contributions onto the new thermodynamic mesh and then synthesizing the mapped kinematic variables on the new kinematic mesh. In this document the map of the thermodynamic variables will be described.
Urogynecologic Surgical Mesh Implants
... be used for urogynecologic procedures, including repair of pelvic organ prolapse (POP) and stress urinary incontinence (SUI). It is ... associated with surgical mesh for transvaginal repair of pelvic organ prolapse 513(e) Proposed Order for Reclassification of Surgical ...
Geometrically Consistent Mesh Modification
Bonito, A.
2010-01-01
A new paradigm of adaptivity is to execute refinement, coarsening, and smoothing of meshes on manifolds with incomplete information about their geometry and yet preserve position and curvature accuracy. We refer to this collectively as geometrically consistent (GC) mesh modification. We discuss the concept of discrete GC, show the failure of naive approaches, and propose and analyze a simple algorithm that is GC and accuracy preserving. © 2010 Society for Industrial and Applied Mathematics.
Botsch, Mario; Pauly, Mark; Alliez, Pierre; Levy, Bruno
2010-01-01
Geometry processing, or mesh processing, is a fast-growing area of research that uses concepts from applied mathematics, computer science, and engineering to design efficient algorithms for the acquisition, reconstruction, analysis, manipulation, simulation, and transmission of complex 3D models. Applications of geometry processing algorithms already cover a wide range of areas from multimedia, entertainment, and classical computer-aided design, to biomedical computing, reverse engineering, and scientific computing. Over the last several years, triangle meshes have become increasingly popular,
To mesh or not to mesh: a review of pelvic organ reconstructive surgery
Dällenbach P
2015-04-01
Full Text Available Patrick Dällenbach Department of Gynecology and Obstetrics, Division of Gynecology, Urogynecology Unit, Geneva University Hospitals, Geneva, Switzerland Abstract: Pelvic organ prolapse (POP is a major health issue with a lifetime risk of undergoing at least one surgical intervention estimated at close to 10%. In the 1990s, the risk of reoperation after primary standard vaginal procedure was estimated to be as high as 30% to 50%. In order to reduce the risk of relapse, gynecological surgeons started to use mesh implants in pelvic organ reconstructive surgery with the emergence of new complications. Recent studies have nevertheless shown that the risk of POP recurrence requiring reoperation is lower than previously estimated, being closer to 10% rather than 30%. The development of mesh surgery – actively promoted by the marketing industry – was tremendous during the past decade, and preceded any studies supporting its benefit for our patients. Randomized trials comparing the use of mesh to native tissue repair in POP surgery have now shown better anatomical but similar functional outcomes, and meshes are associated with more complications, in particular for transvaginal mesh implants. POP is not a life-threatening condition, but a functional problem that impairs quality of life for women. The old adage “primum non nocere” is particularly appropriate when dealing with this condition which requires no treatment when asymptomatic. It is currently admitted that a certain degree of POP is physiological with aging when situated above the landmark of the hymen. Treatment should be individualized and the use of mesh needs to be selective and appropriate. Mesh implants are probably an important tool in pelvic reconstructive surgery, but the ideal implant has yet to be found. The indications for its use still require caution and discernment. This review explores the reasons behind the introduction of mesh augmentation in POP surgery, and aims to
GRChombo: Numerical relativity with adaptive mesh refinement
Clough, Katy; Figueras, Pau; Finkel, Hal; Kunesch, Markus; Lim, Eugene A.; Tunyasuvunakool, Saran
2015-12-01
In this work, we introduce {\\mathtt{GRChombo}}: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation. The code supports non-trivial 'many-boxes-in-many-boxes' mesh hierarchies and massive parallelism through the message passing interface. {\\mathtt{GRChombo}} evolves the Einstein equation using the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. The AMR capability permits the study of a range of new physics which has previously been computationally infeasible in a full 3 + 1 setting, while also significantly simplifying the process of setting up the mesh for these problems. We show that {\\mathtt{GRChombo}} can stably and accurately evolve standard spacetimes such as binary black hole mergers and scalar collapses into black holes, demonstrate the performance characteristics of our code, and discuss various physics problems which stand to benefit from the AMR technique.
3D harmonic phase tracking with anatomical regularization.
Zhou, Yitian; Bernard, Olivier; Saloux, Eric; Manrique, Alain; Allain, Pascal; Makram-Ebeid, Sherif; De Craene, Mathieu
2015-12-01
This paper presents a novel algorithm that extends HARP to handle 3D tagged MRI images. HARP results were regularized by an original regularization framework defined in an anatomical space of coordinates. In the meantime, myocardium incompressibility was integrated in order to correct the radial strain which is reported to be more challenging to recover. Both the tracking and regularization of LV displacements were done on a volumetric mesh to be computationally efficient. Also, a window-weighted regression method was extended to cardiac motion tracking which helps maintain a low complexity even at finer scales. On healthy volunteers, the tracking accuracy was found to be as accurate as the best candidates of a recent benchmark. Strain accuracy was evaluated on synthetic data, showing low bias and strain errors under 5% (excluding outliers) for longitudinal and circumferential strains, while the second and third quartiles of the radial strain errors are in the (-5%,5%) range. In clinical data, strain dispersion was shown to correlate with the extent of transmural fibrosis. Also, reduced deformation values were found inside infarcted segments. PMID:26363844
Algebraic mesh quality metrics
KNUPP,PATRICK
2000-04-24
Quality metrics for structured and unstructured mesh generation are placed within an algebraic framework to form a mathematical theory of mesh quality metrics. The theory, based on the Jacobian and related matrices, provides a means of constructing, classifying, and evaluating mesh quality metrics. The Jacobian matrix is factored into geometrically meaningful parts. A nodally-invariant Jacobian matrix can be defined for simplicial elements using a weight matrix derived from the Jacobian matrix of an ideal reference element. Scale and orientation-invariant algebraic mesh quality metrics are defined. the singular value decomposition is used to study relationships between metrics. Equivalence of the element condition number and mean ratio metrics is proved. Condition number is shown to measure the distance of an element to the set of degenerate elements. Algebraic measures for skew, length ratio, shape, volume, and orientation are defined abstractly, with specific examples given. Combined metrics for shape and volume, shape-volume-orientation are algebraically defined and examples of such metrics are given. Algebraic mesh quality metrics are extended to non-simplical elements. A series of numerical tests verify the theoretical properties of the metrics defined.
Numerical modeling of seismic waves using frequency-adaptive meshes
Hu, Jinyin; Jia, Xiaofeng
2016-08-01
An improved modeling algorithm using frequency-adaptive meshes is applied to meet the computational requirements of all seismic frequency components. It automatically adopts coarse meshes for low-frequency computations and fine meshes for high-frequency computations. The grid intervals are adaptively calculated based on a smooth inversely proportional function of grid size with respect to the frequency. In regular grid-based methods, the uniform mesh or non-uniform mesh is used for frequency-domain wave propagators and it is fixed for all frequencies. A too coarse mesh results in inaccurate high-frequency wavefields and unacceptable numerical dispersion; on the other hand, an overly fine mesh may cause storage and computational overburdens as well as invalid propagation angles of low-frequency wavefields. Experiments on the Padé generalized screen propagator indicate that the Adaptive mesh effectively solves these drawbacks of regular fixed-mesh methods, thus accurately computing the wavefield and its propagation angle in a wide frequency band. Several synthetic examples also demonstrate its feasibility for seismic modeling and migration.
Anatomic Total Shoulder System
Full Text Available GLOBAL AP ANATOMIC TOTAL SHOULDER SYSTEM METHODIST HOSPITAL PHILADELPHIA, PA April 17, 2008 00:00:10 ANNOUNCER: DePuy Orthopedics is continually advancing the standard of orthopedic patient care. In a few ...
Anatomic Total Shoulder System
Full Text Available ... Orthopedics is continually advancing the standard of orthopedic patient care. In a few moments, you'll be ... and version variability which allows adaptability to a patient's unique anatomical makeup. Dr. Gerald R. Williams, Jr., ...
Anatomic Total Shoulder System
Full Text Available ... to a patient's unique anatomical makeup. Dr. Gerald R. Williams, Jr., a shoulder specialist from the Rothman ... That might help. Could you raise the O.R. table, please? 00:28:35 WOMAN: Can you ...
Anatomic Total Shoulder System
Full Text Available GLOBAL AP ANATOMIC TOTAL SHOULDER SYSTEM METHODIST HOSPITAL PHILADELPHIA, PA April 17, 2008 00:00:10 ANNOUNCER: DePuy Orthopedics is continually advancing the standard of orthopedic patient ...
Documentation for MeshKit - Reactor Geometry (&mesh) Generator
Jain, Rajeev [Argonne National Lab. (ANL), Argonne, IL (United States); Mahadevan, Vijay [Argonne National Lab. (ANL), Argonne, IL (United States)
2015-09-30
This report gives documentation for using MeshKit’s Reactor Geometry (and mesh) Generator (RGG) GUI and also briefly documents other algorithms and tools available in MeshKit. RGG is a program designed to aid in modeling and meshing of complex/large hexagonal and rectilinear reactor cores. RGG uses Argonne’s SIGMA interfaces, Qt and VTK to produce an intuitive user interface. By integrating a 3D view of the reactor with the meshing tools and combining them into one user interface, RGG streamlines the task of preparing a simulation mesh and enables real-time feedback that reduces accidental scripting mistakes that could waste hours of meshing. RGG interfaces with MeshKit tools to consolidate the meshing process, meaning that going from model to mesh is as easy as a button click. This report is designed to explain RGG v 2.0 interface and provide users with the knowledge and skills to pilot RGG successfully. Brief documentation of MeshKit source code, tools and other algorithms available are also presented for developers to extend and add new algorithms to MeshKit. RGG tools work in serial and parallel and have been used to model complex reactor core models consisting of conical pins, load pads, several thousands of axially varying material properties of instrumentation pins and other interstices meshes.
Mizoguchi, Tomohiro; Kanai, Satoshi
Along with the rapid growth of industrial X-ray CT scanning systems, it is now possible to non-destructively acquire the entire meshes of assemblies consisting of a set of parts. For the advanced inspections of the assemblies, such as estimation of their assembling errors or examinations of their behaviors in the motions, based on their CT scanned meshes, it is necessary to accurately decompose the mesh and to extract a set of partial meshes each of which correspond to a part. Moreover it is required to create models which can be used for the real-product based simulations. In this paper, we focus on CT scanned meshes of gear assemblies as examples and propose beneficial methods for establishing such advance inspections of the assemblies. We first propose a method that accurately decomposes the mesh into partial meshes each of which corresponds to a gear based on periodicity recognitions. The key idea is first to accurately recognize the periodicity of each gear and then to extract the partial meshes as sets of topologically connected mesh elements where periodicities are valid. Our method can robustly and accurately recognize periodicities from noisy scanned meshes. In contrast to previous methods, our method can deal with single-material CT scanned meshes and can estimate the correct boundaries of neighboring parts with no previous knowledge. Moreover it can efficiently extract the partial meshes from large scanned meshes containing about one million triangles in a few minutes. We also propose a method for creating simulation models which can be used for a gear teeth contact evaluation using extracted partial meshes and their periodicities. Such an evaluation of teeth contacts is one of the most important functions in kinematic simulations of gear assemblies for predicting the power transmission efficiency, noise and vibration. We demonstrate the effectiveness of our method on a variety of artificial and CT scanned meshes.
Badal, Andreu; Kyprianou, Iacovos; Badano, Aldo; Sempau, Josep; Myers, Kyle J.
2007-03-01
X-ray imaging system optimization increases the benefit-to-cost ratio by reducing the radiation dose to the patient while maximizing image quality. We present a new simulation tool for the generation of realistic medical x-ray images for assessment and optimization of complete imaging systems. The Monte Carlo code simulates radiation transport physics using the subroutine package PENELOPE, which accurately simulates the transport of electrons and photons within the typical medical imaging energy range. The new code implements a novel object-oriented geometry package that allows simulations with homogeneous objects of arbitrary shapes described by triangle meshes. The flexibility of this code, which uses the industry standard PLY input-file format, allows the use of detailed anatomical models developed using computer-aided design tools applied to segmented CT and MRI data. The use of triangle meshes highly simplifies the ray-tracing algorithm without reducing the generality of the code, since most surface models can be tessellated into triangles while retaining their geometric details. Our algorithm incorporates an octree spatial data structure to sort the triangles and accelerate the simulation, reaching execution speeds comparable to the original quadric geometry model of PENELOPE. Coronary angiograms were simulated using a tessellated version of the NURBS-based Cardiac-Torso (NCAT) phantom. The phantom models 330 objects, comprised in total of 5 million triangles. The dose received by each organ and the contribution of the different scattering processes to the final image were studied in detail.
Difference schemes on non-uniform mesh and their application
MA Yanwen; GAO Hui; FU Dexun; LI Xinliang
2004-01-01
High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.
Personalized identification of abdominal wall hernia meshes on computed tomography.
Pham, Tuan D; Le, Dinh T P; Xu, Jinwei; Nguyen, Duc T; Martindale, Robert G; Deveney, Clifford W
2014-01-01
An abdominal wall hernia is a protrusion of the intestine through an opening or area of weakness in the abdominal wall. Correct pre-operative identification of abdominal wall hernia meshes could help surgeons adjust the surgical plan to meet the expected difficulty and morbidity of operating through or removing the previous mesh. First, we present herein for the first time the application of image analysis for automated identification of hernia meshes. Second, we discuss the novel development of a new entropy-based image texture feature using geostatistics and indicator kriging. Third, we seek to enhance the hernia mesh identification by combining the new texture feature with the gray-level co-occurrence matrix feature of the image. The two features can characterize complementary information of anatomic details of the abdominal hernia wall and its mesh on computed tomography. Experimental results have demonstrated the effectiveness of the proposed study. The new computational tool has potential for personalized mesh identification which can assist surgeons in the diagnosis and repair of complex abdominal wall hernias. PMID:24184112
Volume-conserving mesh smoothing for front-tracking methods
Among the various direct numerical simulation (DNS) methods dedicated to multiphase flow, the front-tracking methods that use a Lagrangian mesh to describe explicitly the interfaces are generally considered as a very accurate and complex method. In this family of methods, while a fine Lagrangian mesh is desirable for a better representation of the interfacial area, the surface forces and the bubble or droplet volume, one cannot arbitrarily choose the Lagrangian mesh size. Indeed, the Lagrangian mesh displacement algorithm is unstable if the number of Lagrangian degrees of freedom does not match the number of involved Eulerian velocity points. As a consequence, in traditional front-tracking implementations, an accurate description of the interfaces is expensive in terms of Eulerian mesh cells. We demonstrate that a front-tracking interface smoothing (FTIS) method can reduce the constraints on the mesh sizes. It consists in damping the highest spatial frequency components of the Lagrangian mesh to compensate for the lack of Eulerian velocity points. The test case of fundamental proper frequency of a bubble proves the validity of the FTIS method. An example of a 3D-bubble rising shows the interest and the potential applications of the FTIS method. (authors)
Seismic Wave Simulation for Complex Rheologies on Unstructured Meshes
de la Puente, Josep
2008-01-01
The possibility of using accurate numerical methods to simulate seismic wavefields on unstructured meshes for complex rheologies is explored. In particular, the Discontinuous Galerkin (DG) finite element method for seismic wave propagation is extended to the rheological types of viscoelasticity, anisotropy and poroelasticity. First is presented the DG method for the elastic isotropic case on tetrahedral unstructured meshes. Then an extension to viscoelastic wave propagation based upon a Gener...
Mesh generation technology for nuclear reactor simulation; barriers and opportunities
Mesh generation in support of nuclear reactor simulation has much in common with the requirements of other application areas, such as computational fluid dynamics (CFD). Indeed, fluid dynamics analysis of the coolant behavior inside the reactor core is an internal flow problem that requires the resolution of spatial and temporal variations in the flow caused by complex component configurations, fluids/structure interaction, turbulence, and thermal heating of the coolant. Typical concerns of meshing complex geometries; the use of hexahedral vs. tetrahedral elements, element geometric quality, mesh smoothness, use of anisotropic elements in the thermal boundary layer, etc., are all considerations important to the reactor meshing problem. Reactor meshing begins to become more specialized as the need to employ reactor simulation as a predictive design and safety analysis capability grows in importance. First, a predictive capability will require more precise physical models to be included, and these models will need to be supported by a computational science framework that will allow them to be accurately approximated both spatially and temporally during the reactor core analysis. Both the multiphysical nature of the composite reactor model and details of the physics algorithms themselves will place new requirements on the meshing process needed to support multidimensional reactor simulation. This article discusses the current state of meshing technology applied to reactor simulation and examines a set of issues that are important in the generation of high-quality reactor meshes today and in the future
Wang, Xinheng
2008-01-01
Wireless telemedicine using GSM and GPRS technologies can only provide low bandwidth connections, which makes it difficult to transmit images and video. Satellite or 3G wireless transmission provides greater bandwidth, but the running costs are high. Wireless networks (WLANs) appear promising, since they can supply high bandwidth at low cost. However, the WLAN technology has limitations, such as coverage. A new wireless networking technology named the wireless mesh network (WMN) overcomes some of the limitations of the WLAN. A WMN combines the characteristics of both a WLAN and ad hoc networks, thus forming an intelligent, large scale and broadband wireless network. These features are attractive for telemedicine and telecare because of the ability to provide data, voice and video communications over a large area. One successful wireless telemedicine project which uses wireless mesh technology is the Emergency Room Link (ER-LINK) in Tucson, Arizona, USA. There are three key characteristics of a WMN: self-organization, including self-management and self-healing; dynamic changes in network topology; and scalability. What we may now see is a shift from mobile communication and satellite systems for wireless telemedicine to the use of wireless networks based on mesh technology, since the latter are very attractive in terms of cost, reliability and speed. PMID:19047448
Anatomic Total Shoulder System
Full Text Available ... AP ANATOMIC TOTAL SHOULDER SYSTEM METHODIST HOSPITAL PHILADELPHIA, PA April 17, 2008 00:00:10 ANNOUNCER: DePuy ... you don't make a bunch of small passes at the lesser tuberosity and make it a ...
Anatomic Total Shoulder System
Full Text Available ... Anatomic Total Shoulder surgery, which featured the latest innovation in shoulder surgery from DePuy Orthopedics. OR-Live makes it easy for you to learn more. Just click on the "Request Information" button on your webcast screen and open the door to informed medical care. 01:21: ...
Anatomic Total Shoulder System
Full Text Available ... advancing the standard of orthopedic patient care. In a few moments, you'll be able to watch a live global AP anatomic total shoulder surgery from Methodist Hospital in Philadelphia. A revolution in shoulder orthopedics, the Global AP gives ...
Anatomic Total Shoulder System
Full Text Available ... by almost ten years, is shoulders. So by definition, the average shoulder-replacement patient is almost ten ... Anatomic Total Shoulder surgery, which featured the latest innovation in shoulder surgery from DePuy Orthopedics. OR-Live ...
Identification of anatomical terminology in medical text.
Sneiderman, C. A.; Rindflesch, T. C.; Bean, C. A.
1998-01-01
We report on an experiment to use the natural language processing tools being developed in the SPECIALIST system to accurately identify terminology associated with the coronary arteries as expressed in coronary catheterization reports. The ultimate goal is to map from any anatomically-oriented medical text to online images, using the UMLS as an intermediate knowledge source. We describe some of the problems encountered when processing coronary artery terminology and report on the results of a...
Anatomically Plausible Surface Alignment and Reconstruction
Paulsen, Rasmus R.; Larsen, Rasmus
2010-01-01
With the increasing clinical use of 3D surface scanners, there is a need for accurate and reliable algorithms that can produce anatomically plausible surfaces. In this paper, a combined method for surface alignment and reconstruction is proposed. It is based on an implicit surface representation...... energy that has earlier proved to be particularly well suited for human surface scans. The method has been tested on full cranial scans of ten test subjects and on several scans of the outer human ear....
Notes on the Mesh Handler and Mesh Data Conversion
At the outset of the development of the thermal-hydraulic code (THC), efforts have been made to utilize the recent technology of the computational fluid dynamics. Among many of them, the unstructured mesh approach was adopted to alleviate the restriction of the grid handling system. As a natural consequence, a mesh handler (MH) has been developed to manipulate the complex mesh data from the mesh generator. The mesh generator, Gambit, was chosen at the beginning of the development of the code. But a new mesh generator, Pointwise, was introduced to get more flexible mesh generation capability. An open source code, Paraview, was chosen as a post processor, which can handle unstructured as well as structured mesh data. Overall data processing system for THC is shown in Figure-1. There are various file formats to save the mesh data in the permanent storage media. A couple of dozen of file formats are found even in the above mentioned programs. A competent mesh handler should have the capability to import or export mesh data as many as possible formats. But, in reality, there are two aspects that make it difficult to achieve the competence. The first aspect to consider is the time and efforts to program the interface code. And the second aspect, which is even more difficult one, is the fact that many mesh data file formats are proprietary information. In this paper, some experience of the development of the format conversion programs will be presented. File formats involved are Gambit neutral format, Ansys-CFX grid file format, VTK legacy file format, Nastran format and CGNS
Serial and parallel dynamic adaptation of general hybrid meshes
Kavouklis, Christos
The Navier-Stokes equations are a standard mathematical representation of viscous fluid flow. Their numerical solution in three dimensions remains a computationally intensive and challenging task, despite recent advances in computer speed and memory. A strategy to increase accuracy of Navier-Stokes simulations, while maintaining computing resources to a minimum, is local refinement of the associated computational mesh in regions of large solution gradients and coarsening in regions where the solution does not vary appreciably. In this work we consider adaptation of general hybrid meshes for Computational Fluid Dynamics (CFD) applications. Hybrid meshes are composed of four types of elements; hexahedra, prisms, pyramids and tetrahedra, and have been proven a promising technology in accurately resolving fluid flow for complex geometries. The first part of this dissertation is concerned with the design and implementation of a serial scheme for the adaptation of general three dimensional hybrid meshes. We have defined 29 refinement types, for all four kinds of elements. The core of the present adaptation scheme is an iterative algorithm that flags mesh edges for refinement, so that the adapted mesh is conformal. Of primary importance is considered the design of a suitable dynamic data structure that facilitates refinement and coarsening operations and furthermore minimizes memory requirements. A special dynamic list is defined for mesh elements, in contrast with the usual tree structures. It contains only elements of the current adaptation step and minimal information that is utilized to reconstruct parent elements when the mesh is coarsened. In the second part of this work, a new parallel dynamic mesh adaptation and load balancing algorithm for general hybrid meshes is presented. Partitioning of a hybrid mesh reduces to partitioning of the corresponding dual graph. Communication among processors is based on the faces of the interpartition boundary. The distributed
Early fetal anatomical sonography.
Donnelly, Jennifer C
2012-10-01
Over the past decade, prenatal screening and diagnosis has moved from the second into the first trimester, with aneuploidy screening becoming both feasible and effective. With vast improvements in ultrasound technology, sonologists can now image the fetus in greater detail at all gestational ages. In the hands of experienced sonographers, anatomic surveys between 11 and 14 weeks can be carried out with good visualisation rates of many structures. It is important to be familiar with the normal development of the embryo and fetus, and to be aware of the major anatomical landmarks whose absence or presence may be deemed normal or abnormal depending on the gestational age. Some structural abnormalities will nearly always be detected, some will never be and some are potentially detectable depending on a number of factors.
Reference Man anatomical model
Cristy, M.
1994-10-01
The 70-kg Standard Man or Reference Man has been used in physiological models since at least the 1920s to represent adult males. It came into use in radiation protection in the late 1940s and was developed extensively during the 1950s and used by the International Commission on Radiological Protection (ICRP) in its Publication 2 in 1959. The current Reference Man for Purposes of Radiation Protection is a monumental book published in 1975 by the ICRP as ICRP Publication 23. It has a wealth of information useful for radiation dosimetry, including anatomical and physiological data, gross and elemental composition of the body and organs and tissues of the body. The anatomical data includes specified reference values for an adult male and an adult female. Other reference values are primarily for the adult male. The anatomical data include much data on fetuses and children, although reference values are not established. There is an ICRP task group currently working on revising selected parts of the Reference Man document.
Standardized anatomic space for abdominal fat quantification
Tong, Yubing; Udupa, Jayaram K.; Torigian, Drew A.
2014-03-01
The ability to accurately measure subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT) from images is important for improved assessment and management of patients with various conditions such as obesity, diabetes mellitus, obstructive sleep apnea, cardiovascular disease, kidney disease, and degenerative disease. Although imaging and analysis methods to measure the volume of these tissue components have been developed [1, 2], in clinical practice, an estimate of the amount of fat is obtained from just one transverse abdominal CT slice typically acquired at the level of the L4-L5 vertebrae for various reasons including decreased radiation exposure and cost [3-5]. It is generally assumed that such an estimate reliably depicts the burden of fat in the body. This paper sets out to answer two questions related to this issue which have not been addressed in the literature. How does one ensure that the slices used for correlation calculation from different subjects are at the same anatomic location? At what anatomic location do the volumes of SAT and VAT correlate maximally with the corresponding single-slice area measures? To answer these questions, we propose two approaches for slice localization: linear mapping and non-linear mapping which is a novel learning based strategy for mapping slice locations to a standardized anatomic space so that same anatomic slice locations are identified in different subjects. We then study the volume-to-area correlations and determine where they become maximal. We demonstrate on 50 abdominal CT data sets that this mapping achieves significantly improved consistency of anatomic localization compared to current practice. Our results also indicate that maximum correlations are achieved at different anatomic locations for SAT and VAT which are both different from the L4-L5 junction commonly utilized.
VPN Mesh in Industrial Networking
Berndtsson, Andreas
2013-01-01
This thesis report describes the process and present the results gained while evaluating available VPN mesh solutions and equipment for integration into Industrial systems. The task was divided into several sub steps; summarize the previous work done in the VPN mesh area, evaluate the available VPN mesh solutions, verify that the interesting equipment comply with the criteria set by ABB and lastly verify that the equipment can be integrated transparently into already running systems. The resu...
Papadakis, A. P.; Georghiou, G. E.; Metaxas, A. C.
2008-12-01
A new adaptive mesh generator has been developed and used in the analysis of high-pressure gas discharges, such as avalanches and streamers, reducing computational times and computer memory needs significantly. The new adaptive mesh generator developed, uses normalized error indicators, varying from 0 to 1, to guarantee optimal mesh resolution for all carriers involved in the analysis. Furthermore, it uses h- and r-refinement techniques such as mesh jiggling, edge swapping and node addition/removal to develop an element quality improvement algorithm that improves the mesh quality significantly and a fast and accurate algorithm for interpolation between meshes. Finally, the mesh generator is applied in the characterization of the transition from a single electron to the avalanche and streamer discharges in high-voltage, high-pressure gas discharges for dc 1 mm gaps, RF 1 cm point-plane gaps and parallel-plate 40 MHz configurations, in ambient atmospheric air.
Papadakis, A P [Department of Electrical Engineering, Frederick University Cyprus, 7 Y Frederickou Street, Palouriotissa, Nicosia 1036 (Cyprus); Georghiou, G E [Department of Electrical and Computer Engineering, University of Cyprus, 75 Kallipoleos, PO Box 20577, 1678, Nicosia (Cyprus); Metaxas, A C [St John' s College, University of Cambridge, Cambridge, CB2 1TP (United Kingdom)], E-mail: eng.ap@frederick.ac.cy, E-mail: geg@ucy.ac.cy, E-mail: acm33@cam.ac.uk
2008-12-07
A new adaptive mesh generator has been developed and used in the analysis of high-pressure gas discharges, such as avalanches and streamers, reducing computational times and computer memory needs significantly. The new adaptive mesh generator developed, uses normalized error indicators, varying from 0 to 1, to guarantee optimal mesh resolution for all carriers involved in the analysis. Furthermore, it uses h- and r-refinement techniques such as mesh jiggling, edge swapping and node addition/removal to develop an element quality improvement algorithm that improves the mesh quality significantly and a fast and accurate algorithm for interpolation between meshes. Finally, the mesh generator is applied in the characterization of the transition from a single electron to the avalanche and streamer discharges in high-voltage, high-pressure gas discharges for dc 1 mm gaps, RF 1 cm point-plane gaps and parallel-plate 40 MHz configurations, in ambient atmospheric air.
Solving the resonating-group equation on a Lagrange mesh
Hesse, M; Baye, D
2002-01-01
The resonating-group method allows treating reactions in a fully microscopic way. The non-local resonating-group equation can be accurately solved on a Lagrange mesh involving few mesh points. This mesh technique is combined with either the R-matrix method or the Hulthen-Kohn method. The forbidden states can be eliminated by a special treatment. The accuracy of the technique of solution is illustrated on a solvable non-local potential. Phase shifts for the alpha+n and alpha+p scatterings are calculated with both variants of the resonating-group method on a Lagrange mesh and a comparison is performed between them and the equivalent generator-coordinate method.
Planet-disc interaction on a freely moving mesh
Munoz, Diego J; Springel, Volker; Hernquist, Lars
2014-01-01
General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with ...
Mesh quality improvement for SciDAC applications
Accurate and efficient numerical solution of partial differential equations requires well-formed meshes that are non-inverted, smooth, well-shaped, oriented, and size-adapted. The Mesquite mesh quality improvement toolkit is a software library that applies optimization algorithms to create well-formed meshes via node movement. Mesquite can be run standalone using drivers or called directly from an application code. Mesquite can play an essential role in the SLAC accelerator design program as a component in automatic shape optimization software and in manufacturing defect-correction studies to smoothly deform meshes in response to geometric domain deformations guided by the optimization of design parameters. Mesquite has also been applied to problems in fusion, biology, and propellant burn studies
Synthesized Optimization of Triangular Mesh
HU Wenqiang; YANG Wenyu
2006-01-01
Triangular mesh is often used to describe geometric object as computed model in digital manufacture, thus the mesh model with both uniform triangular shape and excellent geometric shape is expected. But in fact, the optimization of triangular shape often is contrary with that of geometric shape. In this paper, one synthesized optimizing algorithm is presented through subdividing triangles to achieve the trade-off solution between the geometric and triangular shape optimization of mesh model. The result mesh with uniform triangular shape and excellent topology are obtained.
Surface meshing with curvature convergence
Li, Huibin
2014-06-01
Surface meshing plays a fundamental role in graphics and visualization. Many geometric processing tasks involve solving geometric PDEs on meshes. The numerical stability, convergence rates and approximation errors are largely determined by the mesh qualities. In practice, Delaunay refinement algorithms offer satisfactory solutions to high quality mesh generations. The theoretical proofs for volume based and surface based Delaunay refinement algorithms have been established, but those for conformal parameterization based ones remain wide open. This work focuses on the curvature measure convergence for the conformal parameterization based Delaunay refinement algorithms. Given a metric surface, the proposed approach triangulates its conformal uniformization domain by the planar Delaunay refinement algorithms, and produces a high quality mesh. We give explicit estimates for the Hausdorff distance, the normal deviation, and the differences in curvature measures between the surface and the mesh. In contrast to the conventional results based on volumetric Delaunay refinement, our stronger estimates are independent of the mesh structure and directly guarantee the convergence of curvature measures. Meanwhile, our result on Gaussian curvature measure is intrinsic to the Riemannian metric and independent of the embedding. In practice, our meshing algorithm is much easier to implement and much more efficient. The experimental results verified our theoretical results and demonstrated the efficiency of the meshing algorithm. © 2014 IEEE.
h-Adaptive Mesh Generation using Electric Field Intensity Value as a Criterion (in Japanese)
Toyonaga, Kiyomi; Cingoski, Vlatko; Kaneda, Kazufumi; Yamashita, Hideo
1994-01-01
Finite mesh divisions are essential to obtain accurate solution of two dimensional electric field analysis. It requires the technical knowledge to generate a suitable fine mesh divisions. In electric field problem, analysts are usually interested in the electric field intensity and its distribution. In order to obtain electric field intensity with high-accuracy, we have developed and adaptive mesh generator using electric field intensity value as a criterion.
Anatomical imaging for radiotherapy
Evans, Philip M [Joint Physics Department, Institute of Cancer Research and Royal Marsden NHS Foundation Trust, Downs Road, Sutton, Surrey SM2 5PT (United Kingdom)], E-mail: phil.evans@icr.ac.uk
2008-06-21
The goal of radiation therapy is to achieve maximal therapeutic benefit expressed in terms of a high probability of local control of disease with minimal side effects. Physically this often equates to the delivery of a high dose of radiation to the tumour or target region whilst maintaining an acceptably low dose to other tissues, particularly those adjacent to the target. Techniques such as intensity modulated radiotherapy (IMRT), stereotactic radiosurgery and computer planned brachytherapy provide the means to calculate the radiation dose delivery to achieve the desired dose distribution. Imaging is an essential tool in all state of the art planning and delivery techniques: (i) to enable planning of the desired treatment, (ii) to verify the treatment is delivered as planned and (iii) to follow-up treatment outcome to monitor that the treatment has had the desired effect. Clinical imaging techniques can be loosely classified into anatomic methods which measure the basic physical characteristics of tissue such as their density and biological imaging techniques which measure functional characteristics such as metabolism. In this review we consider anatomical imaging techniques. Biological imaging is considered in another article. Anatomical imaging is generally used for goals (i) and (ii) above. Computed tomography (CT) has been the mainstay of anatomical treatment planning for many years, enabling some delineation of soft tissue as well as radiation attenuation estimation for dose prediction. Magnetic resonance imaging is fast becoming widespread alongside CT, enabling superior soft-tissue visualization. Traditionally scanning for treatment planning has relied on the use of a single snapshot scan. Recent years have seen the development of techniques such as 4D CT and adaptive radiotherapy (ART). In 4D CT raw data are encoded with phase information and reconstructed to yield a set of scans detailing motion through the breathing, or cardiac, cycle. In ART a set of
Anatomical imaging for radiotherapy
The goal of radiation therapy is to achieve maximal therapeutic benefit expressed in terms of a high probability of local control of disease with minimal side effects. Physically this often equates to the delivery of a high dose of radiation to the tumour or target region whilst maintaining an acceptably low dose to other tissues, particularly those adjacent to the target. Techniques such as intensity modulated radiotherapy (IMRT), stereotactic radiosurgery and computer planned brachytherapy provide the means to calculate the radiation dose delivery to achieve the desired dose distribution. Imaging is an essential tool in all state of the art planning and delivery techniques: (i) to enable planning of the desired treatment, (ii) to verify the treatment is delivered as planned and (iii) to follow-up treatment outcome to monitor that the treatment has had the desired effect. Clinical imaging techniques can be loosely classified into anatomic methods which measure the basic physical characteristics of tissue such as their density and biological imaging techniques which measure functional characteristics such as metabolism. In this review we consider anatomical imaging techniques. Biological imaging is considered in another article. Anatomical imaging is generally used for goals (i) and (ii) above. Computed tomography (CT) has been the mainstay of anatomical treatment planning for many years, enabling some delineation of soft tissue as well as radiation attenuation estimation for dose prediction. Magnetic resonance imaging is fast becoming widespread alongside CT, enabling superior soft-tissue visualization. Traditionally scanning for treatment planning has relied on the use of a single snapshot scan. Recent years have seen the development of techniques such as 4D CT and adaptive radiotherapy (ART). In 4D CT raw data are encoded with phase information and reconstructed to yield a set of scans detailing motion through the breathing, or cardiac, cycle. In ART a set of
Ye, Dexin; Lu, Ling; Joannopoulos, John D; Soljačić, Marin; Ran, Lixin
2016-03-01
A solid material possessing identical electromagnetic properties as air has yet to be found in nature. Such a medium of arbitrary shape would neither reflect nor refract light at any angle of incidence in free space. Here, we introduce nonscattering corrugated metallic wires to construct such a medium. This was accomplished by aligning the dark-state frequencies in multiple scattering channels of a single wire. Analytical solutions, full-wave simulations, and microwave measurement results on 3D printed samples show omnidirectional invisibility in any configuration. This invisible metallic mesh can improve mechanical stability, electrical conduction, and heat dissipation of a system, without disturbing the electromagnetic design. Our approach is simple, robust, and scalable to higher frequencies. PMID:26884208
A discretization of the multigroup PN radiative transfer equation on general meshes
Hermeline, F.
2016-05-01
We propose and study a finite volume method of discrete duality type for discretizing the multigroup PN approximation of radiative transfer equation on general meshes. This method is second order-accurate on a very large variety of meshes, stable under a Courant-Friedrichs-Lewy condition and it preserves naturally the diffusion asymptotic limit.
Rich: Open Source Hydrodynamic Simulation on a Moving Voronoi Mesh
Yalinewich, Almog; Sari, Re'em
2014-01-01
We present here RICH, a state of the art 2D hydrodynamic code based on Godunov's method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time advancement scheme as well as a novel parallelization scheme based on Voronoi tessellation. Using our code we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time advancement scheme is more accurate and robust than AREPO's, when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Llyod iterations) and it effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way, and a sound wave is traveling in the other way (such that relative to the grid the...
3D Mesh Compression and Transmission for Mobile Robotic Applications
Bailin Yang
2016-01-01
Full Text Available Mobile robots are useful for environment exploration and rescue operations. In such applications, it is crucial to accurately analyse and represent an environment, providing appropriate inputs for motion planning in order to support robot navigation and operations. 2D mapping methods are simple but cannot handle multilevel or multistory environments. To address this problem, 3D mapping methods generate structural 3D representations of the robot operating environment and its objects by 3D mesh reconstruction. However, they face the challenge of efficiently transmitting those 3D representations to system modules for 3D mapping, motion planning, and robot operation visualization. This paper proposes a quality-driven mesh compression and transmission method to address this. Our method is efficient, as it compresses a mesh by quantizing its transformed vertices without the need to spend time constructing an a-priori structure over the mesh. A visual distortion function is developed to govern the level of quantization, allowing mesh transmission to be controlled under different network conditions or time constraints. Our experiments demonstrate how the visual quality of a mesh can be manipulated by the visual distortion function.
Transoral Surgery: An Anatomic Study
Rock, Jack P.; Tomecek, Frank J.; Ross, Lawrence
1993-01-01
The transoral approaches have become commonplace in modern neurosurgical practice for treatment of ventral midline lesions of the clivus and upper cervical spine. Although the standard technique of transoral surgery is conceptually simple, anatomic relationships are not so readily appreciated. The present study was undertaken in an effort to define more clearly the midline anatomic relationships as they pertain to the standard transoral and transpalatine operations. The anatomic relationships...
21 CFR 878.3300 - Surgical mesh.
2010-04-01
... GENERAL AND PLASTIC SURGERY DEVICES Prosthetic Devices § 878.3300 Surgical mesh. (a) Identification... acetabular and cement restrictor mesh used during orthopedic surgery. (b) Classification. Class II....
Mesh Adaptation and Shape Optimization on Unstructured Meshes Project
National Aeronautics and Space Administration — In this SBIR CRM proposes to implement the entropy adjoint method for solution adaptive mesh refinement into the Loci/CHEM unstructured flow solver. The scheme will...
Tetrahedral mesh for needle insertion
Syvertsen, Rolf Anders
2007-01-01
This is a Master’s thesis in how to make a tetrahedral mesh for use in a needle insertion simulator. It also describes how it is possible to make the simulator, and how to improve it to make it as realistic as possible. The medical simulator uses a haptic device, a haptic scene graph and a FEM for realistic soft tissue deformation and interaction. In this project a tetrahedral mesh is created from a polygon model, and then the mesh has been loaded into the HaptX haptic scene graph. The object...
Nanowire mesh solar fuels generator
Yang, Peidong; Chan, Candace; Sun, Jianwei; Liu, Bin
2016-05-24
This disclosure provides systems, methods, and apparatus related to a nanowire mesh solar fuels generator. In one aspect, a nanowire mesh solar fuels generator includes (1) a photoanode configured to perform water oxidation and (2) a photocathode configured to perform water reduction. The photocathode is in electrical contact with the photoanode. The photoanode may include a high surface area network of photoanode nanowires. The photocathode may include a high surface area network of photocathode nanowires. In some embodiments, the nanowire mesh solar fuels generator may include an ion conductive polymer infiltrating the photoanode and the photocathode in the region where the photocathode is in electrical contact with the photoanode.
Tangle-Free Mesh Motion for Ablation Simulations
Droba, Justin
2016-01-01
Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and
Shuo Liang
2015-01-01
Full Text Available Background: Although repair augmented with mesh has been proved its priority in anatomical and functional recovery after anterior compartment reconstruction, the data about posterior compartment are scarce. The aim of this study was to compare bowel functional outcome of posterior vaginal compartment repair with and without mesh in patients with pelvic organ prolapse (POP. Methods: This was a prospective, double-blind, clinical pilot study of 22 postmenopausal women with symptomatic POP (overall POP-quantification [POP-Q] Stage III-IV who underwent total pelvic floor reconstruction. Patients were grouped according to the use of mesh for posterior vaginal compartment repair: A mesh group and a nonmesh group. POP-Q stage, the pelvic floor impact questionnaire short form-7 (PFIQ-7 and anorectal manometry were evaluated before and 3 months after surgery. Anatomical success was defined as POP-Q Stage II or less. A t-test was used to compare preoperative with postoperative data in the two groups. Results: Totally, 17 (71% were available for the follow-up. POP-Q measurements improved significantly compared to baseline (P 0.05. Compared with baseline, the nonmesh group exhibited a statistically significant decrease in anal residual pressure, a significant increase in the anorectal pressure difference during bowel movement, and a reduced rate of dyssynergia defecation pattern (P < 0.05. Conclusions: Provided there is sufficient support for the anterior wall and apex of vagina with mesh, posterior compartment repair without mesh may be as effective as repair with mesh for anatomical recovery while providing better anorectal motor function.
Medical Image Processing for Fully Integrated Subject Specific Whole Brain Mesh Generation
Chih-Yang Hsu
2015-05-01
Full Text Available Currently, anatomically consistent segmentation of vascular trees acquired with magnetic resonance imaging requires the use of multiple image processing steps, which, in turn, depend on manual intervention. In effect, segmentation of vascular trees from medical images is time consuming and error prone due to the tortuous geometry and weak signal in small blood vessels. To overcome errors and accelerate the image processing time, we introduce an automatic image processing pipeline for constructing subject specific computational meshes for entire cerebral vasculature, including segmentation of ancillary structures; the grey and white matter, cerebrospinal fluid space, skull, and scalp. To demonstrate the validity of the new pipeline, we segmented the entire intracranial compartment with special attention of the angioarchitecture from magnetic resonance imaging acquired for two healthy volunteers. The raw images were processed through our pipeline for automatic segmentation and mesh generation. Due to partial volume effect and finite resolution, the computational meshes intersect with each other at respective interfaces. To eliminate anatomically inconsistent overlap, we utilized morphological operations to separate the structures with a physiologically sound gap spaces. The resulting meshes exhibit anatomically correct spatial extent and relative positions without intersections. For validation, we computed critical biometrics of the angioarchitecture, the cortical surfaces, ventricular system, and cerebrospinal fluid (CSF spaces and compared against literature values. Volumina and surface areas of the computational mesh were found to be in physiological ranges. In conclusion, we present an automatic image processing pipeline to automate the segmentation of the main intracranial compartments including a subject-specific vascular trees. These computational meshes can be used in 3D immersive visualization for diagnosis, surgery planning with haptics
A Study on the Efficient Mesh Generation for Finite Element Analysis of Electric Machinery
Kim, Jin Tae [Hangul and Computer Co.,Ltd. (Korea); Chung, Tae Kyung [Chungang University (Korea); Kim, Hyeong Seok [Soonchunhyang University (Korea)
1999-04-01
To obtain more accurate result in the finite element analysis for electric machinery, it is important to have a mesh being of good quality. This paper describes a new technique of mesh generation for the finite element method. When the list of points defining the region of analysis is given, an appropriate distribution of interior points is generated first. Secondly the points are connected to form the triangles. Finally the connectivity data are used to reposition the interior points using laplacian smoothing and mesh relaxation technique. In this paper, a mesh searching technique of Lawson which modifies the start mesh is proposed in addition to the above three steps. This algorithm is simple and produces the meshes being of good quality with high speed in comparison with the existing one. (author). 16 refs., 10 figs., 2 tabs.
Mersiline mesh in premaxillary augmentation.
Foda, Hossam M T
2005-01-01
Premaxillary retrusion may distort the aesthetic appearance of the columella, lip, and nasal tip. This defect is characteristically seen in, but not limited to, patients with cleft lip nasal deformity. This study investigated 60 patients presenting with premaxillary deficiencies in which Mersiline mesh was used to augment the premaxilla. All the cases had surgery using the external rhinoplasty technique. Two methods of augmentation with Mersiline mesh were used: the Mersiline roll technique, for the cases with central symmetric deficiencies, and the Mersiline packing technique, for the cases with asymmetric deficiencies. Premaxillary augmentation with Mersiline mesh proved to be simple technically, easy to perform, and not associated with any complications. Periodic follow-up evaluation for a mean period of 32 months (range, 12-98 months) showed that an adequate degree of premaxillary augmentation was maintained with no clinically detectable resorption of the mesh implant. PMID:15959688
Generic Mesh Refinement On GPU
Boubekeur, Tamy; Schlick, Christophe
2005-01-01
International audience Many recent publications have shown that a large variety of computation involved in computer graphics can be moved from the CPU to the GPU, by a clever use of vertex or fragment shaders. Nonetheless there is still one kind of algorithms that is hard to translate from CPU to GPU: mesh refinement techniques. The main reason for this, is that vertex shaders available on current graphics hardware do not allow the generation of additional vertices on a mesh stored in grap...
GENERATION OF IRREGULAR HEXAGONAL MESHES
Vlasov Aleksandr Nikolaevich
2012-07-01
Decomposition is performed in a constructive way and, as option, it involves meshless representation. Further, this mapping method is used to generate the calculation mesh. In this paper, the authors analyze different cases of mapping onto simply connected and bi-connected canonical domains. They represent forward and backward mapping techniques. Their potential application for generation of nonuniform meshes within the framework of the asymptotic homogenization theory is also performed to assess and project effective characteristics of heterogeneous materials (composites.
Image-driven mesh optimization
Lindstrom, P; Turk, G
2001-01-05
We describe a method of improving the appearance of a low vertex count mesh in a manner that is guided by rendered images of the original, detailed mesh. This approach is motivated by the fact that greedy simplification methods often yield meshes that are poorer than what can be represented with a given number of vertices. Our approach relies on edge swaps and vertex teleports to alter the mesh connectivity, and uses the downhill simplex method to simultaneously improve vertex positions and surface attributes. Note that this is not a simplification method--the vertex count remains the same throughout the optimization. At all stages of the optimization the changes are guided by a metric that measures the differences between rendered versions of the original model and the low vertex count mesh. This method creates meshes that are geometrically faithful to the original model. Moreover, the method takes into account more subtle aspects of a model such as surface shading or whether cracks are visible between two interpenetrating parts of the model.
INCISIONAL HERNIA - ONLAY VS SUBLAY MESH HERNIOPLAS T Y
Ravi Kamal Kumar
2015-04-01
Full Text Available BACKGROUND : Incisional hernia is a common surgical problem. Anatomical repair of hernia is now out of vogue. Polypropylene mesh repair has now become accepted. In open mesh repair of incisional hernia cases the site of placement of mesh is still debated. Some surgeo ns favour the onlay repair and others use sublay or retro - rectus plane for deployment of the mesh. AIM: The aim of the study is to examine the pros and cons of both the techniques and find the better one. METHODS : A prospective study was conducted of 37 ca ses of incisional hernia admitted in Govt. General Hospital Guntur from Jan 2012 to Dec 2013. 20 of the cases underwent open mesh repair by onlay method whereas 17 cases underwent open mesh repair by the sublay i . e . retrorectus placement of the mesh. Obser vations were made regarding time taken for both types of repairs, post - operative complications like flap necrosis, wound seroma, wound infecton, postoperative ileus etc., after discharge from the hospital the cases were followed up in the OPD upto Dec . 201 4 and any complications and recurrences were noted. OBSERVATIONS AND RESULTS: Most of the cases (75% were female and the incisional hernias were in the lower abdomen. The time taken for the surgery is more in the sublay group and the postoperative pain score is also more in the sublay group. But the wound complications like wound infection and flap necrosis were more (25 - 30% in the onlay group. Also one case (5% developed recurrence. Though the time taken for the surgery and the skill needed is more fo r the sublay group the wound complications are acceptable in the sublay group. Also there are no recurrences observed in the sublay group. But no statistically significant difference (p<0 . 7 is detected when all the post - operative complications are taken t ogether between the sublay and onlay repair techniques . CONCLUSIONS: Although it can be argued, theoretically and by the wound complication rate, in favour of
Tracking in anatomic pathology.
Pantanowitz, Liron; Mackinnon, Alexander C; Sinard, John H
2013-12-01
Bar code-based tracking solutions, long present in clinical pathology laboratories, have recently made an appearance in anatomic pathology (AP) laboratories. Tracking of AP "assets" (specimens, blocks, slides) can enhance laboratory efficiency, promote patient safety, and improve patient care. Routing of excess clinical material into research laboratories and biorepositories are other avenues that can benefit from tracking of AP assets. Implementing tracking is not as simple as installing software and turning it on. Not all tracking solutions are alike. Careful analysis of laboratory workflow is needed before implementing tracking to assure that this solution will meet the needs of the laboratory. Such analysis will likely uncover practices that may need to be modified before a tracking system can be deployed. Costs that go beyond simply that of purchasing software will be incurred and need to be considered in the budgeting process. Finally, people, not technology, are the key to assuring quality. Tracking will require significant changes in workflow and an overall change in the culture of the laboratory. Preparation, training, buy-in, and accountability of the people involved are crucial to the success of this process. This article reviews the benefits, available technology, underlying principles, and implementation of tracking solutions for the AP and research laboratory. PMID:23634908
Method and system for mesh network embedded devices
Wang, Ray (Inventor)
2009-01-01
A method and system for managing mesh network devices. A mesh network device with integrated features creates an N-way mesh network with a full mesh network topology or a partial mesh network topology.
Optimal Cache-Oblivious Mesh Layouts
Bender, Michael A.; Kuszmaul, Bradley C.; Teng, Shang-Hua; Wang, Kebin
2007-01-01
A mesh is a graph that divides physical space into regularly-shaped regions. Meshes computations form the basis of many applications, e.g. finite-element methods, image rendering, and collision detection. In one important mesh primitive, called a mesh update, each mesh vertex stores a value and repeatedly updates this value based on the values stored in all neighboring vertices. The performance of a mesh update depends on the layout of the mesh in memory. This paper shows how to find a memory...
User Manual for the PROTEUS Mesh Tools
Smith, Micheal A. [Argonne National Lab. (ANL), Argonne, IL (United States); Shemon, Emily R. [Argonne National Lab. (ANL), Argonne, IL (United States)
2015-06-01
This report describes the various mesh tools that are provided with the PROTEUS code giving both descriptions of the input and output. In many cases the examples are provided with a regression test of the mesh tools. The most important mesh tools for any user to consider using are the MT_MeshToMesh.x and the MT_RadialLattice.x codes. The former allows the conversion between most mesh types handled by PROTEUS while the second allows the merging of multiple (assembly) meshes into a radial structured grid. Note that the mesh generation process is recursive in nature and that each input specific for a given mesh tool (such as .axial or .merge) can be used as “mesh” input for any of the mesh tools discussed in this manual.
User Manual for the PROTEUS Mesh Tools
This report describes the various mesh tools that are provided with the PROTEUS code giving both descriptions of the input and output. In many cases the examples are provided with a regression test of the mesh tools. The most important mesh tools for any user to consider using are the MTMeshToMesh.x and the MTRadialLattice.x codes. The former allows the conversion between most mesh types handled by PROTEUS while the second allows the merging of multiple (assembly) meshes into a radial structured grid. Note that the mesh generation process is recursive in nature and that each input specific for a given mesh tool (such as .axial or .merge) can be used as ''mesh'' input for any of the mesh tools discussed in this manual.
An Adaptive Mesh Algorithm: Mesh Structure and Generation
Scannapieco, Anthony J. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2016-06-21
The purpose of Adaptive Mesh Refinement is to minimize spatial errors over the computational space not to minimize the number of computational elements. The additional result of the technique is that it may reduce the number of computational elements needed to retain a given level of spatial accuracy. Adaptive mesh refinement is a computational technique used to dynamically select, over a region of space, a set of computational elements designed to minimize spatial error in the computational model of a physical process. The fundamental idea is to increase the mesh resolution in regions where the physical variables are represented by a broad spectrum of modes in k-space, hence increasing the effective global spectral coverage of those physical variables. In addition, the selection of the spatially distributed elements is done dynamically by cyclically adjusting the mesh to follow the spectral evolution of the system. Over the years three types of AMR schemes have evolved; block, patch and locally refined AMR. In block and patch AMR logical blocks of various grid sizes are overlaid to span the physical space of interest, whereas in locally refined AMR no logical blocks are employed but locally nested mesh levels are used to span the physical space. The distinction between block and patch AMR is that in block AMR the original blocks refine and coarsen entirely in time, whereas in patch AMR the patches change location and zone size with time.
Performance of a mixed-mesh Godunov-based flood inundation model
Kim, B.; Sanders, B. F.; Kim, H.; Famiglietti, J. S.
2011-12-01
Godunov-based finite volume models for solving the shallow-water equations, which are seeing increasing use in hydrology for flood inundation modeling, have almost exclusively adopted either structured meshes of quadrilateral cells or unstructured meshes of triangular cells but have not focused on mixing quadrilateral and triangular cells as is common with finite element models. Triangular meshing is advantageous in complex topography arising from channel junctions, meandering channels, and artificial structures because powerful Delaunay mesh generators easily accommodate internal and external boundary constraints while retaining mesh quality attributes such as area and angle properties. On the other hand, quadrilateral meshing is advantages in the absence of boundary constraints because efficient and low-overhead mesh designs such as Cartesian grids can be used. Motivated by the need for efficient and accurate modeling of flood inundation including channel flows, overbank flows, and overtopping processes, a mixed-mesh version of the BreZo flood inundation model is presented and its ease-of-use and performance in a series of laboratory and field scale test problems is examined in comparison to versions that use either triangular or quadrilateral cells exclusively. We find that mixed meshes are not as easily prepared as triangular meshes for commonly encountered study-site geometries, but execution times and memory requirements for a similar level of accuracy are reduced. The challenges of designing high quality meshes for flood modeling highlights the need for a new class of mesh generators that can scan high resolution topographic data (e.g., lidar) for critical features (e.g., channels and levees) and automate mixed-mesh generation and parameterization.
CCD Photometry of bright stars using objective wire mesh
Obtaining accurate photometry of bright stars from the ground remains problematic due to the danger of overexposing the target and/or the lack of suitable nearby comparison stars. The century-old method of using objective wire mesh to produce multiple stellar images seems promising for the precise CCD photometry of such stars. Furthermore, our tests on β Cep and its comparison star, differing by 5 mag, are very encouraging. Using a CCD camera and a 20 cm telescope with the objective covered by a plastic wire mesh, in poor weather conditions, we obtained differential photometry with a precision of 4.5 mmag per two minute exposure. Our technique is flexible and may be tuned to cover a range as big as 6-8 mag. We discuss the possibility of installing a wire mesh directly in the filter wheel.
Adaptive anisotropic meshing for steady convection-dominated problems
Nguyen, Hoa [Tulane University; Gunzburger, Max [Florida State University; Ju, Lili [University of South Carolina; Burkardt, John [Florida State University
2009-01-01
Obtaining accurate solutions for convection–diffusion equations is challenging due to the presence of layers when convection dominates the diffusion. To solve this problem, we design an adaptive meshing algorithm which optimizes the alignment of anisotropic meshes with the numerical solution. Three main ingredients are used. First, the streamline upwind Petrov–Galerkin method is used to produce a stabilized solution. Second, an adapted metric tensor is computed from the approximate solution. Third, optimized anisotropic meshes are generated from the computed metric tensor by an anisotropic centroidal Voronoi tessellation algorithm. Our algorithm is tested on a variety of two-dimensional examples and the results shows that the algorithm is robust in detecting layers and efficient in avoiding non-physical oscillations in the numerical approximation.
CCD Photometry of bright stars using objective wire mesh
Kamiński, Krzysztof; Zgórz, Marika [Astronomical Observatory Institute, Faculty of Physics, A. Mickiewicz University, Słoneczna 36, 60-286 Poznań (Poland); Schwarzenberg-Czerny, Aleksander, E-mail: chrisk@amu.edu.pl [Copernicus Astronomical Centre, ul. Bartycka 18, PL 00-716 Warsaw (Poland)
2014-06-01
Obtaining accurate photometry of bright stars from the ground remains problematic due to the danger of overexposing the target and/or the lack of suitable nearby comparison stars. The century-old method of using objective wire mesh to produce multiple stellar images seems promising for the precise CCD photometry of such stars. Furthermore, our tests on β Cep and its comparison star, differing by 5 mag, are very encouraging. Using a CCD camera and a 20 cm telescope with the objective covered by a plastic wire mesh, in poor weather conditions, we obtained differential photometry with a precision of 4.5 mmag per two minute exposure. Our technique is flexible and may be tuned to cover a range as big as 6-8 mag. We discuss the possibility of installing a wire mesh directly in the filter wheel.
Relativistic polarizabilities with the Lagrange-mesh method
Filippin, Livio; Baye, Daniel
2016-01-01
Relativistic dipolar to hexadecapolar polarizabilities of the ground state and some excited states of hydrogenic atoms are calculated by using numerically exact energies and wave functions obtained from the Dirac equation with the Lagrange-mesh method. This approach is an approximate variational method taking the form of equations on a grid because of the use of a Gauss quadrature approximation. The partial polarizabilities conserving the absolute value of the quantum number $\\kappa$ are also numerically exact with small numbers of mesh points. The ones where $|\\kappa|$ changes are very accurate when using three different meshes for the initial and final wave functions and for the calculation of matrix elements. The polarizabilities of the $n=2$ excited states of hydrogenic atoms are also studied with a separate treatment of the final states that are degenerate at the nonrelativistic approximation. The method provides high accuracies for polarizabilities of a particle in a Yukawa potential and is applied to a...
Anatomical variations of paranasal sinuses at multislice computed tomography: what to look for
Miranda, Christiana Maia Nobre Rocha de; Maranhao, Carol Pontes de Miranda [Clinica de Medicina Nuclear e Radiologia de Maceio (Medradius), Maceio, AL (Brazil). Setor de Tomografia Computadorizada; Arraes, Fabiana Maia Nobre Rocha [Clinica Sinus, Maceio, AL (Brazil); Padilha, Igor Gomes; Farias, Lucas de Padua Gomes de; Jatoba, Mayara Stephanie de Araujo; Andrade, Anna Carolina Mendonca de; Padilha, Bruno Gomes [Universidade Federal de Alagoas (UFAL), Maceio, AL (Brazil)
2011-07-15
Multislice computed tomography is currently the imaging modality of choice for evaluating paranasal sinuses and adjacent structures. Such a method has been increasingly utilized in the assessment of anatomical variations, allowing their accurate identification with high anatomical details. Some anatomical variations may predispose to sinusal diseases, constituting areas of high risk for injuries and complications during surgical procedures. Therefore, the recognition of such variations is critical in the preoperative evaluation for endoscopic surgery. (author)
Geostrophic balance preserving interpolation in mesh adaptive shallow-water ocean modelling
Maddison, James R; Farrell, Patrick E
2010-01-01
The accurate representation of geostrophic balance is an essential requirement for numerical modelling of geophysical flows. Significant effort is often put into the selection of accurate or optimal balance representation by the discretisation of the fundamental equations. The issue of accurate balance representation is particularly challenging when applying dynamic mesh adaptivity, where there is potential for additional imbalance injection when interpolating to new, optimised meshes. In the context of shallow-water modelling, we present a new method for preservation of geostrophic balance when applying dynamic mesh adaptivity. This approach is based upon interpolation of the Helmholtz decomposition of the Coriolis acceleration. We apply this in combination with a discretisation for which states in geostrophic balance are exactly steady solutions of the linearised equations on an f-plane; this method guarantees that a balanced and steady flow on a donor mesh remains balanced and steady after interpolation on...
A coarse-mesh nodal method-diffusive-mesh finite difference method
Modern nodal methods have been successfully used for conventional light water reactor core analyses where the homogenized, node average cross sections (XSs) and the flux discontinuity factors (DFs) based on equivalence theory can reliably predict core behavior. For other types of cores and other geometries characterized by tightly-coupled, heterogeneous core configurations, the intranodal flux shapes obtained from a homogenized nodal problem may not accurately portray steep flux gradients near fuel assembly interfaces or various reactivity control elements. This may require extreme values of DFs (either very large, very small, or even negative) to achieve a desired solution accuracy. Extreme values of DFs, however, can disrupt the convergence of the iterative methods used to solve for the node average fluxes, and can lead to a difficulty in interpolating adjacent DF values. Several attempts to remedy the problem have been made, but nothing has been satisfactory. A new coarse-mesh nodal scheme called the Diffusive-Mesh Finite Difference (DMFD) technique, as contrasted with the coarse-mesh finite difference (CMFD) technique, has been developed to resolve this problem. This new technique and the development of a few-group, multidimensional kinetics computer program are described in this paper
Connectivity editing for quadrilateral meshes
Peng, Chihan
2011-12-12
We propose new connectivity editing operations for quadrilateral meshes with the unique ability to explicitly control the location, orientation, type, and number of the irregular vertices (valence not equal to four) in the mesh while preserving sharp edges. We provide theoretical analysis on what editing operations are possible and impossible and introduce three fundamental operations to move and re-orient a pair of irregular vertices. We argue that our editing operations are fundamental, because they only change the quad mesh in the smallest possible region and involve the fewest irregular vertices (i.e., two). The irregular vertex movement operations are supplemented by operations for the splitting, merging, canceling, and aligning of irregular vertices. We explain how the proposed highlevel operations are realized through graph-level editing operations such as quad collapses, edge flips, and edge splits. The utility of these mesh editing operations are demonstrated by improving the connectivity of quad meshes generated from state-of-art quadrangulation techniques. © 2011 ACM.
Connectivity editing for quadrilateral meshes
Peng, Chihan
2011-12-01
We propose new connectivity editing operations for quadrilateral meshes with the unique ability to explicitly control the location, orientation, type, and number of the irregular vertices (valence not equal to four) in the mesh while preserving sharp edges. We provide theoretical analysis on what editing operations are possible and impossible and introduce three fundamental operations to move and re-orient a pair of irregular vertices. We argue that our editing operations are fundamental, because they only change the quad mesh in the smallest possible region and involve the fewest irregular vertices (i.e., two). The irregular vertex movement operations are supplemented by operations for the splitting, merging, canceling, and aligning of irregular vertices. We explain how the proposed high-level operations are realized through graph-level editing operations such as quad collapses, edge flips, and edge splits. The utility of these mesh editing operations are demonstrated by improving the connectivity of quad meshes generated from state-of-art quadrangulation techniques.
Mesh Resolution Effect on 3D RANS Turbomachinery Flow Simulations
Yershov, Sergiy
2016-01-01
The paper presents the study of the effect of a mesh refinement on numerical results of 3D RANS computations of turbomachinery flows. The CFD solver F, which based on the second-order accurate ENO scheme, is used in this study. The simplified multigrid algorithm and local time stepping permit decreasing computational time. The flow computations are performed for a number of turbine and compressor cascades and stages. In all flow cases, the successively refined meshes of H-type with an approximate orthogonalization near the solid walls were generated. The results obtained are compared in order to estimate their both mesh convergence and ability to resolve the transonic flow pattern. It is concluded that for thorough studying the fine phenomena of the 3D turbomachinery flows, it makes sense to use the computational meshes with the number of cells from several millions up to several hundred millions per a single turbomachinery blade channel, while for industrial computations, a mesh of about or less than one mil...
Computational mesh generation for vascular structures with deformable surfaces
Computational blood flow and vessel wall mechanics simulations for vascular structures are becoming an important research tool for patient-specific surgical planning and intervention. An important step in the modelling process for patient-specific simulations is the creation of the computational mesh based on the segmented geometry. Most known solutions either require a large amount of manual processing or lead to a substantial difference between the segmented object and the actual computational domain. We have developed a chain of algorithms that lead to a closely related implementation of image segmentation with deformable models and 3D mesh generation. The resulting processing chain is very robust and leads both to an accurate geometrical representation of the vascular structure as well as high quality computational meshes. The chain of algorithms has been tested on a wide variety of shapes. A benchmark comparison of our mesh generation application with five other available meshing applications clearly indicates that the new approach outperforms the existing methods in the majority of cases. (orig.)
Zheng, J.; Zhu, J.; Wang, Z.; Fang, F.; Pain, C. C.; Xiang, J.
2015-06-01
A new anisotropic hr-adaptive mesh technique has been applied to modelling of multiscale transport phenomena, which is based on a discontinuous Galerkin/control volume discretization on unstructured meshes. Over existing air quality models typically based on static-structured grids using a locally nesting technique, the advantage of the anisotropic hr-adaptive model has the ability to adapt the mesh according to the evolving pollutant distribution and flow features. That is, the mesh resolution can be adjusted dynamically to simulate the pollutant transport process accurately and effectively. To illustrate the capability of the anisotropic adaptive unstructured mesh model, three benchmark numerical experiments have been setup for two-dimensional (2-D) transport phenomena. Comparisons have been made between the results obtained using uniform resolution meshes and anisotropic adaptive resolution meshes.
Gosselin, Marie-Christine; Neufeld, Esra; Moser, Heidi; Huber, Eveline; Farcito, Silvia; Gerber, Livia; Jedensjö, Maria; Hilber, Isabel; Di Gennaro, Fabienne; Lloyd, Bryn; Cherubini, Emilio; Szczerba, Dominik; Kainz, Wolfgang; Kuster, Niels
2014-09-01
The Virtual Family computational whole-body anatomical human models were originally developed for electromagnetic (EM) exposure evaluations, in particular to study how absorption of radiofrequency radiation from external sources depends on anatomy. However, the models immediately garnered much broader interest and are now applied by over 300 research groups, many from medical applications research fields. In a first step, the Virtual Family was expanded to the Virtual Population to provide considerably broader population coverage with the inclusion of models of both sexes ranging in age from 5 to 84 years old. Although these models have proven to be invaluable for EM dosimetry, it became evident that significantly enhanced models are needed for reliable effectiveness and safety evaluations of diagnostic and therapeutic applications, including medical implants safety. This paper describes the research and development performed to obtain anatomical models that meet the requirements necessary for medical implant safety assessment applications. These include implementation of quality control procedures, re-segmentation at higher resolution, more-consistent tissue assignments, enhanced surface processing and numerous anatomical refinements. Several tools were developed to enhance the functionality of the models, including discretization tools, posing tools to expand the posture space covered, and multiple morphing tools, e.g., to develop pathological models or variations of existing ones. A comprehensive tissue properties database was compiled to complement the library of models. The results are a set of anatomically independent, accurate, and detailed models with smooth, yet feature-rich and topologically conforming surfaces. The models are therefore suited for the creation of unstructured meshes, and the possible applications of the models are extended to a wider range of solvers and physics. The impact of these improvements is shown for the MRI exposure of an adult
GRChombo : Numerical Relativity with Adaptive Mesh Refinement
Clough, Katy; Finkel, Hal; Kunesch, Markus; Lim, Eugene A; Tunyasuvunakool, Saran
2015-01-01
Numerical relativity has undergone a revolution in the past decade. With a well-understood mathematical formalism, and full control over the gauge modes, it is now entering an era in which the science can be properly explored. In this work, we introduce GRChombo, a new numerical relativity code written to take full advantage of modern parallel computing techniques. GRChombo's features include full adaptive mesh refinement with block structured Berger-Rigoutsos grid generation which supports non-trivial "many-boxes-in-many-boxes" meshing hierarchies, and massive parallelism through the Message Passing Interface (MPI). GRChombo evolves the Einstein equation with the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. We show that GRChombo passes all the standard "Apples-to-Apples" code comparison tests. We also show that it can stably and accurately evolve vacuum black hole spacetimes such as binary black hole mergers, and non-vacuum spacetimes such as scalar collapses into b...
Application of the VOF method based on unstructured quadrilateral mesh
JI Chun-ning; SHI Ying
2008-01-01
To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any complicated boundary treatment or artificial diffusion, this method treated curved boundaries directly by utilizing the inherent merit of unstructured mesh in fitting curves. The PLIC (Piecewise Linear Interface Calculation) method was adopted to obtain a second-order accurate linearized reconstruction approximation and the MLER (Modified Lagrangian-Eulerian Re-map) method was introduced to advect fluid volumes on unstructured mesh. Moreover, an analytical relation for the interface's line constant vs. the volume clipped by the interface was developed so as to improve the method's efficiency. To validate this method, a comprehensive series of large straining advection tests were performed. Numerical results provide convincing evidences for the method's high volume conservative accuracy and second-order shape error convergence rate. Also, a dramatic improvement on computational accuracy over its unstructured triangular mesh counterpart is checked.
Compressive VOF method with skewness correction to capture sharp interfaces on arbitrary meshes
Denner, Fabian; van Wachem, Berend G. M.
2014-12-01
The accurate and efficient modelling of two-phase flows is at present mostly limited to structured, unskewed meshes, due to the additional topological and numerical complexity of arbitrary, unstructured meshes. Compressive VOF methods which discretize the interface advection with algebraic differencing schemes are computationally efficient and inherently applicable to arbitrary meshes. However, compressive VOF methods evidently suffer severely from numerical diffusion on meshes with topological skewness. In this paper we present a compressive VOF method using a state-of-the-art donor-acceptor advection scheme which includes novel modifications to substantially reduce numerical diffusion on arbitrary meshes without adding computational complexity. The new methodology accurately captures evolving interfaces on any arbitrary, non-overlapping mesh and conserves mass within the limits of the applied solver tolerance. A thorough validation of the presented methods is conducted, examining the pure advection of the interface indicator function as well as the application to evolving interfaces with surface tension. Crucially, the results on equidistant Cartesian and arbitrary tetrahedral meshes are shown to be comparable and accurate.
Efficient Packet Forwarding in Mesh Network
Kanrar, Soumen
2012-01-01
Wireless Mesh Network (WMN) is a multi hop low cost, with easy maintenance robust network providing reliable service coverage. WMNs consist of mesh routers and mesh clients. In this architecture, while static mesh routers form the wireless backbone, mesh clients access the network through mesh routers as well as directly meshing with each other. Different from traditional wireless networks, WMN is dynamically self-organized and self-configured. In other words, the nodes in the mesh network automatically establish and maintain network connectivity. Over the years researchers have worked, to reduce the redundancy in broadcasting packet in the mesh network in the wireless domain for providing reliable service coverage, the source node deserves to broadcast or flood the control packets. The redundant control packet consumes the bandwidth of the wireless medium and significantly reduces the average throughput and consequently reduces the overall system performance. In this paper I study the optimization problem in...
Watanabe,Toyohiko
2012-02-01
Full Text Available Polypropylene mesh implants for the correction of pelvic organ prolapse (POP are now available in Japan. We developed an innovative approach for correcting POP by placing polypropylene mesh transvaginally with laparoscopic assistance. From June 2007 through March 2010, sixteen consecutive patients with symptomatic stage 2 or 3 pelvic organ prolapse underwent the laparoscopic-assisted tension-free vaginal mesh procedure at Okayama University Hospital. All patients were evaluated before and at 1, 3, 6, and 12 months after surgery. Female sexual function was also evaluated with the Female Sexual Function Index (FSFI. The procedure was performed successfully without significant complications. Fifteen of 16 patients were considered anatomically cured (93.8% at 12 months postoperatively. One patient with a recurrent stage 3 vaginal vault prolapse required sacral colpopexy six months postoperatively. Total FSFI scores improved significantly from 10.3±1.3 at baseline to 18.0±1.2 at 12 months after surgery. The laparoscopic-assisted trans-vaginal mesh is a safe, effective, and simple procedure for POP repairs. The procedure not only restores anatomic relationships but also improves sexual function.
ANATOMICAL PROPERTIES OF PLANTAGO ARENARIA
Nicoleta IANOVICI; SINITEAN, Adrian; Aurel FAUR
2011-01-01
Psammophytes are marked by a number of adaptations that enable them to exist in the hard environmental conditions of the sand habitats. In this study, the anatomical characteristics of Plantago arenaria were examined. Studies were conducted to assess the diversity of anatomical adaptations of vegetative organs in this taxa. Results are presented with original photographs. The analysis of leaf anatomy in P. arenaria showed that the leaves contained a contained xeromorphic traits. Arbuscular my...
Multiple scale mesh free analysis
Recent developments of mesh free and multi-scale methods and their applications in applied mechanics are surveyed. Three major methodologies are reviewed. First, smoothed particle hydrodynamics (SPH) is discussed as a representative of a non-local kernel, strong form collocation approach. Second, mesh-free Galerkin methods, which have been active research area in recent years, are reviewed. Third, some applications of molecular dynamics (MD) in applied mechanics are discussed. The emphases of this survey are placed on simulations of finite deformations, fracture, shear bands, multi-scale methods, and nano-scale mechanics. Refs. 13 (author)
DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks
Duffell, Paul C
2016-01-01
This work presents the publicly available moving-mesh magnetohydrodynamics code DISCO. DISCO is efficient and accurate at evolving orbital fluid motion in two and three dimensions, especially at high Mach number. DISCO employs a moving-mesh approach utilizing a dynamic cylindrical mesh that can shear azimuthally to follow the orbital motion of the gas. The moving mesh removes diffusive advection errors and allows for longer timesteps than a static grid. Magnetohydrodynamics is implemented in DISCO using an HLLD Riemann solver and a novel constrained transport scheme which is compatible with the mesh motion. DISCO is tested against a wide variety of problems, which are designed to test its stability, accuracy and scalability. In addition, several magnetohydrodynamics tests are performed which demonstrate the accuracy and stability of the new constrained transport approach, including two tests of the magneto-rotational instability (MRI); one testing the linear growth rate and the other following the instability...
A Numerical Study of Blowup in the Harmonic Map Heat Flow Using the MMPDE Moving Mesh Method
Haynes, R.D.; Huang, W.; Zegeling, P.A.
2013-01-01
The numerical solution of the harmonic heat map flow problems with blowup in finite or infinite time is considered using an adaptive moving mesh method. A properly chosen monitor function is derived so that the moving mesh method can be used to simulate blowup and produce accurate blowup profiles wh
Development and verification of unstructured adaptive mesh technique with edge compatibility
In the design study of the large-sized sodium-cooled fast reactor (JSFR), one key issue is suppression of gas entrainment (GE) phenomena at a gas-liquid interface. Therefore, the authors have been developed a high-precision CFD algorithm to evaluate the GE phenomena accurately. The CFD algorithm has been developed on unstructured meshes to establish an accurate modeling of JSFR system. For two-phase interfacial flow simulations, a high-precision volume-of-fluid algorithm is employed. It was confirmed that the developed CFD algorithm could reproduce the GE phenomena in a simple GE experiment. Recently, the authors have been developed an important technique for the simulation of the GE phenomena in JSFR. That is an unstructured adaptive mesh technique which can apply fine cells dynamically to the region where the GE occurs in JSFR. In this paper, as a part of the development, a two-dimensional unstructured adaptive mesh technique is discussed. In the two-dimensional adaptive mesh technique, each cell is refined isotropically to reduce distortions of the mesh. In addition, connection cells are formed to eliminate the edge incompatibility between refined and non-refined cells. The two-dimensional unstructured adaptive mesh technique is verified by solving well-known lid-driven cavity flow problem. As a result, the two-dimensional unstructured adaptive mesh technique succeeds in providing a high-precision solution, even though poor-quality distorted initial mesh is employed. In addition, the simulation error on the two-dimensional unstructured adaptive mesh is much less than the error on the structured mesh with a larger number of cells. (author)
[Sigismund Laskowski and his anatomical preparations technique].
Gryglewski, Ryszard W
2015-01-01
Fixation of the entire bodies or individual organs, and later as well tissues and cellular structures, was and still is often a challenge for anatomists and histologists. Technique that combines extensive knowledge of natural sciences, as well as technical skills, was by those best researchers as Frederik Ruysch, brought to perfection. Preparations, if done with care and talent, are really propelling progress in anatomical studies and determining the quality of education for medical students and young physicians. And as it is true for many of today's medical disciplines and natural sciences, the nineteenth century was in many ways a breaking point for preparatory techniques in the realm of anatomy and histology. Among those who have achieved success, earning notoriety during their lifetime and often going into the annals of European most distinguished scholars were some Polish names: Louis Maurice Hirschfeld, whose preparations of the nervous system earned him well-deserved, international fame, Louis Charles Teichmann, who was the very first so precisely describing the lymphatic system and a creator of unique injection mass, Henry Kadyi, known for his outstanding preparations, especially of vascular system. Henry Frederick Hoyer sen., who was one of the first to use formalin regularly for accurate microscopic preparations, is seen by many as the founder of the Polish histology. In this group of innovators and precursors of modern preparation techniques place should be reserved for Zygmunt (Sigismund) Laskowski, Polish patriot, fighting in January Uprising, later an immigrant, a professor at the university sequentially Paris and Geneva. Acclaimed author of anatomical tables and certainly creator of one of the groundbreaking techniques in anatomical preparations. Based after many years of research on the simple glycerine-phenol mixture achieved excellent results both in fixation of entire bodies and organs or tissues. Quality of those preparations was as high and
Numerical solution of the diffusion equation plays a key role in the study of inertial confinement fusion (ICF). In this paper, based on the global support operator method, a flux-based scheme is proposed. The scheme has local stencil with second-order accuracy both in space and time. For strongly distorted meshes, a procedure of normal direction fix is adopted with proper methods for the computation of corner volume weights, which obtains accurate discretization of the face flux. Numerical experiments show that the scheme can obtain accurate solution for linear problems on non-convex meshes. The method has second-order spatial and temporal accuracy on non-smooth meshes. The method can also preserve the symmetry well and can be extended to the three dimensional unstructured meshes. (authors)
6th International Meshing Roundtable '97
White, D.
1997-09-01
The goal of the 6th International Meshing Roundtable is to bring together researchers and developers from industry, academia, and government labs in a stimulating, open environment for the exchange of technical information related to the meshing process. In the pas~ the Roundtable has enjoyed significant participation born each of these groups from a wide variety of countries. The Roundtable will consist of technical presentations from contributed papers and abstracts, two invited speakers, and two invited panels of experts discussing topics related to the development and use of automatic mesh generation tools. In addition, this year we will feature a "Bring Your Best Mesh" competition and poster session to encourage discussion and participation from a wide variety of mesh generation tool users. The schedule and evening social events are designed to provide numerous opportunities for informal dialog. A proceedings will be published by Sandia National Laboratories and distributed at the Roundtable. In addition, papers of exceptionally high quaIity will be submitted to a special issue of the International Journal of Computational Geometry and Applications. Papers and one page abstracts were sought that present original results on the meshing process. Potential topics include but are got limited to: Unstructured triangular and tetrahedral mesh generation Unstructured quadrilateral and hexahedral mesh generation Automated blocking and structured mesh generation Mixed element meshing Surface mesh generation Geometry decomposition and clean-up techniques Geometry modification techniques related to meshing Adaptive mesh refinement and mesh quality control Mesh visualization Special purpose meshing algorithms for particular applications Theoretical or novel ideas with practical potential Technical presentations from industrial researchers.
Damiani, G R; Riva, D; Pellegrino, A; Gaetani, M; Tafuri, S; Turoli, D; Croce, P; Loverro, G
2016-04-01
117 women with severe pelvic organ prolapse (POP; stage > 2) were enrolled to elucidate a 24-month outcome of POP surgery, using conventional or mesh repair with 3 techniques. 59 patients underwent conventional repair and 58 underwent mesh repair. Two types of mesh were used: a trocar-guided transobturator polypropylene (Avaulta, Bard Inc.) and a porcine dermis mesh (Pelvisoft, Bard Inc.). Women with recurrences, who underwent previous unsuccessful conventional repair, were randomised. Primary outcome was the evaluation of anatomic failures (prolapse stage > 1) in treated and untreated compartments. Anatomic failure was observed in 11 of 58 patients (19%; CI 8.9-29) in the mesh group and in 16 of 59 patients (27.1%; p value = 0.3) in the conventional group. 9 of 11 failures in the mesh group (15.5%; CI 6.2-24.8) were observed in the untreated compartment (de novo recurrences), 14.3% in Pelvisoft and 16.7% in Avaulta arm, while only 1 recurrence in the untreated compartment (1.7%) was observed in the conventional group (odds ratio 10.6, p = 0.03). PMID:26492359
On-line residual capacity estimation for resource allocation in wireless mesh networks
Sarıkaya, Yunus; Sarikaya, Yunus
2008-01-01
Contention-based multi access scheme of 802.11 based wireless mesh networks imposes difficulties in achieving predictable service quality in multi-hop networks. In order to offer effective advanced network services such as flow admission control or load balancing, the residual capacity of the wireless links should be accurately estimated. In this work, we propose and validate an algorithm for the residual bandwidth of wireless mesh network. By collecting transmission statistics from the nearb...
ANATOMICAL PROPERTIES OF PLANTAGO ARENARIA
Nicoleta IANOVICI
2011-01-01
Full Text Available Psammophytes are marked by a number of adaptations that enable them to exist in the hard environmental conditions of the sand habitats. In this study, the anatomical characteristics of Plantago arenaria were examined. Studies were conducted to assess the diversity of anatomical adaptations of vegetative organs in this taxa. Results are presented with original photographs. The analysis of leaf anatomy in P. arenaria showed that the leaves contained a contained xeromorphic traits. Arbuscular mycorrhizal symbiosis seems to be critical for their survival.
Feature-preserving mesh denoising via normal guided quadric error metrics
Yu, Jinze; Wei, Mingqiang; Qin, Jing; Wu, Jianhuang; Heng, Pheng-Ann
2014-11-01
While modern optical and laser 3D scanners can generate high accuracy mesh models, to largely avoid their introducing noise which prohibits practical applications still results in high cost. Thus, optimizing noisy meshes while preserving their geometric details is necessary for production, which still remains as challenging work. In this paper we propose a novel and efficient two-stage feature-preserving mesh denoising framework which can remove noise while preserving fine features of a surface mesh. We improve the capability of feature preservation of our vertex updating scheme by employing an extension of the quadric error metrics (QEM), which can track and minimize updating errors and hence well preserve the overall shape as well as detailed features of a mesh. We further leverage vertex normals to guide the vertex updating process, as the normal field of a mesh reflects the geometry of the underlying surface. In addition, to obtain a more accurate normal field to guide vertex updating, we develop an improved normal filter by integrating advantages of existing filters. Compared with traditional gradient descent based schemes, our method performs better on challenging regions with rich geometric features. Moreover, a local entropy metric is proposed to measure stability of a mesh and the effectiveness of vertex updating algorithms. Qualitative and quantitative experiments demonstrate that our approach can effectively remove noise from noisy meshes while preserving or recovering geometrical features of original objects.
Simple mesh stent placement for treating intracranial aneurysm: progress in research
Treatment of an aneurysm with a mesh stent alone becomes an emerging technique. The mechanism involves the mesh stents, when it crosses the neck of an aneurysm could change the internal circulation and induce the formation of stable thrombus and in turn assists the growth of neointima for anatomical healing of the aneurysmal neck. The mesh stent technique aimed at vascular reconstruction with expected curing effect for aneurysm together as a simple and safe way to keep the patency of the small arterial branches, just contrary to the covered stent. This technique couldn't be carried out practically in wide scale because of immaturity, therefore we give a comprehesive review in the progress of this field. (authors)
[Antique anatomical collections for contemporary museums].
Nesi, Gabriella; Santi, Raffaella
2013-01-01
Anatomy and Pathology Museum collections display a great biological value and offer unique samples for research purposes. Pathological specimens may be investigated by means of modern radiological and molecular biology techniques in order to provide the etiological background of disease, with relevance to present-day knowledge. Meanwhile, historical resources provide epidemiologic data regarding the socio-economic conditions of the resident populations, the more frequently encountered illnesses and dietary habits. These multidisciplinary approaches lead to more accurate diagnoses also allowing new strategies in cataloguing and musealization of anatomical specimens. Further, once these data are gathered, they may constitute the basis of riedited Museum catalogues feasible to be digitalized and displayed via the Web. PMID:25807710
21st International Meshing Roundtable
Weill, Jean-Christophe
2013-01-01
This volume contains the articles presented at the 21st International Meshing Roundtable (IMR) organized, in part, by Sandia National Laboratories and was held on October 7–10, 2012 in San Jose, CA, USA. The first IMR was held in 1992, and the conference series has been held annually since. Each year the IMR brings together researchers, developers, and application experts in a variety of disciplines, from all over the world, to present and discuss ideas on mesh generation and related topics. The technical papers in this volume present theoretical and novel ideas and algorithms with practical potential, as well as technical applications in science and engineering, geometric modeling, computer graphics, and visualization.
22nd International Meshing Roundtable
Staten, Matthew
2014-01-01
This volume contains the articles presented at the 22nd International Meshing Roundtable (IMR) organized, in part, by Sandia National Laboratories and was held on Oct 13-16, 2013 in Orlando, Florida, USA. The first IMR was held in 1992, and the conference series has been held annually since. Each year the IMR brings together researchers, developers, and application experts in a variety of disciplines, from all over the world, to present and discuss ideas on mesh generation and related topics. The technical papers in this volume present theoretical and novel ideas and algorithms with practical potential, as well as technical applications in science and engineering, geometric modeling, computer graphics and visualization.
The moving mesh code Shadowfax
Vandenbroucke, Bert
2016-01-01
We introduce the moving mesh code Shadowfax, which can be used to evolve a mixture of gas, subject to the laws of hydrodynamics and gravity, and any collisionless fluid only subject to gravity, such as cold dark matter or stars. The code is written in C++ and its source code is made available to the scientific community under the GNU Affero General Public License. We outline the algorithm and the design of our implementation, and demonstrate its validity through the results of a set of basic test problems, which are also part of the public version. We also compare Shadowfax with a number of other publicly available codes using different hydrodynamical integration schemes, illustrating the advantages and disadvantages of the moving mesh technique.
Confined helium on Lagrange meshes
Baye, Daniel
2015-01-01
The Lagrange-mesh method has the simplicity of a calculation on a mesh and can have the accuracy of a variational method. It is applied to the study of a confined helium atom. Two types of confinement are considered. Soft confinements by potentials are studied in perimetric coordinates. Hard confinement in impenetrable spherical cavities is studied in a system of rescaled perimetric coordinates varying in [0,1] intervals. Energies and mean values of the distances between electrons and between an electron and the helium nucleus are calculated. A high accuracy of 11 to 15 significant figures is obtained with small computing times. Pressures acting on the confined atom are also computed. For sphere radii smaller than 1, their relative accuracies are better than $10^{-10}$. For larger radii up to 10, they progressively decrease to $10^{-3}$, still improving the best literature results.
ROAMing terrain (Real-time Optimally Adapting Meshes)
Duchaineau, M.; Wolinsky, M.; Sigeti, D.E.; Miller, M.C.; Aldrich, C.; Mineev, M.
1997-07-01
Terrain visualization is a difficult problem for applications requiring accurate images of large datasets at high frame rates, such as flight simulation and ground-based aircraft testing using synthetic sensor stimulation. On current graphics hardware, the problem is to maintain dynamic, view-dependent triangle meshes and texture maps that produce good images at the required frame rate. We present an algorithm for constructing triangle meshes that optimizes flexible view-dependent error metrics, produces guaranteed error bounds, achieves specified triangle counts directly, and uses frame-to-frame coherence to operate at high frame rates for thousands of triangles per frame. Our method, dubbed Real-time Optimally Adapting Meshes (ROAM), uses two priority queues to drive split and merge operations that maintain continuous triangulations built from pre-processed bintree triangles. We introduce two additional performance optimizations: incremental triangle stripping and priority-computation deferral lists. ROAM execution time is proportionate to the number of triangle changes per frame, which is typically a few percent of the output mesh size, hence ROAM performance is insensitive to the resolution and extent of the input terrain. Dynamic terrain and simple vertex morphing are supported.
Bercea, Gheorghe-Teodor; McRae, Andrew T. T.; Ham, David A.; Mitchell, Lawrence; Rathgeber, Florian; Nardi, Luigi; Luporini, Fabio; Kelly, Paul H. J.
2016-01-01
We present a generic algorithm for numbering and then efficiently iterating over the data values attached to an extruded mesh. An extruded mesh is formed by replicating an existing mesh, assumed to be unstructured, to form layers of prismatic cells. Applications of extruded meshes include, but are not limited to, the representation of 3D high aspect ratio domains employed by geophysical finite element simulations. These meshes are structured in the extruded direction. The algorithm presented ...
Adaptive Mesh Refinement in CTH
This paper reports progress on implementing a new capability of adaptive mesh refinement into the Eulerian multimaterial shock- physics code CTH. The adaptivity is block-based with refinement and unrefinement occurring in an isotropic 2:1 manner. The code is designed to run on serial, multiprocessor and massive parallel platforms. An approximate factor of three in memory and performance improvements over comparable resolution non-adaptive calculations has-been demonstrated for a number of problems
Mesh networked unattended ground sensors
Colling, Kent; Calcutt, Wade; Winston, Mark; Jones, Barry
2006-05-01
McQ has developed a family of low cost unattended ground sensors that utilize self-configured, mesh network communications for wireless sensing. Intended for use in an urban environment, the area monitored by the sensor system poses a communication challenge. A discussion into the sensor's communication performance and how it affects sensor installation and the operation of the system once deployed is presented.
Image meshing via hierarchical optimization＊
Hao XIE; Ruo-feng TONGS
2016-01-01
Vector graphic, as a kind of geometric representation of raster images, has many advantages, e.g., definition independence and editing facility. A popular way to convert raster images into vector graphics is image meshing, the aim of which is to find a mesh to represent an image as faithfully as possible. For traditional meshing algorithms, the crux of the problem resides mainly in the high non-linearity and non-smoothness of the objective, which makes it difficult to find a desirable optimal solution. To ameliorate this situation, we present a hierarchical optimization algorithm solving the problem from coarser levels to finer ones, providing initialization for each level with its coarser ascent. To further simplify the problem, the original non-convex problem is converted to a linear least squares one, and thus becomes convex, which makes the problem much easier to solve. A dictionary learning framework is used to combine geometry and topology elegantly. Then an alternating scheme is employed to solve both parts. Experiments show that our algorithm runs fast and achieves better results than existing ones for most images.
Image meshing via hierarchical optimization
Hao XIE; Ruo-feng TONG‡
2016-01-01
Vector graphic, as a kind of geometric representation of raster images, has many advantages, e.g., defi nition independence and editing facility. A popular way to convert raster images into vector graphics is image meshing, the aim of which is to fi nd a mesh to represent an image as faithfully as possible. For traditional meshing algorithms, the crux of the problem resides mainly in the high non-linearity and non-smoothness of the objective, which makes it diﬃcult to fi nd a desirable optimal solution. To ameliorate this situation, we present a hierarchical optimization algorithm solving the problem from coarser levels to fi ner ones, providing initialization for each level with its coarser ascent. To further simplify the problem, the original non-convex problem is converted to a linear least squares one, and thus becomes convex, which makes the problem much easier to solve. A dictionary learning framework is used to combine geometry and topology elegantly. Then an alternating scheme is employed to solve both parts. Experiments show that our algorithm runs fast and achieves better results than existing ones for most images.
Paulo Palma
2010-04-01
Full Text Available PURPOSE: This prospective study was performed to achieve visualization of the reestablishment of anatomy after reconstructive surgery in the different pelvic compartments with non-absorbable radiopaque meshes, providing valuable anatomic information for surgeons implanting meshes. MATERIALS AND METHODS: A total of 30 female patients with stress urinary incontinence (SUI, anterior and posterior vaginal wall prolapse, or both underwent surgical repair using radiopaque meshes after written informed consent. Patients with SUI underwent five different surgeries. Patients with anterior vaginal prolapse underwent a procedure using a combined pre-pubic and transobturator mesh, and those with posterior vaginal prolapse underwent posterior slingplasty. Three-dimensional reconstruction using helical CT was performed four weeks postoperatively. RESULTS: In all cases, the mesh was clearly visualized. Transobturator slings were shown at the midurethra, and the anchoring tails perforated the obturator foramen at the safety region. Mini-slings were in the proper place, and computed angiography revealed that the anchoring system was away from the obturator vessels. In patients undergoing procedure for anterior vaginal prolapse, both pre-pubic armpit and obturator slings were clearly seen and the mesh was in the proper position, supporting the bladder base and occluding the distal part of the urogenital hiatus. Transcoccygeal sacropexy revealed indirectly a well-supported "neo rectovaginal fascia" and the anchoring tails at the level of ischial spines. CONCLUSION: Three-dimensional helical tomography images of the female pelvis using radiopaque meshes have a potential role in improving our understanding of pelvic floor reconstructive surgeries. These radiopaque meshes might be the basis of a new investigative methodology.
SHARP/PRONGHORN Interoperability: Mesh Generation
Avery Bingham; Javier Ortensi
2012-09-01
Progress toward collaboration between the SHARP and MOOSE computational frameworks has been demonstrated through sharing of mesh generation and ensuring mesh compatibility of both tools with MeshKit. MeshKit was used to build a three-dimensional, full-core very high temperature reactor (VHTR) reactor geometry with 120-degree symmetry, which was used to solve a neutron diffusion critical eigenvalue problem in PRONGHORN. PRONGHORN is an application of MOOSE that is capable of solving coupled neutron diffusion, heat conduction, and homogenized flow problems. The results were compared to a solution found on a 120-degree, reflected, three-dimensional VHTR mesh geometry generated by PRONGHORN. The ability to exchange compatible mesh geometries between the two codes is instrumental for future collaboration and interoperability. The results were found to be in good agreement between the two meshes, thus demonstrating the compatibility of the SHARP and MOOSE frameworks. This outcome makes future collaboration possible.
Optimizing the geometrical accuracy of curvilinear meshes
Toulorge, Thomas; Remacle, Jean-François
2015-01-01
This paper presents a method to generate valid high order meshes with optimized geometrical accuracy. The high order meshing procedure starts with a linear mesh, that is subsequently curved without taking care of the validity of the high order elements. An optimization procedure is then used to both untangle invalid elements and optimize the geometrical accuracy of the mesh. Standard measures of the distance between curves are considered to evaluate the geometrical accuracy in planar two-dimensional meshes, but they prove computationally too costly for optimization purposes. A fast estimate of the geometrical accuracy, based on Taylor expansions of the curves, is introduced. An unconstrained optimization procedure based on this estimate is shown to yield significant improvements in the geometrical accuracy of high order meshes, as measured by the standard Haudorff distance between the geometrical model and the mesh. Several examples illustrate the beneficial impact of this method on CFD solutions, with a part...
Adaptive mesh generation for image registration and segmentation
Fogtmann, Mads; Larsen, Rasmus
This paper deals with the problem of generating quality tetrahedral meshes for image registration. From an initial coarse mesh the approach matches the mesh to the image volume by combining red-green subdivision and mesh evolution through mesh-to-image matching regularized with a mesh quality...
Biologic mesh for abdominal wall reconstruction
King KS
2014-11-01
Full Text Available Kathryn S King,1 Frank P Albino,2 Parag Bhanot3 1School of Medicine, Georgetown University Hospital, Washington, DC, USA; 2Department of Plastic Surgery, 3Department of General Surgery, Georgetown University Hospital, Washington, DC, USA Background: Mesh reinforcement significantly decreases rates of recurrence following ventral hernia repair. Historically, biologic mesh was touted as superior in the setting of infection; however, selecting the appropriate mesh for a given clinical scenario is often a matter of debate. The purpose of this review is to highlight a number of the more commonly used biologic mesh products with a review of outcomes from the current literature. Methods: Outcomes following abdominal wall reconstruction using biologic mesh were reviewed for acellular cadaveric human dermis, cross-linked porcine dermis, non-cross-linked porcine dermis, porcine small intestine submucosa, acellular bovine pericardial, and acellular bovine dermal mesh. Studies with rigorous methods, adequate patient samples, and sufficient follow-up were selected for review. Results: Hernia recurrence rates following biologic mesh reinforcement vary widely. Porcine small intestine submucosa and bovine pericardium were associated with the lowest hernia recurrence rates. Porcine cross-linked dermal mesh products resulted in higher rates of adhesion formation and lower rates of tissue incorporation compared to non-cross-linked porcine mesh. Conclusion: Successful ventral hernia repair can be achieved with acceptable complications rates for each of the reviewed mesh products. Biologic meshes have an advantage over synthetic mesh in contaminated wounds but their use may not be cost-effective in all patient populations. Those with and/or at high risk for wound complications may also undergo repair with biologic mesh. Keywords: biologic mesh, ventral hernia repair, acellular dermal matrix
Automatic mesh adaptivity for CADIS and FW-CADIS neutronics modeling of difficult shielding problems
The CADIS and FW-CADIS hybrid Monte Carlo/deterministic techniques dramatically increase the efficiency of neutronics modeling, but their use in the accurate design analysis of very large and geometrically complex nuclear systems has been limited by the large number of processors and memory requirements for their preliminary deterministic calculations and final Monte Carlo calculation. Three mesh adaptivity algorithms were developed to reduce the memory requirements of CADIS and FW-CADIS without sacrificing their efficiency improvement. First, a macro-material approach enhances the fidelity of the deterministic models without changing the mesh. Second, a deterministic mesh refinement algorithm generates meshes that capture as much geometric detail as possible without exceeding a specified maximum number of mesh elements. Finally, a weight window coarsening algorithm de-couples the weight window mesh and energy bins from the mesh and energy group structure of the deterministic calculations in order to remove the memory constraint of the weight window map from the deterministic mesh resolution. The three algorithms were used to enhance an FW-CADIS calculation of the prompt dose rate throughout the ITER experimental facility. Using these algorithms resulted in a 23.3% increase in the number of mesh tally elements in which the dose rates were calculated in a 10-day Monte Carlo calculation and, additionally, increased the efficiency of the Monte Carlo simulation by a factor of at least 3.4. The three algorithms enabled this difficult calculation to be accurately solved using an FW-CADIS simulation on a regular computer cluster, obviating the need for a world-class super computer. (authors)
Low-diffusivity scalar transport using a WENO scheme and dual meshing
Kubrak, B; Herlina, H; Greve, F; Wissink, JG
2013-01-01
Interfacial mass transfer of low-diffusive substances in an unsteady flow environment is marked by a very thin boundary layer at the interface and other regions with steep concentration gradients. A numerical scheme capable of resolving accurately most details of this process is presented. In this scheme, the fifth-order accurate WENO method developed by [13] was implemented on a non-uniform staggered mesh to discretize the scalar convection while for the scalar diffusion a fourth-order accur...
Nonlinear Analyses of Adobe Masonry Walls Reinforced with Fiberglass Mesh
Vincenzo Giamundo
2014-02-01
Full Text Available Adobe constructions were widespread in the ancient world, and earth was one of the most used construction materials in ancient times. Therefore, the preservation of adobe structures, especially against seismic events, is nowadays an important structural issue. Previous experimental tests have shown that the ratio between mortar and brick mechanical properties (i.e., strength, stiffness and elastic modulus influences the global response of the walls in terms of strength and ductility. Accurate analyses are presented in both the case of unreinforced and reinforced with fiberglass mesh when varying the mechanical properties of the materials composing the adobe masonry structure. The main issues and variability in the behavior of seismic resisting walls when varying the mechanical properties are herein highlighted. The aim of the overall research activity is to improve the knowledge about the structural behavior of adobe structural members unreinforced and reinforced with fiberglass mesh inside horizontal mortar joints.
Anisotropic Diffusion in Mesh-Free Numerical Magnetohydrodynamics
Hopkins, Philip F
2016-01-01
We extend recently-developed mesh-free Lagrangian methods for numerical magnetohydrodynamics (MHD) to arbitrary anisotropic diffusion equations, including: passive scalar diffusion, Spitzer-Braginskii conduction and viscosity, cosmic ray diffusion/streaming, anisotropic radiation transport, non-ideal MHD (Ohmic resistivity, ambipolar diffusion, the Hall effect), and turbulent 'eddy diffusion.' We study these as implemented in the code GIZMO for both new meshless finite-volume Godunov schemes (MFM/MFV) as well as smoothed-particle hydrodynamics (SPH). We show the MFM/MFV methods are accurate and stable even with noisy fields and irregular particle arrangements, and recover the correct behavior even in arbitrarily anisotropic cases. They are competitive with state-of-the-art AMR/moving-mesh methods, and can correctly treat anisotropic diffusion-driven instabilities (e.g. the MTI and HBI, Hall MRI). We also develop a new scheme for stabilizing anisotropic tensor-valued fluxes with high-order gradient estimators ...
Jennings Jason
2010-01-01
Full Text Available Laparoscopic inguinal herniorraphy via a transabdominal preperitoneal (TAPP approach using Polypropylene Mesh (Mesh and staples is an accepted technique. Mesh induces a localised inflammatory response that may extend to, and involve, adjacent abdominal and pelvic viscera such as the appendix. We present an interesting case of suspected Mesh-induced appendicitis treated successfully with laparoscopic appendicectomy, without Mesh removal, in an elderly gentleman who presented with symptoms and signs of acute appendicitis 18 months after laparoscopic inguinal hernia repair. Possible mechanisms for Mesh-induced appendicitis are briefly discussed.
The moving mesh code Shadowfax
Vandenbroucke, Bert; De Rijcke, Sven
2016-01-01
We introduce the moving mesh code Shadowfax, which can be used to evolve a mixture of gas, subject to the laws of hydrodynamics and gravity, and any collisionless fluid only subject to gravity, such as cold dark matter or stars. The code is written in C++ and its source code is made available to the scientific community under the GNU Affero General Public License. We outline the algorithm and the design of our implementation, and demonstrate its validity through the results of a set of basic ...
Digital photography in anatomical pathology
Leong F; Leong A
2004-01-01
Digital imaging has made major inroads into the routine practice of anatomical pathology and replaces photographic prints and Kodachromes for reporting and conference purposes. More advanced systems coupled to computers allow greater versatility and speed of turnaround as well as lower costs of incorporating macroscopic and microscopic pictures into pathology reports and publications. Digital images allow transmission to remote sites via the Internet for consultation, quality assurance and ed...
Development of a pregnant woman phantom using polygonal mesh, for dosimetric evaluations
Due to the embryo/fetus radiosensitivity the accurate estimation of the absorbed dose distribution in the abdominal area is an additional problem caused by the exposure of pregnant women to ionizing radiation in medical applications. This paper reports the construction and insertion of a fetal representation in a female geometry by means of 3D modeling techniques. In order to characterize an ECM the Grupo de Dosimetria Numerica (GDN) is using, mainly, simulators emitting gamma sources and voxel phantoms coupled to a MC code. The phantoms are predominantly constructed from stacks of magnetic resonance images (MRI), computed tomography (CT) (obtained from scans of real patients) or from 3D modeling techniques. Due to the difficulty of obtaining medical images of pregnant women, 3D objects in several formats (.obj, .max, .blend, etc.) were acquired for anatomical representation of a non-pregnant adult. To construct a fetal representation, the 3D modeling technique called Poly Modeling (polygon mesh) was used inside of the software Autodesk 3ds Max 2014 (free student version). Information about the radiosensibility of organs included in the abdominal area will be used to fit and use the pregnant phantom in numerical dosimetry. For this, the phantom will be voxelized and the masses of organs of interest will be adjusted according to data provided by International Commission on Radiological Protection (ICRP). Finally, the phantom will be coupled to a MC code creating a MCE that will serve as base for the construction of several other models involving pregnant women submitted to ionizing radiation. (author)
Development of a pregnant woman phantom using polygonal mesh, for dosimetric evaluations
Cabral, Manuela O.M.; Vieira, Jose W., E-mail: manuela.omc@gmail.com [Universidade Federal de Pernambuco (DEN/UFPE), Recife, PE (Brazil). Departamento de Energia Nuclear; Leal Neto, Viriato, E-mail: viriatoleal@yahoo.com.br [Instituto Federal de Educacao, Ciencia e Tecnologia de Pernambuco (IFPE), Recife, PE (Brazil); Lima, Fernando R.A., E-mail: falima@cnen.gov.br [Centro Regional de Ciencias Nucleares do Nordeste (CRCN-NE/CNEN-PE), Recife, PE (Brazil)
2014-07-01
Due to the embryo/fetus radiosensitivity the accurate estimation of the absorbed dose distribution in the abdominal area is an additional problem caused by the exposure of pregnant women to ionizing radiation in medical applications. This paper reports the construction and insertion of a fetal representation in a female geometry by means of 3D modeling techniques. In order to characterize an ECM the Grupo de Dosimetria Numerica (GDN) is using, mainly, simulators emitting gamma sources and voxel phantoms coupled to a MC code. The phantoms are predominantly constructed from stacks of magnetic resonance images (MRI), computed tomography (CT) (obtained from scans of real patients) or from 3D modeling techniques. Due to the difficulty of obtaining medical images of pregnant women, 3D objects in several formats (.obj, .max, .blend, etc.) were acquired for anatomical representation of a non-pregnant adult. To construct a fetal representation, the 3D modeling technique called Poly Modeling (polygon mesh) was used inside of the software Autodesk 3ds Max 2014 (free student version). Information about the radiosensibility of organs included in the abdominal area will be used to fit and use the pregnant phantom in numerical dosimetry. For this, the phantom will be voxelized and the masses of organs of interest will be adjusted according to data provided by International Commission on Radiological Protection (ICRP). Finally, the phantom will be coupled to a MC code creating a MCE that will serve as base for the construction of several other models involving pregnant women submitted to ionizing radiation. (author)
Sierra toolkit computational mesh conceptual model
The Sierra Toolkit computational mesh is a software library intended to support massively parallel multi-physics computations on dynamically changing unstructured meshes. This domain of intended use is inherently complex due to distributed memory parallelism, parallel scalability, heterogeneity of physics, heterogeneous discretization of an unstructured mesh, and runtime adaptation of the mesh. Management of this inherent complexity begins with a conceptual analysis and modeling of this domain of intended use; i.e., development of a domain model. The Sierra Toolkit computational mesh software library is designed and implemented based upon this domain model. Software developers using, maintaining, or extending the Sierra Toolkit computational mesh library must be familiar with the concepts/domain model presented in this report.
Improved AFEM algorithm for bioluminescence tomography based on dual-mesh alternation strategy
Wei Li; Heng Zhao; Xiaochao Qu; Yanbin Hou; Xueli Chen; Duofang Chen; Xiaowei He; Qitan Zhang; Jimin Liang
2012-01-01
Adaptive finite element method (AFEM) is broadly adopted to recover the internal source in biological tissues.In this letter,a novel dual-mesh alternation strategy (dual-mesh AFEM) is developed for bioluminescence tomography.By comprehensively considering the error estimation of the finite element method solution on each mesh,two different adaptive strategies based on the error indicator of the reconstructed source and the photon flux density are used alternately in the process.Combined with the constantly adjusted permissible region in the adaptive process,the new algorithm can achieve a more accurate source location compared with the AFEM in the previous experiments.%Adaptive finite element method (AFEM) is broadly adopted to recover the internal source in biological tissues. In this letter, a novel dual-mesh alternation strategy (dual-mesh AFEM) is developed for biolumi-nescence tomography. By comprehensively considering the error estimation of the finite element method solution on each mesh, two different adaptive strategies based on the error indicator of the reconstructed source and the photon flux density are used alternately in the process. Combined with the constantly adjusted permissible region in the adaptive process, the new algorithm can achieve a more accurate source location compared with the AFEM in the previous experiments.
Towards Perceptual Quality Evaluation of Dynamic Meshes
Torkhani, Fakhri; Wang, Kai; Montanvert, Annick
2011-01-01
In practical applications, it is common that a 3D mesh undergoes some lossy operations. Since the end users of 3D meshes are often human beings, it is thus important to derive metrics that can faithfully assess the perceptual distortions induced by these operations. Like in the case of image quality assessment, metrics based on mesh geometric distances (e.g. Hausdorff distance and root mean squared error) cannot correctly predict the visual quality degradation. Recently, several perceptually-...
Unstructured Polyhedral Mesh Thermal Radiation Diffusion
Palmer, T.S.; Zika, M.R.; Madsen, N.K.
2000-07-27
Unstructured mesh particle transport and diffusion methods are gaining wider acceptance as mesh generation, scientific visualization and linear solvers improve. This paper describes an algorithm that is currently being used in the KULL code at Lawrence Livermore National Laboratory to solve the radiative transfer equations. The algorithm employs a point-centered diffusion discretization on arbitrary polyhedral meshes in 3D. We present the results of a few test problems to illustrate the capabilities of the radiation diffusion module.
Unstructured Polyhedral Mesh Thermal Radiation Diffusion
Unstructured mesh particle transport and diffusion methods are gaining wider acceptance as mesh generation, scientific visualization and linear solvers improve. This paper describes an algorithm that is currently being used in the KULL code at Lawrence Livermore National Laboratory to solve the radiative transfer equations. The algorithm employs a point-centered diffusion discretization on arbitrary polyhedral meshes in 3D. We present the results of a few test problems to illustrate the capabilities of the radiation diffusion module
Association Discovery Protocol for Hybrid Wireless Mesh Networks
Adjih, Cédric; Cho, Song Yean; Jacquet, Philippe
2006-01-01
Wireless mesh networks (WMNs) consist of two kinds of nodes: mesh routers which form the backbones of WMNs and mesh clients which associate with mesh routers to access networks. Because of the discrepancy between mesh routers and mesh clients, WMNs have a hybrid structure. Their hybrid structure presents an opportunity to integrate WMNs with different networks such as wireless LAN, Bluetooth and sensor networks through bridging functions in mesh routers. Because of the ability to integrate va...
Transport of phase space densities through tetrahedral meshes using discrete flow mapping
Bajars, Janis; Sondergaard, Niels; Tanner, Gregor
2016-01-01
Discrete flow mapping was recently introduced as an efficient ray based method determining wave energy distributions in complex built up structures. Wave energy densities are transported along ray trajectories through polygonal mesh elements using a finite dimensional approximation of a ray transfer operator. In this way the method can be viewed as a smoothed ray tracing method defined over meshed surfaces. Many applications require the resolution of wave energy distributions in three-dimensional domains, such as in room acoustics, underwater acoustics and for electromagnetic cavity problems. In this work we extend discrete flow mapping to three-dimensional domains by propagating wave energy densities through tetrahedral meshes. The geometric simplicity of the tetrahedral mesh elements is utilised to efficiently compute the ray transfer operator using a mixture of analytic and spectrally accurate numerical integration. The important issue of how to choose a suitable basis approximation in phase space whilst m...
Quadrilateral mesh fitting that preserves sharp features based on multi-normals for Laplacian energy
Yusuke Imai
2014-04-01
Full Text Available Because the cost of performance testing using actual products is expensive, manufacturers use lower-cost computer-aided design simulations for this function. In this paper, we propose using hexahedral meshes, which are more accurate than tetrahedral meshes, for finite element analysis. We propose automatic hexahedral mesh generation with sharp features to precisely represent the corresponding features of a target shape. Our hexahedral mesh is generated using a voxel-based algorithm. In our previous works, we fit the surface of the voxels to the target surface using Laplacian energy minimization. We used normal vectors in the fitting to preserve sharp features. However, this method could not represent concave sharp features precisely. In this proposal, we improve our previous Laplacian energy minimization by adding a term that depends on multi-normal vectors instead of using normal vectors. Furthermore, we accentuate a convex/concave surface subset to represent concave sharp features.
Converting skeletal structures to quad dominant meshes
Bærentzen, Jakob Andreas; Misztal, Marek Krzysztof; Welnicka, Katarzyna
2012-01-01
We propose the Skeleton to Quad-dominant polygonal Mesh algorithm (SQM), which converts skeletal structures to meshes composed entirely of polar and annular regions. Both types of regions have a regular structure where all faces are quads except for a single ring of triangles at the center of each...... polar region. The algorithm produces high quality meshes which contain irregular vertices only at the poles or where several regions join. It is trivial to produce a stripe parametrization for the output meshes which also lend themselves well to polar subdivision. After an initial description of SQM, we...
MOAB : a mesh-oriented database.
Tautges, Timothy James; Ernst, Corey; Stimpson, Clint; Meyers, Ray J.; Merkley, Karl
2004-04-01
A finite element mesh is used to decompose a continuous domain into a discretized representation. The finite element method solves PDEs on this mesh by modeling complex functions as a set of simple basis functions with coefficients at mesh vertices and prescribed continuity between elements. The mesh is one of the fundamental types of data linking the various tools in the FEA process (mesh generation, analysis, visualization, etc.). Thus, the representation of mesh data and operations on those data play a very important role in FEA-based simulations. MOAB is a component for representing and evaluating mesh data. MOAB can store structured and unstructured mesh, consisting of elements in the finite element 'zoo'. The functional interface to MOAB is simple yet powerful, allowing the representation of many types of metadata commonly found on the mesh. MOAB is optimized for efficiency in space and time, based on access to mesh in chunks rather than through individual entities, while also versatile enough to support individual entity access. The MOAB data model consists of a mesh interface instance, mesh entities (vertices and elements), sets, and tags. Entities are addressed through handles rather than pointers, to allow the underlying representation of an entity to change without changing the handle to that entity. Sets are arbitrary groupings of mesh entities and other sets. Sets also support parent/child relationships as a relation distinct from sets containing other sets. The directed-graph provided by set parent/child relationships is useful for modeling topological relations from a geometric model or other metadata. Tags are named data which can be assigned to the mesh as a whole, individual entities, or sets. Tags are a mechanism for attaching data to individual entities and sets are a mechanism for describing relations between entities; the combination of these two mechanisms is a powerful yet simple interface for representing metadata or application
Mesh Migration Into Urinary Bladder After Open Ventral Herniorrhaphy With Mesh: A Case Report
Su, Yann-Rong; Chan, Pei-Hui
2014-01-01
Ventral hernia repair with mesh products is of increasing popularity. The long-term results of mesh repair of ventral hernia are superior to primary suture repair. However, occasional complications may still present. We report on a 77-year-old man who underwent ventral hernia repair with a mesh 5 years ago with complication of mesh migration into the urinary bladder and enterovesical fistula. The patient presented with lower urinary tract symptoms initially. By urinalysis, persistent hematuri...
MeshLab: an open-source 3D mesh processing system
Cignoni, Paolo; Corsini, Massimiliano; Ranzuglia, Guido
2008-01-01
MeshLab is a free and open-source general-purpose mesh processing system designed to assist in the management of not-so-small, unstructured 3D models that typically occur in the pipeline when processing 3D scanned data in the context of Cultural Heritage. MeshLab provides a set of tools for editing, cleaning, healing, inspecting, rendering and converting the resulting meshes
Pattern recognition of anatomical shapes in CT scans
In medical image processing pattern recognition has become of major value in anatomical analysis and in computer aided information processing. Specifically, pattern recognition techniques simplify software development by means of which clinicians can manipulate anatomical relationships. As part of an overall CT pattern recognition system, a sequential edge tracking routine was devised together with a normalized Fourier descriptor analysis of identified shapes. A collection of shapes were extracted from CT scans of two patients and entered into an anatomic shape dictionary. This dictionary was employed in pattern matching experiments and in three-dimensional anatomical reconstruction. A sequential-edge tracking algorithm of high reliability, consistency, and image invariance, capable of utilizing heuristic and statistical rules, was demonstrated. Tests of pattern matching algorithms based on Fourier descriptors provided rapid and accurate body organ recognition of shapes extracted from de novo images using the shape dictionary. Results indicate that automated contour extraction and object recognition from cross-sectional images of human anatomy can be performed effectively, reliably, and rapidly. This abstract discusses an image processing environment that circumvents manual and subjective shape extraction, by substituting automatic and quantitative shape extraction, pattern matching and object recognition
Mesh Exposure and Associated Risk Factors in Women Undergoing Transvaginal Prolapse Repair with Mesh
Elizabeth A. Frankman
2013-01-01
Full Text Available Objective. To determine frequency, rate, and risk factors associated with mesh exposure in women undergoing transvaginal prolapse repair with polypropylene mesh. Methods. Retrospective chart review was performed for all women who underwent Prolift Pelvic Floor Repair System (Gynecare, Somerville, NJ between September 2005 and September 2008. Multivariable logistic regression was performed to identify risk factors for mesh exposure. Results. 201 women underwent Prolift. Mesh exposure occurred in 12% (24/201. Median time to mesh exposure was 62 days (range: 10–372. When mesh was placed in the anterior compartment, the frequency of mesh exposure was higher than that when mesh was placed in the posterior compartment (8.7% versus 2.9%, P=0.04. Independent risk factors for mesh exposure were diabetes (AOR = 7.7, 95% CI 1.6–37.6; P=0.01 and surgeon (AOR = 7.3, 95% CI 1.9–28.6; P=0.004. Conclusion. Women with diabetes have a 7-fold increased risk for mesh exposure after transvaginal prolapse repair using Prolift. The variable rate of mesh exposure amongst surgeons may be related to technique. The anterior vaginal wall may be at higher risk of mesh exposure as compared to the posterior vaginal wall.
Mesh Geometric Editing Approach Based on Gpu Texture
Guiping Qian
2012-09-01
Full Text Available This paper presents a novel interactive mesh editing approach based on GPU texture mapping. The main feature is that it copies 2D surface geometry information to GPU frame buffer. The planar mesh information is transformed into GPU texture and placed on apposite position of target mesh. 3D information is retrieved after stitching two mesh components from the primitive vertex coordinates. When running real-time mesh cloning operator, our mesh editing approach can copy arbitrary irregular geometric features from source mesh to target mesh. Experimental results indicate that our method can outperform previous related mesh editing techniques.
Adaptive mesh refinement in titanium
Colella, Phillip; Wen, Tong
2005-01-21
In this paper, we evaluate Titanium's usability as a high-level parallel programming language through a case study, where we implement a subset of Chombo's functionality in Titanium. Chombo is a software package applying the Adaptive Mesh Refinement methodology to numerical Partial Differential Equations at the production level. In Chombo, the library approach is used to parallel programming (C++ and Fortran, with MPI), whereas Titanium is a Java dialect designed for high-performance scientific computing. The performance of our implementation is studied and compared with that of Chombo in solving Poisson's equation based on two grid configurations from a real application. Also provided are the counts of lines of code from both sides.
CDCC calculations with the Lagrange-mesh technique
We apply the Lagrange-mesh technique to the Continuum Discretized Coupled Channel (CDCC) theory. The CDCC equations are solved with the R-matrix method, using Lagrange functions as variational basis. The choice of Lagrange functions is shown to be efficient and accurate for elastic scattering as well as for breakup reactions. We describe the general formalism for two-body projectiles, and apply it to the d+58Ni collision at Ed=80 MeV. Various numerical and physical aspects are discussed. Benchmark calculations on elastic scattering and breakup are presented.
Optimal Point Placement for Mesh Smoothing
Amenta, Nina; Bern, Marshall; Eppstein, David
1998-01-01
We study the problem of moving a vertex in an unstructured mesh of triangular, quadrilateral, or tetrahedral elements to optimize the shapes of adjacent elements. We show that many such problems can be solved in linear time using generalized linear programming. We also give efficient algorithms for some mesh smoothing problems that do not fit into the generalized linear programming paradigm.
7th International Meshing Roundtable '98
Eldred, T.J.
1998-10-01
The goal of the 7th International Meshing Roundtable is to bring together researchers and developers from industry, academia, and government labs in a stimulating, open environment for the exchange of technical information related to the meshing process. In the past, the Roundtable has enjoyed significant participation from each of these groups from a wide variety of countries.
Anatomically based lower limb nerve model for electrical stimulation
Soboleva Tanya K
2007-12-01
Full Text Available Abstract Background Functional Electrical Stimulation (FES is a technique that aims to rehabilitate or restore functionality of skeletal muscles using external electrical stimulation. Despite the success achieved within the field of FES, there are still a number of questions that remain unanswered. One way of providing input to the answers is through the use of computational models. Methods This paper describes the development of an anatomically based computer model of the motor neurons in the lower limb of the human leg and shows how it can be used to simulate electrical signal propagation from the beginning of the sciatic nerve to a skeletal muscle. One-dimensional cubic Hermite finite elements were used to represent the major portions of the lower limb nerves. These elements were fit to data that had been digitised using images from the Visible Man project. Nerves smaller than approximately 1 mm could not be seen in the images, and thus a tree-branching algorithm was used to connect the ends of the fitted nerve model to the respective skeletal muscle. To simulate electrical propagation, a previously published mammalian nerve model was implemented and solved on the anatomically based nerve mesh using a finite difference method. The grid points for the finite difference method were derived from the fitted finite element mesh. By adjusting the tree-branching algorithm, it is possible to represent different levels of motor-unit recruitment. Results To illustrate the process of a propagating nerve stimulus to a muscle in detail, the above method was applied to the nerve tree that connects to the human semitendinosus muscle. A conduction velocity of 89.8 m/s was obtained for a 15 μm diameter nerve fibre. This signal was successfully propagated down the motor neurons to a selected group of motor units in the muscle. Conclusion An anatomically and physiologically based model of the posterior motor neurons in the human lower limb was developed. This
Adaptive radial basis function mesh deformation using data reduction
Gillebaart, T.; Blom, D. S.; van Zuijlen, A. H.; Bijl, H.
2016-09-01
Radial Basis Function (RBF) mesh deformation is one of the most robust mesh deformation methods available. Using the greedy (data reduction) method in combination with an explicit boundary correction, results in an efficient method as shown in literature. However, to ensure the method remains robust, two issues are addressed: 1) how to ensure that the set of control points remains an accurate representation of the geometry in time and 2) how to use/automate the explicit boundary correction, while ensuring a high mesh quality. In this paper, we propose an adaptive RBF mesh deformation method, which ensures the set of control points always represents the geometry/displacement up to a certain (user-specified) criteria, by keeping track of the boundary error throughout the simulation and re-selecting when needed. Opposed to the unit displacement and prescribed displacement selection methods, the adaptive method is more robust, user-independent and efficient, for the cases considered. Secondly, the analysis of a single high aspect ratio cell is used to formulate an equation for the correction radius needed, depending on the characteristics of the correction function used, maximum aspect ratio, minimum first cell height and boundary error. Based on the analysis two new radial basis correction functions are derived and proposed. This proposed automated procedure is verified while varying the correction function, Reynolds number (and thus first cell height and aspect ratio) and boundary error. Finally, the parallel efficiency is studied for the two adaptive methods, unit displacement and prescribed displacement for both the CPU as well as the memory formulation with a 2D oscillating and translating airfoil with oscillating flap, a 3D flexible locally deforming tube and deforming wind turbine blade. Generally, the memory formulation requires less work (due to the large amount of work required for evaluating RBF's), but the parallel efficiency reduces due to the limited
Mesh deformation based on artificial neural networks
Stadler, Domen; Kosel, Franc; Čelič, Damjan; Lipej, Andrej
2011-09-01
In the article a new mesh deformation algorithm based on artificial neural networks is introduced. This method is a point-to-point method, meaning that it does not use connectivity information for calculation of the mesh deformation. Two already known point-to-point methods, based on interpolation techniques, are also presented. In contrast to the two known interpolation methods, the new method does not require a summation over all boundary nodes for one displacement calculation. The consequence of this fact is a shorter computational time of mesh deformation, which is proven by different deformation tests. The quality of the deformed meshes with all three deformation methods was also compared. Finally, the generated and the deformed three-dimensional meshes were used in the computational fluid dynamics numerical analysis of a Francis water turbine. A comparison of the analysis results was made to prove the applicability of the new method in every day computation.
Resolution of the Gross-Pitaevskii equation with the imaginary-time method on a Lagrange mesh.
Baye, D; Sparenberg, J-M
2010-11-01
The Lagrange-mesh method is an approximate variational calculation which has the simplicity of a mesh calculation. Combined with the imaginary-time method, it is applied to the iterative resolution of the Gross-Pitaevskii equation. Two variants of a fourth-order factorization of the exponential of the Hamiltonian and two types of mesh (Lagrange-Hermite and Lagrange-sinc) are employed and compared. The accuracy is checked with the help of these comparisons and of the virial theorem. The Lagrange-Hermite mesh provides very accurate results with short computing times for values of the dimensionless parameter of the nonlinear term up to 10⁴. For higher values up to 10⁷, the Lagrange-sinc mesh is more efficient. Examples are given for anisotropic and nonseparable trapping potentials. PMID:21230613
Digital imaging in anatomic pathology.
O'Brien, M J; Sotnikov, A V
1996-10-01
Advances in computer technology continue to bring new innovations to departments of anatomic pathology. This article briefly reviews the present status of digital optical imaging, and explores the directions that this technology may lead over the next several years. Technical requirements for digital microscopic and gross imaging, and the available options for image archival and retrieval are summarized. The advantages of digital images over conventional photography in the conference room, and the usefulness of digital imaging in the frozen section suite and gross room, as an adjunct to surgical signout and as a resource for training and education, are discussed. An approach to the future construction of digital histologic sections and the computer as microscope is described. The digital technologic applications that are now available as components of the surgical pathologist's workstation are enumerated. These include laboratory information systems, computerized voice recognition, and on-line or CD-based literature searching, texts and atlases and, in some departments, on-line image databases. The authors suggest that, in addition to these resources that are already available, tomorrow's surgical pathology workstation will include network-linked digital histologic databases, on-line software for image analysis and 3-D image enhancement, expert systems, and ultimately, advanced pattern recognition capabilities. In conclusion, the authors submit that digital optical imaging is likely to have a significant and positive impact on the future development of anatomic pathology. PMID:8853053
Characteristics of Mesh Wave Impedance in FDTD Non-Uniform Mesh
REN Wu; LIU Bo; GAO Ben-qing
2005-01-01
In order to increase the evaluating precision of mesh reflection wave, the mesh wave impedance(MWI) is extended to the non-uniform mesh in 1-D and 2-D cases for the first time on the basis of the Yee's positional relation for electromagnetic field components. Lots of characteristics are obtained for different mesh sizes and frequencies. Then the reflection coefficient caused by the non-uniform mesh can be calculated according to the theory of equivalent transmission line. By comparing it with that calculated by MWI in the uniform mesh, it is found that the evaluating error can be largely reduced and is in good agreement with that directly computed by FDTD method. And this extension of MWI can be used in the error analysis of complex mesh.
Update on Development of Mesh Generation Algorithms in MeshKit
Jain, Rajeev [Argonne National Lab. (ANL), Argonne, IL (United States); Vanderzee, Evan [Argonne National Lab. (ANL), Argonne, IL (United States); Mahadevan, Vijay [Argonne National Lab. (ANL), Argonne, IL (United States)
2015-09-30
MeshKit uses a graph-based design for coding all its meshing algorithms, which includes the Reactor Geometry (and mesh) Generation (RGG) algorithms. This report highlights the developmental updates of all the algorithms, results and future work. Parallel versions of algorithms, documentation and performance results are reported. RGG GUI design was updated to incorporate new features requested by the users; boundary layer generation and parallel RGG support were added to the GUI. Key contributions to the release, upgrade and maintenance of other SIGMA1 libraries (CGM and MOAB) were made. Several fundamental meshing algorithms for creating a robust parallel meshing pipeline in MeshKit are under development. Results and current status of automated, open-source and high quality nuclear reactor assembly mesh generation algorithms such as trimesher, quadmesher, interval matching and multi-sweeper are reported.
A moving mesh unstaggered constrained transport scheme for magnetohydrodynamics
Mocz, Philip; Pakmor, Rüdiger; Springel, Volker; Vogelsberger, Mark; Marinacci, Federico; Hernquist, Lars
2016-08-01
We present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an unstructured representation of the magnetic vector potential, making the numerical method simple and computationally efficient. The scheme is implemented in the moving mesh code AREPO. We demonstrate the performance of the approach with simulations of driven MHD turbulence, a magnetized disc galaxy, and a cosmological volume with primordial magnetic field. We compare the outcomes of these experiments to those obtained with a previously implemented Powell divergence-cleaning scheme. While CT and the Powell technique yield similar results in idealized test problems, some differences are seen in situations more representative of astrophysical flows. In the turbulence simulations, the Powell cleaning scheme artificially grows the mean magnetic field, while CT maintains this conserved quantity of ideal MHD. In the disc simulation, CT gives slower magnetic field growth rate and saturates to equipartition between the turbulent kinetic energy and magnetic energy, whereas Powell cleaning produces a dynamically dominant magnetic field. Such difference has been observed in adaptive-mesh refinement codes with CT and smoothed-particle hydrodynamics codes with divergence-cleaning. In the cosmological simulation, both approaches give similar magnetic amplification, but Powell exhibits more cell-level noise. CT methods in general are more accurate than divergence-cleaning techniques, and, when coupled to a moving mesh can exploit the advantages of automatic spatial/temporal adaptivity and reduced advection errors, allowing for improved astrophysical MHD simulations.
Constrained CVT Meshes and a Comparison of Triangular Mesh Generators
Nguyen, Hoa [Florida State University; Burkardt, John [Florida State University; Gunzburger, Max [Florida State University; Ju, Lili [University of South Carolina; Saka, Yuki [Florida State University
2009-01-01
Mesh generation in regions in Euclidean space is a central task in computational science, and especially for commonly used numerical methods for the solution of partial differential equations, e.g., finite element and finite volume methods. We focus on the uniform Delaunay triangulation of planar regions and, in particular, on how one selects the positions of the vertices of the triangulation. We discuss a recently developed method, based on the centroidal Voronoi tessellation (CVT) concept, for effecting such triangulations and present two algorithms, including one new one, for CVT-based grid generation. We also compare several methods, including CVT-based methods, for triangulating planar domains. To this end, we define several quantitative measures of the quality of uniform grids. We then generate triangulations of several planar regions, including some having complexities that are representative of what one may encounter in practice. We subject the resulting grids to visual and quantitative comparisons and conclude that all the methods considered produce high-quality uniform grids and that the CVT-based grids are at least as good as any of the others.
Method of modifying a volume mesh using sheet extraction
Borden, Michael J.; Shepherd, Jason F.
2007-02-20
A method and machine-readable medium provide a technique to modify a hexahedral finite element volume mesh using dual generation and sheet extraction. After generating a dual of a volume stack (mesh), a predetermined algorithm may be followed to modify the volume mesh of hexahedral elements. The predetermined algorithm may include the steps of determining a sheet of hexahedral mesh elements, generating nodes for merging, and merging the nodes to delete the sheet of hexahedral mesh elements and modify the volume mesh.
A high-precision simulation method for gas-liquid two-phase flows on unstructured meshes has been developed as a part of numerical studies on a gas entrainment phenomenon in the sodium-cooled fast reactor (JSFR). In this study, a two-dimensional unstructured adaptive mesh algorithm is developed because an adaptive mesh technique is necessary to simulate the local gas entrainment phenomenon accurately in large size JSFR. In a proposed two-dimensional adaptive mesh algorithm, each cell is isotropically subdivided to reduce distortions of the mesh. In addition, a connection cell is formed to eliminate the edge incompatibility between a refined and a non-refined cells. When forming connection cells, patterns of each connection cell is determined by subdivision condisions of neighboring cells. After checking the developed two-dimensional unstructured adaptive mesh manipulations (subdivision and merging of cells and construction of connection cells), the present adaptive mesh algorithm is verified by solving well-kwon driven cavity problem. As the result, the present unstructed adaptive mesh algorithm succeeds in reproducing vortical flow field in the cavity using relatively small cell number. (author)
Adaptive mesh refinement and automatic remeshing in crystal plasticity finite element simulations
In finite element simulations dedicated to the modelling of microstructure evolution, the mesh has to be fine enough to: (i) accurately describe the geometry of the constituents; (ii) capture local strain gradients stemming from the heterogeneity in material properties. In this paper, 3D polycrystalline aggregates are discretized into unstructured meshes and a level set framework is used to represent the grain boundaries. The crystal plasticity finite element method is used to simulate the plastic deformation of these aggregates. A mesh sensitivity analysis based on the deformation energy distribution shows that the predictions are, on average, more sensitive near grain boundaries. An anisotropic mesh refinement strategy based on the level set description is introduced and it is shown that it offers a good compromise between accuracy requirements on the one hand and computation time on the other hand. As the aggregates deform, mesh distortion inevitably occurs and ultimately causes the breakdown of the simulations. An automatic remeshing tool is used to periodically reconstruct the mesh and appropriate transfer of state variables is performed. It is shown that the diffusion related to data transfer is not significant. Finally, remeshing is performed repeatedly in a highly resolved 500 grains polycrystal subjected to about 90% thickness reduction in rolling. The predicted texture is compared with the experimental data and with the predictions of a standard Taylor model
Multi-Dimensional, Compressible Viscous Flow on a Moving Voronoi Mesh
Muñoz, Diego; Marcus, Robert; Vogelsberger, Mark; Hernquist, Lars
2012-01-01
Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow a significant reduction of advection errors and an accurate realization of curved and moving boundary conditions. Here we describe a novel formulation of viscous continuum hydrodynamics that solves the equations of motion on a Voronoi mesh created by a set of mesh-generating points. The points can move in an arbitrary manner, but the most natural motion is that given by the fluid velocity itself, such that the mesh dynamically adjusts to the flow. Owing to the mathematical properties of the Voronoi tessellation, pathological mesh-twisting effects are avoided. Our implementation considers the full Navier-Stokes equat...
Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes
Mavriplis, Dimitri J.; Thomas, James L. (Technical Monitor)
2003-01-01
The accuracy of the least-squares technique for gradient reconstruction on unstructured meshes is examined. While least-squares techniques produce accurate results on arbitrary isotropic unstructured meshes, serious difficulties exist for highly stretched meshes in the presence of surface curvature. In these situations, gradients are typically under-estimated by up to an order of magnitude. For vertex-based discretizations on triangular and quadrilateral meshes, and cell-centered discretizations on quadrilateral meshes, accuracy can be recovered using an inverse distance weighting in the least-squares construction. For cell-centered discretizations on triangles, both the unweighted and weighted least-squares constructions fail to provide suitable gradient estimates for highly stretched curved meshes. Good overall flow solution accuracy can be retained in spite of poor gradient estimates, due to the presence of flow alignment in exactly the same regions where the poor gradient accuracy is observed. However, the use of entropy fixes has the potential for generating large but subtle discretization errors.
Mesh Surgery for Anterior Vaginal Wall Prolapse: A Meta-analysis.
Juliato, Cássia Raquel Teatin; Santos Júnior, Luiz Carlos do; Haddad, Jorge Milhem; Castro, Rodrigo Aquino; Lima, Marcelo; Castro, Edilson Benedito de
2016-07-01
Purpose Pelvic organ prolapse (POP) is a major health issue worldwide, affecting 6-8% of women. The most affected site is the anterior vaginal wall. Multiple procedures and surgical techniques have been used, with or without the use of vaginal meshes, due to common treatment failure, reoperations, and complication rates in some studies. Methods Systematic review of the literature and meta-analysis regarding the use of vaginal mesh in anterior vaginal wall prolapse was performed. A total of 115 papers were retrieved after using the medical subject headings (MESH) terms: 'anterior pelvic organ prolapse OR cystocele AND surgery AND (mesh or colporrhaphy)' in the PubMed database. Exclusion criteria were: follow-up shorter than 1 year, use of biological or absorbable meshes, and inclusion of other vaginal wall prolapses. Studies were put in a data chart by two independent editors; results found in at least two studies were grouped for analysis. Results After the review of the titles by two independent editors, 70 studies were discarded, and after abstract assessment, 18 trials were eligible for full text screening. For final screening and meta-analysis, after applying the Jadad score (> 2), 12 studies were included. Objective cure was greater in the mesh surgery group (odds ratio [OR] = 1,28 [1,07-1,53]), which also had greater blood loss (mean deviation [MD] = 45,98 [9,72-82,25]), longer surgery time (MD = 15,08 [0,48-29,67]), but less prolapse recurrence (OR = 0,22 [01,3-0,38]). Dyspareunia, symptom resolution and reoperation rates were not statistically different between groups. Quality of life (QOL) assessment through the pelvic organ prolapse/urinary incontinence sexual questionnaire (PISQ-12), the pelvic floor distress inventory (PFDI-20), the pelvic floor impact questionnaire (PFIQ-7), and the perceived quality of life scale (PQOL) was not significantly different. Conclusions Anterior vaginal prolapse mesh surgery has greater anatomic
4D cone-beam CT reconstruction using multi-organ meshes for sliding motion modeling
Zhong, Zichun; Gu, Xuejun; Mao, Weihua; Wang, Jing
2016-02-01
A simultaneous motion estimation and image reconstruction (SMEIR) strategy was proposed for 4D cone-beam CT (4D-CBCT) reconstruction and showed excellent results in both phantom and lung cancer patient studies. In the original SMEIR algorithm, the deformation vector field (DVF) was defined on voxel grid and estimated by enforcing a global smoothness regularization term on the motion fields. The objective of this work is to improve the computation efficiency and motion estimation accuracy of SMEIR for 4D-CBCT through developing a multi-organ meshing model. Feature-based adaptive meshes were generated to reduce the number of unknowns in the DVF estimation and accurately capture the organ shapes and motion. Additionally, the discontinuity in the motion fields between different organs during respiration was explicitly considered in the multi-organ mesh model. This will help with the accurate visualization and motion estimation of the tumor on the organ boundaries in 4D-CBCT. To further improve the computational efficiency, a GPU-based parallel implementation was designed. The performance of the proposed algorithm was evaluated on a synthetic sliding motion phantom, a 4D NCAT phantom, and four lung cancer patients. The proposed multi-organ mesh based strategy outperformed the conventional Feldkamp-Davis-Kress, iterative total variation minimization, original SMEIR and single meshing method based on both qualitative and quantitative evaluations.
An adaptive mesh finite volume method for the Euler equations of gas dynamics
Mungkasi, Sudi
2016-06-01
The Euler equations have been used to model gas dynamics for decades. They consist of mathematical equations for the conservation of mass, momentum, and energy of the gas. For a large time value, the solution may contain discontinuities, even when the initial condition is smooth. A standard finite volume numerical method is not able to give accurate solutions to the Euler equations around discontinuities. Therefore we solve the Euler equations using an adaptive mesh finite volume method. In this paper, we present a new construction of the adaptive mesh finite volume method with an efficient computation of the refinement indicator. The adaptive method takes action automatically at around places having inaccurate solutions. Inaccurate solutions are reconstructed to reduce the error by refining the mesh locally up to a certain level. On the other hand, if the solution is already accurate, then the mesh is coarsened up to another certain level to minimize computational efforts. We implement the numerical entropy production as the mesh refinement indicator. As a test problem, we take the Sod shock tube problem. Numerical results show that the adaptive method is more promising than the standard one in solving the Euler equations of gas dynamics.
Anatomic consideration for preventive implantation.
Denissen, H W; Kalk, W; Veldhuis, H A; van Waas, M A
1993-01-01
The aim of preventive implant therapy is to prevent or delay loss of alveolar ridge bone mass. For use in an anatomic study of 60 mandibles, resorption of the alveolar ridge was classified into four preventive stages: (1) after extraction of teeth; (2) after initial resorption; (3) when the ridge has atrophied to a knife-edge shape; and (4) when only basal bone remains. Implantation in stage 3 necessitates removal of the knife-edge ridge to create space for cylindrical implants. Therefore, implantation in stage 2 is advocated to prevent the development of stage 3. The aim of implantation in stage 4 is to prevent total loss of function of the atrophic mandible. PMID:8359876
Polygonal Prism Mesh in the Viscous Layers for the Polyhedral Mesh Generator, PolyGen
Polyhedral mesh has been known to have some benefits over the tetrahedral mesh. Efforts have been made to set up a polyhedral mesh generation system with open source programs SALOME and TetGen. The evaluation has shown that the polyhedral mesh generation system is promising. But it is necessary to extend the capability of the system to handle the viscous layers to be a generalized mesh generator. A brief review to the previous works on the mesh generation for the viscous layers will be made in section 2. Several challenging issues for the polygonal prism mesh generation will be discussed as well. The procedure to generate a polygonal prism mesh will be discussed in detail in section 3. Conclusion will be followed in section 4. A procedure to generate meshes in the viscous layers with PolyGen has been successfully designed. But more efforts have to be exercised to find the best way for the generating meshes for viscous layers. Using the extrusion direction of the STL data will the first of the trials in the near future
Zhang, Fang
2011-02-01
Mesh current collectors made of stainless steel (SS) can be integrated into microbial fuel cell (MFC) cathodes constructed of a reactive carbon black and Pt catalyst mixture and a poly(dimethylsiloxane) (PDMS) diffusion layer. It is shown here that the mesh properties of these cathodes can significantly affect performance. Cathodes made from the coarsest mesh (30-mesh) achieved the highest maximum power of 1616 ± 25 mW m-2 (normalized to cathode projected surface area; 47.1 ± 0.7 W m-3 based on liquid volume), while the finest mesh (120-mesh) had the lowest power density (599 ± 57 mW m-2). Electrochemical impedance spectroscopy showed that charge transfer and diffusion resistances decreased with increasing mesh opening size. In MFC tests, the cathode performance was primarily limited by reaction kinetics, and not mass transfer. Oxygen permeability increased with mesh opening size, accounting for the decreased diffusion resistance. At higher current densities, diffusion became a limiting factor, especially for fine mesh with low oxygen transfer coefficients. These results demonstrate the critical nature of the mesh size used for constructing MFC cathodes. © 2010 Elsevier B.V. All rights reserved.
Engagement of Metal Debris into Gear Mesh
handschuh, Robert F.; Krantz, Timothy L.
2010-01-01
A series of bench-top experiments was conducted to determine the effects of metallic debris being dragged through meshing gear teeth. A test rig that is typically used to conduct contact fatigue experiments was used for these tests. Several sizes of drill material, shim stock and pieces of gear teeth were introduced and then driven through the meshing region. The level of torque required to drive the "chip" through the gear mesh was measured. From the data gathered, chip size sufficient to jam the mechanism can be determined.
Application of mesh network radios to UGS
Calcutt, Wade; Jones, Barry; Roeder, Brent
2008-04-01
During the past five years McQ has been actively pursuing integrating and applying wireless mesh network radios as a communications solution for unattended ground sensor (UGS) systems. This effort has been rewarded with limited levels of success and has ultimately resulted in a corporate position regarding the use of mesh network radios for UGS systems. A discussion into the background of the effort, the challenges of implementing commercial off-the-shelf (COTS) mesh radios with UGSs, the tradeoffs involved, and an overview of the future direction is presented.
Evaluation of Different Meshing Techniques for the Case of a Stented Artery.
Lotfi, Azadeh; Simmons, Anne; Barber, Tracie
2016-03-01
The formation and progression of in-stent restenosis (ISR) in bifurcated vessels may vary depending on the technique used for stenting. This study evaluates the effect of a variety of mesh styles on the accuracy and reliability of computational fluid dynamics (CFD) models in predicting these regions, using an idealized stented nonbifurcated model. The wall shear stress (WSS) and the near-stent recirculating vortices are used as determinants. The meshes comprise unstructured tetrahedral and polyhedral elements. The effects of local refinement, as well as higher-order elements such as prismatic inflation layers and internal hexahedral core, have also been examined. The uncertainty associated with individual mesh style was assessed through verification of calculations using the grid convergence index (GCI) method. The results obtained show that the only condition which allows the reliable comparison of uncertainty estimation between different meshing styles is that the monotonic convergence of grid solutions is in the asymptotic range. Comparisons show the superiority of a flow-adaptive polyhedral mesh over the commonly used adaptive and nonadaptive tetrahedral meshes in terms of resolving the near-stent flow features, GCI value, and prediction of WSS. More accurate estimation of hemodynamic factors was obtained using higher-order elements, such as hexahedral or prismatic grids. Incorporating these higher-order elements, however, was shown to introduce some degrees of numerical diffusion at the transitional area between the two meshes, not necessarily translating into high GCI value. Our data also confirmed the key role of local refinement in improving the performance and accuracy of nonadaptive mesh in predicting flow parameters in models of stented artery. The results of this study can provide a guideline for modeling biofluid domain in complex bifurcated arteries stented in regards to various stenting techniques. PMID:26784359
Numerical simulation of H2/air detonation using unstructured mesh
Togashi, Fumiya; Löhner, Rainald; Tsuboi, Nobuyuki
2009-06-01
To explore the capability of unstructured mesh to simulate detonation wave propagation phenomena, numerical simulation of H2/air detonation using unstructured mesh was conducted. The unstructured mesh has several adv- antages such as easy mesh adaptation and flexibility to the complicated configurations. To examine the resolution dependency of the unstructured mesh, several simulations varying the mesh size were conducted and compared with a computed result using a structured mesh. The results show that the unstructured mesh solution captures the detailed structure of detonation wave, as well as the structured mesh solution. To capture the detailed detonation cell structure, the unstructured mesh simulations required at least twice, ideally 5times the resolution of structured mesh solution.
Mesh Processing in Medical Image Analysis
The following topics are dealt with: mesh processing; medical image analysis; interactive freeform modeling; statistical shape analysis; clinical CT images; statistical surface recovery; automated segmentation; cerebral aneurysms; and real-time particle-based representation....
Shape space exploration of constrained meshes
Yang, Yongliang
2011-01-01
We present a general computational framework to locally characterize any shape space of meshes implicitly prescribed by a collection of non-linear constraints. We computationally access such manifolds, typically of high dimension and co-dimension, through first and second order approximants, namely tangent spaces and quadratically parameterized osculant surfaces. Exploration and navigation of desirable subspaces of the shape space with regard to application specific quality measures are enabled using approximants that are intrinsic to the underlying manifold and directly computable in the parameter space of the osculant surface. We demonstrate our framework on shape spaces of planar quad (PQ) meshes, where each mesh face is constrained to be (nearly) planar, and circular meshes, where each face has a circumcircle. We evaluate our framework for navigation and design exploration on a variety of inputs, while keeping context specific properties such as fairness, proximity to a reference surface, etc.
Shape space exploration of constrained meshes
Yang, Yongliang
2011-12-12
We present a general computational framework to locally characterize any shape space of meshes implicitly prescribed by a collection of non-linear constraints. We computationally access such manifolds, typically of high dimension and co-dimension, through first and second order approximants, namely tangent spaces and quadratically parameterized osculant surfaces. Exploration and navigation of desirable subspaces of the shape space with regard to application specific quality measures are enabled using approximants that are intrinsic to the underlying manifold and directly computable in the parameter space of the osculant surface. We demonstrate our framework on shape spaces of planar quad (PQ) meshes, where each mesh face is constrained to be (nearly) planar, and circular meshes, where each face has a circumcircle. We evaluate our framework for navigation and design exploration on a variety of inputs, while keeping context specific properties such as fairness, proximity to a reference surface, etc. © 2011 ACM.
LR: Compact connectivity representation for triangle meshes
Gurung, T; Luffel, M; Lindstrom, P; Rossignac, J
2011-01-28
We propose LR (Laced Ring) - a simple data structure for representing the connectivity of manifold triangle meshes. LR provides the option to store on average either 1.08 references per triangle or 26.2 bits per triangle. Its construction, from an input mesh that supports constant-time adjacency queries, has linear space and time complexity, and involves ordering most vertices along a nearly-Hamiltonian cycle. LR is best suited for applications that process meshes with fixed connectivity, as any changes to the connectivity require the data structure to be rebuilt. We provide an implementation of the set of standard random-access, constant-time operators for traversing a mesh, and show that LR often saves both space and traversal time over competing representations.
Obtuse triangle suppression in anisotropic meshes
Sun, Feng
2011-12-01
Anisotropic triangle meshes are used for efficient approximation of surfaces and flow data in finite element analysis, and in these applications it is desirable to have as few obtuse triangles as possible to reduce the discretization error. We present a variational approach to suppressing obtuse triangles in anisotropic meshes. Specifically, we introduce a hexagonal Minkowski metric, which is sensitive to triangle orientation, to give a new formulation of the centroidal Voronoi tessellation (CVT) method. Furthermore, we prove several relevant properties of the CVT method with the newly introduced metric. Experiments show that our algorithm produces anisotropic meshes with much fewer obtuse triangles than using existing methods while maintaining mesh anisotropy. © 2011 Elsevier B.V. All rights reserved.
Metal Mesh Filters for Terahertz Receivers Project
National Aeronautics and Space Administration — The technical objective of this SBIR program is to develop and demonstrate metal mesh filters for use in NASA's low noise receivers for terahertz astronomy and...
Bosmans, Lode; Valente, Giordano; Wesseling, Mariska; Van Campen, Anke; De Groote, Friedl; De Schutter, Joris; Jonkers, Ilse
2015-07-16
Scaled generic musculoskeletal models are commonly used to drive dynamic simulations of motions. It is however, acknowledged that not accounting for variability in musculoskeletal geometry and musculotendon parameters may confound the simulation results, even when analysing control subjects. This study documents the three-dimensional anatomical variability of musculotendon origins and insertions of 33 lower limb muscles determined based on magnetic resonance imaging in six subjects. This anatomical variability was compared to the musculotendon point location in a generic musculoskeletal model. Furthermore, the sensitivity of muscle forces during gait, calculated using static optimization, to perturbations of the musculotendon point location was analyzed with a generic model. More specific, a probabilistic approach was used: for each analyzed musculotendon point, the three-dimensional location was re-sampled with a uniform Latin hypercube method within the anatomical variability and the static optimization problem was then re-solved for all perturbations. We found that musculotendon point locations in the generic model showed only variable correspondences with the anatomical variability. The anatomical variability of musculotendon point location did affect the calculated muscle forces: muscles most sensitive to perturbations within the anatomical variability are iliacus and psoas. Perturbation of the gluteus medius anterior, iliacus and psoas induces the largest concomitant changes in muscle forces of the unperturbed muscles. Therefore, when creating subject-specific musculoskeletal models, these attachment points should be defined accurately. In addition, the size of the anatomical variability of the musculotendon point location was not related to the sensitivity of the calculated muscle forces. PMID:25979383
Automatic Assignment of Non-Leaf MeSH Terms to Biomedical Articles.
Kavuluru, Ramakanth; Rios, Anthony
2015-01-01
Assigning labels from a hierarchical vocabulary is a well known special case of multi-label classification, often modeled to maximize micro F1-score. However, building accurate binary classifiers for poorly performing labels in the hierarchy can improve both micro and macro F1-scores. In this paper, we propose and evaluate classification strategies involving descendant node instances to build better binary classifiers for non-leaf labels with the use-case of assigning Medical Subject Headings (MeSH) to biomedical articles. Librarians at the National Library of Medicine tag each biomedical article to be indexed by their PubMed information system with terms from the MeSH terminology, a biomedical conceptual hierarchy with over 27,000 terms. Human indexers look at each article's full text to assign a set of most suitable MeSH terms for indexing it. Several recent automated attempts focused on using the article title and abstract text to identify MeSH terms for the corresponding article. Despite these attempts, it is observed that assigning MeSH terms corresponding to certain non-leaf nodes of the MeSH hierarchy is particularly challenging. Non-leaf nodes are very important as they constitute one third of the total number of MeSH terms. Here, we demonstrate the effectiveness of exploiting training examples of descendant terms of non-leaf nodes in improving the performance of conventional classifiers for the corresponding non-leaf MeSH terms. Specifically, we focus on reducing the false positives (FPs) caused due to descendant instances in traditional classifiers. Our methods are able to achieve a relative improvement of 7.5% in macro-F1 score while also increasing the micro-F1 score by 1.6% for a set of 500 non-leaf terms in the MeSH hierarchy. These results strongly indicate the critical role of incorporating hierarchical information in MeSH term prediction. To our knowledge, our effort is the first to demonstrate the role of hierarchical information in improving
Doing the twist: diagonal meshes are isomorphic to twisted toroidal Meshes
Pearlmutter, Barak A
1996-01-01
We show that a k x n diagonal mesh is isomorphic to a n+k/2 x n+k/2 - nk/2 twisted toroidal mesh, i.e., a network similar to a standard n+k/2 x n-k/2 toroidal mesh, but with opposite handed twists of n-k/2 in the two directions, which results in a loss of (n-k/2)2 nodes.
Gamra: Simple Meshes for Complex Earthquakes
Landry, Walter; Barbot, Sylvain
2016-01-01
The static offsets caused by earthquakes are well described by elastostatic models with a discontinuity in the displacement along the fault. A traditional approach to model this discontinuity is to align the numerical mesh with the fault and solve the equations using finite elements. However, this distorted mesh can be difficult to generate and update. We present a new numerical method, inspired by the Immersed Interface Method, for solving the elastostatic equations with embedded discontinui...
MHD simulations on an unstructured mesh
Two reasons for using an unstructured computational mesh are adaptivity, and alignment with arbitrarily shaped boundaries. Two codes which use finite element discretization on an unstructured mesh are described. FEM3D solves 2D and 3D RMHD using an adaptive grid. MH3D++, which incorporates methods of FEM3D into the MH3D generalized MHD code, can be used with shaped boundaries, which might be 3D
Mesh Processing in Medical Image Analysis
The following topics are dealt with: mesh processing; medical image analysis; interactive freeform modeling; statistical shape analysis; clinical CT images; statistical surface recovery; automated segmentation; cerebral aneurysms; and real-time particle-based representation.......The following topics are dealt with: mesh processing; medical image analysis; interactive freeform modeling; statistical shape analysis; clinical CT images; statistical surface recovery; automated segmentation; cerebral aneurysms; and real-time particle-based representation....
Massively parallel computation on anisotropic meshes
Digonnet, Hugues; Silva, Luisa; Coupez, Thierry
2013-01-01
In this paper, we present developments done to obtain efficient parallel computations on supercomputers up to 8192 cores. While most massively parallel computation are shown using regular grid it is less common to see massively parallel computation using anisotropic adapted unstructured meshes. We will present here two mains components done to reach very large scale calculation up to 10 billions unknowns using a muligrid method over unstructured mesh running on 8192 cores. We firstly focus on...
Brain Morphometry Using Anatomical Magnetic Resonance Imaging
Bansal, Ravi; Gerber, Andrew J.; Peterson, Bradley S.
2008-01-01
The efficacy of anatomical magnetic resonance imaging (MRI) in studying the morphological features of various regions of the brain is described, also providing the steps used in the processing and studying of the images. The ability to correlate these features with several clinical and psychological measures can help in using anatomical MRI to…
Mesh geometry impact on Micromegas performance with an Exchangeable Mesh prototype
Kuger, F.; Bianco, M.; Iengo, P.; Sekhniaidze, G.; Veenhof, R.; Wotschack, J.
2016-07-01
The reconstruction precision of gaseous detectors is limited by losses of primary electrons during signal formation. In addition to common gas related losses, like attachment, Micromegas suffer from electron absorption during its transition through the micro mesh. This study aims for a deepened understanding of electron losses and their dependency on the mesh geometry. It combines experimental results obtained with a novel designed Exchangeable Mesh Micromegas (ExMe) and advanced microscopic-tracking simulations (ANSYS and Garfield++) of electron drift and mesh transition.
Unstructured Mesh Movement and Viscous Mesh Generation for CFD-Based Design Optimization Project
National Aeronautics and Space Administration — The innovations proposed are twofold: 1) a robust unstructured mesh movement method able to handle isotropic (Euler), anisotropic (viscous), mixed element (hybrid)...
How to model wireless mesh networks topology
The specification of network connectivity model or topology is the beginning of design and analysis in Computer Network researches. Wireless Mesh Networks is an autonomic network that is dynamically self-organised, self-configured while the mesh nodes establish automatic connectivity with the adjacent nodes in the relay network of wireless backbone routers. Researches in Wireless Mesh Networks range from node deployment to internetworking issues with sensor, Internet and cellular networks. These researches require modelling of relationships and interactions among nodes including technical characteristics of the links while satisfying the architectural requirements of the physical network. However, the existing topology generators model geographic topologies which constitute different architectures, thus may not be suitable in Wireless Mesh Networks scenarios. The existing methods of topology generation are explored, analysed and parameters for their characterisation are identified. Furthermore, an algorithm for the design of Wireless Mesh Networks topology based on square grid model is proposed in this paper. The performance of the topology generated is also evaluated. This research is particularly important in the generation of a close-to-real topology for ensuring relevance of design to the intended network and validity of results obtained in Wireless Mesh Networks researches
MHD simulations on an unstructured mesh
We describe work on a full MHD code using an unstructured mesh. MH3D++ is an extension of the PPPL MH3D resistive full MHD code. MH3D++ replaces the structured mesh and finite difference / fourier discretization of MH3D with an unstructured mesh and finite element / fourier discretization. Low level routines which perform differential operations, solution of PDEs such as Poisson's equation, and graphics, are encapsulated in C++ objects to isolate the finite element operations from the higher level code. The high level code is the same, whether it is run in structured or unstructured mesh versions. This allows the unstructured mesh version to be benchmarked against the structured mesh version. As a preliminary example, disruptions in DIIID reverse shear equilibria are studied numerically with the MH3D++ code. Numerical equilibria were first produced starting with an EQDSK file containing equilibrium data of a DIII-D L-mode negative central shear discharge. Using these equilibria, the linearized equations are time advanced to get the toroidal mode number n = 1 linear growth rate and eigenmode, which is resistively unstable. The equilibrium and linear mode are used to initialize 3D nonlinear runs. An example shows poloidal slices of 3D pressure surfaces: initially, on the left, and at an intermediate time, on the right
Hydrodynamic simulations on a moving Voronoi mesh
Springel, Volker
2011-01-01
At the heart of any method for computational fluid dynamics lies the question of how the simulated fluid should be discretized. Traditionally, a fixed Eulerian mesh is often employed for this purpose, which in modern schemes may also be adaptively refined during a calculation. Particle-based methods on the other hand discretize the mass instead of the volume, yielding an approximately Lagrangian approach. It is also possible to achieve Lagrangian behavior in mesh-based methods if the mesh is allowed to move with the flow. However, such approaches have often been fraught with substantial problems related to the development of irregularity in the mesh topology. Here we describe a novel scheme that eliminates these weaknesses. It is based on a moving unstructured mesh defined by the Voronoi tessellation of a set of discrete points. The mesh is used to solve the hyperbolic conservation laws of ideal hydrodynamics with a finite volume approach, based on a second-order Godunov scheme with an exact Riemann solver. A...
Discrete differential geometry: The nonplanar quadrilateral mesh
Twining, Carole J.; Marsland, Stephen
2012-06-01
We consider the problem of constructing a discrete differential geometry defined on nonplanar quadrilateral meshes. Physical models on discrete nonflat spaces are of inherent interest, as well as being used in applications such as computation for electromagnetism, fluid mechanics, and image analysis. However, the majority of analysis has focused on triangulated meshes. We consider two approaches: discretizing the tensor calculus, and a discrete mesh version of differential forms. While these two approaches are equivalent in the continuum, we show that this is not true in the discrete case. Nevertheless, we show that it is possible to construct mesh versions of the Levi-Civita connection (and hence the tensorial covariant derivative and the associated covariant exterior derivative), the torsion, and the curvature. We show how discrete analogs of the usual vector integral theorems are constructed in such a way that the appropriate conservation laws hold exactly on the mesh, rather than only as approximations to the continuum limit. We demonstrate the success of our method by constructing a mesh version of classical electromagnetism and discuss how our formalism could be used to deal with other physical models, such as fluids.
Accurate Finite Difference Algorithms
Goodrich, John W.
1996-01-01
Two families of finite difference algorithms for computational aeroacoustics are presented and compared. All of the algorithms are single step explicit methods, they have the same order of accuracy in both space and time, with examples up to eleventh order, and they have multidimensional extensions. One of the algorithm families has spectral like high resolution. Propagation with high order and high resolution algorithms can produce accurate results after O(10(exp 6)) periods of propagation with eight grid points per wavelength.
Recommendation for maximum allowable mesh size for plant combustion analyses with CFD codes
Highlights: ► Used mesh size has to be small enough to resolve all pressure waves relevant for the structural response analyses. ► Maximum allowable mesh size for a combustion pressure load calculation decreases with increasing relevant natural frequency of the structure. ► Maximum allowable mesh size for a combustion pressure load calculation increases with increasing of the speed of the sound in the gas mixture. ► Maximum allowable mesh size can be calculated from the developed analytical formula. - Abstract: The selection of the maximum allowable mesh size for a fluid dynamic calculation with Computational Fluid Dynamic (CFD) codes is essential for the reliability of the results assuming suitable physical and numerical models are used. Calculations with CFD codes are necessary for the assessment of the consequences of pressure loads on containment structures due to possible hydrogen combustion in nuclear power plants in a severe accident and on piping system due to pressure wave propagation in case of a pipe break accident or fast closing of a valve in a pipe with forced flow. CFD simulations of the transport and distribution of the released hydrogen/steam as well as the possible combustion during the transient in the containment require an appropriate mesh size to resolve the relevant phenomena and loads. The determination of the mesh size has to take into account: •adequate delineation of the containment geometry for accurate hydrogen distribution calculations, •sufficient conservative resolution of the combustion phenomena for the determination of pressure wave propagation and pressure loads, •no loss of pressure wave loads with relevant frequencies for the structural response analysis of the containment during the combustion calculation. In this paper, it is found that the accuracy of the calculated pressure wave associated with its frequency depends on the mesh size and a simple and easily useable analytical formula for the determination of
Recommendation for maximum allowable mesh size for plant combustion analyses with CFD codes
Movahed-Shariat-Panahi, M.A., E-mail: Mohammad-Ali.Movahed@areva.com [AREVA GmbH Offenbach (Germany)
2012-12-15
Highlights: Black-Right-Pointing-Pointer Used mesh size has to be small enough to resolve all pressure waves relevant for the structural response analyses. Black-Right-Pointing-Pointer Maximum allowable mesh size for a combustion pressure load calculation decreases with increasing relevant natural frequency of the structure. Black-Right-Pointing-Pointer Maximum allowable mesh size for a combustion pressure load calculation increases with increasing of the speed of the sound in the gas mixture. Black-Right-Pointing-Pointer Maximum allowable mesh size can be calculated from the developed analytical formula. - Abstract: The selection of the maximum allowable mesh size for a fluid dynamic calculation with Computational Fluid Dynamic (CFD) codes is essential for the reliability of the results assuming suitable physical and numerical models are used. Calculations with CFD codes are necessary for the assessment of the consequences of pressure loads on containment structures due to possible hydrogen combustion in nuclear power plants in a severe accident and on piping system due to pressure wave propagation in case of a pipe break accident or fast closing of a valve in a pipe with forced flow. CFD simulations of the transport and distribution of the released hydrogen/steam as well as the possible combustion during the transient in the containment require an appropriate mesh size to resolve the relevant phenomena and loads. The determination of the mesh size has to take into account: Bullet adequate delineation of the containment geometry for accurate hydrogen distribution calculations, Bullet sufficient conservative resolution of the combustion phenomena for the determination of pressure wave propagation and pressure loads, Bullet no loss of pressure wave loads with relevant frequencies for the structural response analysis of the containment during the combustion calculation. In this paper, it is found that the accuracy of the calculated pressure wave associated with its
The existing studies concerning image-free navigated implantation of hip resurfacing arthroplasty are based on analysis of the accuracy of conventional biplane radiography. Studies have shown that these measurements in biplane radiography are imprecise and that precision is improved by use of three-dimensional (3D) computer tomography (CT) scans. To date, the accuracy of image-free navigation devices for hip resurfacing has not been investigated using CT scans, and anteversion accuracy has not been assessed at all. Furthermore, no study has tested the reliability of the navigation software concerning the automatically calculated implant position. The purpose of our study was to analyze the accuracy of varus-valgus and anteversion using an image-free hip resurfacing navigation device. The reliability of the software-calculated implant position was also determined. A total of 32 femoral hip resurfacing components were implanted on embalmed human femurs using an image-free navigation device. In all, 16 prostheses were implanted with the proposed position generated by the navigation software; the 16 prostheses were inserted in an optimized valgus position. A 3D CT scan was undertaken before and after operation. The difference between the measured and planned varus-valgus angle averaged 1 deg (mean±standard deviation (SD): group I, 1 deg±2 deg; group II, 1 deg±1 deg). The mean±SD difference between femoral neck anteversion and anteversion of the implant was 4 deg (group I, 4 deg±4 deg; group II, 4 deg±3 deg). The software-calculated implant position differed 7 deg±8 deg from the measured neck-shaft angle. These measured accuracies did not differ significantly between the two groups. Our study proved the high accuracy of the navigation device concerning the most important biomechanical factor: the varus-valgus angle. The software calculation of the proposed implant position has been shown to be inaccurate and needs improvement. Hence, manual adjustment of the implant position in the software-planning step is frequently required. (author)
Relativistic Vlasov-Maxwell modelling using finite volumes and adaptive mesh refinement
Wettervik, Benjamin Svedung; Siminos, Evangelos; Fülöp, Tünde
2016-01-01
The dynamics of collisionless plasmas can be modelled by the Vlasov-Maxwell system of equations. An Eulerian approach is needed to accurately describe processes that are governed by high energy tails in the distribution function, but is of limited efficiency for high dimensional problems. The use of an adaptive mesh can reduce the scaling of the computational cost with the dimension of the problem. Here, we present a relativistic Eulerian Vlasov-Maxwell solver with block-structured adaptive mesh refinement in one spatial and one momentum dimension. The discretization of the Vlasov equation is based on a high-order finite volume method. A flux corrected transport algorithm is applied to limit spurious oscillations and ensure the physical character of the distribution function. We demonstrate a speed-up by a factor of five, because of the use of an adaptive mesh, in a typical scenario involving laser-plasma interaction in the self-induced transparency regime.
Detecting Translation Errors in CAD Surfaces and Preparing Geometries for Mesh Generation
Petersson, N Anders; Chand, K K
2001-08-27
The authors have developed tools for the efficient preparation of CAD geometries for mesh generation. Geometries are read from IGES files and then maintained in a boundary-representation consisting of a patchwork of trimmed and untrimmed surfaces. Gross errors in the geometry can be identified and removed automatically while a user interface is provided for manipulating the geometry (such as correcting invalid trimming curves or removing unwanted details). Modifying the geometry by adding or deleting surfaces and/or sectioning it by arbitrary planes (e.g. symmetry planes) is also supported. These tools are used for robust and accurate geometry models for initial mesh generation and will be applied to in situ mesh generation requirements of moving and adaptive grid simulations.
Randomized clinical trial of self-gripping mesh versus sutured mesh for Lichtenstein hernia repair
Jorgensen, L N; Sommer, T; Assaadzadeh, S;
2012-01-01
BACKGROUND: Many patients develop discomfort after open repair of a groin hernia. It was hypothesized that suture fixation of the mesh is a cause of these symptoms. METHODS: This patient- and assessor-blinded randomized multicentre clinical trial compared a self-gripping mesh (Parietene Progrip...
Ultrasound appearances after mesh implantation-evidence of mesh contraction or folding?
Švabík, K.; Martan, A.; Mašata, J.; Haddad El, R.; Hubka, P.; Pavlíková, Markéta
2011-01-01
Roč. 22, č. 5 (2011), s. 529-533. ISSN 0937-3462 Grant ostatní: GA MZd(CZ) NR9216 Institutional research plan: CEZ:AV0Z10300504 Keywords : prolift anterior * mesh shrinking * mesh retraction * vaginal ultrasound * vaginal surgery Subject RIV: FK - Gynaecology, Childbirth Impact factor: 1.832, year: 2011
Atlas-Based Automatic Generation of Subject-Specific Finite Element Tongue Meshes.
Bijar, Ahmad; Rohan, Pierre-Yves; Perrier, Pascal; Payan, Yohan
2016-01-01
Generation of subject-specific 3D finite element (FE) models requires the processing of numerous medical images in order to precisely extract geometrical information about subject-specific anatomy. This processing remains extremely challenging. To overcome this difficulty, we present an automatic atlas-based method that generates subject-specific FE meshes via a 3D registration guided by Magnetic Resonance images. The method extracts a 3D transformation by registering the atlas' volume image to the subject's one, and establishes a one-to-one correspondence between the two volumes. The 3D transformation field deforms the atlas' mesh to generate the subject-specific FE mesh. To preserve the quality of the subject-specific mesh, a diffeomorphic non-rigid registration based on B-spline free-form deformations is used, which guarantees a non-folding and one-to-one transformation. Two evaluations of the method are provided. First, a publicly available CT-database is used to assess the capability to accurately capture the complexity of each subject-specific Lung's geometry. Second, FE tongue meshes are generated for two healthy volunteers and two patients suffering from tongue cancer using MR images. It is shown that the method generates an appropriate representation of the subject-specific geometry while preserving the quality of the FE meshes for subsequent FE analysis. To demonstrate the importance of our method in a clinical context, a subject-specific mesh is used to simulate tongue's biomechanical response to the activation of an important tongue muscle, before and after cancer surgery. PMID:26577253
Cartesian anisotropic mesh adaptation for compressible flow
Simulating transient compressible flows involving shock waves presents challenges to the CFD practitioner in terms of the mesh quality required to resolve discontinuities and prevent smearing. This paper discusses a novel two-dimensional Cartesian anisotropic mesh adaptation technique implemented for compressible flow. This technique, developed for laminar flow by Ham, Lien and Strong, is efficient because it refines and coarsens cells using criteria that consider the solution in each of the cardinal directions separately. In this paper the method will be applied to compressible flow. The procedure shows promise in its ability to deliver good quality solutions while achieving computational savings. The convection scheme used is the Advective Upstream Splitting Method (Plus), and the refinement/ coarsening criteria are based on work done by Ham et al. Transient shock wave diffraction over a backward step and shock reflection over a forward step are considered as test cases because they demonstrate that the quality of the solution can be maintained as the mesh is refined and coarsened in time. The data structure is explained in relation to the computational mesh, and the object-oriented design and implementation of the code is presented. Refinement and coarsening algorithms are outlined. Computational savings over uniform and isotropic mesh approaches are shown to be significant. (author)