• Kyungpook Mathematical Journal
• /
• v.47 no.4
• /
• pp.529-548
• /
• 2007

The purpose of this article is to find a relation between the finite difference method and the boundary element method, and propose a new approach deriving a discrete approximation formula as like that of the finite difference method for harmonic functions. We develop a discrete approximation formula on a uniform grid based on the boundary integral formulations. We consider three different boundary integral formulations and derive one discrete approximation formula on the uniform grid for the harmonic function. We show that the proposed discrete approximation formula has the same computational molecules with that of the finite difference formula for the Laplace operator ${\nabla}^2$.

• Proceedings of the KIEE Conference
• /
• 1997.07a
• /
• pp.309-311
• /
• 1997

In this paper, it is proposed that a algorithm which is applicable to the transient analysis by combined use of the Laplace transform and the finite element method. The proposed method removes the time terms using the Laplace transform and then solves the associated equation with the finite element method. The solution which is solved at frequency domain is transformed into time domain by use of the Laplace inversion. To verify proposed algorithm, heat conduction problem is analysed and found a good agreement with analytic solution.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2000.10a
• /
• pp.51-58
• /
• 2000

Dynamic relaxation method is a shape finding analysis method for flexible structures by introducing the dynamic equilibrium equation. However, it is difficult for shape finding to estimate the most appropriate values for the mass and damping on each shape because the values are random one. In this study, the unit mass, the unit damping and the principal direction stiffness are utilized to avoid the random values, and the Newmarks assumption is introduced during the dynamic analysis. By introducing variant time increment method presented, the convergence time is reduced, that is, it can be reduced the total times for analysis.

• International Journal of Safety
• /
• v.9 no.2
• /
• pp.6-10
• /
• 2010

The flash point of liquid solution is one of the most important flammability properties that used in hazard and risk assessments. Minimum flash point behaviour (MFPB) is showed when the flash point of a liquid mixture is below the flash points of the individual components. In this paper, the lower flash points for the flammable binary system, n-decane+n-octanol, were measured by Pensky-Martens closed cup tester. This binary mixture exhibited MFPB. The measured flash points were compared with the values calculated by the Raoult's law and the optimization method using van Laar and UNIQUAC equations. The optimization method were found to be better than those based on the Raoult's law, and successfully estimated MFPB. The opimization method based on the van Laar equation described the experimentally-derived data more effectively than was the case when the prediction model was based upon the UNIQUAC.

• Journal of Ship and Ocean Technology
• /
• v.5 no.1
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.3
• /
• pp.610-625
• /
• 2015

In this paper the wave run-up around a circular column in regular waves is numerically calculated to investigate the applicability of the Modified Marker-Density (MMD) method to prediction of wave run-up around an offshore platform. The MMD method is one of the methods to define the highly nonlinear free surface. The governing equations are the Navier-Stokes equations and the continuity equation which are computed in Cartesian grid system. To validate incident waves generated by numerical simulation, those are compared with the solutions of the Stokes $5^{th}$ order wave theory. The wave run-up simulations are performed varying the steepness and period of incident waves as referred experimental data. The numerical results are compared to the experimental data and the results show good agreements.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.11 no.2
• /
• pp.1644-1650
• /
• 1969

