University:
Widener University
Course:
ESSC 103 | Intriduction To Earth Science
Academic year:

2023
Views:

199

Pages:

1
Author:

Caitlin Mann

Problem set 5 1) Produce the Boussinesq continuity, momentum, and energy equations by starting from the full versions of the equations we derived in class. The purpose of the problem is to motivate you to think critically about scaling arguments, so please justify each simplification with a scaling argument, and comment on the quality of each assumption. 2) Solve the Dirchlet problem, ∇ 2φ = 0 , for a square basin using Matlab and the following boundary conditions: a) Top boundary: Bottom boundary: Left boundary: Right boundary: bctop[1:n]=1 bcbottom[1:n]=2 bcleft[1:n]=3 bcright[1:n]=1 b) Top boundary: Bottom boundary: Left boundary: Right boundary: bctop[1:n]=cos((1:n)/(2*pi)+(t/pi)) bcbottom[1:n]=sin((1:n)/(2*pi)+(t/pi)) bcleft[1:n]=cos((1:n)/(2*pi)+(t/pi)) bcright[1:n]=sin((1:n)/(2*pi)+(t/pi)) For t=1:9 You should make your code general enough that it will support any value of n, but for the work you hand in, let n=20 (hint: KC04 has an entire chapter on potential flow). Turn in a printout of your code (complete with plenty of (relevant) comments), a plot of the potential (see pcolor in Matlab), and a plot of velocity (see quiver in Matlab) for each of the 10 cases given above. For part (b), put the nine plots of potential on a single page and the nine plots of velocity on a single page (see subplot and “orient tall” in Matlab) 3) A hemispherical vessel of radius R has a small rounded orifice of area A at the bottom. Show that the time required to lower the level from h1 to h2 is given by 2π ⎡ 2 1 ⎤ t= R(h13 / 2 − h23 / 2 ) − (h15 / 2 − h25 / 2 )⎥. ⎢ 5 A 2g ⎣ 3 ⎦ [KC04: Ex 4.10]