applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/calcFsiResidual.H

   1 {
   2     const vectorField& solidPatchDisplacement =
   3         DU.boundaryField()[solidPatchID];
   4
   5     vectorField solidZoneDisplacement
   6         (
   7             stressMesh.faceZones()[solidZoneID]().size(),
   8             vector::zero
   9         );
  10
  11     const label solidPatchStart =
  12         stressMesh.boundaryMesh()[solidPatchID].start();
  13
  14     forAll(solidPatchDisplacement, i)
  15     {
  16         solidZoneDisplacement
  17         [
  18             stressMesh.faceZones()[solidZoneID]
  19            .whichFace(solidPatchStart + i)
  20         ] =
  21             solidPatchDisplacement[i];
  22     }
  23
  24     // Parallel data exchange: collect displacement field on all processors
  25     reduce(solidZoneDisplacement, sumOp<vectorField>());
  26
  27     vectorField fluidZoneDisplacement =
  28         interpolatorSolidFluid.faceInterpolate
  29         (
  30             solidZoneDisplacement
  31         );
  32
  33     vectorField fluidPatchDisplacement
  34         (
  35             mesh.boundary()[fluidPatchID].size(),
  36             vector::zero
  37         );
  38
  39     const label fluidPatchStart =
  40         mesh.boundaryMesh()[fluidPatchID].start();
  41
  42     forAll(fluidPatchDisplacement, i)
  43     {
  44         fluidPatchDisplacement[i] =
  45             fluidZoneDisplacement
  46             [
  47                 mesh.faceZones()[fluidZoneID].whichFace(fluidPatchStart + i)
  48             ];
  49     }
  50
  51     primitivePatchInterpolation fluidPatchInterpolator
  52         (
  53             mesh.boundaryMesh()[fluidPatchID]
  54         );
  55
  56     solidPatchPointsDispl =
  57         fluidPatchInterpolator.faceToPointInterpolate
  58         (
  59             fluidPatchDisplacement
  60         );
  61
  62     fsiResidualOld = fsiResidual;
  63
  64     fsiResidual = solidPatchPointsDispl - fluidPatchPointsDispl;
  65
  66 //     maxFsiResidual =
  67 //        gMax
  68 //         (
  69 //             mag(fsiResidual)
  70 //            /(mag(solidPatchPointsDispl) + SMALL)
  71 //         );
  72
  73 //     Info << "Maximal fsi residual: " << maxFsiResidual << endl;
  74
  75     fsiResidualNorm = ::sqrt(gSum(magSqr(fsiResidual)));
  76
  77     if (outerCorr == 1)
  78     {
  79         initialFsiResidualNorm = fsiResidualNorm;
  80     }
  81
  82     fsiResidualNorm /= initialFsiResidualNorm + SMALL;
  83
  84     Info << "Current fsi residual norm: " << fsiResidualNorm << endl;
  85 }