Extensions to OpenGL for CAGD.

Many computer graphic API’s, including OpenGL, emphasize modeling with rectangular patches, which are especially useful in Computer Aided Geomeric Design (CAGD). However, not all shapes are rectangular; some are triangular or more complex. This paper extends the OpenGL library to support the modeling of triangular patches, Coons patches, and Box-splines patches. Compared with the triangular patch created from degenerate rectangular Bezier patch with the existing functions provided by OpenGL, the triangular Bezier patches can be used in certain design situations and allow designers to achieve high-quality results that are less CPU intense and require less storage space. The addition of Coons patches and Box splines to the OpenGL library also give it more functionality. Both patch types give CAGD users more flexibility in designing surfaces. A library for all three patch types was developed as an addition to OpenGL

Efficient and High-Quality Rendering of Higher-Order Geometric Data Representations

Computer-Aided Design (CAD) bezeichnet den Entwurf industrieller Produkte mit Hilfe von virtuellen 3D Modellen. Ein CAD-Modell besteht aus parametrischen Kurven und Flächen, in den meisten Fällen non-uniform rational B-Splines (NURBS). Diese mathematische Beschreibung wird ebenfalls zur Analyse, Optimierung und Präsentation des Modells verwendet. In jeder dieser Entwicklungsphasen wird eine unterschiedliche visuelle Darstellung benötigt, um den entsprechenden Nutzern ein geeignetes Feedback zu geben. Designer bevorzugen beispielsweise illustrative oder realistische Darstellungen, Ingenieure benötigen eine verständliche Visualisierung der Simulationsergebnisse, während eine immersive 3D Darstellung bei einer Benutzbarkeitsanalyse oder der Designauswahl hilfreich sein kann. Die interaktive Darstellung von NURBS-Modellen und -Simulationsdaten ist jedoch aufgrund des hohen Rechenaufwandes und der eingeschränkten Hardwareunterstützung eine große Herausforderung. Diese Arbeit stellt vier neuartige Verfahren vor, welche sich mit der interaktiven Darstellung von NURBS-Modellen und Simulationensdaten befassen. Die vorgestellten Algorithmen nutzen neue Fähigkeiten aktueller Grafikkarten aus, um den Stand der Technik bezüglich Qualität, Effizienz und Darstellungsgeschwindigkeit zu verbessern. Zwei dieser Verfahren befassen sich mit der direkten Darstellung der parametrischen Beschreibung ohne Approximationen oder zeitaufwändige Vorberechnungen. Die dabei vorgestellten Datenstrukturen und Algorithmen ermöglichen die effiziente Unterteilung, Klassifizierung, Tessellierung und Darstellung getrimmter NURBS-Flächen und einen interaktiven Ray-Casting-Algorithmus für die Isoflächenvisualisierung von NURBSbasierten isogeometrischen Analysen. Die weiteren zwei Verfahren beschreiben zum einen das vielseitige Konzept der programmierbaren Transparenz für illustrative und verständliche Visualisierungen tiefenkomplexer CAD-Modelle und zum anderen eine neue hybride Methode zur Reprojektion halbtransparenter und undurchsichtiger Bildinformation für die Beschleunigung der Erzeugung von stereoskopischen Bildpaaren. Die beiden letztgenannten Ansätze basieren auf rasterisierter Geometrie und sind somit ebenfalls für normale Dreiecksmodelle anwendbar, wodurch die Arbeiten auch einen wichtigen Beitrag in den Bereichen der Computergrafik und der virtuellen Realität darstellen. Die Auswertung der Arbeit wurde mit großen, realen NURBS-Datensätzen durchgeführt. Die Resultate zeigen, dass die direkte Darstellung auf Grundlage der parametrischen Beschreibung mit interaktiven Bildwiederholraten und in subpixelgenauer Qualität möglich ist. Die Einführung programmierbarer Transparenz ermöglicht zudem die Umsetzung kollaborativer 3D Interaktionstechniken für die Exploration der Modelle in virtuellenUmgebungen sowie illustrative und verständliche Visualisierungen tiefenkomplexer CAD-Modelle. Die Erzeugung stereoskopischer Bildpaare für die interaktive Visualisierung auf 3D Displays konnte beschleunigt werden. Diese messbare Verbesserung wurde zudem im Rahmen einer Nutzerstudie als wahrnehmbar und vorteilhaft befunden.In computer-aided design (CAD), industrial products are designed using a virtual 3D model. A CAD model typically consists of curves and surfaces in a parametric representation, in most cases, non-uniform rational B-splines (NURBS). The same representation is also used for the analysis, optimization and presentation of the model. In each phase of this process, different visualizations are required to provide an appropriate user feedback. Designers work with illustrative and realistic renderings, engineers need a comprehensible visualization of the simulation results, and usability studies or product presentations benefit from using a 3D display. However, the interactive visualization of NURBS models and corresponding physical simulations is a challenging task because of the computational complexity and the limited graphics hardware support. This thesis proposes four novel rendering approaches that improve the interactive visualization of CAD models and their analysis. The presented algorithms exploit latest graphics hardware capabilities to advance the state-of-the-art in terms of quality, efficiency and performance. In particular, two approaches describe the direct rendering of the parametric representation without precomputed approximations and timeconsuming pre-processing steps. New data structures and algorithms are presented for the efficient partition, classification, tessellation, and rendering of trimmed NURBS surfaces as well as the first direct isosurface ray-casting approach for NURBS-based isogeometric analysis. The other two approaches introduce the versatile concept of programmable order-independent semi-transparency for the illustrative and comprehensible visualization of depth-complex CAD models, and a novel method for the hybrid reprojection of opaque and semi-transparent image information to accelerate stereoscopic rendering. Both approaches are also applicable to standard polygonal geometry which contributes to the computer graphics and virtual reality research communities. The evaluation is based on real-world NURBS-based models and simulation data. The results show that rendering can be performed directly on the underlying parametric representation with interactive frame rates and subpixel-precise image results. The computational costs of additional visualization effects, such as semi-transparency and stereoscopic rendering, are reduced to maintain interactive frame rates. The benefit of this performance gain was confirmed by quantitative measurements and a pilot user study

