• 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 Earthquake Engineering Society of Korea
• /
• v.10 no.4 s.50
• /
• pp.25-33
• /
• 2006

The purpose of this study is to develop the curve fitted equations for practicality and actual calculation during seismic performance evaluation of buried pipelines. Curve fitting for strain curve according to the wavelength of the seismic wave was produced using the non-linear least square method and the equations with the best results was suggested. In addition, a degree and coefficient of polynomial fitting equation needed to use curve fitted equation were identified. Interpreting process during the test of resistance of earthquake of buried pipelines with various end boundary conditions were provided through example questions. The results of this study were used to conduct a dynamic response analysis and a seismic performance evaluation of concrete, steel, and FRP pipes with various end boundary conditions.

• Journal for History of Mathematics
• /
• v.30 no.4
• /
• pp.203-219
• /
• 2017

As is well known, the most important development in the history of Chinese mathematics is materialized in Song-Yuan era through tianyuanshu up to siyuanshu for constructing equations and zengcheng kaifangfa for solving them. There are only two authors in the period, Li Ye and Zhu Shijie who left works dealing with them. They were almost forgotten until the late 18th century in China but Zhu's Suanxue Qimeng(1299) had been a main reference for the Joseon mathematics. Commentary by Luo Shilin on Zhu's Siyuan Yujian(1303) was brought into Joseon in the mid-19th century which induced a great attention to Joseon mathematicians with a thorough understanding of Zhu's tianyuanshu. We discuss the history that Joseon mathematicians succeeded to obtain the mathematical structures of Siyuan Yujian based on the Zhu's tianyuanshu.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.198-198
• /
• 2011

The numerical simulation of dam break problem suffers from several challenges in terms of accuracy, stability, and versatility of the simulation algorithm since the water flow is generally discontinuous and presents abrupt variations. Thus, to obtain stable and accurate solutions, flow models for this purpose require numerical schemes provided with shock-capturing properties, and with the ability to work with flexible two-dimensional meshes. In this context, SU/PG method(Hughes and Brooks, 1979) is excellent candidate for the solution of the dam break problem. The weak formulation of the equations and the discontinuous polynomial basis lead to an accurate representation of bore waves(shocks). Furthermore, the discretization of the domain in finite elements is extremely effective in modeling complex geometries. In this study, a finite element model based on the SU/PG scheme is developed to solve shallow water equations and the model is applied to dam break problem. It is found that the present model accurately captures the bore wave that propagates downstream while spreading laterally and the depression wave that moves upstream. Furthermore, the propagation and formation of water surface profile compared favorably with those obtained by the previously published results.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.08a
• /
• pp.308-312
• /
• 1999

Developed the model equations which calculate roll force, roll power during hot rolling in real time. The variables which mainly effect on the roll force, roll power are shape factor, reduction, roll diameter, roll velocity, strip inlet temperature, carbon content of strip and strip-roll contact friction coefficient. Among these variables roll diameter, roll velocity, inlet temperature, carbon content and friction coefficient can be excluded in interpolated model equation by introducing equation of die force(F'), power(p') of the frictionless uniform plane strain compression which can be calculated without iteration. At the case of coulomb friction coefficient of 0.3, we evaluated coefficient of polynomial equations of {{{{ { F} over {F' } }}}}, {{{{ { Pf} over {Pd }, { Pd} over {P' } }}}} from the result of finite element analysis using interpolation. It was found that the change of values of {{{{ { F} over {F' }, { P} over {P' } }}}} with the friction coefficient tend to straight line which slope depend only on shape factor. With these properties, developed model equations could be extended to other values of coulomb friction coefficient. To verify developed roll force, roll power model equation we compared the results from these model equation with the results from these model equation with the results from finite element analysis in factory process condition.

• Journal of the Society of Disaster Information
• /
• v.10 no.1
• /
• pp.123-140
• /
• 2014

This paper proposes a simple and effective method for determining the dynamic characteristics of revolution shells. This is a weighted residual method in which the collocation points are taken at the roots of orthogonal polynomial. In this paper the collocation method is employed to replace a partical differential eqations by a system of ordinary differential equations in time, and the resulting equations are solved by two different numerical methods of time integration : an implicit method and an explicit method. The proposed approach is formulated in some detail. The versatility and accuracy are illustrated through several numerical examples. The method appears to be relatively easy to set up and gives satisfactory results.

