기초전산유체역학 과제2

소개글

"기초전산유체역학 과제2"에 대한 내용입니다.

목차

1. 개요
2.1. Code (Block Tridiagonal Matrix)
2.2. Code (Gauss-Seidel Iteration)
3. 결과

본문내용

1. 개요
Block Tridiagonal Matrix와 Gauss-Siedel Iteration, 두 가지 방법을 통해 Laplace Equation을 matlab을 통해 수치적으로 풀어낸 뒤, 두 방법으로 구한 data를 3차원 그래프로 시각화하여 유사한 결과 값이 나오는 지에 대해 알아본다.

2.1. Code (Block Tridiagonal Matrix)
%% Block Tridiagonal matrix
clear;

% Initial Condition
Lx=5; Ly=5; % 시작 지점의 좌표는 (0,0)
nx=20; ny=20;

dx=Lx/nx; dy=Ly/ny; % dx=dy

if dx ~= dy
fprintf ('Error : dx is not equal to dy. \n Please edit the value of nx or ny. \n');
pause
end

% Populate the diagonals
A=btridiagonal(nx,ny);

% Boundary condition (bc = bc1 + bc2)
bc = zeros((nx-1)*(ny-1), 1);
for i = 1 : ny-1

참고 자료

없음

압축파일 내 파일목록

assignment 2.hwp