• The Journal of Engineering Geology
• /
• v.33 no.4
• /
• pp.573-586
• /
• 2023

Joints or weak planes can induce anisotropy in the strength and deformability of fractured rock masses. Comprehending this anisotropic behavior is crucial to engineering geology. This study used plaster as a friction material to mold specimens with a single joint. The strength and deformability of the specimens were measured in true triaxial compression tests. The measured results were compared with three-dimensional numerical analysis based on the distinct element method, conducted under identical conditions, to assess the reliability of the modeled values. The numerical results highlight that the principal stress conditions in the field, in conjunction with joint orientations, are crucial factors to the study of the strength and deformability of fractured rock masses. The strength of a transversely isotropic rock mass derived numerically considering changes in the dip angle of the joint notably increases as the intermediate principal stress increases. This increment varies depending on the dip of the joint. Moreover, the interplay between the dip direction of the joint and the two horizontal principal stress directions dictates the strength of the transversely isotropic rock mass. For a rock mass with two joint sets, the set with the steeper dip angle governs the overall strength. If a rock bridge effect occurs owing to the limited continuity of one of the joint sets, the orientation of the set with longer continuity dominates the strength of the entire rock mass. Although conventional three-dimensional failure criteria for fractured rock masses have limited applicability in the field, supplementing them with numerical analysis proves highly beneficial.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2000.04a
• /
• pp.16-16
• /
• 2000

충격파를 포함하는 초음속 유동을 해석하는 수치해법 중에서 많이 사용되어진 것은 엄밀 및 근사 리만 해법과 플럭스 분할 기법들로서 이들은 Euler 방정식에 기반을 두고 선형 또는 비선형파의 상호작용을 풍상 차분법으로 기술하는 방법들이다. 이러한 수치기법들은 과거 광범위하게 사용되어 왔으나 최근 여러 가지 단점이 발견되었다. 이와 같은 문제점을 극복하고자 입자의 통계적인 운동을 기술하는 기체 운동론에 근거하여 BGK 수치기법이 제시되었다. 이는 비충돌 볼츠만 방정식으로부터 입자의 수준에서 플럭스 분할 기법 형태의 풍상차분법을 구현하는 것으로 볼츠만 방정식의 충돌항을 BGK 모델로 대치하고 이것의 적분해로부터 수치 플럭스를 구한다. 이 수치기법은 기존의 리만해법에 비하여 수치적으로나 물리적으로 매우 타당한 성질인 강건성, 정확성, 엔트로피 조건, 양수보존성 등을 가지고 있음이 밝혀졌다. 이와 같은 수치기법을 사용하여 로켓 노즐 내의 아음속, 천이음속, 초음속에서의 유동장 해석을 위한 프로그램을 작성하였다. 시간 적분에 대하여는 정상 상태의 계산을 위하여 내재적 시간 적분 방법을 사용하였으며, 공간 이산화 방법으로는 임의의 제어체적에 대하여 적분형 보존 방정식을 적용하는 유한 체적법을 사용하였다. 초음속 입구 유동과 출구에서 초음속과 저음속 유동의 두가지 경우를 고려하여 얻은 결과를 기존의 연구 결과와 비교하여 본 결과 잘 일치하였다. 입구 유동이 저음속이고 출구 유동이 초음속인 경우에 대하여도 해석결과가 실험결과와 잘 일치하였다. 상대적으로 낮은 온도, 압력 조건과 높은 온도, 압력 조건을 가지는 고체 로켓 모터 노즐 내의 유동을 해석하였다. 이들 해석 결과를 전압, 전온도로 표준화시킨 결과 서로 일치하였으며, 파라서 저온, 저압에서 얻은 결과도 표준화시킬 경우, 고온, 고압에서도 사용될 수 있음을 알 수 있었다.의 영향에 초점을 맞추었다.다고 판단되며 배기 가스 자체에 대기 공기중에 함유되어 있던 습기가 얼어붙는(Icing화) 문제가 발생하기 때문에 배기가스의 Icing을 방지하기 위하여 압축기 끝단에서 공기를 추출하여 배기부분에 송출할 필요성이 있는 것으로 판단되었다. 출구가스의 기체 유동속도가 매우 빠르므로 (100-l10m.sec) 이를 완화하기 위한 디퓨저의 설계가 요구된다고 판단된다. 또 연소기 후방에 물을 주입하는 경우 열교환기 및 기타 부분품에 발생할 수 있는 부식 및 열교환 효율 저하도 간과할 수 없는 문제로 파악되었다. 이러한 기술적 문제가 적절히 해결되는 경우 비활성 가스 제너레이터는 민수용으로는 대형 빌딩, 산림, 유조선 등의 화재에 매우 적절히 사용되어 질 수 있을 뿐 아니라 군사적으로도 군사작전 중 및 공군 기지의 화재 그리고 지하벙커에 설치되어 있는 고급 첨단 군사 장비 등의 화재 뿐 아니라 대간첩작전 등에 효과적으로 활용될 수 있을 것으로 판단된다.가 작으며, 본 연소관에 충전된 RDX/AP계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..7%), 혈액투석, 식도부분절제술 및 위루술·위회장문합술을 시행한 경우가 각 1례(2.9%)씩이었다. 13) 심각한 합병증은 9례(26.5%)에서 보였는데 그중 식도협착증이 6례(17.6%), 급성신부전증 1례(2.9%), 종격동기흉과 폐염이 병발한 경우와 폐염이 각 1례(2.9%)였다. 14) 식도경 시행회수는 1회가 17례(54.8%), 2회가 9례(29.0%), 3회 이상이 5례(16.1%)였다.EX>$IC_{50}$/ 값이 210 $\mu\textrm{g}$<

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.6
• /
• pp.84-93
• /
• 2018

