Discontinuities in numerical radiative transfer

Observations and magnetohydrodynamic simulations of solar and stellar atmospheres reveal an intermittent behavior or steep gradients in physical parameters, such as magnetic field, temperature, and bulk velocities. The numerical solution of the stationary radiative transfer equation is particularly challenging in such situations, because standard numerical methods may perform very inefficiently in the absence of local smoothness. However, a rigorous investigation of the numerical treatment of the radiative transfer equation in discontinuous media is still lacking. The aim of this work is to expose the limitations of standard convergence analyses for this problem and to identify the relevant issues. Moreover, specific numerical tests are performed. These show that discontinuities in the atmospheric physical parameters effectively induce first-order discontinuities in the radiative transfer equation, reducing the accuracy of the solution and thwarting high-order convergence. In addition, a survey of the existing numerical schemes for discontinuous ordinary differential systems and interpolation techniques for discontinuous discrete data is given, evaluating their applicability to the radiative transfer problem

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

A novel parallel algorithm for surface editing and its FPGA implementation

A thesis submitted to the University of Bedfordshire in partial fulfilment of the requirements for the degree of Doctor of PhilosophySurface modelling and editing is one of important subjects in computer graphics. Decades of research in computer graphics has been carried out on both low-level, hardware-related algorithms and high-level, abstract software. Success of computer graphics has been seen in many application areas, such as multimedia, visualisation, virtual reality and the Internet. However, the hardware realisation of OpenGL architecture based on FPGA (field programmable gate array) is beyond the scope of most of computer graphics researches. It is an uncultivated research area where the OpenGL pipeline, from hardware through the whole embedded system (ES) up to applications, is implemented in an FPGA chip. This research proposes a hybrid approach to investigating both software and hardware methods. It aims at bridging the gap between methods of software and hardware, and enhancing the overall performance for computer graphics. It consists of four parts, the construction of an FPGA-based ES, Mesa-OpenGL implementation for FPGA-based ESs, parallel processing, and a novel algorithm for surface modelling and editing. The FPGA-based ES is built up. In addition to the Nios II soft processor and DDR SDRAM memory, it consists of the LCD display device, frame buffers, video pipeline, and algorithm-specified module to support the graphics processing. Since there is no implementation of OpenGL ES available for FPGA-based ESs, a specific OpenGL implementation based on Mesa is carried out. Because of the limited FPGA resources, the implementation adopts the fixed-point arithmetic, which can offer faster computing and lower storage than the floating point arithmetic, and the accuracy satisfying the needs of 3D rendering. Moreover, the implementation includes Bézier-spline curve and surface algorithms to support surface modelling and editing. The pipelined parallelism and co-processors are used to accelerate graphics processing in this research. These two parallelism methods extend the traditional computation parallelism in fine-grained parallel tasks in the FPGA-base ESs. The novel algorithm for surface modelling and editing, called Progressive and Mixing Algorithm (PAMA), is proposed and implemented on FPGA-based ES’s. Compared with two main surface editing methods, subdivision and deformation, the PAMA can eliminate the large storage requirement and computing cost of intermediated processes. With four independent shape parameters, the PAMA can be used to model and edit freely the shape of an open or closed surface that keeps globally the zero-order geometric continuity. The PAMA can be applied independently not only FPGA-based ESs but also other platforms. With the parallel processing, small size, and low costs of computing, storage and power, the FPGA-based ES provides an effective hybrid solution to surface modelling and editing

Iterative geometric design for architecture

This work investigates on computer aided integrated architectural design and production. The aim is to provide integral solutions for the design and the production of geometrically complex free-form architecture. Investigations on computer aided geometric design and integrated manufacturing are carried out with equal importance. This research is considering an integral and interdisciplinary approach, including computer science, mathematics and architecture. Inspired by fractal geometry, the IFS formalism is studied with regards to discrete architectural geometric design. The geometric design method studied provides new shape control possibilities unifying two separate design paradigms of rough and smooth objects. Capable to design fractal geometric figures, the method also covers the generation of classical objects such as conics and NURBS-curves. Close attention has been paid to the design of iterative free-form surfaces, which are composed entirely out of planar elements. A surface method based on projected vector sums is proposed. The resulting geometric figures are expressed in a discrete form and can be easily translated into a coherent set of constructional elements. The studies for translation of the geometrical elements into constructional elements consider integrated manufacturing. Addressing and numbering of the elements by iterative geometric design are investigated and compared to lexicographically ordered addressing systems, in order to provide an adequate data structure for the design, production and assembly of the constructional elements. For the generation of the data describing constructional elements, problems related to thickening and offset meshes are discussed. Once the global geometry of the constructional part has been computed, parameters are defined for generic automated detailing. Hereby the entire description of the constructional elements is completed. These elements are mapped and packed with regards to the coordinate system of a CNC-machine and the properties and the dimensions of the raw material, providing the complete set of workshop plans needed for integrated manufacturing. For automated generation of machine instructions (G-code), machining strategies – depending on the type of machine used, tool and material properties – are elaborated. Finally, the integrated digital design methods studied within the scope of this thesis are tested and verified by the realization of different reduced scale prototypes. The studied applications range from bearing vault structures to fractal and smooth timber panel shell structures. The developed methods have shown to be efficient for the design and the realization of geometrically complex architectural objects. The required planning effort to handle and manipulate the design and the production data has been greatly reduced. Some of the proposed methods have proved to be robust and general enough to be applied on real world applications. Iterative geometric design provides high degree of design possibilities offering an efficient tool for the creation of smooth and rough free form objects. The possibility to incorporate successive folds in free-form objects allows structural applications