% CFD EXPLICIT METHOD FOR A PLATE
% Fluid Mechanics II - Homework 3 - October 2017
% Carlos Fuentes Carrillo - 3� Aerospace Engineering
% Universidad Europea de Madrid
clear all;
close all;
clc;
% Boundary conditions
max_iter = 1000;
mesh_elements = 317;
wall_1 = 0;
wall_2 = 500;
wall_3 = 0;
wall_4 = 100;
% Initial_mesh creation
initial_mesh = zeros(mesh_elements);
initial_mesh(2:mesh_elements - 1,1) = wall_1;
initial_mesh(1,2:mesh_elements - 1) = wall_2;
initial_mesh(2:mesh_elements - 1,mesh_elements) = wall_3;
initial_mesh(mesh_elements,2:mesh_elements - 1) = wall_4;
% solution (i = column, j = row)
iter = 1;
M = initial_mesh;
Aux(:,:,iter) = initial_mesh(:,:);
Ware = zeros(mesh_elements - 2);
while iter <= max_iter
for i = 2:(mesh_elements - 1)
for j = 2:(mesh_elements - 1)
M(j,i) = 0.25*(M(j+1,i) + M(j,i-1) + M(j-1,i) + M(j,i+1));
Ware(j-1,i-1) = M(j,i);
if i==(mesh_elements - 1) && j==(mesh_elements - 1)
iter = iter + 1;
M(2:mesh_elements-1,2:mesh_elements-1) = Ware;
Aux(:,:,iter) = M(:,:);
end
M(j,i) = Aux(j,i,iter);
end
end
end
% Plots
imagesc(Aux(:,:,max_iter))
colorbar