Three-dimensional (3D) numerical simulations have been carried out to study how coolant injection conditions influence the cooling efficiency and projectile ejection performance in a mixture-gas ejection system (or gas-steam launch system). The 3D single-phase computational model was verified using a 1D model constructed with reference to the previous research and then a two-phase flow computation simulating coolant injection on to hot gas was performed using a DPM (Discrete Phase Model). As a result of varying the coolant flow rate and number of injection holes, cooling efficiency was improved when the number of injection holes were increased. In addition, the change of the coalescence frequency and spatial distribution of coolant droplets caused by the injection condition variation resulted in a change of the droplet diameter, affecting the evaporation rate of coolant. The evaporation was found to be a critical factor in the design optimization of the ejection system by suppressing the pressure drop while the temperature decreases inside the breech.

• Economic and Environmental Geology
• /
• v.44 no.2
• /
• pp.161-170
• /
• 2011

The multiphase flow simulator, CHEMPS, was developed based on the fractional flow approach reported in the petroleum engineering literature considering fully three phase flow in physically and chemically heterogeneous media. It is a extension of MPS developed by Suk and Yeh (2008) to include the effect of wettability on the migration of NAPL. The fractional flow approach employs water, total liquid saturation and total pressure as the primary variables. Most existing models are limited to two-phase flow and specific boundary conditions when considering physically heterogeneous media. In addition, these models focused mainly on the water-wet media. However, in a real system, variations in wettability between water-wet and oil-wet media often occur. Furthermore, the wetting of porous media by oil can be heterogeneous, or fractional, rather than uniform due to the heterogeneous nature of the subsurface media and the factors that affect the wettability. Therefore, in this study, the chemically heterogeneous media considering fractional wettability as well as physically heterogeneous media were simulated using CHEMPS. In addition, the general boundary conditions were considered to be a combination of two types of boundaries of individual phases, flux-type and Dirichlet type boundaries.

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

