|
[Sponsors] |
Issues with DES simulation of centrifugal compressor |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 18, 2016, 09:23 |
Issues with DES simulation of centrifugal compressor
|
#1 |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Hi all,
I am modelling the near surge point of a centrifugal compressor with DES in CFX 15. I've attached the images of blending function and strain rate invariant. I am seeking some direction on couple of points mentioned below: 1. As the CFX modelling guide says "In regions where the function is zero, the LES model is used, and the region where its value is one, the RANS model is activated". I can't see any 0 in my model, the lowest is of the order of 10^-2 but not 0. Can anyone share their knowledge if 0 means absolute 0 in this case or a small value? Also, what about the grey region between 0 & 1. Is there any transition function involved or the grey region is to be avoided using grid / timestep refinement? 2. In the case of strain rate invariant, I can see only one scale. The theory guide suggests using reduced timestep, refine grid, use Fsst = F1 or Fsst = 0. I did try reduced TS and grid size but didn't see any more scales in invariant isosurface and don't really understand the Fsst = F1 / 0. Can someone please share some pointers to what am I doing wrong and how can I navigate through this problem. 3. Acoustic courant number is always 999.99, any ideas why? Please correct me if I am wrong, CFX isn't highly dependent on courant number to accurately resolve unsteadiness and I can get away without maintaining the courant number as unity?? In this case, the RMS courant number is 19 and max are 952. Any inputs on this would be great. 4. I am currently using the 4-degree rotation of wheel per timestep as physical TS value and using 8 coefficients loops. I've seen the various post by Glenn and Ansys resources emphasising 3-5 loops to achieve the convergence in each TS by using necessary TS. I did use the lowest TS I could use to (computation power constraint) but still couldn't achieve the convergence in fewer loops. Any direction in here on what should be done would be incredibly useful. I have attached the plots of residuals too. 5. The objective is to model the in-duct pressure fluctuations leading to acoustic spectrum. Any do's and don't from this perspective are highly welcomed. Thank you for contributing. |
|
July 18, 2016, 09:46 |
|
#2 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
I will answer your questions which I think I might be able to help with:
3) Then run a smaller time step and see if it makes a difference. It is best to establish these things yourself using your simulation. 4) 4 degrees per time step sounds very large, especially for a DES/LES model. Are you using double precision numerics? Also you will need to check your mesh is OK. Have you done a mesh sensitivity check? |
|
July 18, 2016, 10:35 |
|
#3 |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Thank you for the inputs and reply, Glenn.
3. I am planning to do a time step independence study but the computational constraints are in the process of sorting out. 4. 4 deg is large but there is published literature with similar timestep size. I am running one model with 2 deg TS too but results aren't substantially different; RMS courant number drops to 12 but acoustic courant number is still 999.99 - Yes, I am using double numerics. - I've done the mesh sensitivity/independence study but it's based on RANS rather than DES (i've attached the mesh study screen shot with the highlighted mesh as the one used in this model). I am not sure if it's the best way to go about it but computational constraints are forcing my hand this way. Although the cfx modelling guide does say that "While the grid resolution requirements are not significantly higher than RANS for simulations, the time resolution imposes high CPU demands." Based on what I've read from your post, you prefer to use adaptative TS based on courant number target, you think I should try that approach instead even CFX isn't highly courant number dependent? Thank you |
|
July 18, 2016, 21:38 |
|
#4 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
Mesh and time step sensitivity of DES/LES will be very different to RANS. You will have to redo it for a DES/LES simulation.
No, I recommend adaptive time steps homing in on 3-5 coeff loops per iteration. |
|
July 19, 2016, 08:32 |
|
#5 | |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Quote:
I've set-up simulation with adaptive timestep and target of 3-5 coefficient loops instead of courant number and also dropped the convergence target for residual RMS from 10^-6 to 10^-4. I don't think it'll possible for me with current resources to do a grid independence study for DES. It would be great if I can do a timestep independence study alone. I'll update you how it goes and thank you again. I appreciate your time. |
||
July 19, 2016, 08:53 |
|
#6 | |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
Quote:
|
||
July 19, 2016, 09:02 |
|
#7 | |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Quote:
I'll measure the acoustic PSD and correlate it with numerical results. I'll only work in the zone of similar trends between two and reason the computational constraints for not working on other regions. Moreover, if i can see positive design effects using the current set-up, I can test to see the actual improvements and that should be enough for thesis ( I hope). |
||
July 19, 2016, 20:47 |
|
#8 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
You don't have to do the grid independence study on the full thing. Either cut out a section of the geometry or use a simpler analogy.
If it takes months to run then it is even more important to do all the checks in advance or you will be wasting months of time! |
|
July 21, 2016, 10:18 |
|
#9 | |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Quote:
1. I ran the adaptive timesteps with target min and max coefficient loops as 3 & 5 respectively. I specified the timestep range from 5 degrees /TS to 0.5 deg/TS. Cfx ran till the least timestep value i.e. 0.5 deg/TS and still wasn't able to achieve the target residuals of 10^-4 and hence, was running with 8 coefficient loops (max specified in solver control). I've attached the relevant screenshots. I don't think it'll be wise to reduce timestep any further, the computational requirements are already huge. I was wondering in your experience what would be the consequences of not following 10^-4 residual levels to the dot? This being said I've come across one paper which specified the residual targets as 10^-3 for near surge points. It's peculiar that not many research groups talk about the convergence criterion used by them in their papers and validates directly with test results. 2. Would cut out section is credible for grid independence study? As in, if I only model impeller section and leave out rest, although I do think immediate inlet and diffuser + volute sees a lot of (dynamic) instability (from steady state results and published work). You think if this will be a valid approach? and also, if you could elucidate on analogy, I am not sure what you meant by using analogy for grid independence. Thank you. |
||
July 21, 2016, 19:37 |
|
#10 | ||
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
Quote:
Quote:
|
|||
July 25, 2016, 06:37 |
|
#11 | |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Quote:
One more thing, I was wondering if you've any directions for implementing anechoic termination/non-reflecting bc in cfx in order to curb the spurious reflections. Thank you. |
||
July 25, 2016, 06:38 |
|
#12 |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
It was in quotes in last reply, my bad.
I extended the adaptive time step study and relaxed the minimum time step limit. To get the initial results, I ran the simulation with frozen (MRF) rotor interface rather than sliding mesh (transient rotor-stator). The residual target to 10^-4 rms is achieved with 0.095 degrees/timestep giving rms courant number of 0.56. I'll run with this timestep and compare its acoustic psd with that of bigger time steps. This, in turn, should reinforce time step independence too. |
|
July 25, 2016, 07:46 |
|
#13 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
Non-reflecting BC: A simple way of doing this which does not always work is to deliberately grossly coarsen your mesh immediately at the boundary. The wave hits the coarse mesh and gets dissipated by the coarse mesh. Use a GGI interface to connect the very coarse mesh to the finer main mesh. This method is crude but sometimes can work.
There is also a beta feature for non-reflective boundary conditions. Enable beta features in CFX-Pre and give it a try as well. |
|
July 28, 2016, 07:40 |
|
#14 |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Hi Glenn, thank you for the information and tips to employ nrbc.
I explored the beta options in cfx but as far as I gather, nrbc can be only employed at the inlet and not on the outlet. Is that so and if yes, what you think is the intent of it? I've come across another method using "beamforming" (published by researchers in CMT Valencia) to resolve the reflection issue. I am reading and try to employ beamforming but haven't had any luck at this point in time. Would you have any idea about beamforming? Thank you |
|
July 28, 2016, 08:01 |
|
#15 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
I have never used the non-reflecting BC in CFX. It is a beta feature anyway so it is not fully supported or documented. So I can't help you there.
I can't help you with beamforming either. If you want to go that way you will have to work it out for yourself. Don't forget my suggestion of the grossly coarse mesh at the boundary - it is simple to do and has given quite good results for several people. It is so simple that surely it is worth a try. |
|
July 28, 2016, 10:20 |
|
#16 |
Member
Sidharath
Join Date: May 2015
Location: UK
Posts: 65
Rep Power: 11 |
Thanks a lot for all your replies Glenn, I really appreciate your time. I should put you as my co-supervisor
I'll surely try the coarse mesh at boundaries and keep you posted about the success of methods and try to validate it with beamforming and experimental results as and when they work out. Thank you once again. |
|
March 21, 2018, 06:30 |
|
#17 |
Member
Tingyun YIN
Join Date: Apr 2017
Posts: 31
Rep Power: 9 |
Hi, Sidharath and Glenn Horrocks,
Currently, I am working on the DES simulation. But I have some problems about the numerical scheme and want to ask for help. When we switch on DES or SAS model in CFX, an additional option occur in the solver control plane. It's called 'CDS blending'. As is well know, different numerical schemes could be applied to different regions when using DES or SAS model: CDS for LES region and Second order for RANS region. So, here are my questions: 1. When we switch on the default DES model in CFX, the default numerical scheme is hybrid scheme? The advection scheme option is available only for RANS region because CDS is used for LES region as default. So, when 'CDS blending' is activated, there are two functions: The first one is to prevent the grid-induced separation (as default, minimum RANS blend= BF1=0, maximum RANS blend= BF2=1) and avoid oscillations due to the central difference scheme (as default, Maximum Courant Number=5, Limiter Exponent=1). 2. If the description of question 1 is correct, so, when we switch on the default SAS model and activate 'CDS blending' function, why minimum RANS blend and maximum RANS blend still occur? As we all know, no grid-induced separation occurs when using SAS model. So, what should I do? Don't activate them? To avoid oscillations of CDS, Maximum Courant Number=5, Limiter Exponent=1 can be activated. 3. If the description of question 1 is wrong, so, when we switch on the default DES model, we have to activate the 'CDS blending' to activate the hybrid scheme? Any comments are appreciated. Thanks in advance. |
|
March 21, 2018, 07:30 |
|
#18 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,867
Rep Power: 144 |
I do not know the answers to your questions so cannot help. All I can suggest is to read the documentation carefully and if you can't work it try ANSYS training and support.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
March 21, 2018, 07:35 |
|
#19 | |
Member
Tingyun YIN
Join Date: Apr 2017
Posts: 31
Rep Power: 9 |
Quote:
Fine. After checking Help Document, I still do not get the answers. That's the reason why I post questions here. Anyway, I will try ANSYS training and support later. Thank you very much! |
||
Tags |
cfx 15, compressor, des turbulence model |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Centrifugal Compressor Efficiency | albertotato | CFX | 3 | May 2, 2020 03:54 |
Steady State Centrifugal Compressor Segment | vigges | OpenFOAM Running, Solving & CFD | 0 | June 16, 2014 05:46 |
Centrifugal compressor modeling by real gas | omidiut | CFX | 3 | May 8, 2013 04:22 |
simulation the centrifugal compressor | layth | STAR-CCM+ | 1 | November 30, 2011 04:17 |
centrifugal compressor | siva appanna | Main CFD Forum | 5 | February 13, 2006 22:07 |