Code covered by the BSD License  

Highlights from
Experimental Investigation on Solar Chimney for room ventilation

image thumbnail

Experimental Investigation on Solar Chimney for room ventilation

by

 

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
% meetashish85@gmail.com,http://www.facebook.com/ashishmeet

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;

Contact us