• Title/Summary/Keyword: Mesh Refinement

Search Result 189, Processing Time 0.025 seconds

Development of 2D Urban Inundation Analysis Model using Adaptive Mesh Refinement Method (메쉬 세분화 기법을 이용한 2차원 침수해석 모형의 개발)

  • Lee, Seung-soo
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2016.05a
    • /
    • pp.93-93
    • /
    • 2016
  • 최근 증가하고 있는 기후변화에 의해 설계빈도를 상회하는 강우의 발생빈도가 증가하고 있으며, 이로 인한 도시유역의 내수범람 피해가 증가하고 있다. 도시유역에서 발생하는 침수 피해의 경우 인적 물적 자원이 집중되어 있는 도시의 특성으로 인해 침수로 인한 직접적 피해 규모가 상당할 뿐만 아니라 침수 발생 후 세균 및 박테리아에 의해 발생하는 수인성 전염병의 유행 등과 같은 2차적 피해 또한 심각한 사회적 비용을 초래할 수 있어 도시유역의 침수 피해를 저감시키기 위한 대책이 절실히 요구되어지고 있다. 도시유역의 침수를 예방하기 위한 대책은 구조적 비구조적 대책으로 구분되어 질 수 있으며 구조적 대책의 경우 침수 피해 예방에 직접적인 효과를 낼 수 있다는 장점이 있으나 대규모 사업예산 및 사업 기간으로 인해 직접적 효과를 보기까지 상대적으로 긴 시간이 필요할 뿐만 아니라 사업 진행 중 대상지역 거주민들의 민원으로 인한 갈등 조정 등으로 인해 사업실행에 어려움을 겪고 있다. 이러한 측면에서 비구조적 대책의 일환인 수치해석을 통한 침수피해 재현 및 침수원인 파악을 통한 구조개선 제안은 구조적 대책의 단점을 보완할 수있는 좋은 대안이 될 수 있다. 도시유역의 경우 비도시유역과 대조적인 차이점으로는 높은 비율의 불투수층, 복잡한 지형, 다수의 인공 구조물 및 배수관망 시스템 등을 들 수 있으며, 침수해석 모형의 정확도를 높이기 위해서는 복잡한 지형의 효율적인 처리가 무엇보다 중요하다. 일반적으로 이용되는 2차원 침수해석 모형들은 직교구조 격자 또는 비구조 격자를 이용하여 지형을 묘사하고 있으며 DEM 자료를 직접 사용하는 직교구조 격자의 경우 지형 데이터 생성이 상대적으로 쉽다는 장점이 있으나 복잡한 지형을 표현하기 위해서는 불필요한 지역까지 높은 해상도를 이용해야 하며 이로 인하여 모의시간이 지나치게 길어지는 문제점이 발생한다. 비구조 격자의 경우 상대적으로 복잡한 도시 유역을 잘 묘사할 수 있다는 장점이 있으나 격자망 생성에 필요한 데이터가 많고 격자망 생성에 지나치게 많은 시간과 노력이 소요된다는 단점이 있다. 따라서 본 연구에서는 위에서 언급한 두 가지 방법의 장점만을 취할 수 있도록 메쉬 세분화 기법을 이용한 2차원 침수해석 모형을 개발 하여 복잡한 지형은 고해상도 메쉬를 이용하여 보다 자세히 묘사하고 상대적으로 복잡하지 않은 지형은 저해상도 메쉬를 이용하여 계산시간을 단축시킬 수 있도록 하였다. 수치해석 기법으로는 엇갈림 격자를 이용하는 Leap-Frog 기법과 유한차분 (Finite difference Method)기법을 이용하였다.

  • PDF

Wind resistance performance of a continuous welding stainless steel roof under static ultimate wind loading with testing and simulation methods

  • Wang, Dayang;Zhao, Zhendong;Ou, Tong;Xin, Zhiyong;Wang, Mingming;Zhang, Yongshan
    • Wind and Structures
    • /
    • v.32 no.1
    • /
    • pp.55-69
    • /
    • 2021
  • Ultrapure ferritic stainless steel provides a new generation of long-span metal roof systems with continuous welding technology, which exhibits many unknown behaviors during wind excitation. This study focuses on the wind-resistant capacity of a new continuous welding stainless steel roof (CWSSR) system. Full-scale testing on the welding joints and the CWSSR system is performed under uniaxial tension and static ultimate wind uplift loadings, respectively. A finite element model is developed with mesh refinement optimization and is further validated with the testing results, which provides a reliable way of investigating the parameter effect on the wind-induced structural responses, namely, the width and thickness of the roof sheeting and welding height. Research results show that the CWSSR system has predominant wind-resistant performance and can bear an ultimate wind uplift loading of 10.4 kPa without observable failures. The welding joints achieve equivalent mechanical behaviors as those of base material is produced with the current of 65 A. Independent structural responses can be found for the roof sheeting of the CWSSR system, and the maximum displacement appears at the middle of the roof sheeting, while the maximum stress appears at the connection supports between the roof sheeting with a significant stress concentration effect. The responses of the CWSSR system are greatly influenced by the width and thickness of the roof sheeting but are less influenced by the welding height.

