JMA: Title Johnson-Mehl-Avrami (JMA)

Description Usage Arguments Value References Examples

Description

simulate a thermogram using JMA theory

Usage

1
2
JMA(A = exp(35), Ea = 120000, q = 50, T0 = -100, T.end = 300,
  npoints = 898, n = 2)

Arguments

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

Value

References

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.

Examples

 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))

takos documentation built on Jan. 13, 2021, 4:11 p.m.

Related to JMA in takos...