• Home
  • Textbooks
  • Fluid Mechanics: Fundamentals and Applications
  • INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS

Fluid Mechanics: Fundamentals and Applications

Yunus Cengel

Chapter 15

INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS - all with Video Answers

Educators


Chapter Questions

Problem 1

A CFD code is used to solve a two-dimensional ( $x$ and $y$ ), incompressible, laminar flow without free surfaces. The fluid is Newtonian. Appropriate boundary conditions are used. List the variables (unknowns) in the problem, and list the corresponding equations to be solved by the computer.

Check back soon!

Problem 2

Write a brief (a few sentences) definition and description of each of the following, and provide example(s) if helpful: (a) computational domain, (b) mesh, (c) transport equation, (d) coupled equations.

Check back soon!

Problem 3

What is the difference between a node and an interval and how are they related to cells? In Fig. P15-3C, how many nodes and how many intervals are on each edge?

Check back soon!

Problem 4

For the two-dimensional computational domain of Fig. P15-3C, with the given node distribution, sketch a simple structured grid using four-sided cells and sketch a simple unstructured grid using three-sided cells. How many cells are in each? Discuss.

Check back soon!

Problem 5

Summarize the eight steps involved in a typical CFD analysis of a steady, laminar flow field.

Check back soon!

Problem 6

Suppose you are using CFD to simulate flow through a duct in which there is a circular cylinder as in Fig. P15-6C. The duct is long, but to save computer resources you choose a computational domain in the vicinity of the cylinder only. Explain why the downstream edge of the computational domain should be further from the cylinder than the upstream edge.

Check back soon!

Problem 7

Write a brief (a few sentences) discussion about the significance of each of the following in regards to an iterative CFD solution: (a) initial conditions, (b) residual, (c) iteration, (d) postprocessing.

Check back soon!

Problem 8

Briefly discuss how each of the following is used by CFD codes to speed up the iteration process: (a) multigridding and (b) artificial time.

Check back soon!

Problem 9

Of the boundary conditions discussed in this chapter, list all the boundary conditions that may be applied to the right edge of the two-dimensional computational domain sketched in Fig. P15-9C. Why can't the other boundary conditions be applied to this edge?

Check back soon!

Problem 10

What is the standard method to test for adequate grid resolution when using CFD?

Check back soon!

Problem 11

What is the difference between a pressure inlet and a velocity inlet boundary condition? Explain why you cannot specify both pressure and velocity at a velocity inlet boundary condition or at a pressure inlet boundary condition.

Check back soon!

Problem 12

An incompressible CFD code is used to simulate the flow of air through a two-dimensional rectangular channel (Fig. P15-12C). The computational domain consists of four blocks, as indicated. Flow enters block 4 from the upper right and exits block 1 to the left as shown. Inlet velocity $V$ is known and outlet pressure $P_{\text {out }}$ is also known. Label the boundary conditions that should be applied to every edge of every block of this computational domain.

Check back soon!

Problem 13

Consider Prob. 15-12C again, except let the boundary condition on the common edge between blocks 1 and 2 be a fan with a specified pressure rise from right to left across the fan. Suppose an incompressible CFD code is run for both cases (with and without the fan). All else being equal, will the pressure at the inlet increase or decrease? Why? What will happen to the velocity at the outlet? Explain.

Check back soon!

Problem 14

List six boundary conditions that are used with CFD to solve incompressible fluid flow problems. For each one, provide a brief description and give an example of how that boundary condition is used.

Check back soon!

Problem 15

A CFD code is used to simulate flow over a twodimensional airfoil at an angle of attack. A portion of the computational domain near the airfoil is outlined in Fig. P15-15 (the computational domain extends well beyond the region outlined by the dashed line). Sketch a coarse structured grid using four-sided cells and sketch a coarse unstructured grid using three-sided cells in the region shown. Be sure to cluster the cells where appropriate. Discuss the advantages and disadvantages of each grid type.

Check back soon!

Problem 16

For the airfoil of Prob. 15-15, sketch a coarse hybrid grid and explain the advantages of such a grid.

Check back soon!

Problem 17

An incompressible CFD code is used to simulate the flow of water through a two-dimensional rectangular channel in which there is a circular cylinder (Fig. P15-17). A timeaveraged turbulent flow solution is generated using a turbulence model. Top-bottom symmetry about the cylinder is assumed. Flow enters from the left and exits to the right as shown. Inlet velocity $V$ is known, and outlet pressure $P_{\text {out }}$ is also known. Generate the blocking for a structured grid using
four-sided blocks, and sketch a coarse grid using four-sided cells, being sure to cluster cells near walls. Also be careful to avoid highly skewed cells. Label the boundary conditions that should be applied to every edge of every block of your computational domain. (Hint: Six to seven blocks are sufficient.)

Check back soon!

Problem 18

An incompressible CFD code is used to simulate the flow of gasoline through a two-dimensional rectangular channel in which there is a large circular settling chamber (Fig. P15-18). Flow enters from the left and exits to the right as shown. A time-averaged turbulent flow solution is generated using a turbulence model. Top-bottom symmetry is assumed. Inlet velocity $V$ is known, and outlet pressure $P_{\text {out }}$ is also known. Generate the blocking for a structured grid using four-sided blocks, and sketch a coarse grid using four-sided cells, being sure to cluster cells near walls. Also be careful to avoid highly skewed cells. Label the boundary conditions that should be applied to every edge of every block of your computational domain.

Check back soon!

Problem 19

Redraw the structured multiblock grid of Fig. $15-12 b$ for the case in which your CFD code can handle only elementary blocks. Renumber all the blocks and indicate how many $i$ - and $j$-intervals are contained in each block. How many elementary blocks do you end up with? Add up all the cells, and verify that the total number of cells does not change.

Check back soon!

Problem 20

Suppose your CFD code can handle nonelementary blocks. Combine as many blocks of Fig. 15-12 $b$ as you can. The only restriction is that in any one block, the number of $i$ intervals and the number of $j$-intervals must be constants. Show that you can create a structured grid with only three nonelementary blocks. Renumber all the blocks and indicate how many $i$ - and $j$-intervals are contained in each block. Add up all the cells and verify that the total number of cells does not change.

Check back soon!

Problem 21

A new heat exchanger is being designed with the goal of mixing the fluid downstream of each stage as thoroughly as possible. Anita comes up with a design whose cross section for one stage is sketched in Fig. P15-21. The geometry extends periodically up and down beyond the region shown here. She uses several dozen rectangular tubes inclined at a high angle of attack to ensure that the flow separates and mixes in the wakes. The performance of this geometry is to be tested using two-dimensional time-averaged CFD simulations with a turbulence model, and the results will be compared to those of competing geometries. Sketch the simplest possible computational domain that can be used to simulate this flow. Label and indicate all boundary conditions on your diagram. Discuss.

Check back soon!

Problem 22

Sketch a coarse structured multiblock grid with foursided elementary blocks and four-sided cells for the computational domain of Prob. 15-21.

Check back soon!

Problem 23

Anita runs a CFD code using the computational domain and grid developed in Probs. 15-21 and 15-22. Unfortunately, the CFD code has a difficult time converging and Anita realizes that there is reverse flow at the outlet (far right edge of the computational domain). Explain why there is reverse flow, and discuss what Anita should do to correct the problem.

Check back soon!

Problem 24

As a follow-up to the heat exchanger design of Prob. 15-21, suppose Anita's design is chosen based on the results
of a preliminary single-stage CFD analysis. Now she is asked to simulate $t$ wo stages of the heat exchanger. The second row of rectangular tubes is staggered and inclined oppositely to that of the first row to promote mixing (Fig. P15-24). The geometry extends periodically up and down beyond the region shown here. Sketch a computational domain that can be used to simulate this flow. Label and indicate all boundary conditions on your diagram. Discuss.

Check back soon!

Problem 25

Sketch a structured multiblock grid with four-sided elementary blocks for the computational domain of Prob. $15-24$. Each block is to have four-sided structured cells, but you do not have to sketch the grid, just the block topology. Try to make all the blocks as rectangular as possible to avoid highly skewed cells in the corners. Assume that the CFD code requires that the node distribution on periodic pairs of edges be identical (the two edges of a periodic pair are "linked" in the grid generation process). Also assume that the CFD code does not allow a block's edges to be split for application of boundary conditions.

Check back soon!

Problem 26

In this exercise, we examine how far away the boundary of the computational domain needs to be when simulating external flow around a body in a free stream. We choose a two-dimensional case for simplicity-flow at speed $V$ over a rectangular block whose length $L$ is 1.5 times its height $D$ (Fig. P15-26a). We assume the flow to be symmetric about the centerline ( $x$-axis), so that we need to model only the upper half of the flow. We set up a semicircular computational domain for the CFD solution, as sketched in Fig. P15-26b. Boundary conditions are shown on all edges. We run several values of outer edge radius $R(5<R / D<$ 500) to determine when the far field boundary is "far enough" away. Run FlowLab, and start template Block_domain.
(a) Calculate the Reynolds number based on the block height D. What is the experimentally measured value of the drag coefficient for this two-dimensional block at this Reynolds number (see Chap. 11)?
(b) Generate CFD solutions for various values of $R / D$. For each case, calculate and record drag coefficient $C_D$. Plot $C_D$ as a function of $R / D$. At what value of $R / D$ does $C_D$ become independent of computational extent to three significant digits of precision? Report a final value of $C_D$, and discuss your results.
(c) Discuss some reasons for the discrepancy between the experimental value of $C_D$ and the value obtained here using CFD.
(d) Plot streamlines for two cases: $R / D=5$ and 500 . Compare and discuss.

Check back soon!

Problem 27

Using the geometry of Prob. 15-26, and the case with $R / D=500$, the goal of this exercise is to check for grid independence. Run FlowLab, and start template Block_mesh. Run various values of grid resolution, and tabulate drag coefficient $C_D$ as a function of the number of cells. Has grid independence been achieved? Report a final value of $C_D$ to three significant digits of precision. Does the final value of drag coefficient agree better with that of this experiment? Discuss.

Check back soon!

Problem 28

In Probs. 15-26 and 15-27, we used air as the fluid in our calculations. In this exercise, we repeat the calculation of drag coefficient, except we use different fluids. We adjust the inlet velocity appropriately such that the calculations are always at the same Reynolds number. Run FlowLab, and start template Block_fluid. Compare the value of $C_D$ for all three cases (air, water, and kerosene) and discuss.

Check back soon!

Problem 29

Experiments on two-dimensional rectangular blocks in an incompressible free-stream flow reveal that the drag coefficient is independent of Reynolds number for Re greater than about $10^4$. In this exercise, we examine if CFD calculations are able to predict the same independence of $C_D$ on Re. Run FlowLab, and start template Block_Reynolds. Calculate and record $C_D$ for several values of Re. Discuss.

Check back soon!

Problem 30

In Probs. 15-26 through 15-29, the $k-\varepsilon$ turbulence model is used. The goal of this exercise is to see how sensitive the drag coefficient is to our choice of turbulence model and to see if a different turbulence model yields better agreement with experiment. Run FlowLab, and start template Block_turbulence_model. Run the simulation with all the available turbulence models. For each case, record $C_D$. Which one gives the best agreement with experiment? Discuss.

Check back soon!

Problem 31

Experimental drag coefficient data are available for two-dimensional blocks of various shapes in external flow. In this exercise, we use CFD to compare the drag coefficient of
rectangular blocks with $L D D$ ranging from 0.1 to 3.0 (Fig. P15-31). The computational domain is a semicircle similar to that sketched in Fig. P15-26b; we assume steady, incompressible, turbulent flow with symmetry about the $x$-axis. Run FlowLab, and start template Block_length.
(a) Run the CFD simulation for various values of $L / D$ between 0.1 and 3.0. Record the drag coefficient for each case, and plot $C_D$ as a function of $L D D$. Compare to experimentally obtained data on the same plot. Discuss.
(b) For each case, plot streamlines near the block and in its wake region. Use these streamlines to help explain the trend in the plot of $C_D$ versus $L / D$.
(c) Discuss possible reasons for the discrepancy between CFD calculations and experimental data and suggest a remedy.

Check back soon!

Problem 32

Repeat Prob. 15-26 for the case of axisymmetric flow over a blunt-faced cylinder (Fig. P15-32), using FlowLab template Block_axisymmetric. The grids and all the parameters are the same as those in Prob. 15-26, except the symmetry boundary condition is changed to "axis," and the flow solver is axisymmetric about the $x$-axis. In addition to the questions listed in Prob. 15-26, compare the two-dimensional and axisymmetric cases. Which one requires a greater extent of the far field boundary? Which one has better agreement with experiment? Discuss. (Note: The reference area for $C_D$ in the axisymmetric case is the frontal area $A=\pi D^2 / 4$.)

Check back soon!

Problem 33

Air flows through a conical diffuser in an axisymmetric wind tunnel (Fig. P15-33a-drawing not to scale). $\theta$ is the diffuser half-angle (the total angle of the diffuser is equal to $2 \theta$ ). The inlet and outlet diameters are $D_1=0.50 \mathrm{~m}$ and $D_2=1.0 \mathrm{~m}$, respectively, and $\theta=20^{\circ}$. The inlet velocity is nearly uniform at $V=10.0 \mathrm{~m} / \mathrm{s}$. The axial distance upstream of the diffuser is $L_1=1.50 \mathrm{~m}$, and the axial distance from the start of the diffuser to the outlet is $L_2$ $=8.00 \mathrm{~m}$. We set up a computational domain for a CFD solution, as sketched in Fig. P15-33b. Since the flow is axisymmetric and steady in the mean, we model only one two-dimensional slice as shown, with the bottom edge of the domain specified as an axis. The goal of this exercise is to test for grid independence. Run FlowLab, and start template Diffuser_mesh.
(a) Generate CFD solutions for several grid resolutions. Plot streamlines in the diffuser section for each case. At what grid resolution does the streamline pattern appear to be grid independent? Describe the flow field for each case and discuss.
(b) For each case, calculate and record pressure difference $\Delta P=P_{\text {in }}-P_{\text {out }}$. At what grid resolution is the $\Delta P$ grid independent (to three significant digits of precision)? Plot $\Delta P$ as a function of number of cells. Discuss your results.

Check back soon!

Problem 34

Repeat Prob. 15-33 for the finest resolution case, but with the "pressure outlet" boundary condition changed to an "outflow" boundary condition instead, using FlowLab template Diffuser_outflow. Record $\Delta P$ and compare with the result of Prob. 15-33 for the same grid resolution. Also compare the pressure distribution at the outlet for the case with the pressure outlet boundary condition and the case with the outflow boundary condition. Discuss.

Check back soon!

Problem 35

Barbara is designing a conical diffuser for the axisymmetric wind tunnel of Prob. 15-33. She needs to achieve at least 40 Pa of pressure recovery through the diffuser, while keeping the diffuser length as small as possible. Barbara decides to use CFD to compare the performance of diffusers of various half-angles ( $5^{\circ} \leq \theta \leq 90^{\circ}$ ) (see Fig. P15-33 for the definition of $\theta$ and other parameters in the problem). In all cases, the diameter doubles through the dif-fuser-the inlet and outlet diameters are $D_1=0.50 \mathrm{~m}$ and $D_2$ $=1.0 \mathrm{~m}$, respectively. The inlet velocity is nearly uniform at $V=10.0 \mathrm{~m} / \mathrm{s}$. The axial distance upstream of the diffuser is $L_1=1.50 \mathrm{~m}$, and the axial distance from the start of the diffuser to the outlet is $L_2=8.00 \mathrm{~m}$. (The overall length of the computational domain is 9.50 m in all cases.)

Run FlowLab with template Diffuser_angle. In addition to the axis and wall boundary conditions labeled in Fig. P15-33, the inlet is specified as a velocity inlet and the outlet is specified as a pressure outlet with $P_{\mathrm{cot}}=0$ gage pressure for all cases. The fluid is air at default conditions, and turbulent flow is assumed.
(a) Generate CFD solutions for half-angle $\theta=5,7.5,10$, $12.5,15,17.5,20,25,30,45,60$, and $90^{\circ}$. Plot streamlines for each case. Describe how the flow field changes with the diffuser half-angle, paying particular attention to flow separation on the diffuser wall. How small must $\theta$ be to avoid flow separation?
(b) For each case, calculate and record $\Delta P=P_{\text {in }}-P_{\text {out }}$. Plot $\Delta P$ as a function of $\theta$ and discuss your results. What is the maximum value of $\theta$ that achieves Barbara's design objectives?

Check back soon!

Problem 36

Consider the diffuser of Prob. 15-35 with $\theta=90^{\circ}$ (sudden expansion). In this exercise, we test whether the grid is fine enough by performing a grid independence check. Run FlowLab, and start template Expansion_mesh. Run the CFD code for several levels of grid refinement. Calculate and record $\Delta P$ for each case. Discuss.

Check back soon!

Problem 37

Water flows through a sudden contraction in a small round tube (Fig. P15-37a). The tube diameters are $D_1$ $=8.0 \mathrm{~mm}$ and $D_2=2.0 \mathrm{~mm}$. The inlet velocity is nearly uniform at $V=0.050 \mathrm{~m} / \mathrm{s}$, and the flow is laminar. Shane wants to predict the pressure difference from the inlet $\left(x=-L_1\right)$ to the axial location of the sudden contraction $(x=0)$. He sets up the computational domain sketched in Fig. P15-37b. Since the flow is axisymmetric and steady, Shane models only one slice, as shown, with the bottom edge of the domain specified as an axis. In addition to the boundary conditions labeled in Fig. P15-37b, the inlet is specified as a velocity inlet, and the outlet is specified as a pressure outlet with $P_{\text {out }}=0$ gage. What Shane does not know is how far he needs to extend the domain downstream of the contraction in order for the flow field to be simulated accurately upstream of the contraction. (He has no interest in the flow downstream of the contraction.) In other words, he does not know how long to make $L_{\text {extend }}$. Run FlowLab, and start template Contraction_domain.
(a) Generate solutions for $L_{\text {eutend }} / D_2=0.25,0.5,0.75,1.0$, $1.25,1.5,2.0,2.5$, and 3.0. How big must $L_{\text {extend }} / D_2$ be in order to avoid reverse flow at the pressure outlet? Explain. Plot streamlines near the sudden contraction to help explain your results.
(b) For each case, record gage pressures $P_{\text {in }}$ and $P_1$, and calculate $\Delta P=P_{\text {in }}-P_1$. How big must $L_{\text {extend }} / D_2$ be in order for $\Delta P$ to become independent of $L_{\text {extend }}$ (to three significant digits of precision)?
(c) Plot inlet gage pressure $P_{\text {in }}$ as a function of $L_{\text {eutend }} / D_2$. Discuss and explain the trend. Based on all your results taken collectively, which value of $L_{\text {extend }} / D_2$ would you recommend to Shane?

Check back soon!

Problem 38

Consider the sudden contraction of Prob. 15-37 (Fig. P15-37). Suppose Shane were to disregard the downstream extension entirely ( $L_{\text {extend }} / D_2=0$ ). Run FlowLab, and start template Contraction_zerolength. Iterate to convergence. Is there reverse flow? Explain. Plot streamlines near the outlet, and compare with those of Prob. 15-37. Discuss. Calculate $\Delta P=P_{\text {in }}-P_{\text {out }}$, and calculate the percentage error in $\Delta P$ under these conditions, compared to the converged value of Prob. 15-37. Discuss.

