
[Sponsors] 
Residual Convergence Issues with Overset Meshing? 

LinkBack  Thread Tools  Search this Thread  Display Modes 
May 1, 2017, 19:00 
Residual Convergence Issues with Overset Meshing?

#1 
New Member
Dan
Join Date: Aug 2016
Location: Montreal, Canada
Posts: 7
Rep Power: 5 
Hi all,
I am working on a simulation of a vertical axis wind turbine operating at low TSR and thus in the dynamic stall regime. I am using an overset region around each blade and polygonal elements everywhere. I am having issues with getting the residuals of momentum in particular to drop in an expected manner  i.e. in the range of a 3 order of magnitude drop in a reasonable number of inner iterations. So far I've been primarily using the segregated flow solver, but the residual trend for both momentum quantities is often very irregular and oscillatory (continuity is fine and behaves as expected). To try and solve the problem I've been doing some testing whereby I use the same model setup but look at a static airfoil case with an angle of attack of 15 degrees (see attached mesh image) and Reynolds of approximately 1.1E5. I am using 960 nodes on the airfoil with a thick boundary layer and my timestep should be sufficiently small, being the equivalent of 2000 steps/cycle under the turbine's operating conditions. Under the exact same static setup, I've done a comparison between using the segregated flow solver and coupled implicit one. I've attached an example output of residuals from both of these (turbulent quantities are not displayed because they are always orders of magnitude smaller). Basically, the segregated solver can reach a low residual (say 1E4), but this takes almost 100 inner iterations to accomplish. Conversely, the residuals of the coupled solver completely bottom out and flatline after approximately 50 inner iterations (I set a maximum number of inner iterations of 200 here just to investigate). As my mesh and timestep size should be sufficiently refined for this static airfoil case, I am trying to find the reason for this bizarre residual behaviour and am thinking that it may have something to do with the overset interpolation? But I am really not sure and very much welcome all suggestions. As I imagine someone may propose, I do also monitor the instantaneous forces on the blades as a means of convergence during my actual full turbine simulations. However, having the residuals not converging within timesteps is unsettling and is definitely something I'd like to address. Thanks for any help! Cheers, Dan 

May 2, 2017, 21:43 

#2 
Senior Member
Join Date: Nov 2010
Location: USA
Posts: 1,201
Rep Power: 20 
Why do you need to use overset for this? There's no way to do a sliding mesh?
Do the forces converge before the residuals drop that far? I don't really see any issues with those residual curves, they're just slow. Have you turned the URFs up? If your mesh is good quality you should be able to increase them a lot over the default values. 

May 3, 2017, 16:04 

#3  
New Member
Dan
Join Date: Aug 2016
Location: Montreal, Canada
Posts: 7
Rep Power: 5 
Quote:
Thanks for the reply, I appreciate it. You are indeed correct in that I could just use a sliding mesh for this. I should have explained further: the traditional VAWT here is only the validation case for my project, the actual turbine I will be modelling uses a unique blade path trajectory and thus the overset meshing greatly simplifies the motion definition. Since posting this I have managed to more or less fix the problem I believe (with the help of my postdoc friend ). Still using the Implicit Coupled model, we added in the Coupled Energy model option and changed the fluid model from Constant Density to Ideal Gas. Despite the fact that at this low of a Mach number compressibility effects are of course negligible, this change in model/solver formulation had a huge effect on the residuals trend (see attached plot). Now the residuals fall smoothly below to 1E4 in 4060 inner iterations, at which point the forces are also converged sufficiently. I think the URFs are only applicable to the Segregated solver for pressure and velocity, and after this new discovery with the Coupled solver I'll be sticking with it for now. But thanks for the suggestion! Cheers, Dan 

May 3, 2017, 17:16 

#4 
Senior Member
Join Date: Nov 2010
Location: USA
Posts: 1,201
Rep Power: 20 
Still, 4060 iterations per timestep is very slow. Turning on compressibility making the coupled solver easier to converge isn't much of a surprise, it lets the system be a little less stiff.
What CFL are you using for the solver? 

May 3, 2017, 19:06 

#5 
New Member
Dan
Join Date: Aug 2016
Location: Montreal, Canada
Posts: 7
Rep Power: 5 
If you're referring to using CFL to control the timestep of the transient run, I set my timestep directly as 7.7E5 s, which is the timestep I use in my actual turbine simulations to provide me with 2000 steps per rotor cycle, or equivalently 0.18 degrees of rotor rotation per step.
I see in the Solvers>Coupled Implicit node within the STAR simulation tree there is a Courant Number field, and I have simply left this at the default of 50. Since I specify my timestep manually, I am not sure what this Courant Number specification actually does to be honest. Do you have an explanation? Also, I should mention that the bottleneck on the number of inner iterations now seems to be the force convergence criteria, which I have set quite high for the time being at essentially 0.01% difference between successive inner iteration values. The continuity and momentum residual criteria of 1E4 are usually satisfied 1020 iterations before this. Thanks for the help and discussion. Cheers, Dan Last edited by Dan709; May 3, 2017 at 19:14. Reason: Forgot to include note about residuals 

Tags 
convergence, overset mesh, residuals, vawt, wind turbine 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Moving mesh  Niklas Wikstrom (Wikstrom)  OpenFOAM Running, Solving & CFD  122  June 15, 2014 07:20 
Unstabil Simulation with chtMultiRegionFoam  mbay101  OpenFOAM Running, Solving & CFD  13  December 28, 2013 14:12 
calculation stops after few time steps  sivakumar  OpenFOAM Running, Solving & CFD  7  March 17, 2013 07:37 
SLTS+rhoPisoFoam: what is rDeltaT???  nileshjrane  OpenFOAM Running, Solving & CFD  4  February 25, 2013 05:13 
Differences between serial and parallel runs  carsten  OpenFOAM Bugs  11  September 12, 2008 12:16 