Research article

Prevalence of goiter among children aged 8–10 in Binh Dinh province, Vietnam in 2016–2017

  • Received: 13 February 2019 Accepted: 17 May 2019 Published: 23 May 2019
  • Objective: The study was conducted to estimate the goiter prevalence, and the median urine iodine concentrations among schoolchildren aged 8–10 in Binh Dinh province, Vietnam. Methods: A school-based cross-sectional survey was carried out from May 2016 to May 2017. A multistage, proportional-to-population-size sampling method with 30 clusters was used. The children were examined by palpation for the presence or absence of goiter based on the criteria of the World Health Organization (WHO), the United Nations Children’s Fund (UNICEF), International Council for the Control of Iodine Deficiency (ICCIDD); urinary iodine was determined in microplates by a modification of the Sandell-Kolthoff reaction. The Chi-square test was used to compare prevalences, and the Chi-square test for trend was employed to assess the trend of goiter prevalence and urine iodine levels by age and economic-social areas. Results: 1800 pupils from 8 to 10 years old including 900 males and 900 females were examined and 300 among them were tested for the urinary iodine concentration (UIC). The prevalence of goiter among schoolchildren was 6.6%. The prevalence of goiter tended to increase in areas with disadvantaged conditions, among which the urban areas occupied the lowest prevalence (5%) while the mountainous areas and Midland took the highest (8.8%) (the p-value of 0.0193). The median UIC of the study group was 159.9 µg/L; the 25th and 75th percentile value was 103 µg/L and 230.2 µg/L, respectively. Conclusion: According to the WHO/UNICEF/ICCIDD classification, the goiter prevalence indicated that some regions of Binh Dinh province appeared to be slightly affected by iodine deficiency. These have characterized an important public health challenge, highlighting the need to eliminate iodine deficiency disorders in these areas.

    Citation: Truong Quang Dat, Le Nguyen Huong Giang, Pham Van Bao, Nguyen Thi Hong Tuyen. Prevalence of goiter among children aged 8–10 in Binh Dinh province, Vietnam in 2016–2017[J]. AIMS Public Health, 2019, 6(2): 184-194. doi: 10.3934/publichealth.2019.2.184

    Related Papers:

    [1] Lianbing She, Nan Liu, Xin Li, Renhai Wang . Three types of weak pullback attractors for lattice pseudo-parabolic equations driven by locally Lipschitz noise. Electronic Research Archive, 2021, 29(5): 3097-3119. doi: 10.3934/era.2021028
    [2] Xiu Ye, Shangyou Zhang, Peng Zhu . A weak Galerkin finite element method for nonlinear conservation laws. Electronic Research Archive, 2021, 29(1): 1897-1923. doi: 10.3934/era.2020097
    [3] Guifen Liu, Wenqiang Zhao . Regularity of Wong-Zakai approximation for non-autonomous stochastic quasi-linear parabolic equation on $ {\mathbb{R}}^N $. Electronic Research Archive, 2021, 29(6): 3655-3686. doi: 10.3934/era.2021056
    [4] Jiwei Jia, Young-Ju Lee, Yue Feng, Zichan Wang, Zhongshu Zhao . Hybridized weak Galerkin finite element methods for Brinkman equations. Electronic Research Archive, 2021, 29(3): 2489-2516. doi: 10.3934/era.2020126
    [5] Jun Pan, Yuelong Tang . Two-grid $ H^1 $-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
    [6] 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
    [7] 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
    [8] Noelia Bazarra, José R. Fernández, Ramón Quintanilla . Numerical analysis of a problem in micropolar thermoviscoelasticity. Electronic Research Archive, 2022, 30(2): 683-700. doi: 10.3934/era.2022036
    [9] Yaojun Ye, Qianqian Zhu . Existence and nonexistence of global solutions for logarithmic hyperbolic equation. Electronic Research Archive, 2022, 30(3): 1035-1051. doi: 10.3934/era.2022054
    [10] 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
  • Objective: The study was conducted to estimate the goiter prevalence, and the median urine iodine concentrations among schoolchildren aged 8–10 in Binh Dinh province, Vietnam. Methods: A school-based cross-sectional survey was carried out from May 2016 to May 2017. A multistage, proportional-to-population-size sampling method with 30 clusters was used. The children were examined by palpation for the presence or absence of goiter based on the criteria of the World Health Organization (WHO), the United Nations Children’s Fund (UNICEF), International Council for the Control of Iodine Deficiency (ICCIDD); urinary iodine was determined in microplates by a modification of the Sandell-Kolthoff reaction. The Chi-square test was used to compare prevalences, and the Chi-square test for trend was employed to assess the trend of goiter prevalence and urine iodine levels by age and economic-social areas. Results: 1800 pupils from 8 to 10 years old including 900 males and 900 females were examined and 300 among them were tested for the urinary iodine concentration (UIC). The prevalence of goiter among schoolchildren was 6.6%. The prevalence of goiter tended to increase in areas with disadvantaged conditions, among which the urban areas occupied the lowest prevalence (5%) while the mountainous areas and Midland took the highest (8.8%) (the p-value of 0.0193). The median UIC of the study group was 159.9 µg/L; the 25th and 75th percentile value was 103 µg/L and 230.2 µg/L, respectively. Conclusion: According to the WHO/UNICEF/ICCIDD classification, the goiter prevalence indicated that some regions of Binh Dinh province appeared to be slightly affected by iodine deficiency. These have characterized an important public health challenge, highlighting the need to eliminate iodine deficiency disorders in these areas.


    In the recent few centuries, partial differential equations have provided powerful scientific tools for describing many mathematical models. However, with the development of science and the further understanding of natural phenomena, the classical differential equations can not explain random phenomena in the nature and other fields well. In this context, the stochastic partial differential equations (SPDEs for short) are brought in as the models for depicting a variety of random phenomena [1]. The general form of such equations is usually formulated as

    $ du=(Au+F(u))dt+G(u)dW(t),u(0)=u0.
    $

    Here, $ H $ is a Hilbert space, and $ u(t) $ is an $ H $-valued random process. We denote by $ u_0\in H $ the initial value. $ A $ is a linear, self-adjoint, positive definite, not necessarily bounded operator with a compact inverse and densely defined in a subspace of $ H $. $ F $ and $ G $ are usually nonlinear operators on $ H $. $ W(t) $ is an $ H $-valued $ Q $-Wiener process defined in a filtered probability space $ (\Omega, \mathcal{F}, P, \mathcal{F}_{t}) $.

    In this paper, we mainly concern the following SPDE with additive noise,

    $ du+(Au+F(u))dt=GdW(t),in D,0tT,u=0,on D,0tT,u(0)=u0,in D,
    $
    (1.1)

    where $ \mathcal{D}\subset\mathbb R^{2} $ is a polygonal domain and the linear operator $ G $ is independent of $ u $.

    Up to now, different kinds of numerical methods have been applied to solving the form of the SPDE (1.1), such as finite difference methods, finite element methods, discontinuous Galerkin methods, WG methods, etc. In [2], the author presents a finite difference method for stochastic nonlinear second-order boundary-value problems (BVPs) driven by additive noises, and proves that the finite difference solution converges to the solution to the original stochastic BVP at $ O(h) $ in the mean-square sense. The stochastic Allen-Cahn equation with additive noise is discretized by means of a spectral Galerkin method in space and a tamed version of the exponential Euler method in time [3]. A discontinuous Galerkin method is applied in [4] for stochastic differential equations driven by additive noise, and convergence analysis is provided. In [5], the authors analyze strong approximation errors of both finite element semi-discretization and spatio-temporal full discretization for the stochastic Allen-Cahn equation driven by additive noise in space dimension $ d\leq 3 $. A stochastic analogue of the local discontinuous Galerkin method is constructed for a stochastic two-point boundary-value problem driven by an additive white noise [6]. In [7], the authors adopt the Argyris finite elements to solve (1.1) with $ A = -\Delta^{2} $ and obtain the optimal order $ h^{\beta} $ of error estimates with $ \beta > 0 $. In [8,9], the linear version of (1.1) is investigated via the WG methods and the optimal order estimates in the sense of strong convergence are derived.

    In general, WG methods, firstly proposed by Wang and Ye [10], are newly developed numerical techniques for solving partial differential equations. The essence of this method is the use of weak finite element functions and their weak derivatives computed with a framework that mimics the distribution or generalized functions. Since the method was put forward, it has been applied/extended to different kinds of partial differential equations, such as biharmonic equations [11], Stokes equations [13], linear elasticity equations [22], poroelasticity problems [14], parabolic problems [15,16], SPDEs [8,9], etc.

    In this paper, we adopt the WG method with a parameter-free stabilization term for solving the SPDE (1.1) with $ A = -\Delta $. The main characteristic of this method is that the WG finite element space consists of discontinuous functions, which allows the WG method applied to the general polygonal or polyhedral meshes. This characteristic makes the WG method efficient and highly flexible. The optimal order for strongly convergent error estimates in $ L_2 $-norm is studied based on the established semi-discrete WG scheme. As far as we know, this paper is the first to apply the ideas of WG to nonlinear stochastic models for error analysis.

    This paper is organized as follows. In Section 2, we provide several definitions and assumptions as the preliminaries for the theoretical analysis. Section 3 introduces the details of WG method and sets up the semi-discrete WG scheme of the nonlinear stochastic model. Several error estimates for the related deterministic problem are supplied in Section 4, which is helpful for the derivation of our later main result. And finally, in Section 5, we derive the optimal order for strongly convergent error estimates in $ L_2 $-norm based on the established semi-discrete WG scheme.

    In this section, we introduce several definitions and assumptions as a preparation for the later theoretical analysis.

    Recall that $ \mathcal{D}\subset\mathbb R^{2} $ is a polygonal domain. Unless particularly stated, in this paper, we shall use the standard notations for Sobolev spaces and their associated norms [17]. Let $ H = L_2(\mathcal{D}) $ whose inner product and norm are denoted by $ (\cdot, \cdot) $ and $ \parallel\cdot\parallel $, respectively. Denote by $ H^s = H^s(\mathcal{D}) $ with norm $ \|\cdot\|_s $. Let $ H_0 $ and $ H_0^s $ be the subspaces of $ H $ and $ H^s $, respectively, the elements of which vanish on the boundary $ \partial\mathcal{D} $.

    Denote by $ Q: H\rightarrow H $ a linear self-adjoint operator with eigenvalues $ \gamma_{i} > 0 $ ($ i = 1, 2, \cdots $) and corresponding normalized eigenfunctions $ e_i\in H $ ($ i = 1, 2, \cdots $). Then $ e_i $ ($ i = 1, 2, \cdots $) form a family of completely orthonormal bases of the space $ H $. We further assume:

    ($ \bf{H1} $): The operator $ Q $ is bounded and positive definite.

    ($ \bf{H2} $): The operator $ Q $ has bounded trace, i.e., $ Tr(Q) = \sum\limits_{i = 1}^{\infty}\gamma_{i} < +\infty $.

    Since $ W(t) $ is a $ Q $-Wiener process defined on a given filtrated probability space $ (\Omega, \mathcal{F}, P, \mathcal{F}_{t}) $, $ W(t) $ is $ H $-valued. Denote by $ \beta_{i}(t) $ ($ i = 1, 2, ... $) a family of Brownian motions in which the elements are independently and identically distributed. Then the Wiener process $ W(t) $ can be written in the form of its Fourier expansion [18]:

    $ W(t)=i=1γ1/2ieiβi(t).
    $

    Next, we define several operator spaces. Let $ L(Q^{1/2}(H), H) $ be the space of bounded linear operators from $ Q^{1/2}(H) $ to $ H $ and denote by $ L_{2}^{0}(Q^{1/2}(H), H) $ a subspace of $ L(Q^{1/2}(H), H) $ satisfying:

    $ L02(Q1/2(H),H)={ϕL(Q1/2(H),H):l=1ϕQ1/2el2<}.
    $

    Let $ L(H) $ be the linear bounded operator space from $ H $ to $ H $. Then we define a Hilbert-Schmidt operator space $ L_{HS}(H)\subset L(H) $, i.e.,

    $ LHS(H)={ΦL(H):i=1Φei2<}
    $

    with norm

    $ ΦHS=(i=1Φei2)1/2,ΦLHS(H).
    $

    It is not hard to see that for any $ \psi\in L_{2}^{0}(Q^{1/2}(H), H) $, the operator $ \psi Q^{1/2}\in L_{HS}(H) $. Let $ {\bf{E}} $ represent the standard mathematical expectation. Then, for any $ \psi\in L_{2}^{0}(Q^{1/2}(H), H) $, the following isometry equation holds.

    $ Et0ψ(s)dW(s)2=t0ψ(s)Q1/22HSds.
    $
    (2.1)

    Additionally, define the space $ L_{2}(\Omega; H) $,

    $ L2(Ω;H)={v:Ωv2dP(ω)<}
    $
    (2.2)

    with norm

    $ vL2(Ω;H)=(Ωv2dP(ω))1/2,vL2(Ω;H).
    $

    Similarly, we can define $ L_{\infty}(\Omega; H) $.

    Let $ A = -\Delta $. Assume the operator $ A $ defined on $ H_{0}^{2} $ has real eigenvalues $ \lambda_{i} > 0 $ ($ i = 1, 2, ... $) with corresponding eigenfunctions $ g_{i}\in H_{0}^{2}\subset H $ ($ i = 1, 2, ... $). Then for all $ v\in H_{0}^{2} $, we have

    $ Av=i=1λi(v,gi)gi.
    $

    For $ s > 0 $, we define the space $ \dot{H}^{s} $:

    $ ˙Hs={vH:∥As/2v∥=(i=1λsi(v,gi)2)1/2<}
    $

    with norm $ \mid\cdot\mid_{s} = \parallel A^{s/2}\cdot\parallel $. Similar to (2.2), we can define $ L_{2}(\Omega; \dot{H}^{s}) $.

    The lemma below provides the relationship between $ \dot{H}^{s} $ and $ H_0^s $.

    Lemma 2.1. [19,Lemma 3.1] For any $ s > 0 $, we have

    $ ˙Hs={vHs:Δjv=0onD,j<s/2}.
    $

    Moreover, $ \mid\cdot\mid_{s} $ is equivalent to $ \parallel\cdot\parallel_{s} $ in $ \dot{H}^{s} $, where

    $ vs={Δpv,ifs=2p,(Δpv),ifs=2p+1.
    $
    (2.3)

    Now we write three assumptions for the SPDE (1.1) we consider.

    ($ \bf{H3} $): The initial value $ u_{0}\in L_{2}(\Omega; \dot{H}^{2}) $.

    ($ \bf{H4} $): The operator $ F:H_{0}\rightarrow \dot H^{1} $ satisfies:

    $ F(u1)F(u2)1Cu1u2,u1,u2H0,
    $
    (2.4)

    where $ C $ is a positive constant. In this paper, the letter $ C $ denotes a generic positive constant which may be different at different occurrences. Furthermore, if we take $ u_{2}\equiv 0 $, then

    $ F(u1)1C(u1u2+F(u2)1)C(u1+1),u1H0.
    $
    (2.5)

    ($ \bf{H5} $): The operator $ G\in L_{2}^{0}(Q^{1/2}(H), H) $ satisfies

    $ A1/2GQ1/2HS<.
    $
    (2.6)

    Let $ E(t) = e^{-tA} $. Then (1.1) admits a unique mild solution of the form [18]:

    $ u(t)=E(t)u0t0E(ts)F(u(s))ds+t0E(ts)GdW(s).
    $
    (2.7)

    From [20], the following two lemmas hold true.

    Lemma 2.2. If the assumptions $ \bf{H3}-\bf{H5} $ hold, $ u $ is the mild solution of (1.1). Then for $ 0\le t_{1}\le t_{2}\le T $, there exists a constant $ C $ such that

    $ u(t1)u(t2)L2(Ω;H)Ct1t21/2.
    $
    (2.8)

    Lemma 2.3. If the assumptions $ \bf{H3}-\bf{H5} $ hold, $ u $ is the mild solution of (1.1). Then

    $ sups[0,T]E[u(s)2]<.
    $
    (2.9)

    In this section, the details of WG method is introduced, and the semi-discrete WG scheme for the SPDE (1.1) is established.

    Let $ \mathcal{T}_{h} $ be a regular partition of the domain $ \mathcal{D} $ satisfying the shape regularity requirements A1-A4 in [21]. The boundary of each element $ K\in\mathcal{T}_{h} $ is denoted by $ \partial K $. Let $ h_{K} $ be the diameter of $ K $ and $ h = \max\limits_{K\in\mathcal{T}_{h}} h_{K} $.

    We define the weak function space:

    $ W(K){v={v0,vb};v0L2(K),vbL2(K)}.
    $

    Here $ v_{0} $ is the value of $ v $ in $ K $ and $ v_{b} $ is the value of $ v $ on the boundary $ \partial K $.

    With an inclusion map $ i : H^{1}(K)\rightarrow W(K) $ [21],

    $ i(v)={vK,vK}, vH1(K),
    $

    we can embed $ H^{1}(K) $ into the weak function space $ W(K) $.

    For any $ u, v\in W(K) $, denote by $ (u, v)_{K} $ the standard $ L^{2} $-inner product in $ K $ and by $ \parallel\cdot\parallel_{K} $ the corresponding norm. Similarly, $ \langle u, v\rangle_{\partial K} $ represents the $ L^{2} $-inner product on $ \partial K $, and the norm is notated by $ \parallel\cdot\parallel_{\partial K} $.

    Let $ r $ be a non-negative integer. For each $ K\in\mathcal{T}_{h} $, let $ P_{r+1}(K) $ and $ P_{r+1}(\partial K) $ be the sets of polynomials with degree no more than $ r+1 $ in $ K $ and on $ \partial K $, respectively. Define a discrete weak function space $ W_{r+1}(K)\subset W(K) $, i.e.,

    $ Wr+1(K){v={v0,vb};v0Pr+1(K),vbPr+1(K)}.
    $

    For each $ K\in\mathcal{T}_{h} $, let $ V(K, r) = [P_{r}(K)]^{2} $. Define the discrete weak gradient operator $ \nabla_{d}:W_{r+1}(K)\rightarrow V(K, r) $, such that for any $ v = \{v_{0}, v_{b}\}\in W_{r+1}(K) $, $ \nabla_{d}v\in V(K, r) $ is the unique vector-valued polynomial satisfying:

    $ (dv,q)K(v0,(q))K+vb,(qn)K,qV(K,r).
    $
    (3.1)

    Now we extend the definition of the weak function space $ W(K) $ and the discrete weak function space $ W_{r+1}(K) $ to the whole domain $ \mathcal{D} $. Denote by $ W(\mathcal{D}) = \{v:v\mid_{K}\in W(K)\} $ the weak function space defined on the domain $ \mathcal{D} $. Let $ S_{h}(r+1)\subset W(\mathcal{D}) $ be the the discrete weak function space satisfying:

    $ Sh(r+1){v={v0,vb};vKWr+1(K),KTh},
    $

    and denote by $ S_{h}^{0}(r+1) $ the subspace of $ S_{h}(r+1) $ with vanishing values on the boundary $ \partial\mathcal{D} $. If no confusion occurs, we use respectively $ S_{h} $ to denote $ S_{h}(r+1) $, and use $ S_{h}^{0} $ to denote $ S_{h}^{0}(r+1) $ throughout this paper.

    We also define the following global vector-valued polynomial space:

    $ V(r){q:qKV(K,r)}.
    $

    Then we extend the definition of the weak gradient operator $ \nabla_{d} $ from each $ K\in\mathcal{T}_{h} $ to the whole domain $ \mathcal{D} $. That is to say, for any $ v\in S_{h}(r+1) $, define the operator $ \nabla_{d}: S_{h}(r+1)\rightarrow V(r) $, such that

    $ dvKd(vK).
    $
    (3.2)

    Next, we bring in several locally defined projection operators which are helpful for our theoretical analysis. For each $ K\in\mathcal{T}_{h} $, the $ L_2 $-projection operators $ Q_0:L_{2}(K)\rightarrow P_{r+1}(K) $ and $ Q_b:L_{2}(\partial K)\rightarrow P_{r+1}(\partial K) $ are defined piecewisely in $ K $ and on $ \partial K $, respectively. Define $ Q_{h}u = \{Q_{0}u_{0}, Q_{b}u_{b}\}:W(\mathcal{D})\rightarrow S_{h}(r+1) $ satisfying:

    $ QhK:W(K)L2Wr+1(K),KTh.
    $

    For each $ K\in\mathcal{T}_{h} $, let $ R_{h}:[L_{2}(\mathcal{D})]^{2}\rightarrow V(r) $ be the $ L_2 $-projection operator defined by

    $ RhK:[L2(K)]2L2V(K,r),KTh.
    $
    (3.3)

    Now we introduce a bilinear form as follows: For any $ u_{h}, v_{h}\in S_{h} $,

    $ as(uh,vh)=KTh(duh,dvh)K+KThh1Ku0ub,v0vbK,uh,vhSh.
    $

    Take the following elliptic problem into account:

    $ Δu=f,inD,u=0,onD.
    $
    (3.4)

    The variational formulation for (3.4) is to find $ u\in H_{0}^{1}(\mathcal{D}) $ such that

    $ (u,v)=(f,v),vH10(D).
    $
    (3.5)

    The WG scheme for (3.4) is to find $ u_{h}\in S_{h}^{0} $ such that

    $ as(uh,vh)=(f,vh)=(Qhf,vh),vhS0h.
    $
    (3.6)

    Let $ A_{h}:S_{h}^{0}\rightarrow S_{h}^{0} $ be the operator satisfying:

    $ (Ahuh,vh)=as(uh,vh),vhS0h.
    $

    Here $ A_{h} $ is a linear, self-adjoint, symmetric, positive definite operator. Then, the numerical scheme (3.6) is equivalent to finding $ u_{h}\in S_{h}^{0} $ such that

    $ Ahuh=Qhf.
    $
    (3.7)

    Based on [21,Theorem 8.2] and [8,(4.7)], we have the following lemma about the error estimates of the elliptic problem (3.4) with WG method.

    Lemma 3.1. Assuming $ u\in H^{2}(\mathcal{D}) $ and $ u_h\in S_{h} $ are the solutions of (3.4) and (3.7), respectively. Then there exists a positive constant $ C $ which depends only on the domain $ \mathcal{D} $, such that

    $ Qhuuh∥≤Ch2f.
    $
    (3.8)

    Define two operators $ G = A^{-1} $ and $ G_{h} = A_{h}^{-1} $. $ G: H_{0}\rightarrow H_{0}^{1} $ and $ G_{h}: S_{h}^{0}\rightarrow S_{h}^{0} $ are the solution operators of (3.4) and (3.7), respectively. In other words, $ u = Gf $ and $ u_{h} = G_{h}Q_{h}f $ are respectively the solutions of (3.4) and (3.7). It is easy to see that $ G $ and $ G_h $ are linear, symmetric, positive definite operators. Then (3.8) can be written as

    $ (QhGGhQh)f∥≤Ch2f.
    $
    (3.9)

    Now we approximate the stochastic problem (1.1) with WG method. The semi-discrete WG scheme for (1.1) is to find an $ H- $valued random process $ u_{h}(\cdot, t)\in S_{h}^{0} $ with $ u_{h}(0) = Q_{h}u_{0} $, such that for any $ v_{h}\in S_{h}^{0} $ and $ 0\le t\le T $,

    $ (uh(t),vh)(uh(0),vh)+t0(Ahuh,vh)ds+t0(QhF(uh),vh)ds=(t0QhGdW(s),vh).
    $
    (3.10)

    In fact, (3.10) is equivalent to

    $ duh+Ahuhdt+QhF(uh)dt=QhGdW,in D, 0tT,uh(0)=Qhu0,in D.
    $
    (3.11)

    Let $ E_{h}(t) = e^{-tA_{h}} $, $ t\geq 0 $. It is obvious that the equation (3.11) has a mild solution

    $ uh(t)=Eh(t)Qhu0t0Eh(ts)QhF(uh(s))ds+t0Eh(ts)QhGdW(s).
    $
    (3.12)

    In this section, we provide several error estimates with respect to the related deterministic problem, which are used in the subsequent section.

    Lemma 4.1. [12,Lemma 3.2] For any $ \alpha $, $ \beta\in\mathbb{R} $ and $ l\ge 0 $, we have

    $ DltE(t)vβCt(βα)/2lvα, t>0, 2l+βα,
    $
    (4.1)

    and

    $ t0sαDltE(s)v2βdsCv22l+βα1, t>0, α0,
    $
    (4.2)

    where $ D_{t}^{l} $ is the $ l $-th derivative with respect to $ t $.

    Consider the following two problems:

    $ ut+Au=0,u(0)=u0,
    $
    (4.3)

    and

    $ uh,t+Ahuh=0,uh(0)=Qhu0.
    $
    (4.4)

    Then $ u = E(t)u_{0} $, $ u_{h} = E_{h}(t)Q_{h}u_{0} $ are the solutions of (4.3) and (4.4), respectively. Making use of the forms of $ u $ and $ u_{h} $, the error equation is shown as follows:

    $ Ghet+e=ρ,
    $
    (4.5)

    where $ e(t) = u_{h}(t)-Q_{h}u(t) $, $ \rho = (Q_{h}G-G_{h}Q_{h})u_{t} $ and $ e(0) = u_{h}(0)-Q_{h}u(0) = 0 $.

    Indeed, by virtue of (4.3) and (4.4), we supply

    $ Ghet+e=(Ghuh,t+uh)(GhQhut+Qhu)=Gh(uh,t+Ahuh)(GhQhut+Qhu)=(QhGGhQh)ut.
    $

    According to the error equation (4.5), it follows from [8] that

    $ t0e2dst0ρ2ds,
    $
    (4.6)

    and

    $ te214te2+4tρ2+t0(2sρt∥∥e+e2+2ρ∥∥e)ds.
    $
    (4.7)

    Furthermore, we can derive the following two lemmas.

    Lemma 4.2. For any $ t\in[0, T] $, assume that $ e(t)\in S_{h}^{0} $, $ e(0) = 0 $ and (4.5) holds. Then there exists a positive constant $ C $, such that

    $ e(t)∥≤C(sup0stsρt(s)+sup0stρ(s)),t0.
    $
    (4.8)

    Proof. Due to (4.7), (4.6) and the mean inequality, we obtain

    $ e(t)2C(sup0sts2ρt(s)2+sup0stρ(s)2),t0,
    $

    which completes the proof.

    Lemma 4.3. Under the same assumptions of Lemma 4.2, for any fixed positive number $ \epsilon $, there exists a positive constant $ C_{\epsilon} $ depending on $ \epsilon $, such that

    $ e(t)∥≤(ϵsup0stsρt(s)+Cϵsup0stρ(s)),t0.
    $
    (4.9)

    Proof. The proof is similar to the one in Lemma 4.2. Noticing that

    $ 2sρt∥∥e∥≤(ϵ2s2ρt2+1ϵ2e2),
    $
    (4.10)

    together with (4.7) and (4.6), we finish the proof.

    Notate $ F_{h}(t) = E_{h}(t)Q_{h}-Q_{h}E(t) $. Then we render the following estimate results.

    Lemma 4.4. [8,Lemma 4.4] For $ v\in\dot{H}^{2} $, then we have

    $ FhvL([0,T];H)=sup0tTFhv∥≤Ch2v2.
    $
    (4.11)

    If $ v\in\dot{H}^{1} $ and $ 0\le t\le T $, then

    $ FhvL2([0,t];H)=(t0Fhv2ds)1/2Ch2v1,
    $
    (4.12)

    where $ C $ is a positive constant only depending on the domain $ \mathcal{D} $.

    Theorem 4.1. Assuming that $ u_{0}\in \dot{H}^{1} $ and $ t > 0 $, then we have

    $ t0Fhu0ds∥≤Ch2u0,
    $
    (4.13)

    where $ C $ is a positive constant only depending on the domain $ \mathcal{D} $.

    Proof. Denote by $ \tilde{e}(t) $ and $ \tilde{\rho}(t) $ the integral $ \int_{0}^{t}e(s)ds $ and $ \int_{0}^{t}\rho(s)ds $, respectively. Together with $ e(0) = 0 $, we acquire

    $ Gh˜et+˜e=Ghe(t)+t0e(s)ds=t0(Ghet+e)ds=˜ρ.
    $

    It is easy to check $ \tilde{e}(0) = 0 $. According to (4.8), we provide

    $ ˜e(t)∥≤C(sup0stsρ(s)+sup0st˜ρ(s)).
    $
    (4.14)

    It follows from Lemma 4.1 and (3.9) that

    $ sρ(s)∥=s(GhQhQhG)ut(s)∥≤Ch2sut(s)∥≤Ch2u0,
    $
    (4.15)

    and

    $ ˜ρ(s)∥=∥s0(GhQhQhG)ut(τ)dτ∥=∥(GhQhQhG)(u(s)u0)∥≤Ch2u0.
    $
    (4.16)

    Here we notice that $ E_{h}(s) = e^{-sA_{h}} $ is a bounded operator. Then the proof is completed.

    Theorem 4.2. Under the same assumptions of Theorem 4.1, we supply

    $ Fhu0∥≤Ch2t1u0,
    $
    (4.17)

    where $ C $ is a positive constant only depending on the domain $ \mathcal{D} $.

    Proof. Define $ \hat{e}(t) = te $, then by (4.5),

    $ Ghˆet+ˆe=Ghe+Gh(tet)+te=Ghe+tρ.
    $

    We denote $ \chi(t) = G_{h}e(t)+t\rho(t) $. It is easy to check $ \hat{e}(0) = 0 $. Due to (4.9), for any $ \epsilon > 0 $, we present

    $ te(t)∥=∥ˆe(t)∥≤ϵsup0stsχt(s)+Cϵsup0stχ(s).
    $

    For the two parts on the right side of the estimate above, we have

    $ χ(s)∥≤sρ(s)+Ghe(s),
    $

    and it follows from (4.5) that

    $ sχt(s)s2ρt(s)+sρ(s)+sGhet(s)s2ρt(s)+sρ(s)+s(ρ(s)+e(s))=s2ρt(s)+2sρ(s)+ˆe(s).
    $

    Let $ \epsilon = 1/2 $, then

    $ ˆe(t)∥≤(1/2)sup0stˆe(s)+C(sup0sts2ρt(s)+sup0stsρ(s)+sup0stGhe(s)).
    $

    Choose $ 0\le s_{0}\le t $ such that $ \parallel\hat{e}(s_{0})\parallel = \sup\limits_{0\le s\le t}\parallel\hat{e}(s)\parallel $, then

    $ ˆe(t)∥≤∥ˆe(s0)∥≤(1/2)ˆe(s0)+C(sup0sts2ρt(s)+sup0stsρ(s)+sup0stGhe(s)).
    $

    Hence

    $ te(t)∥=∥ˆe(t)∥≤∥ˆe(s0)∥≤C(sup0sts2ρt(s)+sup0stsρ(s)+sup0stGhe(s)).
    $
    (4.18)

    Now we estimate $ \parallel G_{h}e(s)\parallel $. Let $ \tilde{e}(t) = \int_{0}^{t}e(s)ds $, then

    $ Ghe+˜e=Gh˜et+˜e=˜ρ.
    $

    By virtue of (4.14), we obtain

    $ Ghe(t)∥≤∥˜e(t)+˜ρ(t)∥≤C(sup0stsρ(s)+sup0st˜ρ(s)).
    $
    (4.19)

    Combining (4.18) with (4.19), we have

    $ e(t)∥≤Ct1(sup0sts2ρt(s)+sup0stsρ(s)+sup0st˜ρ(s)).
    $

    Because of Lemma 4.1 and (3.9), we render

    $ s2ρt(s)∥=s2(GhQhQhG)utt(s)∥≤Ch2s2utt(s)∥≤Ch2u0.
    $

    With the help of (4.15) and (4.16), we complete the proof.

    In this section, we derive the optimal order for strongly convergent error estimates between the mild solution (2.7) of the SPDE (1.1) and its semi-discrete WG approximation (3.12) in $ L_{2} $-norm.

    The next lemma plays a very important role in getting our main result.

    Lemma 5.1. [20] For all $ C_{1}, C_{2}\ge 0, \alpha > 0 $, $ t\in[0, T] $, let $ \phi:[0, T]\rightarrow \mathbb{R} $ be a nonnegative and continuous function. If

    $ ϕ(t)=C1+C2t0(ts)α1ϕ(s)ds,
    $
    (5.1)

    then there exists a constant $ C = C(C_{2}, T, \alpha) $ such that for all $ t\in[0, T] $,

    $ ϕ(t)CC1.
    $
    (5.2)

    Theorem 5.1. Let $ u $ and $ u_{h} $ be the mild solutions of (1.1) and (3.11), respectively. If the assumptions $ \bf{H3}-\bf{H5} $ in Section 2 hold, then there exists a constant $ C $ depending on the domain $ \mathcal{D} $ and the upper bound of time $ T $ such that

    $ uh(t)Qhu(t)L2(Ω;H)Ch2.
    $

    Proof. Making use of (2.7) and (3.12), we write

    $ uh(t)Qhu(t)=Fh(t)u0+t0(QhE(ts)F(u(s))Eh(ts)QhF(uh(s))ds)ds+t0Fh(ts)GdW(s):=F1+F2+F3.
    $

    We record the three parts on the right side as $ F_{1} $, $ F_{2} $ and $ F_{3} $, respectively, and estimate each part one by one. Firstly, we can easily obtain $ \parallel F_{1}\parallel_{L_{2}(\Omega; H)}\le h^{2}\parallel u_{0}\parallel_{L_{2}(\Omega; \dot{H}^{2})} $ by (4.11).

    $ F_{2} $ can be written as a combination of three parts:

    $ F2=t0QhE(ts)F(u(s))dst0Eh(ts)QhF(uh(s))ds=t0Eh(ts)Qh(F(u(s))F(uh(s)))dst0Fh(F(u(s))F(u(t)))dst0Fh(F(u(t)))ds=I1I2I3.
    $

    Next, we estimate $ I_{1} $, $ I_{2} $ and $ I_{3} $, respectively. $ E_{h}(t-s) $ and $ Q_{h} $ are both bounded operators. Then together with (2.4), we find

    $ I1L2(Ω;H)Ct0uh(s)Qhu(s)L2(Ω;H)ds.
    $

    By Theorem 4.2, Lemma 2.1, (2.4) and (2.8), we have

    $ I2L2(Ω;H)t0Fh(F(u(s))F(u(t)))L2(Ω;H)dsCh2t0(ts)1F(u(s))F(u(t))L2(Ω;H)dsCh2t0(ts)1F(u(s))F(u(t))L2(Ω;˙H1)dsCh2t0(ts)1u(s)u(t)L2(Ω;H)dsCh2t0(ts)1/2dsCt1/2h2CT1/2h2.
    $

    From Theorem 4.1, Lemma 2.1, (2.5) and (2.9), we supply

    $ I3L2(Ω;H)Ch2F(u(t))L2(Ω;H)Ch2F(u(t))L2(Ω;˙H1)Ch2(u(t)L2(Ω;H)+1)Ch2(sups[0,T]u(s)L2(Ω;H)+1).
    $

    Thus,

    $ F2L2(Ω;H)Ch2+Ct0uh(s)Qhu(s)L2(Ω;H)ds.
    $

    Now we consider $ F_{3} $. With the help of (2.1), we provide

    $ F32L2(Ω;H)=Et0Fh(ts)GdW(s)2=t0Fh(ts)GQ1/22HSds=l=1t0Fh(ts)GQ1/2el2ds.
    $

    From (4.12), it follows that

    $ F32L2(Ω;H)Cl=1h4gQ1/2el21=Ch4A1/2gQ1/22HS.
    $

    Hence,

    $ uh(t)Qhu(t)L2(Ω;H)F1L2(Ω;H)+F2L2(Ω;H)+F3L2(Ω;H)=Ch2+Ct0uh(s)Qhu(s)L2(Ω;H)ds.
    $

    By Lemma 5.1 with $ \alpha = 1 $ and $ \phi(t) = \parallel u_{h}(t)-Q_{h}u(t)\parallel_{L_{2}(\Omega; H)} $, the proof is completed.

    This work was sponsored by the Research Foundation for Beijing University of Technology New Faculty Grant No. 006000514122516.

    All authors declare no conflicts of interest in this paper.


    Acknowledgments



    The research team sincerely thanks Binh Dinh Province People's Committee, the teachers and pupils at the research sites for their assistance in the implementation of this work.
    Financial support and sponsorship: Binh Dinh Province People's Committee.

    Conflict of interest



    All authors declare no conflicts of interest in this paper.

    [1] Gunnarsdottir I, Dahl L (2012) Iodine intake in human nutrition: a systematic literature review. Food Nutr Res 56.
    [2] World Health Organization (2004) Iodine status worldwide: WHO Global Database on iodine deficiency. World Health Organization, Switzerland.
    [3] World Health Organization, Micronutrient deficiencies. Available from: https://www.who.int/nutrition/topics/idd/en/.
    [4] Andersson M, Takkouche B, Egli I, et al. (2005) Current global iodine status and progress over the last decade towards the elimination of iodine deficiency. Bull World Health Organ 83: 518–525.
    [5] UNICEF–WHO Joint Committee on Health Policy (1994) World Summit for Children-Mid Decade Goal: Iodine Deficiency Disorders. Geneva, United Nations Children's Fund, World Health Organization; (JCHPSS/94/2.7).
    [6] United Nations Children's Fund (2013) State of the World's Children 2014: Every Child Counts. UNICEF. Available from: http://data.unicef.org/nutrition/iodine.html.
    [7] Andersson M, Karumbunathan V, Zimmermann MB (2012) Global iodine status in 2011 and trends over the past decade. J Nutr 142: 744–750. doi: 10.3945/jn.111.149393
    [8] Prime Minister of Viet Nam (1994) Decision No. 481/TTg of the Prime Minister dated September 1994 on organizing and calling upon all people to buy and use iodized salt.
    [9] Hospital of Endocrinology (2012) Report on Iodine Deficiency Disorder Preventative Activities in Viet Nam.
    [10] Codling K, Quang NV, Phong L, et al (2015) The Rise and Fall of Universal Salt Iodization in Vietnam: Lessons Learned for Designing Sustainable Food Fortification Programs with a Public Health Impact. Food Nutr Bull 36: 441–454. doi: 10.1177/0379572115616039
    [11] Binh Dinh Statistics Office (2016) Binh Dinh Statistical Yearbook 2015, Statistical publishing house 2016.
    [12] Binh Dinh Center for Prevention of Malaria and Endocrine Diseases (2016) Report on the program to prevent iodine deficiency disorders in 2011–2015.
    [13] WHO/UNICEF/ICCIDD (2007) Assessment of iodine deficiency disorders and monitoring their elimination. A Guide for Programme Managers 2007.
    [14] Cambell MJ, Swinsco TDV (2009) Statistics at square one,11th ed. UK: John Wiley & Sons Ltd, 86–95.
    [15] DeMaeyer EM, Lowenstein FW, Thilly CH (1979) The control of endemic goiter, World Health Organization.
    [16] Luc DL (2007) Survey of goiter among students in the last 3 grades of primary school in Dau Tieng district in Binh Duong province in 2005. Full report on scientific topics of the 3rd National Scientific Conference on Endocrinology and Metabolism, Medical Publishing House, Hanoi, 190–199.
    [17] Ta VB, Hoang KH (2007) Study on evaluating the rate of goiter and median urinary iodine in 8 to 10-year-old students in 7 ecological regions in Vietnam in 2005. Report on the full text of scientific projects of the 3rd National Scientific Conference on Endocrinology and Metabolism. Medical Publishing House, 165–175.
    [18] Hoang S (2009) Situation and some factors related to simple goiter prevalence of pupils aged from 8 to 12 in Nam Giang District, Quang Nam Province. Hue Univ J Sci 55: 149–156.
    [19] Tigabu EK, Bekele B, Dachew BA (2017) Prevalence of goiter and associated factors among schoolchildren in northeast Ethiopia. Epidemiol Health 39: e2017055. doi: 10.4178/epih.e2017055
    [20] Mardani M, Gholaami F, Rezapour P, et al. (2015) The Prevalence of Endemic Goiter and Its Relationship to Urine Iodine in Primary School Age Children in Khoramabad. J Health Sci Surveill Sys 3: 71–75.
    [21] Ahmed FBM, Ahmed IEBM (2015) Goiter Incidence Among Basic Schools Children in Shendi Area-Northern Sudan. Sci J Clin Med 4: 117–120. doi: 10.11648/j.sjcm.20150406.11
    [22] Kamath R, Bhat V, Rao RS, et al. (2009) Prevalence of goitre among school children in Belgaum district. Indian J Pediatr 76: 825–828. doi: 10.1007/s12098-009-0101-6
    [23] Nguyen TH, Nguyen TD (2002) Endemic goiter status in one district of Bac Bo Delta. Vietnam J Physiol 6: 17–22.
    [24] Zein AZ, Al-Haithamy S, Obadi Q, et al. (2009) The epidemiology of iodine deficiency disorders in Yemen. Public Health Nutr 3: 245–252.
    [25] Chandra AK, Bhattacharjee A, Malik T, et al. (2008) Goiter prevalence and iodine nutritional status of school Children in a Sub-Himalyan Tarai region of Eastern Uttar Pradesh. India Pediatr 45: 469–474.
    [26] Sharma U, Sharma JP, Sharma A, et al. (2015) Prevalence of goitre among school going children in urban area of Dehradun. Int J Res Med Sci 3: 198–200.
    [27] Gupta RK, Langer B, Raina KS, et al. (2016) Goiter prevalence in school-going children: A cross-sectional study in two border districts of sub-Himalayan Jammu and Kashmir. J Family Med Prim Care 5: 825–828. doi: 10.4103/2249-4863.201150
    [28] Delange F (1994) The disorders induced by iodine deficiency. Thyroid 4: 107–128. doi: 10.1089/thy.1994.4.107
    [29] Truong QD, Le NHG, Nguyen TTL, et al. (2018) The prevalence of malnutrition based on anthropometry among primary schoolchildren in Binh Dinh province, Vietnam in 2016. AIMS Public Health 5: 203–216. doi: 10.3934/publichealth.2018.3.203
    [30] Jimba M, Aitken IW, Joshi AB, et al. (2006) A challenge for monitoring iodine deficiency disorders in rural Nepal. Trop Doct 37:106–107.
    [31] Laillou A, Sophonneary P, Kuong K, et al. (2016) Low Urinary Iodine Concentration among Mothers and Children in Cambodia. Nutrients 8: 172. doi: 10.3390/nu8040172
    [32] Zou Y, Luo X, Ding G, et al. (2014) An assessment of iodine nutritional status and thyroid hormone levels in children aged 8–10 year living in Zhejiang province, China: a cross-sectional study. Eur J Pediatr 173: 929–934. doi: 10.1007/s00431-014-2273-y
    [33] Alsanosy RM, Gaffar AM, Khalafalla HE, et al. (2012) Current iodine nutrition status and progress toward elimination of iodine deficiency disorders in Jazan, Saudi Arabia. BMC Public Health 12: 1006. doi: 10.1186/1471-2458-12-1006
    [34] Leung AM, Braverman LE (2014) Consequences of excess iodine. Nat Rev Endocrinol 10: 136–142. doi: 10.1038/nrendo.2013.251
    [35] Katagiri R, Yuan XY, Kobayashi S, et al. (2017) Effect of excess iodine intake on thyroid diseases in different populations: A systematic review and meta-analyses including observational studies. PLoS One 12: e0173722. doi: 10.1371/journal.pone.0173722
    [36] Manjunath B, Suman G, Hemanth T, et al. (2016) Prevalence and factors associated with goiter among 6–12-year-old children in a Rural area of Karnataka in South India. Biol Trace Elem Res 69: 22–26.
  • This article has been cited by:

    1. Shanshan Gu, Shimin Chai, Chenguang Zhou, Jinhui Zhou, Weak Galerkin finite element method for linear poroelasticity problems, 2023, 190, 01689274, 200, 10.1016/j.apnum.2023.04.015
  • Reader Comments
  • © 2019 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(6353) PDF downloads(1505) Cited by(1)

Figures and Tables

Tables(5)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog