## BALANCING AND STUDYING EFFECTS OF A CHEMICAL COMBUSTION EQUATION IN MATLAB

STOICHIOMETRIC COMBUSTION

DESCRIPTION:

Stoichiometric combustion defines the ideal combustion that needs to take place for a given chemical equation. This is actually the general form of equation for alkanes.

The general form might differ for each of the Carbon family: Alkanes, Alkenes and Alkynes as shown below,

FOR ALKANES, C will have 4 single bonds.

Alkanes: CnH2n+2

FOR ALKENES, C will have 2 double bonds.

Alkenes: CnH2n

FOR ALKYNES, C will have a triple bond.

Alkynes: CnH2n−2

Based on this concept, we should be balancing the equation.

The Stoichiometric co-efficient refers to the co-efficient of moles needed in a chemical equation to require balance through out. So, to find the stoichiometric co-efficient of 'ar' given here, we need to balance the whole equation with respect to the elements present in the reactant and the product sides.

SOLVING PROCEDURE FOR ALKANES:

1. Number of moles of C in LHS and RHS:
n = a

2.  Number of moles of H in LHS and RHS:
⦁ Alkanes: 2n+2 = 2b

3. Number of moles of N in LHS and RHS:
3.76ar = 2c

4. Number of moles of O in LHS and RHS:
2ar = 2a+b

Solving for ar value:

2*ar = (2*a) + b

2*ar = (2*n) + (2*n +1)

ar = n + (n+1)/2

ar = ((3*n)+1)/2

SOLVING PROCEDURE FOR ALKENES:

CnH2n + ar(O2 + 3.76N2) = a(CO2) + b(H2O) + c(N2)

1. Number of moles of C in LHS and RHS:
n = a

2.  Number of moles of H in LHS and RHS:
⦁ Alkenes: 2n = 2b

3. Number of moles of N in LHS and RHS:
3.76ar = 2c

4. Number of moles of O in LHS and RHS:
2ar = 2a+b

Solving for ar value:

2*ar = (2*a) + b

2*ar = (2*n) + (n)

ar = (3*n)/2

SOLVING PROCEDURE FOR ALKYNES:

CnH2n-2 + ar(O2 + 3.76N2) = a(CO2) + b(H2O) + c(N2)

1. Number of moles of C in LHS and RHS:
n = a

2.  Number of moles of H in LHS and RHS:
⦁ Alkynes: 2n-2 = 2b

3. Number of moles of N in LHS and RHS:
3.76ar = 2c

4. Number of moles of O in LHS and RHS:
2ar = 2a+b

Solving for ar value:

2*ar = (2*a) + b

2*ar = (2*n) + (n-1)

ar = n + ((n-1)/2)

ar = ((3*n)-1)/2

SOLVED FORMULA:

Alkanes: ar = [(3/2) x n] + 0.5

Alkenes: ar = [3/2 x n]

Alkynes: ar = [(3/2) x n] - 0.5

CODE:

close all
clear all
clc

% a = n;
% b = n+1;
% c = ar;
n = linspace(1,10,10)
% for the given equation in the challenge

%assingning ar value for alkane, alkene and alkyne
ar_alkane = (3*n+1)/2;
ar_alkene = 3*n/2;
ar_alkyne = (3*n-1)/2;

figure(1)
plot(n,ar_alkane,'b','linewidth',1.5)
hold on
plot(n,ar_alkene,'g','linewidth',1.5)
plot(n,ar_alkyne,'r','linewidth',1.5)
legend('Alkane','Alkene','Alkyne','location','northwest')
xlabel('no of moles')
ylabel('stoichiometric co-eff')
title('Stoichiometric Combustion')

PLOT:

### DETERMINATION AND STUDY OF EIGEN VALUE AND SPECTRAL RADIUS BY ITERATIVE METHODS Yokesh R · 2019-05-29 14:50:23

MAIN CODE: clear all close all clc A=[5 1 2 ; -3 9 4; 1 2 -7]; B=[10 ;-14; 33]; u = [0 1 2; 0 0 4; 0 0 0]; l = [0 0 0; -3 0 0; 1 2 0]; di = [5 0 0; 0 9 0; 0 0 -7]; mag = [0.25,0.5,0.75,1,1.5,2]; for i = 1 : length(mag) d = di*mag(i); Read more

### STUDY ON BEHAVIOR OF ITERATIVE TECHNIQUES FOR TWO DIMENSIONAL CONDUCTION DISTRIBUTION Yokesh R · 2019-05-25 10:39:05

DESCRIPTION: In this code, we are going to solve the two dimesnional convection equation defining a temperture distribution. The equation can be solved by 2 methods. 1. Implicit method 2. Explicit method In explicit method, we will be solving the equation sequencial Read more

### TAYLOR TABLE MODELLING FOR FOURTH-ORDER-APPROXIMATION OF SECOND ORDER DERIVATIVE Yokesh R · 2019-05-23 13:49:10

DESCRIPTION: Using the central and skewed scheme methods initially, we have to develop the Taylor table, which should be converted into a matrix and solved for the constants. These constants help in solving the second order differential equation. Finally, the error res Read more

### STUDY ON 1-D CONVECTION DISTRIBUTION FOR DIFFERENT TIME-STEPS Yokesh R · 2019-05-23 06:51:17

DESCRIPTION: In this code, we are going to compare the results of the 1D convection distribution equation describing the velocity for a range of time-step values. As a result, we will be getting the plots on the final velocity profile on time marching and computational Read more

### DISCRETIZATION OF FIRST ORDER DERIVATIVE USING FIRST, SECOND AND FOURTH ORDER APPX Yokesh R · 2019-05-22 11:24:56

DESCRIPTION:      In this code, we are going to discretize the function for the first-order derivative using first, second and fourth order approximation techniques and studying the error bounced by every method by a bar chart.   CODE: clear all Read more

### DISCRETIZATION OF A FIRST ORDER DERIVATIVE USING A RANGE OF dx TERMS Yokesh R · 2019-05-22 10:39:19

DESCRIPTION:      Here in this code, we are using a range of dx terms for discretizing the first order derivative and studying the results of the dx vs error plot  CODE: clear all close all clc x = pi/3; dx = linspace(pi/4,pi/4000,30); %y Read more

### STUDY ON EFFECT OF GRID POINTS IN 1D-CONVECTION DISTRIBUTION Yokesh R · 2019-05-10 07:29:29

DESCRIPTION:         In this report, we will be studying the effects of the grid points in the first order convection equation describing one dimensional velocity distribution.   frac{delu}{delt} = c*frac{delu}{delx}   CODE: cl Read more