Next Article in Journal
Global Well-Posedness and Analyticity of Generalized Porous Medium Equation in Fourier-Besov-Morrey Spaces with Variable Exponent
Next Article in Special Issue
Study of the Boundary Value Problems for Nonlinear Wave Equations on Domains with a Complex Structure of the Boundary and Prehistory
Previous Article in Journal
Neutrosophic Completion Technique for Incomplete Higher-Order AHP Comparison Matrices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discontinuous Galerkin Isogeometric Analysis of Convection Problem on Surface

1
Department of Mathematics, Beijing Institute of Technology, Beijing 100081, China
2
State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, China Meteorological Administration, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Submission received: 20 January 2021 / Revised: 24 February 2021 / Accepted: 25 February 2021 / Published: 28 February 2021
(This article belongs to the Special Issue Advanced Methods in Computational Mathematical Physics)

Abstract

:
The objective of this work is to study finite element methods for approximating the solution of convection equations on surfaces embedded in R 3 . We propose the discontinuous Galerkin (DG) isogeometric analysis (IgA) formulation to solve convection problems on implicitly defined surfaces. Three numerical experiments shows that the numerical scheme converges with the optimal convergence order.

1. Introduction

Surface partial differential equations (SPDEs) arise in natural sciences and applied areas, such as minimal surface equation, Willmore flow, transport of surfactants along interfaces in multiphase fluids [1], and lipid interactions in cell membranes [2]. In this paper, we consider the DG-IgA methods for the following model problem:
Γ · ( β u ) + γ u = f ,
where Γ · = div Γ denotes the surface divergence, Γ is a smooth two-dimensional surface without boundary embedded in R 3 . β = ( b 1 , b 2 , b 3 ) is the advective velocity with b i W 1 ( Γ ) , γ L ( Γ ) is reaction coefficient. We call this model problem the convection equation on surface Γ . The problem (1) has a unique solution u H ( d i v , Γ ) satisfying Γ u d σ = 0 , where d σ is surface measure. Moreover, for source term f L 2 ( Γ ) we have Γ f d σ = 0 .
To ensure the existence and uniqueness of a solution u H 1 ( Γ ) to the model problem (1.1), we adopt the following hypothesis: There exists a positive constant γ 0 such that
γ ( x ) + 1 2 Γ · β ( x ) = ( c 0 ( x ) ) 2 γ 0 , f o r a l l x Γ .
Finite element methods of SPDEs have been extensively developed. G. Dziuk [3] first used the finite element method for solving the surface elliptic problem. A. Bonito and J. E. Pasciak [4] designed and analyzed multigrid algorithms for the Laplace-Beltrami operator on a smooth and closed surface. DG methods for solving the first-order hyperbolic problems were defined in [5,6,7,8]. Since DG methods have well-known stability and local conservation in numerical application of partial differential equation problems [9,10] and capture solution blow-ups, it is natural to extend the DG approximation to PDEs on surfaces or manifolds. Ref. [11] extended the discontinuous Galerkin (DG) methods to the first-order hyperbolic and advection-dominated problems on surfaces.
IgA has been introduced in [12] as a new tool for solving numerically PDEs with the complicated geometric domains, in particular, surfaces. DG-IgA methods have been applied for approximating solutions of elliptic problems on surfaces [13,14]. Ref. [15] developed DG-IgA numerical schemes for solving problems on segmentations with gaps. Ref. [16] studied NURBS-based isogeometric analysis for the computation of flows about rotating components. In [17], a new stabilized symmetric Nitsche method was proposed for enforcement of Dirichlet boundary conditions for elliptic problems in Cut-IgA. The remainder of the paper is organized as follows: First, we introduce some preliminaries about IgA and discrete NURBS finite space V ( h ) . Then we derive the DG-IgA scheme of convection problem (1.1). Finally, we present three numerical experiments to illustrate the discrete formulation.

2. Model Problem and Discretization

2.1. Differential Operators on Surfaces

