• Proceedings of the Korea Water Resources Association Conference
• /
• 1989.07a
• /
• pp.113-114
• /
• 1989

원자력, 화력발전소, 제철소 등의 다양한 임해 공업시설로부터 방출되는 냉각용 온배수는 연안일대 해수의 열균형을 파괴시켜 생태계의 보존 및 환경관리에 좋지 못한 영향을 야기케 된다. 이러한 영향은 해안 또는 만내의 수중온도를 전반적으로 높일뿐 아니라 가동중단시 갑자기 수온을 떨어뜨려 해양생물상에 피해를 줄수도 있다. 또한 온배수에 의하여 온도가 상승된 해수가 취수구를 통하여 재순환되어 냉각기능의 부진을 초래하게되면 발전효율 또는 기계가동율을 저하시키게 된다. 이러한 측면에서 온배수의 확산에 대한 정성, 정량적인 예측은 환경영향평가, 취.배수구 설계조건의 산정 등에 매우 중요한 문제라 하겠다. 본 연구는 정지수역으로 유입하는 3차원 정상류 표면온배수 해석모형의 개발로서 개발된 모형의 수치실험을 통하여 온배수 확산의 물리적 특성을 규명한다. 지배방정식에 나타나는 Reynolds 응력항($)과 온도유동 프럭스항($)의 해석에서 필요한 난류모델은 k-$\varepsilon$ 모형에 난류 평균자승 온도유동($) 및 그 감쇄방정식을 추가한 4-방정식 모델로서 구성하였다. 아울러 3차원 정상류 모형에서 야기되는 타원형 방정식을 포물형 방정식의 형태로 전환하여 효과적으로 해석할 수 있도록 모델의 특성을 정리하였다. 본 모델의 검증을 위하여 Lal 및 Rajaratnam(1977)의 물리적 실험값과 비교해본 결과 온배수 거동의 물리적현상이 잘 일치하였다. 또한 McGuirk 및 Rodi(1979)에 의해 개발된 2-방정식 k-$\varepsilon$ 난류모형의 해석결과에 대하여 비교분석을 실시하였다.

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

동해를 전파하는 지진해일은 다른 지역에서 발생하는 지진해일과 비교하였을 때 상대적으로 파장이 짧고, 파장에 비해 먼거리를 전파한다. 따라서 지진해일이 전파할 시 물리적인 분산효과가 매우 중요하다. 그러므로 동해에서 발생하는 지진해일을 수치모의 할 때는 분산효과가 충분히 고려될 수 있는 선형 Boussinesq 방정식을 사용한다. 그러나 이를 직접 풀 경우에는 상당히 많은 시간이 소비되며 효율적이지 못하다. 이와 같은 단점을 극복하기 위해 기존의 연구에서는 leap-frog 기법을 사용하여 선형 천수방정식을 차분할 때 발생하는 수치분산항을 선형 Boussinesq 방정식의 분산항과 같은 형태를 가질 수 있도록 분산보정계수를 사용하여 수치모의를 수행하였다. 하지만 이때 사용된 지배방정식은 수심이 일정하다는 가정을 이용하여 유도된 것이므로, 실제 경사가 있는 지형을 통과할 때의 수치모의 결과는 정확하다고 할 수 없다. 본 연구에서는 이를 극복하기 위하여 바닥 지형이 1차원으로 변한다는 가정으로 새로운 지배방정식을 유도하였으며, 수심변화로 인해 새로 발생하는 항을 기존의 분산보정기법에 추가하였다. 또한 수심이 변화는 지형을 통과하는 지진해일의 분산효과가 충분히 고려되는지 확인하기 위하여 Gaussian hump를 이용하여 가상 지진해일을 원형 천퇴지형에 통과시키는 수치모의를 수행하였다. 결과의 비교를 위한 정확해가 없으므로, 비선형 Boussinesq 방정식을 직접 차분하여 푸는 FUNWAVE를 이용하여 동일한 조건으로 수치모의를 수행하였다. 수치모의 시 중심선에 4개의 가상 gage를 설치하였으며, 이를 통해 각각의 수치모의 실험에 대한 자유수면 변위를 관찰하여 비교하였다. 수치모의 결과에 대한 비교를 통하여 기존의 분산보정기법에 비해 본 연구에서 제안한 새로운 수치기법이 분산효과를 비교적 잘 반영하는 것으로 나타났으며, 비교적으로 실제 지형에 적용하였을 때 정확도 향상의 가능성이 높다고 판단하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.4
• /
• pp.181-189
• /
• 2000

A weakly nonlinear mild-slope equation has been derived directly from the continuity equation with the aid of the Galerkin's method. The equation is combined with the momentum equations defined at the mean water level. A single component model has also been obtained in terms of the surface displacement. The linearized form is completely identical with the time-dependent mild-slope equation proposed by Smith and Sprinks(1975). For the verification purposes of the present nonlinear model, the degenerate forms were compared with Airy(1845)'s non-dispersive nonlinear wave equation, classical Boussinesq equation, andsecond¬order permanent Stokes waves. In this study, the present nonlinear wave equations are discretized by the approximate factorization techniques so that a tridiagonal matrix solver is used for each direction. Through the comparison with physical experiments, nonlinear wave model capacity was examined and the overall agreement was obtained.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.4
• /
• pp.355-364
• /
• 2011