Check back soon!

Problem 39

In this exercise, we apply different back pressures to the sudden contraction of Prob. 15-37 (Fig. P15-37), for the case with $L_{\text {extend }} / D_2=2.0$. Run FlowLab, and start template Contraction_pressure. Set the pressure boundary condition at
the outlet to $P_{\text {out }}=-50,000 \mathrm{~Pa}$ gage (about $1 / 2 \mathrm{~atm}$ below atmospheric pressure). Record $P_{\text {in }}$ and $P_1$, and calculate $\Delta P$ $=P_{\text {in }}-P_1$. Repeat for $P_{\text {out }}=0 \mathrm{~Pa}$ gage and $P_{\text {oat }}=50,000$ Pa gage. Discuss your results.

Check back soon!

Problem 40

Consider the sudden contraction of Prob. 15-37, but this time with turbulent rather than laminar flow. The dimensions shown in Fig. P15-37 are scaled proportionally by a factor of 100 everywhere so that $D_1=0.80 \mathrm{~m}$ and $D_2$ $=0.20 \mathrm{~m}$. The inlet velocity is also increased to $V=1.0 \mathrm{~m} / \mathrm{s}$. A 10 percent turbulence intensity is specified at the inlet. The outlet pressure is fixed at zero gage pressure for all cases. Run FlowLab, and start template Contraction_turbulent.
(a) Calculate the Reynolds numbers of flow through the large tube and the small tube for Prob. 15-37 and also for this problem. Are our assumptions of laminar versus turbulent flow reasonable for these problems?
(b) Generate CFD solutions for $L_{\text {ettend }} / D_2=0.25,0.5,0.75$, $1.0,1.25,1.5$, and 2.0 . How big must $L_{\text {extend }} / D_2$ be in order to avoid reverse flow at the pressure outlet? Plot streamlines for the case in which $L_{\text {extend }} / D_2=0.75$ and compare to the corresponding streamlines of Prob. 15-37 (laminar flow). Discuss.
(c) For each case, record gage pressures $P_{\text {in }}$ and $P_1$, and calculate $\Delta P=P_{\text {in }}-P_1$. How big must $L_{\text {extend }} / D_2$ be in order for $\Delta P$ to become independent of $L_{\text {extend }}$ (to three significant digits of precision)?

Check back soon!

Problem 41

Run FlowLab, and start template Contraction_outflow. The conditions are identical to Prob. 15-40 for the case with $L_{\text {extend }} / D_2=0.75$, but with the "pressure outlet" boundary condition changed to an "outflow" boundary condition instead. Record $P_{\text {in }}$ and $P_1$, calculate $\Delta P=P_{\text {in }}-P_1$, and compare with the result of Prob. 15-40 for the same geometry. Discuss.

Check back soon!

Problem 42

Run FlowLab with template Contraction_2d. This is identical to the sudden contraction of Prob. 15-40, but the flow is two-dimensional instead of axisymmetric. (Note that the "axis" boundary condition is replaced by "symmetry.") As previously, the outlet pressure is set to zero gage pressure.
(a) Generate CFD solutions for $L_{\text {eatend }} / D_2=0.25,0.5,0.75$, $1.0,1.25,1.5,2.0,3.0$, and 4.0. How big must $L_{\text {extend }} / D_2$ be in order to avoid reverse flow at the pressure outlet? Plot streamlines for the case in which $L_{\text {extend }} / D_2=0.75$, and compare to the corresponding streamlines of Prob. 15-40 (axisymmetric flow). Discuss.
(b) For each case, record gage pressures $P_{\text {in }}$ and $P_1$, and calculate $\Delta P=P_{\text {in }}-P_1$. How big must $L_{\text {extend }} / D_2$ be in order for $\Delta P$ to become independent of $L_{\text {extend }}$ (to three significant digits of precision)?

Check back soon!

Problem 43

Air flows through a "jog" in a rectangular channel (Fig. P15-43a, not to scale). The channel dimension is $D_1$ $=1.0 \mathrm{~m}$ everywhere, and it is wide enough (into the page of Fig. P15-43) that the flow can be considered two-dimensional. The inlet velocity is nearly uniform at $V=1.0 \mathrm{~m} / \mathrm{s}$. The distance upstream of the jog is $L_1=5.0 \mathrm{~m}$, the overall jog length is $L_j=3.0 \mathrm{~m}$, and the distance from the end of the jog to the outlet is $L_2=10.0 \mathrm{~m}$. We set up a computational domain for a CFD solution, as sketched in Fig. P15-43b. In addition to the wall boundary conditions labeled in Fig. P15-43b, the inlet is specified as a velocity inlet and the outlet is specified as a pressure outlet with $P_{\text {out }}=0$ gage pressure. The fluid is air at default conditions, and turbulent flow is assumed. The goal of this exercise is to test for grid independence in this flow field. Run FlowLab with template Jog_turbulent_mesh.
(a) Generate CFD solutions for various levels of grid resolution. Plot streamlines in the region of the jog for each case. At what grid resolution does the streamline pattern appear to be grid independent? Discuss.
(b) For each case, calculate and record $\Delta P=P_{\text {in }}-P_{\text {out }}$. At what grid resolution is $\Delta P$ grid independent (to three significant digits of precision)? Plot $\Delta P$ as a function of the number of cells. Discuss your results.

Check back soon!

Problem 44

Repeat Prob. 15-43, but for laminar flow, using Jog_laminar_mesh as the FlowLab template. The jog is identical in shape, but scaled down by a factor of 1000 compared to that of Prob. 15-43 (the channel width is $D_1=1.0 \mathrm{~mm}$ everywhere). The inlet velocity is nearly uniform at $V$ $=0.10 \mathrm{~m} / \mathrm{s}$, and the fluid is changed to water at room temperature. Discuss your results.

Check back soon!

Problem 45

Repeat Prob. 15-44, but for laminar flow at a higher Reynolds number, using FlowLab template Jog_high_Re.

Everything is identical to Prob. 15-44, except the inlet velocity is increased from to $V=0.10$ to $1.0 \mathrm{~m} / \mathrm{s}$. Compare results and the Reynolds numbers for the two cases and discuss.

Check back soon!

Problem 46

