• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.413-420
• /
• 1994

Dynamics of the river plume is a very complicated non-linear problem with the free boundary changing in time and space. Mixing with the ambient water through the boundary makes the problem more complicated. In this paper we reduced 3-dimensional problem into 1-dimensional one by using the integral analysis method. Basic equations have been integrated over the lateral and vertical variations. For these integrations we adopted the well-established assumption that the flow-axis component of plume velocity and the density difference of the plume with the ambient water have Gaussian distributions in directions which are perpendicular to the flow-axis of the plume. We also used the result of our previous study on the lateral spreading velocity of the plume derived under the same assumption. And entrainment was included as a mixing process. The resultant 1-dimensional equations were solved by Runge-Kutta numerical method. Consequently, comparatively easy method of numerical analysis is presented for the 3-dimensional river plume. The method can also be used for the analysis of the thermal plume of cooling water of power plants.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.1A
• /
• pp.11-18
• /
• 2012

This paper deals with the numerical determination of the stress intensity factors of cracked aluminum plates under the mixed mode of $K_I$ and $K_{II}$ in glass-epoxy fiber reinforced composites. For the stress intensity factors, two different models are reviewed such as VCCT and two-step extension method. The p-convergent partial layerwise model is adopted to determine the fracture parameters in terms of energy release rates and stress intensity factors. The p-convergent approach is based on the concept of subparametric element. In assumed displacement field, strain-displacement relations and 3-D constitutive equations of a layer are obtained by combination of 2-D and 1-D higher-order shape functions. In the elements, Lobatto shape functions and Gauss-Lobatto technique are employed to interpolate displacement fields and to implement numerical quadrature. Using the models and techniques considered, effects of composite laminate configuration according to inclined angles and adhesive properties on the performance of bonded composite patch are investigated. In addition to these, the out-of-plane bending effect has been investigated across the thickness of patch repaired laminate plates due to the change of neutral axis. The present model provides accuracy and simplicity in terms of stress intensity factors, stress distribution, number of degrees of freedom, and energy release rates as compared with previous works in literatures.

• 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.

• The Korean Journal of Applied Statistics
• /
• v.23 no.5
• /
• pp.871-882
• /
• 2010

In animal tumorigenicity data, the occurrence time of tumor is not observed because the existence of a tumor is examined only at either time of natural death or time of sacrifice for the animal. A three-state model (Health-Tumor onset-Death) is widely used to model the incomplete data. In this paper, we employed a frailty effect into the three-state model to incorporate the dependency of death on tumor occurrence when the time of natural death works as an informative censoring against the tumor onset time. For the inference of parameters, then the EM algorithm is considered in order to deal with missing quantities of tumor onset time and random frailty. The proposed method is applied to the bladder tumor data taken from Lindsey and Ryan (1993, 1994) and a simulation study is performed to show the behavior of the proposed estimators.

• 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.

• The Korean Journal of Nuclear Medicine
• /
• v.32 no.3
• /
• pp.225-237
• /
• 1998

Purpose: We investigated the statistical methods to compose the functional brain map of human working memory and the principal factors that have an effect on the methods for localization. Materials and Methods: Repeated PET scans with successive four tasks, which consist of one control and three different activation tasks, were performed on six right-handed normal volunteers for 2 minutes after bolus injections of 925 MBq $H_2^{15}O$ at the intervals of 30 minutes. Image data were analyzed using SPM96 (Statistical Parametric Mapping) implemented with Matlab (Mathworks Inc., U.S.A.). Images from the same subject were spatially registered and were normalized using linear and nonlinear transformation methods. Significant difference between control and each activation state was estimated at every voxel based on the general linear model. Differences of global counts were removed using analysis of covariance (ANCOVA) with global activity as covariate. Using the mean and variance for each condition which was adjusted using ANCOVA, t-statistics was performed on every voxel To interpret the results more easily, t-values were transformed to the standard Gaussian distribution (Z-score). Results: All the subjects carried out the activation and control tests successfully. Average rate of correct answers was 95%. The numbers of activated blobs were 4 for verbal memory I, 9 for verbal memory II, 9 for visual memory, and 6 for conjunctive activation of these three tasks. The verbal working memory activates predominantly left-sided structures, and the visual memory activates the right hemisphere. Conclusion: We conclude that rCBF PET imaging and statistical parametric mapping method were useful in the localization of the brain regions for verbal and visual working memory.