우수 관거 시스템에서 맨홀의 설치 시 연결관 내부와 맨홀의 내부는 여러 가지 수리학적 조건이 다르므로 수두손실의 발생이 필연적일 수밖에 없다. 현재 계획 또는 설계단계에서 수행되고 있는 관거 시설의 수리계산에는 연결관 내에서의 마찰손실만을 고려하여 설계를 수행하고 있으며, 맨홀에서의 수두 손실은 거의 대부분 고려되지 않고 있다. 단지 맨홀에서의 수두손실을 저감하기 위하여 하수도시설기준(환경부, 2005)상의 단차 및 인버트 규정만 있을 뿐, 우수 관거 설계에 직접적으로 필요한 적절한 맨홀의 손실계수가 제시되지 않고 있는 실정이다. 국외에서는 축소 수리 모형을 이용한 실험과 수치해석 기법 등을 이용하여 맨홀에서의 손실계수를 산정하는 연구가 꾸준히 진행되어 왔으나 국내에서는 맨홀의 손실계수 산정에 관한 연구가 미흡한 실정이며, 더욱이 맨홀에서의 손실계수 산정을 위한 상사성 적용에 관한 연구는 전무한 실정이다. 그러므로 본 연구에서는 맨홀의 축척 변화에 따른 손실계수의 변화를 분석하기 위하여 하수도시설기준(환경부, 2005)의 특 1호(사각형) 맨홀을 각각 1/2과 1/5로 축소 제작하고, 수리실험 장치를 제작하였다. Froude 상사 법칙을 적용하여 1/2의 축소 모형의 실험 조건을 1/5 축소 모형의 값으로 환산하였으며, 각 축소 모형에 대한 수리 실험을 실시하였다. 과부하된 맨홀의 손실계수를 예측하는데 Froude 상사법칙의 사용 가능성을 확인 하였으며, 1/2 축소 모형과 1/5 축소 모형에서 산정된 손실계수 값이 0.45로 일치하고 있으므로 우수 관거 시스템의 맨홀 설계 시, 축소 수리 모형실험에서 산정된 손실계수의 직접적인 적용이 가능하다고 판단된다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.1-7
• /
• 2011

In this paper, based on an isogeometric approach, we have developed a shape design optimization method for plane elasticity problems subjected to design-dependent loads. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry. In an isogeometric analysis, however, the geometric properties are already embedded in the B-spline basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. The solution space for the response analysis can be represented in terms of the same NURBS basis functions to represent the geometry, which yields a precise analysis model that exactly represents the normal and curvature depending on the applied loads. A continuum-based isogeometric adjoint sensitivity is extensively derived for the plane elasticity problems under the design-dependent loads. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.404-404
• /
• 2011

최근 기상변동성 증가와 극치수문사상의 발생빈도 증가로 인한 기상재해가 빈번하게 일어나고 있다. 이러한 기상현상으로 인한 재해의 예방을 위해서 사전에 위험을 인지하고 그 규모를 예측할 수 있는 여러 기법들이 기상레이더 또는 수치예보자료 등을 이용하여 개발 및 적용되고 있다. 이 과정에서 해결해야 할 여러 문제점들이 있는데, 우선 수치예보자료 또는 기상레이더자료를 종관기상관측소 및 자동기상관측지점의 지상관측 강수량과 연계하여 평가하는 과정이 필요하고, 현재시점에 형성되어 있는 강우장의 공간 이동 예측 기법이 확보되어야 할 것이다. 전북지역은 게릴라성 집중호우가 빈번한 산악형 강수와 산지유역의 급한 하천경사가 맞물려 인명 및 재산피해가 매년 발생하고 있으며, 과거 돌발홍수가 발생한 사례가 있어 이상기후 및 기후변화로 인한 홍수 위험도가 커질 것으로 전망되고 있다. 본 연구는 전라북도의 기상재해 예측모형 개발을 위한 사전 분석과정으로 전라북도지역에서 관측된 기존의 대규모 강수사상을 이용한 강수사상의 특성 분류 및 관측소간 공간상관성을 분석하는데 목적을 두고 있다. 강수사상의 특성분류를 통해 강수 발생형태에 따른 기상학적 영향인자, 강수의 발생량 및 이동특성 예측의 정도를 향상시킬 수 있으며, 분류 기법으로 SVM(support vector machine)을 이용한 자동분류를 적용한다. 또한 관측소간 공간상관성 분석을 위하여 각 관측소 강수량간의 조건부 확률을 이용한다. 예로써 부안관측소에 강수가 발 생했을 때, 부안관측소의 강수량 조건에 의한 전주관측소 강수량 확률을 다음과 같이 구성할 수 있다. �揚滑斂�수량�咀刮활�수량��. 공간상관성 분석과정에서 관측소간 강수 이동시간에 따른 강수 발생 시간의 차이 또한 고려하며, 과거 기상관측 자료의 분석을 통해 전라북도지역의 관측소간 강수발생의 공간적 상관성을 규명하고, 단기예측 모델 개발을 위한 기초자료로 활용할 수 있을 것이다. 또한, 기후변화시나리오에 의한 미래 강수량의 지역적 상세화 과정에도 본 연구를 통한 결과를 이용할 수 있을 것이라 판단된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.44-44
• /
• 2015

