Code covered by the BSD License

# Experimental Investigation on Solar Chimney for room ventilation

### Ashish Meshram (Meet) (view profile)

14 May 2013 (Updated )

This project implements experimental Investigation on Solar chimney for room ventilation

MatCrea(Coeff,Flux,Wall,Tf,Ta,Ls,z,W,d)
```function [A B Area M] = MatCrea(Coeff,Flux,Wall,Tf,Ta,Ls,z,W,d)
% This function creates the matrix i.e. the coefficients of temperature and
% right hand side vector along with that some dimensional attributes.
%
% INPUT:
%           Coeff = Structure containing various heat transfer
%                   coefficients
%           Flux = Structure containing solar radiation heat flux (W/m^2)
%           Wall = Structure containing properties of Wall
% OUTPUT:
%           A = coefficients of temperature
%           B = right hand side vector
%           Area = dimensional attributes (m^2)
%           M = Mass flow rate (kg/s)

% Implemented by ASHISH MESHRAM

Lw = Ls - z/2;
%---area of wall and glass (m^2)
Area.Aw = Lw*W;
Area.Ag = Ls*W;
%---cross sectional area of inlet and outlet to air flow channel (m^2)
Area.Ai = z*W;
Area.Ao = d*W;
Area.Ar = Area.Ao/Area.Ai;

M = (0.57*Wall.Rof*Area.Ao/((1 + Area.Ar^2)^(1/2))*((2*9.81*Ls*(Tf-Ta)/Ta))^(1/2));

%---Creating Temperature Coefficients
A(1,1) = Coeff.hg*Area.Ag + Coeff.hrwg*Area.Aw + Coeff.Ut*Area.Ag;
A(1,2) = -Coeff.hg*Area.Ag;
A(1,3) = -Coeff.hrwg*Area.Aw;
A(2,1) = Coeff.hg*Area.Ag;
A(2,2) = -(Coeff.hg*Area.Ag + Coeff.hw*Area.Aw + M*Wall.Cf/0.75);
A(2,3) = Coeff.hw*Area.Aw;
A(3,1) = -Coeff.hrwg*Area.Aw;
A(3,2) = -Coeff.hw*Area.Aw;
A(3,3) = Coeff.hw*Area.Aw + Coeff.hrwg*Area.Aw + Coeff.Ub*Area.Aw;

%---Creating Right Hand Side Vector
B(1,1) = Flux.S1*Area.Ag + Coeff.Ut*Area.Ag*Ta;
B(2,1) = -(M*Wall.Cf/0.75)*Ta;
B(3,1) = Flux.S2*Area.Aw + Coeff.Ub*Area.Aw*Ta;```