An isogeometric analysis for elliptic homogenization problems

A novel and efficient approach which is based on the framework of isogeometric analysis for elliptic homogenization problems is proposed. These problems possess highly oscillating coefficients leading to extremely high computational expenses while using traditional finite element methods. The isogeometric analysis heterogeneous multiscale method (IGA-HMM) investigated in this paper is regarded as an alternative approach to the standard Finite Element Heterogeneous Multiscale Method (FE-HMM) which is currently an effective framework to solve these problems. The method utilizes non-uniform rational B-splines (NURBS) in both macro and micro levels instead of standard Lagrange basis. Beside the ability to describe exactly the geometry, it tremendously facilitates high-order macroscopic/microscopic discretizations thanks to the flexibility of refinement and degree elevation with an arbitrary continuity level provided by NURBS basis functions. A priori error estimates of the discretization error coming from macro and micro meshes and optimal micro refinement strategies for macro/micro NURBS basis functions of arbitrary orders are derived. Numerical results show the excellent performance of the proposed method

Solving boundary value problems via the Nyström method using spline Gauss rules

We propose to use spline Gauss quadrature rules for solving boundary value problems (BVPs) using the Nyström method. When solving BVPs, one converts the corresponding partial differential equation inside a domain into the Fredholm integral equation of the second kind on the boundary in the sense of boundary integral equation (BIE). The Fredholm integral equation is then solved using the Nyström method, which involves the use of a particular quadrature rule, thus, converting the BIE problem to a linear system. We demonstrate this concept on the 2D Laplace problem over domains with smooth boundary as well as domains containing corners. We validate our approach on benchmark examples and the results indicate that, for a fixed number of quadrature points (i.e., the same computational effort), the spline Gauss quadratures return an approximation that is by one to two orders of magnitude more accurate compared to the solution obtained by traditional polynomial Gauss counterparts

Solving Boundary Value Problems Via the Nyström Method Using Spline Gauss Rules

We propose to use spline Gauss quadrature rules for solving boundary value problems (BVPs) using the Nyström method. When solving BVPs, one converts the corresponding partial differential equation inside a domain into the Fredholm integral equation of the second kind on the boundary in the sense of boundary integral equation (BIE). The Fredholm integral equation is then solved using the Nyström method, which involves a use of a particular quadrature rule, thus, converting the BIE problem to a linear system. We demonstrate this concept on the 2D Laplace problem over domains with smooth boundary as well as domains containing corners. We validate our approach on benchmark examples and the results indicate that, for a fixed number of quadrature points (i.e., the same computational effort), the spline Gauss quadratures return an approximation that is by one to two orders of magnitude more accurate compared to the solution obtained by traditional polynomial Gauss counterparts

Numerical Solutions of Sixth Order Linear and Nonlinear Boundary Value Problems

The aim of paper is to find the numerical solutions of sixth order linear and nonlinear differential equations with two point boundary conditions. The well known Galerkin method with Bernstein and modified Legendre polynomials as basis functions is exploited. In this method, the basis functions are transformed into a new set of basis functions, which satisfy the homogeneous form of Dirichlet boundary conditions. A rigorous matrix formulation is derived for solving the sixth order BVPs. Several numerical examples are considered to verify the efficiency and implementation of the proposed method. The numerical results are compared with both the exact solutions and the results of the other methods available in the literature. The comparison shows that the performance of the present method is more efficient and yields better results

Smooth Subdivision Surfaces: Mesh Blending and Local Interpolation

