Smooth parametric surfaces and n-sided patches

The theory of 'geometric continuity' within the subject of CAGD is reviewed. In particular, we are concerned with how parametric surface patches for CAGD can be pieced together to form a smooth Ck surface. The theory is applied to the problem of filling an n-sided hole occurring within a smooth rectangular patch complex. A number of solutions to this problem are surveyed

PARALLEL √3-SUBDIVISION with ANIMATION in CONSIDERATION of GEOMETRIC COMPLEXITY

We look at the broader field of geometric subdivision and the emerging field of parallel computing for the purpose of creating higher visual fidelity at an efficient pace. Primarily, we present a parallel algorithm for √3-Subdivision. When considering animation, we find that it is possible to do subdivision by providing only one variable input, with the rest being considered static. This reduces the amount of data transfer required to continually update a subdividing mesh. We can support recursive subdivision by applying the technique in passes. As a basis for analysis, we look at performance in an OpenCL implementation that utilizes a local graphics processing unit (GPU) and a parallel CPU. By overcoming current hardware limitations, we present an environment where general GPU computation of √3-Subdivision can be practical

Center Vortex Model for the Infrared Sector of SU(3) Yang-Mills Theory - Vortex Free Energy

The vortex free energy is studied in the random vortex world-surface model of the infrared sector of SU(3) Yang-Mills theory. The free energy of a center vortex extending into two spatial directions, which is introduced into Yang-Mills configurations when acting with the 't Hooft loop operator, is verified to furnish an order parameter for the deconfinement phase transition. It is shown to exhibit a weak discontinuity at the critical temperature, corresponding to the weak first order character of the transition.Comment: 13 pages, 2 figures containing 3 eps file

Computerized Analysis of Magnetic Resonance Images to Study Cerebral Anatomy in Developing Neonates

The study of cerebral anatomy in developing neonates is of great importance for the understanding of brain development during the early period of life. This dissertation therefore focuses on three challenges in the modelling of cerebral anatomy in neonates during brain development. The methods that have been developed all use Magnetic Resonance Images (MRI) as source data. To facilitate study of vascular development in the neonatal period, a set of image analysis algorithms are developed to automatically extract and model cerebral vessel trees. The whole process consists of cerebral vessel tracking from automatically placed seed points, vessel tree generation, and vasculature registration and matching. These algorithms have been tested on clinical Time-of- Flight (TOF) MR angiographic datasets. To facilitate study of the neonatal cortex a complete cerebral cortex segmentation and reconstruction pipeline has been developed. Segmentation of the neonatal cortex is not effectively done by existing algorithms designed for the adult brain because the contrast between grey and white matter is reversed. This causes pixels containing tissue mixtures to be incorrectly labelled by conventional methods. The neonatal cortical segmentation method that has been developed is based on a novel expectation-maximization (EM) method with explicit correction for mislabelled partial volume voxels. Based on the resulting cortical segmentation, an implicit surface evolution technique is adopted for the reconstruction of the cortex in neonates. The performance of the method is investigated by performing a detailed landmark study. To facilitate study of cortical development, a cortical surface registration algorithm for aligning the cortical surface is developed. The method first inflates extracted cortical surfaces and then performs a non-rigid surface registration using free-form deformations (FFDs) to remove residual alignment. Validation experiments using data labelled by an expert observer demonstrate that the method can capture local changes and follow the growth of specific sulcus

Kirchhoff-Love shell representation and analysis using triangle configuration B-splines

This paper presents the application of triangle configuration B-splines (TCB-splines) for representing and analyzing the Kirchhoff-Love shell in the context of isogeometric analysis (IGA). The Kirchhoff-Love shell formulation requires global $C^1$-continuous basis functions. The nonuniform rational B-spline (NURBS)-based IGA has been extensively used for developing Kirchhoff-Love shell elements. However, shells with complex geometries inevitably need multiple patches and trimming techniques, where stitching patches with high continuity is a challenge. On the other hand, due to their unstructured nature, TCB-splines can accommodate general polygonal domains, have local refinement, and are flexible to model complex geometries with $C^1$ continuity, which naturally fit into the Kirchhoff-Love shell formulation with complex geometries. Therefore, we propose to use TCB-splines as basis functions for geometric representation and solution approximation. We apply our method to both linear and nonlinear benchmark shell problems, where the accuracy and robustness are validated. The applicability of the proposed approach to shell analysis is further exemplified by performing geometrically nonlinear Kirchhoff-Love shell simulations of a pipe junction and a front bumper represented by a single patch of TCB-splines

