V. Problems of 5th and 11th Chapter (Viscous Flow and Pipe Flow)

5.1 Derive the equation of the velocity profile in the pipe for the case of laminar and turbulent flow using the contents of Brebbia textbook (층류인 경우 관에서의 유속분포와 수두손실에 대한 식(Poiseulli 식)을 유도하라. 난류인 경우 수두손실에 대한 식(Daycy-Weisbach식)을 유도하라.) 마찰력(전단력) = 압력에 의한 힘 = 수두손실(압력손실) 위의 식을 정리하면 다음과 같다. or

<관로 내부에서의 유속, 압력, 전단응력> <관류의 전단응력 및 유속 분포> <관로 내부에서의 유속, 압력, 전단응력> <관류의 전단응력 및 유속 분포>

뉴턴의 제1점성법칙은 다음과 같다. 위의 식을 에 대입하면 다음과 같다 뉴턴의 제1점성법칙은 다음과 같다. 위의 식을 에 대입하면 다음과 같다. 위의 유속에 대한 식을 관의 내경에 대하여 적분하면 다음과 같다. 위의 식의 적분상수 C는 다음과 같이 구할 수 있다. at , 최대유속은 r=0인 점에서 발생하며, 다음과 같다.

평균유속은 다음과 같이 유속의 유한단면적을 곱한 다음 전체를 적분하여 전체 단면적으로 나누어준다 평균유속은 다음과 같이 유속의 유한단면적을 곱한 다음 전체를 적분하여 전체 단면적으로 나누어준다. 다음과 같이 평균유속은 최대유속의 1/2이다. 마찰에 의한 에너지 손실 : 유속 - 압력의 구배도 – 마찰력 압력손실 - 마찰 - 에너지 손실 위식을 정리하면 다음과 같은 압력손실을 구할 수 있다.

관류의 동수반경은 다음과 같다. 동수반경 = 단면적 / 윤변 압력손실(수두손실)은 다음과 같다 관류의 동수반경은 다음과 같다. 동수반경 = 단면적 / 윤변 압력손실(수두손실)은 다음과 같다. 위의 식을 Poiseulli식이라 한다. 층류의 에너지 손실식은 점성계수, 유속, 관의 내경, 길이에 관계가 있다. 정수장에서 배수관망 설계, 상수관망 설계는 다음과 같다. 1. 수요 : 원하는 유량 : 20만 * 0.5톤/일/인 = 10만톤/일 필요한 공급 유량 Q 2. 배수관망의 조건에 따라 필요한 관의 P, L이 결정된다. 3. V, D는 위의 조건으로부터 구한다.

점성법칙이 유효하게 적용되기 위해서는, 작은 내경의 관로의 점성유체에 대한 층류이어야 한다 * 점성법칙이 유효하게 적용되기 위해서는, 작은 내경의 관로의 점성유체에 대한 층류이어야 한다. 난류 <층류와 난류의 유속 분포도>

층류의 경우에는 난류인 경우에는 실험 결과로부터, 여기서, 층류 : Re < 2000 난류 : Re > 4000 천이영역 : 2000<Re<4000

5. 2 Derive the equations of the head loss for the problem 5 5.2 Derive the equations of the head loss for the problem 5.1 such as Darcy-Weisbach equation. 난류 유동 : Darcy-Weisbach 공식 관류유동에 있어서 난류인 경우의 수두손실은 Darcy-Weisbach식으로 다음과 같다. 여기서, 는 마찰상수이다. 개수로 유동에서의 Darcy-Weisbach식은 다음과 같다. 관의 형상에 따른 국부 수두손실은 다음의 식으로 평가된다.

압력과 마찰력의 평형 관계로부터 다음의 식을 구성할 수 있다 압력과 마찰력의 평형 관계로부터 다음의 식을 구성할 수 있다. 마찰력(전단력) = 압력에 의한 힘 = 수두손실(압력손실) 위의 식을 정리하면 다음과 같다. or 난류에 대한 실험에 의하면 마찰력은 다음과 같다. 층류에 대한 뉴턴의 제1점성법칙은 다음과 같다. 따라서, 수두손실은 다음과 같다.

Osborne Reynolds는은 유동특성을 규명하기 위하여 다음과 같은 Reynolds 번호를 고안하였다 Osborne Reynolds는은 유동특성을 규명하기 위하여 다음과 같은 Reynolds 번호를 고안하였다. : 동점성계수 위의 식에서 알수 있드시 관내의 유체거동은 에 따라 결정된다. 유속과 수두손실과의 관계는 다음 그림과 같다. Figure 3.4 유속과 수두손실

층류 Re < 2000 Laminar Flow 천이영역 2000 < Re < 4000 Transition Region 난류 Re > 4000 Turbulent Flow 층류에 대해서는 다음과 같이 Poiseuille의 식을 사용한다. 마찰력 -> 압력감소 -> 수두손실 -> 에너지 손실 평형개념 : 마찰력-압력감소 -> 관류에서의 유속의 식 레이놀즈 번호 : 관성력 / 점성력 -> 층류, 난류, 천이영역

5.3 Derive the equations of the velocity from the head loss equation for the problem 5.2 such as Manning's, Chezy's, and Hagen-William's formula. 1) Manning 공식 for turbulent regions of rough pipes, or for canals. R = 동수반경 = 윤변/단면적 = perimeter/area n = 조도계수 = roughness coefficient : hydraulic radius ( for a circular pipe)

2) Hazen-Williams 공식 천이영역에 대해서는 다음과 같이 Hazen-Williams 공식을 적용한다 2) Hazen-Williams 공식 천이영역에 대해서는 다음과 같이 Hazen-Williams 공식을 적용한다. for SI system →