In this study, the details of WSA(wavelet series analysis) have been demonstrated to solve the 4th-order elliptic differential equation. It is clear to solve the 2nd-order elliptic differential equation with the basis function of Hat wavelet series that is used in the previous study existed in $H^1$-space. However, it is difficult to solve the 4th order differential equation with same basis function of Hat wavelet series because of insufficient differentiability and integrability. To overcome this problem, the linear equations in terms of moment and deflection have been formulated and solved sequentially that are similar to extension of Elastic Load Method and Moment Area Method in some senses. Also, the differences and common points between the proposed method and the meshless method are discussed in the procedure of WSA formulation. As we expect, it is easy to ascertain that the more terms of Hat wavelet series are used, the better numerical solutions are improved. Also the solutions obtained by WSA have been compared with the conventional FEM solutions in case of Euler beam problems with stress singularity.

• Fire Protection Technology
• /
• s.45
• /
• pp.4-4
• /
• 2008

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.4
• /
• pp.5-11
• /
• 2018

In conventional analytical solutions for rainfall-induced soil slope stability, the Green-Ampt (1911) equation for estimating the saturation depth and the Skempton & DeLory (1957) equation for calculating the infinite slope shallow failure were compared with the numerical analysis to confirm the error. In the simple evaluation of the reason of soil slope instability due to rainfall using the conventional equations, there are many errors and, overestimation or underestimation of the calculation results. In this study, the equation consisting of the results obtained from infiltration analysis on unsaturated soil slope is proposed by applying the average range of the strength parameters of the granite weathered soils, and its reliability is verified by comparing with the numerical analysis results. The developed equation can be used easily in various fields for the estimation of slope safety factor by checking the rainfall duration and saturation depth.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.4
• /
• pp.191-197
• /
• 2018

This paper details the development of an isotropic elastic-damage constitutive model for polymeric foam based on irreversible thermodynamics to consider the growth and coalescence of voids. The constitutive equations describe the material behavior sustaining unilateral damage. To facilitate finite element analysis, the material properties for specific types of polymeric foams are applied to the developed model; the model is then implemented in ABAQUS as a user-defined material subroutine. To validate the developed damage model, the simulated results are compared to the results of a series of tensile tests on various polymeric foams. The proposed damage model can be utilized to further research on continuum damage mechanics and finite element analysis of polymeric foams in computational engineering.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.3D
• /
• pp.425-433
• /
• 2006

The presented work provided a predictive equation for dynamic modulus of hot mix asphalt, which showed higher reliability and more simplicity. Lots of test result by UTM at laboratory has been used to develop the precise predictive equation. Evaluation of dynamic modulus for 13mm and 19mm surface course and 25mm of base course of hot mix asphalt with granite aggregate and two asphalt binders (AP-3 and AP-5) were carried out. Superpave Level 1 Mix Design with gyrator compactor was adopted to determine the optimum asphalt binder content (OAC) and the measured ranges of OAC were between 5.1% and 5.4% for surface HMA, and around 4.2% for base HMA. The dynamic modulus and phase angle were determined by testing on UTM, with 5 different testing temperature (-10, 5, 20, 40, & $55^{\circ}C$) and 5 different loading frequencies (0.05, 0.1, 1, 10, 25 Hz). Using the measured dynamic modulus and phase angle, the input parameters of Sigmoidal function equation to represent the master curve were determined and these will be adopted in FEM analysis for asphalt pavements. The effect of each parameter for equation has been compared. Due to the limitation of laboratory tests, the reliability of predictive equation for dynamic modulus is around 80%.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.2
• /
• pp.85-93
• /
• 2024

In this paper, we propose a method for establishing a state-space equation model for the motion analysis of floating structures subjected to wave loads, by applying system-identification techniques. Traditionally, the motion of floating structures has been analyzed in the time domain by integrating the Cummins equation over time, which utilizes a convolution integral term to account for the effects of the retardation function. State-space equation models have been studied as a way to efficiently solve floating-motion equations in the time domain. The proposed approach outlines a procedure to derive the target transfer function for the load-displacement input/output relationship in the frequency domain and subsequently determine the state-space equation that closely approximates it. To obtain the state-space equation, the method employs the N4SID system-identification method and an optimization approach that treats the coefficients of the numerator and denominator polynomials as design variables. To illustrate the effectiveness of the proposed method, we applied it to the analysis of a single-degree-of-freedom model and the motion of a six-degree-of-freedom barge. Our findings demonstrate that the presented state-space equation model aligns well with the existing analysis results in both the frequency and time domains. Notably, the method ensures computational accuracy in the time-domain analysis while significantly reducing the calculation time.