MATLAB for Chemical engineer

Transcription

1 University of Baghdad College of Engineering Department of Chemical Engineering MATLAB for Chemical engineer Basic and Applications Lecture No. 9 Dr. Mahmood Khazzal Hummadi Samar Kareem 2016

2 Lecture No. 9 Looping statements and descision statements: The conditions are given below. The column of the left hand side is to compare two values and those on the right hand side are to add multiple conditions: 1. LOOPING: Looping or iterative functions are extremely useful in engineering problem solving. 2 MKH

3 1. For statements loop a specific number of times, and keep track of each iteration with an incrementing index variable. A loop statement allows us to execute a statement or group of statements multiple times and following is the general form of a loop statement in most of the programming languages. For example, preallocate a 10-element vector, and calculate five values: x = ones(1,10); for n = 2:6 x(n) = 2 * x(n - 1); 2. While statements loop as long as a condition remains true. For example, find the first integer n for which factorial(n) is a 100-digit number: n = 1; nfactorial = 1; while nfactorial < 1e100 n = n + 1; nfactorial = nfactorial * n; Each loop requires the keyword. 3 MKH

4 It is a good idea to indent the loops for readability, especially when they are nested (that is, when one loop contains another loop): A = zeros(5,100); for m = 1:5 for n = 1:100 A(m, n) = 1/(m + n - 1); You can programmatically exit a loop using a break statement, or skip to the next iteration of a loop using a continue statement. Exit Loop before Expression Is False. Example: Sum a sequence of random numbers until the next random number is greater than an upper limit. Then, exit the loop using a break statement. limit = 0.8; s = 0; while 1 tmp = rand; if tmp > limit break s = s + tmp; 4 MKH

5 2. Decision Decision making structures require that the programmer should specify one or more conditions to be evaluated or tested by the program, along with a statement or statements to be executed if the condition is determined to be true, and optionally, other statements to be executed if the condition is determined to be false. 1. if... statement If the expression evaluates to true, then the block of code inside the if statement will be executed. If the expression evaluates to false, then the first set of code after the statement will be executed. Example: Create a script file and type the following code : a= 10; % check the condition using if statement if a < 20 % if condition is true then print the following fprintf('a is less than 20\n' ); fprintf('value of a is : %d\n', a); 5 MKH

6 2. if...else... statement An if statement can be followed by an optional else statement, which executes when the expression is false. If the boolean expression evaluates to true, then the if block of code will be executed, otherwise else block of code will be executed. If the boolean expression evaluates to true, then the if block of code will be executed, otherwise else block of code will be executed. Example Create a script file and type the following code a = 100; % check the boolean condition if a < 20 else % if condition is true then print the following fprintf('a is less than 20\n' ); % if condition is false then print the following fprintf('a is not less than 20\n' ); fprintf('value of a is : %d\n', a); 3. if...elseif...elseif...else... Statements. An if statement can be followed by one (or more) optional elseif... and anelse statement, which is very useful to test various conditions. When using if... elseif...else statements, are few points to keep in mind: there An if can have zero or one else's and it must come after any elseif's. An if can have zero to many elseif's and they must come before the else. Once an else if succeeds, none of the remaining elseif's or else's will be tested. if <expression 1> % Executes when the expression 1 is true <statement(s)> 6 MKH

7 elseif <expression 2> % Executes when the boolean expression 2 is true <statement(s)> Elseif <expression 3> % Executes when the boolean expression 3 is true <statement(s)> else % executes when the none of the above condition is true <statement(s)> Example: a = 100; %check the boolean condition if a == 10 % if condition is true then print the following fprintf('value of a is 10\n' ); elseif( a == 20 ) % if else if condition is true fprintf('value of a is 20\n' ); elseif a == 30 else % if else if condition is true fprintf('value of a is 30\n' ); % if none of the conditions is true ' fprintf('none of the values are matching\n'); fprintf('exact value of a is: %d\n', a ); 4. The Nested if Statements It is always legal in MATLAB to nest if-else statements which means you can use one if or elseif statement inside another if or elseif statement(s). Syntax 7 MKH

8 The syntax for a nested if statement is as follows if <expression 1> % Executes when the boolean expression 1 is true if <expression 2> % Executes when the boolean expression 2 is true You can nest elseif...else in the similar way as you have nested if statement. Example Create a script file and type the following code in it a = 100; b = 200; % check the boolean condition if( a == 100 ) % if condition is true then check the following if( b == 200 ) % if condition is true then print the following fprintf('value of a is 100 and b is 200\n' ); fprintf('exact value of a is : %d\n', a ); fprintf('exact value of b is : %d\n', b ); Answer Value of a is 100 and b is 200 Exact value of a is : 100 Exact value of b is : The switch Statement A switch block conditionally executes one set of statements from several choices. Each choice is covered by a case statement. An evaluated switch_expression is a scalar or string. An evaluated case_expression is a scalar, a string or a cell array of scalars or strings. 8 MKH

9 The switch block tests each case until one of the cases is true. A case is true when For numbers, eq(case_expression,switch_expression). For strings, strcmp(case_expression,switch_expression). For objects that support theeq(case_expression,switch_expression). For a cell array case_expression, at least one of the elements of the cell array matches switch_expression, as defined above for numbers, strings and objects. When a case is true, MATLAB executes the corresponding statements and then exits the switch block. The otherwise block is optional and executes only when no case is true. Syntax The syntax of switch statement in MATLAB is switch <switch_expression> case <case_expression> <statements> case <case_expression> <statements> otherwise <statements> Example Create a script file and type the following code in it grade = 'B'; switch(grade) case 'A' fprintf('excellent!\n' ); 9 MKH