Let us consider a two-dimensional surface Γ R 3 defined in the physical space R 3 . Assume that the surface is characterized by a geometrical mapping from a parameter space Γ ^ R 2 . Let X C 2 ( Γ ^ , R 3 ) be a local parameterization of θ Γ ^ , where the vector-valued independent variable θ = ( θ 1 , θ 2 ) is called parametric coordinate. Indeed, X defines the geometrical mapping as follows:
X : Γ ^ Γ R 3 , θ X ( θ ) .
To define surface gradient, we introduce the Jacobi matrix J ( θ ) = ( J i j ) 3 × 2 as
J ( θ ) = X θ 1 , X θ 2 .
So, the metric tensor G ( θ ) = ( g i j ) 2 × 2 is represented by
g i j ( θ ) = J ( θ ) T J ( θ ) = X ( θ ) θ i X ( θ ) θ j , i , j = 1 , 2 .
The inverse of G ( θ ) is denoted by G 1 ( θ ) = ( g i j ) 2 × 2 . The tangential or surface gradient is given by
Γ u ( X ( θ ) ) = G 1 ( θ ) d u ( X ( θ ) ) d θ .
So, we can deduce the tangential or surface divergence. For υ is a C 1 vector field on Γ , we have
Γ · υ = 1 | G | θ 1 , θ 2 ( | G | G 1 J T υ ) .
Remark that the vector field υ may not be tangential to surface Γ . However, the 2-dimensional vector G 1 J T υ is the projection of υ in parametric coordinate θ , which is a tangent vector field of Γ . So, the definition of surface divergence is well-posed.
Let L 2 ( Γ ) denote the usual L 2 space on the surface Γ with norm
f L 2 ( Γ ) = Γ f 2 d σ 1 2 .
Furthermore, we use standard notations W p m ( Γ ) and H m ( Γ ) for Sobolev spaces on Γ with norm and semi-norm
f W p m ( Γ ) p = k = 0 m f W p k ( Γ ) p a n d f W p m ( Γ ) = D Γ m f L p ( Γ ) , f H m ( Γ ) 2 = k = 0 m f H k ( Γ ) 2 a n d f H k ( Γ ) = D Γ k f L 2 ( Γ ) .
To present the weak form of problem (1), we introduce the Stokes theorem on surface Γ . Let υ be a C 1 vector field on surface Γ , q C 1 ( Γ ) , then the following identity holds:
Γ ( υ · Γ q + q Γ · υ ) d σ = Γ q υ · n d s .

2.2. Isogeometric Analysis

To apply the IgA methodology of the problem (1), the domain Γ is partitioned into some closed sub-patches such that
Γ = i = 1 N Γ i w i t h Γ i Γ j = , i f i j .
We denote the set of sub-patches as S ( Γ ) = { Γ i } i = 1 N .
Without loss of generality, we simply assume a parametric domain D ^ of unit length, i.e., D ^ = [ 0 , 1 ] 2 . For each patch Γ i , we associate the knot vectors V i on D ^ , which defines a partition T ^ h , D ^ i = { K ^ m } m = 1 M i , where K ^ m are 2-dimensional closed elements. We refer to T ^ h , D ^ i as the parametric mesh of patch Γ i and denote T ^ h = i = 1 N T ^ h , D ^ i . Any patch Γ i can be represented by a parametrization map as follows
Φ i : D ^ Γ i , Φ i ( x ^ ) = j = 1 N i C j ( i ) B ^ j ( i ) ( x ^ ) = x Γ i ,
where C j ( i ) are the control points and B ^ j ( i ) are B-spline basis functions or NURBS of degree k [18], N i denotes the number of basis functions B ^ j ( i ) on patch Γ i .
We obtain the partition T h , Γ i = { K m } m = 1 M i of the patch Γ i , whose vertices are the images of the vertices of the corresponding parametric partition T ^ h , D ^ i by the map Φ i . Now we can construct the partition of the domain Γ , denoted by T h as follows,
T h ( Γ ) = i = 1 N T h , Γ i .
The set of all the edges of partition T h is denoted by E as follows,
E = K 1 , K 2 T h ( K 1 K 2 ) .
We denote the faces of all patches as F defined by
F = i { e E : e Γ i } .
Let h K and h e denote the size of element K T h and the length of the face e E respectively. h Γ i = max K Γ i { h K } is the mesh size of patch Γ i . The global mesh size of partition T h is defined as
h = max 1 i N { h Γ i } .
We assume that the shape of the elements is regular and quasi-uniform, i.e., h e h and h K h Γ i .
Next, we define the broken space on the physical domain Γ associated with T h by using the introduced push-forward function Ψ i ( x ) = Φ i 1 ( x ) = x ^ , for any u ^ H s ( D ^ ) ,
u ( x ) = u ^ ( Ψ i ( x ) ) , x Γ i .
Now we can define the broken Sobolev space H k ( T h ) :
H k ( T h ) = { v L 2 ( Γ ) : v | Γ i H k ( Γ i ) , i = 1 , 2 , , N } .
To apply DG method to problem (1), we introduce the jump [ u h ] of u h on F i j . For u h H 1 ( T h ) , using the notation u h i = u h | Γ i , we define
[ u h ] = u h i · n i + u h j · n j , o n F i j ,
where n i denotes the unit normal vector of Γ i on F i j pointing exterior to Γ j .