3) Chezy's 단면적이 일정한 유로를 비압축성유체가 정상, 등속의 난류로 흐를 때 벽에서의 전단응력은 평균속도의 제곱에 비례한다. (5.6.1) 여기서 는 무차원계수이다. 개수로와 비원형폐수로에서는 전단응력이 全表面(전표면)에 걸쳐 일정하지가 않다. 이런 경우에는 로서 벽에서의 평균전단응력을 사용한다. 非圓形流路(비원형유로)에서 나타나는 부차유동(secondary flow)은 壁剪斷應力(벽전단응력)을 균등화시키려는 작용을 한다.

그림 5. 15 水路(수로)의 검사체적에 작용하는 축방향 힘 그림 5 그림 5.15 水路(수로)의 검사체적에 작용하는 축방향 힘 그림 5.15는 개수로 또는 폐수로를 흐르는 정상, 둥속유동을 표시한 그림이다. 개수로에 대해서는 친과 f2가 같다. 그리고 유동이 일어나는 것은 퍼텐셜에너지의 감소 에 기인한다. 폐수로에서는 유동을 야기 시키는 에너지는 퍼텐셜에너지의 감소 이외에 압력강하 도 있다. 원관 내에서 연직하향으로 흐르는 경우 압력 는 유동에 따라 증가할 것이다.

그러나 壁面(벽면) 剪 斷應力(전단응력)을 극복하고 유동을 계속할 수 있는 에너지를 하려면 퍼텐설에너지의 감소 가 보다 커야 한다. 손실을 유용에너지의 감소와 결부시켜 에너지방정식인 식 (3.10.1)을 다음과 같이 나타낼 수 있다. 속도수두 은 동일하므로 다음과 같이 나타낼 수 있다. 손실 (5.6.2)

등속도유동이라고 가정하였으므로, 방향에 적용한 선형운동량 방정식인 식 (3. 11. 2)는 로 된다 등속도유동이라고 가정하였으므로, 방향에 적용한 선형운동량 방정식인 식 (3.11.2)는 로 된다. 여기서 P는 水路(수로)의 濕潤長(습윤장)(wetted perimeter)이다. 다시 말해서 유체와 접하고 있는 (자유액체표면을 제외) 벽 둘레 길이이다. L sin θ = z1 - z2이므로 (5.6.3) 이다. 식 (5.6.2), 식 (5.6.3)으로부터 식 (5.6.1)을 사용하면 손실 (5.6.4)

이 식에서 로 놓았다. R을 유로의 수력반지름(hydraulic radius)이라 한다. R은 문제를 다롤 때 대단히 유용하다 이 식에서 로 놓았다. R을 유로의 수력반지름(hydraulic radius)이라 한다. R은 문제를 다롤 때 대단히 유용하다. 원관의 경우 이다. 식 (5.6.4)에서 손실항의 단위는 m · N/N또는 ft ·lb/lb이구 마찰에 의한 水頭損失(수두손실)(head loss due to friction)이라는 의미에서 로 표시한다. 水路(수로)의 단위길이 당 유체의 단위중량당의 손실을 S 로 정의하면 (5.6.5) 이고, V에 관하여 풀면 (5.6.6) 로 된다.

계수또는 C는 실험에 의하여 결정된다. 이 식이 Chezy의 공식이다 계수또는 C는 실험에 의하여 결정된다. 이 식이 Chezy의 공식이다. Chezy계수 C는 원천적으로 어느 정도의 치수를 갖는 水路(수로)나 壁面狀態(벽면상태)에 대하여 상수로 간주하였다. 그러나 현재는 C에 관한 여러 종류의 공식들이 일반적으로 사용된다.

5. 4 Explain the concept of Reynold's number and hydraulic radius 5.4 Explain the concept of Reynold's number and hydraulic radius. 1)reynold's number 영국의 유체역학자 O.레이놀즈가 발견한 것으로 움직이는 유체 내에 물체를 놓거나 유체가 관속을 흐를 때 난류와 층류의 경계가 되는 값을 말한다. 물체가 정지유체 속을 진행하는 경우 그 물체의 길이 또는 관지름을 L, 진행속도 또는 평균속도를 U,ρ을 유체의 밀도라 할 때 R=ρUL/μ이며 R이 레이놀즈 수이다.

2)동수반경 흐르는 물의 단면을, 물길 옆으로 스며든 물의 길이로 나눈 값 2)동수반경 흐르는 물의 단면을, 물길 옆으로 스며든 물의 길이로 나눈 값. 동수반경 = 단면적 / 윤변 Reynolds수 는 점성력에 대한 관성력의 비이다. 임계 Reynolds수는 管內流動(관내유동), 경계층유동 또는 유체 속에 잠긴 물체 둘레의 유동이 층류유동이냐 혹은 난류유동 이냐와 같이, 유동형태를 구별하는 판정값이다. 압축성유동에서는 Reynolds수보다 Mach수가 일반적으로 더 중요한 의미를 갖는다.

Reynolds수 는 점성력에 대한 관성력의 비이다 Reynolds수 는 점성력에 대한 관성력의 비이다. 임계 Reynolds수는 管內流動(관내유동), 경계층유동 또는 유체 속에 잠긴 물체 둘레의 유동이 층류유동이냐 혹은 난류유동 이냐와 같이, 유동형태를 구별하는 판정 값이다. 압축성 유동에서는 Reynolds수보다 Mach수가 일반적으로 더 중요한 의미를 갖는다. 층류유동은 유체입자가 층 또는 막 안에서, 한 층이 인접한 층 위를 원활하게 미끄러지도록 운동하는 유동을 말한다. 그러므로 층류유동에서의 운동량수송은 층과 층 사이에서 분자적 운동량교환만으로 이루어진다. 불안정성, 즉 난류로 전향하려는 경향은 인접층간의 상대운동을 저지하려는 粘性剪斷力(점성전단력)에 의하여 억제된다. 그러나 난류유동은 격렬한 橫斷운동량수송(transverse interchange of momentum)을 가지면서 유체 입자가 매우 산만한 운동을 하는 유동이다. 유동의 성질, 다시 말해서 층류냐 혹은 난류냐, 그리고 층류화 경향에 미치는 난류의 상대적 중요성을 나타내는 척도는 Reynolds 수에 의하여 정량적으로 표시된다.

