• Journal of the Korean Society of Propulsion Engineers
• /
• v.23 no.4
• /
• pp.81-91
• /
• 2019

The 3-dimensional grain burn-back analysis is performed using the face offsetting method for calculating the solid rocket motor performance. The grain burning configuration analysis is a moving surface problem that calculates the regression of the burning surface. In the previous study, various moving interface analysis methods were applied for the grain burn-back analysis, but the results were imperfect. In this study, a 3-dimensional grain burn-back analysis module is developed using the face offsetting method, which combines the advantages of the existing moving interface analysis methods to increase the accuracy and robustness. As a result, the face offsetting method is proved to be efficient for the grain burn-back analysis.

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

해안 대수층은 해수와 담수가 공존하는 지역으로 상대적으로 밀도가 큰 해수가 대수층의 담수 아래에 쐐기형태로 존재하게 된다. 이러한 쐐기형태의 해수와 담수의 경계면은 압력경도의 평형에 의해 경계면이 유지되며, 해수면 또는 지하수위가 변동할 경우 해수-담수 경계면의 균형이 무너짐과 더불어 압력경도의 평행이 이루어질 때 까지 해수-담수 경계면의 이동이 계속 진행된다. 수위 변화의 주요 원인으로는 지구온난화 및 기후변화로 인한 지속적인 해수면 상승과 도서지역의 인구증가 및 산업화로 인한 무분별한 지하수의 사용 등에 의한 지하수위 저하 등을 꼽을 수 있다. 이와 같은 원인으로 해안 및 도서지역에서는 해안 대수층의 해수침투거리가 증가하여 지하수 이용에 큰 어려움을 겪고 있다. 이에 해안 대수층의 해수침투 범위 및 거리를 추정하기 위한 많은 연구들이 다양한 분야에서 지속해서 이루어지고 있지만, 서로 밀도가 다른 해수와 담수가 공존하는 해안 대수층 내의 수리특성을 명확히 파악하기에는 아직까지 미흡한 점들이 많다. 과거에는 Darcy의 법칙 및 Ghyben-Herzberg 식에 근거한 이론적인 연구들이 주로 이루어졌고, 근래에 현장관측이나 수리모형실험이 국내 외적으로 수행되고 있으나, 모든 영역의 지하수의 특성을 조사하는 것이 사실상 불가능하다. 이에 최근에는 컴퓨터 성능의 비약적인 발전과 더불어 다양한 수치해석방법에 의한 수치모델들이 개발되어 시뮬레이션에 적용되고 있다. 하지만 거의 대부분의 수치모델은 해안 대수층 수리특성을 투수계수에 의존하고 있을 뿐, 대수층 내부의 해수-담수에 의한 밀도류의 유동특성을 전혀 고려하지 못한 채 정수압에 근거한 해수-담수 경계면에 대해 모의하고 있는 정도이다. 따라서 본 연구에서는 해안 대수층 내부의 유동현상을 투수계수에 의존하는 방법에서 탈피하여 대수층 매체의 입경, 공극, 형상 등을 고려할 수 있을 뿐만 아니라, 염분 및 온도차에 의한 밀도류를 해석할 수 있는 강비선형 수치모델을 개발하여 해수침투 현상을 직접 모의한다. 나아가 대부분의 이전 연구들에서 간과하고 있는 해안지역의 대표적 물리력인 파랑과 조석의 영향이 해안 대수층의 해수침투에 미치는 영향, 해안 대수층의 지하수위 및 해수면의 수위차에 의한 해수침투 특성 그리고 이를 제어 할 수 있는 새로운 대응기술을 제안하는 것을 목적으로 한다.

• Solar Energy
• /
• v.17 no.3
• /
• pp.13-21
• /
• 1997

In this study, the theoretical equation of thermal storage efficiency was established to applied long term hot water storage system. The, effective thermal diffusivity and storage efficiency were, measured through the experiment to predict the degree of mixture in thermal storage tank. The effective thermal diffusivity was inversely preportional to the storage efficiency. The most effective storage efficiency was obtained under condition of low flow rate and using the perforated distributor.

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

