THE COPPERBELT UNIVERSITY SCHOOL OF MINES AND MINERAL TENC! DEPARTMENT OF CHEMICAL ENGINE 2022/23 SESSIONAL EXAMINATION COURSE: NUMERICAL METHODS CODE: MA 311 INSTRUCTIONS: 1. The paper contains 7 questions. 2. Question 1 is COMPULSORY. Answer it plus any other 4 questions. 3. Total marks for the paper is 100. Each question carries equal marks. 4. Relevant additional information will be provided. 5. Scientific calculators can be used 6. Time allowed is 3 hours.c equations} & Question 1 is compulsory [Linear algebrai For the following set of linear equations: 4x, -23y - 3x, + 6x, #12 6x, #78, + 6.5%, -6x = 65 6.2514 + 5.5%, =16 al? x, +75 422x, +15.5%, 2x, (i) write the equations in the form aug A= [A|b|E] (79 use the Gauss-Jordan method to compute the solution of the system of equations (i) find A“, (iv) show that AA =T. [20] Question 2 [Regression] The data below is to be fitted to the following equation y a 0 1 2 ¥. 1 2 2 a 1 (i) Calculate co, ¢1 and c2 (show what summations need to be calculated). (#/) How do you find co, cy and ca if is set equal to zer0? (iti) Ifthe desired equation were y=a,xe"™ , how would you use the least squares to find a; and a» ? [20]Question 3 [Non-linear equations] For turbulent flow of fluid in a hydraulically smooth pipe, Prandtl’s universal resistance law relates the friction factor, fand the Reynolds number, Re, according to the following relationship: 1 4d 3 (o- ff 0.40+4.0in(ReJ FX OMG WF n( e i) ig 5-8 ~0-86F 16 (i) Definex =, and rearrange the above equation in the standard form: f (x) = 0, Compute f for Re = 10° using (ii) Bisection method; /"' = 0.002, f°! = 0.01 (ii) false position method, f'!! = 0.002, f = 0.01. Terminate your iterations when the relative error, &, < 5%. Report the friction factor obtained at the end of each iteration. [20] ~ O. Clee a, ; = Ba-AGq3 ~ F° C9CF M165 y Question 4 [Interpolation] ‘The following data were measured precisely: 4 - x ee t 2 2.1 2.2 27 8 3.4 Zz 6 7.752 10.256 36.576 66 125.168 (a) Use Newton interpolating polynomials to determine z at = 2.5. Make sure that you order your points to attain the most accurate results. What do your results tell you regarding the order of the polynomial used to generate the data? (b) Use a third-order Lagrange interpolating polynomial to determine y at r= 2.5. [20]Question 5 [Numerical differentiation) de Chemical reactions often follow the model W ic reaction rate, and # = reaction order, Given values of ¢ and de“dt, ke and # can aay linear regression of the logarithm of this equation toa{- a nloge } ‘Use this approach along with the following data to estimate k and m: .* where ¢ = concentration, ¢ = time, k = be evaluated by a 7 10 20 30 40 50 60 c 3.5 25 L8 13 09 0.6 [20] Question 6 [Numerical Integration] The dimensionless temperature of a fluid under steady, fully developed Jat cylindrical pipe with walls heated electrically is given by @=—4Z—1" E 4 wninar flow in a q 7 4 where r is.the 24 ~ dimensionless radius. The cup mixing (bulk) dimensionless temperature is given by ule (i) Determine @, analytically (ii) generate values of @ at Z= 0.4 and r= 0, 0.2, 0.4, 0.6, 0.8,1.0 (iti) use @ combination of Simpson’s 1/3 and 3/8 rules to estimate & numerically (iv) determine the absolute relative true error, |é,| . [20] vit %y yor Kg dey Se OY £4 ere 5 Yah -6 we oor - 0 G4 4b XeQuestion 7 [ODE-IVP] ‘Use (a) Euler’s and (b) the fourth-order RK method to solve ca , Da -2y+4e' ao over the range ‘= 0 to Lusing a step size of 0.2 with y (0) = 2 and z (0) = 4. [20]