An immersed boundary hierarchical B-spline method for flexoelectricity

© 2019. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/This paper develops a computational framework with unfitted meshes to solve linear piezoelectricity and flexoelectricity electromechanical boundary value problems including strain gradient elasticity at infinitesimal strains. The high-order nature of the coupled PDE system is addressed by a sufficiently smooth hierarchical B-spline approximation on a background Cartesian mesh. The domain of interest is embedded into the background mesh and discretized in an unfitted fashion. The immersed boundary approach allows us to use B-splines on arbitrary domain shapes, regardless of their geometrical complexity, and could be directly extended, for instance, to shape and topology optimization. The domain boundary is represented by NURBS, and exactly integrated by means of the NEFEM mapping. Local adaptivity is achieved by hierarchical refinement of B-spline basis, which are efficiently evaluated and integrated thanks to their piecewise polynomial definition. Nitsche's formulation is derived to weakly enforce essential boundary conditions, accounting also for the non-local conditions on the non-smooth portions of the domain boundary (i.e. edges in 3D or corners in 2D) arising from Mindlin's strain gradient elasticity theory. Boundary conditions modeling sensing electrodes are formulated and enforced following the same approach. Optimal error convergence rates are reported using high-order B-spline approximations. The method is verified against available analytical solutions and well-known benchmarks from the literature.Peer ReviewedPostprint (author's final draft

Subdivision Surface based One-Piece Representation

Subdivision surfaces are capable of modeling and representing complex shapes of arbi-trary topology. However, methods on how to build the control mesh of a complex surfaceare not studied much. Currently, most meshes of complicated objects come from trian-gulation and simplification of raster scanned data points, like the Stanford 3D ScanningRepository. This approach is costly and leads to very dense meshes.Subdivision surface based one-piece representation means to represent the final objectin a design process with only one subdivision surface, no matter how complicated theobject\u27s topology or shape. Hence the number of parts in the final representation isalways one.In this dissertation we present necessary mathematical theories and geometric algo-rithms to support subdivision surface based one-piece representation. First, an explicitparametrization method is presented for exact evaluation of Catmull-Clark subdivisionsurfaces. Based on it, two approaches are proposed for constructing the one-piece rep-resentation of a given object with arbitrary topology. One approach is to construct theone-piece representation by using the interpolation technique. Interpolation is a naturalway to build models, but the fairness of the interpolating surface is a big concern inprevious methods. With similarity based interpolation technique, we can obtain bet-ter modeling results with less undesired artifacts and undulations. Another approachis through performing Boolean operations. Up to this point, accurate Boolean oper-ations over subdivision surfaces are not approached yet in the literature. We presenta robust and error controllable Boolean operation method which results in a one-piecerepresentation. Because one-piece representations resulting from the above two methodsare usually dense, error controllable simplification of one-piece representations is needed.Two methods are presented for this purpose: adaptive tessellation and multiresolutionanalysis. Both methods can significantly reduce the complexity of a one-piece represen-tation and while having accurate error estimation.A system that performs subdivision surface based one-piece representation was im-plemented and a lot of examples have been tested. All the examples show that our ap-proaches can obtain very good subdivision based one-piece representation results. Eventhough our methods are based on Catmull-Clark subdivision scheme, we believe they canbe adapted to other subdivision schemes as well with small modifications

In-vivo heterogeneous functional and residual strains in human aortic valve leaflets

Residual and physiological functional strains in soft tissues are known to play an important role in modulating organ stress distributions. Yet, no known comprehensive information on residual strains exist, or non-invasive techniques to quantify in-vivo deformations for the aortic valve (AV) leaflets. Herein we present a completely non-invasive approach for determining heterogeneous strains – both functional and residual – in semilunar valves and apply it to normal human AV leaflets. Transesophageal 3D echocardiographic (3DE) images of the AV were acquired from open-heart transplant patients, with each AV leaflet excised after heart explant and then imaged in a flattened configuration ex-vivo. Using an established spline parameterization of both 3DE segmentations and digitized ex-vivo images (Aggarwal et al., 2014), surface strains were calculated for deformation between the ex-vivo and three in-vivo configurations: fully open, just-coapted, and fully-loaded. Results indicated that leaflet area increased by an average of 20% from the ex-vivo to in-vivo open states, with a highly heterogeneous strain field. The increase in area from open to just-coapted state was the highest at an average of 25%, while that from just-coapted to fully-loaded remained almost unaltered. Going from the ex-vivo to in-vivo mid-systole configurations, the leaflet area near the basal attachment shrank slightly, whereas the free edge expanded by ~10%. This was accompanied by a 10° −20° shear along the circumferential-radial direction. Moreover, the principal stretches aligned approximately with the circumferential and radial directions for all cases, with the highest stretch being along the radial direction. Collectively, these results indicated that even though the AV did not support any measurable pressure gradient in the just-coapted state, the leaflets were significantly pre-strained with respect to the excised state. Furthermore, the collagen fibers of the leaflet were almost fully recruited in the just-coapted state, making the leaflet very stiff with marginal deformation under full pressure. Lastly, the deformation was always higher in the radial direction and lower along the circumferential one, the latter direction made stiffer by the preferential alignment of collagen fibers. These results provide significant insight into the distribution of residual strains and the in-vivo strains encountered during valve opening and closing in AV leaflets, and will form an important component of the tool that can evaluate valve׳s functional properties in a non-invasive manner

The effect of geometrical configurations on flows in idealised and realistic vascular geometries

This thesis reports the use of computational fluid dynamics (CFD) to investigate geometrical effects on flows in idealised non-branching double curved geometries (Study A) and in realistic distal anastomoses where geometries have been determined in vivo using magnetic resonance imaging (Study B). The purpose of this research is to improve understanding of the effects of geometrical configurations, especially curvature and non-planarity, on steady flow in idealised non-branching double curved geometries typical of arteries such as the aortic arch, the right coronary artery or the femoral arteries and on pulsatile flow in realistic distal anastomosis geometries. It is explained that the further knowledge gained from these idealised geometries can be useful to understand flows in anatomically correct geometries in order to optimise the design of end-to-side bypass graft vessels in clinical surgery. In the Study A, three-dimensional computations of steady flows in planar and non-planar double bends with 8 = 0.25 (curvature ratio) at Reynolds numbers of 125 and 500 were performed using the Navier-Stokes solver called Nektar that is based on spectral/hp element methods. The numerical haemodynamics analysis is presented in terms of the various mechanical factors which primarily involve axial velocity, transverse flows, vorticity, coherent vertical structure and wall shear stress. From these results, we can anticipate the wall shear stress distributions and secondary flow patterns in various double bend geometries with different non-planarity at low Reynolds numbers. Non-planarity plays a significant role on the wall shear stress distribution and the mixing properties of flow in the double bends, both of which are believed to be important factors for the patency of bypass grafts. In the Study B, Three sets of MRI data from patients undergoing tunnelled or superficial femoral bypass surgery were processed to give input velocity waveforms and geometries. From the latter and using the same Navier Stokes solver, detailed patient-specific pulsatile haemodynamics were calculated. Wall shear stress and velocity are influenced by the anatomical geometry, and wall shear stress distributions in pulsatile flow are compared with those in steady flow. The correlation is discussed between the haemodynamics and the remodelling of the vessel following bypass surgery determined in follow-up studies carried out a few months or a few years later. The results suggest that better designs should be possible for bypass grafts and indicate how inflow conditions can affect the flow field. The findings imply that proper anastomotical configurations might, induce haemodynamics environments that may prevent cardiovascular disorders or delay the progression of vascular disease.Open acces

Optimal Surface Fitting of Point Clouds Using Local Refinement : Application to GIS Data

This open access book provides insights into the novel Locally Refined B-spline (LR B-spline) surface format, which is suited for representing terrain and seabed data in a compact way. It provides an alternative to the well know raster and triangulated surface representations. An LR B-spline surface has an overall smooth behavior and allows the modeling of local details with only a limited growth in data volume. In regions where many data points belong to the same smooth area, LR B-splines allow a very lean representation of the shape by locally adapting the resolution of the spline space to the size and local shape variations of the region. The iterative method can be modified to improve the accuracy in particular domains of a point cloud. The use of statistical information criterion can help determining the optimal threshold, the number of iterations to perform as well as some parameters of the underlying mathematical functions (degree of the splines, parameter representation). The resulting surfaces are well suited for analysis and computing secondary information such as contour curves and minimum and maximum points. Also deformation analysis are potential applications of fitting point clouds with LR B-splines