The Construction of Optimized High-Order Surface Meshes by Energy-Minimization

Despite the increasing popularity of high-order methods in computational fluid dynamics, their application to practical problems still remains challenging. In order to exploit the advantages of high-order methods with geometrically complex computational domains, coarse curved meshes are necessary, i.e. high-order representations of the geometry. This dissertation presents a strategy for the generation of curved high-order surface meshes. The mesh generation method combines least-squares fitting with energy functionals, which approximate physical bending and stretching energies, in an incremental energy-minimizing fitting strategy. Since the energy weighting is reduced in each increment, the resulting surface representation features high accuracy. Nevertheless, the beneficial influence of the energy-minimization is retained. The presented method aims at enabling the utilization of the superior convergence properties of high-order methods by facilitating the construction of coarser meshes, while ensuring accuracy by allowing an arbitrary choice of geometric approximation order. Results show surface meshes of remarkable quality, even for very coarse meshes representing complex domains, e.g. blood vessels

Software for Modeling Ultrasound Breast Cancer Imaging

Computer-based models are increasingly used in biomedical imaging research to clarify links between anatomical structure, imaging physics, and the information content of medical images. A few three-dimensional breast tissue software models have been developed for mammography simulations to optimize current mammography systems or to test novel systems. It would be beneficial in the development of ultrasound breast imaging to have a similar computational model for simulation. A three-dimensional breast anatomy model with the lobular ducts, periductal and intralobular loose fibrous tissue, interlobular dense fibrous tissue, fat, and skin has been implemented. The parenchymal density of the model can be varied from about 20 to 75% to represent a range of clinically relevant densities. The anatomical model was used as a foundation for a three-dimensional breast tumour model. The tumour model was designed to mimic the ultrasound appearance of features used in tumour classification. Simulated two-dimensional ultrasound images were synthesized from the models using a first-order k-space propagation simulator. Similar to clinical ultrasound images, the simulated images of normal breast tissue exhibited non-Rayleigh speckle in regions of interest consisting of primarily fatty, primarily fibroglandular, and mixed tissue types. The simulated images of tumours reproduced several shape and margin features used in breast tumour diagnosis. The ultrasound wavefront distortion produced in simulations using the anatomical model was evaluated and a second method of modeling wavefront distortion was also proposed in which 10 to 12 irregularly shaped, strongly scattering inclusions were iii superimposed on multiple parallel time-shift screens to create the screen-inclusion model. Simulations of planar pulsed wave propagation through the two proposed models, a conventional parallel time-shift screen model, and digitized breast tissue specimens were compared. The anatomical model and screen-inclusion model were able to produce arrival-time fluctuation and energy-level fluctuation characteristics comparable to the digitized tissue specimens that the parallel-screen model was unable to reproduce. This software is expected to be valuable for imaging simulations that require accurate and detailed representation of the ultrasound characteristics of breast tumours

Compact and Intuitive Airfoil Parameterization Method through Physics-aware Variational Autoencoder

Airfoil shape optimization plays a critical role in the design of high-performance aircraft. However, the high-dimensional nature of airfoil representation causes the challenging problem known as the "curse of dimensionality". To overcome this problem, numerous airfoil parameterization methods have been developed, which can be broadly classified as polynomial-based and data-driven approaches. Each of these methods has desirable characteristics such as flexibility, parsimony, feasibility, and intuitiveness, but a single approach that encompasses all of these attributes has yet to be found. For example, polynomial-based methods struggle to balance parsimony and flexibility, while data-driven methods lack in feasibility and intuitiveness. In recent years, generative models, such as generative adversarial networks and variational autoencoders, have shown promising potential in airfoil parameterization. However, these models still face challenges related to intuitiveness due to their black-box nature. To address this issue, we developed a novel airfoil parameterization method using physics-aware variational autoencoder. The proposed method not only explicitly separates the generation of thickness and camber distributions to produce smooth and non-intersecting airfoils, thereby improving feasibility, but it also directly aligns its latent dimensions with geometric features of the airfoil, significantly enhancing intuitiveness. Finally, extensive comparative studies were performed to demonstrate the effectiveness of our approach.Comment: 33 pages, 19 figures, preprint for journal submissio

