수치해석 Homework 환경공학과 20061469 안 영 수.

Presentation on theme: "수치해석 Homework 환경공학과 20061469 안 영 수."— Presentation transcript:

1 수치해석 Homework 환경공학과 안 영 수

2 Part II I. 유한차분법

3 I.1 물질이동식을 해석하기 위한 범용적 유한차분법 알고리즘을 유도하라.

4

5 I.2 교과서에서 제시된 일반적 유한차분법에 대한 해를 Fortran 및 Basic 프로그램으로 개발한 전산모형과 Excel을 이용한 Excel 모형으로 구하여라. 확산계수 및 유속은 교과서에 있는 값을 사용하여라. 모델링 값을 비교 분석하여 그림에 나타내어라.

6 그림 2.4 M1.FOR 모형의 계산 결과 (유속에 의한 이동이 지배적인 경우)

7 그림 2.5 M1.FOR 모형의 계산 결과 (확산 및 유속 이동이 존재하는 경우)

8 그림 2.6 M1.FOR 모형의 계산 결과 (확산에 의한 이동이 지배적인 경우)

9 그림 2.7 M1-Excel 모형의 운영

10 그림 2.8 M1-Excel 모형의 계산 결과

11 I.3 편미분방정식의 특성식을 설명하고 이 특성식의 근의 수에 따라 편미분방정식을 분류하는 방법을 서술하고, 공학 및 자연과학 문제에 있어서 지배방정식을 이러한 분류된 편미분방정식의 어느 범주에 속하는지 서술하라.

12 편미분 방정식은 특성식에 따라 쌍곡형 (hyperbolic), 포물형(parabolic), 타원형(elliptic) 으로 나누어진다.
편미분방정식은 특성식의 속성에 따라 다음과 같 이 분류된다.

13 위의 편미분방정식중 쌍곡선형인 경우에는 수치 해석 불안정성과 오차가 심하기 때문에 격자망을 작게 사용하거나, 특성법을 사용하여야 한다. 특 성법은 다음과 같이 설명된다. 연속방정식에 전방 차분법을 적용하면 다음과 같다. 위의 전방차분 아날로그에 대하여 특성법에서는 주변수가 유속을 따라 이동한다고 가정한다. 즉. 시간과 공간격자의 크기를 특성선상의 유속에 따 라 다음과 같이 결정한다.

14 따라서, 다음의 식으로 전개된다. 즉, 주변수가 격자망을 따라서 유속에 의하여 이 동하는 것을 알 수 있다. 물질이동식에 GCA의 알고리즘을 적용하기 위해 서는 유속항은 중앙차분법을, 확산항은 특성법을 적용한다. 즉 분산이 유속에 따라 이동된 격자에 서 일어난 것으로 가정한다. 따라서, 분산항은 다 음과 같이 평가된다.

15 일반적으로 지하수 흐름의 지하수 유동 방정식 및 오염물 이동에 관계된 물질이동방정식은 포물형의 편미분 방정식 유형에 속한다
일반적으로 지하수 흐름의 지하수 유동 방정식 및 오염물 이동에 관계된 물질이동방정식은 포물형의 편미분 방정식 유형에 속한다. 이러한 일차원 물질이동방정식에 대하여 여러 수치해석 기법을 적용하여 일차원 모형을 개발하였 다. Excel을 사용하여 모형을 개발함으로서 수치해석 알고 리즘에 대하여 계산결과를 상세히 확인할 수 있었다. 이러 한 계산결과는 저자가 기존에 개발한 BASIC 및 FORTRAN을 이용한 수치해석모형과 비교하여 수치해석 기법의 타당성을 검증하였다. 또한 모든 수치해석 기법의 결과는 이론해와 비교하였다. 수치해석 알고리즘으로 유 한차분법 방법중 범용적 Crank-Nicholson 해법 (GCN : Generalized Crank-Nicholson Method)과 범용적 특성 평 균법 (GCA : Generalized Characteristic Averaging Method)를 사용하였다.

16 I.4 GCA 방법에 의한 물질이동식의 수치해석 기법을 설명하고, 해를 Fortran 및 Basic 프로그램으로 개발한 전산모형과 Excel을 이용한 Excel 모형으로 구하여라. 확산계수 및 유속은 교과서에 있는 값을 사용하여라. 모델링 값을 비교 분석하여 그림에 나타내어라.

17 GCA방법에 의한 알고리즘은 다음과 같이 특성법 (Characteristic Method), 중앙차분평균법 (Centered Difference (Averaging) Method), 범용 적(Generalized) Crank Nicholson 법의 결합에 의 해 유도된다.