Subdivision surfaces are widely used in computer graphics and animation. Catmull-Clark subdivision (CCS) is one of the most popular subdivision schemes. It is capable of modeling and representing complex shape of arbitrary topology. Polar surface, working on a triangle-quad mixed mesh structure, is proposed to solve the inherent ripple problem of Catmull-Clark subdivision surface (CCSS). CCSS is known to be C1 continuous at extraordinary points. In this work, we present a G2 scheme at CCS extraordinary points. The work is done by revising CCS subdivision step with Extraordinary-Points-Avoidance model together with mesh blending technique which selects guiding control points from a set of regular sub-meshes (named dominative control meshes) iteratively at each subdivision level. A similar mesh blending technique is applied to Polar extraordinary faces of Polar surface as well. Both CCS and Polar subdivision schemes are approximating. Traditionally, one can obtain a CCS limit surface to interpolate given data mesh by iteratively solving a global linear system. In this work, we present a universal interpolating scheme for all quad subdivision surfaces, called Bezier Crust. Bezier Crust is a specially selected bi-quintic Bezier surface patch. With Bezier Crust, one can obtain a high quality interpolating surface on CCSS by parametrically adding CCSS and Bezier Crust. We also show that with a triangle/quad conversion process one can apply Bezier Crust on Polar surfaces as well. We further show that Bezier Crust can be used to generate hollowed 3D objects for applications in rapid prototyping. An alternative interpolating approach specifically designed for CCSS is developed. This new scheme, called One-Step Bi-cubic Interpolation, uses bicubic patches only. With lower degree polynomial, this scheme is appropriate for interpolating large-scale data sets. In sum, this work presents our research on improving surface smoothness at extraordinary points of both CCS and Polar surfaces and present two local interpolating approaches on approximating subdivision schemes. All examples included in this work show that the results of our research works on subdivision surfaces are of high quality and appropriate for high precision engineering and graphics usage

Smooth path planning with Pythagorean-hodoghraph spline curves geometric design and motion control

This thesis addresses two significative problems regarding autonomous systems, namely path and trajectory planning. Path planning deals with finding a suitable path from a start to a goal position by exploiting a given representation of the environment. Trajectory planning schemes govern the motion along the path by generating appropriate reference (path) points. We propose a two-step approach for the construction of planar smooth collision-free navigation paths. Obstacle avoidance techniques that rely on classical data structures are initially considered for the identification of piecewise linear paths that do not intersect with the obstacles of a given scenario. In the second step of the scheme we rely on spline interpolation algorithms with tension parameters to provide a smooth planar control strategy. In particular, we consider Pythagorean\u2013hodograph (PH) curves, since they provide an exact computation of fundamental geometric quantities. The vertices of the previously produced piecewise linear paths are interpolated by using a G1 or G2 interpolation scheme with tension based on PH splines. In both cases, a strategy based on the asymptotic analysis of the interpolation scheme is developed in order to get an automatic selection of the tension parameters. To completely describe the motion along the path we present a configurable trajectory planning strategy for the offline definition of time-dependent C2 piece-wise quintic feedrates. When PH spline curves are considered, the corresponding accurate and efficient CNC interpolator algorithms can be exploited

Smooth path planning with Pythagorean-hodoghraph spline curves geometric design and motion control

This thesis addresses two significative problems regarding autonomous systems, namely path and trajectory planning. Path planning deals with finding a suitable path from a start to a goal position by exploiting a given representation of the environment. Trajectory planning schemes govern the motion along the path by generating appropriate reference (path) points. We propose a two-step approach for the construction of planar smooth collision-free navigation paths. Obstacle avoidance techniques that rely on classical data structures are initially considered for the identification of piecewise linear paths that do not intersect with the obstacles of a given scenario. In the second step of the scheme we rely on spline interpolation algorithms with tension parameters to provide a smooth planar control strategy. In particular, we consider Pythagorean–hodograph (PH) curves, since they provide an exact computation of fundamental geometric quantities. The vertices of the previously produced piecewise linear paths are interpolated by using a G1 or G2 interpolation scheme with tension based on PH splines. In both cases, a strategy based on the asymptotic analysis of the interpolation scheme is developed in order to get an automatic selection of the tension parameters. To completely describe the motion along the path we present a configurable trajectory planning strategy for the offline definition of time-dependent C2 piece-wise quintic feedrates. When PH spline curves are considered, the corresponding accurate and efficient CNC interpolator algorithms can be exploited

Quantifying Initial Condition and Parametric Uncertainties in a Nonlinear Aeroelastic System with an Efficient Stochastic Algorithm