다음을 만족할 때 두 유동은 力學的(역학적) 相似(상사)가 이루어진다고 말한다. 1. 기하학적으로 相似(상사)하다 다음을 만족할 때 두 유동은 力學的(역학적) 相似(상사)가 이루어진다고 말한다. 1. 기하학적으로 相似(상사)하다. 즉 대응 線型次元比(선형차원비)가 일정하다. 2. 대응하는 유선(流線)들이 기하학적으로 상사(相似)하다. 또는 대응되는 모든 점에서 일정한 壓力比(압력비)를 갖는다. 2개의 기하학적 相似流動(상사유동)을 생각할 때, Reynolds는 유동을 기술하는 일반 미분방정식이 같으면 力學的(역학적) 相似(상사)도 성립한다는 것을 추론하여 냈다.

한 조의 방정식에 포함되어있는 질량, 길이 및 시간의 단위를 바꾸고, 이들 변형된 방정식이 원래의 방정식과 일치되도록 만족하여야 할 조건을 결정하면, 이 양쪽 경우에 대하여 무차원수 가 같은 값을 가져야 한다는 것을 Reynolds는 증명하였다. 여기서 는 특성속도, 은 특성길이, 는 질량밀도 그리고 는 점성계수이다. 이 무차원수를 Reynolds 수라 R로 쓴다. (5.1.1)

무차원수의 물리적 의미를 결정하기 위하여 Reynolds는 그림 5 무차원수의 물리적 의미를 결정하기 위하여 Reynolds는 그림 5.1에 도시한 바와 같이 유리관을 흐르는 물의 유동에 대한 실험을 하였다. 유리관을 수평으로 설치하고 한 쪽 끝은 물탱크에 연결하고, 다른 쪽 끝에 밸브를 결합하였다. 관의 上端(상단)에는 매끈한 나팔관을 부착하여 물의 교란을 막았고, 나팔관 前面(전면)의 한 점에서 가느다란 染料線(염료선)이 분사(噴射)되도록 染料噴出口(염료분출구)를 설치하였다. Reynolds는 특성속도로 평균속도 V , 특성길이로서 관의 지름 D를 택하여 로 하였다.

유량이 적을 때 염료선은 관속에서 직선을 그리며 흘러 유동이 층류임을 보여준다 유량이 적을 때 염료선은 관속에서 직선을 그리며 흘러 유동이 층류임을 보여준다. 유량이 증가함에 따라 Reynolds 수는 증가한다. 왜냐하면 는 상수이다. 는 유량에 비례하여 증가하기 때문이다. 유량이 계속 증가하면 염료의 사선(絲線)이 흔들거리다가 드디어 갑자기 파괴되고 관 전체로 확산되어 버린다. 이러게 되면 유동은 격렬한 운동량 교환을 일으키면서 난류로 변환되어 질서정연한 층류유동은 완전히 붕괴되어 버린다. 세밀한 실험을 통해 Reynolds는 부터 난류가 발생됨을 알아냈다. 그림 5.1 Reynolds 실험장치

그 후 한 연구자가 Reynolds의 실험 장치를 그대로 사용하여 실험 전 여러 날 동안 탱크의 물을 정지 상태로 유지하고, 물이나 실험장치의 진동을 세심하게 배제시키고 실험한 결과 의 값을 얻기도 하였다. 이들 값들 上臨界(상임계) Reynolds수(Reynolds upper critical number)라 한다. 그러나 통상의 배관은 불규칙하여 훨씬 작은 Reynolds수에서 난류를 발생시키므로 상임계 Reynolds수는 실용적 의미가 없다.

Reynolds는 난류상태로 흐르는 유리관 속의 물의 속도를 점차 줄여가면서 Reynolds수를 계산한 결과 R이 2,000 이하로 되는 속도에서는 항상 층류가 됨을 발견하였다. 이 값을 관 유동에 대한 하임계(下臨界) Reynolds수(Reynolds lower critical number)라 하고, 이 값은 매우 중요한 의미를 갖는다. 통상의 배관에서 층류유동이 난류로 변하는 Reynolds수는 2,000~4,000의 범위이다. 그러나 통상층류와 난류를 구분하는 임계값으로서 을 가정한다. 층류에서 손실은 속도에 비례하지만, 난류에서는 속도의 1.7~2.0 에 비례한다.

오늘날 사용되는 Reynolds 수는 원형단면의 直管(직관)에 대한 Reynolds 수 이외에도 여러 가지가 있다 오늘날 사용되는 Reynolds 수는 원형단면의 直管(직관)에 대한 Reynolds 수 이외에도 여러 가지가 있다. 예를 들면, 유체 속에서 운동하는 구(球)에 대해서는 가 사용된다. 여기서 는 구(球)의 속도,D 는 구(球)의 지름 그리고 와 는 유체밀도와 점성계수이다.

비압축성유체의 주어진 유동성질은 그의 Reynolds 수의 크기에 의하여 특정지어진다 비압축성유체의 주어진 유동성질은 그의 Reynolds 수의 크기에 의하여 특정지어진다. 이 크다는 것은 분자의 어느 한 항 또는 모든 항이 분모에 비하여 상대적으로 크다는 것이다. 이것은 유체의 공간이 넓은 경우, 고속(高速)유동의 경우, 밀도가 큰 경우, 점성계수가 극히 작은 경우 또는 이들 경우들이 조합되어 일어나는 경우에 해당한다. 분자는 관성력(inertial force), 다시 말해서 유체의 가속 또는 감속으로 나타나는 힘에 관련한다. 분모는 점성전단력(粘性剪斷力)의 원인이 되는 항이다. 따라서 Reynolds 수는 점성력에 대한 관성력의 비로 생각할 수 있다. 큰 Reynolds 수는 손실이 속도의 제곱에 비례하는 격렬한 난류를 의미한다.