2.3. DG-IgA Discretization

Next, we introduce the finite element space associated with the partition T h . In general, the basis function B j ( i ) ( x ) of IgA on the patch Γ i are pushed forward from the basis function B ^ j ( i ) ( x ^ ) of the parametric domain D ^ by considering a composition with Ψ i = Φ i 1 , i.e., for any basis function B ^ j ( i ) ,
B j ( i ) ( x ) = B ^ j ( i ) ( Ψ i ( x ) ) .
The DG finite element space V h is defined by
V h = u h L 2 ( T h ) : u h | Γ i Span B j ( i ) , j = 1 , , N i , i = 1 , , N .
The DG-IgA approximation is formulated as follows: Find u h V h such that
a h ( u h , v h ) = F ( v h ) , v h V h ,
where
a h ( u h , v h ) = K T h K ( β · Γ v h + γ v h ) u h d σ + e E e { u h } [ β v h ] d s + e η e 2 [ u h ] [ v h ] d s , F ( v h ) = K T h K f v h d σ .
where η e is the stability parameter, defined as
η e = 1 2 ( | β ( i ) · n i | + | β ( j ) · n j | ) , β ( i ) = β ( j ) , | β · n | , o t h e r w i s e ,
for e = K i K j .
The bilinear map a h ( u h , v h ) can be modified. Indeed, by using the continuity of β and the function of V h in each patch, i.e., [ v h ] | K \ Γ i = 0 and [ β v h ] | K \ Γ i = 0 for any v h V h , we have
e E e [ β v h ] { u h } d s + e η e 2 [ u h ] [ v h ] d s = F i j F F i j [ β v h ] { u h } d s + F i j η e 2 [ u h ] [ v h ] d s .
Similarly, we can define the outflow faces F + . So, the bilinear map a h ( u h , v h ) can be modified as follows
a h ( u h , v h ) = K T h K ( β · Γ v h + γ v h ) u h d σ + F i j F F i j [ β v h ] { u h } d s + F i j η e 2 [ u h ] [ v h ] d s .
We can derive the discrete coercivity, stability and consistency of above numerical scheme with a similar technique in [19]. Due to the discontinuity of β , it needs a few skills to prove the discrete coercivity. So there exists a unique discrete solution u h V h satisfying DG-IgA scheme (5). Here we briefly show the a priori error result without detailed proof. We introduce the dual or adjoint weak form: Find ω H s ( T h ) such that
K T h ( ω , · ( β v ) + γ v ) K F i j F ( [ β v ] , { Π h ω } ) F i j + F i j F ( η e 2 [ v ] , [ ω ] ) F i j = ( u h Π h u , v ) ,
for all v V ( h ) . With the above adjoint weak form, inverse inequality and the approximation of the interpolant operator, we can prove the following result.
Theorem 1.
Let u and u h denote the solutions of (1) and (5), respectively. Assume that u | K H k + 1 ( K ) , ω H k + 1 ( K ) , K T h . Then we have
u u h L 2 ( Γ ) C h k + 1 u H k + 1 ( Γ ) .