Shape analysis of the corpus callosum of autistic and normal subjects in neuroimaging.

Early detection of human disease in today’s society can have an enormous impact on the severity of the disease that is manifested. Disease such as Autism and Dyslexia, which have no current cure or proven mechanism as to how they develop, can often have an adverse physical and physiological impact on the lifestyle of a human being. Although these disease are not fully curable, the severity handicaps that accompany them can be significantly reduced with the proper therapy, and thus the earlier that the disease is detected the faster therapy can be administered. The research in this thesis is an attempt at studying discriminatory shape measures of some brain structures that are known to carry changes from autistics to normal individuals. The focus will be on the corpus callosum. There has been considerable research done on the brain scans (MRI, CT) of autistic individuals vs. control (normal) individuals to observe any noticeable discrepancies through statistical analysis. The most common and powerful tool to analyze structures of the brain, once a specific region has been segmented, is using Registration to match like structures and record their error. The ICP algorithm (Iterative Closest Point) is commonly used to accomplish this task. Many techniques such as level sets and statistical methods can be used for segmentation. The Corpus Callosum (CC) and the cortical surface of the brain are currently where most Autism analysis is performed. It has been observed that the gyrification of the cortical surface is different in the two groups, and size as well as shape of the CC. An analysis approach for autism MRI is quite extensive and involves many steps. This thesis is limited to examination of shape measures of the CC that lend discrimination ability to distinguish between normal and autistic individuals from T1-weigheted MRI scans. We will examine two approaches for shape analysis, based on the traditional Fourier Descriptors (FD) method and shape registration (SR) using the procrustes technique. MRI scans of 22 autistic and 16 normal individuals are used to test the approaches developed in this thesis. We show that both FD and SR may be used to extract features to discriminate between the two populations with accuracy levels over 80% up to 100% depending on the technique

Non-Uniform Rational B-Splines and Rational Bezier Triangles for Isogeometric Analysis of Structural Applications

Isogeometric Analysis (IGA) is a major advancement in computational analysis that bridges the gap between a computer-aided design (CAD) model, which is typically constructed using Non-Uniform Rational B-splines (NURBS), and a computational model that traditionally uses Lagrange polynomials to represent the geometry and solution variables. In IGA, the same shape functions that are used in CAD are employed for analysis. The direct manipulation of CAD data eliminates approximation errors that emanate from the process of converting the geometry from CAD to Finite Element Analysis (FEA). As a result, IGA allows the exact geometry to be represented at the coarsest level and maintained throughout the analysis process. While IGA was initially introduced to streamline the design and analysis process, this dissertation shows that IGA can also provide improved computational results for complex and highly nonlinear problems in structural mechanics. This dissertation addresses various problems in structural mechanics in the context of IGA, with the use of NURBS and rational Bézier triangles for the description of the parametric and physical spaces. The approaches considered here show that a number of important properties (e.g., high-order smoothness, geometric exactness, reduced number of degrees of freedom, and increased flexibility in discretization) can be achieved, leading to improved numerical solutions. Specifically, using B-splines and a layer-based discretization, a distributed plasticity isogeometric frame model is formulated to capture the spread of plasticity in large-deformation frames. The modeling approach includes an adaptive analysis where the structure of interest is initially modeled with coarse mesh and knots are inserted based on the yielding information at the quadrature points. It is demonstrated that improvement on efficiency and convergence rates is attained. With NURBS, an isogeometric rotation-free multi-layered plate formulation is developed based on a layerwise deformation theory. The derivation assumes a separate displacement field expansion within each layer, and considers transverse displacement component as C0-continuous at dissimilar material interfaces, which is enforced via knot repetition. The separate integration of the in-plane and through-thickness directions allows to capture the complete 3D stresses in a 2D setting. The proposed method is used to predict the behavior of advanced materials such as laminated composites, and the results show advantages in efficiency and accuracy. To increase the flexibility in discretizing complex geometries, rational Bézier triangles for domain triangulation is studied. They are further coupled with a Delaunay-based feature-preserving discretization algorithm for static bending and free vibration analysis of Kirchhoff plates. Lagrange multipliers are employed to explicitly impose high-order continuity constraints and the augmented system is solved iteratively without increasing the matrix size. The resulting discretization is geometrically exact, admits small geometric features, and constitutes C1-continuity. The feature-preserving rational Bézier triangles are further applied to smeared damage modeling of quasi-brittle materials. Due to the ability of Lagrange multipliers to raise global continuity to any desired order, the implicit fourth- and sixth-order gradient damage models are analyzed. The inclusion of higher-order terms in the nonlocal Taylor expansion improves solution accuracy. A local refinement algorithm that resolves marked regions with high resolution while keeping the resulting mesh conforming and well-conditioned is also utilized to improve efficiency. The outcome is a unified modeling framework where the feature-preserving discretization is able to capture the damage initiation and early-stage propagation, and the local refinement technique can then be applied to adaptively refine the mesh in the direction of damage propagation.PHDCivil EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/147668/1/ningliu_1.pd