18 1)특성법 편미분방정식은 특성식의 속성에 따라 다음과 같이 분류된다.

19 따라서, 다음의 식으로 전개된다. 즉, 주변수가 격자망을 따라서 유속에 의하여 이 동하는 것을 알 수 있다. 물질이동식에 GCA의 알고리즘을 적용하기 위해 서는 유속항은 중앙차분법을, 확산항은 특성법을 적용한다. 즉 분산이 유속에 따라 이동된 격자에 서 일어난 것으로 가정한다. 따라서, 분산항은 다 음과 같이 평가된다.

20 2) 중앙차분(평균)법 (Centered Difference (Averaging) Method)
물질이동식중 유속에 의한 항은 중앙차분법으로 평가된다. 즉, 시간도함수는 공간절점에 대하여 평균치를 취하고, 공간도함수는 시간절점에 대하 여 평균치를 취한다.

21 위의 식을 정리하면 다음과 같다. 여기서 , , 이다. 만약 평균을 취하지 않는 다면, 위의 식은 다음과 같다 (전방차분특성법). 인 경우, 두 평균 차분 특성법이나 전방 차분 특성법 모두 다음과 같이 유속에 따라 주변수를 추적하는 알고리즘이 되므로 특성법을 입자추적 법이라고도 한다.

22 물질이동식의 수치해석기법 유한요소법은 수치해의 오차를 최소화하도록 해 를 구하는 방법이다. 즉, 수치 해에서는 좌변과 우 변이 다르므로, 좌변과 우변의 차이를 잔차라고 정의한다. 물질이동식의 잔차의 미분운영함수 L(C)는 다음과 같이 정의된다.

23 위의 잔차를 해석하려는 전체 공간에 대하여 최소 화하기 위하여 각 격자 혹은 요소에 대한 가중화 된 잔차를 다 더한 후에 잔차의 합을 0이 되게끔 식을 구성한다. 시간에 대한 미분식은 일반적인 차분법을 적용하므로 잔차는 공간에 대해서만 해 석한다. 물질이동식의 가중잔차의 최소해를 구하 기 위한 가중잔차식은 다음과 같다.

24 위의 방법을 가중잔차법(Weighted Residual Method)이라고 하며, 유 한요소법의 근본원리가 된다
위의 방법을 가중잔차법(Weighted Residual Method)이라고 하며, 유 한요소법의 근본원리가 된다. 즉, 공간 및 시간영역에 대하여 격자화 된 각 계산점에서 수치해의 오차가 최소화되도록, 각 격자점의 수치해 에 가중치를 곱하여 합한 다음, 전체 오차가 0이 되게끔 알고리즘을 설정하는 방법인 것이다. 가중잔차식은 부분적분과 Green의 정리를 이용하여, 확산의 이차 미분항이 일차 미분항으로 변환되고, 유출경 계조건과 결합된다. 일반적으로 컴퓨터를 사용한 미분 혹은 편미분 방정식의 해법에 있어 서 수학적 혹은 이론적인 엄밀해와는 달리 컴퓨터는 공간적, 시간적인 전체 문제영역을 연속적으로는 생각할 수가 없기 때문에 오직 이산화 혹은 격자환된 점(요소내의 절점)에서의 변화에 대해서만 계산이 가 능하다. 따라서, 어떤 수치해석 기법(유한차분법, 유한요소법 등)을 사용하던 간에 지배방정식에 관계되는 모든 주변수, 매개변수, 독립 변수, 자료 등을 이산화하여야 한다. 유한요소법에서는 다음과 같이 대표적인 변수, 파라미터 등을 기저함수(Basis Function) 혹은 형상함 수(Shape Function)를 이용하여 이산화한다.

25 기저함수는 공간영역의 격자화에만 관계되므로 현상자체의 과정에는 의존하지 않으므로, 각 요소 의 형상에 의해서만 결정된다
기저함수는 공간영역의 격자화에만 관계되므로 현상자체의 과정에는 의존하지 않으므로, 각 요소 의 형상에 의해서만 결정된다. 따라서, 기저함수 는 계산시간을 크게 감소시킬 수 있도록, 격자점 의 좌표계만 주어지면 프로그램 운영 시 초기에 1 회 평가된다. 즉, 물질이동식의 물리적, 화학적, 생 물학적 기작을 평가하기 이전에 입력 자료로서 주 어지는 격자망의 구성 방법 혹은 좌표계로서 평가 되는 것이다. 기저함수는 모든 Gauss 지점에서 구 해진 후, 요소행렬들의 적분이 수행될 때 조합된 다.

26 감사합니다.