저수지나 하천 사면에서 발생하는 산사태와 토석류는 저수지와 하천 수체에 충격을 가한다. 이로 인해 발생하는 수면 충격파는 전파되어 반대편 제방으로 파의 처오름 또는 댐 제체위로의 물넘이로 큰 피해를 줄 수 있다. 최근 외국에서는 2차원 충격파 생성 및 전파의 기본 과정을 구명하기 위한 실험적 연구가 이뤄지고 있으며, 이들 연구들은 충격파의 발생과 전파, 사면활동 물질과 수체의 상호작용 그리고 자유 수면과 유속분표의 발달에 대한 자세한 관측 자료를 제시하고 있다. 아울러 충격파에 영향을 주는 지배 매개변수를 제시하고 있다. 하지만, 이러한 실험적 연구의 최근 진보에도 불구하고, 이들 지배 매개변수를 고려한 충격파 지배공식들은 대상 지역의 복잡한 바닥 지형이나, 평면적 지형 변화를 단순한 추정치로만 고려하게 된다. 따라서 복잡한 지형조건에서 토석류와 수체의 상호작용과 수면 충격파의 전파를 합리적으로 해석하는 데는 한계가 있다. 이 경우 수치모델링 기법을 대안으로 적용할 수 있으나, 수치모델링은 수면에서 충격파의 전파와 수중에서 토석류의 전파를 동시에 모의해야 하고, 뉴턴 유체와 비뉴턴 유체의 특성을 동시에 고려해야하므로 수치해석 연구자들에게는 하나의 큰 도전사항이다. 이 연구는 경계면 포착기법을 이용한 계산유체동력학 기법을 이용하여 사면활동과 이로 인한 정지 수역에서의 충격파의 발생 및 전파를 재현하기 위한 수치 모델링 기법을 개발하는 것이 목적이다. 사면활동과 수면의 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축 기법을 이용한다. 지배 방정식은 비압축성(incompressible) 질량 보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 큰와 모의 계열의 난류 모델링 기법을 적용하여 충격파의 전파와 붕괴에 대한 난류의 영향을 고려하였다. 토석류는 비뉴턴 흐름저항 관계식을 적용하여 그 흐름특성을 재현하였다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. 외국의 연구자들이 관측하여 제시한 길이 11 m 그리고 폭 0.5 m의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.

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

