Processing math: 100%
Research article Special Issues

Assessment of health status of tree trunks using ground penetrating radar tomography

  • Typical tree health assessment methods are destructive. Non-destructive Ground Penetrating Radar (GPR) technique can provide a diagnostic tool for assessing the health status of live tree trunks based on internal dielectric permittivity distribution. Typical GPR acquisition technique—common-offset—is not helpful in providing robust and high-resolution quantitative results. In the current work we evaluate the capabilities of GPR tomography on locating tree-decays in a number of different tree species, imaging the interval structure of a healthy tree and quantitative estimation of moisture content (MC) based on distribution of dielectric permittivity, directly related to MC. The measurements described in this work were made on the trunks of live trees of different species in different conditions: a "healthy" English oak (Quercus robur), a "dry" Siberian fir (Picea obovata), a Horse chestnut (Castanea dentata) and a European aspen (Populus tremula) with rot inside. The results of the suggested approach were confirmed by resistography. Different parts of the trunk (bark, core, sapwood), as well as healthy and affected areas differ in moisture content, so the method of GPR tomography allowed us to see both the structure of the trunk of a healthy tree, and the presence and dimensions of defects.

    Citation: Maria Sudakova, Evgeniya Terentyeva, Alexey Kalashnikov. Assessment of health status of tree trunks using ground penetrating radar tomography[J]. AIMS Geosciences, 2021, 7(2): 162-179. doi: 10.3934/geosci.2021010

    Related Papers:

    [1] Yue Feng, Yujie Liu, Ruishu Wang, Shangyou Zhang . A conforming discontinuous Galerkin finite element method on rectangular partitions. Electronic Research Archive, 2021, 29(3): 2375-2389. doi: 10.3934/era.2020120
    [2] Chunmei Wang . Simplified weak Galerkin finite element methods for biharmonic equations on non-convex polytopal meshes. Electronic Research Archive, 2025, 33(3): 1523-1540. doi: 10.3934/era.2025072
    [3] Leilei Wei, Xiaojing Wei, Bo Tang . Numerical analysis of variable-order fractional KdV-Burgers-Kuramoto equation. Electronic Research Archive, 2022, 30(4): 1263-1281. doi: 10.3934/era.2022066
    [4] Guanrong Li, Yanping Chen, Yunqing Huang . A hybridized weak Galerkin finite element scheme for general second-order elliptic problems. Electronic Research Archive, 2020, 28(2): 821-836. doi: 10.3934/era.2020042
    [5] Victor Ginting . An adjoint-based a posteriori analysis of numerical approximation of Richards equation. Electronic Research Archive, 2021, 29(5): 3405-3427. doi: 10.3934/era.2021045
    [6] Jun Pan, Yuelong Tang . Two-grid H1-Galerkin mixed finite elements combined with L1 scheme for nonlinear time fractional parabolic equations. Electronic Research Archive, 2023, 31(12): 7207-7223. doi: 10.3934/era.2023365
    [7] Hongze Zhu, Chenguang Zhou, Nana Sun . A weak Galerkin method for nonlinear stochastic parabolic partial differential equations with additive noise. Electronic Research Archive, 2022, 30(6): 2321-2334. doi: 10.3934/era.2022118
    [8] Zexuan Liu, Zhiyuan Sun, Jerry Zhijian Yang . A numerical study of superconvergence of the discontinuous Galerkin method by patch reconstruction. Electronic Research Archive, 2020, 28(4): 1487-1501. doi: 10.3934/era.2020078
    [9] Suayip Toprakseven, Seza Dinibutun . A weak Galerkin finite element method for parabolic singularly perturbed convection-diffusion equations on layer-adapted meshes. Electronic Research Archive, 2024, 32(8): 5033-5066. doi: 10.3934/era.2024232
    [10] Bin Wang, Lin Mu . Viscosity robust weak Galerkin finite element methods for Stokes problems. Electronic Research Archive, 2021, 29(1): 1881-1895. doi: 10.3934/era.2020096
  • Typical tree health assessment methods are destructive. Non-destructive Ground Penetrating Radar (GPR) technique can provide a diagnostic tool for assessing the health status of live tree trunks based on internal dielectric permittivity distribution. Typical GPR acquisition technique—common-offset—is not helpful in providing robust and high-resolution quantitative results. In the current work we evaluate the capabilities of GPR tomography on locating tree-decays in a number of different tree species, imaging the interval structure of a healthy tree and quantitative estimation of moisture content (MC) based on distribution of dielectric permittivity, directly related to MC. The measurements described in this work were made on the trunks of live trees of different species in different conditions: a "healthy" English oak (Quercus robur), a "dry" Siberian fir (Picea obovata), a Horse chestnut (Castanea dentata) and a European aspen (Populus tremula) with rot inside. The results of the suggested approach were confirmed by resistography. Different parts of the trunk (bark, core, sapwood), as well as healthy and affected areas differ in moisture content, so the method of GPR tomography allowed us to see both the structure of the trunk of a healthy tree, and the presence and dimensions of defects.



    For simplicity, we consider Poisson equation with a Dirichlet boundary condition as our model problem.

    Δu=f,inΩ, (1)
    u=g,onΩ, (2)

    where Ω is a bounded polygonal domain in R2.

    Using integration by parts, we can get the variational form: find uH1(Ω) satisfying u=gonΩ and

    (u,v)=(f,v),vH10(Ω). (3)

    Various finite element methods have been introduced to solve the Poisson equations (1)-(2), such as the Galerkin finite element methods (FEMs)[2, 3], the mixed FEMs [15] and the finite volume methods (FVMs) [6], etc. The FVMs emphasis on the local conservation property and discretize equations by asking the solution satisfying the flux conservation on a dual mesh consisting of control volumes. The mixed FEMs is another category method that based on the variable u and a flux variable usually written as p.

    The classical conforming finite element method obtains numerical approximate results by constructing a finite-dimensional subspace of H10(Ω). The finite element scheme has the same form with the variational form (3): find uhVhH1(Ω) satisfying uh=IhgonΩ and

    (uh,vh)=(f,vh),vhV0h, (4)

    where V0h is a subspace of Vh that satisfying vh=0 on Ω and Ih is the kth order Lagrange interpolation operator. The FE method is a popular and easy-to-implement numerical scheme, however, it is less flexible in constructing elements and generating meshes. These limitations are mainly due to the strong continuity requirements of functions in Vh. One solution to circumvent these limitations is using discontinuous approximations. Since the 1970th, many new finite element methods with discontinuous approximations have been developed, including the early proposed DG methods [1], local discontinuous Galerkin (LDG) methods [8], interior penalty discontinuous Galerkin (IPDG) methods [9], and the recently developed hybridizable discontinuous Galerkin (HDG) methods [7], mimetic finite differences method [10], virtual element (VE) method [4], weak Galerkin (WG) method [19, 20] and references therein.

    One obvious disadvantage of discontinuous finite element methods is their rather complex formulations which are often necessary to ensure connections of discontinuous solutions across element boundaries. For example, the IPDG methods add parameter depending interior penalty terms. Besides additional programming complexity, one often has difficulties in finding optimal values for the penalty parameters and corresponding efficient solvers. Most recently, Zhang and Ye [21] developed a discontinuous finite element method that has an ultra simple weak formulation on triangular/tetrahedal meshes. The corresponding numerical scheme can be written as: find uh˜Vh satisfying uh=IhgonΩ and

    (wuh,wvh)=(f,vh),vhV0h, (5)

    where ˜Vh is the DG finite element space and w is the weak gradient operator. The notion of weak gradient was first introduced by Wang and Ye in the weak Galerkin (WG) methods [19, 20]. The WG methods allow the use of totally discontinuous functions and provides stable numerical schemes that are parameter-independent and free of locking [17] in some applications. Another key feature in the WG methods is it can be used for arbitrary polygonal meshes. The WG finite element method has been rapidly developed and applied to other problems, including the Stokes and Navier-Stokes equations [11, 18], the biharmonic [14, 13] and elasticity equations [12, 17], div-curl systems and the Maxwell's equations and parabolic problem [23], etc. The introduction of the weak gradient operator in the conforming DG methods makes the scheme (5) maintain the simple formulation of conforming finite element method while have the flexibility of using discontinuous approximations. Hence, the programming complexity of this conforming DG scheme is significantly reduced. Furthermore, the scheme results in a simple symmetric and positive definite system.

    Following the work in [21, 22], we propose a new conforming DG finite element method on rectangular partitions in this work. It can be obtained from the conforming formulation simply by replacing by w and enforcing the boundary condition strongly. The simplicity of the conforming DG formulation will ease the complexity for implementation of DG methods. We note that the conforming DG method in [21] is based on triangular/tetrahedal meshes. Then in [22], the method is extended to work on general polytopal meshes by raising the degree of polynomials used to compute weak gradient.

    In this paper, we keep the same finite element space as DG method, replace the boundary function with the average of the inner function, and use the weak gradient arising from local Raviart-Thomas (RT) elements [5] to approximate the classic gradient. Moreover, the derivation process in this paper is based on rectangular RT elements [16]. Error estimates of optimal order are established for the corresponding conforming DG approximation in both a discrete H1 norm and the L2 norm. Numerical verifications have been performed on different kinds of quadrangle finite element space. In particular, super-convergence phenomenon have been observed for Q0 elements.

    The rest of this paper is organized as follows: In Section 2, we shall present the conforming DG finite element scheme for the Poisson equation on rectangular partitions. Section 3 is devoted to a discussion of the stability and solvability of the new method. In Section 4, we shall prepare ourselves for error estimates by deriving some identities. Error estimates of optimal order in H1 and L2 norm are established in Section 5. In Section 6, we present some numerical results to illustrate the theory derived in earlier sections. Finally in section 7, we conclude our major contributions in this article.

    Throughout this paper, we adopt the standard definition of Sobolev space Hs(Ω). For any given open bounded domain KΩ, (,)s,K,s,K, and ||s,K are used to denote the inner product, norm and semi-norm, respectively. The space H0(K) coincides with L2(K), and the subscripts K in the inner product, norm, and semi-norm can be dropped in the case of K=Ω. In particular, the function space H10(Ω) is defined as

    H10(Ω)={vH1(Ω):v|Ω=0},

    and the space H(div,Ω) is defined as the set of vector-valued functions q, which together with their divergence are square integrable, i.e.

    H(div,Ω)={q[L2(Ω)]d:qL2(Ω)}.

    Assume that the domain Ω is of polygonal type and is partitioned into non-overlapping rectangles Th={T}. For each TTh, denote by T0 its interior and T its boundary. Denote by Eh={e} the set of all edges in Th, and E0h=EhΩ the set of all interior edges in Th. For each TTh and eEh, denote by hT and he the diameter of T and e, respectively. h=maxTThhT is the meshsize of Th.

    For any interior edge eE0h, let T1 and T2 be two rectangles sharing e, we define the average {} and the jump [[]] on e for a scalar-valued function v by

    {v}=12(v|T1+v|T2),[[v]]=v|T1n1+v|T2n2, (6)

    where v|Ti,i=1,2 is the trace of v on Ti, n</italic>1 and n</italic>2 are the two unit outward normal vectors on e, associated with T1 and T2, respectively. If e is a boundary edge, we define

    {v}=v|eand[[v]]=v|en. (7)

    We define a discontinuous finite element space

    Vh={vL2(Ω):v|TQk(T),TTh}, (8)

    and its subspace

    V0h={vVh:v=0onΩ}, (9)

    where Qk(T),k1 denotes the set of polynomials with regard to quadrilateral elements. The weak gradient for a scalar-valued function vVh is defined by the following definition

    Definition 2.1. For a given TTh and a function vVh, the discrete weak gradient dvRTk(T) on T is defined as the unique polynomial such that

    (dv,q)T:=(v,q)T+{v},qnT,qRTk(T), (10)

    where n is the unit outward normal on T, RTk(T)=[Qk(T)]2+xQk(T), and {v} is defined in (6) and (7).

    The weak gradient operator d as defined in (10) is a local operator computed at each element. It can be extended to any function vVh by taking weak gradient locally on each element T. More precisely, the weak gradient of any vVh is defined element-by-element as follows:

    (dv)|T=d(v|T).

    We introduce the following bilinear form:

    a(v,w)=(dv,dw),

    the conforming DG algorithm to solve the problems (1) - (2) is given by

    Conforming DG algorithm 1. Find uhVh satisfying uh=IhgonΩ and

    a(uh,vh)=(f,vh),vhV0h, (11)

    where Ih is the kth order Lagrange interpolation.

    We will prove the existence and uniqueness of the solution of equation (11). Firstly, we present the following two useful inequalities to derive the forthcoming analysis.

    Lemma 3.1 (trace inequality). Let T be an element of the finite element partition Th, and e is an edge or face which is part of T. For any function φH1(T), the following trace inequality holds true (see [20] for details):

    φ2eC(h1Tφ2T+hTφ2T), (12)

    where C is a constant independent of h.

    Lemma 3.2 (inverse inequality). Let Th be a finite element partition of Ω that is shape regualr. Assume that Th satisfies all the assumptions A1-A4 in [20]. Then, for any piecewise polynomial function φ of degree n on Th, there exists a constant C=C(n) such that

    φTC(n)h1TφT,TTh. (13)

    Then, we define the following semi-norms in the discontinuous finite element space Vh

    |||v|||2=a(v,v)=TThdv2T, (14)
    v21,h=TThv2T+eE0hh1e[[v]]2e. (15)

    We have the equivalence between the semi-norms |||v||| and v1,h, and it is proved in the following lemma.

    Lemma 3.3. For any vVh, the following equivalence holds true

    C1v1,h|||v|||C2v1,h, (16)

    where C1 and C2 are two constants independent of h.

    Proof. It follows from the definition of dv, integration by parts, the trace inequality, and the inverse inequality that

    dv2T1=(dv,dv)T1=(v,dv)T1+{v}n,dvT1=(v,dv)T1(v{v})n,dvT1vT1dvT1+(v{v})nT1dvT1dvT1(vT1+h12T1(v{v})nT1). (17)

    For any eT1, e=T1T2, we have

    (v{v})|en1=v|T1n112(v|T1+v|T2)n1=12(v|T1n1+v|T2n2)=12[[v]]e.

    Then we can get

    (v{v})n2T112eT1[[v]]2e. (18)

    Substituting (18) into (17) gives

    dv2T1C2dvT1(vT1+eT1h12e[[v]]e),

    this completes the proof of the right-hand of (16).

    To prove the left-hand of (16), we consider the subspace of RTk(T) for any TTh

    D(k,T):={qRTk(T):qn=0onT}.

    Note that D(k,T) is a dual space of [Qk1(T)]2 [13]. Thus, for any v[Qk1(T)]2, we have

    vT=supqD(k,T)(v,q)TqT. (19)

    Using the integration by parts, Cauchy-Schwarz inequality, the definition of D(k,T) and dv, we get

    (v,q)T=(v,q)T+v,qnT=(dv,q)T{v},qnT=(dv,q)TdvTqT,

    where we have used the fact that q</italic>n|T=0 in the definition of D(k,T). Combining the above result with (19), one has

    vTdvT. (20)

    We define the space De(k,T) as the set of all q</italic>RTk(T) such that all degrees of freedom, except those for q</italic>n|e, vanish. Note that De(k,T) is a dual space of [Qk(e)]2 [13]. Thus, we know

    [[v]]e=supqDe(k,T)[[v]],qneqne. (21)

    Following the integration by parts and the definition of d, we can derive that

    (dv,q)T=(v,q)Tv,qne+{v},qne.

    Together with (20), we obtain

    |[[v]],qne|=2|(dv,q)T(v,q)T|2|(dv,q)T|+2|(v,q)T|C(dvTqT+vTqT)CdvTqT.

    Substituting the above inequality into (21), by the scaling argument [13], for such qDe(k,T), we have qTh12qne, then

    [[v]]eCdvTqTqneCh12dvT. (22)

    Combining (20) and (22) gives a proof of the left-hand of (16).

    Lemma 3.4. The semi-norm |||||| defined in (14) is a norm in V0h.

    Proof. We shall only verify the positivity property for ||||||. To this end, assume |||v|||=0 for some vV0h. By Lemma 3.3, it follows that v1,h=0 for all TTh, which means that v=0 for all elements TTh and [[v]]=0 for all edges eE0h. We can derive from v=0 for all TTh that v is a constant in each T. [[v]]=0 on each eE0h implies v is a continuous function. This two conclusions and v=0 on Ω show that v=0, which completes the proof of the lemma.

    The above two lemmas imply the well posedness of the scheme (11). We prove the existence and uniqueness of solution of the conforming DG method in Theorem 3.1.

    Theorem 3.1. The conforming DG scheme (11) has and only has one solution.

    Proof. To prove the scheme (11) is uniquely solvable, it suffices to verify that the homogeneous equation has zero as its unique solution. To this end, let uhVh be the solution of the numerical scheme 11 with homogeneous data f=0, g=0. Letting vh=uh, we obtain

    a(uh,uh)=0,

    which leads to uh=0 by using Lemma 3.4. This completes the proof of the theorem.

    In this section, we will derive an error equation which will be used for the error estimates. For any q</italic>H(div,Ω), we assume that there exist an interpolation operator Πh satisfying ΠhqH(div,Ω)RTk(T) on each element TTh and

    (q,v)T=(Πhq,v)T,vQk(T). (23)

    For any wH1+k(Ω) with k1, from Lemma 7.3 in [20], we have the estimate of Πh as follows.

    Πh(w)wChkw1+k. (24)

    Moreover, it is easy to verify the following property holds true.

    Lemma 4.1. For any qH(div,Ω),

    TTh(q,v)T=TTh(Πhq,dv)T,vV0h. (25)

    Proof. ΠhqH(div,Ω) implies that Πhq is continuous across each interior edge. Since vV0h, we know that {v}=v=0 on Ω. Then

    TTh{v},ΠhqnT=0. (26)

    By the definition of Πh and d and the equation (26), we have

    TTh(q,v)T=TTh(Πhq,v)T=TTh(Πhq,v)T+TTh{v},ΠhqnT=TTh(Πhq,dv)T.

    This completes the proof of the lemma.

    Before establishing the error equation, we define a continuous finite element subspace of Vh as follows

    ˜Vh={vH1(Ω):v|TQk(T),TTh}. (27)

    so as a subspace of ˜Vh

    ˜V0h:={v˜Vh:v|Ω=0}. (28)

    Lemma 4.2. For any v˜Vh, we have

    dv=v.

    Proof. By the definition of d and integration by parts, for any q</italic>RTk(T), we have

    (dv,q)T=(v,q)T+{v},qnT=(v,q)T+v,qnT=(v,q)T,

    which gives

    (dvv,q)T=0,qRTk(T).

    Letting q be dvv in the above equation yields dvv=0, which completes the proof of the lemma.

    Let eh=Ihuuh, where Ih is the kth order Lagrange interpolation, uHk+1(Ω) with k1 is the exact solution of the Poisson equations (1) - (2), and uhVh is the numerical solution of the scheme (11). The following estimate of the Lagrange interpolation operator Ih holds true.

    IhuuChk+1uk+1, (29)
    IhuuChkuk+1. (30)

    It is obvious that ehV0h and Ihu˜Vh. We have the following lemma:

    Lemma 4.3. Denote eh=Ihuuh the error of conforming DG method arising from (11). For any vhV0h, we have

    a(eh,vh)=lu(vh), (31)

    where

    lu(vh)=TTh(IhuΠhu,dvh). (32)

    Proof. Since Ihu˜Vh, we have dIhu=Ihu. Using the property (25), we can derive

    TTh(dIhu,dvh)T=TTh(Ihu,dvh)T=TTh(IhuΠhu+Πhu,dvh)T=TTh(IhuΠhu,dvh)T+TTh(Πhu,dvh)T=lu(vh)TTh(u,vh)T=lu(vh)+(f,vh).

    By the definition of the scheme (11), we have

    TTh(dIhuduh,dvh)T=lu(vh).

    This completes the proof of the lemma.

    The goal of this section is to derive the error estimates in H1 and L2 norms for the conforming DG solution uh.

    Theorem 5.1. Let uHk+1(Ω) with k1 be the exact solution of the Poisson equation (1) - (2), and uhVh be the numerical solution of the scheme (11). Let eh=Ihuuh, there exists a constant C independent of h such that

    |||eh|||Chk|u|k+1. (33)

    Proof. Letting vh=eh in (31), and by the definition of ||||||, we have

    |||eh|||2=lu(eh). (34)

    From the Cauchy-Schwarz inequality, the triangle inequality, the definition of ||||||, (24), and (30), we arrive at

    lu(vh)=TTh(IhuΠh(u),dvh)TTThIhuΠh(u)TdvhT(TThIhuΠh(u)2T)12(TThdvh2T)12=(TThIhuu+uΠh(u)2T)12|||vh|||(TThIhuu2T+uΠh(u)2T)12|||vh|||Chk|u|k+1|||vh|||.

    Then, we have

    lu(eh)Chk|u|k+1|||eh|||. (35)

    Substituting (35) to (34), we obtain

    |||eh|||2Chk|u|k+1|||eh|||,

    which completes the proof of the lemma.

    It is obvious that ˜V0hV0h. Let ˜uh˜Vh be the finite element solution for the problem (1)-(2) which satisfies ˜uh=Ihg on Ω and

    (˜uh,v)=(f,v),v˜V0h. (36)

    For any v˜V0h˜Vh, we have dv=v, i.e.

    (duh˜uh,v)=0,v˜V0h. (37)

    In the rest of this section, we derive an optimal order error estimate for the conforming DG approximation (11) in L2 norm by adopting the duality argument. To this end, we consider the following dual problem that seeks ΦH10(Ω) satisfying

    (Φ)=uh˜uh,inΩ. (38)

    Assume that the dual problem satisfies H2-regularity, which means the following priori estimate holds true

    Φ2Cuh˜uh. (39)

    In the following of this paper, we note εh=uh˜uh for simplicity.

    Theorem 5.2. Assume uHk+1(Ω) with k1 is the exact solution of the Poisson equation (1) - (2), and uhVh is the numerical solution obtained with the scheme (11). Furthermore, assume that (39) holds true. Then, there exists a constant C independent of h such that

    uuhChk+1|u|k+1. (40)

    Proof. First, we shall derive the optimal order for εh in L2 norm. Consider the corresponding conforming DG scheme defined in (11) and let ΦhV0h be the solution satisfying

    a(Φh,v)=(εh,v),vV0h. (41)

    Since IhΦ˜Vh, it follows from (37) that

    (duh˜uh,IhΦ)=0,dIhΦ=IhΦ,

    which gives

    (duh˜uh,dIhΦ)=0. (42)

    Setting v=εh in (41), then by the definition of εh and (42), we have

    εh2=a(Φh,εh)=TTh(dΦh,dεh)T=TTh(d(ΦhIhΦ),duh˜uh)T|||ΦhIhΦ|||(|||uhIhu|||+(Ihu˜uh)).

    Then, by the Cauchy-Schwarz inequality, (33) and (39), we obtain

    εh2Ch|Φ|2hk|u|k+1Chk+1|u|k+1εh,

    which gives

    εhChk+1|u|k+1. (43)

    Combining the error estimate of finite element solution, the triangle inequality and (43) yields (40), which completes the proof of the theorem.

    In this section, we shall present some numerical results for the conforming discontinuous Galerkin method analyzed in the previous sections.

    We solve the following Poisson equation on the unit square domain Ω=(0,1)×(0,1),

    Δu=2π2sin(πx)sin(πy)in Ω (44)
    u=0on Ω. (45)

    The exact solution of the above problem is u=sin(πx)sin(πy). Uniform square grids as shown in Figure 1 are used for computation.

    Figure 1.  The first three grids used in the computation.

    We first use the Pk conforming discontinuous Galerkin spaces (8) to compute the test case (44)-(45), where Pk denotes the set of polynomials of 2 variables of degree less than or equal to k. The weak gradient is computed locally using rectangular RTk polynomials. The errors and the order of convergence of the conforming DG approximations are listed in Table 1. Optimal order of convergence is achieved in every case, which is consistent with our theory. In particular, a superconvergence of order O(h2) was observed in the discrete H1 norm for P0 elements. Furthermore, the results obtained with P0 elements seems to be slightly better than that obtained with P1 elements.

    Table 1.  Error profiles and convergence rates for test case (44)-(45) obtained with uniform grids and Pk conforming DG spaces.
    level uhQhu0 rate |||uhQhu||| rate #Dof
    by P0 conforming discontinuous Galerkin elements
    6 0.1996E-02 1.97 0.8887E-02 1.98 1024
    7 0.5013E-03 1.99 0.2228E-02 2.00 4096
    8 0.1255E-03 2.00 0.5574E-03 2.00 16384
    by P1 conforming discontinuous Galerkin elements
    6 0.2427E-02 1.97 0.1027E+00 1.02 3072
    7 0.6100E-03 1.99 0.5105E-01 1.01 12288
    8 0.1527E-03 2.00 0.2546E-01 1.00 49152
    by P2 conforming discontinuous Galerkin elements
    5 0.1533E-03 3.00 0.2042E-01 2.03 1536
    6 0.1915E-04 3.00 0.5061E-02 2.01 6144
    7 0.2394E-05 3.00 0.1260E-02 2.01 24576
    by P3 conforming discontinuous Galerkin elements
    5 0.7959E-05 4.00 0.1965E-02 3.00 2560
    6 0.4971E-06 4.00 0.2451E-03 3.00 10240
    7 0.3140E-07 3.98 0.3059E-04 3.00 40960
    by P4 conforming discontinuous Galerkin elements
    4 0.1055E-04 4.97 0.1421E-02 4.05 960
    5 0.3314E-06 4.99 0.8735E-04 4.02 3840
    6 0.1057E-07 4.97 0.5417E-05 4.01 15360
    by P5 conforming discontinuous Galerkin elements
    2 0.2835E-02 6.24 0.1450E+00 5.49 84
    3 0.4532E-04 5.97 0.4718E-02 4.94 336
    4 0.7115E-06 5.99 0.1478E-03 5.00 1344

     | Show Table
    DownLoad: CSV

    The same test case is also computed using the Qk conforming DG finite element space, where Qk denotes the set of polynomials of 2 variables defined on Ω, and for each variable, the degree of the variable is at most k. Table 2 illustrates the numerical performance of the corresponding conforming DG scheme. It can be seen from numerical computing that, in this case, the results obtained with the Q1 element are more accurate than those obtained with Q0(=P0) elements (see Table 1). All numerical results converge at the corresponding optimal order, which is consistent with the theory.

    Table 2.  Error profiles and convergence rates for test case (44)-(45) obtained with uniform grids and Qk conforming DG spaces.
    level uhQhu0 rate |||uhQhu||| rate #Dof
    by Q1 conforming discontinuous Galerkin elements
    6 0.4006E-03 1.99 0.2389E-02 1.99 4096
    7 0.1003E-03 2.00 0.5982E-03 2.00 16384
    8 0.2510E-04 2.00 0.1496E-03 2.00 65536
    by Q2 conforming discontinuous Galerkin elements
    6 0.2360E-04 2.99 0.3186E-02 1.99 9216
    7 0.2953E-05 3.00 0.7976E-03 2.00 36864
    8 0.3692E-06 3.00 0.1995E-03 2.00 147456
    by Q3 conforming discontinuous Galerkin elements
    5 0.1413E-04 4.08 0.1650E-02 2.97 4096
    6 0.8676E-06 4.03 0.2072E-03 2.99 16384
    7 0.5398E-07 4.01 0.2593E-04 3.00 65536
    by Q4 conforming discontinuous Galerkin elements
    3 0.2226E-02 4.59 0.5414E-01 3.52 400
    4 0.9610E-04 4.53 0.3723E-02 3.86 1600
    5 0.3279E-05 4.87 0.2392E-03 3.96 6400

     | Show Table
    DownLoad: CSV

    To test the superconvergence of P0 DG element, we solve the following 2nd order elliptic equation on the unit square domain Ω=(0,1)×(0,1),

    Δu+u=fin Ωu=0on Ω,

    where f is chosen so that the exact solution is not symmetric,

    u=(xx2)(yy3). (46)

    Uniform square grids as shown in Figure 1 are used for numerical computation. The numerical results are listed in Table 3. Surprising, for this problem, the H1-like norm of error superconverges at 1.5 order, and the L2 error has one order of superconvergence. But we do not yet know if such a superconvergence exists in general.

    Table 3.  Error profiles and convergence rates for test case (46) obtained with uniform grids and P0 conforming DG spaces.
    level uhQhu0 rate |||uhQhu||| rate #Dof
    by P0 conforming discontinuous Galerkin elements
    3 0.8265E-02 1.06 0.4577E-01 1.14 16
    4 0.2772E-02 1.58 0.1732E-01 1.40 64
    5 0.7965E-03 1.80 0.6331E-02 1.45 256
    6 0.2142E-03 1.90 0.2290E-02 1.47 1024
    7 0.5564E-04 1.94 0.8213E-03 1.48 4096
    8 0.1419E-04 1.97 0.2928E-03 1.49 16384

     | Show Table
    DownLoad: CSV

    To test further the superconvergence of P0 DG element, we solve the following 2nd order elliptic equations on the unit square domain Ω=(0,1)×(0,1),

    (au)=fin Ωu=0on Ω,

    where a=1+x+y and f is chosen so that the exact solution is not symmetric,

    u=(xx3)(y2y3). (47)

    Uniform square grids as shown in Figure 1 are used for computation. The numerical results are listed in Table 4. Surprising, again, the H1-like norm of error superconverges at 1.5 order, and the L2 error has one order of superconvergence for this problem.

    Table 4.  Error profiles and convergence rates for test case (47) obtained with uniform grids and P0 conforming DG spaces.
    level uhQhu0 rate |||uhQhu||| rate #Dof
    by P0 conforming discontinuous Galerkin elements
    3 0.4929E-02 0.97 0.5371E-01 0.80 16
    4 0.1917E-02 1.36 0.2401E-01 1.16 64
    5 0.6004E-03 1.67 0.9407E-02 1.35 256
    6 0.1682E-03 1.84 0.3507E-02 1.42 1024
    7 0.4457E-04 1.92 0.1275E-02 1.46 4096
    8 0.1148E-04 1.96 0.4576E-03 1.48 16384

     | Show Table
    DownLoad: CSV

    In this paper, we establish a new numerical approximation scheme based on the rectangular partition to solve second order elliptic equation. We derived the numerical scheme and then proved the optimal order of convergence of the error estimates in L2 and H1 norms of the conforming DG method. Numerical experiments are then present to verify the theoretical analysis, and all numerical results converging at the corresponding optimal order. Comparing with existing numerical methods, the confoming DG method has the following two characteristics: 1. The formulation is relatively simple. The stabilizer s(,) is no longer needed, and the boundary function ub is omitted, which is replaced by the average of internal function u0; 2. The projection operator Qh used in the traditional WG method is replaced by the Lagrange interpolation operator Ih, which makes the theoretical analysis much easier. As can be seen from the numerical examples in Section 6, this method reduces the programming complexity while ensuring the optimal order of convergence.



    [1] Johnstone D, Moore G, Tausz M, et al. (2010) The measurement of wood decay in landscape trees. Arboric Urban For 36: 121-127.
    [2] Sambuelli L, Socco LV, Godio A, et al. (2003) Ultrasonic, electric and radar measurements for living trees assessment. B Geofis Teor Appl 44: 253-279.
    [3] Rabe C, Ferner D, Fink S, et al. (2004) Detection of decay in trees with stress waves and interpretation of acoustic tomograms. Arboricultural J 28: 3-19. doi: 10.1080/03071375.2004.9747399
    [4] Göcke L, Rust S, Weihs U, et al. (2007) Combining sonic and electrical impedance tomography for the non-destructive testing of trees. Proceedings of 15th international symposium on nondestructive testing of wood. Forest products society. Madison, 31-42.
    [5] Leong EC, Burcham DC, Yok-King F (2012) A purposeful classification of tree decay detection tools. Arboricultural J Int J Urban For 34: 91-115. doi: 10.1080/03071375.2012.701430
    [6] Mattheck CG, Breloer H (1994) Field guide for visual tree assessment (VTA). Arboricultural J Int J Urban For 18: 1-23. doi: 10.1080/03071375.1994.9746995
    [7] Hayes E (2001) Evaluating Tree Defects, Safe trees, 2 Eds., Rochester, 34.
    [8] Luley CL (2005) Wood decay fungi common to living urban trees in the Northeast and Central United States, Naples: Urban Forestry, 61.
    [9] Rinn F, Schweingruber FH, Schiir E (1996) RESISTOGRAPH and X-ray density charts of wood comparative evaluation of drill resistance profiles and X-ray density charts of different wood species. Holzforschung 50: 303-311. doi: 10.1515/hfsg.1996.50.4.303
    [10] Martin T, Günther T (2013) Complex Resistivity Tomography (CRT) for fungus detection on standing oak trees. Eur J For Res 132 : 765-776. doi: 10.1007/s10342-013-0711-4
    [11] Gilbert EA, Smiley ET (2004) Picus sonic tomography for the quantification of decay in white oak (Quercus alba) and hickory (Carya spp.). J Arboric 30: 277-280.
    [12] Wessoly L (1995) Fracture diagnosis of trees. Part 3: Boring is no way for reliable fracture diagnosis. Stadt Und Grün 9: 635-640.
    [13] Hagrey SA (2007) Geophysical imaging of root-zone, trunk, and moisture heterogeneity. J Exp Bot 58: 839-854. doi: 10.1093/jxb/erl237
    [14] Bassuk N, Grabosky J, Mucciardi A, et al. (2011) Ground-penetrating radar accurately locates tree roots in two soil media under pavement. Arboric Urban For 37: 160-166.
    [15] Zhu SP, Huang CL, Su Y, et al. (2014) 3D Ground penetrating radar to detect tree roots and estimate root biomass in the field. Remote Sens 6: 5754-5773. doi: 10.3390/rs6065754
    [16] Satriani A, Loperte A, Proto M, et al. (2010) Building damage caused by tree roots: laboratory experiments of GPR and ERT surveys. Adv Geosci 24: 133-137. doi: 10.5194/adgeo-24-133-2010
    [17] Lorenzo H, Pérez-Gracia V, Novo A, et al. (2010) Forestry applications of ground-penetrating radar. For Syst 19: 5-17.
    [18] Muller W (2003) Timber girder inspection using ground penetrating radar. OR Insight 45: 809-812.
    [19] Pyakurel S (2009) 2D and 3D GPR imaging of wood and fiber reinforced polymer composites. West Virginia University, Morgantown
    [20] Halabe UB, Agrawal S, Gopalakrishnan B (2009) Non-destructive evaluation of wooden logs using ground penetrating radar. Nondestr Test Eval 24: 329-346. doi: 10.1080/10589750802474344
    [21] Colla C (2010) GPR of a timber structural element. Proceedings of 13th international conference on ground penetrating radar, Lecce.
    [22] Rodriguez-Abad I, et al. (2010) Non-destructive methodologies for the evaluation of moisture content in sawn timber structures: ground-penetrating radar and ultrasound techniques. Near Surf Geophys 8: 475-482. doi: 10.3997/1873-0604.2010048
    [23] Jezova J, Lambot S (2019) Influence of bark surface roughness on tree trunk radar inspection. Ground Penetrating Radar 2: 1-25. Available from: https://doi.org/10.26376/GPR2019001.
    [24] Nicolotti G, Socco LV, Martinis R, et al. (2003) Application and comparison of three tomographic techniques for detection of decay in trees. J Arboric 29: 66-78.
    [25] Sudakova MS, Kalashnikov AU, Terentieva EB (2017) GPR tomography as applied to delineation of voids. Construction Unique Build Struct 60: 57.
    [26] Sihvola A (2000) Mixing rules with complex dielectric coefficients. Subsurf Sens Technol Appl 1: 393-415. doi: 10.1023/A:1026511515005
    [27] Huisman JA, Hubbard SS, Redman JD, et al. (2003) Measuring soil water content with ground penetrating radar: a review. Vadose Zone J 2: 476-491. doi: 10.2136/vzj2003.4760
    [28] Hartley I, Hamza MF (2016) Wood: Moisture Content, Hygroscopicity, and Sorption. In: Saleem Hashmi (Editor in chief), Reference Module in Materials Science and Materials Engineering, Oxford: Elsevier.
    [29] Torgovnikov G (1993) Dielectric properties of wood and wood-based materials. Berlin, Heidelberg, New-York: Springer-Verlag, 196.
    [30] Sahin H, Ay N (2004) Dielectric properties of hardwood species at microwave frequencies. J Wood Sci 50: 375-380. doi: 10.1007/s10086-003-0575-1
    [31] Sbartaï, ZM, Lataste JF (2008) Evaluation non-destructive de l'humidité du bois par mesures radar et résistivité électrique. Rapport interne, Université bordeaux 1, US2B, Ghymac.
    [32] Jol HM (2009) Ground penetrating radar theory and applications, Oxford: Elsevier, 523.
    [33] Kabir MF, Khalid KB, Daud WM, et al. (1997) Dielectric properties of rubber wood at microwave frequencies measured with an open-ended coaxial line. Wood Fiber Sci 29: 319-324.
    [34] Afzal MT, Colpits B, Galik K (2003) Dielectric Properties of Softwood Species Measured with an Open-ended Coaxial Probe. In proceedings 8th International IUFRO Wood Drying Conference, Brasvo, Romania, August, 24-29, 2003; 110-115.
    [35] Koubaa A, Perré P, Hutcheon RM, et al. (2008) Complex Dielectric Properties of the Sapwood of Aspen, White Birch, Yellow Birch, and Sugar Maple. Drying Technol 26: 568-578. doi: 10.1080/07373930801944762
    [36] Topp GC, Davis JL, Annan AP (1980) Electromagnetic determination of soil water content: measurements in coaxial transmission lines. Water Resour Res 16: 574-582. doi: 10.1029/WR016i003p00574
    [37] Lin G, Lv JX, Wen J (2014) Research on the relationship between moisture content and the dielectric constant of the tree trunk by the radar wave. Comput Modell New Technol 18: 1171-1175.
    [38] Razafindratsima S, Sbartaï ZM, Demontoux F (2017) Permittivity measurement of wood material over a wide range of moisture content. Wood Sci Technol 51: 1421-1431. doi: 10.1007/s00226-017-0935-4
    [39] Acuna L, Basterra LA, Casado M, et al. (2011) Application of resistograph to obtain the density and to differentiate wood species. Mater Construcc 61: 451-464. doi: 10.3989/mc.2010.57610
  • This article has been cited by:

    1. Xiu Ye, Shangyou Zhang, A weak divergence CDG method for the biharmonic equation on triangular and tetrahedral meshes, 2022, 178, 01689274, 155, 10.1016/j.apnum.2022.03.017
    2. Jun Zhou, Da Xu, Wenlin Qiu, Leijie Qiao, An accurate, robust, and efficient weak Galerkin finite element scheme with graded meshes for the time-fractional quasi-linear diffusion equation, 2022, 124, 08981221, 188, 10.1016/j.camwa.2022.08.022
    3. Xiu Ye, Shangyou Zhang, A conforming discontinuous Galerkin finite element method for the Stokes problem on polytopal meshes, 2021, 93, 0271-2091, 1913, 10.1002/fld.4959
    4. Xiu Ye, Shangyou Zhang, Constructing a CDG Finite Element with Order Two Superconvergence on Rectangular Meshes, 2023, 2096-6385, 10.1007/s42967-023-00330-5
    5. Yan Yang, Xiu Ye, Shangyou Zhang, A pressure-robust stabilizer-free WG finite element method for the Stokes equations on simplicial grids, 2024, 32, 2688-1594, 3413, 10.3934/era.2024158
    6. Xiu Ye, Shangyou Zhang, A superconvergent CDG finite element for the Poisson equation on polytopal meshes, 2023, 0044-2267, 10.1002/zamm.202300521
    7. Xiu Ye, Shangyou Zhang, Two-Order Superconvergent CDG Finite Element Method for the Heat Equation on Triangular and Tetrahedral Meshes, 2024, 2096-6385, 10.1007/s42967-024-00444-4
    8. Xiu Ye, Shangyou Zhang, Order two superconvergence of the CDG finite elements for non-self adjoint and indefinite elliptic equations, 2024, 50, 1019-7168, 10.1007/s10444-023-10100-9
    9. Fuchang Huo, Weilong Mo, Yulin Zhang, A locking-free conforming discontinuous Galerkin finite element method for linear elasticity problems, 2025, 465, 03770427, 116582, 10.1016/j.cam.2025.116582
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4272) PDF downloads(259) Cited by(5)

Figures and Tables

Figures(7)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog