• 한국지하수토양환경학회지:지하수토양환경
• /
• 제13권2호
• /
• pp.54-61
• /
• 2008

수평 혹은 경사 형태 특수 정호 양수량에 대한 시공간적 수위 강하를 지하수 수치 모델링을 활용하여, 평가하였다. 지하수 수치 모델링은 상용 프로그램인 FEFLOW(version 5.1)의 1차원 선형 불연속 특징 요소를 활용하여 수행되었으며, 수치해의 검증을 위해 Zhan과 Zlotnik(2002)이 제안한 연속된 점 형태 배출원 배열 방식 준 해석해와 비교하였다. 비교 검증 결과, 수치해와 준해석해는 최대 수위 강하가 나타나는 양수 최인접부를 제외하고는 거의 일치한 형태를 보여주었다. 검증된 수치적 방법을 이용하여, 강변여과 방식 취수가 검토되는 현장에 대한 수위강하를 정량적으로 평가할 수 있었다.

• Structural Engineering and Mechanics
• /
• 제8권5호
• /
• pp.513-529
• /
• 1999

A concept of hierarchical modeling, the newest modeling technology, has been introduced in early 1990's. This new technology has a great potential to advance the capabilities of current computational mechanics. A first step to implement this concept is to construct hierarchical models, a family of mathematical models sequentially connected by a key parameter of the problem under consideration and have different levels in modeling accuracy, and to investigate characteristics in their numerical simulation aspects. Among representative model problems to explore this concept are elastic structures such as beam-, arch-, plate- and shell-like structures because the mechanical behavior through the thickness can be approximated with sequential accuracy by varying the order of thickness polynomials in the displacement or stress fields. But, in the numerical, analysis of hierarchical models, two kinds of errors prevail, the modeling error and the numerical approximation error. To ensure numerical simulation quality, an accurate estimation of these two errors is definitely essential. Here, a local a posteriori error estimator for elastic structures with thin domain such as plate- and shell-like structures is derived using the element residuals and the flux balancing technique. This method guarantees upper bounds for the global error, and also provides accurate local error indicators for two types of errors, in the energy norm. Compared to the classical error estimators using the flux averaging technique, this shows considerably reliable and accurate effectivity indices. To illustrate the theoretical results and to verify the validity of the proposed error estimator, representative numerical examples are provided.

• 지질학회지
• /
• 제54권6호
• /
• pp.683-694
• /
• 2018

섭입은 지구의 물질 및 에너지 순환에서 중요한 역할을 담당할 뿐만 아니라 우리의 삶에 밀접한 지질 현상인 지진과 호화산을 발생시키므로 많은 연구가 이루어져 왔다. 그 중에서 맨틀 내부의 섭입해양판처럼 우리가 직접 관찰할 수 없는 곳에서 발생하는 지질 현상에 대한 정량적 연구에 컴퓨터 수치모델링이 널리 이용되어 왔다. 이 연구에서는 다양한 연구진들에 의해 사용되고 있는 섭입대 수치모델링의 벤치마크를 수행하였다. 섭입대 수치모델링을 위하여 유한요소법 기반 상용 소프트웨어인 콤솔 멀티피직스를 사용하였으며 계산된 결과는 과거 수행된 벤치마크 결과와 잘 일치하였다.

• 방사성폐기물학회지
• /
• 제20권4호
• /
• pp.429-453
• /
• 2022

Numerical modeling and scenario composition are needed to characterize the geological environment of the disposal site and analyze the long-term evolution of natural barriers. In this study, processes and features of the hydro-mechanical behavior of natural barriers were categorized and represented using the interrelation matrix proposed by SKB and Posiva. A hydro-mechanical coupled model was evaluated for analyzing stress field changes and fracture zone re-activation. The processes corresponding to long-term evolution and the hydro-mechanical mechanisms that may accompany critical processes were identified. Consequently, practical numerical methods could be considered for these geological engineering issues. A case study using a numerical method for the stability analysis of an underground disposal system was performed. Critical stress distribution regime problems were analyzed numerically by considering the strata's movement. Another case focused on the equivalent continuum domain composition under the upscaling process in fractured rocks. Numerical methods and case studies were reviewed, confirming that an appropriate and optimized modeling technique is essential for studying the stress state and geological history of the Korean Peninsula. Considering the environments of potential disposal sites in Korea, selecting the optimal application method that effectively simulates fractured rocks should be prioritized.

• Geomechanics and Engineering
• /
• 제33권4호
• /
• pp.353-367
• /
• 2023

The construction of a protective embankment is a suitable strategy to stop and control high-energy rock blocks' impacts during the rockfall phenomenon. In this paper, based on the discrete element numerical method, by modeling an existing embankment reinforced with geogrid, its stability status under the impact of a rock block with two types of low and high kinetic energy, namely 2402 and 4180 kJ, respectively, has been investigated. The modeling results show that the use of geogrid has caused the displacement in the front and back of the embankment to decrease by more than 30%. In this case, the reinforced embankment has stopped the rock block earlier. The displacements obtained from the DEM modeling are compared with the displacements measured from an actual practical experiment to evaluate the results' validity. Comparison between the results shows that the displacement values are close together, while the maximum percentage error in previous studies by an analytical method and the finite element method was 76.4% and 36.6%, respectively. Therefore, the obtained results indicate the discrete numerical method's high ability compared to other numerical and analytical methods to simulate and design the geogrid-reinforced soil embankment under natural disasters such as rockfall with a minor error.

• 한국해양공학회지
• /
• 제37권2호
• /
• pp.68-79
• /
• 2023

