Explore BrainMass
Share

Explore BrainMass

    Matlab -- Developing 22 Independent Material balances and solving in Matlab

    This content was COPIED from BrainMass.com - View the original, and get the already-completed solution here!

    Please see the attached file for the fully formatted problem(s).
    Find composition and flow rates from mass percent data.

    Here's what the 22 independent material balances should supposedly look like:

    1. W1 - W2 -W3=0
    2. H1 - H2 - H3 = 0
    3. S1 - S2 - S3 = 0
    4. H2 - H4- H5 = 0
    5. W2 - W5 = 0
    6. S2-S5 = 0
    7. W5-W7=0
    8. S5-S7 = 0
    9. H5-H7=0
    10. G6-G7=0
    11. W3+W7-W8=0
    12. H3+H7-H8=0
    13. G7-G8=0
    14. S3-S8+S7=0
    15. -0.685H2+H4=0
    16. W1+H1+S1=7500
    17. -0.42H5 + G6=0
    18. 0.38W1-0.62H1-0.62S1=0
    19. -0.05W1-0.05H1+0.95S1=0
    20. 0.38W2-0.68H2-0.68S2=0
    21. -0.05W2 -0.05H2+0.95S2 =0
    22. 0.021W8+0.021H8+0.021S8-0.979G8=0

    And here are the answers (composition and flow rates? Which is which?)
    x =

    1.0e+003 *

    4.6500
    2.2759
    2.3741
    2.2759
    2.2759
    4.6500
    2.4750
    1.0945
    1.3805
    0.7497
    0.3448
    0.3448
    1.7253
    0.3750
    0.1774
    0.1976
    0.1774
    0.1774
    0.3750
    0.1448
    0.1448
    0.1448

    And here are the variables set up in Excel. Matrix A will be all the coefficients on the left-hand side of the material balances (columns W1-G8) and Matrix B will be the right-hand side of the material balances (the column called 'Ans'). Then we solve for x by x=inv(a)*(b).

    W1 W2 W3 W5 W7 W8 H1 H2 H3 H4 H5 H7 H8 S1 S2 S3 S5 S7 S8 G6 G7 G8 Ans
    1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 -1 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 1 0 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0
    0 1 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 -1 0 0 0 0 0 0
    0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0 0
    0 0 1 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 1 0 0 1 -1 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 -1 0 0 0 0
    0 0 0 0 0 0 0 -0.685 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0
    1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 7500
    0 0 0 0 0 0 0 0 0 0 -0.42 0 0 0 0 0 0 0 0 1 0 0 0
    0.38 0 0 0 0 0 -0.62 0 0 0 0 0 0 -0.62 0 0 0 0 0 0 0 0 0
    -0.05 0 0 0 0 0 -0.05 0 0 0 0 0 0 0.95 0 0 0 0 0 0 0 0 0
    0 0.38 0 0 0 0 0 -0.68 0 0 0 0 0 0 -0.68 0 0 0 0 0 0 0 0
    0 -0.05 0 0 0 0 0 -0.05 0 0 0 0 0 0 0.95 0 0 0 0 0 0 0 0
    0 0 0 0 0 0.021 0 0 0 0 0 0 0.021 0 0 0 0 0 0.021 0 0 -0.979 0

    The Matlab Program
    %Matrix-gotta love it

    clear
    clc

    % variable portion of the matrix
    a=[1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 -1 0 0 0 0 0 0
    0 0 0 0 0 0 0 1 0 -1 -1 0 0 0 0 0 0 0 0 0 0 0
    0 1 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 -1 0 0 0 0 0
    0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0
    0 0 1 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 1 0 0 1 -1 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 -1 0 0 0
    0 0 0 0 0 0 0 -0.685 0 1 0 0 0 0 0 0 0 0 0 0 0 0
    1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 -0.42 0 0 0 0 0 0 0 0 1 0 0
    0.38 0 0 0 0 0 -0.62 0 0 0 0 0 0 -0.62 0 0 0 0 0 0 0 0
    -0.05 0 0 0 0 0 -0.05 0 0 0 0 0 0 0.95 0 0 0 0 0 0 0 0
    0 0.38 0 0 0 0 0 -0.68 0 0 0 0 0 0 -0.68 0 0 0 0 0 0 0
    0 -0.05 0 0 0 0 0 -0.05 0 0 0 0 0 0 0.95 0 0 0 0 0 0 0
    0 0 0 0 0 0.021 0 0 0 0 0 0 0.021 0 0 0 0 0 0.021 0 0 -0.979]

    % answer portion of the matrix
    b=[0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    7500
    0
    0
    0
    0
    0
    0
    ]

    x=inv(a)*(b)

    Test

    a =

    Columns 1 through 7

    1.0000 -1.0000 -1.0000 0 0 0 0
    0 0 0 0 0 0 1.0000
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 1.0000 0 -1.0000 0 0 0
    0 0 0 0 0 0 0
    0 0 0 1.0000 -1.0000 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 1.0000 0 1.0000 -1.0000 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    1.0000 0 0 0 0 0 1.0000
    0 0 0 0 0 0 0
    0.3800 0 0 0 0 0 -0.6200
    -0.0500 0 0 0 0 0 -0.0500
    0 0.3800 0 0 0 0 0
    0 -0.0500 0 0 0 0 0
    0 0 0 0 0 0.0210 0

    Columns 8 through 14

    0 0 0 0 0 0 0
    -1.0000 -1.0000 0 0 0 0 0
    0 0 0 0 0 0 1.0000
    1.0000 0 -1.0000 -1.0000 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 1.0000 -1.0000 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 1.0000 0 0 1.0000 -1.0000 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    -0.6850 0 1.0000 0 0 0 0
    0 0 0 0 0 0 1.0000
    0 0 0 -0.4200 0 0 0
    0 0 0 0 0 0 -0.6200
    0 0 0 0 0 0 0.9500
    -0.6800 0 0 0 0 0 0
    -0.0500 0 0 0 0 0 0
    0 0 0 0 0 0.0210 0

    Columns 15 through 21

    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    -1.0000 -1.0000 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    1.0000 0 -1.0000 0 0 0 0
    0 0 0 0 0 0 0
    0 0 1.0000 -1.0000 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 1.0000 -1.0000
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 1.0000
    0 1.0000 0 1.0000 -1.0000 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    0 0 0 0 0 1.0000 0
    0 0 0 0 0 0 0
    0 0 0 0 0 0 0
    -0.6800 0 0 0 0 0 0
    0.9500 0 0 0 0 0 0
    0 0 0 0 0.0210 0 0

    Column 22

    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    -1.0000
    0
    0
    0
    0
    0
    0
    0
    0
    -0.9790

    b =

    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    0
    7500
    0
    0
    0
    0
    0
    0

    x =

    1.0e+003 *

    4.6500
    2.2759
    2.3741
    2.2759
    2.2759
    4.6500
    2.4750
    1.0945
    1.3805
    0.7497
    0.3448
    0.3448
    1.7253
    0.3750
    0.1774
    0.1976
    0.1774
    0.1774
    0.3750
    0.1448
    0.1448
    0.1448

    © BrainMass Inc. brainmass.com October 9, 2019, 3:39 pm ad1c9bdddf
    https://brainmass.com/engineering/materials-engineering/matlab-developing-independent-material-balances-solving-matlab-9739

    Attachments

    Solution Preview

    Please see the attachment.

    To answer your questions:
    The values in x that you have are

    W1

    W2

    W3

    W5

    W7

    ...

    G6

    G7

    G8

    ]

    which are the entries in the heading of your table except the last column. All these entries are in kg/hr.

    How do you solve this problem:

    You have 4 component flowrates per stream. You have 8 streams so you have a maximum of 32 unknowns. You know that streams 1,2,and 3 do not have glue, thus you set the flowrate of this component in these streams to zero and you are left with 29 unknowns. You also know that stream 4 is pure water so you set the flowrates of the other components to zero. You are left with 26 unknowns. You also know that stream 5 does not have glue so one more unknown out, you're left with 25 unknowns. You know that stream 6 is pure glue so you have three other flowrates as zero, thus you are left with 22 ...

    Solution Summary

    The composition and flow rate of a feed stream is calculated using Matlab. The 22 independent material balances and solving in Matlab are developed.

    $2.19