3. Numerical Experiments

In this section, we present some numerical experiments of convection problems on surface. Numerical examples are presented for a sphere and a quarter of a cylinder.

3.1. Numerical Experiment 1

We first consider the model problem (1) on the unit sphere
Γ · ( β u ) + γ u = f , in Γ
subject to the compatibility condition Γ f d σ = 0 , where β = [ 1 , 1 , 1 ] T and γ = 1 . We select the source function f such that the exact solution is u = sin ( θ ) sin ( ϕ ) , where ( ϕ , θ ) are the spherical coordinates.
We divide the unit sphere Γ into 6 patches. For each patch, the knot vector is taken as [ 0 , 0 , 0 , 1 , 1 , 1 ] × [ 0 , 0 , 0 , 1 , 1 , 1 ] to represent the geometry of each patch. We generate the mesh by refining the parameter element of each patch, whose mesh size is denoted by h. We show the patches and the uniform meshes of the sphere for h = 1 / 4 in Figure 1(Left).
The numerical L 2 errors and convergence results are given in Table 1 for k = 1 , 2 and 3, respectively. In Figure 2, we present the convergence histories of L 2 errors. These results show that the IgA-DG method yields O ( h k + 1 ) convergent solution. We present the numerical solution of convection problem (1) on the sphere in Figure 1(Right).

3.2. Numerical Experiment 2

Here we continue to consider the model problem on the surface of torus
Γ · ( β u ) + γ u = f , in Γ ,
subject to the compatibility condition Γ f d σ = 0 , where β = [ 1 , 1 , 1 ] T and γ = 1 . The torus is the surface
Γ = { ( x , y ) ( 3 , 3 ) 2 , z ( 1 , 1 ) : r 2 = z 2 + ( x 2 + y 2 R 2 ) }
with r = 1 and R = 2 . We take coordinates ( ϕ , θ ) as
ϕ = arctan y x , θ = arctan z x 2 + y 2 R ,
and select the source function f such that the exact solution is u = sin ( θ ) sin ( 3 ϕ ) .
We divide the torus into 8 patches. For each patch, we take the knot vector as [ 0 , 0 , 0 , 1 / 2 , 1 / 2 , 1 , 1 , 1 ] × [ 0 , 0 , 0 , 1 , 1 , 1 ] to give the geometrical representation. We plot the patches and the uniform meshes of the torus for h = 1 / 4 in Figure 3(Left).
The numerical L 2 errors and convergence rates of this problem for k = 1 , 2 and 3 are shown in Table 2. Table 3 indicate that the rates are also O ( h k + 1 ) for L 2 norm. Figure 4 shows the convergence history of errors. Finally, we plot the numerical solution of convection problem (1) on torus in Figure 3(Right).

3.3. Numerical Experiment 3

