• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.354-359
• /
• 2013

일반적으로 건축물의 설계시 풍동 실험을 통한 풍환경의 평가를 수행하고 있으며, 이는 환경 영향 평가법에서 정한 건축 사업 시행 시 수반되어야 할 자연환경, 생활환경 그리고 사회경제환경의 영향 평가의 일환으로 실시되고 있다. 그러나, 풍동 실험의 경우 여러 가지 현실적 제약조건으로 설계와 실험의 피드백 (Feedback)이 원활하지 못하며, 특히 대상 건축물이 공장과 같이 대기 오염원이 되는 경우 실험은 더욱 어려운 형편이다. 이에 대한 보완책으로 전산 유체 역학을 이용한 건축물의 풍압 해석에 의한 풍하중 추정이나 인접 지형-지물의 영향을 고려한 건축물 주위의 풍환경 평가가 있다. 전산 모사에 의해 풍동 실험의 미비점을 보완하고, 보다 상세한 정보를 확보함으로써 건축물의 구조적 안전성의 증대와 환경 피해 감소를 기할 수 있다. 그러나 복잡한 지형-지물이나 건축물 주위의 풍환경에 대한 전산 모사는 주로 두 가지의 기술적 어려움을 수반하게 되다. 그 중 하나는 고정 경계면을 이루는 형상의 복잡성으로 인해 기존에 많이 이용하고 있는 Body-fitted 격자계를 이용하는 경우, 격자 생성 과정이 매우 복잡하고 어려울 뿐 만 아니라 생성된 격자가 주로 비정렬 (unstructured) 특성을 갖게 되어 수치해석 과정의 효율을 저하시키는 요인이 되며, 격자의 형상도 수치해석의 수렴성을 저하시키는 예가 많다. 다른 어려움으로 풍환경은 전형적인 난류 유동장으로서 난류의 전산 해석은 아직도 해결하지 못한 부분이 많다는 점이다. 이에 본 논문에서는 복잡한 지형-지물이나 건축물의 풍하중과 풍환경의 전산 모사 기술 확보를 위하여 수행중인 연구의 일환으로 물체 형상의 기하학적 복잡성의 극복을 위한 가상경계법 (Immersed Boundary Method)과 난류 유동장의 물리적 엄밀성을 높이기 위한 다와동 모사 (Large Eddy Simulation)을 이용한 물체 형상과 무관한 유동장 해석 기술 개발에 대하여 다루고자 한다. 먼저 최근에 유동 해석에 이용되는 방법인 가상경계법(IBM)은 물체를 포함한 전체 전산 영역을 직교 좌표계에 의해 이산화하고, 유동장내 존재하는 물체의 표면에서의 점착 조건을 만족시키기 위하여 지배 방정식에 적절한 외력을 추가로 고려하는 방법이다. 본 연구에서는 가상경계법을 이용하여 경계층에 위치한 건물 형상의 각진 물체 주위 사이에 형성되는 공동 내부의 비정상 유속 및 압력에 대한 전산 해석을 수행하고, 풍상측 전면에 형성되는 경계층에 의한 영향을 분석하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.655-659
• /
• 2008

본 연구에서는 3차원 수치모의를 통하여 횡방향 하상형상 및 격자형 이차흐름 구조의 생성 메커니즘을 분석하였다. 이를 위해 곡선좌표계에서의 지배방정식을 구성하고, 난류 폐합을 위해 Speziale(1987)가 제안한 비등방성 k-$\varepsilon$모형을 이용하였다. 또한 Exner 방정식을 이용하여 시간에 따른 하상변동을 예측하였다. 그 결과 바닥과 측벽 사이에서 발생되는 바닥 이차흐름의 하향류에 의해 측벽 부근부터 하상이 침식되고, 침식된 유사량은 이차흐름의 횡방향 유속에 의해 이동되어 퇴적되어, 결국 횡방향으로 연속적인 언덕 및 저면과 같인 하상형상이 생성되는 것으로 나타났다. 또한 시간에 따른 이차흐름 및 바닥 전단력, 하상고의 변화에 대해 살펴보았다.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.19-30
• /
• 2000

