Description Usage Arguments Value References Examples
simulate a thermogram using JMA theory
1 2 |
A |
pre exponential parameters (1/s) |
Ea |
Activation energy (J/mol) |
q |
= rate of analyis (K/min) |
T0 |
= starting temperature of the simulated thermogram expressed in K |
T.end |
= ending temperature of the simulated thermogram expressed in K |
npoints |
desired number of points of the simulate thermogram |
n |
numerical parameter required by JMA model |
T.C = temperature in Celsius
fi=dα/dt * q
α
time.s = time in second
dα dT
1. Vyazovkin S, Chrissafis K, Di Lorenzo ML, et al. ICTAC Kinetics Committee recommendations for collecting experimental thermal analysis data for kinetic computations. Thermochim Acta. 2014;590:1-23. doi:10.1016/j.tca.2014.05.036.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | data <- JMA(A = exp(35),Ea = 120000,q = 50,T0 = -100,T.end = 300,npoints=898,n=2)
require(data.table)
#choose the rates for the simulation of the thermograms
rates=c(0.5,1,2,5,10,20,50)
#first serie of thermograms for all the chosen rate
a<-lapply(rates, function(x) JMA(A=exp(35),Ea=120000,T0=0,T.end=300,q=x,npoints=5000,n=2))
#setup column names
a<-lapply(seq(1,length(a)), function(x) data.table(a[[x]]$time.s,
a[[x]]$T.C, a[[x]]$dadT, rates[[x]]))
lapply(seq(1,length(a)), function(x) setnames(a[[x]],
c("time.seconds","temperature.s","heat.flow","rates") ) )
#create a plot using the function thermo
amaxH <- max(sapply(a, function(x) max(x$heat.flow))) # calculate the max
plot(c(0,300),c(0,amaxH),mytitle="dataset A 120/60 0.66/0.33",
ylab="ExothermicHeatFlow", xlab="Temperature")
lapply(a, function(x) lines(x$temperature.s,x$heat.flow,lwd=3))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.