Question about implement IQN-ILS algorithm in icoFsiElasticNonLinULSolidFoam
icoFsiElasticNonLinULSolidFoam is available in foam-extend-3.0.
After reading the papar“performance of a new partitioned procedure versus a monolithic in fluid-structure interaction”. http://web.mit.edu/kjb/www/Principal...nteraction.pdf I try to implement IQN-ILS algorithm in icoFsiElasticNonLinULSolidFoam according to paper.Blow is my code. 1.Add two line into createInterfaceFields.H Code:
SLList<scalarRectangularMatrix> deltaInput; Code:
else if (couplingScheme == "IQN-ILS") http://b173.photo.store.qq.com/psb?/...KAAgAAAAADAHI! http://user.qzone.qq.com/489509928http://user.qzone.qq.com/489509928/4 But I change the code line from:fluidPatchPointsDispl[id/dim][i]+=(temp[id+i][0]-vVector[id+i][0]); to:fluidPatchPointsDispl[id/dim][i]+=(temp[id+i][0]); Things get a little better,but it still can't meet convergence criteria. [IMG]http://b170.photo.store.qq.com/psb?/V11dkwCX1wKa4a/GjfsQ*slhvGFcYuS46SR.znc5g*L*kXjlU4qiJ.g7og!/b/dEguYGVAKAAA&bo=1wKAAgAAAAADAHI![/IMG] |
Hey,
did you find a solution? I work with this benchmark too and I am really interessted in implementing the new algorithm in the HronTurekFsi Case. I would be very grateful for your help. Best regards, Stephie |
Sorry,
But actually its quite excited to see "image can not be cited from QQ without permission", um. In my primary language. |
I am currently working with this case, Stephie, keep us posted about your findings! :D Are you using the FSI-PFS-2A setting? I am using M. Breuers data as reference, what data are you using? Could you describe your case setting? Best Regards /Johan
FYI: IQN-ILS is implemented in icoFsiElasticNonLinULSolidFoam. |
Hello Dears
i am using this solver for my case. in the tutorial case of this solver the properties of fluid and solid set to fluid density = 1 kg/m3 kinemaic viscosity = 10e-3 solid density = 10 kg/m3 what is the meaning of this properties.? when i set my real properties for blood and vessel blood density = 1000 kg/m3 vessel density = 1200 kg/m3 my solution crashed by " complex eigenvalue detected for tensor ...." error. best regards |
1 Attachment(s)
Hi I know this is an old post. But could your give me a link to the implementation of IQN-ILS in icoFsiElasticNonLinULSolidFoam. I am not able to find it in foam extend 3.1 as that portion is commented.
|
Quote:
I have implemented my own since the provided use explicit Q, i prefer Q free implementation using householder method. |
All times are GMT -4. The time now is 15:45. |