Consider compressible flow of air through an axisymmetric converging-diverging nozzle (Fig. P15-46), in which the inviscid flow approximation is applied. The inlet conditions are fixed ( $P_{0, \text { inlet }}=220 \mathrm{kPa}, P_{\text {inlet }}=210 \mathrm{kPa}$, and $T_{0, \text { inlet }}=300 \mathrm{~K}$ ), but the back pressure $P_b$ can be varied. Run FlowLab, using template Nozzle_axisymmetric. Do several cases, with back pressure ranging from 100 to 219 kPa . For each case, calculate the mass flow rate ( $\mathrm{kg} / \mathrm{s}$ ) through the nozzle, and plot $\dot{m}$ as a function of $P_b / P_{0, \text { inlet }}$. Explain your results.

Check back soon!

Problem 47

Run FlowLab with template Nozzle_axisymmetric (Prob. 15-46). For the case in which $P_b=100 \mathrm{kPa}\left(P_b / P_{0, \text { inlet }}\right.$ $=0.455$ ), plot pressure and the Mach number contours to verify that a normal shock is present near the outlet of the computational domain. Generate a plot of average Mach number Ma and average pressure ratio $P / P_{0, \text { inlet }}$ across several cross sections of the domain, as in Fig. 15-76. Point out the location of the normal shock, and compare the CFD results to one-dimensional compressible flow theory. Repeat for $P_b$ $=215 \mathrm{kPa}\left(P_b / P_{0, \text { inlet }}=0.977\right)$. Explain.

Check back soon!

Problem 48

Run FlowLab with template Nozzle_2d, which is the same as Prob. 15-46, except the flow is two-dimensional instead of axisymmetric. Note that the "axis" boundary condition is also changed to "symmetry." Compare your results and discuss the similarities and differences.

Check back soon!

Problem 49

Consider flow over a simplified, two-dimensional model of an automobile (Fig. P15-49). The inlet conditions are fixed at $V=60.0 \mathrm{mi} / \mathrm{h}(26.8 \mathrm{~m} / \mathrm{s})$, with 10 percent turbulence intensity. The standard $k-\varepsilon$ turbulence model is used. Run FlowLab with template Automobile_drag. Vary the shape of the rear end of the car, and record the drag coefficient for each shape. Also plot velocity vectors in the vicinity of the rear end for each case. Compare and discuss. Which case gives the lowest drag coefficient? Why?

Check back soon!

Problem 50

In this exercise, we examine the effect of the location of the upper symmetry boundary condition of Prob. 15-49. Run FlowLab with template Automobile_domain for several values of $H / h$ (Fig. P15-49). Plot the calculated value of $C_D$ as a function of $H / h$. At what value of $H / h$ does $C_D$ level off? In other words, how far away must the upper symmetry boundary be in order to have negligible influence on the calculated value of drag coefficient? Discuss.

Check back soon!

Problem 51

Run FlowLab, and start template Automobile_turbulence_model. In this exercise, we examine the effect of turbulence model on the calculation of drag on a simplified, twodimensional model of a car (Fig. P15-49). Run all the available turbulence models. For each case, record $C_D$. Is there much variation in the calculated values of $C_D$ ? Which one is correct? Discuss.

Check back soon!

Problem 52

Run FlowLab, and start template Automobile_3d. In this exercise, we compare the drag coefficient for a fully three-dimensional automobile to that predicted by the twodimensional approximation of Prob. 15-49. Note that the solution takes a long time to converge and requires a significant amount of computer resources. Therefore, the converged solution is already available in this template. Observe the three-dimensional pathlines around the car by rotating the view. Calculate the drag coefficient. Is it larger or smaller than the two-dimensional prediction? Discuss.

Check back soon!

Problem 53

Run FlowLab, and start template Pipe_laminar_ developing. In this exercise, we study laminar flow in the entrance region of a round pipe (Fig. P15-53, not to scale). Because of the axisymmetry, the computational domain consists of one slice (light blue region in Fig. P15-53). Calculate
the flow at various values of the Reynolds number Re, where Re is based on pipe diameter and average speed through the pipe. For each case, study the velocity profiles at several axial locations down the pipe, and estimate the entrance length in each case. Also plot the pressure distribution along the pipe axis for each case. Estimate the end of the entrance region as the location where the pressure begins to drop linearly with $x$. Compare your results with those obtained from the velocity profiles, and also with theory, $L_{\ell} / D \cong 0.06 \mathrm{Re}$. Discuss.

Check back soon!

Problem 54

Run FlowLab, and start template Pipe_turbulent_ developing. In this exercise, we study turbulent flow in the entrance region of a round pipe (Fig. P15-53). Calculate the flow at several values of the Reynolds number. For each case, study the velocity profiles at several axial locations down the pipe, and estimate the entrance length in each case. Also plot the pressure distribution along the pipe axis for each case. Estimate the end of the entrance region as the location where the pressure begins to drop linearly with $x$. Compare your results with those obtained from the velocity profiles, and also with the empirical approximation, $L_e / D \cong 4.4 \mathrm{Re}^{1 / 6}$. Compare your results to those of the laminar flow of Prob. 15-53. Discuss. Which flow regime, laminar or turbulent, has the longer entrance length? Why?

Check back soon!

Problem 55

Consider fully developed, laminar pipe flow (Fig. P15-55). In this exercise, we are not concerned about entrance effects. Instead, we want to analyze the fully developed flow downstream of the entrance region. Because of the axisymmetry, the computational domain consists of one slice (light blue region). The velocity profile at the inlet boundary is set to be the same as that at the outlet boundary, but a pressure drop from $x=0$ to $L$ is imposed to simulate fully developed flow. Run FlowLab with template Pipe_laminar_developed. The template is set up such that the outlet velocity profile gets fed into the inlet. In other words, the inlet and outlet are periodic boundary conditions, but with an imposed pressure drop. Run several cases corresponding to various values of the Reynolds number. For each case, look at velocity profiles to confirm that the flow is fully developed. Calculate and plot Darcy friction factor $f$ as a function of Re , and compare with the theoretical value for laminar flow, $f=64 / \mathrm{Re}$. Discuss the agreement between CFD and theory.

Check back soon!

Problem 56

Repeat Prob. 15-55, except for fully developed turbulent flow through a smooth-walled pipe. Use FlowLab template Pipe_turbulent_developed. Calculate and plot the Darcy friction factor $f$ as a function of Re. Compare $f$ with that predicted in Chap. 8 for fully developed turbulent pipe flow through a smooth pipe. Discuss.

Check back soon!

Problem 57

In Prob. 15-56, we considered fully developed turbulent flow through a smooth pipe. In this exercise, we examine fully developed turbulent flow through a rough pipe. Run FlowLab with template Pipe_turbulent_rough. Run several cases, each with a different value of normalized pipe roughness, $\varepsilon / D$, but at the same Reynolds number. Calculate and tabulate Darcy friction factor $f$ as a function of normalized roughness parameter $\varepsilon / D$. Compare $f$ with that predicted by the Colebrook equation for fully developed turbulent pipe flow in rough pipes. Discuss.

