• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2020년도 학술발표회
• /
• pp.88-88
• /
• 2020

수심평균 2차원 혼합모형은 하천환경에서 다양한 용존성 오염물질의 혼합현상을 모의하기 위해 널리 활용되어왔다. 2차원 혼합모형에서 분산계수는 하천의 전단 흐름에 의해 야기되는 오염물질의 퍼짐 현상을 표현하는 중요한 인자로서 작용하기 때문에 정교한 오염물질 혼합거동을 모의하기 위해서는 적합한 분산계수를 산정하는 것이 필수적이다. 분산계수를 실험적으로 산정하는 방법으로는 크게 모멘트법과 추적법으로 나뉘며, 비정상상태의 혼합거동에 대해 종방향 및 횡방향 분산계수를 동시에 산정할 수 있는 방법은 추적법 계열의 2차원 유관추적법(2D STRP)이 유일하다. 본 연구에서는 하천에 유입된 오염물질의 2차원 혼합해석을 위한 수치모형을 개발하였으며, 개발된 모형의 수치해를 바탕으로 다양한 Peclet 수의 범위에 대해 기존연구에서 제시된 2D STRP의 적용범위 및 성능을 정량적으로 분석하였다. 분석된 정보를 바탕으로 기존 2D STRP의 한계를 극복하기 위한 개선된 2차원 유관추적법(2D STRP-i)을 개발하고, 사행하천을 모형화한 실규모 하천실험시설에서 검증하였다. 기존 2D STRP의 성능평가 결과, Peclet 수가 낮은 조건일수록 농도분포의 예측 정확도가 감소하는 경향을 보였으며, 하안 경계에 도달하는 농도가 증가할수록 부정확한 결과를 초래하는 것으로 나타났다. 본 연구에서는 기존 2D STRP의 한계를 보완하여 더욱 정확한 분산계수를 산정하고자 하안 경계면 조건을 고려한 2차원 유관추적법(2D STRP-i)을 개발하였다. 2D STRP-i는 직교-곡선좌표계 기반의 2차원 이송-분산 방정식을 바탕으로 횡방향 유속분포 및 하안 경계조건을 고려할 수 있도록 개선되었다. 2D STRP-i는 공간적으로 상이한 이송효과 및 하안경계 조건을 적절히 반영함으로써 농도분포의 예측 정확도를 개선 시키는 것으로 평가되었으며, 하안경계면에서 농도가 증가하는 구간에서 기존 2D STRP의 결과와 비교하여 더욱 정확한 농도분포 및 분산계수를 제공하는 것으로 밝혀졌다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2011년도 학술발표회
• /
• pp.238-238
• /
• 2011

복단면 하도는 홍수소통을 위한 상부 하도와 생태계 서식처 강화와 유사이송능력 개선, 경제적인 유지유량의 확보 등을 도모하는 자연적인 저수로 하도로 구성된다. 그러므로 갈수기와 홍수기의 유량차가 큰 하천에서, 복단면 하도는 하천 공학, 환경 생태 관점에서 다양한 이점을 제공한다. 그러나 복단면 하도에서의 흐름저항, 조도계수, 통수능, 수위-유량관계, 횡방향 유속분포, 그리고 하상전단응력 분포는 단단면 하도와 상이하며, 중요한 차이점은 주수로부와 홍수터부 사이의 유속차로 인하여 발생하는 활발한 운동량 교환에 의해 추가적인 전단층이 생성되는 점이다. 주수로와 홍수터 사이에서 발생된 경계전단력(Apparent Shear Force, ASF)은 하천의 전반적인 통수능을 감소시켜 홍수의 하도내 저류효과를 증가시킬 수도 있다. 또한 주수로 및 홍수터가 분담할 수 있는 유량은 경계전단력에 민감하기 때문에, 홍수터에 수목 식재, 구조물 설치 등을 계획할 경우 정량적인 항력저항의 산정을 위하여 매우 중요하다. 현재까지 복단면에서 발생하는 경계전단응력을 추정하는 다양한 방법들이 개발되었으며, 각각의 방법으로 부터 경계전단력을 정량적으로 예측할 수 있다. 그러나 대부분의 방법이 실험에 근거한 경험적 방법이므로 각 식의 개발에 사용된 자료에만 적합할 수 있는 한계가 있으며, 균일한 수로의 기하학적 변수(전단면과 주수로의 하폭비, 홍수터와 주수로의 수심비, 홍수터와 주수로의 조도비, 주수로의 하폭대 수심비 등)에 의한 식으로 구성되어 수치모형에서 직접 활용하기에는 한계가 있다. 따라서 본 연구에서는 복단면 경계전단력 산정을 위하여 개발된 다양한 연구결과들을 비교 검토하고, 기존문헌 또는 웹상에서 가용한 수리실험 자료들을 이용하여 각 방법에서 계산된 유량과 실측 유량을 비교함으로써 각 방법이 지닌 한계를 분석하였다. 본 연구에서 검토된 5가지 방법은 Knight and Demetriou (1983)(이하 KD), Wormleaton and Merret (1990)(이하 WM), Cristodoulou (1992), Bousmar and Zech(1999) (이하 BZ) 그리고 Moreta and Martin-Vide (2010)(이하 MM)이다. BZ 방법을 제외하고 나머지 방법들은 전체 유량의 비교에서 5% 이내의 오차를 나타내었다. 그러나 나머지 4가지 방법 중 주수로부 유속의 비교에 있어서는 소규모 수로 실험자료만 이용한 KD 방법이 5% 이상의 오차를 나타내어 정확도가 떨어지는 것으로 나타났다. 따라서 비교적 단순한 형태의 Cristodoulou (1992) 방법이 적용하고자 하는 수로의 규모와 무관하게 복단면 유량예측에 적합할 것으로 판단된다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2019년도 학술발표회
• /
• pp.93-93
• /
• 2019