A finite-volume method is developed to solve turbulent flows around modern commercial hull forms with bow and stern bulbs. The RANS equations are solved. The cell-centered finite-volume method employs QUICK and central difference scheme for convective and diffusive flux discretization, respectively. The SIMPLEC method is adopted for the velocity-pressure coupling. The developed numerical methods are applied to calculate turbulent flow around KRISO 3600TEU container ship. Surface meshes are generated into five blocks: bow and stern bulbs, overhang, fore and afterbody. 3-D field grid system with O-H topology is generated using elliptic grid generation method. Surface friction lines and wake distribution at propeller plane is compared with experiment. The calculated results show that the present method can be used to predict flow around a modern commercial hull forms with bulbs.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1999.04a
• /
• pp.11-11
• /
• 1999

액체 램제트 연소기는 흡입공기와 분무, 혼합 그리고 이에 따른 연소 등 일련의 과정에 따라 다수의 복잡한 현상들이 상호 밀접하게 관련되어 있다. 본 연구에서는 액체 램제트 연소기내의 유동특성을 파악하기 위해서 2차원 및 3차원 연소기 형상에 대해서 수치적 실험을 수행하였으며, 격자구성은 연소기에 공기를 공급하고 연료를 분무하는 공기 유입관 영역과 연소실 영역, 그리고 출구 대기 영역으로 나누어 독자적으로 격자를 생성시켰다. 2차원과 3차원 유동해석을 비교하였고 분무모델의 적용에 따른 연소특성 및 분사위치에 따른 연소특성을 비교하였다. 유동해석 결과 2차원과 3차원의 유동특성은 달랐으며, 분무모델을 적용해야 정확한 연소 유동 현상을 예측할 수 있음을 알 수 있었다. 그리고 유입관의 안쪽에 연료의 분사위치를 준 경우가 연소의 안정화에 필요한 재순환영역으로의 연료의 혼합이 잘 되어 유입관 바깥쪽에 연료를 분사시키는 것보다 좋은 분사위치임을 알 수 있었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.222-222
• /
• 2023

오픈 소스 유체역학 소프트웨어인 OpenFOAM은 다양한 유체 흐름에 적용 가능한 프로그램들로 구성되어 있다. 이 중 interFoam은 밀도가 다른 두 유체(i.e., 물, 공기) 간의 경계를 추적하는 기법을 기반으로 한 프로그램으로, 파랑의 거동 모의에 주로 쓰이고 있다. 파생형 프로그램으로는 동격자(dynamic mesh) 및 중첩 격자 기법(overset grid method)을 interFoam에 추가한 overInterDyMFoam이 있다. 두 기법을 사용하면 각각 여러 영역에서 유체흐름과 다중 물체 간의 상호작용을 효율적으로 모의할 수 있다. 본 연구에서는 overInterDyMFoam을 사용하여 두 개 수문의 개방 움직임을 구현하고 생성된 파랑이 포말대(swash zone)에 접근하였을 때의 흐름 특성을 조사하였다. 수치모형실험 결과 수문 개방 속도가 댐 붕괴 파랑 흐름 전파속도에 영향을 미치는 사실을 발견하였다. 또한, 처오름과 처내림의 상호작용에 의한 난류 운동 특성을 조사하기 위해 수문 개방시간 간격을 0초~3초로 설정하였다. 수치모형실험 결과는 수리모형실험의 수면 변동 시계열과 속도 시계열 결과와 비교하여 모형의 정확성이 검증되었다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.8
• /
• pp.1007-1013
• /
• 2010

A new local refinement scheme for spline finite element method has been proposed; this scheme involves the use of hierarchical B-spline. NURBS has been widely used in CAD; however, the local refinement of NURBS is difficult due to its tensor-product property. In this study, we attempted to use hierarchical B-splines as local refinement strategy in spline FEM. The regions of high gradients are overlapped by hierarchically-created local meshes. Knot vectors and control points in local meshes are extracted from global meshes, and they are refined using specific schemes. Proper compatibility conditions are imposed between global and local meshes. The effectiveness of the proposed method is verified on the basis of numerical results. Further, it is shown that by using a proposed local refinement scheme, the accuracy of the solution can be improved and it could be higher than that of the solution of a conventional spline FEM with relatively lower degrees of freedom.

