By Dr. A. C. Thompson

Similar languages & tools books

Threaded Interpretive Languages: Their Design and Implementation

Introduces participants possessing microcomputers or minicomputers with minimum peripherals to the layout and implementation of a threaded interpreter as an method of constructing a customary, nonstandard programming language.

Event Management And Best Practices (IBM Redbooks)

Publication through IBM Redbooks

Extra info for Basic Hydrodynamics

Sample text

16 Flow round a cylinder with circulation The force in Hne with the stream will be zero because there are equal p terms at equal y values on the front and rear faces. 18) This remarkably simple result also holds for any other shape in potential flow. 20) V= resultant velocity. The force on a body is then found by integrating the pressures over the surface, so where C = curve of body shape Θ — angle between normal to body surface and x-direction ds = distance along the surface. Where the points on the body are specified, these equations can be expressed as where NP = number of points.

The x, y coordinates of each point on the circle are found in lines 1020 to 1040. These are then transformed to new x, y values in lines 1050 to 1080. At line 1100 the routine returns to line 160. (8) Subroutine 2000 plots the new shape. Scale factor SF, jc-offset XS and y-offset YS are input in lines 2010 to 2030. Lines 2040 to 2050 set up the screen for plotting. In lines 2060 to 2240 the points are plotted. Each x and y value is scaled and offset to fit in the 360 X 160 picture space, in lines 2070 to 2080.

01, reveals the shape of the obstacle more clearly. 0 is added at x= 10, we see that the flow round a finite body is modelled. Increasing the strengths of the sources and bringing them closer together produces a more circular shape. Program SPLOT 100 110 111 112 113 120 130 131 DIM S(21,21),QS(10),XS<10) INPUT "NO. ";J2 INPUT U N0. GRID RQUS? ";U1 INPUT "NO. ";M1 FOR H = 1 TO H1 . INPUT "Q/2PI? ";QS(M) H INPUT "X? 5708 1 295 S(J,K) = S(J,K) + QS(N) * T1 300 NEXT M 310 NEXT K: NEXT J 350 60SÜB 1000 352 INPUT " TO REPLÛT ENTER 1-■ M ;NR 354 IF NR = 1 GOTO 350 500 STOP 1000 REM PLOTS STREAMLINES 1010 HOHE 1020 VTAB 2 4 1030 HGR 1040 INPUT "NO.