난류는 점성작용에 의하여 기계적 에너지를 비가역량으로 급격히 전환시키는 극히 많은 작은 와(渦)로 구성된 미세규모(fine scale)일 수도 있고 거대한 와동(渦動)이나 강(江)에서의 소용돌이, 대기 중에서의 돌풍과 같이 큰 규모(large scale)일 수도 있다. 큰 와동은 보다 작은 와동들을 유발하고, 차례로 미세규모의 난류를 생성해낸다. 난류유동은 매끈한 유동, 경우에 따라서는 등속도 유동에 부차유동(副次流動)을 중첩시킨 유동으로 생각할 수 있다. 미세규모의 난류운동은 고진동수(高振動數)를 갖는 속도의 작은 변동을 갖는다. 이 변동량의 제곱평균제곱근과 변동의 부호가 바뀌는 주파수는 난류의 정도를 나타내는 정량적인 척도이다. 일반적으로 난류의 세기는 Reynolds 수가 증가함에 따라 커진다.

중간정도의R 값에서는 점성과 관성의 영향이 모두 중요하다. 그리고 점성계수가 변하면 속도분포와 유동저항도 변한다 중간정도의R 값에서는 점성과 관성의 영향이 모두 중요하다. 그리고 점성계수가 변하면 속도분포와 유동저항도 변한다. 같은 값의 R에 대하여 기하학적으로 상사(相似)한 閉流路系(폐유로계)는(예컨대 상사비가 1대 2인 경우) 속도수두에 대한 손실비가 같다. Reynolds 수는 한 유체에 대한 실험결과를 가지고 다른 종류의 유체가 흐르는 相似(상사)한 유동에 대하여 특성을 예측할 수 있는 수단이 된다.

hydraulic radius: 흐르는 물의 단면을, 물길 옆으로 스며든 물의 길이로 나눈 값 hydraulic radius: 흐르는 물의 단면을, 물길 옆으로 스며든 물의 길이로 나눈 값. 유로의 수력반지름(hydraulic radius)이라 한다. 문제를 다롤 때 대단히 유용하다. 원관의 경우 이다.

5.5 Explain the computational algorithm of the pipe network using the concept as shown in the section of 4.3-4.5 in the Brebbia's textbook. 상수관망내의 모든 요소와 절점이 정의되면, 각 요소에 대하여 수두손실과 유량과의 관계를 정의하여야 한다. 다음의 I 요소에 대한 관경 D, 조도계수 C, 길이 L을 나타내었다. 각 요소별 절점의 연결도와 유량의 정의

여기서, k, j는 절점번호, i는 요소번호를 나타낸다 여기서, k, j는 절점번호, i는 요소번호를 나타낸다. i 요소로 연결되는 각 절점 k, j의 수두차이는 , 요소내 유량은 이다. 수두손실로부터 유량을 평가하기 위하여 층류에 대한 Poiseuille의 공식을 사용하면 다음과 같다(마찰력과 압력). 여기서, 계수는 유체와 요소의 특성에 따라 결정된다. 위의 식은 선형방정식으로 간단하지만 정확하지는 않다.

천이영역 및 난류에 대해 Hazen-Williams 공식을 사용하면 다음과 같다 (비선형 문제)

요소행렬식(각 관의 절점별 유량과 수두와의 관계) 절점의 유량은 다음의 식과 같이 자기중심적으로 평가된다 * 요소행렬식(각 관의 절점별 유량과 수두와의 관계) 절점의 유량은 다음의 식과 같이 자기중심적으로 평가된다. 확산, 지하수 유동의 개념도 마찬가지이다. 위의 식을 수두차이로 표시하기 위하여, 절점별 유량은 다음과 같이 절점에서 요소로 유입되는 경우에는 “+”, 요소에서 절점으로 유출되는 경우에는 “-” 부호를 가진다(요소 중심적인 관점).

위의 식을 수두차이 에 대한 식으로 정리하면 다음과 같다. 행렬과 벸터로 표시하면 다음과 같다 위의 식을 수두차이 에 대한 식으로 정리하면 다음과 같다. 행렬과 벸터로 표시하면 다음과 같다. 여기서, 는 요소내 절점별 유속 벸터, 는 요소별 특성 행렬, 는 요소별 절점 수두 벸터이다.

전체 시스템 방정식 관망 해석은 전체 절점에 대한 유량 및 수두를 계산하는 것이 주 알고리즘이다 *전체 시스템 방정식 관망 해석은 전체 절점에 대한 유량 및 수두를 계산하는 것이 주 알고리즘이다. 따라서, 전체 절점에 대하여 각 요소별 행렬식을 합해 주어야 한다. 이때, 각 절점별로 연속방정식 즉, 유량의 물질수지가 유지되어야 한다. 따라서, 관망해석은 각 관로(요소)별로 에너지손실과 흐름을 해석하고(왜냐하면 관로 내부에서는 실제 흐름에 의해 마찰이 발생하고, 이러한 마찰과 흐름의 관계를 평형 상태에 도달했다고 가정), 전체 절점에 대한 물질수지를 고려하는 것(여러 관로에서 유입되는 유량이 평형상태에 도달함)이라 할 수 있다. 이러한 관로의 유동 해석은 상수, 하수, 기타 가스유동도 마찬가지이다. 각 절점의 물질수질를 고려한 각 절점별 연속방정식은 다음과 같다. 여기서, 는 j 절점에서의 소비나 공급유량이다.

다음과 같은 상수관망을 고려하여 위의 알고리즘을 나타내보자 다음과 같은 상수관망을 고려하여 위의 알고리즘을 나타내보자. 각 절점별 연속방정식을 평가하면 전체 시스템 방정식을 구성할 수 있다. 각 절점에 대한 연속방정식을 평가하기 위하여 5 절점에 대한 연속방정식은 다음과 같다. 요소 시작절점, 최종절점, 1 2 3 4 5 6 7 8 9

