Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization

In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries consisting of a set of B-spline curves. Instead of forming the computational domain by a simple boundary, planar domains with high genus and more complex boundary curves are considered. Firstly, some pre-processing operations including B\'ezier extraction and subdivision are performed on each boundary curve in order to generate a high-quality planar parameterization; then a robust planar domain partition framework is proposed to construct high-quality patch-meshing results with few singularities from the discrete boundary formed by connecting the end points of the resulting boundary segments. After the topology information generation of quadrilateral decomposition, the optimal placement of interior B\'ezier curves corresponding to the interior edges of the quadrangulation is constructed by a global optimization method to achieve a patch-partition with high quality. Finally, after the imposition of C1=G1-continuity constraints on the interface of neighboring B\'ezier patches with respect to each quad in the quadrangulation, the high-quality B\'ezier patch parameterization is obtained by a C1-constrained local optimization method to achieve uniform and orthogonal iso-parametric structures while keeping the continuity conditions between patches. The efficiency and robustness of the proposed method are demonstrated by several examples which are compared to results obtained by the skeleton-based parameterization approach

Analysis-suitable volume parameterization of multi-block computational domain in isogeometric applications

International audienceParameterization of computational domain is a key step in isogeometric analysis just as mesh generation is in finite element analysis. In this paper, we study the volume parameterization problem of multi-block computational domain in isogeometric version, i.e, how to generate analysis-suitable parameterization of the multi-block computational domain bounded by B-spline surfaces. Firstly, we show how to find good volume parameterization of single-block computational domain by solving a constraint optimization problem, in which the constraint condition is the injectivity sufficient conditions of B-spline volume parametrization, and the optimization term is the minimization of quadratic energy functions related to the first and second derivatives of B-spline volume parameterization. By using this method, the resulted volume parameterization has no self-intersections, and the isoparametric structure has good uniformity and orthogonality. Then we extend this method to the multi-block case, in which the continuity condition between the neighbor B-spline volume should be added to the constraint term. The effectiveness of the proposed method is illustrated by several examples based on three-dimensional heat conduction problem

Optimal analysis-aware parameterization of computational domain in 3D isogeometric analysis

International audienceIn isogeometric analysis framework, computational domain is exactly described using the same representation as that employed in the CAD process. For a CAD object, we can construct various computational domain with same shape but with different parameterization. One basic requirement is that the resulting parameterization should have no self-intersections. In this paper, a linear and easy-to-check sufficient condition for injectivity of trivariate B-spline parameterization is proposed. By an example of 3D thermal conduction problem, we show that different parameterization of computational domain has different impact on the simulation result and efficiency in isogeometric analysis. For problems with exact solutions, we propose a shape optimization method to obtain optimal parameterization of computational domain. The proposed injective condition is used to check the injectivity of initial trivariate B-spline parameterization constructed by discrete Coons volume method, which is the generalization of discrete Coons patch method. Several examples and comparisons are presented to show the effectiveness of the proposed method. Compared with the initial parameterization during refinement, the optimal parameterization can achieve the same accuracy but with less degrees of freedom

Constructing analysis-suitable parameterization of computational domain from CAD boundary by variational harmonic method

In isogeometric anlaysis, parameterization of computational domain has great effects as mesh generation in finite element analysis. In this paper, based on the concept of harmonic mapping from the computational domain to parametric domain, a variational harmonic approach is proposed to construct analysis-suitable parameterization of computational domain from CAD boundary for 2D and 3D isogeometric applications. Different from the previous elliptic mesh generation method in finite element analysis, the proposed method focus on isogeometric version, and converts the elliptic PDE into a nonlinear optimization problem, in which a regular term is integrated into the optimization formulation to achieve more uniform and orthogonal isoparametric structure near convex (concave) parts of the boundary. Several examples are presented to show the efficiency of the proposed method in 2D and 3D isogeometric analysis

Parametrization of computational domain in isogeometric analysis: methods and comparison

International audienceParameterization of computational domain plays an important role in isogeometric analysis as mesh generation in finite element analysis. In this paper, we investigate this problem in the 2D case, i.e, how to parametrize the computational domains by planar B-spline surface from the given CAD objects (four boundary planar B-spline curves). Firstly, two kinds of sufficient conditions for injective B-spline parameterization are derived with respect to the control points. Then we show how to find good parameterization of computational domain by solving a constraint optimization problem, in which the constraint condition is the injectivity sufficient conditions of planar B-spline parametrization, and the optimization term is the minimization of quadratic energy functions related to the first and second derivatives of planar B-spline parameterization. By using this method, the resulted parameterization has no self-intersections, and the isoparametric net has good uniformity and orthogonality. After introducing a posteriori error estimation for isogeometric analysis, we propose $r$-refinement method to optimize the parameterization by repositioning the inner control points such that the estimated error is minimized. Several examples are tested on isogeometric heat conduction problem to show the effectiveness of the proposed methods and the impact of the parameterization on the quality of the approximation solution. Comparison examples with known exact solutions are also presented

Hermite type Spline spaces over rectangular meshes with complex topological structures

International audienceMotivated by the magneto hydrodynamic (MHD) simulation for Tokamaks with Isogeometric analysis, we present a new type of splines defined over a rectangular mesh with arbitrary topology, which are piecewise polynomial functions of bidegree (d,d) and C^r parameter continuity. In particular, We compute their dimension and exhibit basis functions called Hermite bases for bicubic spline spaces. We investigate their potential applications for solving partial differential equations (PDEs) over a complex physical domain in the framework of Isogeometric analysis. In particular, weanalyze the property of approximation of these spline spaces for the L2-norm. Despite the fact that the basis functions are singular at extraordinary vertices, we show that the optimal approximation order and numerical convergence rates are reached by setting a proper parameterization

Exact conversion from Bézier tetrahedra to Bézier hexahedra

International audienceModeling and computing of trivariate parametric volumes is an important research topic in the field of three-dimensional isogeo-metric analysis. In this paper, we propose two kinds of exact conversion approaches from Bézier tetrahedra to Bézier hexahedra with the same degree by reparametrization technique. In the first method, a Bézier tetrahedron is converted into a degenerate Bézier hexahedron, and in the second approach, a non-degenerate Bézier tetrahedron is converted into four non-degenerate Bézier hexahedra. For the proposed methods, explicit formulas are given to compute the control points of the resulting tensor-product Bézier hexahedra. Furthermore, in the second method, we prove that tetrahedral spline solids with C k-continuity can be converted into a set of tensor-product Bézier volumes with G k-continuity. The proposed methods can be used for the volumetric data exchange problems between different trivariate spline representations in CAD/CAE. Several experimental results are presented to show the effectiveness of the proposed methods