A new method for constructing triangular patches is presented. The triangular patch that interpolates the given boundary curves and cross-boundary slopes is constructed by a basic approximation operator plus an additi...A new method for constructing triangular patches is presented. The triangular patch that interpolates the given boundary curves and cross-boundary slopes is constructed by a basic approximation operator plus an additional interpolation operator. The basic approximation operator is constructed by a polynomial surface of degree five which approximates the given interpolation conditions. The additional interpolation operator is formed by the side-vertex method. The basic and the additional operators have different roles in constructing the triangular patch: the first one makes the triangular patch approximate the given interpolation conditions with a polynomial approximation precision of degree five, while the second one makes it satisfy the given interpolation conditions. The triangular patch reproduces polynomial surfaces of degree five. Comparison results of the new method with the other two methods are included.展开更多
Mixed triangular spectral element method using nodal basis on unstructured meshes is investigated in this paper.The method is based on equivalent first order system of the elliptic problem and rectangle-triangle trans...Mixed triangular spectral element method using nodal basis on unstructured meshes is investigated in this paper.The method is based on equivalent first order system of the elliptic problem and rectangle-triangle transforms.It fully enjoys the ten-sorial structure and flexibility in handling complex domains by using nodal basis and unstructured triangular mesh.Different from the usual Galerkin formulation,the mixed form is particularly advantageous in this context,since it can avoid the singularity in-duced by the rectangle-triangle transform in the calculation of the matrices,and does not require the evaluation of the stiffness matrix.An hp a priori error estimate is pres-ented for the proposed method.The implementation details and some numerical exam-ples are provided to validate the accuracy and flexibility of the method.展开更多
基金This work was supported by the National Natural Science Foundation of China(Grant No.601 73052)Shandong Province Natural Science Foundation(Grant No.Z2001G01)Doctoral Program of High Education of China(Grant No.20020422030).
文摘A new method for constructing triangular patches is presented. The triangular patch that interpolates the given boundary curves and cross-boundary slopes is constructed by a basic approximation operator plus an additional interpolation operator. The basic approximation operator is constructed by a polynomial surface of degree five which approximates the given interpolation conditions. The additional interpolation operator is formed by the side-vertex method. The basic and the additional operators have different roles in constructing the triangular patch: the first one makes the triangular patch approximate the given interpolation conditions with a polynomial approximation precision of degree five, while the second one makes it satisfy the given interpolation conditions. The triangular patch reproduces polynomial surfaces of degree five. Comparison results of the new method with the other two methods are included.
基金The first and second authors gratefully acknowledge the financial support provided by NSFC(grant 11771137)。
文摘Mixed triangular spectral element method using nodal basis on unstructured meshes is investigated in this paper.The method is based on equivalent first order system of the elliptic problem and rectangle-triangle transforms.It fully enjoys the ten-sorial structure and flexibility in handling complex domains by using nodal basis and unstructured triangular mesh.Different from the usual Galerkin formulation,the mixed form is particularly advantageous in this context,since it can avoid the singularity in-duced by the rectangle-triangle transform in the calculation of the matrices,and does not require the evaluation of the stiffness matrix.An hp a priori error estimate is pres-ented for the proposed method.The implementation details and some numerical exam-ples are provided to validate the accuracy and flexibility of the method.