• 한국수자원학회논문집
• /
• 제41권3호
• /
• pp.243-249
• /
• 2008

본 연구에서는 지표흐름 습지의 수질관리를 위하여 습지 내 오염물질을 적절히 모의할 수 있는 연속교반탱크반응기(CSTR) 모델을 구축하였다. 이 모델은 4차 Runge-Kutta법을 사용하고, 실측치와 계산치의 차이를 최소화 하는 최적화 기법으로 해를 구하며, 미국 EPA 습지 데이터베이스에 수록된 자료중 수질 및 수리자료가 충분하며, 분석이 용이한 하나의 수생대를 갖는 습지를 선택하여, 수질항목중 관심도가 높은 BOD, TSS (Total Suspended Solid)의 모의에 적용하였다. 습지의 체적이 일정할 경우와 체적 및 유량이 일정할 경우 각각에 대해 반응속도상수만을 고려할 경우, 물질흐름양만을 고려할 경우 및 반응속도상수와 물질흐름양 모두를 고려할 경우로 나누어 모의를 실시하였다. 모의 결과 기존의 반응속도상수만을 고려한 경우보다 반응속도상수 및 물질흐름양 모두를 고려하는 것이 실제현상을 잘 나타냈으며, 습지의 TSS 침강속도는 $0.3{\sim}0.4\;m/d$의 값을 나타내는 것으로 파악되었다. 이 모델은 습지의 수질관리에 적절히 활용할 수 있을 것으로 판단된다.

• Geomechanics and Engineering
• /
• 제13권4호
• /
• pp.585-600
• /
• 2017

This study deals with simple solutions for a spherical or circular opening excavated in elastic-brittle plastic rock mass compatible with a linear Mohr-Coulomb (M-C) or a nonlinear Hoek-Brown (H-B) yield criterion. Based on total strain approach, the closed-form solutions of stresses and displacement are derived simultaneously for circular and spherical openings using original H-B and M-C yield criteria. Two simple numerical procedures are proposed for the solution of generalized H-B and M-C yield criteria. Based on incremental approach, the similarity solution is derived for circular and spherical openings using generalized H-B and M-C yield criteria. The classical Runge-Kutta method is used to integrate the first-order ordinary differential equations. Using three data sets for M-C and H-B models, the results of the radial displacements, the spreading of the plastic radius with decreasing pressure, and the radial and circumferential stresses in the plastic region are compared. Excellent agreement among the solutions is obtained for all cases of spherical and circular openings. The importance of the use of proper initial values in the similarity solution is discussed.

• 한국공간구조학회논문집
• /
• 제20권2호
• /
• pp.59-68
• /
• 2020

In this paper, the dynamic snapping of the 3-free-nodes spatial truss model was studied. A governing equation was derived considering geometric nonlinearity, and a model with various conditions was analyzed using the fourth order Runge-Kutta method. The dynamic buckling phenomenon was observed in consideration of sensitive changes to the force mode and the initial condition. In addition, the critical load level was analyzed. According to the results of the study, the level of critical buckling load elevated when the shape parameter was high. Parallelly, the same result was caused by the damping term. The sensitive asymmetrical changes showed complex orbits in the phase space, and the critical load level was also becoming lowly. In addition, as the value of damping constant was high, the level of critical load also increases. In particular, the larger the damping constant, the faster it converges to the equilibrium point, and the occurrence of snapping was suppressed.

• 한국지구과학회지
• /
• 제39권4호
• /
• pp.317-326
• /
• 2018

Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.

• 한국가시화정보학회지
• /
• 제12권1호
• /
• pp.13-17
• /
• 2014

The growth of micro bubble has been simulated under the variation of ambient pressure. The Rayleigh-Plesset equation governs the dynamic growth and collapse of a bubble according to pressure and temperature conditions. The Rayleigh-Plesset equation was solved by 4th-order Runge-Kutta method for wide range of pressure variations. As numerical parameters, the pressure difference between initial and final pressures, and the temporal pressure gradient are changed. The results show that the pressure difference has little effect on the growth rate of the micro bubble in the inertia controlled growth region. On the other hand, the growth rate increases linearly with the increase of the pressure gradient.

• 한국소음진동공학회논문집
• /
• 제15권1호
• /
• pp.29-38
• /
• 2005

