• The Korean Journal of Rheology
• /
• v.3 no.1
• /
• pp.47-55
• /
• 1991

본 연구에서는 구형기포가 미분형의 upper convected Maxwell 모델을 따르는 유체 내에서 수축할 때의 현상을 이론적으로 해석하였다. 수치해법으로는 Lagrangian 좌표계에서 지배방정식을 유도 사용함으로써 자연스럽게 자유표면을 추적하는 동시에 압력변수도 반복 에 의하지 않고 직접적인 방법으로 계산할 수 있는 Galerkin-유한요소법을 개발사용하였다. 본연구의 결과 유체의 탄성은 변형초기에 충분히 발달치 못하기 때문에 수축을 가속화시키 지만 수축 후기에는 지연시킴을 알 수 있었다. 수축의 속도는 적분형 Maxwell 유체에서 보 다 빠른 것을 알수 있었는데 이는 유체의 정지이력에 의한 것으로 판단되었다, 또한 Maxwell 유체내에서 기포가 수축할 경우 탄성에 의한 반동현상이 나타나며 반동이 점성에 의하여 감쇄될 때의 진폭과 주기는 비례함을 보였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.3
• /
• pp.165-173
• /
• 1987

A numerical model is presented that predicts directly the wave angle and height at every point on the grid. The governing equations used are conservation of waves equation and conservation of energy equation which are derived from the basic linear potential equations by means of an asympotic approximation. Finite difference methods are used to solve the governing equations and the solution is obtained for a finite number of rectilinear grid cells that comprise the domain of interest. Model results are compared with the results obtained from wave ray methods and it shows no significant differences between two results. The model is especially efficient for modeling large areas of coastline with arbitrary bathymetry, and therefore it is anticipated to be used in many coastal engineering problems such as littoral drift problems.

• Water for future
• /
• v.18 no.1
• /
• pp.85-94
• /
• 1985

A method is developed for the separation of the major base flow in a river hydrograph combining the numerical techniques and the empirical methods. The linearized Boussinesq equation and the storage function are used to obtain the base flow recession. The shape of base flow curve made by the recharge of the groundwater table aquifer resulting from rainfall in determined by the Singh and Stall's graphical method, and the continuous from for the curve is approximated by the multiple and polynomial regression. this procedure was successfully tested for the separation of base flow and the establishment of hydrograph in a natural watershed. It was found that the direct numerical method applied to the homogeneous linear second order ordinary differential equation system is not suited to obtain the recession curve, and the case that the loss is generated in the partially penetrating stream can not be solved by the method of this study.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.31 no.2
• /
• pp.202-209
• /
• 1998

A three-dimensional numerical method based on the Green's integral equation is developed to predict the motion response and drift force in waves for the ocean monitoring facilities. In this method, we use source and doublet distribution, and triangular and rectangular eliments. To eliminate the irregular frequency phenomenon, the method of improved integral equation is applied and the time-mean drift force is calculated by the method of direct pressure integration over the body surface. To conform the validity of the present numerical method, some calculations for the floating sphere are performed and it is shown that the present method provides sufficiently reliable results. As a calculation example for the real facilities, the motion response and the drift force of the vertical cylinder type ocean monitoring buoy with 2.6 m diameter and 3,77 m draft are calculated and discussed. The obtained results of motion response can be used to determine the shape and dimension of the buoy to reduce the motion response, and other data such as the effect of motion reduction due to a damper can be predictable through these motion calculations. Also, the calculation results of drift force can be used in the design procedure of mooring system to predict the maximum wave load acting on the mooring system. The present method has, in principle, no restriction in the application to the arbitrary shape facilities. So, this method can be a robust tool for the design, installation, and operation of various kinds of the floating-type ocean monitoring facilities.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.35 no.6
• /
• pp.1237-1246
• /
• 2015

A novel technique was proposed to calculate fluxes accurately by separation of flow area into a part of step face which is dominated by flow resistance of it and an upper part which is relatively less affected by the step face in analyzing shallow-water flows over discontinuous topography. This technique gives fairly good agreement with exact solutions, 3D simulations, and experimental results. It has been possible to directly analyze shallow-water flows over discontinuous topography by the technique developed in this study. It is expected to apply the developed technique to accurate evaluation of overflows over weirs or retaining walls (riverside roads) and areas flooded by the inundation in the city covered in discontinuous topography.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1990.11a
• /
• pp.69-72
• /
• 1990

최근의 전자및 광학기기 분야에 있어서의 준부신 발전은 다면경 가공기나 초정밀 절삭, 연삭기와 같은 초정밀 가공기계의 개발과 실용화에 힘입은 바 크다. 이러한 초정밀 가공기의 성능을 좌우하는 핵심 요소로서 주축계를 들 수 있으며, 비교적 소형 경량의 공작물을 가공하는 기계의 주축용 베어링으로는 볼 베어링이나 오일 베어링을 대신하여 공기베어링이 점차 널리 사용되고 있다. 일반적으로 주축으로 사용되는 베어링은 원통형 레이디얼 베어링과 원판형 스러스트 베어링이 결합된 형식이 주류를 이루나 이러한 베어링은 스러스트 판과 축의 직가곧 가공오차가 존재하기 때문에 가공하기는 쉬우나 회전시에 의의 영향에 의해 회전 정밀도 유지가 어렵다는 단점을 지니고 있다. 이러한 단점을 보완하기 위하여 사용되는 베어링에는 원추형(conical) 베어링과 구면형(spherical) 베어링이 쓰이고 있다. 이러한 원추형 베어링과 구면형 베어링은 가공오차를 베어링과 축의 현압 연마로써 없애줄 수 있으며 베어링이 축방향 하중과 경방향 하중을 동시에 지지하여 줌으로써 기계 전체의 부피를 줄이고 회전 정밀도를 향상시켜 주는 것으로 알려져 있다. 그러나 구면의 베어링 간극을 정확히 가공하기 어려운 단점이 있어 축과 베어링을 현압연마하여 가공한 후에 두부품을 중심선상에서 분리시키므로써 요구되는 간극을 얻을 수 있는 원추형 베어링이 많이 쓰이고 있다. 본 연구에서는 직접 수치 해법을 이요하여 원추형 베어링의 유막내의 압력 분포를 계산하고 이 합력인 하중지지 용량이 축방향 하중과 경방향 하중을 지지하는 특성을 이론적으로 검토하여 외부 가압 원추형 베어링으 특성수를 파악하여 설계자료를 제시하고자 한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.1088-1092
• /
• 2005