Even though submerged breakwater reduces incident wave energy, it redistributes the coastal area's wave-induced current, sediment transport, and morphological change. This study examines the coastal hydrodynamics and the morphological response of a wave-dominated beach with submerged breakwaters installed through field observation and quasi-3D numerical modeling. The pre-and post-storm bathymetry, water level, and offshore wave under storm forcing were collected in Bongpo Beach on the East coast of Korea and used to analyze the coastal hydrodynamic response. Four vertically equidistant layers were used in the numerical simulation, and the wave-induced current was examined using quasi-3D numerical modeling. The shore normal incident wave (east-northeast) generated strong cross-shore and longshore currents toward the hinterland of the submerged breakwater. However, the oblique incident wave (east-southeast) induced the southeastward longshore current and the sedimentation in the northeast area of the beach. The results suggested that the incident wave direction is a significant factor in determining the current and sediment transport patterns in the presence of the submerged breakwaters. Moreover, the quasi-3D numerical modeling is more appropriate for estimating the wave transformation, current, and sediment transport pattern in the coastal area with the submerged breakwater.

• Geomechanics and Engineering
• /
• 제31권3호
• /
• pp.319-328
• /
• 2022

Ensuring the integrity of a country's infrastructure is necessary to protect surrounding communities in case of disaster. Embankment dam systems across the US are an essential component of infrastructure, referred to as lifeline structures. Embankment dams are crucial to the survival of life and if these structures were to fail, it is imperative that states be prepared. Southern California is particularly concerned with the stability of embankment dams due to the frequent seismic activity that occurs in the state. The purpose of this study was to create a numerical model of an existing embankment dam simulated under seismic loads using previously recorded data. The embankment dam that was studied in Los Angeles, California was outfitted with accelerometers provided by the California Strong Motion Instrumentation Program that have recorded strong motion data for decades and was processed by the Center for Engineering Strong Motion Data to be used in future engineering applications. The accelerometer data was then used to verify the numerical model that was created using finite element modeling software RS2. The results from this study showed Puddingstone Dam's simulated response was consistent with that experienced during previous earthquakes and therefore validated the predicted behavior from the numerical model. The study also identified areas of weakness and instability on the dam that posed the greatest risk for its failure. Following this study, the numerical model can now be used to predict the dam's response to future earthquakes, develop plans for its remediation, and for emergency response in case of disaster.

• 한국도로학회논문집
• /
• 제16권2호
• /
• pp.43-52
• /
• 2014

PURPOSES : To solve problems in current compaction control DCPT(Dynamic Cone Penetrometer Test), highly correlated with various testing methods, simple, and economic is being applied. However, it、s hard to utilize DCPT results due to the few numerical analyses for DCPT have been performed and the lack of data accumulation. Therefore, this study tried to verify the validation of numerical modeling for DCPT by comparing and analyzing the results of numerical analyses with field tests. METHODS: The ground elastic modulus and PR(Penetration Rate) value were estimated by using PFC(Particle Flow Code) 3D program based on the discrete element method. Those values were compared and analyzed with the result of field tests. Also, back analysis was conducted to describe ground elastic modulus of field tests. RESULTS : Relative errors of PR value between the numerical analyses and field tests were calculated to be comparatively low. Also, the relationship between elastic modulus and PR value turned out to be similar. CONCLUSIONS : Numerical modeling of DCPT is considered to be suitable for describing field tests by carrying out numerical analysis using PFC 3D program.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권3호
• /
• pp.385-392
• /
• 2018

This paper presents 2D numerical modeling to calculate ship-ice interactions that occur when an icebreaker advances into ice-covered water. The numerical model calculates repeated icebreaking of an ice plate and removal of small ice floes. The icebreaking of the ice plate is calculated using a ship-ice contact detection technique and fluid-structural interaction of ice plate bending behavior. The ship-ice interactions in small ice floes are calculated using a physically based modeling with 3DOF rigid body equations. The ice plate is broken in crushing, bending, and splitting mode. The ice floes drift by wind or current and by the force induced by the ship-ice interaction. The time history of ice force and ice floe distribution when an icebreaker advances into the ice-covered water are obtained numerically. Numerical results demonstrate that the time history of ice force and distribution of ice floes (ice channel width) depend on the ice floe size, ship motion and ice drifting by wind or current. It is shown that the numerical model of ship maneuvering in realistic ice conditions is necessary to obtain precise information about the ship in ice-covered water. The proposed numerical model can be useful to provide data of a ship operating in ice-covered water.

• 지질공학
• /
• 제26권4호
• /
• pp.461-471
• /
• 2016

수주지열정의 지하수 유동 및 지열 이동 현상을 모사할 수 있는 수치 모델링 기법을 제시하기 위하여 열-수리학적 모델인 TOUGH2를 기반으로 한 모듈을 개발하였고 두 유형의 수주지열정에서의 냉방 및 난방 운영을 수치 모의하였다. 수치 모델링 결과들은 수주지열정의 냉난방 운영이 주변 지층 및 지하수와의 열교환 및 유체 혼합에 크게 영향을 받으며, 특히 PVC 또는 PE 파이프 내부와 수주지열정 내부에서 서로 반대방향으로 유동하는 유입수와 유출수 간의 열교환이 중요한 요소임을 보여준다. 또한 본 연구에서 개발 및 제안된 수치 모델링 기법이 수주지열정의 내부 구조 및 유체 유동과 열 이동 현상을 합리적으로 모사할 수 있음을 보여준다. 이러한 수치 모델링 기법은 수주지열정의 설계, 시공 및 운영 단계에서 열교환 성능을 정량적으로 분석할 때 유용하게 활용될 수 있을 것으로 기대된다.