기초전산유체역학 과제2
- 최초 등록일
- 2021.05.17
- 최종 저작일
- 2018.04
- 8페이지/ 압축파일
- 가격 1,000원
소개글
"기초전산유체역학 과제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