하도에서의 홍수추적시 유역내의 지류 혹은 지표 및 지하수등은 추적대상이 되는 하도구간내에서 측방유입량의 되어 유출수문곡선의 첨두유량, 첨두시간, 수문곡선의 형태등에 영향을 주므로 정확한 산정이 필요하며, 직접유출수문곡선에서 측방유입량은 지표유출에 의해 발생하므로 강우발생시 유역에서 하도까지 걸리는 도달시간의 산정이나 측방유입속도의 결정이 필요하다. 기존의 강우-유출 수문모형은 지표수흐름의 복잡한 메카니즘 및 수리특성을 규명하는데 어려움이 있다. 본 연구에서는 관측유입수문곡선 및 유출수문곡선을 이용하여 측방유입량을 산정하고, 하도구간으로 유입되는 기지의 측방유입량으로부터 수리학적 홍수추적을 위한 지배방정식인 Saint-Venant방정식의 수치해법중 하나인 양해법에 diffusing scheme을 적용하였다. 또한 하도 전구간에 동일한 측방유입속도로 유입될 경우와 하도중심을 기점으로 상류부와 하류부로 구분하여 두 구간의 측방유입속도가 다른 두가지 경우에 대해 측방유입속도를 역추정하였으며, 계산 유출 수문곡선과 관측 유출수문곡선을 비교$\cdot$분석함으로써 구성한 홍수추적모형에 대한 정확성과 타당성, 적용 가능성등을 검증하고자 하였다.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.19-27
• /
• 2016

The flow analyses through a porous hydraulic fractures is among the most important tasks in recently developed shale reservoirs but is rendered difficult by non-Darcy effects and geometric changes in the hydraulic fractures during production. In this study, several Computational Fluid Dynamics(CFD) models of hydraulic fractures, with a simple shape such as that of parallel plates, filled with proppants were built. Direct Numerical Simulation(DNS) analyses were then carried out to examine the flow loss characteristics of the fractures. The hydraulic diameters for the simulation models were calculated using the DNS results, and then they were compared with the results from Kozeny's definition of hydraulic diameter which is most widely used in the flow analysis field. Also, the characteristic parameters based on both hydraulic diameters were estimated for the investigation of the flow loss variation features. Consequently, it was checked in this study that the hydraulic diameter based on Kozeny's definition is not accordant to the results from the DNS analyses, and the case using the CFD results exhibits f Re robustness like general pipe flows, whereas the other case using Kozeny's definition doesn't. Ultimately, it is expected that discoveries reported in this study would help further porous flow analyses such as hydraulic fracture flows.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.3
• /
• pp.91-105
• /
• 1992

For the analysis of multistory frames with sidesway, no adequate procedures can be found in the classical methods of structural analysis. Even well-known procedures such as the slope-deflection method and the moment distribution method may not be effective tools since those methods require a multiple of computational labor and/or yield results of approximate values. In this study, a direct method is developed and proposed for the analysis of multistory frames with sidesway, which is due to the lateral loads, asymmetry of the structure itself, or asymmetry of vertical loadings. The proposed method is to obtain simple forms of equations derived by a mathematical formulation of the moment distribution procedure combined with successive correction concept. Numerical illustrations show that the results obtained by the proposed method agree well with those by rigorous ones. Undoubtedly, this newly developed method can be applied more easily for the analysis of structural frames without joint translation as well as continuous beams.

• Geophysics and Geophysical Exploration
• /
• v.26 no.3
• /
• pp.114-125
• /
• 2023

The physics-informed neural network (PINN) has been proposed to overcome the limitations of various numerical methods used to solve partial differential equations (PDEs) and the drawbacks of purely data-driven machine learning. The PINN directly applies PDEs to the construction of the loss function, introducing physical constraints to machine learning training. This technique can also be applied to wave equation modeling. However, to solve the wave equation using the PINN, second-order differentiations with respect to input data must be performed during neural network training, and the resulting wavefields contain complex dynamical phenomena, requiring careful strategies. This tutorial elucidates the fundamental concepts of the PINN and discusses considerations for wave equation modeling using the PINN approach. These considerations include spatial coordinate normalization, the selection of activation functions, and strategies for incorporating physics loss. Our experimental results demonstrated that normalizing the spatial coordinates of the training data leads to a more accurate reflection of initial conditions in neural network training for wave equation modeling. Furthermore, the characteristics of various functions were compared to select an appropriate activation function for wavefield prediction using neural networks. These comparisons focused on their differentiation with respect to input data and their convergence properties. Finally, the results of two scenarios for incorporating physics loss into the loss function during neural network training were compared. Through numerical experiments, a curriculum-based learning strategy, applying physics loss after the initial training steps, was more effective than utilizing physics loss from the early training steps. In addition, the effectiveness of the PINN technique was confirmed by comparing these results with those of training without any use of physics loss.