• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.359-366
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• Journal of the Korean Wood Science and Technology
• /
• v.35 no.3
• /
• pp.1-9
• /
• 2007

To analyze the behavior of multiple bolted connection with analytical method, there was a major problem in terms of the number of bolt as two. Multiple bolted connection is wrong in other way, that is impossible to analyze the behavior of connection with various types. So the numerical analysis model was developed to analyze the behavior of connection with various types, and to overcome the limitation of number of bolt. The optimum values for analysis the behavior of multiple bolted connection with numerical methods were obtained as, r=0.35, Kx=137.5, and Ky=257.4. The numerical analysis method was developed in this study showed high efficiency for the existing methods in the behavior analysis of connection.

• Structural Engineering and Mechanics
• /
• v.49 no.5
• /
• pp.555-568
• /
• 2014

In this article, stability of composite conical shells subjected to dynamic external pressure is investigated by numerical and experimental methods. In experimental tests, cross-ply glass woven fabrics were selected for manufacturing of specimens. Hand-layup method was employed for fabricating the glass-epoxy composite shells. A test-setup that includes pressure vessel and data acquisition system was designed. Also, numerical analyses are performed. In these analyses, effect of actual geometrical imperfections of experimental specimens on the numerical results is investigated. For introducing the imperfections to the numerical models, linear eigen-value buckling analyses were employed. The buckling modes are multiplied by very small numbers that are derived from measurement of actual specimens. Finally, results are compared together while a good agreement between results of imperfect numerical analyses and experimental tests is observed.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.5-9
• /
• 2010

A numerical wind tunnel simulation is performed in order to predict wind loads acting on a building. The aim of the present study is to suggest a guideline for the numerical wind tunnel analysis, which could provide more detail wind load distributions compared to the wind code and expensive wind tunnel experiments. To validate the present numerical simulation, wind-induced loads on a 6 m cube model is predicted. Atmospheric boundary layer is used as a inlet boundary condition. Various effect of numerical methods are investigated such as size of computational domain, grid density, turbulence model and discretization scheme. The appropriate procedure for the numerical wind tunnel analysis is suggested through the present study.

• Journal of Industrial Technology
• /
• v.21 no.A
• /
• pp.7-14
• /
• 2001

In building structure, the story height can be minimized by providing openings in beams which serves for the utility equipments passing through. The dead space in false ceiling thus put to economical use in the form of a substantial reduction in materials and construction cost. In the case of steel structure, there is no critical risk in the structural strength because of reinforcing methods of stiffness and steel plate but in the case of reinforced concrete structure, proper provision should be made in designing these openings, otherwise there is a risk that these opening will possibly weaken the structural strength of the building frame to a critical degree. In this paper, for the numerical analysis of the reinforced concrete beams with circular opening in the web, expecting stress concentration of the circular opening, reinforcing methods were studied. Twenty test pieces with each different reinforcing methods were tested and their resisting forces were defined. From the numerical analysis and test results, the followings were founded;(1)high shear stress distributed around the openings reduce the shearing strength, (2)from the numerical analysis, the maximum tensile stress occurred at opening nodes 1,7, these phenomena were agreed with the test results, (3)reinforcing method around openings have to carried out for stopping diagonal cracks, and (4)both, by steel plate, and wire mesh, are effective reinforcing methods.

• International Journal of Ocean System Engineering
• /
• v.1 no.3
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.

• The Journal of Engineering Geology
• /
• v.33 no.3
• /
• pp.439-459
• /
• 2023

Numerical analysis has been adopted in nearly all modern scientific and engineering fields due to the rapid and ongoing evolution of computational technology, with the number of grid or mesh points in a given data field also increasing. Some values must be extracted from large data fields to evaluate and supplement numerical analysis results and observational data, thereby highlighting the need for a fast and effective interpolation approach. The quadrilateral irregular network (QIN) proposed in this study is a fast and reliable interpolation method that is capable of sufficiently satisfying these demands. A comparative sensitivity analysis is first performed using known test functions to assess the accuracy and computational requirements of QIN relative to conventional interpolation methods. These same interpolation methods are then employed to produce simple numerical model results for a real-world comparison. Unlike conventional interpolation methods, QIN can obtain reliable results with a guaranteed degree of accuracy since there is no need to determine the optimal parameter values. Furthermore, QIN is a computationally efficient method compared with conventional interpolation methods that require the entire data space to be evaluated during interpolation, even if only a subset of the data space requires interpolation.

• Journal of applied mathematics & informatics
• /
• v.26 no.1_2
• /
• pp.135-149
• /
• 2008

Many deterministic systems are described by Ordinary differential equations and can often be improved by including stochastic effects, but numerical methods for solving stochastic differential equations(SDEs) are required, and work in this area is far less advanced than for deterministic differential equations. In this paper,first we follow [7] to describe Runge-Kutta methods with order 2 from Taylor approximations in the weak sense and present two well known Runge-Kutta methods, RK2-TO and RK2-PL. Then we obtain a new 3-stage explicit Runge-Kutta with order 2 in weak sense and compare the numerical results among these three methods.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.325-343
• /
• 2012

We present different numerical methods for solving the shallow shelf equations with basal drag (SSAB). An alternative approach of splitting the SSAB equation into a Laplacian and diagonal shift operator is discussed with respect to the underlying eigenvalue problem. First, we solve the equations using standard methods. Then, the coupled equations are decomposed into operators for membranes stresses, basal shear stress and driving stress. Applying reasonable parameter values, we demonstrate that the operator of the membrane stresses is much stiffer than the operator of the basal shear stress. Here, we could apply a new splitting method, which alternates between the iteration on the membrane-stress operator and the basal-shear operator, with a more frequent iteration on the operator of the membrane stresses. We show that this splitting accelerates and stabilize the computational performance of the numerical method, although an appropriate choice of the standard method used to solve for all operators in one step speeds up the scheme as well.

• Journal of applied mathematics & informatics
• /
• v.33 no.1_2
• /
• pp.61-76
• /
• 2015

In this paper, a fifth order extension of Potra's third order iterative method is proposed for solving nonlinear equations. A convergence theorem along with the error bounds is established. The method takes three functions and one derivative evaluations giving its efficiency index equals to 1.495. Some numerical examples are also solved and the results obtained are compared with some other existing fifth order methods. Next, the interval extension of both third and fifth order Potra's method are developed by using the concepts of interval analysis. Convergence analysis of these methods are discussed to establish their third and fifth orders respectively. A number of numerical examples are worked out using INTLAB in order to demonstrate the efficacy of the methods. The results of the proposed methods are compared with the results of the interval Newton method.