Reference no: EM13708859
Question: You need to get flow around cylinders in a channel by creating a grid and specifying the intial flow velocity.
1. Consider the 4 × 2m channel to be covered by an x-y grid of finite spacing ? = ?x = ?y. The aim is to determine the stream function at each grid point on the grid. Select a grid spacing such that the height of the channel is eg. 50? and the width 75?.
2. Use meshgrid to define arrays x arr and y arr that contain the x and y co-ordinates of all points in the grid. Define an array Psi in which to store the estimated stream function values at each grid point.
3. Put Psi = Uo * y arr this will generate a uniform stream function at each grid point.
4. Now accumulate the contributions from each cylinder, eg. for i = 1:no_of_cyls x_norm = x_arr - x_c(i); y_norm = y_arr - y_c(i); Psi = Psi - Uo .* radii(i)^2 .* y_norm ./(x_norm.^2 + y_norm.^2); end
5. Specify a starting guess for the stream function by evaluating the analytical estimate at each point on the grid.
6. Calculate the corrected stream function Ψc = Ψ1+Ψ2+Ψ3+Ψ4 4 for each grid point (xi , yj ) inside the domain.
7. Set the stream function of all points coinciding with the boundary streamline of each cylinder as the average of all such boundary values.
8. The top and bottom edges of the domain are streamlines, so the stream function must be constant on these edges.
9. At inlet and outlet the stream function of the boundary points should be reset to their original free-stream values with no interior obstacles (uniform flow U0y) 21
10. Every 50 iterations visualise the stream function using the Matlab functions contour(Psi) and mesh(Psi). Velocity vectors can be plotted using the quiver command, you will need to differentiate the stream function to obtain them.
11. Iterate: repeat from step 6 until changes at each iteration are below some small level (eg 0.01%)
As its a 2d grid the mass is accumulated as the y value increases.
An inital stream function is calculated by adding up 4 points and taking the mean.
the stream function is then corrected until the change each time is smaller than a given value
The problem i have is that not all the 6 cylinders are showing up on my final contour plot
|
Atmosphere is connected to a major artery
: A vertical tube open to the atmosphere is connected to a major artery in the arm of a person. To what maximum height does the blood rise in the tube? Take SG blood = 1.04, and the systolic blood pressure of the individual to be 120 mm-Hg.
|
|
Find the net upward force on the inflated balloon
: A rubber balloon has a mass of 1500 g when it is deflated, i.e., empty. When the balloon is filled up with Helium gas (He = 0.18 kg/m3 ), it takes the form of a sphere with a 4 m diameter. (V-sphere= (4/3)(pi)(r)^3) Find the net upward force on the i..
|
|
Change in internal energy and enthalpy of oxygen
: Oxygen at 450K is heated to 850K. Assume oxygen as an ideal gas. Determine the change in internal energy and enthalpy of oxygen per unit mass, using: Data for line over u and line over h from the oxygen table
|
|
What is the work done on or by the steam
: A piston-cylinder contains 0.1 kg of superheated steam at 180 ¯ C and 5.0 bars. The piston is frictionless and freely floating. Heat is added until the final volume is twice the initial volume.
|
|
Creating a grid and specifying the intial flow velocity
: Need to get flow around cylinders in a channel by creating a grid and specifying the intial flow velocity.
|
|
Radius of moons orbit relative to the earth
: What is the total change in velocity required to go from a 200 km altitude circular, Earth orbit to circular orbit with a radius of 384,000 km (the mean radius of moons orbit relative to the earth, but ignore the moons presence}?
|
|
What are the fin efficiency and effectiveness
: Aluminium fins of triangular profile are attached to a plane wall whose surface temperature is 250 oC. The fin base thickness is 2 mm, and its length is 6 mm. The system is in ambient air at a temperature of 20 oC, and the surface convection coeffici..
|
|
Expression for the rate of heat loss to surrounding fluid
: A spine attached to a wall maintained at a uniform temperature Tb has the shape of a circular cone with base diameter b, height L and a concave parabolic profile, and it is exposed to a fluid at a uniform temperature T?.
|
|
Quality of the refrigerant and the pressure at final state
: A well-insulated piston-cylinder device initially contains 5kg of saturated mixture of refrigerant-134a at -12°C, in which 2.55kg of the refrigerant is in the liquid phase. An electric switch is turned on, and a 110-V source supplies current to a res..
|