통상 방조제 배수갑문의 방류능력은 위어공식 또는 오리피스공식을 이용하여 산정하는데, 이 경우 지형특성, 배수갑문 형상에 따른 흐름의 간섭현상, 유입부와 유출부의 형상 등이 방류량에 미치는 영향을 고려하지 못한다. 본 연구에서는 도류벽, 배수문, 물받이 등 방조제 배수갑문의 형상과 배치가 방류량에 미치는 영향을 상용 프로그램인 FLOW-3D를 이용하여 정량적으로 해석하였다. 이를 통해 배수갑문의 방류능력과 유황을 개선할 수 있는 방안을 도출하였다. 본 연구에서는 시화조력발전소를 대상으로 배수갑문의 방류능력 개선에 3차원 수치모의가 효과적으로 적용될 수 있음을 보였다. 본 연구에서 도출한 주요 결론은 다음과 같다. 1) 유출부측의 물받이길이를 40 m 증가시킴에 따라 전체 방류량은 계획안에 비해 약 $10\% 증가하는 것으로 나타났다. 2) 물받이 끝과 원지반의 연결부 사면을 1:1에서 1:5의 완경사로 변화시킴에 따라 전체 방류량은 약 $2\%$ 증가하는 것으로 나타났다. 3) 배수문과 수차발전 구조물 사이의 유선형 연결구조물을 제거함에 따라 전체 방류량은 약 $3\%$ 증가하는 것으로 나타났다. 4) 도류벽의 접근각도를 $10^{\circ}$ 감소시키거나 증가시킴에 따라 전체 방류량은 약 $5\% 감소 또는 증가하는 것으로 나타났다. 본 연구는 배수갑문의 설계시 방류능력 개선을 위해서는 수리학적 검토가 필요하며 수치모형실험이 수리모형실험과 더불어 유용한 해석도구로 이용될 수 있음을 보인 것으로, 이후 관련 구조물의 설계시 참고자료로 이용 가능할 것으로 사료된다.다. 실험 결과, Escarameia와 May가 제안한 공식을 더 확장하여 적용할 수 있는 실험 공식으로 개선하였으며 다양한 조건에 대한 실험을 수행하여 보다 정밀한 공식으로 개선할 수 있었다.$10,924m^3/s$ 및 $10,075m^3/s$로서 실험 I의 $2,757m^3/s$에 비해 통수능이 많이 개선되었음을 알 수 있다.함을 알 수 있다. 상수관로 설계 기준에서는 관로내 수압을 $1.5\~4.0kg/cm^2$으로 나타내고 있는데 $6kg/cm^2$보다 과수압을 나타내는 경우가 $100\%$로 밸브를 개방하였을 때보다 $60\%,\;80\%$ 개방하였을 때가 더 빈번히 발생하고 있으므로 대상지역의 밸브 개폐는 $100\%$ 개방하는 것이 선계기준에 적합한 것으로 나타났다. 밸브 개폐에 따른 수압 변화를 모의한 결과 밸브 개폐도를 적절히 유지하여 필요수량의 확보 및 누수방지대책에 활용할 수 있을 것으로 판단된다.8R(mm)(r^2=0.84)$로 지수적으로 증가하는 경향을 나타내었다. 유거수량은 토성별로 양토를 1.0으로 기준할 때 사양토가 0.86으로 가장 작았고, 식양토 1.09, 식토 1.15로 평가되어 침투수에 비해 토성별 차이가 크게 나타났다. 이는 토성이 세립질일 수록 유거수의 저항이 작기 때문으로 생각된다. 경사에 따라서는 경사도가 증가할수록 증가하였으며 $10\% 경사일 때를 기준으로 $Ro(mm)=Ro_{10}{\times}0.797{\times}e^{-0.021s(\

• Ecology and Resilient Infrastructure
• /
• v.8 no.4
• /
• pp.245-252
• /
• 2021

A turbidity current in a river and a lake occurs due to diverse nutrient loading including suspended sediment in sediment runoff, which affects water withdrawal and river environments. We developed one dimensional time-variant numerical model based on Python for the Nagdonggang mainstream. We examined the numerical stability and the applicability of the model by performing the simulation of quasi-steady flow in non-flooding for three cases, which are different according to the point and the amount of turbidity inflows in the Nagdonggang upstream and a tributary. The result was reasonable in the respect of the conservation of matter. The model will facilitate to simulate a large river if we can secure the data of turbidity variations in a target river reach or measured points in a field.