Multi-material topology optimization for crack problems based on eXtended isogeometric analysis

  • Banh, Thanh T.;Lee, Jaehong;Kang, Joowon;Lee, Dongkyu
    • Steel and Composite Structures
    • /
    • v.37 no.6
    • /
    • pp.663-678
    • /
    • 2020
  • This paper proposes a novel topology optimization method generating multiple materials for external linear plane crack structures based on the combination of IsoGeometric Analysis (IGA) and eXtended Finite Element Method (X-FEM). A so-called eXtended IsoGeometric Analysis (X-IGA) is derived for a mechanical description of a strong discontinuity state's continuous boundaries through the inherited special properties of X-FEM. In X-IGA, control points and patches play the same role with nodes and sub-domains in the finite element method. While being similar to X-FEM, enrichment functions are added to finite element approximation without any mesh generation. The geometry of structures based on basic functions of Non-Uniform Rational B-Splines (NURBS) provides accurate and reliable results. Moreover, the basis function to define the geometry becomes a systematic p-refinement to control the field approximation order without altering the geometry or its parameterization. The accuracy of analytical solutions of X-IGA for the crack problem, which is superior to a conventional X-FEM, guarantees the reliability of the optimal multi-material retrofitting against external cracks through using topology optimization. Topology optimization is applied to the minimal compliance design of two-dimensional plane linear cracked structures retrofitted by multiple distinct materials to prevent the propagation of the present crack pattern. The alternating active-phase algorithm with optimality criteria-based algorithms is employed to update design variables of element densities. Numerical results under different lengths, positions, and angles of given cracks verify the proposed method's efficiency and feasibility in using X-IGA compared to a conventional X-FEM.

Wintertime Extreme Storm Waves in the East Sea: Estimation of Extreme Storm Waves and Wave-Structure Interaction Study in the Fushiki Port, Toyama Bay (동해의 동계 극한 폭풍파랑: 토야마만 후시키항의 극한 폭풍파랑 추산 및 파랑 · 구조물 상호작용 연구)

  • Lee, Han Soo;Komaguchi, Tomoaki;Yamamoto, Atsushi;Hara, Masanori
    • Journal of Korean Society of Coastal and Ocean Engineers
    • /
    • v.25 no.5
    • /
    • pp.335-347
    • /
    • 2013
  • In February 2008, high storm waves due to a developed atmospheric low pressure system propagating from the west off Hokkaido, Japan, to the south and southwest throughout the East Sea (ES) caused extensive damages along the central coast of Japan and along the east coast of Korea. This study consists of two parts. In the first part, we estimate extreme storm wave characteristics in the Toyama Bay where heavy coastal damages occurred, using a non-hydrostatic meteorological model and a spectral wave model by considering the extreme conditions for two factors for wind wave growth, such as wind intensity and duration. The estimated extreme significant wave height and corresponding wave period were 6.78 m and 18.28 sec, respectively, at the Fushiki Toyama. In the second part, we perform numerical experiments on wave-structure interaction in the Fushiki Port, Toyama Bay, where the long North-Breakwater was heavily damaged by the storm waves in February 2008. The experiments are conducted using a non-linear shallow-water equation model with adaptive mesh refinement (AMR) and wet-dry scheme. The estimated extreme storm waves of 6.78 m and 18.28 sec are used for incident wave profile. The results show that the Fushiki Port would be overtopped and flooded by extreme storm waves if the North-Breakwater does not function properly after being damaged. Also the storm waves would overtop seawalls and sidewalls of the Manyou Pier behind the North-Breakwater. The results also depict that refined meshes by AMR method with wet-dry scheme applied capture the coastline and coastal structure well while keeping the computational load efficiently.

An evaluation of wall functions for RANS computation of turbulent flows (난류 흐름의 RANS 수치모의를 위한 벽함수 성능 평가)

  • Yoo, Donggeun;Paik, Joongcheol
    • Journal of Korea Water Resources Association
    • /
    • v.53 no.1
    • /
    • pp.1-13
    • /
    • 2020
  • The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.

Development of Structure Analysis Program for Jointed Concrete Pavement Applying Load Discretization Algorithm (하중변환 알고리듬을 적용한 줄눈 콘크리트 포장해석 프로그램 개발)

  • Yun, Tae-Young;Kim, Ji-Won;Cho, Yoon-Ho
    • International Journal of Highway Engineering
    • /
    • v.5 no.4 s.18
    • /
    • pp.1-11
    • /
    • 2003
  • Recently, the new pavement design method considering Korean environment and the specification for improving performance of pavement are being developed in Korea. The Jointed Concrete Pavement Program Applying Load Discretization Algorithm (called HEART-JCP) is one of the results of Korea Pavement Research Project in Korea. HEART-JCP program is developed to analyze various loading condition using the load discretization algorithm without mesh refinement. In addition, it can be modified easily into multi-purpose concrete pavement nidyses program because of the modularized structure characteristic of HEART-JCP. The program consists of basic program part and load discretization part. In basic program part, the displacement and stress are computed in the concrete slab, sub-layer, and dowel bar, which are modeled with plate/shell element, spring element and beam element. In load discretization program part, load discretization algorithm that was used for the continuum solid element is modified to analyze the model with plate and shell element. The program can analyze the distributed load, concentrated load, thermal load and body load using load discretization algorithm. From the result of verification and sensitivity study, it was known that the loading position, the magnitude of load, and the thickness of slab were the major factors of concrete pavement behavior as expected. Since the result of the model developed is similar to the results of Westergaard solution and ILLISLAB, the program can be used to estimate the behavior of jointed concrete pavement reasonably.

  • PDF