• Journal of Korean Society for Geospatial Information Science
• /
• v.1 no.2 s.2
• /
• pp.83-96
• /
• 1993

The orthophto map is seen as the form of picture with the uniform reduced scale as the current terrain map. Thus it provides a reasonable feeling of scene and is easy to be interpreted. Furthermore, digital orthophoto is currently used as the basic terrain information data of the Geo-Spatial Information System(GSIS). Therefore, the orhtophoto map has high potential use as a future terrain map. This paper studies the method of producing orthophoto map by using the digital satellite imagery data taken from SPOT satellite of France. The production of orthophoto map requires the process of generating orthophoto imagery with digital elevation model, which process is called digital differential rectification. As the final accuracy of orthophoto map depends on that of digital elevation model, the precise and efficient production method of digital elevation model should be preceded. This study investigated the method of producing digital elevation model directly from SPOT satellite imagery and generated ortho-image by resampling the original SPOT imagery through digital differential rectification. Finally, Simple orthophoto map was made by overlaying the ortho-image and the contour map from digital elevation model.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1634-1638
• /
• 2008

하천의 2차원 흐름 해석, 유사이동 해석, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있다. 유체에 대한 수치해석 기법으로 전통적으로 가장 많이 사용되고 있는 방법은 유한차분법이지만, 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 다양한 형태의 하천 해석에는 더욱 적합할 것이다. 본 연구에서는 비구조적 요소망을 advanced front method를 이용하여 생성해 보았다. Advanced front method는 해석하고자 하는 영역에 적절한 절점들을 생성한 후 삼각 요소망을 구성하는 grid based advanced front method와 절점들을 생성하지 않고 해석 영역에 삼각 요소를 바로 구성하는 element based advanced front method로 나눌 수 있다. Grid based advanced front method에서 해석 영역에 적절한 절점을 생성하는 방법으로는 일반적인 격자 구조의 절점 생성 방법을 적용하였으며 경계와의 거리가 가까운 절점은 생성되지 않으며, 삼각 요소를 구성할 때에는 두 개의 인접 절점을 비교하여 최적의 삼각 요소를 구성하게 된다. 단 두 개의 인접 절점만을 비교함으로서 비교적 빠른 시간 안에 최적의 삼각 요소망을 구성할 수 있다. 삼각 요소망을 생성한 후에는 Laplacian smoothing을 이용하여 삼각 요소망의 형질을 개선하였다. Element based advanced front method는 외부 경계에서부터 시작된 Front가 내부 영역으로 확대되어지며 각 Front에서 적절한 절점을 직접 생성하여 바로 삼각 요소를 구성하게 된다. Front에서 생성된 절점은 인접 절점들이 있는지 검색하여 인접 절점이 있다면 생성된 절점은 삭제되어지며 인접 절점이 삼각 요소를 위한 나머지 한 점으로 채택되어진다. Front는 외부 경계와 교차되어지지 않아야 하며 또한 연속된 Front를 효율적으로 관리하기 위해 list 자료 구조를 활용하였다.

• Korean Journal of Remote Sensing
• /
• v.20 no.2
• /
• pp.91-102
• /
• 2004

This paper reflects the estimation of using the EOC(Electro-optical Camera) images supporting GIUH(geomorphological instantaneous unit hydrograph) approach. We have analyzed GIUH in its density and frequency distribution by creating a DEM(digital elevation model) for the sub basin produced from the EOC images and examined topographical and hydrological application possibility of the EOC images. In this process, we have topographical basin characteristic analysis that use the remote sensing technique analyzing the DEM creation process of the EOC stereo images by studying the basic topographical hydrology analysis about abstraction technique since it is flirty complex and is more time-consuming than other method. we executed statistical analysis of a basin size and river length using the frequency function after divided lattice spacing applied have to the sub river basin from the image data and the digital map into 10m intervals ranging from 10m to 100m. After comparing and examining the peak and time to peak of the GIUH, we proceeded with a comparative analysis by lattice concerning the topographical divergence rate, area ratio, length ratio. Accumulating the peak and time to peak of the GIUH is altered to non-linear form in accordance to lattice dimension as well as basin factor. It was proved that the lattice dimension is one of the important factors about the peak and time to peak of the GIUH.