위식의 유량을 평가하기 위하여 5 절점에 연결된 요소를 다음과 같이 살펴보자 위식의 유량을 평가하기 위하여 5 절점에 연결된 요소를 다음과 같이 살펴보자. 각 요소의 유량은 다음과 같이 수두차이에 의해 결정된다. 따라서, 5 절점으로 유입되거나, 유출되는 유량은 다음과 같다.

위의 3개의 식을 합하면 다음과 같다. 위의 과정이 5 절점에 대한 각 요소별 절점 유량을 합해주는 과정으로 전체 절점에 대한 합산 과정은 전체 절점에 대하여 연결된 요소의 절점의 유량을 합쳐주는 과정을 의미한다. 관망이외에 일반적인 유한요소법의 수치해석 방법에 있어서, 전체 절점에 대한 행렬식을 구성하기 위해서는 각 요소별 절점에 대한 식을 합쳐주어야 한다. 일반적으로, 전체 요소를 결합하기 위하여 다음과 같은 체계적으로 합치는 방법을 제시하였다.

(i) 유한요소법과 마찬가지로 각 요소별 절점의 연결도표를 작성한다. (ii) 양의 유량은 시작절점에서 최종절점으로 흐른다 (i) 유한요소법과 마찬가지로 각 요소별 절점의 연결도표를 작성한다. (ii) 양의 유량은 시작절점에서 최종절점으로 흐른다. (iii) 각 요소에 대하여 계수를 계산한다. (iv) 해당되는 절점에 알맞게 합쳐주기 위하여 각 요소의 절점의 연결도를 고려하여 전체 요소를 합쳐준다. 요소 시작절점, 최종절점, 1 2 3 4 5 6 7 8 9