식생이 있는 개수로 흐름의 분석이 중요한 이유는 식생이 흐름 특성과 유사이송 및 하상변동에 영향을 주어 복잡한 흐름 현상을 만들어 내기 때문이다. 식생 밀도, 식생 길이, 식생이 식재된 형태에 따라 흐름의 현상이 달라지며, 흐름 특성의 변화로 인해 하상 및 주변 하천구조물에 영향을 줄 수 있다. 하천의 식생은 대부분 폭과 길이가 제한된 모습으로 존재하는 경우가 많으며, 식생이 존재하는 영역에서는 식생항력에 의해 유속이 감소하고 식생이 없는 영역에서 유속이 증가한다. 식생항력에 의해 변화된 흐름 특성은 일정 거리를 지나면 평형상태에 도달하게 되며, 식생의 경계면에서는 전단층이 발달하여 모멘텀의 교환이 활발하게 발생하는 현상을 보인다. 식생 흐름의 특성을 분석하기 위해 수치모형으로는 식생 흐름의 강한 비등방성을 모의할 수 있는 비등방성 수치모형이 적합하다. 본 연구에서는 부분 식생이 식재된 개수로에서 RSM (Reynolds stress model)을 이용하여 식생과 식생 주변의 흐름 특성 및 난류량을 분석하였다. 흐름의 변화에 따라 크게 3 구간으로 나누어 흐름 특성을 파악할 수 있으며, 식생의 밀도와 식생 영역의 폭의 길이에 따라 각 구간의 길이가 달라지는 것을 확인하였다.

• 자동차안전학회지
• /
• 제4권1호
• /
• pp.23-26
• /
• 2012

본 논문은 듀얼클러치 변속기 변속 메커니즘 개발을 위한 동역학적 시뮬레이션 해석 및 피로수명 해석에 관한 연구이다. 캠을 이용한 쉬프트 포크의 이송으로 기어를 변속하는 개념을 갖는 쉬프트 메커니즘의 작동성 및 설계안 검증을 위해 캠의 회전 속도, 접촉면의 경계조건 및 싱크로나이저 슬리브에서의 저항을 반영한 모델링을 구성하여 동역학적 시뮬레이션을 실시하였다. 각 파트의 응력 및 변형률 발생을 해석하여 설계 개선안을 도출하고, 개선안이 반영된 최종 시작모델의 시뮬레이션 상 최대 응력 발생 시점에서 시프트 포크의 각 접촉면 별 외력 작용을 분석하여 이를 적용한 피로 해석을 진행하였다.

• 한국지반환경공학회 논문집
• /
• 제13권10호
• /
• pp.55-62
• /
• 2012

본 연구에서는 개별요소해석 프로그램인 $PFC^{3D}$를 이용하여 대입경 조립재료의 특정 입도분포를 구현하는 알고리즘을 개발하였다. 기존의 입자 형상이나 입도를 구현하기 위해 사용되었던 clump logic 또는 cluster logic을 사용하지 않으며 요소의 경계면 이탈 현상과 경계면 파괴 등을 방지할 수 있는 초기 개별요소 모델링 기법을 고안하였다. 최종적으로 대입경 조립재료에 대한 대형 삼축압축시험을 수치 모델링하고 실내시험 결과와 비교하였다. 해석 결과, 실제 시료의 입도분포와 매우 흡사한 분포의 개별요소를 생성할 수 있었고, 적정 미시물성치 산정 과정(calibration)을 통해 다양한 구속응력 조건에 대한 대입경 조립재료의 특정 입도분포하에서의 전체적인 수치 모델링이 가능하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2021년도 학술발표회
• /
• pp.14-14
• /
• 2021

