• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1993.04a
• /
• pp.19-26
• /
• 1993

The p-version crack model based on integrals of Legendre polynomial and virtual crack extension method is proposed with its potential for application to stress intensity factor computations in linear elastic fracture mechanics. The main advantage of this model is that the data preparation effort is minimal because only a small number of elements are used and the high accuracy and the rapid rate of convergence can be achieved in the vicinity of crack tip. There are two important findings from this study. Firstly, the limit value, the strain energy of the exact solution can be estimated with successive three p-version approximations by ascertaining the approximations is entered the asymptotic range. Secondly, the rate of convergence of p-version model is almost twice that of h-version model on the basis of uniform or quasiuniform mesh refinement for the cracked panel problem subjected tension.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1143-1156
• /
• 2017

This paper describes recent development activities of the GENESIS code, which is a transport code for heterogeneous three-dimensional geometry, focusing on applications to reactor core analysis. For the treatment of anisotropic scattering, the concept of the simplified Pn method is introduced in order to reduce storage of flux moments. The accuracy of the present method is verified through a benchmark problem. Next, the iteration stability of the GENESIS code for the highly voided condition, which would appear in a severe accident (e.g., design extension) conditions, is discussed. The efficiencies of the coarse mesh finite difference and generalized coarse mesh rebalance acceleration methods are verified with various stabilization techniques. Use of the effective diffusion coefficient and the artificial grid diffusion coefficients are found to be effective to stabilize the acceleration calculation in highly voided conditions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.3
• /
• pp.211-236
• /
• 2019

This paper is devoted to the study and investigation of the Runge-Kutta discontinuous Galerkin method for a system of differential equations consisting of two hyperbolic conservation laws. The numerical coupling flux which is used at a given interface (x = 0) is the upwind flux. Moreover, in the linear case, we derive optimal convergence rates in the $L_2$-norm, showing an error estimate of order ${\mathcal{O}}(h^{k+1})$ in domains where the exact solution is smooth; here h is the mesh width and k is the degree of the (orthogonal Legendre) polynomial functions spanning the finite element subspace. The underlying temporal discretization scheme in time is the third-order total variation diminishing Runge-Kutta scheme. We justify the advantages of the Runge-Kutta discontinuous Galerkin method in a series of numerical examples.

• Journal of the Korean earth science society
• /
• v.34 no.5
• /
• pp.393-401
• /
• 2013

Spherical harmonics power spectrum of the geopotential field of Gaussian-bell type on the sphere was investigated using integral formula that is associated with Legendre polynomials. The geopotential field of Gaussian-bell type is defined as a function of sine of angular distance from the bell's center in order to guarantee the continuity on the global domain. Since the integral-formula associated with the Legendre polynomials was represented with infinite series of polynomial, an estimation method was developed to make the procedure computationally efficient while preserving the accuracy. The spherical harmonics power spectrum was shown to vary significantly depending on the scale parameter of the Gaussian bell. Due to the accurate procedure of the new method, the power (degree variance) spanning over orders that were far higher than machine roundoff was well explored. When the scale parameter (or width) of the Gaussian bell is large, the spectrum drops sharply with the total wavenumber. On the other hand, in case of small scale parameter the spectrum tends to be flat, showing very slow decaying with the total wavenumber. The accuracy of the new method was compared with theoretical values for various scale parameters. The new method was found advantageous over discrete numerical methods, such as Gaussian quadrature and Fourier method, in that it can produce the power spectrum with accuracy and computational efficiency for all range of total wavenumber. The results of present study help to determine the allowable maximum scale parameter of the geopotential field when a Gaussian-bell type is adopted as a localized function.

• Asian-Australasian Journal of Animal Sciences
• /
• v.30 no.10
• /
• pp.1382-1387
• /
• 2017

Objective: During the last decade, genetic evaluation of dairy cows using longitudinal data (test day milk yield or 305-day milk yield) using random regression method has been officially adopted in several countries. The objectives of this study were to estimate covariance functions for genetic and permanent environmental effects and to obtain genetic parameters of 305-day milk yield over seven parities. Methods: Data including 60,279 total 305-day milk yield of 17,309 Iranian Holstein dairy cows in 7 parities calved between 20 to 140 months between 2004 and 2011. Residual variances were modeled by homogeneous and step functions with 7 and 10 classes. Results: The results showed that a third order polynomial for additive genetic and permanent environmental effects plus a step function with 10 classes for the residual variance was the most adequate and parsimonious model to describe the covariance structure of the data. Heritability estimates obtained by this model varied from 0.17 to 0.28. The performance of this model was better than repeatability model. Moreover, 10 classes of residual variance produce the more accurate result than 7 classes or homogeneous residual effect. Conclusion: A quadratic Legendre polynomial for additive genetic and permanent environmental effects with 10 step function residual classes are sufficient to produce a parsimonious model that explained the change in 305-day milk yield over consecutive parities of Iranian Holstein cows.

