applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/setInterfaceDisplacement.H

   1 Info << "\nTime = " << runTime.timeName()
   2     << ", iteration: " << outerCorr << endl;
   3
   4 if (outerCorr < 3 || couplingScheme == "FixedRelaxation")
   5 {
   6     Info << "Current fsi under-relaxation factor: "
   7         << fsiRelaxationFactor << endl;
   8
   9     fluidPatchPointsDisplOld = fluidPatchPointsDispl;
  10
  11     fluidPatchPointsDispl += fsiRelaxationFactor*fsiResidual;
  12 }
  13 else
  14 {
  15     if (couplingScheme == "Aitken")
  16     {
  17         fsiRelaxationFactor =
  18            -fsiRelaxationFactor
  19            *(
  20                 gSum
  21                 (
  22                     fsiResidualOld
  23                    &(fsiResidual - fsiResidualOld)
  24                 )
  25                /(
  26                     gSum
  27                     (
  28                         (fsiResidual - fsiResidualOld)
  29                        &(fsiResidual - fsiResidualOld)
  30                     )
  31                 )
  32            );
  33
  34         fsiRelaxationFactor = mag(fsiRelaxationFactor);
  35
  36         Info << "Current fsi under-relaxation factor (Aitken): "
  37             << fsiRelaxationFactor << endl;
  38
  39         fluidPatchPointsDisplOld = fluidPatchPointsDispl;
  40
  41         fluidPatchPointsDispl +=
  42             fsiRelaxationFactor*fsiResidual;
  43     }
  44 //     else if (couplingScheme == "IQN-ILS")
  45 //     {
  46
  47 //     }
  48 }