CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Aitken adaptive under-relaxation for FSI

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 3, 2010, 09:27
Default Aitken adaptive under-relaxation for FSI
  #1
Member
 
Wolfgang W.
Join Date: Nov 2009
Location: Switzerland
Posts: 57
Rep Power: 8
WiWo is on a distinguished road
Hello everyone,

I'm currently using icoFsiFoam extended with an outer iteration loop as a solver for strongly coupled FSI systems (blood - vessel interaction in my case). By setting some under-relaxation factors I could get the solver to handle problems even close to a density ratio (solid/fluid) of 1. Anyway, this was at enormous expense of computational time and I learned that one can do much better by implementing adaptive under-relaxation. Aitken's under-relaxation algorithm was mentioned several times in this context (Tukovic ; Degroote et al, Computers and Structures 86 (2008) 2224 - 2234 ; Kuettler et al, Comput Mech (2008) 43:61-72).

I implemented the algorithm, which appears to be quite simple using Usolid (the displacement of the solid mesh) as parameter.
Still, the algorithm is yielding obviously wrong results ranging from negative values to values > 1 for the relaxation factor. I checked my code many times and then re-implemented it two times from scratch using analytically equivalent formulations of Aitken's algorithm. All three versions yield the exact same (wrong) results.

Looking at the formulation of Aitken's algorithm more closely, I see no reason why result of the process should remain in the interval [0, 1] anyway. It seems to me that the previous value of the relaxation factor is essentially multiplied by a term that can have values ranging from -infinity to +infinity. Of course I'm not assuming that something is wrong in an algorithm as widely used as this one ... but it seems that I'm clearly missing a point here.

So, my question is if anybody has any experience using adaptive under-relaxation or even the Aitken algorithm and could provide me with some deeper insight. Moreover, i was wondering if maybe somebody has already implemented this algorithm for an FSI coupling scheme. Maybe I'm just working on the totally wrong parameter (Usolid)?

I would greatly appreciate your help on this issue!

Best,
Wolfgang
WiWo is offline   Reply With Quote

Old   February 4, 2010, 08:22
Default
  #2
New Member
 
Thomas Gallinger
Join Date: Mar 2009
Posts: 28
Rep Power: 8
thomas is on a distinguished road
Dear Wolfgang,

Aitken relaxation is a form of fixed-point iteration working at the interface.

If I understand right, you are working with the whole displacement field of the solid region. But you have to use the interface displacement field as input for Aitken's formula.

Hope this helps!

Best regards
Thomas
thomas is offline   Reply With Quote

Old   February 22, 2010, 08:50
Default
  #3
Member
 
Wolfgang W.
Join Date: Nov 2009
Location: Switzerland
Posts: 57
Rep Power: 8
WiWo is on a distinguished road
Hi Thomas,

Thank you for your kind support! ... and sorry for my delayed response. My time was recently consumed entirely by other urgent business thus I could just return to the Aitken problem today.
I just had a first run with the Aitken algoritm working only on the interface region. It seems to perform better than the previous version but after some iterations the relaxation factor still starts to play crazy (jumping outside the interval [0 1]).

What i did exactly is implementing the formula for the Aitken algorithm as given by Tukovic for the inner loop of a strongly coupled FSI solver in order to increase the efficiency. Anyway, looking at the formula i don't see a mathematical reason why it should yield results between 0 and 1 - could you explain that a little bit more closely?

Thanks for your time and efforts!

Best regards,
Wolfgang
WiWo is offline   Reply With Quote

Old   February 22, 2010, 10:47
Default
  #4
New Member
 
Thomas Gallinger
Join Date: Mar 2009
Posts: 28
Rep Power: 8
thomas is on a distinguished road
Dear Wolfgang,

I do not know, if the relaxation parameter you get by using Aitken's formula is restricted to values between 0 and 1.

Maybe the parameter is calculated right, but something with your field solvers is wrong. I would suggest to first test your solver using constant underrelaxation and a relaxation parameter of maybe 0.1. If the problem converges, you can replace the relaxation parameter then with Aitken's method.

As a starting point for the implementation I suggest the PhD thesis of Mok. You can find it at http://www.ibb.uni-stuttgart.de/index.html - Publikationen. It gives a nice version of Aitkens formula.

Best regards
Thomas
thomas is offline   Reply With Quote

Old   August 5, 2013, 23:04
Default
  #5
New Member
 
Vincent
Join Date: Aug 2013
Location: Singapore
Posts: 7
Rep Power: 4
VincentChai is on a distinguished road
Dear Wolfgang,

May I know is the icoFsiFoam extended you used is available in the openFOAM 1.6 ext? or can you advise me how to make the icoFsiFoam become a strongly coupled FSI solver?
Thanks

Regards,
Vincent
VincentChai is offline   Reply With Quote

Old   January 4, 2016, 02:49
Default
  #6
New Member
 
anand sudhi
Join Date: Sep 2015
Posts: 16
Rep Power: 2
anandsudhi is on a distinguished road
Hi Wolfgang,

I am working on a similar problem with strong coupling and I have the same behavior.
I would like to know if you have found a solution to your problem or is there no definite interval restriction for the relaxation factors of Aitken algorithm.

Thank you
anandsudhi is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Relaxation and convergence sammi Phoenics 0 March 20, 2008 04:32
CFX-10 - Adaptive Time Stepping James Date CFX 3 October 17, 2005 22:15
relaxation parameter rvndr Main CFD Forum 0 May 5, 2004 12:34
Question on adjusting relaxation factor CFD Rookie Main CFD Forum 3 January 26, 2004 15:37
Relaxation factor Moon CD-adapco 1 June 13, 2003 11:13


All times are GMT -4. The time now is 02:21.