개수로에서는 반드시 자유수면이 있으며, 따라서 물과 공기와의 마찰은 관수로에 비해 상대적으로 매우 작고, 개수로의 전단응력 분포는 관수로와 달리 근본적으로 비대칭이다. 따라서 전단응력은 수로 바닥이나 측면에서만 작용하게 된다. 이러한 평균 전단응력 개념은 흐름에 의해 경계면 구성재료가 이동하는 이동상 수리학에서 유사이송 능력을 해석하는 기준이 되며, 경계면의 전단응력은 힘으로 표시하여 통상 소류력이라 한다. 이러한 복잡한 유체거동은 하천시설물 설계, 시공 및 관리에 있어서 구성재료의 보호능력에 따라 예상하지 못한 조건에서 쉽게 파손될 수 있다. 국내 하천의 경우 한계유속과 한계소류력에 의해 하천설계에 적용되고 있다. 한계 유속의 경우 간단한 수식에 의해 산정될 수 있지만 실제 하천의 보호능력을 대표하기는 힘들기 때문에 한계소류력이 동시에 고려되어야 한다. 한계소류력은 개수로 흐름에서 복잡하게 발생하는 이차류나, 난류 특성에 의해 산정하거나 예측하기는 매우 어렵다. 한계 소류력 뿐만 아니라 하천을 구성하는 재질의 조도계수 역시 균일하지 못하고 매우 예측하기 어렵다. 따라서 본 연구에서는 이러한 복잡한 양상을 나타내는 수리학적 요소에 대해 표준화된 실험수로에서 실험을 통해 평가하고, 체계화된 설계 지침이 되고자 연구를 진행하였다. 본 연구에서는 자연하천과 유사한 조건의 경사를 가지는 경사수로와 경사의 영향에서 자유롭게 평가를 진행하고자 기존 연구를 바탕으로 제작된 소류력 측정장치를 이용하였다. 하천의 설계나 평가에 적용되는 평균 소류력 개념은 복잡한 난류흐름에서 평가지표로써 대표하기 힘들기 때문에 유사 하천환경의 바닥에서 발생하는 소류력을 직접 측정하였다. 본 연구에 사용된 장치는 난류유속 u', v'을 이용하여 Reynolds stress산정하여 Total shear stress를 추정하는 기법을 사용하여 검증하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2004년도 학술발표회
• /
• pp.537-542
• /
• 2004

오염물질의 이송${\cdot}$확산 과정을 해석하기 위하여 계산효율이 높은 3차원 퍼프모형을 개발하였다. 본 연구에서 개발된 퍼프모형은 추적방법에 따라 전방추적모형과 후방추적모형으로 나눌 수 있으며, 이 두 가지 추적방법은 계산효율과 수치오차에 있어서 상이한 특성을 나타내었다. 전방추적 퍼프모형은 일정한 시간간격을 가지므로 정상상태의 연속오염원의 경우에 각각의 퍼프들이 동일한 질량을 갖는다. 그러므로 전방추적 퍼프 모형은 각 퍼프들간의 중첩정도가 일정하지 않다. 이에 관한 오차분석을 수행하기 위하여 본 연구에서는 무차원 퍼프중첩계수를 정의하였다. 퍼프중첩계수란 퍼프의 크기에 대하여 퍼프중심간의 거리가 떨어진 정도를 나타내는 무차원수로서 너무 작은 경우에는 정확도가 떨어지고 너무 큰 경우에는 계산효율이 감소한다. 전방 추적 퍼프모형의 경우, 중첩계수가 작은 초기구간에는 정확도가 떨어지며, 시간이 지남에 따라 중첩계수가 필요이상으로 증가하여 계산효율이 떨어지는 것으로 나타났다. 이에 비해 일정한 중첩정도를 갖는 후방추적 퍼프모형의 경우에는 전 영역에 걸쳐서 정확도와 계산효율이 높은 것으로 나타났다. 하지만 일정한 시간간격 마다 농도장을 계산하고자 할 때, 전방추적법은 단 한번의 전체계산을 통하여 수행가능하지만 후방추적법의 경우에는 매 출력시간마다 초기시점까지 반복해서 계산해야하는 단점이 있다. 경계처리에 있어서 입자추적모형과 상이한 방법을 사용하는 퍼프모형은 폐경계에서는 입자추적모형과 동일한 결과를 나타내지만 개경계에서는 다른 결과를 나타내었다. 또한 오염원이 임의의 공간적 분포를 갖는 경우, 퍼프모형은 입자추적모형보다는 적은 수의 퍼프를 사용할 수 있지만 이에 따른 경계면에서의 수치오차를 발생하였다. 본 연구에서는 개발된 퍼프모형을 검증하고 장${\cdot}$단점을 분석하기 위하여 흐름이 일정한 경우와 전단흐름의 경우에 대하여 이송${\cdot}$확산 수치모의를 수행하였으면, 이를 각각의 경우의 해석해 결과와 비교${\cdot}$분석하였다. 후방추적 퍼프모형은 전방추적 퍼프모형에 비하여 사용된 퍼프수와 관계없이 작은 오차를 발생하였으며, 전체적으로 퍼프 모형이 입자모형보다는 훨씬 적은 수의 계산을 통해서도 작은 오차를 나타낼 수 있다는 것을 알 수 있었다. 그러나 Gaussian 분포를 갖는 퍼프모형은 전단흐름에서의 긴 유선형 농도분포를 모의할 수 없었고, 이에 관한 오차는 전단계수가 증가함에 따라 비선형적으로 증가하였다. 향후, 보다 다양한 흐름영역에서 장${\cdot}$단점 분석 및 오차해석을 수행한 후에 각각의 Lagrangian 모형의 장점만을 갖는 모형결합 방법을 제시할 수 있을 것으로 판단된다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2012년도 제38회 춘계학술대회논문집
• /
• pp.304-308
• /
• 2012

