fluid_mechanics
708 Appendix A ■ Computational Fluid Dynamics and FlowLab y U a + r θ j Δy (a) + Δ x i x (b) ψ i, j + 1 Δy ψ i – 1, j ψi, j ψ i + 1, j Δ x Δ x Δy ψ i, j – 1 (c) F I G U R E EA.2 storage may be required. Considerable work has gone into forming appropriate grids 1Ref. 62. We consider a grid that is uniformly spaced in the x and y directions, as shown in Fig. EA.2b. As shown in Eq. 6.112, the exact solution to Eq. 1 1in terms of polar coordinates r, u rather than Cartesian coordinates x, y2 is c Ur 11 a 2 r 2 2 sin u. The finite difference solution approximates these stream function values at a discrete 1finite2 number of locations 1the grid points2 as c i, j , where the i and j indices refer to the corresponding x i and y j locations. The derivatives of c can be approximated as follows: and This particular approximation is called a forward-difference approximation. Other approximations are possible. By similar reasoning, it is possible to show that the second derivatives of c can be written as follows: and 0c 0x 1 ¢x 1c i1, j c i, j 2 0c 0y 1 ¢y 1c i, j1 c i, j 2 0 2 c 0x 2 1 1¢x2 2 1c i1, j 2c i, j c i1, j 2 0 2 c 0y 2 1 1¢y2 2 1c i, j1 2c i, j c i, j1 2 (2) (3) Thus, by combining Eqs. 1, 2, and 3 we obtain 0 2 c 0x 2 02 c 0y 2 1 1¢x2 2 1c i1, j c i1, j 2 1 1¢y2 2 1c i, j1 Equation 4 can be solved for the stream function at and to give c i, j c i, j1 2 2 a 1 1¢x2 2 1 1¢y2 2b c i, j 0 1 31¢y2 2 1c 231¢x2 2 1¢y2 2 i1, j c i1, j 2 4 1¢x2 2 1c i, j1 c i, j1 24 c i, j Note that the value of depends on the values of the stream function at neighboring grid points on either side and above and below the point of interest 1see Eq. 5 and Fig. EA. 2c2. To solve the problem 1either exactly or by the finite difference technique2 it is necessary to specify boundary conditions for points located on the boundary of the flow domain 1see Section 6.6.32. For example, we may specify that c 0 on the lower boundary of the domain 1see Fig. EA.2b2 and c C, a constant, on the upper boundary of the domain. Appropriate boundary conditions on the two vertical ends of the flow domain can also be specified. Thus, for points interior to the boundary Eq. 5 is valid; similar equations or specified values of c i, j are valid for boundary points. The result is an equal number of equations and unknowns, c i, j , one for every grid point. For this problem, these equations represent a set of linear algebraic equations for c i, j , the solution x i y j (4) (5)
A.6 Methodology 709 of which provides the finite difference approximation for the stream function at discrete grid points in the flow field. Streamlines 1lines of constant c2 can be obtained by interpolating values of c i, j between the grid points and “connecting the dots” of c constant. The velocity field can be obtained from the derivatives of the stream function according to Eq. 6.74. That is, u 0c 0y 1 ¢y 1c i, j1 c i, j 2 and v 0c 0x 1 ¢x 1c i1, j c i, j 2 Further details of the finite difference technique can be found in standard references on the topic 1Refs. 5, 7, 82. Also, see the completely solved viscous flow CFD problem in Section A6. The preceding two examples are rather simple because the governing equations are not too complex. A finite difference solution of the more complicated, nonlinear Navier–Stokes equation (Eq. 6.127) requires considerably more effort and insight and larger and faster computers. A typical finite difference grid for a more complex flow, the flow past a turbine blade, is shown in Fig. A.5. Note that the mesh is much finer in regions where large gradients are to be expected (i.e., near the leading and trailing edges of the blade) and more coarse away from the blade. F I G U R E A.5 Finite difference grid for flow past a turbine blade. (From Ref. 9, used by permission.) A.6 Methodology In general, most applications of CFD take the same basic approach. Some of the differences include problem complexity, available computer resources, available expertise in CFD, and whether a commercially available CFD package is used, or a problem-specific CFD algorithm is developed. In today’s market, there are many commercial CFD codes available to solve a wide variety of problems. However, if the intent is to conduct a thorough investigation of a specific fluid flow problem such as in a research environment, it is possible that taking the time to develop a problem-specific algorithm may be most efficient in the long run. The features common to most CFD applications can be summarized in the flow chart shown in Fig. A.6. A complete, detailed CFD solution for a viscous flow obtained by using the steps summarized in the flow chart can be accessed from the book’s website at www.wiley.com/college/munson. CFD Methodology Physics Grid Discretize Solve Analyze Problem Geometry Discretization Method Algorithm Development Verification & Validation Governing Equations Structured or Unstructured Accuracy Steady/ Unsteady Postprocess Values Models Special Requirements Implicit or Explicit Run Simulation Visualize Flow Field Assumptions & Simplifications Convergence Interpret Results F I G U R E A.6 Flow chart of general CFD methodology.
- Page 682 and 683: 658 Chapter 12 ■ Turbomachines Th
- Page 684 and 685: 660 Chapter 12 ■ Turbomachines 50
- Page 686 and 687: 662 Chapter 12 ■ Turbomachines Si
- Page 688 and 689: 664 Chapter 12 ■ Turbomachines (2
- Page 690 and 691: 666 Chapter 12 ■ Turbomachines cu
- Page 692 and 693: 668 Chapter 12 ■ Turbomachines SO
- Page 694 and 695: 670 Chapter 12 ■ Turbomachines h
- Page 696 and 697: 672 Chapter 12 ■ Turbomachines He
- Page 698 and 699: 674 Chapter 12 ■ Turbomachines Ro
- Page 700 and 701: V 1 W 1 = W 2 U 676 Chapter 12 ■
- Page 702 and 703: 678 Chapter 12 ■ Turbomachines E
- Page 704 and 705: 680 Chapter 12 ■ Turbomachines U
- Page 706 and 707: 682 Chapter 12 ■ Turbomachines V1
- Page 708 and 709: 684 Chapter 12 ■ Turbomachines Fo
- Page 710 and 711: 686 Chapter 12 ■ Turbomachines Co
- Page 712 and 713: 688 Chapter 12 ■ Turbomachines Co
- Page 714 and 715: 690 Chapter 12 ■ Turbomachines us
- Page 716 and 717: 692 Chapter 12 ■ Turbomachines es
- Page 718 and 719: 694 Chapter 12 ■ Turbomachines 1
- Page 720 and 721: 696 Chapter 12 ■ Turbomachines 55
- Page 722 and 723: 698 Chapter 12 ■ Turbomachines Se
- Page 724 and 725: 700 Chapter 12 ■ Turbomachines Co
- Page 726 and 727: 702 Appendix A ■ Computational Fl
- Page 728 and 729: 704 Appendix A ■ Computational Fl
- Page 730 and 731: 706 Appendix A ■ Computational Fl
- Page 734 and 735: 710 Appendix A ■ Computational Fl
- Page 736 and 737: 712 Appendix A ■ Computational Fl
- Page 738 and 739: A ppendix B Physical Properties of
- Page 740 and 741: 716 Appendix B ■ Physical Propert
- Page 742 and 743: 718 Appendix B ■ Physical Propert
- Page 744 and 745: 720 Appendix C ■ Properties of th
- Page 746 and 747: 722 Appendix D ■ Compressible Flo
- Page 748 and 749: 724 Appendix D ■ Compressible Flo
- Page 751 and 752: Answers to Selected Even-Numbered H
- Page 753 and 754: Answers to Selected Even-Numbered H
- Page 755 and 756: Answers to Selected Even-Numbered H
- Page 757 and 758: I ndex Absolute pressure, 13, 48 Ab
- Page 759 and 760: Index I-3 guidelines for selection,
- Page 761 and 762: Index I-5 hydrostatic: on curved su
- Page 763 and 764: Index I-7 piezometer tube, 50, 105,
- Page 765 and 766: Index I-9 Pressure force, 40 Pressu
- Page 767 and 768: Index I-11 Stress field, 277 Strouh
- Page 769: Index I-13 Weber, M., 29, 350 Weber
- Page 772 and 773: V4.8 Jupiter red spot V4.9 Streamli
- Page 774: V10.3 Water strider V10.13 Triangul
- Page 781 and 782: ■ TABLE 1.4 Conversion Factors fr
A.6 Methodology 709<br />
of which provides the finite difference approximation for the<br />
stream function at discrete grid points in the flow field. Streamlines<br />
1lines of constant c2 can be obtained by interpolating values<br />
of c i, j between the grid points and “connecting the dots” of<br />
c constant. The velocity field can be obtained from the derivatives<br />
of the stream function according to Eq. 6.74. That is,<br />
u 0c<br />
0y 1 ¢y 1c i, j1 c i, j 2<br />
and<br />
v 0c<br />
0x 1 ¢x 1c i1, j c i, j 2<br />
Further details of the finite difference technique can be found in<br />
standard references on the topic 1Refs. 5, 7, 82. Also, see the completely<br />
solved viscous flow CFD problem in Section A6.<br />
The preceding two examples are rather simple because the governing equations are not too<br />
complex. A finite difference solution of the more complicated, nonlinear Navier–Stokes equation<br />
(Eq. 6.127) requires considerably more effort and insight and larger and faster computers. A typical<br />
finite difference grid for a more complex flow, the flow past a turbine blade, is shown in Fig.<br />
A.5. Note that the mesh is much finer in regions where large gradients are to be expected (i.e.,<br />
near the leading and trailing edges of the blade) and more coarse away from the blade.<br />
F I G U R E A.5 Finite difference<br />
grid for flow past a turbine blade. (From Ref. 9,<br />
used by permission.)<br />
A.6 Methodology<br />
In general, most applications of CFD take the same basic approach. Some of the differences<br />
include problem complexity, available computer resources, available expertise in CFD, and<br />
whether a commercially available CFD package is used, or a problem-specific CFD algorithm is<br />
developed. In today’s market, there are many commercial CFD codes available to solve a wide<br />
variety of problems. However, if the intent is to conduct a thorough investigation of a specific<br />
<strong>fluid</strong> flow problem such as in a research environment, it is possible that taking the time to develop<br />
a problem-specific algorithm may be most efficient in the long run. The features common to most<br />
CFD applications can be summarized in the flow chart shown in Fig. A.6. A complete, detailed<br />
CFD solution for a viscous flow obtained by using the steps summarized in the flow chart can<br />
be accessed from the book’s website at www.wiley.com/college/munson.<br />
CFD Methodology<br />
Physics Grid<br />
Discretize Solve Analyze<br />
Problem<br />
Geometry<br />
Discretization<br />
Method<br />
Algorithm<br />
Development<br />
Verification<br />
& Validation<br />
Governing<br />
Equations<br />
Structured or<br />
Unstructured<br />
Accuracy<br />
Steady/<br />
Unsteady<br />
Postprocess<br />
Values<br />
Models<br />
Special<br />
Requirements<br />
Implicit or<br />
Explicit<br />
Run Simulation<br />
Visualize<br />
Flow Field<br />
Assumptions &<br />
Simplifications<br />
Convergence<br />
Interpret<br />
Results<br />
F I G U R E A.6<br />
Flow chart of general CFD methodology.