FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvOptions;
}

electron-beam-heating
{
  type            scalarCodedSource;
  active          yes;
  name            heating;

  scalarCodedSourceCoeffs
  {
    selectionMode   cellZone;
    cellZone        heate;
    fields          ( T );

    codeInclude     #{ 
    #include "IFstream.H"
    #include "OFstream.H"
    #};

    codeCorrect     #{ #};

    codeAddSup
    #{

      const scalar rhom = readScalar(coeffs_.lookup("rho"));
      const scalar Cp = readScalar(coeffs_.lookup("Cp"));
      const scalar q0 = readScalar(coeffs_.lookup("power"));
      const scalar z0 = readScalar(coeffs_.lookup("z0"));
      const scalar x0 = readScalar(coeffs_.lookup("x0"));
      const scalar zte = readScalar(coeffs_.lookup("zte"));
      const scalar re = readScalar(coeffs_.lookup("re"));
      const word zone_name(coeffs_.lookup("target"));
      scalarField& src = eqn.source();
      const auto& cell_zone = mesh_.cellZones()[zone_name];
        
        Info<<"cell_zone:"<<endl;
        Info<<cell_zone<<endl;
        ofstream test;
        test.open("test.txt");

      forAll(cell_zone, i) 
      {

        const label cell_i = cell_zone[i];

        Info<<"cell_i"<<endl;
        Info<<cell_i<<endl;


        const scalar z = mesh_.C()[cell_i].y();
        const scalar x = mesh_.C()[cell_i].x();
        const scalar xa =mag( x - x0);
        const scalar za =mag( z - z0);
        const scalar pre1 = (xa*xa)/(re*re);

        Info<<"Abstand in z"<<endl;
        Info<<za<<endl;
        Info<<"z-Koordinate"<<endl;
        Info<<z<<endl;
        Info<<"Faktor von e"<<endl;
        Info<<pre1<< endl; 

         const scalar qvo = q0 * (-1.4 * ((za * za) / (zte * zte)) + 1.4 * (za/zte) + 0.65);
        Info<<"qvo"<<endl;
        Info<<qvo<<endl;

        const scalar q = qvo*exp(-pre1); // i.e. q0*exp(-distance)

        Info<<"q"<<endl;
        Info<<q<<endl;

        src[cell_i] += -q/(rhom*Cp);
        Info<<"q/rho/cp"<<endl;
        Info<<src[cell_i]<<endl;

        //test<<"cell_i:"<<" "<<cell_i<<";"<<" Abstand in z von Brennpunkt:"<<za<<";"<<" z-Koordinate der Zelle:"<<z<<";"<<" Faktor von e:"<<pre1<<";"<<" Schwächung von Q in z-Richtung:"<<qvo<<";"<<" Abgeschwächtes Q an bei der Zelle:"<<q<<";"<<" Vorwert für die Temperatur Q/Cp/rho:"<<src[cell_i]<<";"<<"\n"<<endl;
       
       }

    #};

    codeSetValue    #{ #};

    code
    #{
      $codeAddSup
    #};
  }

  heatingCoeffs
  {
    target    heate;  // cellZone name
    power     0.1;   // W
    z0        1e-3;     // m
    x0        0.5e-3; 
    rho       10070;   // kg/m3
    Cp        120;    // J/kg/K
    zte       239.4381122e-6;
    re        0.5e-3;
  }
}
// vi: set ft=foam et sw=2 ts=2 sts=2: