Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Upload

Welcome to Scribd!

  • 3 views

    Assignment 1

    Uploaded by

    Faisal Mumtaz
    The document describes solving a non-linear equation using MATLAB for two cases of a CSTR reactor in series. For case 1 with 5 reactors, it defines initial guesses, vectors for variables and the Jacobian, and iterates to solve the system of equations. For case 2 with 10 reactors, it similarly defines the initial guesses, vectors, Jacobian and iterates to solve the larger system of equations.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd

    Assignment 1

    Uploaded by

    Faisal Mumtaz
    3 views3 pages
    The document describes solving a non-linear equation using MATLAB for two cases of a CSTR reactor in series. For case 1 with 5 reactors, it defines initial guesses, vectors for variables and the Jacobian, and iterates to solve the system of equations. For case 2 with 10 reactors, it similarly defines the initial guesses, vectors, Jacobian and iterates to solve the larger system of equations.

    Original Title

    assignment-1

    Copyright

    © © All Rights Reserved

    Available Formats

    DOCX, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document describes solving a non-linear equation using MATLAB for two cases of a CSTR reactor in series. For case 1 with 5 reactors, it defines initial guesses, vectors for variables and the Jacobian, and iterates to solve the system of equations. For case 2 with 10 reactors, it similarly defines the initial guesses, vectors, Jacobian and iterates to solve the larger system of equations.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    3 views3 pages

    Assignment 1

    Uploaded by

    Faisal Mumtaz
    The document describes solving a non-linear equation using MATLAB for two cases of a CSTR reactor in series. For case 1 with 5 reactors, it defines initial guesses, vectors for variables and the Jacobian, and iterates to solve the system of equations. For case 2 with 10 reactors, it similarly defines the initial guesses, vectors, Jacobian and iterates to solve the larger system of equations.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    of 3

    Assignment-1:

    CSTR reactor in series, solve non-linear equation using Matlab


    Prepared by: Faisal Mumtaz

    Case 1: for n=5 :-

    % Defining Initial guess:

    >> a0=5;
    >> a1(1)=4;
    >> a2(1)=3;
    >> a3(1)=2;
    >> a4(1)=1;
    >> beta=5;
    >> a5=0.5;

    >> %Defining two vectors of x, y & jacobian:

    >> x=[a1(1); a2(1); a3(1); a4(1); beta(1);];


    >> jacobian=[2*beta(1)*a1(1)+1 0 0 0 a1(1)^2; -1 2*beta(1)*a2(1)+1 0 0 a2(1)^2; 0 -1
    2*beta(1)*a3(1)+1 0 a3(1)^2 ; 0 0 -1 2*beta(1)*a4(1)+1 a4(1)^2; 0 0 0 -1 a5^2];
    >> y=[beta(1)*a1(1)^2-(a0-a1(1)); beta(1)*a2(1)^2-(a1(1)-a2(1)); beta(1)*a3(1)^2-(a2(1)-a3(1));
    beta(1)*a4(1)^2-(a3(1)-a4(1)); beta(1)*a5^2-(a4(1)-a5)];

    >> % initial for iterative term j as 1;

    >> j=1;

    >> while norm(jacobian^-1*y) > 1e-2


    x = [a1(j); a2(j); a3(j); a4(j); beta(j)]
    jacobian = [2*beta(j)*a1(j)+1 0 0 0 a1(j)^2; -1 2*beta(j)*a2(j)+1 0 0 a2(j)^2; 0 -1 2*beta(j)*a3(j)+1 0
    a3(j)^2 ; 0 0 -1 2*beta(j)*a4(j)+1 a4(j)^2; 0 0 0 -1 a5^2];
    y=[beta(j)*a1(j)^2-(a0-a1(j)); beta(j)*a2(j)^2-(a1(j)-a2(j)); beta(j)*a3(j)^2-(a2(j)-a3(j)); beta(j)*a4(j)^2-
    (a3(j)-a4(j)); beta(j)*a5^2-(a4(j)-a5)];
    x = x - jacobian^-1*y;
    a1(j+1) = x(1);
    a2(j+1) = x(2);
    a3(j+1) = x(3);
    a4(j+1) = x(4);
    beta(j+1) = x(5);
    j=j+1;
    end
    Assignment-1:
    CSTR reactor in series, solve non-linear equation using Matlab
    Prepared by: Faisal Mumtaz

    Case 2: for n=10 :-

    % Defining Initial guess:

    >> a0=5;
    >> a1(1)=4.5;
    >> a2(1)=4;
    >> a3(1)=3.5;
    >> a4(1)=3;
    >> a5(1)=2.5;
    >> a6(1)=2;
    >> a7(1)=1.5;
    >> a8(1)=1;
    >> a9(1)=0.8;
    >> a10=0.5;
    >> beta=5;

    >> %Defining two vectors of x, y & jacobian:

    >> x=[a1(1); a2(1); a3(1); a4(1); a5(1); a6(1); a7(1); a8(1); a9(1); beta(1)];
    >> jacobian=[
    2*beta(1)*a1(1)+1 0 0 0 0 0 0 0 0 a1(1)^2;
    -1 2*beta(1)*a2(1)+1 0 0 0 0 0 0 0 a2(1)^2;
    0 -1 2*beta(1)*a3(1)+1 0 0 0 0 0 0 a3(1)^2;
    0 0 -1 2*beta(1)*a4(1)+1 0 0 0 0 0 a4(1)^2;
    0 0 0 -1 2*beta(1)*a5(1)+1 0 0 0 0 a5(1)^2;
    0 0 0 0 -1 2*beta(1)*a6(1)+1 0 0 0 a6(1)^2;
    0 0 0 0 0 -1 2*beta(1)*a7(1)+1 0 0 a7(1)^2;
    0 0 0 0 0 0 -1 2*beta(1)*a8(1)+1 0 a8(1)^2;
    0 0 0 0 0 0 0 -1 2*beta(1)*a9(1)+1 a9(1)^2;
    0 0 0 0 0 0 0 0 -1 a10^2;
    ];
    >> y=[

    beta(1)*a1(1)^2-(a0-a1(1));
    beta(1)*a2(1)^2-(a1(1)-a2(1));
    beta(1)*a3(1)^2-(a2(1)-a3(1));
    beta(1)*a4(1)^2-(a3(1)-a4(1));
    beta(1)*a5(1)^2-(a4(1)-a5(1));
    beta(1)*a6(1)^2-(a5(1)-a6(1));
    beta(1)*a7(1)^2-(a6(1)-a7(1));
    beta(1)*a8(1)^2-(a7(1)-a8(1));
    beta(1)*a9(1)^2-(a8(1)-a9(1));
    beta(1)*a10^2-(a9(1)-a10)];
    Assignment-1:
    CSTR reactor in series, solve non-linear equation using Matlab
    Prepared by: Faisal Mumtaz

    >> % initial for iterative term j as 1;

    >> j=1;

    >> while norm(jacobian^-1*y) > 1e-2


    x = [a1(j); a2(j); a3(j); a4(j); a5(j); a6(j); a7(j); a8(j); a9(j); beta(j)]

    jacobian = [
    2*beta(j)*a1(j)+1 0 0 0 0 0 0 0 0 a1(j)^2;
    -1 2*beta(j)*a2(j)+1 0 0 0 0 0 0 0 a2(j)^2;
    0 -1 2*beta(j)*a3(j)+1 0 0 0 0 0 0 a3(j)^2;
    0 0 -1 2*beta(j)*a4(j)+1 0 0 0 0 0 a4(j)^2;
    0 0 0 -1 2*beta(j)*a5(j)+1 0 0 0 0 a5(j)^2;
    0 0 0 0 -1 2*beta(j)*a6(j)+1 0 0 0 a6(j)^2;
    0 0 0 0 0 -1 2*beta(j)*a7(j)+1 0 0 a7(j)^2;
    0 0 0 0 0 0 -1 2*beta(j)*a8(j)+1 0 a8(j)^2;
    0 0 0 0 0 0 0 -1 2*beta(j)*a9(j)+1 a9(j)^2;
    0 0 0 0 0 0 0 0 -1 a10^2;
    ]

    y=[
    beta(j)*a1(j)^2-(a0-a1(j));
    beta(j)*a2(j)^2-(a1(j)-a2(j));
    beta(j)*a3(j)^2-(a2(j)-a3(j));
    beta(j)*a4(j)^2-(a3(j)-a4(j));
    beta(j)*a5(j)^2-(a4(j)-a5(j));
    beta(j)*a6(j)^2-(a5(j)-a6(j));
    beta(j)*a7(j)^2-(a6(j)-a7(j));
    beta(j)*a8(j)^2-(a7(j)-a8(j));
    beta(j)*a9(j)^2-(a8(j)-a9(j));
    beta(j)*a10^2-(a9(j)-a10)];

    x = x - jacobian^-1*y;
    a1(j+1) = x(1);
    a2(j+1) = x(2);
    a3(j+1) = x(3);
    a4(j+1) = x(4);
    a5(j+1) = x(5);
    a6(j+1) = x(6);
    a7(j+1) = x(7);
    a8(j+1) = x(8);
    a9(j+1) = x(9);
    beta(j+1) = x(10);
    j=j+1;
    end

    You might also like