Next, we solve on the surface of a square of the cylinder the model problem
Γ · ( β u ) + γ u = f , in Γ , u = g , o n Γ ,
where γ = 1 . The domain Γ is the surface of a quarter of the cylinder, shown in Figure 5. In contrast to the case of sphere, this problem needs the boundary condition Γ determined by β . Then, writing ( r = 1 , θ , z ) to denote the system of cylindrical coordinate, we impose an appropriate boundary condition g for u and source function f so that the exact solution is u = sin ( z θ ) .
We consider the model problem on the surface of a square of the cylinder with continuous and discontinuous coefficient β .
(1) We take continuous coefficient as β = [ 1 , 1 , 1 ] T . We divide the cylinder into 9 patches. For each patch, we take the knot vector as [ 0 , 0 , 0 , 1 , 1 , 1 ] × [ 0 , 0 , 1 , 1 ] to give the geometrical representation. Similarly, we plot the patches and the uniform meshes of the cylinder for h = 1 / 4 in Figure 5(Left).
We present numerical L 2 errors and convergence rates of this problem for k = 1 , 2 and 3 in Table 3. Table 3 indicate that the rates are also O ( h k + 1 ) for L 2 norm. Figure 6 shows the convergence history of errors. Finally, we plot the numerical solution of convection problem (1) on cylinder in Figure 5(Right).
(2) We consider the discontinuous advective velocity β as
β ( 1 ) = β ( 4 ) = [ 1 ; 1 ; 1 ] , β ( 2 ) = β ( 3 ) = [ 2 ; 4 ; 3 ] ,
where the index ( i ) is the patch number. Observe that the source function is discontinuous across the patches according to the choice of β .
According to the discontinuity of β , we divide the cylinder into 4 patches. For each patch, we take the knot vector as [ 0 , 0 , 0 , 1 , 1 , 1 ] × [ 0 , 0 , 1 , 1 ] to give the geometrical representation. Similarly, we plot the patches and the uniform meshes of the cylinder for h = 1 / 4 in Figure 7(Left), where patch 1 and patch 4 are painted blue, and patch 2, 3 are painted yellow.
The numerical L 2 errors and convergence rates of this problem for k = 1 , 2 and 3 are shown in Table 4. Table 4 indicate that the rates are also O ( h k + 1 ) for L 2 norm. Figure 8 shows the convergence history of errors. Finally, we plot the numerical solution of convection problem (1) on cylinder in Figure 7(Right).

4. Conclusions

In this paper, we present the new penalty discontinuous Galerkin (DG) isogeometric analysis(IgA) methods to solve convection problems with continuous or discontinuous coefficient on implicitly defined surfaces. For further purpose, it is worthy studying the stability and error analysis of this method and more practical problems on surfaces.

Author Contributions

