3D boundary recovery by constrained Delaunay tetrahedralization

Three-dimensional boundary recovery is a fundamental problem in mesh generation. In this paper, we propose a practical algorithm for solving this problem. Our algorithm is based on the construction of a {\it constrained Delaunay tetrahedralization} (CDT) for a set of constraints (segments and facets). The algorithm adds additional points (so-called Steiner points) on segments only. The Steiner points are chosen in such a way that the resulting subsegments are Delaunay and their lengths are not unnecessarily short. It is theoretically guaranteed that the facets can be recovered without using Steiner points. The complexity of this algorithm is analyzed. The proposed algorithm has been implemented. Its performance is reported through various application examples

Drift-diffusion models for innovative semiconductor devices and their numerical solution

We present charge transport models for novel semiconductor devices which may include ionic species as well as their thermodynamically consistent finite volume discretization

Full Envelope Control of Nonlinear Plants with Parameter Uncertainty by Fuzzy Controller Scheduling

A full envelope controller synthesis technique is developed for multiple-input single-output (MISO) nonlinear systems with structured parameter uncertainty. The technique maximizes the controller\u27s valid region of operation, while guaranteeing pre-specified transient performance. The resulting controller does not require on-line adaptation, estimation, prediction or model identification. Fuzzy Logic (FL) is used to smoothly schedule independently designed point controllers over the operational envelope and parameter space of the system\u27s model. These point controllers are synthesized using techniques chosen by the designer, thus allowing an unprecedented amount of design freedom. By using established control theory for the point controllers, the resulting nonlinear dynamic controller is able to handle the dynamics of complex systems which can not otherwise be addressed by Fuzzy Logic Control. An analytical solution for parameters describing the membership functions allows the optimization to yield the location of point designs: both quantifying the controller\u27s coverage, and eliminating the need of extensive hand tuning of these parameters. The net result is a decrease in the number of point designs required. Geometric primitives used in the solution all have multi-dimensional interpretations (convex hull, ellipsoid, Voronoi-Delaunay diagrams) which allow for scheduling on n-dimensions, including uncertainty due to nonlinearities and parameter variation. Since many multiple-input multiple-output (MIMO) controller design techniques are accomplished by solving several MISO problems, this work bridges the gap to full envelope control of MIMO nonlinear systems with parameter variation

A Unified Framework for Parallel Anisotropic Mesh Adaptation

Finite-element methods are a critical component of the design and analysis procedures of many (bio-)engineering applications. Mesh adaptation is one of the most crucial components since it discretizes the physics of the application at a relatively low cost to the solver. Highly scalable parallel mesh adaptation methods for High-Performance Computing (HPC) are essential to meet the ever-growing demand for higher fidelity simulations. Moreover, the continuous growth of the complexity of the HPC systems requires a systematic approach to exploit their full potential. Anisotropic mesh adaptation captures features of the solution at multiple scales while, minimizing the required number of elements. However, it also introduces new challenges on top of mesh generation. Also, the increased complexity of the targeted cases requires departing from traditional surface-constrained approaches to utilizing CAD (Computer-Aided Design) kernels. Alongside the functionality requirements, is the need of taking advantage of the ubiquitous multi-core machines. More importantly, the parallel implementation needs to handle the ever-increasing complexity of the mesh adaptation code. In this work, we develop a parallel mesh adaptation method that utilizes a metric-based approach for generating anisotropic meshes. Moreover, we enhance our method by interfacing with a CAD kernel, thus enabling its use on complex geometries. We evaluate our method both with fixed-resolution benchmarks and within a simulation pipeline, where the resolution of the discretization increases incrementally. With the Telescopic Approach for scalable mesh generation as a guide, we propose a parallel method at the node (multi-core) for mesh adaptation that is expected to scale up efficiently to the upcoming exascale machines. To facilitate an effective implementation, we introduce an abstract layer between the application and the runtime system that enables the use of task-based parallelism for concurrent mesh operations. Our evaluation indicates results comparable to state-of-the-art methods for fixed-resolution meshes both in terms of performance and quality. The integration with an adaptive pipeline offers promising results for the capability of the proposed method to function as part of an adaptive simulation. Moreover, our abstract tasking layer allows the separation of different aspects of the implementation without any impact on the functionality of the method

Delaunay triangulation in R3 on the GPU

Ph.DDOCTOR OF PHILOSOPH

Computational Multiscale Methods

Computational Multiscale Methods play an important role in many modern computer simulations in material sciences with different time scales and different scales in space. Besides various computational challenges, the meeting brought together various applications from many disciplines and scientists from various scientific communities

Resonant Inelastic X-ray Scattering Studies of Elementary Excitations

In the past decade, Resonant Inelastic X-ray Scattering (RIXS) has made remarkable progress as a spectroscopic technique. This is a direct result of the availability of high-brilliance synchrotron X-ray radiation sources and of advanced photon detection instrumentation. The technique's unique capability to probe elementary excitations in complex materials by measuring their energy-, momentum-, and polarization-dependence has brought RIXS to the forefront of experimental photon science. We review both the experimental and theoretical RIXS investigations of the past decade, focusing on those determining the low-energy charge, spin, orbital and lattice excitations of solids. We present the fundamentals of RIXS as an experimental method and then review the theoretical state of affairs, its recent developments and discuss the different (approximate) methods to compute the dynamical RIXS response. The last decade's body of experimental RIXS data and its interpretation is surveyed, with an emphasis on RIXS studies of correlated electron systems, especially transition metal compounds. Finally, we discuss the promise that RIXS holds for the near future, particularly in view of the advent of x-ray laser photon sources.Comment: Review, 67 pages, 44 figure