
[Sponsors] 
February 6, 2009, 12:40 
Hello. I'm pretty new to OF bu

#1 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
Hello. I'm pretty new to OF but I've been browsing the forums for some time now. I am doing some research on an airship and I want to get lift, drag, moment coefficients for it in various configurations.
I've borrowed Mr. Barnaby's vwt setup and I'm trying to apply it to my own models. I have been trying to run simpleFoam on the model "solarship.stl" in the directory, but was consistently getting really unreasonable results, with drag orders of magnitude larger than lift. I then made a model of a 3d airfoil (e197) section with high aspect ratio and the results again are consistently wrong compared to XFOIL (even with 2D>3D errors expected) OF ended up stabilizing with a Cd = 0.2908 and CL = 0.0754 but XFOIL gives roughly Re = 150000 Cl = 0.2741 Cd = 0.01423 So I'm obviously doing something wrong. I was wondering if someone would have a look at my setup and let me know if there's something obvious and stupid that I have wrong? Thanks for the help vwtairship.tar.gz 

February 6, 2009, 12:44 
PS. Sorry, the stl of the airs

#2 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
PS. Sorry, the stl of the airship was too large to post, but the airfoil is present.
The problem was present with both models so perhaps it's an error in my setup. 

February 9, 2009, 01:51 
Hi Kyle
Well just a suggesti

#3 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 17 
Hi Kyle
Well just a suggestion... reduce grid size...check might need to run for longer no. iterations... Just suggestion might not be correct... Regards sachin 

February 9, 2009, 02:17 
Hello Kyle,
I see that you us

#4 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 20 
Hello Kyle,
I see that you used a low Re ke turbulence model. The first paper I found on the internet regarding turbulence models comparison says that a high Re turbulence model with one equation gives a better prediction of the transition point than a high Re ke or a low Re ke (like Launder and Sharma). As far as I know, OpenFOAM has such a model implemented (Spalart Allmaras), why don't you give it a try? I hope this is helpful, Dragos 

February 9, 2009, 02:41 
Hello
Dragos thanks for that

#5 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 17 
Hello
Dragos thanks for that paper but to what i read i didnt find it opposed high keps model but combines high keps model with one equation model Correct me if i am wrong... does OpenFOAM have Spalart Allmaras combined with high keps model....if so Please can you give me some details about it? regards, Sachin 

February 9, 2009, 02:56 
No, you're right. They are usi

#6 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 20 
No, you're right. They are using a high Re ke in the bulk flow and a one equation model near wall.
I misread the paper, and I don't think there is a similar implementation in the current OpenFOAM release. However, my point was that one equation turbulence models (like Spalart Allmaras) usually gives better results for airfoils like bodies, than two equation models (or at least fluent manual claims that, if I remember correctly). Dragos 

February 9, 2009, 03:15 
thanks dragos...
was about it

#7 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 17 
thanks dragos...
was about it implement in my case Sachin 

February 10, 2009, 23:51 
Thanks very much for the input

#8 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
Thanks very much for the input, everyone. I did realize that the problem was in the turbulence because the viscous forces were much higher than they should've been. Because my test case should not need turbulence I ended up disabling it.
Afterwards, using simpleFoam without turbulence I was able to get results very close to the predicted values, but I am having trouble achieving convergence. I also cannot get convergence on my airship model. Any suggestions for improving the convergence speed in simpleFoam? The results stabilize around 40 iterations and the residuals decrease for the next 100 or 200 and then start to drift and eventually the solution explodes. Thanks for the help! 

February 13, 2009, 02:55 
So I am still unable to get an

#9 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
So I am still unable to get any decent convergence. The residuals get down on the order of 10^3 and will hover there until the end of my 1000 iterations.
I'll attach my fvSolutions and fvSchemes and a log of the most recent run. I used blockmesh at a resolution of (25,25,25) for a test section of (+25,28,20) and then used snappyHexMesh to create a mesh around the airship. Here is the checkMesh output: /**\  =========    \ / F ield  OpenFOAM: The Open Source CFD Toolbox   \ / O peration  Version: 1.5   \ / A nd  Web: http://www.OpenFOAM.org   \/ M anipulation   \**/ Exec : checkMesh Date : Feb 13 2009 Time : 01:52:56 Host : Orion PID : 12753 Case : /home/kyle/OpenFOAM/kyle1.5/run/solarShip nProcs : 1 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = constant Time = constant Mesh stats points: 266591 faces: 596209 internal faces: 502073 cells: 165657 boundary patches: 4 point zones: 0 face zones: 0 cell zones: 1 Number of cells of each type: hexahedra: 132556 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 33101 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Topological cell zipup check OK. Face vertices OK. Faceface connectivity OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface inlet 625 676 ok (not multiply connected) outlet 625 676 ok (not multiply connected) channelWalls 2500 2600 ok (not multiply connected) car_Mesh 90386 92000 ok (closed singly connected surface) Checking geometry... Domain bounding box: (25 28 20) (25 28 20) Boundary openness (1.8302e17 7.9149e18 1.07047e16) OK. Max cell openness = 1.85863e16 OK. Max aspect ratio = 1.4 OK. Minumum face area = 0.003125. Maximum face area = 4.48. Face area magnitudes OK. Min volume = 0.00021875. Max volume = 7.168. Total volume = 111076. Cell volumes OK. Mesh nonorthogonality Max: 32.0305 average: 14.0889 Nonorthogonality check OK. Face pyramids OK. Max skewness = 1 OK. All angles in faces OK. Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1 All face flatness OK. Mesh OK. And the files: fvSolution fvSchemes logtruncatedrelaxed.txt Any help getting this to converge is most appreciated. Thanks! 