Check back soon!

Problem 58

Consider the laminar boundary layer developing over a flat plate (Fig. P15-58). Run FlowLab with template Plate_laminar. The inlet velocity and length are chosen such that the Reynolds number at the end of the plate, $\mathrm{Re}_L$ $=\rho V L / \mu$, is approximately $1 \times 10^5$, just on the verge of transition toward turbulence. From your CFD results, calculate the following, and compare to theory: (a) the boundary layer profile shape at $x=L$ (compare to the Blasius profile), (b) boundary layer thickness $\delta$ as a function of $x$, and (c) drag coefficient on the plate.

Check back soon!

Problem 59

Repeat Prob. 15-58, but for turbulent flow on a smooth flat plate. Use FlowLab template Plate_turbulent. The Reynolds number at the end of the plate is approximately $1 \times 10^7$ for this case-well beyond the transition region.

Check back soon!

Problem 60

Run FlowLab, and start template Plate_turbulence_models. In this exercise, we examine the effect of the turbulence model on the calculation of the drag coefficient on a flat plate (Fig. P15-58). Run for each of the available turbulence models. For each case, record $C_D$. Is there much variation in the calculated values of $C_D$ ? Which turbulence model yields the most correct value of drag coefficient? Discuss.

Check back soon!

Problem 61

Consider laminar flow on a smooth heated flat plate (Fig. P15-61). Run FlowLab template Plate_laminar_temperature for two fluids: air and water. The inlet velocity is adjusted such that the Reynolds number for the air and water cases are approximately equal. Compare the 99 percent temperature thickness at the end of the plate to the 99 percent velocity thickness. Discuss your results. (Hint: What is the Prandtl number of air and of water?)

Check back soon!

Problem 62

Repeat Prob. 15-61, except for turbulent flow on a smooth heated flat plate (Fig. P15-61). Use FlowLab template Plate_turbulent_temperature. Discuss the differences between the laminar and turbulent calculations. Specifically, which regime (laminar or turbulent) produces the largest variation between 99 percent temperature thickness and 99 percent velocity thickness? Explain.

Check back soon!

Problem 63

Consider turbulent flow of water through a smooth, $90^{\circ}$, flanged elbow in a round pipe (Fig. P15-63). Because of symmetry, only half of the pipe is modeled; the center plane is specified as a "symmetry" boundary condition. The pipe walls are smooth. The inlet velocity and pipe diameter are chosen to yield a Reynolds number of 20,000 . For the first (default) case, the standard $k-\varepsilon$ turbulence model is used. Run FlowLab with template Elbow. This is a three-dimensional calculation, so expect significantly longer run times. The average pressure is calculated across several cross sections of the pipe: upstream of the elbow, in the elbow, and downstream of the elbow (sections A-A, B-B, etc., in Fig. P15-63). Plot average pressure as a function of axial distance along the pipe. Where does most of the pressure drop occur-in the pipe section upstream of the elbow, in the elbow itself, immediately downstream of the elbow, or in the pipe section downstream of the elbow? Discuss.

Check back soon!

Problem 64

Run FlowLab with template Elbow, again using the standard $k-\varepsilon$ turbulence model. In this exercise, we study velocity vectors in the plane of several cross sections along the pipe. Compare the velocity vectors at a section upstream of the elbow, at a section in the elbow, and at several sections downstream of the elbow. At which locations do you observe counter-rotating eddies? How does the strength of the counter-rotating eddies change with downstream distance? Discuss. Explain why many manufacturers of pipe flowmeters recommend that their flowmeter be installed at least 10 or 20 pipe diameters downstream of an elbow.

Check back soon!

Problem 65

Run FlowLab with template Elbow, again using the standard $k-\varepsilon$ turbulence model. In this exercise, we calculate the minor loss coefficient $K_L$ for the elbow of Prob. 15-63. In order to do so, we compare the pressure drop calculated through the pipe with the elbow to that through a straight pipe of the same overall length, and with identical inlet and outlet conditions. Calculate the pressure drop from inlet to outlet for both geometries. To calculate $K_L$ for the elbow, subtract $\Delta P$ of the straight pipe from $\Delta P$ of the pipe with the elbow. The difference thus represents the pressure drop due to the elbow alone. From this pressure drop and the average velocity through the pipe, calculate minor loss coefficient $K_L$, and compare to the value given in Chap. 8 for a smooth, $90^{\circ}$, flanged elbow.

Check back soon!

Problem 66

In this exercise, we examine the effect of the turbulence model on the calculation of the minor loss coefficient of a pipe elbow (Fig. P15-63). Using FlowLab template Elbow, repeat Prob. 15-65, but with various turbulence models. For each case, calculate $K_L$. Is there much variation in the calculated values of $K_L$ ? Which turbulence model yields the most correct value, compared with the empirical result of Chap. 8? The Spallart-Allmaras model is the simplest, while the Reynolds stress model is the most complicated of the four. Do the calculated results improve with turbulence model complexity? Discuss.

Check back soon!

Problem 67

Consider flow over a two-dimensional airfoil of chord length $L_c$ at an angle of attack $\alpha$ in a flow of freestream speed $V$ with density $\rho$ and viscosity $\mu$. Angle $\alpha$ is measured relative to the free-stream flow direction. (Fig. P15-67). In this exercise, we calculate the nondimensional lift and drag coefficients $C_L$ and $C_D$ that correspond to lift
and drag forces $F_L$ and $F_D$, respectively. Free-stream velocity and chord length are chosen such that the Reynolds number based on $V$ and $L_c$ is $1 \times 10^7$ (turbulent boundary layer over nearly the entire airfoil). Run FlowLab with template Airfoil_angle at several values of $\alpha$, ranging from -2 to $20^{\circ}$. For each case, calculate $C_L$ and $C_D$. Plot $C_L$ and $C_D$ as functions of $\alpha$. At approximately what angle of attack does this airfoil stall?

Check back soon!

Problem 68

In this problem, we study the effect of Reynolds number on the lift and drag coefficients of an airfoil at various angles of attack. Note that the airfoil used here is of a different shape than that used in Problem 15-67. Run FlowLab with template "Airfoil_Reynolds." For the case with Reynolds number equal to $3 \times 10^6$, calculate and plot $C_L$ and $C_D$ as functions of $\alpha$, with $\alpha$ ranging from -2 to $24^{\circ}$. What is the stall angle for this case? Repeat for $\operatorname{Re}=6 \times 10^6$. Compare the two results and discuss the effect of Reynolds number on lift and drag of this airfoil.

