Develop a MATLAB (or programming language of your choice) codeto show the flow field for lifting cylinder flow using thefollowing specifications.
Set your free-stream velocity, cylinder radius, and vortexstrength in the code. There are a total of four cases for which youwill run your code (or have all four in one script). You will onlybe changing the vortex strength for the different cases. Therefore,set ?∞ = 1 and ? = 1 for all cases. Set your vortex strength to thefollowing values for each case, Case 1: ? = 0, Case 2: ? = 8, Case3: ? = 16, Case 4: Find the value for vortex strength that producesone stagnation point on the surface of the cylinder.
Set V = R = 1; Set Γ = value desired for specific case; Use meshgrid to create my domain, but do this in polar coordinates;Calculate x and y coordinates based on mesh grid; Calculate streamfunction as a function of r and theta; Calculate Vr and Vtheta as afunction of r and theta; Convert Vr and Vtheta to u and v; Plotcontours of stream function; Plot velocity field using quiverfunction;