본 연구는 에틸렌-공기 혼합물로 채워져 있는 변형 가능한 구리 관에서의 초음속 화염 전파를 수치적으로 살펴보았다. 탄화수소의 화염 전파를 해석하기 위하여 지배방정식으로 Navier-Stokes 방정식과 Arrhenius 형태의 1단계 화학 반응식을 활용하였으며 변형 가능한 관을 해석하기 위하여 Inviscid Euler 방정식을 활용하였다. 또한, 두 물질 간 경계면 추적을 위하여 Level-set 기법을, 경계값 결정을 위하여 ghost fluid 기법을 사용하였다. 이러한 수치적 기법을 바탕으로 관의 변화에 따른 초음속 화염 내 팽창파의 전파 및 그에 따른 간섭 현상을 밀도 및 속도 변화를 통해 확인하였으며 초음속 화염 전파에도 안전성이 확보되는 최소 관 두께를 예측할 수 있는 수치적 기반을 마련하였다.

• 한국수자원학회논문집
• /
• 제42권10호
• /
• pp.891-896
• /
• 2009

준설 시 발생하는 부유물의 이송 확산 과정을 해석하기 위하여 퍼프모형을 개발하였다. 본 연구에서 개발된 퍼프 모형은 추적방법에 따라 전방추적모형과 후방추적모형으로 나눌 수 있으며, 이 두 가지 추적방법은 계산효율과 수치 오차에 있어서 상이한 특성을 나타내었다. 경계처리에 있어서 입자추적모형과 상이한 방법을 사용하는 퍼프모형은 폐경계에서는 입자추적모형과 동일한 결과를 나타내지만 개경계에서는 다른 결과를 나타내었다. 또한 오염원이 임의의 공간적 분포를 갖는 경우, 퍼프모형은 입자추적모형보다는 적은 수의 퍼프를 사용할 수 있지만 이에 따른 경계면에서의 수치오차를 발생하였다. 흐름이 일정한 경우와 전단흐름의 경우에 대하여 이송 확산 수치모의를 수행하였으며, 이를 각각의 경우의 해석해 결과와 비교 분석하였다. 후방추적 퍼프모형은 전방추적 퍼프모형에 비하여 사용된 퍼프수와 관계없이 작은 오차를 발생하였으며, 전체적으로 퍼프모형이 입자모형보다는 훨씬 적은 수의 계산을 통해서도 작은 오차를 나타낼 수 있다는 것을 알 수 있었다. 그러나 Gaussian 분포를 갖는 퍼프모형은 전단흐름에서의 긴 유선형 농도분포를 모의할 수 없었고, 이에 관한 오차는 전단계수가 증가함에 따라 비선형적으로 증가하였다. 이와 같이 퍼프모형은 다양한 수환경에 적용할 경우, 뛰어난 효율성에 비해 정확도가 다소 감소하는 경향이 있지만, 입자추적모형과의 연계 모의 등을 통해 준설지점 인근의 근역에서의 오염원 해석에 사용될 수 있다.

• 대한조선학회논문집
• /
• 제47권3호
• /
• pp.279-290
• /
• 2010

In this paper, a VOF method called RHRIC (refined high resolution intertace capturing) is introduced for solving the motion of the free surface and applied to the simulation of the advection of rigid interiaces of different shapes and a 20 dam-break problem, which are typical benchmark test cases. The numerical results for the interface advection cases are compared to the analytic solutions, while the available experimental data and other numerical results of various free surface methods for the dam-break problem are provided for the validation of the proposed VOF method. The same simulations were also carried out using the original HRIC scheme and a modified HRIC scheme called MHRIC for comparison. Although the RHRIC uses a simple order scheme, a basis of the original HRIC scheme, lower than the third-order ULTIMATE-QUICKEST used by the MHRIC, it provides an improved accuracy over the two previous HRIC methods.