Check back soon!

Problem 69

In this exercise, we examine the effect of grid resolution on the calculation of airfoil stall (flow separation) for the airfoil of Problem 15-67 at $\alpha=15^{\circ}$ and $\operatorname{Re}=1 \times 10^7$. Run FlowLab with template Airfoil_mesh. Run for several levels of grid resolution. For each case, calculate $C_L$ and $C_D$. How does grid resolution affect the stall angle? Has grid independence been achieved?

Check back soon!

Problem 70

Consider creeping flow produced by the body of a microorganism swimming through water, represented here as a simple $2 \times 1$ ellipsoid (Fig. P15-70, not to scale). Applied boundary conditions are shown for each edge in parentheses. The flow is laminar, and the default values of $V$ and $L$ are chosen such that the Reynolds number $\operatorname{Re}=\rho V L / \mu$ is equal to 0.20 . Run FlowLab with template Creep_domain. Vary the computational domain radius from $R / L=3$ to 2000 . For each case, calculate the drag coefficient $C_D$ on the body. How large of a computational domain is required for the drag coefficient to level off (far field boundary conditions no longer have significant influence)? Discuss. For the largest computational domain case ( $R / L=2000$ ), plot velocity vectors along a vertical line coincident with the $y$-axis. Compare to the velocity profile we would expect at very high Reynolds numbers. Discuss.

Check back soon!

Problem 71

Run FlowLab with template Creep_Reynolds. In this exercise, the Reynolds number is varied from 0.1 to 100 for flow over an ellipsoid (Fig. P15-70). Plot $C_D$ as a function of Re , and compare velocity profiles along the $y$-axis as Re increases above the creeping flow regime. Discuss.

Check back soon!

Problem 72

Consider the two-dimensional wye of Fig. P15-72. Dimensions are in meters, and the drawing is not to scale. Incompressible flow enters from the left, and splits into two parts. Generate three coarse grids, with identical node distributions on all edges of the computational domain: (a) structured multiblock grid, (b) unstructured triangular grid, and (c) unstructured quadrilateral grid. Compare the number of cells in each case and comment about the quality of the grid in each case.

Check back soon!

Problem 73

Choose one of the grids generated in Prob. 15-72, and run a CFD solution for laminar flow of air with a uniform inlet velocity of $0.02 \mathrm{~m} / \mathrm{s}$. Set the outlet pressure at both outlets to the same value, and calculate the pressure drop through the wye. Also calculate the percentage of the inlet flow that goes out of each branch. Generate a plot of streamlines.

Check back soon!

Problem 74

Repeat Prob. 15-73, except for turbulent flow of air with a uniform inlet velocity of $10.0 \mathrm{~m} / \mathrm{s}$. In addition, set the turbulence intensity at the inlet to 10 percent with a turbulent length scale of 0.5 m . Use the $k-\varepsilon$ turbulence model with wall functions. Set the outlet pressure at both outlets to the same value, and calculate the pressure drop through the wye. Also calculate the percentage of the inlet flow that goes out of each branch. Generate a plot of streamlines. Compare results with those of laminar flow (Prob. 15-73).

Check back soon!

Problem 75

Generate a computational domain to study the laminar boundary layer growing on a flat plate at $\mathrm{Re}=10,000$. Generate a very coarse mesh, and then continually refine the mesh until the solution becomes grid independent. Discuss.

Check back soon!

Problem 76

Repeat Prob. 15-75, except for a turbulent boundary layer at $\operatorname{Re}=10^6$. Discuss.

Check back soon!

Problem 77

Generate a computational domain to study ventilation in a room (Fig. P15-77). Specifically, generate a rectangular room with a velocity inlet in the ceiling to model the supply air, and a pressure outlet in the ceiling to model the return air. You may make a two-dimensional approximation for simplicity (the room is infinitely long in the direction normal to the page in Fig. P15-77). Use a structured rectangular grid. Plot streamlines and velocity vectors. Discuss.

Check back soon!

Problem 78

Repeat Prob. 15-77, except use an unstructured triangular grid, keeping everything else the same. Do you get the same results as those of Prob. 15-77? Compare and discuss.

Check back soon!

Problem 79

Repeat Prob. 15-77, except move the supply and/or return vents to various locations in the ceiling. Compare and discuss.

Check back soon!

Problem 80

Choose one of the room geometries of Probs. 15-77 and $15-79$, and add the energy equation to the calculations. In particular, model a room with air-conditioning, by specifying the supply air as $\operatorname{cool}\left(T=18^{\circ} \mathrm{C}\right)$, while the walls, floor, and ceiling are warm ( $T=26^{\circ} \mathrm{C}$ ). Adjust the supply air speed until the average temperature in the room is as close as possible to $22^{\circ} \mathrm{C}$. How much ventilation (in terms of number of room air volume changes per hour) is required to cool this room to an average temperature of $22^{\circ} \mathrm{C}$ ? Discuss.

Check back soon!

Problem 81

Repeat Prob. 15-80, except create a three-dimensional room, with an air supply and an air return in the ceiling. Compare the two-dimensional results of Prob. 15-80 with the more realistic three-dimensional results of this problem. Discuss.

Check back soon!

Problem 82

Generate a computational domain to study compressible flow of air through a converging nozzle with atmospheric pressure at the nozzle exit (Fig. P15-82). The nozzle walls may be approximated as inviscid (zero shear stress). Run several cases with various values of inlet pressure. How much inlet pressure is required to choke the flow? What happens if the inlet pressure is higher than this value? Discuss.

Check back soon!

Problem 83

Repeat Prob. 15-82, except remove the inviscid flow approximation. Instead, let the flow be turbulent, with smooth, no-slip walls. Compare your results to those of Prob. 15-82. What is the major effect of friction in this problem? Discus.

Check back soon!

Problem 84

Generate a computational domain to study incompressible, laminar flow over a two-dimensional streamlined body (Fig. P15-84). Generate various body shapes, and calculate the drag coefficient for each shape. What is the smallest value of $C_D$ that you can achieve? (Note: For fun, this problem can be turned into a contest between students. Who can generate the lowest-drag body shape?)

Check back soon!

Problem 85

Repeat Prob. 15-84, except for an axisymmetric, rather than a two-dimensional, body. Compare to the twodimensional case. Which has the lower drag coefficient? Discuss.

Check back soon!

Problem 86

