emis_evap: Estimation of evaporative emissions

Description Usage Arguments Value References Examples

View source: R/emis_evap.R

Description

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

Usage

1
2
3
emis_evap(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)

Arguments

veh

Total number of vehicles by age of use. If is a lsit 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 celcius degrees

nd3

Number of days with temperature between 10 and 25 celcius degrees

nd2

Number of days with temperature between 0 and 15 celcius degrees

nd1

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

hs_nd4

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

hs_nd3

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

hs_nd2

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

hs_nd1

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

rl_nd4

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

rl_nd3

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

rl_nd2

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

rl_nd1

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

d_nd4

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

d_nd3

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

d_nd2

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

d_nd1

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

Value

dataframe of emission estimation in grams/days

References

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
{
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_evap(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_evap(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))
}

ibarraespinosa/vein documentation built on Dec. 5, 2018, 4:28 a.m.