• Journal of Welding and Joining
• /
• v.30 no.6
• /
• pp.74-79
• /
• 2012

It depends on the joint configuration, dimensions and constraints of the joint whether the residual stress at the root of single-sided butt weld is tensile or not. Therefore, recommendation is generally made that high R ratio should be used in the fatigue test of welded joint in order to prevent excessively long life caused by compressive residual stress. In this research, the residual stress profile in butt weld joint was obtained through compliance method, using successive extension of a slot and measurement of the variation of strain during the slot extension. The residual stress profile was firstly assumed to be the linear summation of Legendre polynomials up to 9th order excluding 0th and 1st order. Strain variation on the surface was measured while the slot was being extended by cutting to find out the 8 unknown coefficients of each polynomial term. The cut was made by the electric discharge machine. It was concluded that the residual stress near the surface is positive valued, however, it turned into the negative value as soon as it passed through 2 or 3 mm of the depth.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2681-2692
• /
• 1996

The finite-volume method for radiation in a three-dimensional non-orthogonal gas turbine combustion chamber with absorbing, emitting and anisotropically scattering medium is presented. The governing radiative transfer equation and its discretization equation using the step scheme are examined, while geometric relations which transform the Cartesian coordinate to a general body-fitted coordinate are provided to close the finite-volume formulation. The scattering phase function is modeled by a Legendre polynomial series. After a benchmark solution for three-dimensional rectangular combustor is obtained to validate the present formulation, a problem in three-dimensional non-orthogonal gas turbine combustor is investigated by changing such parameters as scattering albedo, scattering phase function and optical thickness. Heat flux in case of isotropic scattering is the same as that of non-scattering with specified heat generation in the medium. Forward scattering is found to produce higher radiative heat flux at hot and cold wall than backward scattering and optical thickness is also shown to play an important role in the problem. Results show that finite-volume method for radiation works well in orthogonal and non-orthogonal systems.

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.57-66
• /
• 1993

The p-version crack model based on integrals of Legendre polynomial and virtual crack extension method is proposed with its potential for application to stress intensity factor computations in linear elastic fracture mechanics. The main advantage of this model is that the data preparation effort is minimal because only a small number of elements are used and high accuracy and the rapid convergence can be achieved in the vicinity of crack tip. There are two important findings from this study. Firstly, the limit value, the strain energy of the exact solution, can be estimated with successive three p-version approximations by ascertaining that the approximations enter the asymptotic range. Secondly, the rate of convergence of p-version model is almost twice that of h-version model on the basis of uniform or quasiuniform mesh refinement for the cracked panel problem subjected to tension.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.151-160
• /
• 1993

This paper is concerned with formulations of the hierarchical $C^{\circ}$-plate element on the basis of Reissner-Mindlin plate theory. On reason for the development of the aforementioned element based on Integrals of Legendre shape functions is that it is still difficult to construct elements based on h-version concepts which are accurate and stable against the shear locking effects. An adaptive mesh refinement and selective p-distribution of the polynomial degree using hp-version of the finite element method are proposed to verify the superior convergence and algorithmic efficiency with the help of the simply supported L-shaped plate problems.

• Asian-Australasian Journal of Animal Sciences
• /
• v.17 no.7
• /
• pp.914-918
• /
• 2004

The random regression model methodology was applied into the estimation of genetic parameters for body weights in Chinese Simmental cattle to replace the traditional multiple trait models. The variance components were estimated using Gibbs sampling procedure on Bayesion theory. The data were extracted for Chinese Simmental cattle born during 1980 to 2000 from 6 national breeding farms, where records from 3 months to 36 months were only used in this study. A 3 orders Legendre polynomial was defined as the submodel to describe the general law of that body weight changing with months of age in population. The heritabilities of body weights from 3 months to 36 months varied between 0.31 and 0.48, where the heritabilities from 3 months to 12 months slightly decreased with months of age but ones from 13 months to 36 months increased with months of age. Specially, the heritabilities at eighteenth and twenty-fourth month of age were 0.33 and 0.36, respectively, which were slightly greater than 0.30 and 0.31 from multiple trait models. In addition, the genetic and phenotypic correlations between body weights at different month ages were also obtained using regression model.