Repeat Prob. 15-85, except for turbulent, rather than laminar, flow. Compare to the laminar case. Which has the lower drag coefficient? Discuss.

Check back soon!

Problem 87

Generate a computational domain to study Mach waves in a two-dimensional supersonic channel (Fig. P15-87). Specifically, the domain should consist of a simple rectangular channel with a supersonic inlet $(\mathrm{Ma}=2.0)$, and with a very small bump on the lower wall. Using air with the inviscid flow approximation, generate a Mach wave, as sketched. Measure the Mach angle, and compare with theory (Chap. 12). Also discuss what happens when the Mach wave hits the opposite wall. Does it disappear, or does it reflect, and if so, what is the reflection angle? Discuss.

Check back soon!

Problem 88

Repeat Prob. 15-87, except for several values of the Mach number, ranging from 1.10 to 3.0 . Plot the calculated Mach angle as a function of Mach number and compare to the theoretical Mach angle (Chap. 12). Discuss.

Check back soon!

Problem 89

For each statement, choose whether the statement is true or false, and discuss your answer briefly.
(a) The physical validity of a CFD solution always improves as the grid is refined.
(b) The $x$-component of the Navier-Stokes equation is an example of a transport equation.
(c) For the same number of nodes in a two-dimensional mesh, a structured grid typically has fewer cells than an unstructured triangular grid.
(d) A time-averaged turbulent flow CFD solution is only as good as the turbulence model used in the calculations.

Check back soon!

Problem 90

In Prob. 15-18 we take advantage of top-bottom symmetry when constructing our computational domain and grid. Why can't we also take advantage of the right-left symmetry in this exercise? Repeat the discussion for the case of potential flow.

Check back soon!

Problem 91

Gerry creates the computational domain sketched in Fig. P15-91C to simulate flow through a sudden contraction in a two-dimensional duct. He is interested in the timeaveraged pressure drop (minor loss coefficient) created by the sudden contraction. Gerry generates a grid and calculates the flow with a CFD code, assuming steady, turbulent, incompressible flow (with a turbulence model).
(a) Discuss one way that Gerry could improve his computational domain and grid so that he would get the same results in approximately half the computer time.
(b) There may be a fundamental flaw in how Gerry has set up his computational domain. What is it? Discuss what should be different about Gerry's setup.

Check back soon!

Problem 92

Think about modern high-speed, large-memory computer systems. What feature of such computers lends itself nicely to the solution of CFD problems using a multiblock grid with approximately equal numbers of cells in each individual block? Discuss.

Check back soon!

Problem 93

What is the difference between multigridding and multiblocking? Discuss how each may be used to speed up a CFD calculation. Can these two be applied together?

Check back soon!

Problem 94

Suppose you have a fairly complex geometry and a CFD code that can handle unstructured grids with triangular cells. Your grid generation code can create an unstructured grid very quickly. Give some reasons why it might be wiser to take the time to create a multiblock structured grid instead. In other words, is it worth the effort? Discuss.

Check back soon!

Problem 95

Generate a computational domain and grid, and calculate flow through the single-stage heat exchanger of Prob. $15-21$, with the heating elements set at a $45^{\circ}$ angle of attack with respect to horizontal. Set the inlet air temperature to $20^{\circ} \mathrm{C}$, and the wall temperature of the heating elements to $120^{\circ} \mathrm{C}$. Calculate the average air temperature at the outlet.

Check back soon!

Problem 96

Repeat the calculations of Prob. 15-95 for several angles of attack of the heating elements, from 0 (horizontal) to $90^{\circ}$ (vertical). Use identical inlet conditions and wall conditions for each case. Which angle of attack provides the most heat transfer to the air? Specifically, which angle of attack yields the highest average outlet temperature?

Check back soon!

Problem 97

Generate a computational domain and grid, and calculate flow through the two-stage heat exchanger of Prob. 15-24, with the heating elements of the first stage set at a $45^{\circ}$ angle of attack with respect to horizontal, and those of the second stage set to an angle of attack of $-45^{\circ}$. Set the inlet air temperature to $20^{\circ} \mathrm{C}$, and the wall temperature of the heating elements to $120^{\circ} \mathrm{C}$. Calculate the average air temperature at the outlet.

Check back soon!

Problem 98

Repeat the calculations of Prob. 15-97 for several angles of attack of the heating elements, from 0 (horizontal)
to $90^{\circ}$ (vertical). Use identical inlet conditions and wall conditions for each case. Note that the second stage of heating elements should always be set to an angle of attack that is the negative of that of the first stage. Which angle of attack provides the most heat transfer to the air? Specifically, which angle of attack yields the highest average outlet temperature? Is this the same angle as calculated for the single-stage heat exchanger of Prob. 15-96? Discuss.

Check back soon!

Problem 99

Generate a computational domain and grid, and calculate stationary turbulent flow over a spinning circular cylinder (Fig. P15-99). In which direction is the side force on the body-up or down? Explain. Plot streamlines in the flow. Where is the upstream stagnation point?

Check back soon!

Problem 100

For the spinning cylinder of Fig. P15-99, generate a dimensionless parameter for rotational speed relative to free-stream speed (combine variables $\omega, D$, and $V$ into a nondimensional Pi group). Repeat the calculations of Prob. $15-99$ for several values of angular velocity $\omega$. Use identical inlet conditions for each case. Plot lift and drag coefficients as functions of your dimensionless parameter. Discuss.

Check back soon!

Problem 101

Consider the flow of air into a two-dimensional slot along the floor of a large room, where the floor is coincident with the $x$-axis (Fig. P15-101). Generate an appropriate computational domain and grid. Using the inviscid flow approximation, calculate vertical velocity component $y$ as a function of distance away from the slot along the $y$-axis. Compare with the potential flow results of Chap. 10 for flow into a line sink. Discuss.

Check back soon!

Problem 102

For the slot flow of Prob. 15-101, change to laminar flow instead of inviscid flow, and recompute the flow field. Compare your results to the inviscid flow case and to the potential flow case of Chap. 10. Plot contours of vorticity. Where is the irrotational flow approximation appropriate? Discuss.

Check back soon!

Problem 103

Generate a computational domain and grid, and calculate the flow of air into a two-dimensional vacuum cleaner inlet (Fig. P15-103), using the inviscid flow approximation. Compare your results with those predicted in Chap. 10 for potential flow. Discuss.

Check back soon!

Problem 104

For the vacuum cleaner of Prob. 15-103, change to laminar flow instead of inviscid flow, and recompute the flow field. Compare your results to the inviscid flow case and to the potential flow case of Chap. 10. Discuss.

Check back soon!