applications/solvers/surfaceTracking/bubbleInterTrackFoam/solveBulkSurfactant.H

   1
   2 if(CPtr)
   3 {
   4     volScalarField& C = *CPtr;
   5
   6     const dimensionedScalar& D
   7         = interface.surfactant().surfactBulkDiffusion();
   8
   9     scalar ka = interface.surfactant().surfactAdsorptionCoeff().value();
  10     scalar kb = interface.surfactant().surfactDesorptionCoeff().value();
  11     scalar CsInf = interface.surfactant().surfactSaturatedConc().value();
  12
  13     scalarField& Cs =
  14         interface.surfactantConcentration().internalField();
  15
  16     scalarField CP =
  17         C.boundaryField()[interface.aPatchID()].patchInternalField();
  18
  19     const scalarField& Cfs = C.boundaryField()[interface.aPatchID()];
  20
  21     scalarField dn = 1.0/mesh.boundary()[interface.aPatchID()].deltaCoeffs();
  22
  23     if
  24     (
  25         C.boundaryField()[interface.aPatchID()].type()
  26      == fixedGradientFvPatchScalarField::typeName
  27     )
  28     {
  29         fixedGradientFvPatchScalarField& CA =
  30             refCast<fixedGradientFvPatchScalarField>
  31             (
  32                 C.boundaryField()[interface.aPatchID()]
  33             );
  34
  35         CA.gradient() = (ka*kb*Cs - ka*(CsInf-Cs)*Cfs)/D.value();
  36     }
  37     else
  38     {
  39         FatalErrorIn(args.executable())
  40             << "Bulk concentration boundary condition "
  41                 << "at the freeSurface patch is not "
  42                 << fixedGradientFvPatchScalarField::typeName
  43                 << exit(FatalError);
  44     }
  45
  46     solve
  47     (
  48         fvm::ddt(C)
  49       + fvm::div(phi - fvc::meshPhi(rho,U), C, "div(phi,C)")
  50       - fvm::laplacian(D, C, "laplacian(D,C)")
  51     );
  52
  53     CP = C.boundaryField()[interface.aPatchID()].patchInternalField();
  54
  55 //     Info << gMax(CP) << ", "<< gAverage(CP) << endl;
  56 }