|
[Sponsors] |
Why do oscilations appear when reducing the time step? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 28, 2015, 07:35 |
Why do oscilations appear when reducing the time step?
|
#1 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
Hi,
I am simulating a flow past a circular cylinder, and according to references, the fluid flow establishes a oscilating pattern. Using a certain time step, there is a deviation (5%) in the frequency of my simulation and the claimed frequency in references. For this reason I have reduced the time step in order to obtain better accuracy. But, when reducing the time step of the simulation, several small oscilations appear additionally to the main oscilation expected. Why are these small oscialtions appearing? Which is the physics behind these oscilations? I was wondering if there could be an error in the implementation of the algorithm. Best regards, Hector. |
|
May 28, 2015, 08:17 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
Quote:
On the other hand, reducing the computational time step means an increasing in the Nyquist cut-off pi/dt, therefore you velocity field can show more frequencies. But this is relevant in DNS/LES formulation, what are you using? |
||
May 28, 2015, 09:35 |
|
#3 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
Hi Filippo,
Thanks for your reply. I do not use any turbulence model. The simulation is for laminar flow. No DNS or LES formulation is used. Then, according to your statement, it could be an error in the integration method. The algorithm contains a relaxation factor for the integration method. Right now, I have set this relaxation factor to 0.6666, which means that 0.666 times of the velocity at time n is take into account and 0.333 times of the velocity at time n+1 is take into account. Based on your suggestion, the problem could be in the value used for this relaxation factor. Best regards, Hector. |
|
May 28, 2015, 11:12 |
|
#4 |
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 |
What is your Reynolds number?
|
|
May 28, 2015, 11:15 |
|
#5 |
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 |
Oh, and given your time step, what is your maximum CFL number?
|
|
May 28, 2015, 11:20 |
|
#6 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
||
May 28, 2015, 12:09 |
|
#7 |
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 |
My thought was that he might be using an implicit method.
Edit: Or something with dual time stepping, in which case the max CFL for the outer loop is the number I'm interested in. |
|
May 28, 2015, 12:47 |
|
#8 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
The Reynolds number is 100 (Which I consider is laminar).
The maximum time step for the algorithm that fulfills the CFL condition is 0.001 (both convective and diffusive limitation), and initially I am using 0.0005 as time step for the simulation. With this value, no small oscilations happens. When using 0.0001 (even more reduced value) as time step for the simulation is when the small osicilations appear. |
|
May 28, 2015, 12:53 |
|
#9 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
I am using an implicit method in the equation of the pressure. The value of the relaxation factor for the pressure equation I am using is 0.666.
The type of algortithm I am using is the so called Characteristics Based Split algorithm (CBS), that uses a split scheme with a correction step once the pressure is calculated. So, both relaxation factors I am using: One for the velocity (explicit equation): 0.6666 Another for the pressure (implicit equation): 0.6666 |
|
May 28, 2015, 12:59 |
|
#10 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
could you post some results showing the comparison for both time step?
A further issue to check is the threshold for the iterative solvers, it could be useful to diminuish it according to the time step |
|
May 28, 2015, 14:37 |
|
#11 | |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34 |
Quote:
Also is it a collocated solver? Because pressure velocity may decouple at small time steps and could give you oscillations or checkerboarding. |
||
May 28, 2015, 15:29 |
|
#12 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
I am attaching the results I am getting:
- VelocityYDirectionTimeStep1 When using a time step of 0.0005 - VelocityYDirectionTimeStep2 When using a time step of 0.0002 As you can see, when reducing the time step, some oscilations appear in between the main oscilations of the pattern flow. |
|
May 28, 2015, 15:40 |
|
#13 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
Try reducing all thresholds in the iterative solvers.
|
|
May 28, 2015, 17:36 |
|
#14 |
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 |
My area of knowledge is centered around implicit compressible time marching methods. So I'm not sure if that carries over to what you're doing. However, a CFL of 0.5 may be too large to insure accurate time marching for a method with an explicit component. It may be stable but possibly not accurate. One may need to go for dual time stepping, fully implicit, higher order time discretization, or lower CFL to achieve higher time accuracy.
That being said, an Re of 100 is low, so unless your grid is really coarse, I don't see what what the problem is. Hopefully one of the other posts above helps you out. Good luck, BTW, I am assuming that the boundary conditions are implemented correctly and that grid is not being subdivided. If the grid is being subdivided and distributed to more than one processor then the B.C.s at the grid edges may be lagged. In which case dual time stepping might be required. But, I'm not sure if that applies to your algorithm. |
|
May 29, 2015, 07:41 |
|
#15 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
Quote:
the shape of oscillations makes me suppose that you should try to work on the residuals, try two order of magnitude lower. Do you check similar problem by refining the spatial grid? |
||
May 29, 2015, 09:56 |
|
#16 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
This problem does not appear in a coarse grid.
It has appeared when I have used a finer grid. I understand that your statement "to work on the residuals," refers to the tolerance (threshold) of the iterative solver, doesn't it? |
|
May 29, 2015, 10:21 |
|
#17 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
but you have simultaneously decreased both the time step and the grid sizes?
|
|
May 29, 2015, 10:40 |
|
#18 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
Yes I have drecreased both time step and grid size simultaneously, since the shorter the grid the more reduced the time step should be used. Otherwise the solution does not converge.
|
|
May 29, 2015, 11:12 |
|
#19 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
||
May 29, 2015, 16:16 |
|
#20 |
Senior Member
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 243
Rep Power: 16 |
I have kept constant the dt/h ratio between grids.
|
|
Tags |
oscilations |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
negative min.(alpha1) value interFoam | Arjun Jayakumar | OpenFOAM | 11 | December 21, 2019 10:59 |
High Courant Number @ icoFoam | Artex85 | OpenFOAM Running, Solving & CFD | 11 | February 16, 2017 13:40 |
Rapidly decreasing deltaT for interDyMFoam | chrisb2244 | OpenFOAM Running, Solving & CFD | 3 | July 1, 2014 16:40 |
InterFoam negative alpha | karasa03 | OpenFOAM | 7 | December 12, 2013 03:41 |
calculation diverge after continue to run | zhajingjing | OpenFOAM | 0 | April 28, 2010 04:35 |