There is a growing interest in understanding how uncertainties in flight conditions and structural parameters affect the character of a limit cycle oscillation (LCO) response, leading to failure of an aeroelastic system. Uncertainty quantification of a stochastic system (parametric uncertainty) with stochastic inputs (initial condition uncertainty) has traditionally been analyzed with Monte Carlo simulations (MCS). Probability density functions (PDF) of the LCO response are obtained from the MCS to estimate the probability of failure. A candidate approach to efficiently estimate the PDF of an LCO response is the stochastic projection method. The objective of this research is to extend the stochastic projection method to include the construction of B-spline surfaces in the stochastic domain. The multivariate B-spline problem is solved to estimate the LCO response surface. An MCS is performed on this response surface to estimate the PDF of the LCO response. The probability of failure is then computed from the PDF. This method is applied to the problem of estimating the PDF of a subcritical LCO response of a nonlinear airfoil in inviscid transonic flow. The stochastic algorithm provides a conservative estimate of the probability of failure of this aeroelastic system two orders of magnitude more efficiently than performing an MCS on the governing equations

Theory and applications of bijective barycentric mappings

Barycentric coordinates provide a convenient way to represent a point inside a triangle as a convex combination of the triangle's vertices, and to linearly interpolate data given at these vertices. Due to their favourable properties, they are commonly applied in geometric modelling, finite element methods, computer graphics, and many other fields. In some of these applications it is desirable to extend the concept of barycentric coordinates from triangles to polygons. Several variants of such generalized barycentric coordinates have been proposed in recent years. An important application of barycentric coordinates consists of barycentric mappings, which allow to naturally warp a source polygon to a corresponding target polygon, or more generally, to create mappings between closed curves or polyhedra. The principal practical application is image warping, which takes as input a control polygon drawn around an image and smoothly warps the image by moving the polygon vertices. A required property of image warping is to avoid fold-overs in the resulting image. The problem of fold-overs is a manifestation of a larger problem related to the lack of bijectivity of the barycentric mapping. Unfortunately, bijectivity of such barycentric mappings can only be guaranteed for the special case of warping between convex polygons or by triangulating the domain and hence renouncing smoothness. In fact, for any barycentric coordinates, it is always possible to construct a pair of polygons such that the barycentric mapping is not bijective. In the first part of this thesis we illustrate three methods to achieve bijective mappings. The first method is based on the intuition that, if two polygons are sufficiently close, then the mapping is close to the identity and hence bijective. This suggests to ``split'' the mapping into several intermediate mappings and to create a composite barycentric mapping which is guaranteed to be bijective between arbitrary polygons, polyhedra, or closed planar curves. We provide theoretical bounds on the bijectivity of the composite mapping related to the norm of the gradient of the coordinates. The fact that the bound depends on the gradient implies that these bounds exist only if the gradient of the coordinates is bounded. We focus on mean value coordinates and analyse the behaviour of their directional derivatives and gradient at the vertices of a polygon. The composition of barycentric mappings for closed planar curves leads to the problem of blending between two planar curves. We suggest to solve it by linearly interpolating the signed curvature and then reconstructing the intermediate curve from the interpolated curvature values. However, when both input curves are closed, this strategy can lead to open intermediate curves. We present a new algorithm for solving this problem, which finds the closed curve whose curvature is closest to the interpolated values. Our method relies on the definition of a suitable metric for measuring the distance between two planar curves and an appropriate discretization of the signed curvature functions. The second method to construct smooth bijective mappings with prescribed behaviour along the domain boundary exploits the properties of harmonic maps. These maps can be approximated in different ways, and we discuss their respective advantages and disadvantages. We further present a simple procedure for reducing their distortion and demonstrate the effectiveness of our approach by providing examples. The last method relies on a reformulation of complex barycentric mappings, which allows us to modify the ``speed'' along the edges to create complex bijective mappings. We provide some initial results and an optimization procedure which creates complex bijective maps. In the second part we provide two main applications of bijective mapping. The first one is in the context of finite elements simulations, where the discretization of the computational domain plays a central role. In the standard discretization, the domain is triangulated with a mesh and its boundary is approximated by a polygon. We present an approach which combines parametric finite elements with smooth bijective mappings, leaving the choice of approximation spaces free. This approach allows to represent arbitrarily complex geometries on coarse meshes with curved edges, regardless of the domain boundary complexity. The main idea is to use a bijective mapping for automatically warping the volume of a simple parametrization domain to the complex computational domain, thus creating a curved mesh of the latter. The second application addresses the meshing problem and the possibility to solve finite element simulations on polygonal meshes. In this context we present several methods to discretize the bijective mapping to create polygonal and piece-wise polynomial meshes