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

Question about implement IQN-ILS algorithm in icoFsiElasticNonLinULSolidFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 15, 2014, 09:02
Question Question about implement IQN-ILS algorithm in icoFsiElasticNonLinULSolidFoam
  #1
New Member
 
Jinbin Ge
Join Date: Mar 2013
Posts: 6
Rep Power: 4
gejiabin is on a distinguished road
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;
SLList<scalarRectangularMatrix> deltaOutput;
2.The main code is in setInterfaceDisplacement.H
Code:
  else if (couplingScheme == "IQN-ILS")
    {
     Info<< "IQN-ILS iteration: "<<endl;
     //if(outerCorr>20)
     //{
     // deltaInput.removeHead();
     // deltaOutput.removeHead();
    // }
     //-------------------------------------------------------------------------
     int numColumn(0);
     int dim=2;//3
     int numDegreeOfFreedom( dim*mesh.boundaryMesh()[fluidPatchID].nPoints() );
     scalarRectangularMatrix vMatrix(numDegreeOfFreedom,deltaInput.size()+1,0);
     scalarRectangularMatrix wMatrix(numDegreeOfFreedom,deltaInput.size()+1,0);
     scalarRectangularMatrix vVector(numDegreeOfFreedom,1,0);
     scalarRectangularMatrix temp(numDegreeOfFreedom,1,0);
     
     //1.Assemble vMatrix  
     for(int id=0;id<numDegreeOfFreedom;id=id+dim)
      for(int i=0;i<dim;i++)
       temp[id+i][0]=(fsiResidualOld[id/dim][i]-fsiResidual[id/dim][i]);
     //
     forAllIter(SLList<scalarRectangularMatrix>,deltaInput,iter)
     { 
       (*iter)=(*iter)+temp;
       for(int i=0;i<numDegreeOfFreedom;i++)    
        vMatrix[i][numColumn]=(*iter)[i][0];
       numColumn+=1;
     } 
     //
     deltaInput.append(temp);
     for(int i=0;i<numDegreeOfFreedom;i++) 
      vMatrix[i][numColumn]=temp[i][0];
     
     //2.Assemble vVecotr
     for(int id=0;id<numDegreeOfFreedom;id=id+dim)
      for(int i=0;i<dim;i++)
       vVector[id+i][0]=-fsiResidual[id/dim][i];
     
     //3.Assemble wMatrix 
     numColumn=0;
     for(int id=0;id<numDegreeOfFreedom;id=id+dim)
      for(int i=0;i<dim;i++)
       temp[id+i][0]=(fluidPatchPointsDisplOld[id/dim][i]-fluidPatchPointsDispl[id/dim][i]);
     //
     forAllIter(SLList<scalarRectangularMatrix>,deltaOutput,iter)
     {  
       (*iter)=(*iter)+temp;
       for(int i=0;i<numDegreeOfFreedom;i++)    
        wMatrix[i][numColumn]=(*iter)[i][0];
       numColumn+=1;
     } 
     //
     deltaOutput.append(temp);
     for(int i=0;i<numDegreeOfFreedom;i++) 
      wMatrix[i][numColumn]=temp[i][0];
       //4.Update fluidPatchPointsDispl
     fluidPatchPointsDisplOld = fluidPatchPointsDispl;
     multiply(temp,wMatrix,SVDinv(vMatrix),vVector);
     for(int id=0;id<numDegreeOfFreedom;id=id+dim)
      for(int i=0;i<dim;i++)
       fluidPatchPointsDispl[id/dim][i]+=(temp[id+i][0]-vVector[id+i][0]) ;
     
      Info<<"end:"<<endl;
     //-------------------------------------------------------------------------
     }
}
But the new solver can't convergence for the benchmark proposed in the Hron & Turek paper.









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]

Last edited by gejiabin; April 17, 2014 at 03:31.
gejiabin is offline   Reply With Quote

Old   June 28, 2015, 08:05
Default
  #2
Member
 
Stephanie
Join Date: Feb 2015
Location: Magdeburg, Germany
Posts: 62
Rep Power: 2
stephie is on a distinguished road
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
stephie is offline   Reply With Quote

Old   June 28, 2015, 12:00
Default
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 676
Rep Power: 8
sharonyue is on a distinguished road
Sorry,

But actually its quite excited to see "image can not be cited from QQ without permission", um. In my primary language.
__________________
Im the translator of OpenFOAM User Guide Chinese Edition! But always newbie on CFD.
Never too old to learn CFD. Worship these CFD experts.
sharonyue is offline   Reply With Quote

Old   July 1, 2015, 05:44
Default
  #4
Member
 
Johan Lorentzon
Join Date: Mar 2009
Location: Lunds University, Sweden
Posts: 71
Rep Power: 14
pi06jl6 will become famous soon enough
I am currently working with this case, Stephie, keep us posted about your findings! 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.
pi06jl6 is offline   Reply With Quote

Reply

Tags
elasticnonlinulsolidfoam, icofsi, iqn-ils

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
On the pimpleFoam solution algorithm vkrastev OpenFOAM 1 February 1, 2011 09:25
Question on Pressure and Velocity Fields with SIMPLE algorithm Phiper Main CFD Forum 0 November 18, 2010 19:29
Question regarding the F2F particle tracking algorithm implementation. andrewryan OpenFOAM 2 October 22, 2009 01:55
Question to SIMPLE algorithm by Patankar's book. Rafal Main CFD Forum 2 December 21, 2005 04:48
SIMPLER Algorithm question Erik Main CFD Forum 1 May 23, 2004 03:57


All times are GMT -4. The time now is 09:57.