Learn R Programming

vein (version 1.1.3)

emis_evap2: Estimation of evaporative emissions 2

Description

emis_evap performs the estimation of evaporative emissions from EMEP/EEA emission guidelines with Tier 2.

Usage

emis_evap2(
  veh,
  name,
  size,
  fuel,
  aged,
  nd4,
  nd3,
  nd2,
  nd1,
  hs_nd4,
  hs_nd3,
  hs_nd2,
  hs_nd1,
  rl_nd4,
  rl_nd3,
  rl_nd2,
  rl_nd1,
  d_nd4,
  d_nd3,
  d_nd2,
  d_nd1
)

Value

dataframe of emission estimation in grams/days

Arguments

veh

Total number of vehicles by age of use. If is a list of 'Vehicles' data-frames, it will sum the columns of the eight element of the list representing the 8th hour. It was chosen this hour because it is morning rush hour but the user can adapt the data to this function

name

Character of type of vehicle

size

Character of size of vehicle

fuel

Character of fuel of vehicle

aged

Age distribution vector. E.g.: 1:40

nd4

Number of days with temperature between 20 and 35 Celsius degrees

nd3

Number of days with temperature between 10 and 25 Celsius degrees

nd2

Number of days with temperature between 0 and 15 Celsius degrees

nd1

Number of days with temperature between -5 and 10 Celsius degrees

hs_nd4

average daily hot-soak evaporative emissions for days with temperature between 20 and 35 Celsius degrees

hs_nd3

average daily hot-soak evaporative emissions for days with temperature between 10 and 25 Celsius degrees

hs_nd2

average daily hot-soak evaporative emissions for days with temperature between 0 and 15 Celsius degrees

hs_nd1

average daily hot-soak evaporative emissions for days with temperature between -5 and 10 Celsius degrees

rl_nd4

average daily running losses evaporative emissions for days with temperature between 20 and 35 Celsius degrees

rl_nd3

average daily running losses evaporative emissions for days with temperature between 10 and 25 Celsius degrees

rl_nd2

average daily running losses evaporative emissions for days with temperature between 0 and 15 Celsius degrees

rl_nd1

average daily running losses evaporative emissions for days with temperature between -5 and 10 Celsius degrees

d_nd4

average daily diurnal evaporative emissions for days with temperature between 20 and 35 Celsius degrees

d_nd3

average daily diurnal evaporative emissions for days with temperature between 10 and 25 Celsius degrees

d_nd2

average daily diurnal evaporative emissions for days with temperature between 0 and 15 Celsius degrees

d_nd1

average daily diurnal evaporative emissions for days with temperature between -5 and 10 Celsius degrees

References

Mellios G and Ntziachristos 2016. Gasoline evaporation. In: EEA, EMEP. EEA air pollutant emission inventory guidebook-2009. European Environment Agency, Copenhagen, 2009

Examples

Run this code
if (FALSE) {
data(net)
PC_G <- c(33491,22340,24818,31808,46458,28574,24856,28972,37818,49050,87923,
          133833,138441,142682,171029,151048,115228,98664,126444,101027,
          84771,55864,36306,21079,20138,17439, 7854,2215,656,1262,476,512,
          1181, 4991, 3711, 5653, 7039, 5839, 4257,3824, 3068)
veh <- data.frame(PC_G = PC_G)
pc1 <- my_age(x = net$ldv, y = PC_G, name = "PC")
ef1 <- ef_evap(ef = "erhotc",v = "PC", cc = "<=1400", dt = "0_15", ca = "no")
dfe <- emis_evap2(veh = pc1,
                 name = "PC",
                 size = "<=1400",
                 fuel = "G",
                 aged = 1:ncol(pc1),
                 nd4 = 10,
                 nd3 = 4,
                 nd2 = 2,
                 nd1 = 1,
                 hs_nd4 = ef1*1:ncol(pc1),
                 hs_nd3 = ef1*1:ncol(pc1),
                 hs_nd2 = ef1*1:ncol(pc1),
                 hs_nd1 = ef1*1:ncol(pc1),
                 d_nd4 = ef1*1:ncol(pc1),
                 d_nd3 = ef1*1:ncol(pc1),
                 d_nd2 = ef1*1:ncol(pc1),
                 d_nd1 = ef1*1:ncol(pc1),
                 rl_nd4 = ef1*1:ncol(pc1),
                 rl_nd3 = ef1*1:ncol(pc1),
                 rl_nd2 = ef1*1:ncol(pc1),
                 rl_nd1 = ef1*1:ncol(pc1))
lpc <- list(pc1, pc1, pc1, pc1,
            pc1, pc1, pc1, pc1)
dfe <- emis_evap2(veh = lpc,
                 name = "PC",
                 size = "<=1400",
                 fuel = "G",
                 aged = 1:ncol(pc1),
                 nd4 = 10,
                 nd3 = 4,
                 nd2 = 2,
                 nd1 = 1,
                 hs_nd4 = ef1*1:ncol(pc1),
                 hs_nd3 = ef1*1:ncol(pc1),
                 hs_nd2 = ef1*1:ncol(pc1),
                 hs_nd1 = ef1*1:ncol(pc1),
                 d_nd4 = ef1*1:ncol(pc1),
                 d_nd3 = ef1*1:ncol(pc1),
                 d_nd2 = ef1*1:ncol(pc1),
                 d_nd1 = ef1*1:ncol(pc1),
                 rl_nd4 = ef1*1:ncol(pc1),
                 rl_nd3 = ef1*1:ncol(pc1),
                 rl_nd2 = ef1*1:ncol(pc1),
                 rl_nd1 = ef1*1:ncol(pc1))
}

Run the code above in your browser using DataLab