Non convergence for a cyclic Geometry/values for Spalart-Allmaras model
I am trying to simulate 3D air flow inside the vanes of a closed ventilated rotor.
I am using SRFSimpleFoam/Mixer as my sample case, and was able to run a successful initial simulation with my cyclic geometry, with slight modifications in the files from the tutorial case. However,when I tried to check the convergence for my initial variables, p,Urel,epsilon, and omega, I hardly see any convergence in the results. Initially I checked the mesh, and mesh seems fine. Skewness:(0.79) and non-orthogonality Max:56.2308, avergae:16.9368 Now, I as am trying different settings in the system/fvSchemesand system/fvSchems, I plan to take it to to next step by moving to Spalart-Allmaras model(Using as template:http://rapidof.com/centrifugal-pump-simulation/). After looking through a lot of forums, I was able to find the variables like nu,nut and nuTilda, however, I still cant find the values for DnuTildaEff,DepsilonEff,DkEff and nuEff. ... I have looked through a lot of tutorials, to find any files I could use as a template(I plan to not create the input file from scratch , since i am new, and I done know what class types should I use for my boundary layers)..... Does anybody by any chance know, where can I find a sample tutorial files with those values, or if you could forsee any other issue that could cause divergence? -------------------------------------------------------------------------------- I am bit new to OpenFoam, so I might have missed some documentation |
I found the tutorials for turbulence models in the $FOAM_SRC folder.
However, they only contain the .H and .C files. Is there any way i could figure out what is the class types of my input variables. Is there anywhere, where I can get a simple turbulent 1D or 2 equation tutorial models? And the second question about improving the convergence of my model still remains. |
Solved the first issue
Only the second issue about imporving the covergence remains I am playing around with fvSchemes currently, is there any other thing that could be wrong |
Arpit, I don't know if that is suitable for your problem, but I had a problem with cyclic as well. I was trying to get a simple periodic pipe or channel flow running and could not get any convergence, no matter what I set in fvSchemes and fvSolution. This was driving me insane. I figured out that I can do exactly what I want with "mappedPatch" boundaries, namely directing the outflow of one boundary to the inflow of a second boundary. For the case "piece of pipe" with periodic boundaries, I modified the "boundary" file in the "constant\polyMesh" directory:
Code:
/*--------------------------------*- C++ -*----------------------------------*\ In your \0 boundary files you must set mapped as type: Code:
/*--------------------------------*- C++ -*----------------------------------*\ Cheers, Philipp |
Hi Philipp
Thank you, taking your time out to review my posts. I am currently facing a bigger problem than convergence. Although , my residual is 0.002 > (desired based on forum recommendations= e-4),I am not too worried about it. However, my results from the openfoam are completely off, and I feel I have inputted wrong boundary condition somewhere . If you could please take some more time, and review my case, where I went, I would really appreciate it. I tried your settings with the type mappedFace(Settings 2), but nothing has changed. Please follow the below link, which describes my case files and results https://www.dropbox.com/sh/js5zmei0e...EhRP0OwDa?dl=0 MY CASE I am trying to do a CFD on a car rotor, especially the air flow inside the vanes of a ventilated rotor. I am rotating the rotor anticlockwise, such that it acts as a pump, and pushes the air out. I am trying to model it as a stationary rotor, and moving fluid around it. There is no forced flow, everything is induced by the angular velocity of the rotor. Thanks a lot for the help .... , unfortunately , I am still stuck |
1 Attachment(s)
Attached is an excel sheetAttachment 42747 detailing my boundary conditions.
I have tried changing the direction of the rotation(angular velocity) in SRFProperties, but my flow is still flowing inwards as compared to outwards(out of the vane like an impeller). |
Arpit when I look into the "settings_1" fvSchemes I see very bad schemes for the first run.
Try this: Code:
divSchemes |
Hi Philipp
I changed my fvSchemes settings based on what you recommended and also increased my number of iterations until my residuals were mostly steady, however my results are still not correct. I am modelling this blade as an impeller blade, and no matter what what rotation I specify my velocity vectors are pointing inwards(towards the center) as compared to outwards(as it should). You can see a few pictures of my results in "Settings1_Results " at the same link I think there is something else wrong, which I cant figure out. I feel my setup is somewhere wrong, and does not capture the necessary physics. I dont how,or what is wrong |
Hey guys
I am really stuck on simulation and would appreciate any help or suggestions you could provide. I am trying a CFD simulation of a ventilated car rotor. I am trying to run it as am impeller , such that the flow of the air is outwards(away from the center). However , my simulation results indicates the direction of velocity vectors to be inwards(I tried reversing the direction of my rotation, but the simulated flow is still inwards). I am trying to use the frozen rotor approach where,my cyclic geometry is stationary and there is air all around it(which is moving based on a rotating frame of reference) , which .Thus, I believe there is something wrong with setting up my problem.My original post can be found here http://www.cfd-online.com/Forums/ope...tml#post569114 Link to my result and setup folder is-https://www.dropbox.com/sh/js5zmei0e703ptm/AAAwMqsDmycfVFMwEhRP0OwDa?dl=0 Thanks for all the help you could provide. Arpit |
Arpit, your numerical settings are not good.
You stop at t=0.3 and don't see what happens after that. You should (at the beginning) always set an incredibly high endTime, such as 1e9 and let the simulation run until all residuals settled down. Then you will see, that your fvSchemes are not able to obtain a convergent solution. Try these settings for a safe run: Code:
divSchemes |
1 Attachment(s)
Thanks Phillip
I changed my boundary conditions and increased my number of iterations and solution worked. At least the cyclic geometry is working as an impeller(throwing air away from the center) as it is supposed to. Thanks a lot. However, if you happen to know how to create streamlines from inlet(a boundary patch) as compared to a point, let me know. I followed up this link - http://www.cfd-online.com/Forums/ope...seed-type.html , but am having trouble doing the step number one. Since my inlet is a curved plane, it is hard for me to specify an accurately enough. Currently i am using a point in middle of my inlet as a point source seed and playing around with seed radius to cover my geometry. My current streamline file is as shown(i wished I could specify the direction of the streamline flow and change it from horizontal to along the curved geometry) Attachment 42951 Thanks a lot, btw for your responses |
I never did this, sorry.
|
Cool.
Hopefully, somebody else might be might be able to help me out. Thanks a lot, though |
All times are GMT -4. The time now is 20:20. |