Conceptualization and methodology, C.X.; formal analysis, software and validation, L.W.; visualization, X.Y.; writing—review and editing, H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Basic Research and Development Program of China grant number 2017YFC1502201, NSFC grant number 11672032, NSFC grant number 10871218 and NSFC grant number 61873036.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Milliken, W.J.; Stone, H.A.; Leal, L.G. The effect of surfactant on transient motion of newtonian drops. Phys. Fluids A 1993, 5, 69–79. [Google Scholar] [CrossRef] [Green Version]
  2. Elliott, C.M.; Stinner, B. Modeling and computation of two phase geometric biomembranes using surface finite elements. J. Comput. Phys. 2010, 229, 6585–6612. [Google Scholar] [CrossRef] [Green Version]
  3. Dziuk, G. Finite elements for the Beltrami operator on arbitrary surfaces. Partial differential equations and calculus of variations. Lect. Notes Math. 1988, 1357, 142–155. [Google Scholar]
  4. Bonito, A.; Pasciak, J.E. Convergence analysis of variational and non-variational multi-grid algorithms for the Laplace-Beltrami operator. Math. Comp. 2012, 81, 1263–1288. [Google Scholar] [CrossRef]
  5. Houston, P.; Schwab, C.; Süli, E. Discontinuous hp-finite element methods for advection-diffusion-reaction problems. SIAM J. Numer. Anal. 2002, 39, 2133–2163. [Google Scholar] [CrossRef]
  6. Xiong, C.; Li, Y. Error analysis for optimal control problem governed by convection diffusion equations: DG method. J. Comput. Appl. Math. 2011, 235, 3163–3177. [Google Scholar] [CrossRef] [Green Version]
  7. Xiong, C.; Luo, F.; Ma, X.; Li, Y. A priori error analysis for optimal distributed control problem governed by the first order linear hyperbolic equation: hp-streamline diffusion discontinuous galerkin method. J. Numer. Math. 2016, 24, 125–134. [Google Scholar] [CrossRef]
  8. Xiong, C.; Li, Y. A posteriori error estimators for optimal distributed control governed by the first-order linear hyperbolic equation: DG method. NMPDEs 2011, 27, 491–506. [Google Scholar] [CrossRef]
  9. Xiong, C.; Becker, R.; Luo, F.; Ma, X. A priori and a posteriori error analysis for the mixed discontinuous Galerkin finite element approximations of the biharmonic problems. NMPDEs 2017, 33, 318–353. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, L.; Xiong, C.; Wu, H.; Luo, F. A priori and a posteriori error analysis for discontinuous Galerkin finite element approximations of biharmonic eigenvalue problems. Adv. Comput. Math. 2013, 45, 2623–2646. [Google Scholar] [CrossRef]
  11. Dedner, A.; Madhavan, P. Discontinuous Galerkin Methods for Hyperbolic and Advection-Dominated Problems on Surfaces. 2015. Available online: https://arxiv.org/pdf/1505.06752.pdf (accessed on 20 December 2020).
  12. Hughes, T.J.R.; Cottrell, J.A.; Bazilevs, T. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Comput. Methods Appl. Mech. Engrg. 2005, 194, 4135–4195. [Google Scholar] [CrossRef] [Green Version]
  13. Langer, U.; Toulopoulos, I. Analysis of multipatch discontinuous Galerkin IgA approximations to elliptic boundary value problems. Comput. Vis. Sci. 2015, 17, 217–233. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, F.; Xu, Y.; Chen, F. Discontinuous Galerkin methods for isogeometric analysis for elliptic equations on surfaces. Commun. Math. Stat. 2014, 2, 431–461. [Google Scholar] [CrossRef] [Green Version]
  15. Hofer, C.; Langer, U.; Toulopoulos, I. Discontinuous Galerkin Isogeometric Analysis of Elliptic Diffusion Problems on Segmentations with Gaps. Comput. Math. Appl. 2015, 72, 1811–1827. [Google Scholar] [CrossRef]
  16. Bazilevs, Y.; Hughes, T.J.R. NURBS-based isogeometric analysis for the computation of flows about rotating components. Comput. Mech. 2008, 43, 143–150. [Google Scholar] [CrossRef]
  17. Elfverson, D.; Larson, M.G.; Larsso, K. A new least squares stabilized Nitsche method for cut isogeometric analysis. Comput. Methods Appl. Mech. Eng. 2019, 349, 1–16. [Google Scholar] [CrossRef] [Green Version]
  18. Schumaker, L.L. Spline Functions: Basic Theory, 3rd ed.; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  19. Brezzi, F.; Marini, L.D.; Süli, E. Discontinuous Galerkin methods for first-order hyperbolic problems. Math. Model. Methods Appl. Sci. 2004, 14, 1893–1903. [Google Scholar] [CrossRef]
