Consider a fluid bounded by two parallel plates extended to infinity such that no end effects are encountered. The walls and the fluid are initially at rest. The lower wall is suddenly accelerated in the x-direction. The coordinate system is such that the lower wall coincides with the xz-plane to which the y-axis is perpendicular. The spacing between the two plates is denoted by h.
The equation for this problem is
Wheren is the kinematic viscosity of the fluid. It is required to compute the velocity profile u = u(y,t). The problem has the following initial and boundary conditions
Initial Condition:
Boundary Condition:
The fluid has a kinematic viscosity of 0.000217 m2/s and h = 40 mm. Your solution should be run until t = 1.08 seconds. Various time steps are to be used to investigate the numerical schemes and the effect of the time step on stability and accuracy.
Write a code (or multiple codes) to solve this problem for the following schemes and conditions.
a) FTCS Explicit with time steps 0.002 and 0.00233 sec
c) FTCS Implicit with time steps 0.002 and 0.01 sec
If j = 1 at the lower surface and a step size of Dy = 0.001 m is used, then j at the upper surface will be 41. Note that n = 1 corresponds to t = 0.
For each case of each method, plot the velocity profile for four difference times on the same graph: 0.27 sec, 0.54 sec, 0.81sec and 1.08 sec. (You will have 4 plots in all, each with four curves.) Remember to use proper plotting techniques, including title, axis labels, units, legend, and legibility. Submit your plots in hardcopy, and a copy of your code both in hardcopy and electronically on Canvas.