The track pull-in behavior analysis in an optical disk drive (ODD) using plane phase and the evaluations for effecting parameters of it are discussed. Track pull-in, track capture procedure to do track following control, is a key factor to increase data transfer rate. First, the relative velocity between the beam spot of an optical pick-up and the target track of an optical disk is analyzed during the track pull-in procedure. In this process, it is showed that the track error signal has nonlinear characteristics which are depending on the time. Second, Runge-Kutta method to solve the nonlinear equation is applied to find the track pull-in behavior, and some optimal parameters to get stable and fast pull-in condition are obtained. Then, the phase plane analysis for track pull-in procedure is presented. Finally, some comments for the simulated results are discussed briefly.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 추계학술대회논문집
• /
• pp.1266-1271
• /
• 2007

General coil spring has linearity. However, disc spring has non-linearity so that using this non-linearity disc spring can be designed to do shock-absorbing in cases we need because shock response also has non-linearity. By changing the shape and stacking number, it is satisfactory with response of displacement, velocity and acceleration. Conventionally, disc spring was used to control the vibration against huge load and limited space. However, it is limitedly used because of difficulty of the designing guidance. Therefore, disc spring is needed to study further in order to apply it widely. Response of disc spring is compared to response of coil spring by changing ho/t radio with computer simulation and the usage of disc spring is increased through analysis of effect of design factors. The purpose of this paper is that the shock response of disc spring is calculated through numerical simulation and effect of ho/t and stiffness is analyzed to broad usage so that design factor of disc spring is presented.

• 대한기계학회논문집
• /
• 제17권8호
• /
• pp.2122-2130
• /
• 1993

A numerical and experimental study has been performed on the flow in a shallow rectangular tank accompanying a vortex shedding. The model is composed of a rectangular tank with a vertical plate with a length half the width of the tank. The tank is subject to a horizontal sinusoidal oscillation. The numerical analysis shows that the pattern of vortex shedding changes considerably when the Reynolds number $R_e$ is varied from 500 to 7500. It is symmetric for $R_e$ <1500 and asymmetric for $R_e$ > 1500. The kinetic energies of the right-hand and left-hand sides of the vertical plate are used to quantify the degree of the asymmetry. Experimental visualization is carried out at $R_e$ = 3876 and 52000. The development of the streamline pattern at $R_e$ = 3876 is in closer agreement with the numerical result at $R_e$ = 1000 than that at $R_e$ =3876. The asymmetric pattern is observed at $R_e$ = 52000.

• Journal of Mechanical Science and Technology
• /
• 제14권6호
• /
• pp.690-698
• /
• 2000

The laminar boundary layer flow and heat transfer of anisotropic fluids in the vicinity of a wedge have been examined with constant surface temperature. The similarity variables found by Falkner and Skan are employed to reduce the stream wise-dependence in the coupled nonlinear boundary layer equations. The numerical solutions are presented using the fourth-order Runge - Kutta method and the distribution of velocity, micro-rotation, shear and couple stresses and temperature across the boundary layer are plotted. These results are also compared with the corresponding flow problems for Newtonian fluid over wedges. It is found that for a constant wedge angle, the skin friction coefficient is lower for micropolar fluid, as compared to Newtonian fluid. For the case of the constant material parameter K, however, the magnitude of velocity for anisotropic fluid is greater than that of Newtonian fluid. The numerical results also show that for a constant wedge angle with a given Prandtl number, Pr = I, the effect of increasing values of K results in increasing thermal boundary layer thickness for anisotropic fluid, as compared with Newtonian fluid. For the case of the constant material parameter K, however, the heat transfer rate for anisotropic fluid is lower than that of Newtonian fluid.

• 수산해양기술연구
• /
• 제42권3호
• /
• pp.127-133
• /
• 2006

An estimation of the headline height of a bottom trammel net set across under uniform current was achieved numerically from a differential equations describing the forces of the net and compared with the measured value in a flume tank experiment. The analysis on the shape of the bottom trammel net with the headline free was based on the equilibrium equation of the bottom gill net which was modified and slack of the trammel net was varied with net depth as shown in the tank experiment. The differential equations were solved by a forth-order Runge-Kutta method. The estimated headline heights with varied slack was found to be closer than that with constant slack when compared with the actual values.