• Journal of the Korean Society of Safety
• /
• v.33 no.4
• /
• pp.21-29
• /
• 2018

In many companies handling flammable liquids, explosion-proof electrical equipment have been installed according to the Korean Industrial Standards (KS C IEC 60079-10-1). In these standards, hazardous area for explosive gas atmospheres has to be classified by the evaluation of the evaporation rate of flammable liquid leakage. The evaporation rate is an important factor to determine the zones classification and hazardous area distance. However, there is no systematic method or rule for the estimation of evaporation rate in these standards and the first principle equations of a evaporation rate are very difficult. Thus, it is really hard for industrial workplaces to employ these equations. Thus, this problem can trigger inaccurate results for evaluating evaporation range. In this study, empirical models for estimating an evaporation rate of flammable liquid have been developed to tackle this problem. Throughout the sensitivity analysis of the first principle equations, it can be found that main factors for the evaporation rate are wind speed and temperature and empirical models have to be nonlinear. Polynomial regression is employed to build empirical models. Methanol, benzene, para-xylene and toluene are selected as case studies to verify the accuracy of empirical models.

• Environmental and Resource Economics Review
• /
• v.13 no.4
• /
• pp.717-734
• /
• 2004

This paper examined the statistical goodness-of-fit tests for biological growth model in bioeconomic analysis. Some authors estimated usually growth function for fish in the world. However, few studies have estimated growth equations for the bivalve species. Thus, this paper studied the common functional forms of fitting growth equations for cham scallops considering environmental factors and production styles. The following functional forms are considered: linear, log-reciprocal, double log, polynomial and linear with interactions. Results of fitting these various functional forms with real data are compared and evaluated using standard statistical goodness-of-fit tests. Results also indicate that log-reciprocal function is statistically the best fit to the real data. Therefore, the log-reciprocal function is decided the best function describing cham scallop biological growth and hence might be useful for economic evaluation(i.e., optimal harvesting time).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.493-499
• /
• 2007

The Taylor expansion expresses a differentiable function and its coefficients provide good approximations for the given function and its derivatives. In this study, m-th order Taylor Polynomial is constructed and the coefficients are computed by the Moving Least Squares method. The coefficients are applied to the governing partial differential equation for solid problems including crack problems. The discrete system of difference equations are set up based on the concept of point collocation. The developed method effectively overcomes the shortcomings of the finite difference method which is dependent of the grid structure and has no approximation function, and the Galerkin-based meshfree method which involves time-consuming integration of weak form and differentiation of the shape function and cumbersome treatment of essential boundary.

• Korea-Australia Rheology Journal
• /
• v.16 no.4
• /
• pp.183-191
• /
• 2004

High Deborah or Weissenberg number problems in viscoelastic flow modeling have been known formidably difficult even in the inertialess limit. There exists almost no result that shows satisfactory accuracy and proper mesh convergence at the same time. However recently, quite a breakthrough seems to have been made in this field of computational rheology. So called matrix-logarithm (here we name it tensor-logarithm) formulation of the viscoelastic constitutive equations originally written in terms of the conformation tensor has been suggested by Fattal and Kupferman (2004) and its finite element implementation has been first presented by Hulsen (2004). Both the works have reported almost unbounded convergence limit in solving two benchmark problems. This new formulation incorporates proper polynomial interpolations of the log­arithm for the variables that exhibit steep exponential dependence near stagnation points, and it also strictly preserves the positive definiteness of the conformation tensor. In this study, we present an alternative pro­cedure for deriving the tensor-logarithmic representation of the differential constitutive equations and pro­vide a numerical example with the Leonov model in 4:1 planar contraction flows. Dramatic improvement of the computational algorithm with stable convergence has been demonstrated and it seems that there exists appropriate mesh convergence even though this conclusion requires further study. It is thought that this new formalism will work only for a few differential constitutive equations proven globally stable. Thus the math­ematical stability criteria perhaps play an important role on the choice and development of the suitable con­stitutive equations. In this respect, the Leonov viscoelastic model is quite feasible and becomes more essential since it has been proven globally stable and it offers the simplest form in the tensor-logarithmic formulation.