Questions and Answers of Numerical Methods With Chemical Engineering Applications

Solve the following ordinary differential equation dy dx = Inx, y(1)=2
Solve Eq. (1.2.70) subject to an initial condition c(x, 0) = 1 and homogeneous boundary conditions at x = 0 and x = L. Әс at О a2c = D— 3х2 (1.2.70)
Solve the following ordinary differential equation y = 3xy + xy(3) = 1
Write a MATLAB function that determines the product n S= = П j=jo (1.4.1)
Solve the following ordinary differential equation dy dx =y-2, y(0) = 2
Write aMATLAB program that computes the product in Eq. (1.4.1) until it exceeds a particular value smax. S= п П j=jo (1.4.1)
Solve the following ordinary differential equation y" + 3y + 2y = 0, y(0) = 1, y(2) = 2
Write a MATLAB program that computes n terms in the sum n - Σ' sinj| j=1 S= (1.4.3)
Solve the following ordinary differential equation y"-2y, y(0) = 2, y(0) = 0 =
Solve the convection–diffusion–reaction problemsubject to no-flux at x = 0,and a concentration c = c0 at x = L. Before solving this problem, you should convert it into dimensionless form, wherein
Solve the following ordinary differential equation d²y dx² dy - 2- dx = 0, y(0) = 2, dy -(0) = 1 dx
Solve the following ordinary differential equation d'y dx² = 4y+1, y(0) = 1, dy dx -(0) = 0
Solve the following ordinary differential equation for the function y(x)subject to y(0) = 0 and y(1) = 1. y"+y=(1-²) sin лx
Solve the system of ordinary differential equationssubject to x(0) = 0 and y(0) = 1. dx dt dy dt = 3x + 2y || = 4x + 2y =
Determine the trajectory of the kinetic systemwith initial conditions A(0) = 1 and B(0) = 2. Does the system reach steady state? dA dt dB dt = 2A- 3B -4A + B
Use separation of variables to obtain an eigenfunction solution of the formfor the unsteady-diffusion equationsubject to the constant concentration boundary condition c(0, t) = 0 and the reaction–
Consider the reversible reaction A ⇌ B with a forward rate constant of 3 and a reverse rate constant of 2. If we start with 1 mole of A and 1 mole of B, compute the number of moles of each species
Use separation of variables to solve the unsteady-diffusion equationsubject to no flux at the left boundaryand zero concentration c(1, t) = 0 on the right boundary. The initial condition is c(x, 0) =
Consider the partial differential equationsubject to ∂y/∂x = 0 at x = 0 and y = 0 at x = 1. When you solve this problem using separation of variables, what are the eigenvalues? ду at а у zre
Solvesubject to the initial condition T(x, 0) = 1 and the boundary conditions T(0, t) = 0 and T(1, t) = 0. Report your answer in terms of the non-zero Fourier modes. aT at || а т ах2
Consider a box of length L = 1 that is initially in contact with a heat reservoir on the left at T = 1 and a heat reservoir on the right with T = 0. When this system reaches equilibrium, it will have
Consider a slab of thickness 2W that is initially at a temperature TH. The slab is in contact with a cold reservoir at some temperature TC < TH. The slab has a thermal conductivity k and the heat
What is the output of the following MATLAB code? 1 2 3 4 5 6 7 counter = 3; 8 stop = 0; 9 10 11 12 N 13 function out = 15 16 17 for i = 4:-1:1 x(i) = 1^2; end while stop == 0; end problem1_36 if x
In numerical methods, there are often many ways to make a calculation that give the same answer but have vastly different efficiencies. In this problem, we are going to consider a common calculation
Express the 103493.234 as a 16 bit word in base 10.
Using the LU decomposition in Problem 2.32, determine the value of x2 for Ax = b if the forcing function isData from Problem 2.32:In LU decomposition, determine the missing value of L for
What happens if we switch the order of the equations in Example 2.8?Example 2.8Solve the systemby Jacobi’s method with an initial guess (0, 0). 2x1 + x₂ = 2 x12x₂ = -2 (2.13.7)
Use LU decomposition to solve the system of equations 2x13x2 + x3 = 7 1x₁x22x3 = -2 3x1 + x2x3 = 0
Use LU decomposition with Doolittle’s method to solve the set of equations given by 1 2 1 -1 1 2 -2 1 2 3 3 5-59 1 X11 X12 X13 X21 X22 X23 X31 X32 X33 X41 X42 X43 2 9 6 15 1 12 -2 10 15 42 -2 3
Use LU decomposition (Doolitle’s method) to solve the system of equations 4x1 -2x1 1x1 + 3x2 X3 4x2 + 5x3 + 2x2 + 6x3 - - 1 2 23 3
Use LU decomposition to compute the inverse of 2 -2 4 4 1 14 3 1
Determine the solution tousing(a) Gauss elimination without pivoting(b) LU decomposition. Show the values of U, L, and y. 2 2 3 xX1 BAHA X2 X3 457 246 5 11 8
Solve the following problem(a) Using Cramer’s rule(b) Using Gauss elimination with partial pivoting(c) Using LU decomposition. [²2][3]-[2] y 12
For the matrixwhich definition of the norm gives the largest value and which definition gives the smallest value? You only need to consider the L1, Le, and L∞ norms. A = 6 1 2-2 1
Solve the following system of equationsusing(a) Naive Gauss elimination(b) Gauss elimination with pivoting(c) LU decomposition. 6 [][][] y 15 Z 6 1 2 3 4 5 6 1 3 2
Consider the matrixDetermine the values of x such that the Euclidean norm of the matrix is larger than the 2-norm. A = [2, 4] -1
What is the condition number on the system below using the Euclidean norm? 3x + 2y + z = 5 x-y + 2z = 5 x - 3z = -1
Compute the L1 norm for the matrix A || 2 -4 2 6-2 8 10 -13 9
Compute the L1, L2, and L∞ norms for the vector x = [10 3 –4 –1 5].
For the matrixwhich definition of the norm gives the largest value and which definition gives the smallest value? You only need to consider the L1, Le, and L∞ norms. A || 2 1 -5 11 -3 2 8 3 1 2 2
Compute the L1, Le, and L∞ norms for the matrix 1 2 3 2 -5 4 6 -3 1-4 5 2 216 10
Use the 1-norm to determine the condition number of A: || 1 2 1 1 12 402 1 22
Compute the condition number forusing the Euclidean norm. 1 1 12 13 0 1 1 2
Compute the condition number for the matrixusing the Euclidean norm. A = [ 1 ² ] 12 01
Compute the 1-norm, Euclidean norm, spectral norm, and infinity norm of 3 - [22] 1 A =
Consider the 2 × 2 system of equationsWhat is the value of the 1-norm of ||Ax(1) − b||1 after one iteration of Jacobi’s method with the initial guess x(0)1 = 1, x(0)2 = 0? 2x1 + x₂ = 4 x1 +
Perform the first two iterations of Jacobi’s method for the system of equationsusing an initial guess x = 0. 6 -2 1 2 10-3 2 3 8 -6 10 10 13 -2 9 X1 X2 X3 X4 || -9 9 17 13
Perform the first two iterations of Gauss–Seidel for the system of equationsusing an initial guess x = 0. 6 1 03 -28 2 3 16 10 -2 2 1 3 9 2 X3 X4 || -9 9 17 13
Determine x(1)2 for the system of equationsusing Gauss–Seidel with an initial guess of x1 = 0, x2 = 0, x3 = 0. 3 24 X1 [CHO 1 30 X2 24 1 X3 6 3 4
Show that the matrix A in Problem 2.58 is diagonally dominant.Data from Problem 2.58Perform the first two iterations of Jacobi’s method for the system of equationsusing an initial guess x = 0.
Is the matrix from Problem 2.46 diagonally dominant?Problem 2.46Compute the L1 norm for the matrix A || 2 -4 2 6-2 8 10 -13 9
Perform linear regression for the following data set. X 1.1 1.6 2.0 2.3 2.8 у 2.4 2.4 3.8 5.2 6.0 7.9
For numerical calculations, the computer is limited by two major constraints: the amount of memory available and the speed of the processor. Let’s examine how these constraints affect our ability
If I want to solve a 200 × 200 matrix with a bandwidth of 5 and p = q, what is the ratio of the time required for full Gauss elimination relative to using a banded Gauss elimination solver?
Let’s assume we have a computer that can solve a 1000 × 1000 system in 2 seconds using Gauss elimination. If this is in fact a banded system with p = q = 2 and I can solve the same problem in five
Let’s assume we have a computer that can solve a 100 × 100 system in 1 second using Gauss elimination. Estimate the time required to solve the following problems on the same computer. You must
What are the magnitudes of the eigenvalues of 1 A=[27] 1
What are the criteria for the number of iterations such that Jacobi’s method is faster than Gauss elimination for solving a linear problem that is not banded?
Consider the matrix(a) What is the determinant of the matrix A?(b) What is the bandwidth of the matrix A?(c) What is the 1-norm of the matrix A?(d) How many eigenvalues do you expect to have for A?
Answer the following questions about this 12 × 12 matrix:(a) What is the bandwidth?(b) Assuming this system has a solution for a given vector b, are you guaranteed to have convergence using
What numerical method is implemented by the code listed below (be specific). 1 2 3 4 5 6 7 8 9 10 11 function out = problem2_75 n= size (A, 1); for k=1:n-1 for i=k+1:n m=A (i, k)/A (k, k); for
Consider the following program:What system of equations is being solved by this method? What method is used for the solution? Will this program converge to a solution? To save you time, we have
Answer the following questions about this MATLAB program. You can assume that the norm computed by MATLAB is the 1-norm. (We are using the 1-norm to simplify the problem, not because MATLAB actually
Answer the following questions about this program:(a) What mathematical problem does this program solve?(b) What numerical method does this program use?(c) What is the initial condition used to start
Consider the process flow diagram shown in Fig. 2.14. We can write this out as the system of equationsThe first four equations are the mass balances; the next eight equations are the process
In this problem, you are going to investigate the work required to solve a matrix system using Gauss–Seidel. The matrix system is defined by the entriesand the forcing vector entries areYour
A convenient way to visualize the evolution of a two-dimensional system is to make a phase plane plot. We will do more of this type of calculation in Chapter 5. For now, we would like to see how the
We want to explore the accuracy of MATLAB’s solution to the following set of equationsfor values of ϵ from 10−1 to 10−15. For this problem, you are allowed to use the slash command to solve
In this problem, we want to figure out how the rate of convergence to the successive over-relaxation (SOR) solution of a linear system of equations depends on the relaxation parameter for w = 0.1,
We saw that the time t required to solve a banded matrix of size n × n with a bandwidth p = q scales like t ∼ np2. In this problem, we will explore the ability of MATLAB’s internal solver to
Show thatis a linear equation. f(x) = x (3.1.9)
Show thatis a nonlinear function. f(x)=x² (3.1.13)
Consider the Newton’s method calculation of the roots of the nonlinear equation f (x) = x3. If the initial guess is x(0) = 2, what is x(1)?
For the nonlinear functions that follow, provide the Newton iteration formula that you would use to find the roots of f (x) = 0: f(x) f(x) f(x) = 32x5 - 64x +31 e-3x X COS X - X
For the nonlinear functions that follow, provide the Newton iteration formula that you would use to find the roots of f (x) = 0: (a) f(x)=x²-3x - 2 (b) f(x) = xe-t (c) f(x) = (x - 1) In(x)
Construct a fixed-point method to find the root of the functionusing an initial guess x(0) = 1. f(x) = ex-x (3.2.13)
Consider using Newton’s method to find the real root ofwith an initial guess x(0) = 2. What is x(1)? Do you see any potential problem here? f(x)= Inx+2
Use Picard’s method to find the roots offor the initial guesses −0.1, 0.6, 1.99, and 2.01. f(x)=x²-3x+2=0 (3.3.2)
Consider the system of nonlinear equationsA root of this system isShow that the convergence of Newton–Raphson is independent of the value of ϵ for an initial guessfor any real number ϵ. R = xx2 +
Use Newton’s method to find the root of Eq. (3.2.13) using an initial guess x(0) = 1. f(x) = ex-x (3.2.13)
For what values of x(0) do you expect to have a problem with the Newton’s method solution to compute the roots of f(x) = tan x?
Determine the number of steps required to find the root of f (x) = x2 by Newton’s method as a function of the initial guess x0 if the goal is to reach a tolerance |f (x∗)| ≤ ϵ.
This problem involves computing the liquid–vapor P–x–y diagram from the two parameter Margulies model. In this model, the activity coefficients are given byandThe parameters α and β are
Recall the derivation of the Newton iteration formula discussed in the lecture; it was based on a Taylor series expansion of the function f(x) around an estimate of the root x(k−1), truncated after
Write a MATLAB program that uses Newton–Raphson to compute the solution to the following system of equationsfor an initial guess of x = y = z = 0.2. How many iterations are required such that
Consider the simple equation f (x) = x3 − x2. Do the following analysis for this equation:(a) Determine the three roots by hand.(b) Perform the first five iterations by hand for Picard’s method
In this problem, you will use your knowledge of solving systems of nonlinear equations and continuation methods to determine both the composition and the temperature distribution in an ideal staged
We want to explore the ability of Newton’s method to find roots to the equation sin(πx) = 0. You should already know the roots of this equation. Write a MATLAB program that determines the roots of
Consider the system of equationsFind all of the real roots of this system of equations by hand. Write aMATLAB program that demonstrates quadratic convergence towards each of the real roots of the
In this problem, you are going to look at how the numerical accuracy of the derivative affects the rate of convergence in Newton’s method. Consider the calculation of the root ofwith an initial
For a packed bed, the Ergun equation provides a relationship between the pressure drop per unit length and the properties of the bed,where η is the fluid viscosity, v0 is the superficial velocity,
Consider the partial differential equationsubject to the constant concentration boundary condition c(0, t) = 0 and the reaction–diffusion boundary conditionand the initial condition c(x, 0) = 1.
A vapor–liquid mixture of ethylbenzene and toluene has a partial pressure of 250 mm Hg of ethylbenzene and 343 mm Hg of toluene. Write a MATLAB program that will compute the composition of the
The overall goal of this problem is to compute the P–V and P–T equilibrium diagramsfor a single component fluid described by the van derWaals equation of state. Let us recall the key things we
This problem involves the analysis of a counter-current heat exchanger. The design equation for a counter-current heat exchanger iswhere Q is the heat transferred between streams, A is the area of
Consider non-ideal liquid phases where the activity coefficient is given by the one-parameter modelThe saturation pressure of species 1 is 800 mm Hg, and the saturation pressure of species 2 is 1000
Let’s look at another way to compute phase diagrams from an equation of state without using departure functions but still requiring a numerical solution. The van der Waals equation of state can be
Consider the implicit Euler integration of the unsteady mass balance equations with time step h corresponding to the batch reactionwhich we can write aswith a reaction rate k. What are the equations
Consider the nonlinear ordinary differential equationwith y(2) = 2. Estimate y(2.5)(a) Using the explicit Euler’s method with h = 0.1(b) Using the fourth-order Runge–Kutta method with h = 0.5. dy
Write a program to solve the system of equationsusing implicit Euler. The initial conditions are y1(0) = 1 and y2(0) = 2. dy1 = yiy dt dy2 dt = (y₁ - y2)² (4.6.43) (4.6.44)

Showing 100 - 200 of 376