• Korean Journal of Mathematics
• /
• v.22 no.3
• /
• pp.429-441
• /
• 2014

In this work we study two a posteriori error estimators of hierarchical type for lowest-order mixed finite element methods. One estimator is computed by solving a global defect problem based on the splitting of the lowest-order Brezzi-Douglas-Marini space, and the other estimator is locally computable by applying the standard localization to the first estimator. We establish the reliability and efficiency of both estimators by comparing them with the standard residual estimator. In addition, it is shown that the error estimator based on the global defect problem is asymptotically exact under suitable conditions.

• Structural Engineering and Mechanics
• /
• v.8 no.5
• /
• pp.513-529
• /
• 1999

A concept of hierarchical modeling, the newest modeling technology, has been introduced in early 1990's. This new technology has a great potential to advance the capabilities of current computational mechanics. A first step to implement this concept is to construct hierarchical models, a family of mathematical models sequentially connected by a key parameter of the problem under consideration and have different levels in modeling accuracy, and to investigate characteristics in their numerical simulation aspects. Among representative model problems to explore this concept are elastic structures such as beam-, arch-, plate- and shell-like structures because the mechanical behavior through the thickness can be approximated with sequential accuracy by varying the order of thickness polynomials in the displacement or stress fields. But, in the numerical, analysis of hierarchical models, two kinds of errors prevail, the modeling error and the numerical approximation error. To ensure numerical simulation quality, an accurate estimation of these two errors is definitely essential. Here, a local a posteriori error estimator for elastic structures with thin domain such as plate- and shell-like structures is derived using the element residuals and the flux balancing technique. This method guarantees upper bounds for the global error, and also provides accurate local error indicators for two types of errors, in the energy norm. Compared to the classical error estimators using the flux averaging technique, this shows considerably reliable and accurate effectivity indices. To illustrate the theoretical results and to verify the validity of the proposed error estimator, representative numerical examples are provided.

• Journal of Theoretical and Applied Mechanics
• /
• v.3 no.1
• /
• pp.16-33
• /
• 2002

A concept of hierarchical modeling, the newest modeling technology. has been introduced early In 1990. This nu technology has a goat potential to advance the capabilities of current computational mechanics. A first step to Implement this concept is to construct hierarchical models, a family of mathematical models which are sequentially connected by a key parameter of the problem under consideration and have different levels in modeling accuracy, and to investigate characteristics In their numerical simulation aspects. Among representative model problems to explore this concept are elastic structures such as beam-, arch-. plate- and shell-like structures because the mechanical behavior through the thickness can be approximated with sequential accuracy by varying the order of thickness polynomials in the displacement or stress fields. But, in the numerical analysis of hierarchical models, two kinds of errors prevail: the modeling error and the numerical approximation errors. To ensure numerical simulation quality, an accurate estimation of these two errors Is definitely essential. Here, a local a posteriori error estimator for elastic structures with thin domain such as plate- and shell-like structures Is derived using element residuals and flux balancing technique. This method guarantees upper bounds for the global error, and also provides accurate local error Indicators for two types of errors, in the energy norm. Comparing to the classical error estimators using flux averaging technique, this shows considerably reliable and accurate effectivity indices. To illustrate the theoretical results and to verify the validity of the proposed error estimator, representative numerical examples are provided.

• Journal of the Korean Statistical Society
• /
• v.27 no.1
• /
• pp.1-10
• /
• 1998

In the problem of estimating a vector of unknown regression coefficients under the sum of squared error losses in a hierarchical linear model, we propose the hierarchical Bayes estimator of a vector of unknown regression coefficients in a hierarchical linear model, and then prove the admissibility of this estimator using Blyth's (196\51) method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.139-170
• /
• 2013

We extend a p-hierarchical decomposition of the second degree finite element space of N$\acute{e}$d$\acute{e}$lec for tetrahedral meshes in three dimensions given in [1] to meshes with hexahedral elements, and derive p-hierarchical decompositions of the second degree finite element space of Raviart-Thomas in two dimensions for triangular and quadrilateral meshes. After having proved stability of these subspace decompositions and requiring certain saturation assumptions to hold, we construct a local a posteriori error estimator for fem and bem coupling of a time-harmonic electromagnetic eddy current problem in $\mathbb{R}^3$. We perform some numerical tests to underline reliability and efficiency of the estimator and test its usefulness in an adaptive refinement scheme.

• Journal of the Korean Data and Information Science Society
• /
• v.20 no.5
• /
• pp.895-903
• /
• 2009

The penalized partial likelihood based on restricted maximum likelihood method has been widely used for the inference of frailty models. However, the standard-error estimate for frailty parameter estimator can be downwardly biased. In this paper we show that such underestimation can be corrected by using hierarchical likelihood. In particular, the hierarchical likelihood gives a statistically efficient procedure for various random-effect models including frailty models. The proposed method is illustrated via a numerical example and simulation study. The simulation results demonstrate that the corrected standard-error estimate largely improves such bias.

• Communications for Statistical Applications and Methods
• /
• v.3 no.2
• /
• pp.205-216
• /
• 1996

Consider the problem of estimating a multivariate normal mean with an unknown covarience matrix under a weighted sum of squared error losses. We first provide hierarchical Bayes estimators which shrink the usual (maximum liklihood, uniformly minimum variance unbiased) estimator towards a regression surface and then prove the admissibility of these estimators using Blyth's (1951) method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.481-486
• /
• 2007

The Zienkiewicz-Zhu(Z/Z) error estimate is slightly modified for the hierarchical p-refinement, and is then applied to L-shaped plates subjected to bending to demonstrate its effectiveness. An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the superconvergent patch recovery(SPR) technique. The modified Z/Z error estimate p-refinement is different from the conventional approach because the high order shape functions based on integrals of Legendre polynomials are used to interpolate displacements within an element, on the other hand, the same order of basis function based on Pascal's triangle tree is also used to interpolate recovered stresses. The least-square method is used to fit a polynomial to the stresses computed at the sampling points. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly or selectively. It is noted that the error decreases rapidly with an increase in the number of degrees of freedom and the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.533-541
• /
• 2007

The Zienkiewicz-Zhu(Z/Z) error estimate is slightly modified for the hierarchical p-refinement, and is then applied to L-shaped plates subjected to bending to demonstrate its effectiveness. An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the superconvergent patch recovery(SPR) technique. The modified Z/Z error estimate p-refinement is different from the conventional approach because the high order shape functions based on integrals of Legendre polynomials are used to interpolate displacements within an element, on the other hand, the same order of basis function based on Pascal's triangle tree is also used to interpolate recovered stresses. The least-square method is used to fit a polynomial to the stresses computed at the sampling points. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly or selectively. It is noted that the error decreases rapidly with an increase in the number of degrees of freedom and the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.4 s.74
• /
• pp.429-440
• /
• 2006

This paper comprises two specific objectives. The first is to examine the applicability of ordinary kriging interpolation(OK) to the p-adaptivity of the finite element method that is based on variogram modeling. The second objective Is to present the adaptive procedure by the hierarchical p-refinement in conjunction with a posteriori error estimator using the modified S.P.R. (superconvergent patch recovery) method. The ordinary kriging method that is one of weighted interpolation techniques is applied to obtain the estimated exact solution from the stress data at the Gauss points. The weight factor is determined by experimental and theoretical variograms for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. To verify the performance of the modified S.P.R. method, the new error estimator based on limit value has been proposed. The validity of the proposed approach has been tested with the help of some benchmark problems of linear elastic fracture mechanics such as a centrally cracked panel, a single edged crack, and a double edged crack.