In the design of general hydrological structures, it is well know that the design flood is of importance in the design of those structures. As the design flood is estimated using the design storm, the design storm is defined by the rainfall intensity itself. Though I had studied and reported many times the reports about the rainfall-intensity in my country, poorly I did not study the long-period variation of the intensity through each section in my country before. But now, in the basin area of the Han river and the Keum river, the self-recorded rainfall charts of the single storms, which are mostly above rainfall amount of 30mm and data of about 4500 with the 150 stationyear, were analyzed, And then, the intensity formula of the hourly unit is estimated using the period from 10 minutes to 5 days. The method to analyze and estimate them, and the final results will be summarized as mentioned below: (i) At first I intended to select out the homogeneous watersheds of three, one in the Han river and two in the Keum river. But I would select the northern and the sourthern river basins, and westward from Koan station, in the basins of the Han river. Also I would select the upstream area, and the downstream area including the watershed of Chungioo, Kongjoo, Chupungryung, and the Mt. Sock, in the basins of the Keum river. Finally, I could find that there couldn't in the Keum river basin. So, I decided out and analyze only river basins of the Han river with limitation mentioned above. (ii) The statistical method to select out the homogenous watersheds is the test of homogeneous variance, and it is estimated from the following equation: $$X_{k1}^2=[{\Sigma}(n_i-1)log\bar{S^2}-\Sigma(n_i-1)log\bar{S^2}]{\times}loge$$ (iii) Actually, each homogeneous watershed has individually its own intensity formula, But I would express them as the actual amount, because the equation of intensity variance is experiential and theoretical equation of the variance. Therefore the caluating equation is actually more convenient in the actual uses. (iv) This report is one of the series for me to give the basis to the actual designs. The cost for this study is provided by the Ministry of Construction. And the designs of the hydrological structures in the watersheds with limitation mentioned above may be concerned with and based upon this report.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• International Journal of Safety
• /
• v.1 no.1
• /
• pp.36-42
• /
• 2002

This paper is aimed at predicting the life of rubber insulating gloves under normal operating stresses from relatively rapid test performed at higher stresses. Specimens of rubber insulating gloves are subject to multiple stress conditions, i.e. combined electrical and thermal stresses. Two modes of electrical stress, step voltage stress and constant voltage stress are used in specimen aging. There are two types of test for electrical stress in this experiment: the one is Breakdown Voltage (BDV) test under step voltage stress and thermal stress and the other is lifetime test under constant voltage stress and temperature stress. The ac breakdown voltage defined as the break-down point of insulation that leakage current excesses a limit value, l0mA in this experiment, is determined. Because the very high variability of aging data requires the application of statistical model, Weibull distribution is used to represent the failure times as the straight line on Weibull probability paper. Weibull parameters are deter-mined by three statistical methods i.e. maximum likelihood method, graphical method and least squares method, which employ SAS package, Weibull probability paper and FORTRAN, respectively. Two chosen models for predicting the life under simultaneous electrical and thermal stresses are inverse power model and exponential model. And the constants of life equation for multistress aging are calculated using numerical method, such as Gauss Jordan method etc.. The completion of life equation enables to estimate the life at normal stress based on the data collected from accelerated aging test. Also the comparison of the calculated lifetimes between the inverse power model and the exponential model is carried out. And the lifetimes calculated by three statistical methods with lower voltage than test voltage are compared. The results obtained from the suggested experimental method are presented and discussed.

• Journal of Korea Soil Environment Society
• /
• v.5 no.1
• /
• pp.3-12
• /
• 2000

The multi-region model, to describe preferential flow, is an equation representing solute transport in soils by dividing soil into numerous pore groups and using the hydraulic properties of the soil. As the model partial differential equation (PDE) is solved numerically with finite difference methods. a modified equivalent partial differential equation(MEPDE) of the partial differential equation of the multi-region model is derived to analyze the accuracy and consistency of the solution of the model PDE and the Von Neumann method is used to analyze the stability of the finite difference scheme. The evaluation obtained from the MEPDE indicated that the finite difference scheme was found to be consistent with the model PDE and had the second order accuracy The stability analysis is performed to analyze the model PDE with the amplification ratio and the phase lag using the Von Neumann method. The amplification ratio of the finite difference scheme gave non-dissipative results with various Peclet numbers and yielded the most high values as the Peclet number was one. The phase lag showed that the frequency component of the finite difference scheme lagged the true solution. From the result of the stability analysis for the model PDE, it is analyzed that the model domain should be discretized in the range of Pe < 1.0 and Cr < 2.0 to obtain the more accurate solution.