다음에 각 요소행렬이 합쳐지는 예를 나타내었다. (전체 절점에 대한 합산과정은 각 요소에 대해서 계산해주고 합쳐주는 과정이다 다음에 각 요소행렬이 합쳐지는 예를 나타내었다.(전체 절점에 대한 합산과정은 각 요소에 대해서 계산해주고 합쳐주는 과정이다. 즉, 전체 요소에 대한 요소행렬식을 평가하여 합쳐주면 계산이 완료됨.) 위의 요소행렬을 더하는 방법은 i 요소내 k, j 절점의 경우 다음과 같이 일반화할 수 있다. - (k,k)와 (j,j) 위치에 계수를 더해준다. - (k,j)와 (j,k) 위치에 계수를 더해준다.

위의 과정을 Fortran으로 프로그램밍하면 다음과 같다. c c. initialize global matrix 위의 과정을 Fortran으로 프로그램밍하면 다음과 같다. c c... initialize global matrix ... c npmax : total number of nodes = 6 c glob : global matrix after assembling c do i=1,npmax do j=1,npmax glob(i,j) = 0 enddo enddo c c... assemble each element matrix for overall element ...

c nelmax : total nuber of elements : 9 c do nel=1, nelmax c c c nelmax : total nuber of elements : 9 c do nel=1, nelmax c c... read information about element connectivity ... c nodemax : total number of node in each elment : 2 c c read(nodeele(nel,ie), ie=1, nodemax) c... c... compute each element of element matrix at each node ... element 1st node 2nd node 1 2 3

c. c. compute each element of element matrix at each node c... c... compute each element of element matrix at each node ... c ig : global nodal number of ieth node in element nel c jg : global nodal number of jeth node in element nel c do ie=1, nodemax do je=1, nodemax ig=nodeele(nel,ie) jg=nodeele(nel,je) c c... compute each element in k matrix ...

c ek= rk(nel). (-1). (ie+je) c c. assemble each element in k matrix c ek= rk(nel)*(-1)**(ie+je) c c... assemble each element in k matrix ... c glob(ig,jg) = glob(ig,jg)+ek c enddo enddo enddo 합쳐진 최종방정식은 다음과 같이 행렬과 벸터를 이용하여 간단히 나타낼 수 있다. 여기서, 는 관망의 특성 행렬, 는 관망의 수두 벸터, 는 관망 소비 벸터이다.

5.6 Determine the discharge through the system as shown below by using the concept of the equivalent pipe where, =0.5, =300m, =600mm, = 2mm, =240 m, =1m, =0.3mm, , and H=6m.

[예제 11. 4] 그림 11. 7에서 , , , , , , , 및 이다. 이 시스템을 흐르는 유량을 구하라 [예제 11.4] 그림 11.7에서 , , , , , , , 및 이다. 이 시스템을 흐르는 유량을 구하라. <풀이> 에너지방정식으로부터 정리하면

= 0.0033, = 0.0003의 값과 그림 5.21로부터 의 값을 완전난류 범위에 대해서 가정하면, 이 값들로 을 구함으로써 = 2.848 m/s, = 1.025 m/s 그림 5.21에서 = 0.0265, = 0.0168, 에 대해 다시 풀면 = 2.819m/s 그리고 이다.

5.7 Determine flow through each pipe and the pressure at B as shown below for a total flow of 12 cfs. Where, =3000ft, =1ft, =0.00 ft; =2000ft, =8in, =0.0001ft; =4000ft, =16in, =0.0008ft; =2.00 slugs/, =0.00003/s, =8psi, =100ft, =80ft. [예제 11.6] 그림 11.9에서 = 3,000ft, = 1ft, = 0.001ft; = 2,000ft, = 8in, = 0.0001ft; = 4,000ft, = 16in, = 0.0008ft; = 2.0slug/ft3 , = 0.00003 ft2 /s, = 80psi, = 100ft, = 80ft. 총 유량 12cfs에 대해 각 파이프를 통과하는 유량과 B에서의 압력을 구하라.

유동이 파이프들 사이로 나뉘어졌다가 다시 합쳐지도록, 2개 또는 그 이상의 파이프로 연결된 것을 병렬파이프(parallel-pipe)시스템이라 한다. 직렬파이프 시스템에서는 같은 유체유동이 모든 파이프를 통과하면서 수두손실이 합산되나, 병렬파이프 시스템에서의 수두손실은 어느 파이프에서나 같으며 유량이 합산된다. 병렬파이프 시스템을 해석하기 위해서, 국부손실은 등가길이로서 각 파이프의 길이에 더해진다고 가정한다. 그림 11.9로부터 만족되어야 하는 조건은 (11.5.1) 이다. 여기서 , 는 점 와 의 높이이고, 는 진입파이프나 출구파이프를 통과하는 송출량이다.

이때 두 가지 유형의 문제가 발생한다. (1) A와 B 에서 수력구배선의 높이를 알고 유량Q 를 구하는 문제 이때 두 가지 유형의 문제가 발생한다. (1) A와 B 에서 수력구배선의 높이를 알고 유량Q 를 구하는 문제. (2)Q 를 알고 수두손실과 유동의 분포량을 구하는 문제, 파이프의 크기, 유체의 성질 및 조도는 알려져 있다고 가정한다 첫 번째 유형은, 수두손실이 수력구배선의 강하에 해당하기 때문에 유량에 대한 단순한 파이프문제의 해가 된다. 이 유량들을 더하면 총 유량이 된다. 두 번째 유형의 문제는 어느 한 파이프의 수두손실 뿐 아니라 유량도 알려져 있지 않기 때문에 더 복잡하다. 따라서 아래와 같이 풀어나가는 것이 좋다.

1. 파이프 1을 통한 유량 를 가정한다. 2. 가정된 유량을 사용하여 를 구한다. 3. 를 이용하여 , 를 구한다. 4 1. 파이프 1을 통한 유량 를 가정한다. 2. 가정된 유량을 사용하여 를 구한다. 3. 를 이용하여 , 를 구한다. 4. 공통의 수두손실에 대한 3개의 유량을 가지고, 문제에서 주어진 Q가 , , 의 같은 비율로 다음과 같이 나뉘어진다고 가정한다. 5. 계산된 , , 에 대해 , , 를 계산하여 이 유량의 타당성을 검토한다

이 과정은 많은 수의 파이프 시스템에도 적용된다 이 과정은 많은 수의 파이프 시스템에도 적용된다. 시스템을 통과하는 전체유동 중 지름, 길이 및 조도를 기초로 하여 파이프 1의 통과비율을 추정하여 를 현명하게 선택함으로써 식 (11.5.2)에 의해 몇 퍼센트의 허용오차내의 값을 얻을 수 있다. 이 값들은 마찰계수의 정확도 범위 내에서 잘 일치한다. <풀이> = 3cfs라고 가정한다. 그러면 = 3.82, = 3.82(1/0.00003)= 127,000, = 0.001, = 0.022이고 다음과 같다. 파이프 2에 대해서는

그리고 = 0. 00015이므로, = 0. 020으로 가정한다. 그러면 = 4. 01ft/s, = 4. 01(⅔)(1/0 그리고 = 0.00015이므로, = 0.020으로 가정한다. 그러면 = 4.01ft/s, = 4.01(⅔)(1/0.00003) = 89,000, = 0.019, = 4.11ft/s, = 1.44 cfs. 파이프 3에 대해서는, 그리고 = 0.0006이므로, = 0.020으로 가정한다. 그러면 = 4.01ft/s, = 4.01(1.333/0.00003) = 178,000, = 0.020, = 5.60 cfs으로 된다. 가정된 조건에 대한 총 유량은 따라서

h1 , h2 및 h3 의 값을 검토하면 는 대략 0. 018과 0. 019사이의 중간값이다. 만일 0 h1 , h2 및 h3 의 값을 검토하면 는 대략 0.018과 0.019사이의 중간값이다. 만일 0.018이 선택됐다면 = 20.4ft가 될 것이다. 를 구하기 위해 또는 여기서 평균 수두손실이 사용되었으므로 다음과 같이 된다.

5.8 Find the discharges for water at 20°C with the following pipe data and reservoir elevations as shown below. Where, 간단한 분기파이프 시스템이 나타나 있다. 이 문제에서 구하고자 하는 것은 수조의 높이가 주어졌을 때 각 파이프를 통과하는 유량이다. 파이프의 크기, 형태 및 유체의 성질을 알고 있다고 가정한다. 각 파이프에 대해 Darcy-Weisbach 방정식과 연속방정식이 만족되어야 한다.

연결점 에서 출입하는 유량은 같다. 유동은 가장 높은 수조에서 가장 낮은 수조로 흐르므로 연속방정식은 다음과 같다 연결점 에서 출입하는 유량은 같다. 유동은 가장 높은 수조에서 가장 낮은 수조로 흐르므로 연속방정식은 다음과 같다. 또는 연결점에서의 수력구배선의 높이가 중간에 위치한 수조의 높이보다 위에 있으면 유동은 그곳으로 흘러갈 것이나 만일 수력구배선의 높이가 그보다 낮은 곳에 위치하면 유동은 그곳에서 흘러나올 것이다. 국부손실은 등가길이로서 표시되고 파이프의 실제길이에 더해진다. 이 문제에 대한 해는 먼저 연결점에서의 수력구배선의 높이를 가정하고, , 및 를 계산하여, 이를 연속방정식에 대입함으로써 얻어질 수 있다. 만일 연결점에 흘러들어가는 유동이 너무 크다면, 유입을 줄이고 유출을 크게 하도록 구배선의 높이를 더 높이 가정한다.

[예제 11. 8] 그림 11. 11에서 다음과 같은 값들을 가질 메, 20〫 C 물에 대한 유량을 구하라

이므로 유입이 유출보다 다음 양만큼 더 크다. 라고 가정하면, 유입이 아직도 0. 029m3/s만큼 크다

5.9 Explain the algorithm of the Hardy Cross Method for the pipe network analysis. 관망문제를 해석적으로 해결하는 것은 현실적이 되지 못하므로 반복근사해법이 이용된다. Hardy Cross방법 (3)은 모든 연결점에서 연속방정식이 만족되도록 각 파이프에 대해 유동을 가정하는 방법이다. 그리고 각 회로에서의 유동에 대한 보정값을 차례로 계산하고, 이 값을 사용해서 회로들의 평형상태를 더욱 향상시킨다. 국부손실은 각 파이프의 등가길이로서 포함된다. 지수방정식은 보통 의 형태로 사용되며, 여기서 r은 식 (11.1.1)에서 이다. r의 값은 각각 관로에서 일정한 값이며 (Darcy-Weisbach의 식이 사용되지 않는 한), 루프-균형(loop-balancing) 과정에 앞서서 결정된다. 보정함은 아래와 같이 얻어진다.

임의 파이프에 대해 초기유량 Q0가 가정될 때, (11. 7. 1) 이다. 여기서 Q는 보정된 유량이고 는 보정량이다 임의 파이프에 대해 초기유량 Q0가 가정될 때, (11.7.1) 이다. 여기서 Q는 보정된 유량이고 는 보정량이다. 이때 각각의 파이프에 대해 만일 가 에 Q 비해 작다면 두 번째 항 이외의 모든 항은 무시될 수 있다. 지금 하나의 회로에 대해 이다.

는 회로내의 모든 파이프에서 동일하기 때문에 로부터 밖으로 나왔고, 절대값 기호는 회로에 따른 합의 방향을 고려하기 위해 부가되었다. 마지막 식을 풀면 관망중의 각 회로에 대한 를 구할 수 있다. (11.7.2) 를 식 (11.7.1)에 따라 회로내의 각 파이프에 적용할 때 방향이 매우 중요하다. 즉, 유동이 시계방향일 때는 더해 주고 반시계 방향일 때는 빼준다.

계산은 아래와 같은 단계를 밟는다. 1. 관망을 잘 살펴보아서, 연속식을 만족시키는 가장 좋은 유동분포를 가정한다. 2 계산은 아래와 같은 단계를 밟는다. 1. 관망을 잘 살펴보아서, 연속식을 만족시키는 가장 좋은 유동분포를 가정한다. 2. 하나의 기본회로내의 각 파이프에 대해서 순 수두손실을 계산하여 이들의 합 을 구한다. 또한 그 회로의 를 계산한다. 식 (11.7.2)에서 음의 비율로써 보정량을 구하고, 유량을 보정하기 위해서 이것을 회로내의 유동에 대수적으로 더해준다. 3. 다른 기본회로들에 대해 적용하여 2의 보정과정을 반복한다. 이를 모든 기본회로에 대해 계속한다. 4. 보정량(he)이 임의의 적은 양이 되기까지 2와 3의 과정을 여러 번 반복한다.

값은 분자와 분모 양쪽에 나타난다. 따라서 유량분포를 구하기 위해서는 실제의 r에 비례하는 값을 이용할 수 있다

관망 여러 회로들로부터 하나의 주어진 출구로 유동하도록 상호연결된 파이프를 파이프관망이라 한다 관망 여러 회로들로부터 하나의 주어진 출구로 유동하도록 상호연결된 파이프를 파이프관망이라 한다. 이것은 전기관망과 많은 점에서 유사하다. 이에 관한 문제는 일반적으로 복잡하며 기본회로들이 유동조건을 만족할 때까지 차례로 균형을 이루도록 시키는데, 시행 해법을 필요로 한다. 파이프관망에서는 아래 조건들이 반드시 만족되어야 한다. 1. 각 회로에서의 압력강하의 대수합은 0이어야 한다. (에너지방정식) 2. 각 연결부로 들어가는 유동과 나가는 유동은 같아야 한다. (연속방정식) 3. 각 파이프에 대한 Darcy-Weisbach방정식 혹은 이에 등가한 지수마찰공식이 만족되어야 한다. 즉, 수두손 실과 유량 사이의 적절한 관계가 각 파이프에서 유지 되어야 한다. (수두손실식 : 압력=마찰력, 운동방정식)

유체의 유동 고체의 움직임 질량보전 : 연속방정식 힘보전 : 운동방정식 에너지보전 : 에너지방정식 그림 11 유체의 유동 고체의 움직임 질량보전 : 연속방정식 힘보전 : 운동방정식 에너지보전 : 에너지방정식 그림 11.15 관망 첫 번째 조건은 회로의 어느 두 점, 예를 들면 그림 11.15의 A와 G사이의 압력강하는 파이프 AG를 통과하든지 흑은 AFEDG를 통과하든지 모두 같아야만 한다는 것이다. 두 번째 조건은 연속방정식이다.

5.10 The distribution of flow through the network as shown below is desired for the inflows and outflows as given. For simplicity n has been given the value 2.0. Implement the Hardy Cross Method for this problem. (a)

그림 11. 16(a)에는 가정된 유동분포가 표시되어 있다. 위쪽 좌측에는 아래쪽의 1번회로에 대한 의 항이 계산되어 있다 그림 11.16(a)에는 가정된 유동분포가 표시되어 있다. 위쪽 좌측에는 아래쪽의 1번회로에 대한 의 항이 계산되어 있다. 그 옆에는 같은 회로에 대한 를 계산한 것이다. 그림의 우측 위편에는 제2회로에 대한 것이 같은 형식으로 나타나 있다. 첫 번 단계에서 보정된 유량은 맨 위쪽 수평 관에서 15 + 11.06 = 26.06이 결정되고, 대각선 파이프는 35 + (-21.17) + (-11.06) = 2.77이다. 그림 11.16(b)에는 이와 같이 한 번 보정된 유량이 제시되어 있고 그림 11.16(c)는 4번째 보정된 값이 제시되어 있다.

그림 11.16에서 보여 지는 것과 같은 매우 단순한 관망의 경우는, 약 15개의 기억 용량과 100개의 프로그램 단계를 갖는 계산기가 있으면 수동으로 계산될 수 있다. 앞의 예보다 큰 관망에 대해서나 또는 복수의 수조나, 공급펌프, 승압펌프 둥을 포함하는 관망에 대해서는 Hardy Cross의 루프-균형해법에 의해 디지털 컴퓨터로 수치 해석할 수 있도록 프로그래밍 할 수 있다. 이에 대한 프로그램은 다음 절에서 취급된다.

그림 11.16 단순관망에서의 유동에 대한 해 (a) (b)

보다 일반적인 방법들[4-7]이 있으나, 이들은 기본적으로 Hardy Cross루프균형법, 또는 격자점균형(node balancing)법에 근거를 두고 있다. 이 시스템을 보다 일반적인방법으로 다루려고 하는 경우에는 한 쌍의 연립방정식으로 나타내고 그들을 Newton-Raphson법으로 푸는 것이 보통이다. 몇몇의 프로그램해법[5,6]에는 연결점압력과 유량에 및 붙여 파이프 크기 또는 조도가 미지수로 취급되어 있으므로 설계수단으로서 매우유용하다.

5.11 The program in Fig. 11.18 in the text book is used to solve the network problem displayed as below. The pump data are as follows: Implement the pipe network modeling for this problem by using the suggested program of Hardy Cross Method. 그림 5.1 0 0.03 0.06 0.09 30 29 26 20

<풀이> 그림 5.1과 같이 압력수두의 높이가 여러 곳에서 고정되어 있는 시스템의 경우는 수조에서의 미지의 유입량과 유출량을 고려하기 위해, 그리고 균형을 이를 때까지 연속식을 만족시키기 위해 가요소(pseudo elements)를 도입한다. 이의 요소에 의해 각 쌍의 고정압력수준 사이를 연결하면 가상적인 또는 인위적인 회로가 형성된다. 이의 요소들은 유동을 수반하지는 않지만, 수조들의 높이차와 같은 크기만큼의 수력구매선의수준을 강하시킨다. 만일의 요소에서 수두강하가 가정된 양의 방향으로 양의 값을 갖는다면 루프 3(그림 5.1)의 보정량은 다음과 같다. (11.8.1)

이 보정량은 단지 파이프 1과 4에만 적용된다. 만일 실제 관로가 가상루프 내에 추가로 존재한다면, 각각은 루프균형의 반복에 따라서 적절하게 조절될 것이다. 시스템내의 펌프는 그것을 통과하는 유량에 대응하는 펌프의 수두상승과 같은 크기의 음의 수두손실을 나타내는 유동요소로서 취급된다. 그림 5.1에서 요소 8의 펌프 양정-유량 곡선은 3차방정식으로 표시될 수 있다. 이때 는 점프의 차단양정이다. 루프 4의 보정량은 (11.8.2) 이 보정량은 루프 중의 파이프 5와 펌프 8에 적용된다. 펌프를 가진 관망을 만족스럽게 균형 시키기 위해서는, 양정-유량곡선의 기울기가 0보다 작거나 같아야만 한다.

모든 파이프에 대해 Hazen-Williams의 관로계수가 100이다. 그림 5. 2는 입력데이터를 나타내고, 그림 5 모든 파이프에 대해 Hazen-Williams의 관로계수가 100이다. 그림 5.2는 입력데이터를 나타내고, 그림 5.3은 이 문제의 컴퓨터 출력이다. DATA EXAMPLE 11.12 DATA SI,30,.001,.000001,100. DATA 5,HW DATA 1 , .12,600.,0.3,.0 DATA 2,.03,300., .15, .0 DATA 3,.0,500.,.6,.0 DATA 4,.03,400.,.3,.0 DATA 5,.03,30 0.,.3,.0 DATA 2,PS DATA 6,15.0 DATA 7,18.0 DATA 1,PU DATA 8,.06,.03,30.,29.,26.,20. DATA 16.IND DATA 3,2,1,-3,3,4,-5,3,3,6,-4,-1,3,5,7,8 DATA 1,NODES DATA 5,117. DATA 9,IX DATA 5,8,4,2,2,1,1,4,3 그림 5.2 그림5.1의 Hardy Cross 프로그램의 데이터

문제5. 11 SI UNITS SPEC. VISCOSITY IN M^2/S=. 000001 DESIRED TOLERANCE= 문제5.11 SI UNITS SPEC..VISCOSITY IN M^2/S= .000001 DESIRED TOLERANCE= .001 NO. OF ITERATIONS=30 PIPE Q(CFS OR M^3/S) L(FT OR M) D(FT OF M) HW C OR EFS 1 0.12000 600.00000 0.30000 100.00000 2 0.03000 300.00000 0.15000 100.00000 3 0.00000 500.00000 0.60000 100.00000 4 0.03000 400.00000 0.30000 100.00000 5 0.03000 300.00000 0.30000 100.00000 6 RESERVOIR ELEV. DIFFERENCE= 15 7 RESERVOIR ELEV. DIFFERENCE= 18 8 PUMP CURVE, DQ= .03 H= 30 29 26 20 COEF. IN PUMP EQ.=30 -11.11112 -555.5555 -6172.841 IND= 3 2 1 -3 3 4 -5 3 3 6 -4 -1 3 5 7 8

ITERATION NO. 1 SUM OF FLOW CORRECTION= 0. 1385 ITERATION NO ITERATION NO. 1 SUM OF FLOW CORRECTION= 0.1385 ITERATION NO. 2 SUM OF FLOW CORRECTION= 0.1040 ITERATION NO. 3 SUM OF FLOW CORRECTION= 0.0372 ITERATION NO. 4 SUM OF FLOW CORRECTION= 0.0034 ITERATION NO. 5 SUM OF FLOW CORRECTION= 0.0006 ELEMENT FLOW 1 0.143 2 -0.034 3 0.027 4 0.080 5 0.094 8 0.087 IX= 5 8 4 2 2 1 1 4 3

JUNCTION HEAD 1 137.811 2 150.044 3 135.044 4 137.797 5 117.000 그림 5.3 문제 5.11의 프로그램 출력