Modelling and Simulation Mid
Modelling and Simulation Mid
Please answer all questions. Use calculators to compute desired values wherever needed. You may use MatLab
software to perform necessary computing. Should you use MatLab, you must report the commands given, inputs
supplied, and outputs obtained.
1. (20 points) Consider a square plate ABCD as shown in the figure below. The plate has been discretized
according to the grid drawn in the figure. The temperatures on the edges are maintained at 30 oC, 25oC, 20oC,
and 15oC respectively on AC, CD, DB, and BA. (Ignore the discontinuity at the corner points.) Assume that the
temperature at a gridpoint inside the boundaries is an average of the temperatures at the the nearest grid-points.
Find the temperature values at the internal gridpoints (from point E to point T).
2. (20 points) Consider a gas obeying the following van der Waals equation of state (relating pressure P,
temperature T, and specific volume V) :
P
a
v2
v−b =RT
where a = 3.592 atm liter2/mol2, and b= 0.04267 liter/mol. R is 0.082054 liter atm /(Mol.K). Using Newton
Raphson method, find out the specific volume of the gas at T = 315.6 K and P = 12 atm. Show your work.
3. (25 points) The fully developed velocity field in a pressure driven flow between two plates moving at differnt
velocities as shown in the figure below, is governed by the following differential equation:
2
P d v
= 2
x dy
P
where v is the velocity along x direction, the pressure gradient along x direction is a constant, and μ is
x
the viscosity of the fluid. The plate velocities are v(y = 0) = 0 at the bottom plate and v(y = 0.1 m) = 1 cm/s at
P
the top plate. The pressure gradient is = 1.5 Pa/m, and the viscosity is μ = 0.1 Poise. By applying finite
x
difference approximation to the derivatives, set up the matrix equation whose solution gives the velocity profile
v(y). Obtain the velocity profile by discretizing the y axis with Δy = 0.01 m. Assume no slip boundary
condition.
4. (25 points) Consider reaction and one dimensional diffusion of a gas species A in a thin rectangular slab of a
catalyst of thickness t = 1 cm, as shown below in the figure.
a) Write the steady state mass balance equation for the diffusing species A inside the slab.
t
b) At the solid-gas interfaces located at x=± , the concentrations are maintained as follows: CA(x = t/2) =
2
5 mM, and CA(x = -t/2) = 1 mM. Values of diffusion coefficient and rate constant are respectively D =
5⨯10-2 cm2/s and k = 1 min-1. Divide the slab in six sections, and numerically obtain the concentration
values at the intermediate grid points.
Suppose you want to solve it by assuming that the approximate solution has the following polynomial form:
2 3 4 5
Y x=a0 a 1 xa2 x a 3 x a4 x a5 x , where the coefficients ais are yet unknon.
a) By requiring that the presumed solution satisfies boundary conditions at the terminal points, and the given
ODE at the terminal and intermediate points, set up the linear system in terms of the unknown coefficients
ais.
b) Solve this linear system and report the y values at the intermediate points.