이번년부터 석사생활을 하고 있습니다. 저희 실험실이 열유체역학을 기반으로 여러개의 팀으로 이루어져 있습니다.

주 특기는 mems기술이라 미세입자를 만드는 것입니다.

그래서 CFD나 코딩 관련해서 아는 분이 한분밖에 없으신데 지금 박사 졸업 논문 작성중이라서 많이 바쁘십니다.

석사생이지만 아는 분이 없으셔서 돈을 주고 인강을 신청 했습니다. 

교수님께서 저를 열유체 해석팀으로 키우고 싶으셔서 전산유체역학 수업을 추천하셔서 듣고 있습니다.

하지만, 이론적으로는 알겠는데 아무리 책을 보고 검색을 해봐도 도저히 감이 안잡힙니다. 

어떻게 해야 하는지 알려주세요. 

clear all; close all; clc;

%**************************************************************************

% Spatial and Temporal discretization

%**************************************************************************

Lx = 1.0; % Width

Nx = 31;

x = linspace(0,Lx,Nx);

dx = x(2) - x(1);

Ly = 1.0; % Width

Ny = 31;

y = linspace(0,Ly,Ny);

dy = y(2) - y(1);

maxIter = 500;

%**************************************************************************

% Properties of pure silver

%**************************************************************************

Rho = 1000;               % kg/m3

Cp = 2000;               % j/kg-k

k = 1;                % w/m-k

alpha = k/(Rho*Cp);         % m2/s

%**************************************************************************

% Initial Condition

%**************************************************************************

T  = zeros(Nx,Ny);

%**************************************************************************

% Boundary Condition

%**************************************************************************

T_Left  =   50;

T_Bottom  =   50;

T_Top   =   100;

T_Right   =   50;

T(:,1) = 30; % Left

T(:,end) = 30; % Right

T(1,:) = 100; % Top

T(end,:) = 30; % Bottom

%% Solving

for k=1:maxIter

    Told = T;

    for j=2:Ny-1

        for i=2:Nx-1

           T(i,j) = (T(i-1,j)+T(i+1,j)+T(i,j+1)+T(i,j-1))/4;

        end

    end

    numErr = abs (T - Told);

    if numErr < 1e-6

        disp(k)

        break

    end

end  

이건 FD로 한 코딩이구요. 박사행님이 도와주셔가지구요. 근대 이것도 맞는지 모르겠네요. 그리고 여러가지 iteration method를 사용하여 비교해봐라고 하는데 그건 어떻게 해야 하나요? method는 알지만 매트랩에서 코딩으로 하는 수식을 입력을 하기가 어렵네요.

그리고 이거 첨부파일로 올려드린 1번 문제에 대한 코딩인데 이까지는 하겠는데 여기서 막히네요.

이건 방정식에 대해 코딩을 할려고 하는데 쌍곡선이 들어가니까 어떻게 코딩을 해야 할지 막막합니다ㅠㅠ

clc ; clear all;

format long

L = 1.0;

W = 1.0;

x = [0:0.05:1];

for n = 1:100

    for i = 1:31 % x-dir

    for j = 1:31 % y-dir

        theta(n) = (((-1)^(n+1)+1)/n)*sin(n*pi*x(i)/L)*(

    end

end

S = sum(theta)

for i=1:4

    for j=1:4

        disp(j*(i-1)+j)

    end

end  

그리고 1번과 2번 문제에 대한 전반적으로 어떻게 접근해야 하는지 가르쳐주세요ㅜㅜ. 

교수님께 찾아가도 가르쳐 주시길 않네요. 개인적으로 공부해야 하는거라며ㅠㅠㅠ 

방법은 사비를 들어서 인강을 신청하는 방법밖에 없었습니다. 

지금 패키지로 매트랩 모든 강의를 신청하고 일주일내로 밤을 새서라도 들을 예정입니다. 

감사합니다.