본 논문에서는 해석영역의 형상변화, 즉 로켓노즐 내열재의 삭마로 인해 야기된 경계 면의 이동을 고려하면서 2차원 비정상 비등방성 재료의 열전달 문제를 해석할 수 있는 수치해석에 대해 기술하였다. 수치해석 알고리즘은 유한요소법이며 열해석시 경계면 이동으로 인한 격자계의 절점 좌표점이 계산과정 동안에 이동하는 변형 가능한 유한요소격자(transformable finite-element grid)를 사용하였다. 본 수치해석기법의 타당성 입증을 위해 극심한 열하중이 부여된 조건하에서 엄밀해가 존재하는 비정상 축대칭문제 및 고체로켓 노즐내열재에 대해 열해석을 수행하였으며, 그 결과 수치해는 엄밀해 또는 실험치와 잘 일치함을 보이었다. 여러 가지 복합재로 구성된 내삽노즐 또는 외삽노즐에 대해서도 안정된 수치해를 얻을 수 있었다. 아울러 노즐목삽입재로 탄소-탄소 복합재를 적용한 고체추진기관 내삽노즐을 해석모델로 택하여 열전달 해석을 수행하고 해석결과를 분석하였다. 노즐의 표면산화반응에 대한 열반응상수, 즉 Arrhenius 형태로 표시된 식에서 pre-exponential factor 및 activation energy 변화가 탄소-탄소 복합재 및 탄소-페놀릭 복합재의 삭마량에 미치는 영향에 대해서 고찰하였다.

• Proceedings of the KSEG Conference
• /
• 2002.04a
• /
• pp.193-200
• /
• 2002

울산광역시 울주군 삼남면 가천리와 상천리 일원에 발달하는 양산단층대 중남부의 발달특성과 제4기단층을 기재한다. 이 지역에는 중생대 백악기의 퇴적암과 화강암의 경계부 근처에 양산단층대의 주단층대와 부단층대로 확인되는 대규모의 단층파쇄대가 북북동-서남서 내지 남-북의 주향에 거의 수직으로 발달하고 있다. 단층조선은 수평에 가까우며 단층대내의 구조에 의하면, 주로 우향의 주향이동운동이 우세하다. 한편, 상천리와 가천리에는 이들 기반암과 제4기의 하성 사력층의 경계부 부근에서 제4기단층이 2조 발달하고 있다. 이들은 가천 제1단층과 가천 제2단층으로 기존의 양산단층대 일부가 제4기에 재활동한 것으로, 북북동 방향의 주향에 동측으로 고각의 경사를 보인다 단층조선은 거의 수평이며, 제4기 역들이 단층끌림에 의해 배열된 상태나 단층엽리내의 구조에 의하면 우향의 주향이동성운동이 우세하다.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.1
• /
• pp.171-181
• /
• 2013

The radial mapping rule in bounding surface model was experimentally investigated by analyzing the drained stress probe tests on Chicago clays. The experimental data obtained from 10 drained stress probe paths were analyzed to calculate the directions of the plastic strain increments. The anisotropic bounding surface model was adopted to represent a bounding yield surface which resides in the pre-consolidation yield stress of undisturbed clays. The projection origins were estimated by finding the interceptions of the straight lines passing through the current stress point and the imaginary yield stress point on the bounding surface. The results show that the projection origin is not fixed at a point but moves toward the direction of the stress probe path after it is established around the initial stress point.

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

담수가 해수에서 흘러드는 하구에서는 성층이 관측되며 이것은 난류의 미세구조를 변화시키는 주요 원인으로 작용한다. 이러한 성층화 현상은 하구 내 부유사의 군집인 하구 최대혼탁수(Estuarine Turbidity Maximum, ETM)의 형성에 영향을 주게 된다. 본 연구는 성층의 하구 최대 혼탁수 생성 메커니즘에 관심을 두고 수치모델링을 활용한 미세 난류의 부유사 거동 분석에 초점을 두었다. 성층과 전단응력 사이의 난류 혼합을 대표하는 유동인 켈빈-헬름홀츠 불안정성(Kelvin-Helmholtz Instability)을 도입하고 성층 경계면 근처에서 부유사의 이송을 높은 레이놀즈수(Reynolds number) 유동에서 RANS(Reynolds-averaged Navier-Stokes Simulation)보다 다양한 규모의 에너지 획득이 가능하여 미세 난류 구조 재현에 장점을 갖는 Large-eddy Simulation(LES)를 활용하여 모사하였다. 여기에서, 부유사는 주위 유동의 물리적 특성 변화에 영향을 미치지 않는 Passive scalar로 가정하였으며 $6^{th}$-order Lagrangian 다항식 보간법을 적용하여 입자의 이동 속도를 계산하고 이를 시간에 대해 적분함으로써 이동 궤적을 추적하였다. 수치 모델 결과 Lock-exchange 유동 내에서 켈빈-헬름홀츠 불안정성이 발생함에 따라 경계면 주위에 위치한 부유사가 billow 내에서 트랩핑(trapping)되는 것을 보여주어 KH-billow 혹은 braids 내의 미세 난류에 의한 영향이 확인되었다. 본 연구에서는 LES를 활용하여 성층류 및 성층류 내의 부유사 혼합을 모사하여 난류의 정도에 따른 이동 궤적의 차이에 대해서 분석함으로써 성층의 난류 강도 저하에 따른 부유사의 군집으로의 영향에 대해 서술한다.

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