Semi-Automatic Method for Constructing 2D and 3D Indoor GIS Maps based on Point Clouds from Terrestrial LiDAR (지상 라이다의 점군 데이터를 이용한 2차원 및 3차원 실내 GIS 도면 반자동 구축 기법 개발)

  • Hong, Sung Chul;Jung, Jae Hoon;Kim, Sang Min;Hong, Seung Hwan;Heo, Joon
    • Journal of Korean Society for Geospatial Information Science
    • /
    • v.21 no.2
    • /
    • pp.99-105
    • /
    • 2013
  • In rapidly developing urban areas that include high-rise, large, and complex buildings, indoor and outdoor maps in GIS become a basis for utilizing and sharing information pertaining to various aspects of the real world. Although an indoor mapping has gained much attentions, research efforts are mostly in 2D and 3D modeling of terrain and buildings. Therefore, to facilitate fast and accurate construction of indoor GIS, this paper proposes a semi-automatic method consisting of preprocessing, 2D mapping, and 3D mapping stages. The preprocessing is designed to estimate heights of building interiors and to identify noise data from point clouds. In the 2D mapping, a floor map is extracted with a tracing grid and a refinement method. In the 3D mapping, a 3D wireframe model is created with heights from the preprocessing stage. 3D mesh data converted from noise data is combined with the 3D wireframe model for detail modeling. The proposed method was applied to point clouds depicting a hallway in a building. Experiment results indicate that the proposed method can be utilized to construct 2D and 3D maps for indoor GIS.

Stereo Image-based 3D Modelling Algorithm through Efficient Extraction of Depth Feature (효율적인 깊이 특징 추출을 이용한 스테레오 영상 기반의 3차원 모델링 기법)

  • Ha, Young-Su;Lee, Heng-Suk;Han, Kyu-Phil
    • Journal of KIISE:Computer Systems and Theory
    • /
    • v.32 no.10
    • /
    • pp.520-529
    • /
    • 2005
  • A feature-based 3D modeling algorithm is presented in this paper. Since conventional methods use depth-based techniques, they need much time for the image matching to extract depth information. Even feature-based methods have less computation load than that of depth-based ones, the calculation of modeling error about whole pixels within a triangle is needed in feature-based algorithms. It also increase the computation time. Therefore, the proposed algorithm consists of three phases, which are an initial 3D model generation, model evaluation, and model refinement phases, in order to acquire an efficient 3D model. Intensity gradients and incremental Delaunay triangulation are used in the Initial model generation. In this phase, a morphological edge operator is adopted for a fast edge filtering, and the incremental Delaunay triangulation is modified to decrease the computation time by avoiding the calculation errors of whole pixels and selecting a vertex at the near of the centroid within the previous triangle. After the model generation, sparse vertices are matched, then the faces are evaluated with the size, approximation error, and disparity fluctuation of the face in evaluation stage. Thereafter, the faces which have a large error are selectively refined into smaller faces. Experimental results showed that the proposed algorithm could acquire an adaptive model with less modeling errors for both smooth and abrupt areas and could remarkably reduce the model acquisition time.

Applicability Evaluation of Flood Inundation Analysis using Quadtree Grid-based Model (쿼드트리 격자기반 모형의 홍수범람해석 적용성 평가)

  • Lee, Dae Eop;An, Hyun Uk;Lee, Gi Ha;Jung, Kwan Sue
    • Journal of Korea Water Resources Association
    • /
    • v.46 no.6
    • /
    • pp.655-666
    • /
    • 2013
  • Lately, intensity and frequency of natural disasters such as flood are increasing because of abnormal climate. Casualties and property damages due to large-scale floods such as Typhoon Rusa in 2002 and Typhoon Maemi in 2003 rapidly increased, and these show the limits of the existing disaster prevention measures and flood forecasting systems regarding irregular climate changes. In order to efficiently respond to extraordinary flood, it is important to provide effective countermeasures through an inundation model that can accurately simulate flood inundation patterns. However, the existing flood inundation analysis model has problems such as excessive take of analysis time and accuracy of the analyzed results. Therefore, this study conducted a flood inundation analysis by using the Gerris flow solver that uses quadtree grid, targeting the Baeksan Levee in the Nakdong River Basin that collapsed because of a concentrated torrential rainfall in August, 2002. Through comparisons with the FLUMEN model that uses unstructured grid among the existing flood inundation models and the actual flooded areas, it determined the applicability and efficiency of the quadtree grid-based flood inundation model of the Gerris flow solver.