<XDFCONSTANT uniqueid="0x0">
    <title>KRKTE</title>
    <description>Conversion of relative fuel mass (rk) to effective injector on time (te)</description>
    <CATEGORYMEM index="0" category="29" />
    <EMBEDDEDDATA mmedtypeflags="0x02" mmedaddress="0x1C2BC" mmedelementsizebits="16" />
    <units>cc/min</units>
    <outputtype>2</outputtype>
    <decimalpl>0</decimalpl>
    <datatype>0</datatype>
    <unittype>0</unittype>
    <DALINK index="0" />
    <MATH equation="307432.000000 * X">
      <VAR id="X" />
    </MATH>
  </XDFCONSTANT>

#include <stdio.h>
#include <math.h>

// Injector information : From Datasheet
//    Static flow rate @ 3.00 BAR, medium N-Heptane: 213,2 g/min
// { 43.5 PSI (3.00 BAR)   = 267 cc/min    gMin = fuelNHeptane*ccMin }    

void main(void)
{
   double i;
   double barConversion         = 14.5038;                 // conversion of bar   : psi
   double lbhrConversion        = 0.095200;                // conversion of lb/hr : cc/min

   double rho0Air               = 1.293;                   // Air Density (@ 0C and 1013 hPa) [g/dm^3]

   double engineCC              = 3586.0;	               // Engine Displacement             [dm^3]
   double numCyl                = 8;			           // Number of cylinders
   double Vhzyl                 = (engineCC/numCyl)/1000;  // Cylinder Displacement           [dm^3]

   double ccMin                 = 267;                     // cc/min:  fuel injector cc/min @ given psi in Injector Datasheet
   double origFuelPressurePsi   = 43.5;                    // psi:     fuel injector cc/min @ given psi in Injector Datasheet
   double fuelPressurePsi       = 50.0;                    // *actual* fuel injector cc/min pressure used in this application (check workshop manual or measure it)

   double fuelPressureBar       = fuelPressurePsi / barConversion;
   double fuelInj_ccMin_at_1psi = ccMin/origFuelPressurePsi;
   double final_ccMin           = fuelPressurePsi * fuelInj_ccMin_at_1psi;    

   double Lst                   = 14.7;		               // afr for stoich
   double McuClockSpeed			= 40.0;					   // Frequency (MHz)
   double McuClockAdjust		= 0.00000041676;
   double McuAdaption           = McuClockSpeed*McuClockAdjust;   
   double Normmk                = McuAdaption;			   // 0.00001667 = (min/ms Adaptation of units @ 40MHz )

   double fuelNHeptane          = 0.75;	                   // 0.694, 0.71 or 0.75?
   double nHeptaneCorrection    = 1.05;
   double Qstat                 = fuelNHeptane*final_ccMin;

   double start_psi             = 43.5;
   double stop_psi              = 60.0;
   double psi_increment         = 0.1;

   printf("Cubic centimeters per minute flow rate unit: cc/min\n");
   printf("engineCC : %2.2f (Engine Displacement [dm^3])\n",engineCC);
   printf("numCyl   : %2.2f (Number of cylinders)\n",numCyl);
   printf("Vhzyl    : %2.2f (Cylinder Displacement [dm^3])\n\n",Vhzyl);
   printf("\tGiven Spec Static flow rate @ %2.2f PSI ( %2.2f BAR) is %2.2f cc/min (%2.2f lb/hr) : spec of injectors given\n",origFuelPressurePsi, origFuelPressurePsi / barConversion, ccMin, ccMin*lbhrConversion );
   printf("\tCalculated Static flow rate @ %2.2f PSI ( %2.2f BAR) is %2.2f cc/min (%2.2f lb/hr) : for *this* application (derived from given spec)\n\n--\n",fuelPressurePsi,fuelPressureBar, final_ccMin, final_ccMin*lbhrConversion);

   // display calculated KRKTE @ different fuel pressures 
   for(i=start_psi;i<stop_psi;i+= psi_increment)       // range of psi to show conversions for
   {
       printf("\tDerived  : %2.2f cc/min (%2.2f LB/HR) @ %2.1f PSI (%2.2f BAR) : ", (i * fuelInj_ccMin_at_1psi),(i * fuelInj_ccMin_at_1psi)*lbhrConversion,i, i/barConversion);
       final_ccMin    = i * fuelInj_ccMin_at_1psi;
       Qstat          = fuelNHeptane*final_ccMin;
	   printf("KRKTE = %2.4f ms/%%\n", (rho0Air*Vhzyl ) /	(100 * Lst * Normmk * nHeptaneCorrection * Qstat ) );
   }

}