자연하천의 주요한 특징 중 하나인 하천의 사행은 직선수로에서 예측되는 유속분포를 왜곡시키며 매우 복잡한 흐름구조를 형성한다. 이는 하상 경계면에서 발생하는 전단응력 분포의 변화를 야기하는데 하상 경계면에서의 전단응력은 다양한 경험적 관계에 의존하는 유사이동의 한계 소류력 산정 및 오염물질 거동해석의 분산계수 산정에 많은 영향을 미치게 된다. 물리적인 관측을 통한 하상 경계면에서의 전단응력의 관측은 다소 제한적이며 많은 비용을 요구한다. 따라서 하상 경계면에서 발생하는 전단응력의 경우 수심의 20% 이하의 연직 유속분포를 벽법칙에 적용하여 추정하는 방법이 주로 이루어지고 있다. 벽법칙을 이용한 하상 경계면의 전단응력을 계산하는 경우 대수중복층의 유속 분포 $u/u^*=(1/{\kappa})ln(zu^*/{\nu})+B$에서 무차원상수 ${\kappa}$와 B의 적절한 추정이 요구되어 진다. 일반적으로 무차원상수 ${\kappa}$와 B는 수리학적으로 매끄러운 벽면에서 대략 ${\kappa}=0.4$, B=5.5로서 경험적으로 이용되고 있다. 본 연구에서는 직선수로 및 다양한 사행수로의 3차원 흐름장 모의를 수행하여 벽법칙의 대수 중복층을 따르는 주흐름 방향 유속의 연직분포를 비교하였다. 수치모의 소프트웨어로서 Linux 기반의 OpenFOAM이 사용되었으며 모델의 검증을 위해 Chang(1971)에 의해 수행 된 사행수로에서의 유속장 관측 결과와 비교하였고 수치모의 결과가 실험 관측치와 잘 일치하는 것으로 판단되었다. 수치모의에 적용 된 사행수로의 형상은 Hey(1976)에 의해 제안 된 사행하천의 지형학적 인자들 간에 관계를 이용하여 사행도 1.03에서 2.42까지 총 7개의 사행수로 지형을 생성하였다. 사행도의 변화에 따라 만곡부 정점에서 대수중복층 구간의 주흐름 방향 유속의 연직분포를 비교한 결과, 본 연구에서 생성 된 모든 사행수로에서 대수중복층 구간의 무차원 유속 $u^+$와 무차원 거리 $z^+$가 로그 분포를 따르는 것으로 나타났으나 경험적으로 사용되었던 무차원상수 B의 경우 사행도가 증가 할수록 대수적으로 감소하는 경향이 나타났다. 본 연구에서는 이러한 관계가 무차원 상수 B값에 미치는 영향을 반영하여 수리학적으로 매끄러운 벽면에서 적용이 가능한 수정된 대수중복층 식을 제시하고자 한다.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.18-27
• /
• 2002

In order to improve the accuracy of the solution near the boundary in an analysis of viscous flow around an arbitrary boundary which move and be deformed using an Eulerian concept, a level-set based grid deformation method is introduced to concentrate grid points near the boundary. This paper presents a new monitor function which can easily control the level of the concentration of grid points along the boundary. Computations for steady flow around a semi-circular cylinder mounted on the bottom of the flow domain were carried out to check the improvement of the solution using the adaptive grid system with an immersed boundary method. The present numerical results show a good agreement with the solutions obtained by a body fitted grid system and more accurate solutions than those computed with non-adaptive grid system. For the validation of mechanical usefulness of the present method, an expanded analysis of flow around multi-body fixed in the flow domain was carried out. Finally, the present moving adaptive grid method was applied to a two-dimensional bubble rise problem. The computed results show well adapted grid points around the boundary of the bubble at every time and a good agreement with the result calculated by fixed grid system.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.10
• /
• pp.26-34
• /
• 2005

In this paper, a novel moving least squares interface welding method is proposed to carry out the coupled analysis of whole model composed of independently modeled finite element substructures with nodal mismatching interfaces. To verify the validity, and efficiency of the proposed interface welding method, various numerical examples are worked out including patch tests, convergence tests, and examples of coupled analyses of the structural systems with mismatching substructures. From the numerical tests, it is confirmed that one can efficiently carry out the coupled analysis of whole model composed of mismatching finite element substructures through the proposed method without any remeshing or any additional unknown.