10 case 'B' fprintf('well done\n' ); case 'C' fprintf('well done\n' ); case 'D' fprintf('you passed\n' ); case 'F' fprintf('better try again\n' ); otherwise fprintf('invalid grade\n' ); When you run the file, it displays Well done 6. The Nested switch Statements It is possible to have a switch as part of the statement sequence of an outer switch. Even if the case constants of the inner and outer switch contain common values, no conflicts will arise. Syntax The syntax for a nested switch statement is as follows switch(ch1) case 'A' fprintf('this A is part of outer switch'); switch(ch2) case 'A' fprintf('this A is part of inner switch' ); case 'B' fprintf('this B is part of inner switch' ); 10 MKH

11 case 'B' fprintf('this B is part of outer switch' ); Example Create a script file and type the following code in it a = 100; b = 200; switch(a) case 100 fprintf('this is part of outer switch %d\n', a ); switch(b) case 200 fprintf('this is part of inner switch %d\n', a ); fprintf('exact value of a is : %d\n', a ); fprintf('exact value of b is : %d\n', b ); When you run the file, it displays This is part of outer switch 100 This is part of inner switch 100 Exact value of a is : 100 Exact value of b is : MKH

12 Problems Q1) Find the value of y. Use function. A) B) Q2)Write the code required to determine the value of Q. The arrays Z and Y are generated as following: Z i =i+1, Y i =i 2 -i.use Subroutine Q N i0 Y i N! N i0 Z! i Q3)Solve the second order algebraic equation with one variable ax 2 + bx + c =0. The solution. Note use Function. b x b2 4ac 2a Q4) Calculate the kinetic and potential energies of an object at different velocities and heights respectively.. Use velocities and heights ranged from 1 to 1000 m ( step 1). Q5) Make an array for the numbers from 0 to 100, then make a second one represents the results of multiplying each sequential two numbers in the first array. Q6) The one dimension arrays A and B,and the number for the each array are 20 value. calculate the summation of array A and the summation of array B and find the y value, where A(i)=i and B(i)= i^2, where the value of y is calculate as following : Where MB is the average of B array. 12 MKH

13 Q7) Design windows application to find the densities of different six gases ( O 2,H 2, N 2, CO 2, CH 4,C 3 H 8, C 2 H 2 ) at temperatures ranged from 0 to 300 o C (step 10 ). The user should enter pressure in atm. Q8)Use loop to find the bubble point of ternary system (Ethanol 40 mol%, Water 20 mol% and Benzene 40 mol%). Knowing that the vapor pressure for three components are calculated by: Ethanol Po e=exp( /(t-50.5)) Water Po w=exp( /(t-46.13)) Benzene Po b=exp( /(t-52.36)) Where Ki= P o i /Pt,Pt=760,yi =Kixi, At Bubble point Σyi=ΣKi xi =1 Q9)Given that the vapor pressure of methyl chloride at K is bar, write a code to calculate the molar volume of saturated vapor at these conditions using Redlich/Kwong equation. Knowing; a= *r 2 T c 2.5 /Pc b= *rt c /P c V i+1 =(RT/P)+b- (a*(v i -b))/(t 1/2 PV i (V i +b)) R=83.14, Tc= k, Pc= 66.8 bar Q10) A simple force balance on a spherical particle reaching terminal velocity in a fluid is given by; Where: V t : Terminal velocity in m/s g : Acceleration of gravity p p : Particle density D p : The diameter of the spherical particle in m C D : Dimensionless drag coefficient. The drag coefficient on a spherical particle at terminal velocity varies with Reynolds number (Re) as followings: C D =24/Re for Re< 0.1 C D =24*(1+0.14* Re^0.7)/Re for 0.1=<Re=< 1000 CD=0.44 for 1000<Re=< CD=0.19-8*10^4/Re for <Re Q11)Write a computer program to calculate molar volume and compressibility factor for gaseous ammonia at a pressure values (1-30 atm) by using of the Redlich-Kwong equation of state. The equations and variables are listed below. R= gas constant (R = atm L/g-mol K),Tc= the critical temperature (405.5 K for ammonia) Pc= the critical pressure (111.3 atm for ammonia) Compressibility factor is given by Z= PV/RT 13 MKH

14 Q12)Calculate pressure drop (N/m 2 ) across a pipe of diameter d (m) and length L (m). The fluid properties are velocity u (m/s), density ρ (kg/m 3 ), and viscosity µ (kg/m.s). du Re For laminar flow: Re < 2100: 8 j f Re For turbulent flow: 10 5 > Re > 2100: j f 1/ 4 Re L 2 P 4 j f u d Q13) write the code required to calculate the density in g/lit of 50 mol% propane (C 3 H 8 )-50 mol% isobutene (i-c 4 H 8 ) gaseous mixture.the pressure in atm, temperature in K and critical constants Pc and Tc for propane and iso-butane.(note: Use function ). M p z z i ri 0.2T mix mix c ri M yi p, T p y i p y z PM z RT i i ri i T T ci ri lit. atm R gmol. k Comp. Pc, atm Tc, K M, g/gmol C 4 H C 3 H Q14) Estimate the heat capacity of difference gases (O 2,H 2, N 2, CO 2, CH 4,C 3 H 8, C 2 H 2 ) at different temperatures ( from 0 to 1000 o C step 10 )., Cp J/mol K, T deg K. component a b c d O E E E-08 H E E E-10 N E E E-08 CO E E E-09 CH E E E-08 C 3 H E E E-09 C 2 H E E E MKH