Figure 1. Sphere Case. (Left): The patches and meshes; (Right): Numerical solution.
Figure 1. Sphere Case. (Left): The patches and meshes; (Right): Numerical solution.
Mathematics 09 00497 g001
Figure 2. Error convergence result for sphere case.
Figure 2. Error convergence result for sphere case.
Mathematics 09 00497 g002
Figure 3. Torus Case. (Left): The patches and meshes; (Right): Numerical solution.
Figure 3. Torus Case. (Left): The patches and meshes; (Right): Numerical solution.
Mathematics 09 00497 g003
Figure 4. Error convergence result for torus case.
Figure 4. Error convergence result for torus case.
Mathematics 09 00497 g004
Figure 5. Cylinder Case. (Left): The patches and meshes; (Right): Numerical solution.
Figure 5. Cylinder Case. (Left): The patches and meshes; (Right): Numerical solution.
Mathematics 09 00497 g005
Figure 6. Error convergence result for cylinder case.
Figure 6. Error convergence result for cylinder case.
Mathematics 09 00497 g006
Figure 7. Cylinder Case. (Left): The patches and meshes; (Right): Numerical solution.
Figure 7. Cylinder Case. (Left): The patches and meshes; (Right): Numerical solution.
Mathematics 09 00497 g007
Figure 8. Error convergence result for cylinder case.
Figure 8. Error convergence result for cylinder case.
Mathematics 09 00497 g008
Table 1. Errors and convergence order of sphere case for k = 1 , 2 and 3.
Table 1. Errors and convergence order of sphere case for k = 1 , 2 and 3.
Degree kMesh Size hErrorConvergence Order
11/29.934 × 10 2
1/41.926 × 10 2 2.366
1/84.561 × 10 3 2.078
1/161.128 × 10 3 2.016
1/322.813 × 10 4 2.004
21/21.580 × 10 2
1/41.656 × 10 3 3.255
1/81.799 × 10 4 3.202
1/162.202 × 10 5 3.031
1/322.748 × 10 6 3.002
31/22.840 × 10 3
1/43.006 × 10 4 3.240
1/81.210 × 10 5 4.635
1/166.826 × 10 7 4.148
1/324.150 × 10 8 4.040
Table 2. Errors and convergence orders for the torus for k = 1 , 2 and 3.
Table 2. Errors and convergence orders for the torus for k = 1 , 2 and 3.
Degree kMesh Size hErrorConvergence Order
11/22.2825 × 10 0
1/43.5096 × 10 1 2.7013
1/87.2088 × 10 2 2.2835
1/161.7199 × 10 2 2.0674
1/324.2620 × 10 3 2.0127
21/29.4780 × 10 1
1/41.2046 × 10 1 2.9760
1/89.5761 × 10 3 3.6530
1/161.1486 × 10 3 3.0595
1/321.4448 × 10 5 2.9910
31/20.3400 × 10 1
1/45.4673 × 10 2 2.6368
1/89.7914 × 10 4 5.8032
1/164.6835 × 10 5 4.3859
1/322.7729 × 10 6 4.0781
Table 3. Errors and convergence order of cylinder case for k = 1 , 2 and 3.
Table 3. Errors and convergence order of cylinder case for k = 1 , 2 and 3.
Degree kMesh Size hErrorConvergence Order
11/27.832 × 10 3
1/41.912 × 10 3 2.0344
1/84.748 × 10 4 2.0097
1/161.185 × 10 4 2.0023
1/322.962 × 10 5 2.0004
21/24.5649 × 10 4
1/45.7310 × 10 5 2.9937
1/87.2113 × 10 6 2.9905
1/169.0619 × 10 7 2.9924
1/321.1364 × 10 7 2.9954
31/21.7089 × 10 5
1/41.5120 × 10 6 3.4986
1/81.0073 × 10 7 3.9079
1/166.4755 × 10 9 3.9594
1/324.1048 × 10 10 3.9796
Table 4. Errors and convergence orders for the cylinder for k = 1 , 2 and 3.
Table 4. Errors and convergence orders for the cylinder for k = 1 , 2 and 3.
Degree kMesh Size hErrorConvergence Order
11/21.774 × 10 2
1/44.277 × 10 3 2.0531
1/81.063 × 10 3 2.0085
1/162.656 × 10 4 2.0006
1/326.640 × 10 5 2.0002
21/21.7036 × 10 3
1/41.9959 × 10 4 3.0935
1/82.4770 × 10 5 3.0104
1/163.1024 × 10 6 2.9971
1/323.8885 × 10 7 2.9961
31/21.0065 × 10 4
1/48.1947 × 10 6 3.6186
1/85.3064 × 10 7 3.9489
1/163.3827 × 10 8 3.9715
1/322.1381 × 10 9 3.9838
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, L.; Xiong, C.; Yuan, X.; Wu, H. Discontinuous Galerkin Isogeometric Analysis of Convection Problem on Surface. Mathematics 2021, 9, 497. https://0-doi-org.brum.beds.ac.uk/10.3390/math9050497

AMA Style

Wang L, Xiong C, Yuan X, Wu H. Discontinuous Galerkin Isogeometric Analysis of Convection Problem on Surface. Mathematics. 2021; 9(5):497. https://0-doi-org.brum.beds.ac.uk/10.3390/math9050497

Chicago/Turabian Style

Wang, Liang, Chunguang Xiong, Xinpeng Yuan, and Huibin Wu. 2021. "Discontinuous Galerkin Isogeometric Analysis of Convection Problem on Surface" Mathematics 9, no. 5: 497. https://0-doi-org.brum.beds.ac.uk/10.3390/math9050497

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop