applications/solvers/solidMechanics/viscoElasticSolidFoam/createFluid.H

   1
   2     fileName casePath(args.additionalArgs()[0]);
   3     fileName fluidRootDir = casePath.path();
   4     fileName fluidCaseDir = casePath.name();
   5
   6     Info<< "Fluid: " << fluidRootDir << " " << fluidCaseDir << endl;
   7
   8     Time fluidRunTime
   9     (
  10         Time::controlDictName,
  11         fluidRootDir,
  12         fluidCaseDir
  13     );
  14
  15     fluidRunTime.setDeltaT(runTime.deltaT().value());
  16
  17     fvMesh fluidMesh
  18     (
  19         IOobject
  20         (
  21             fvMesh::defaultRegion,
  22             fluidRunTime.timeName(),
  23             fluidRunTime
  24         )
  25     );
  26
  27     label fluidInterfacePatchID =
  28         fluidMesh.boundaryMesh().findPatchID("interface");
  29
  30     if(fluidInterfacePatchID == -1)
  31     {
  32         FatalErrorIn(args.executable())
  33             << "Can't find patch with name interface" << abort(FatalError);
  34     }
  35
  36     patchToPatchInterpolation interpolator
  37     (
  38         fluidMesh.boundaryMesh()[fluidInterfacePatchID],
  39         mesh.boundaryMesh()[interfacePatchID],
  40         intersection::VISIBLE
  41     );
  42
  43     scalarField error =
  44         mag
  45         (
  46             mesh.C().boundaryField()[interfacePatchID]
  47           - interpolator.faceInterpolate(fluidMesh.C().boundaryField()[fluidInterfacePatchID])
  48         );
  49
  50     Info << "Patch-to-patch interpolation error, max: " << max(error) << ", avg: " << average(error) << endl;