# Unexpected results for 2D cylinder flow

 Register Blogs Members List Search Today's Posts Mark Forums Read

March 6, 2019, 17:32
Unexpected results for 2D cylinder flow
#1
New Member

D. de Vries
Join Date: Feb 2019
Posts: 3
Rep Power: 4
I am trying to learn how to use SU2. As a test, I am trying to solve the 2D, incompressible, inviscid flow around a cylinder using the Euler equations. I have created a simple triangular mesh using Gmsh (see attachment.), and have changed the CFG file from the Inviscid Hydrofoil example on the SU2 website to work with my mesh instead (see attachment). I am then running the case in parallel using the command
Code:
`parallel_computation.py -f cylinder.cfg -n 16`
. The simulation exits after the convergence criteria are met. The inlet flow conditions are 10 m/s in the positive x-direction, so I expect a symetric flow field and no residual force in the y-direction. However, the flow field is not symmetric and there is a residual force in the negative y-direction. Note also that the velocity of the free-stream is not 10 m/s. I have attached two images depicting a set of streamlines and the pressure coefficient and velocity magnitude respectively. Can anyone tell me what I'm doing wrong here?
Attached Images
 flowFieldAndPressureCoefficient.jpg (54.9 KB, 18 views) flowFieldAndVelocityMagnitude.jpg (72.0 KB, 14 views)
Attached Files
 cylinder.zip (3.1 KB, 6 views)

 March 10, 2019, 21:12 #2 New Member   Alex Join Date: Jul 2018 Posts: 8 Rep Power: 4 Glancing over the .cfg file you posted, here are a few things that I would change/adjust: -It looks like you are specifying a CFL of 1e6 on the first iteration, then switching to an adaptive one in the interval [1, 100] for all other iterations (you can monitor this value in the history file). -You are specifying a convergence criterion of 10 orders of magnitude reduction, but you also have the minimum value of the residual set to -6 (log10). Perhaps your simulation is exiting prematurely; you can try decreasing RESIDUAL_MINVAL to see if the same setup will still converge. -I would also consider turning off multi-grid until you are achieving appropriate results. I’m not sure these modifications will solve your problem, but it may help point out other issues. As for the velocity magnitude not equaling 10 m/s: The values are non-dimensionalized with respect to the initial values by default. You can alter this by using the “INC_NONDIM= ” option. Reference values will be output at the beginning of the simulation, or you can view them in the forces_breakdown.dat file. Hope this helps, Alex

 Tags cylinder, inviscid flow, test case

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post aakie OpenFOAM Running, Solving & CFD 3 November 13, 2018 04:39 Scabbard OpenFOAM Running, Solving & CFD 37 March 21, 2016 16:16 Wyrold Main CFD Forum 0 October 22, 2015 08:48 raf1111 OpenFOAM 1 December 16, 2013 09:45 ozzythewise Main CFD Forum 8 June 13, 2012 06:24

All times are GMT -4. The time now is 13:57.