February 16, 2009, 15:21 
Sorry to keep posting, but I a

#10 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
Sorry to keep posting, but I am still completely unable to get any convergence on these test cases. I've tried increasing the mesh resolution significantly (up to ~1.5m cells)
I've tried initializing the flow field with potentialFoam. This seems to bring the residuals down to 10^3 in fewer iterations but it still will not go to convergence. I also tried playing around with the relaxation factors and orthogonality correctors, which did not seem to help much. I just want to get the forces for an airship in incompressible, viscid, steady, nonturbulent flow. Any suggestions are most appreciated! 

February 16, 2009, 15:31 
Also, I just want to verify th

#11 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
Also, I just want to verify that these boundary conditions make sense for a "virtual wind tunnel"
For 0/U: FoamFile { version 2.0; format ascii; root ""; case "3d"; instance "0"; local ""; class volVectorField; object U; } dimensions [0 1 1 0 0 0 0]; internalField uniform (12.803 0 2.257); boundaryField { inlet { type fixedValue; value uniform (12.803 0 2.257); } outlet { type zeroGradient; } channelWalls { type fixedValue; value uniform (0 0 0); } car_Mesh { type fixedValue; value uniform (0 0 0); } } and for 0/p: FoamFile { version 2.0; format ascii; root ""; case "3d"; instance "0"; local ""; class volScalarField; object p; } dimensions [0 2 2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } channelWalls { type zeroGradient; } car_Mesh { type zeroGradient; } } As you can see, I wanted the ship to be flying at an angle of attack of 10deg at V(magnitude)=13m/s. Is this the proper way to set up an angle of attack flight? Because when I look at the flow in paraView after some iterations it seems like the flow levels off and decreases the angle of attack. Also, when I initialize the flow with potentialFoam it "levels off" the flow. Is this a boundary condition issue? Thanks! 

February 16, 2009, 18:57 
Hello again. I think there is

#12 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
Hello again. I think there is something wrong with my boundary conditions. After running another simulation I am looking at how the flow develops with the iterations, and it strangely produces a region of flow reversal in the tunnel. Check out the images:
Now it's pretty obvious that this is not going to reach a steady state solution, or at least it will not be what I want. Any ideas how I can fix this so that it actually represents my model flying at an angle of attack? Thanks! 

February 16, 2009, 19:01 
Sorry. Images didn't get uploa

#13 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
Sorry. Images didn't get uploaded.


February 17, 2009, 00:47 
Hello
I am not very sure...j

#14 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 17 
Hello
I am not very sure...just a thought if channel wall is the outer wall of wind tunnel...it should be considered that it plays no role in your flow field...maybe 0 0 0 for channel wall is not correct....try with zeroGradient...which just puts the value whatever is in the cell next to it Hope this helps... Sachin 

February 17, 2009, 04:21 
Hello Hyle,
I had a problem s

#15 
Senior Member
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 23 
Hello Hyle,
I had a problem similar to yours a couple of days ago... The solution was to impose in the upper, lower, front and rear surface the same bc you have at the inlet. This will simulate better a freestream condition. For some more information, have a look here! Good luck! Maddalena 

February 17, 2009, 06:16 
Thanks very much for the input

#16 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
Thanks very much for the input. It did end up being a problem of the boundary conditions.
I ended up setting the front, top, back and bottom all to the same constant velocity, and had zero pressure gradient across all these faces. This has given me a flow like this: Now I think I will be able to get most of the data I need. I have two questions though 1) Is it possible to record multiple moments on a patch? I want to get Mx, My, Mz all referenced at the same point, and obviously don't want to run the simulation three times to do so. 2) I want to find the dynamic stability derivatives and to do this I will need to accelerate the flow field. Is there a way to do this? Any suggestions on setting up a time varying flow field? Thanks so much for your help! 

February 18, 2009, 06:19 
hello,
1. not much idea...but

#17 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 17 
hello,
1. not much idea...but guess should be possible 2. time varying flow field ...change the solver to transient solver sachin 

February 19, 2009, 03:02 
So far I am just going to reco

#18 
New Member
Kyle Dyroff
Join Date: Mar 2009
Posts: 10
Rep Power: 17 
So far I am just going to record the force data and go back to calculate the other moments later. Not too worried about that.
I am however, worried about setting up a time dependent flow field. I need to find the stability derivatives of the aircraft and so, for example, Cxq will be how the x force on the aircraft varies with respect to rolling rate. Or Cxu(dot) is how the xforce varies with respect to acceleration of the flow. If anyone can point me in the direction of how to do this I'd be much appreciative. I feel like I need to set up some flow source that has a variation in time, but have no idea how to do this. Thanks, Kyle 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
reasonable mixing rate for C2H6 Oxygen mixture  N. Schiepel  CFX  0  September 12, 2008 05:57 
Unable to Initialising to M=1.71  ravi  FLUENT  2  September 14, 2007 03:37 
reasonable result for air flow in city buildings?  George  FLUENT  0  August 21, 2006 21:36 
Unable to block  Tang  CFX  2  October 7, 2005 10:12 
This result reasonable help me  zou_mo  OpenFOAM Running, Solving & CFD  4  August 16, 2005 05:49 