In this paper, an ASPEN plus version 10.0 (Aspen Technology Inc., Cambridge, MA, USA) process simulator is utilised to simulate the ISHSE process and DST process. This simulator is selected since it enables the generation of accurate mass and energy balance data in a timely fashion. This simulation software also facilitates the computation of the flow rates, compositions, and physical characteristics of the input and output streams which will constitute the basis for subsequent energetic and environmental performance assessments. The processes of DST and ISHSE for biodiesel production to be investigated are illustrated in
Figure 1 and
Figure 2 respectively. In the alternative pathways investigated, the biodiesel production capacity is fixed to a production capacity of approximately 6120 t per y with the plants operating for 7200 h per y. In
Figure 1 and
Figure 2, the ‘dashed’ blocks represent the overall processes and the solid blocks denote the individual processes.
Figure 1 shows that the overall DST biodiesel production process involves a supercritical transesterification reaction of triglycerides in the presence of methanol at a temperature of 280 °C and a pressure 28 MPa. At the conclusion of the transesterification reaction the excess and unreacted methanol is recovered via vaporisation from the FAME and glycerol products. Further purification of the FAME product to remove the glycerol and any remaining methanol impurity is subsequently achieved via the distillation process shown in
Figure 1 using a distillation column.
Figure 2 shows that the overall ISHSE biodiesel production process involves an initial water-subcritical hydrolysis of triglycerides at a temperature of 270 °C and a pressure 7 MPa and then the generated fatty acids are separated via decantation, with residual polar aqueous phase containing glycerol and water hydrolysis being left. The separated fatty acids are esterified in the presence of methanol under supercritical conditions of 250 °C and 8 MPa of methanol. At the conclusion of the supercritical esterification reaction, the FAME product is purified by separating the FAMEs from the unreacted methanol and water via vaporisation of methanol and water impurities. Recovery of methanol from the residual methanol-water mixture is undertaken using the distillation columns as shown in
Figure 2. The following mass and energy balance equations at steady state were used and solved using ASPEN plus for each of the individual processes of the DST and ISHSE biodiesel production processes shown in
Figure 1 and
Figure 2,
where,
is the net heat transfer flow rate in kJ/s into the system in the control volume of the individual process, per s;
hi is the specific enthalpy of the
ith stream, in kJ/kg;
is the rate of net work done by the system in the control volume of the individual process in kJ/s; and
is the mass flow rate of the
ith stream, in kg/s.
To minimise the cooling and heating duties of both high temperature and high pressure biodiesel production processes, Aspen Energy Analyser
® (Aspen Technology Inc., Cambridge, MA, USA) is used to perform heat integration and to design a possible heat exchanger network based on classic pinch analysis methods. These pinch analysis methods incorporate the assessment of heat flow rate inequalities and stream splitting rules via tick-off heuristics [
20]. In utilising Aspen Energy Analyser
® for the reduction of heating and the cooling duties, the procedure proposed by Sadhukha et al. [
21] is utilised. Necessary thermodynamic data highlighting the supply temperature and target temperature of relevant streams are extracted from simulation output sheet generated by ASPEN plus. Using these data the composite curve is constructed to aid in the determination of the minimum energy consumption target since the overlap of heat availability (hot composite curve) and the heat requirement (cold composite curve) provides an indication of the maximum possible process heat recovery [
22]. The determination of the possible heat recovery enables the estimation of the remaining heating requirement (called the minimum hot Q(H)
min) and cooling requirements (called the minimum cold Q(C)
min). Further optimisation of heat exchanger design to minimise economic cost are outside the scope of this research as process optimisation approaches are not considered in this study. For simplicity we have adopted the minimum allowable temperature difference (ΔT
min) at ‘pinch point’, in a heat exchanger to be 10 K [
23]. All calculations are performed using the Aspen Energy Analyser and results reported directly. Finally, the heat exchanger network is fully solved using the Aspen Energy Analyser
® with all the process streams satisfied based on the adopted minimum allowable temperature difference.
2.1.1. Calculation of Physical and Thermodynamic Properties
The NRTL-Redlich-Kwong property method in ASPEN plus was utilised in predicting the properties of the chemical species in the liquid and vapour phases. This is because the NRTL activity coefficient model was shown to be sufficient in predicting the vapour-liquid equilibria of chemical species under moderate conditions while the Redlich-Kwong model adequately predicts the vapour-liquid equilibria of the chemical species under supercritical and subcritical conditions [
24]. The triglyceride molecule in the reaction was modelled as triolein in this paper, in accordance with the experimental data of subcritical lipid hydrolysis of triolein and the kinetic models of the supercritical esterification of oleic acid and supercritical transesterification of triolein for producing methyl oleate [
6,
25,
26]. This methyl oleate product sufficiently models biodiesel since its major fuel properties satisfy existing European biodiesel standard (EN 14214) and American biodiesel standard (ASTM D 6751) as presented in
Table 1.
In this work, the values of the critical temperature, critical pressure and acentric factor for triolein calculated by Tang et al. [
33] were used to define triolein in the component database of the ASPEN plus. The boiling temperature of triolein was specified to be 548.3 °C [
34] rather than 846 °C available in ASPEN plus since previous studies have established that the boiling point temperature for triolein present in the ASPEN plus database is overestimated [
34].
The vapour pressure was calculated as follows [
35],
The thermodynamic properties of other chemicals such as glycerol, methanol, oleic acids, methyl oleate and water were obtained from the component database of the ASPEN plus. For simplicity perfect mixing of the reacting species has also been assumed.
2.1.2. Kinetic Modelling of Subcritical Lipid Hydrolysis and Supercritical Fatty Acid Esterification
In modelling the subcritical lipid hydrolysis reaction the FA yield data reported by Minami and Saka [
25] for subcritical triglyceride (TAG) hydrolysis of lipids undertaken at temperatures ranging from 250 °C to 320 °C and pressures ranging from 7 MPa to 20 MPa in a batch-wise reactor were analysed. In the experiments by Minami and Saka, excess water was utilised such that the molar ratio of water to TAG was 53.8:1 and the initial concentrations of the TAG and water were 0.56 mole/L and 27.75 mole/L, respectively. The experimental data obtained by Minami and Saka have been presented in
Figure 3.
The subcritical TAG hydrolysis was modelled as occurring according to the following reaction equation,
where, Gly represents the glycerol and FA represents the fatty acid.
This study recognises that TAG hydrolysis reaction have also been reported as occuring via a stage wise pathway under moderate reaction conditions with the intermediate molecules of diglyceride and mono glyceride shown to be formed during different stages of the process. Crucially however, under subcritical conditions such intermediate steps can be ignored with authors such as Kocsisová et al. [
36] and Sturzenegge and Stum [
37] demonstrating experimentally the validity of an alternative pseudo-homogenous First order irreversible kinetic relation provided the conditions of excess water and subcritical reaction conditions are maintained. According to these authors under conditions of excess moisture the rate of the TAG hydrolysis under subcritical conditions can be presented as follows,
Such that the integration of Equation (5) provides a relationship between the fractional conversion,
X, of the TAG, as follows,
or,
with the reaction rate constant (
khy) expressed by the following relation,
In Equation (8), khy, Ahy, Ehy, R and T represent the approximate overall rate constant, in s−1; pre-exponential constant, in s−1; activation energy, in kJ/kmol, universal gas constant, in 8.314 kJ/kmol·K and temperature in K, respectively.
Thus, a plot of Ln [1 − X] against time is a straight line and the slope of the line equal to value of −khy for a set of experimental data.
Using the data shown in
Figure 3, we have evaluated the conversion factor at different temperatures and the results are in
Figure 4. Applying Equation (7), fitted straight lines representing the relations between the conversion factor and reacting time at different temperatures are also presented in
Figure 4. From
Figure 4, it can be seen that all the straight lines have a correlation coefficient (
R2) value greater than 0.9. Since the
R2 values determined are all greater than the accepted minimum
R2 value of 0.7 specified in modelling work [
38], it can be suggested that the kinetic relationship represented by the linear plot is sufficient to predict the yield of FAs generated during subcritical hydrolysis and can be employed to achieve the aims of this paper. This implies therefore that the kinetic parameters (
Ahy and
Ehy), can be easily estimated by plotting Ln
khy, which is obtained from the graph slopes in
Figure 4 against 1/
T such that the slope and the intercept will give the
Ehy/
R and the
Ahy respectively. The plot of Ln
khy against 1/
T for temperatures ranging from 250 °C to 320 °C has therefore been presented in
Figure 5.
Figure 5 shows that the relationship between
khy and
T is as follows,
or
The kinetic relation presented in Equation (10) was therefore utilised in modelling the subcritical hydrolysis of TAGs in the present study.
The supercritical esterification of the FAs generated from the subcritical hydrolysis of the TAG in the ISHSE process is expressed by the following chemical reaction,
In this paper, the supercritical esterification of FAs was modelled as occurring via irreversible First order kinetics, following the experimental analysis by Jin et al. [
6] as follows,
Possible side reactions which are associated with the thermal cracking, isomerisation and decarboxylation of FAMEs will not be considered in this study since the mechanisms and reaction equations of these side reactions are presently unknown [
39].
The rate constant of the supercritical esterification reaction (
) at the operating temperature, molar ratio of methanol to FA and pressure of 260 °C, 20:1 and 8.1 MPa was reported, and is given by [
6],
where
khy is the rate constant of supercritical esterification in s
−1, the constant number 0.136 is the pre-exponential constant in s
−1, the constant value 21,980 is the activation energy in kJ/kmol,
R is the universal gas constant, 8.314 kJ/kmol·K and
T is the temperature in
K.
The supercritical transesterification of TAGs using methanol in the DST process has also been modelled in this paper as an irreversible a first order reaction, following the experimental analysis undertaken by He et al. [
26] as follows,
such that,
The rate constant of the supercritical transesterification reaction at the operating temperature, molar ratio of methanol to TAG and pressure of 280 °C, 42:1 and 28 MPa was reported as follows [
26],
where
ksp is the rate constant of supercritical transesterification, in s
−1, the constant value 141.796 is the pre-exponential constant, in s
−1, the constant value 56,000 is the activation energy, in kJ/kmole,
R is the universal gas constant, 8.314 kJ/kmole.
K and
T is the temperature, in
K.
2.1.3. Potential Environmental Impact Assessment
The potential environmental impacts (PEIs) of the alternative catalyst-free biodiesel production processes were investigated using a comprehensive environmental impact assessment (EIA) tool, the waste reduction (WAR) algorithm (WAR version 1.0.17, United States Environmental Protection Agency, Washington, DC, USA). The WAR impact assessment tool has been identified as the preferred impact assessment tool when it is necessary to compare anticipated environmental impacts associated with only the chemical production processing facility [
40] and has been used extensively in previous work [
41,
42]. Considering that the objectives of this study does not encompass the determination of suitable waste minimisation opportunities [
43], assessment of environmental risk (environmental impact minimization and the environmental fate and risk assessment tool) [
44] or the assessment of impacts associated with raw material acquisition, product distribution or product recycle (Ecoindicator-99 for life cycle analysis tool) [
45], the utilisation of these aforementioned alternative environmental assessment tools is pointless. The sufficiency of the utilisation of the WAR algorithm in achieving the research objectives is therefore justified.
WAR algorithm is a powerful environmental impact assessment tool that utilises the mass balance data describing the flow rate and composition of the streams entering into and exiting a chemical process and energy balance data describing the rate of energy consumption from fossil-based sources to calculate the potential impact of the a process if the chemicals utilised in the process are directly emitted to the environment [
45]. This environmental impact assessment tool was developed in response to the need for rapid assessments of environmental impacts during the conceptual and design stages of manufacturing processes [
46]. The WAR algorithm enables the determination of the conceptual PEI of a process via an approach analogous to the classic mass or energy balances, since the theoretical PEI of an a chemical process is assumed to be conserved. Therefore for steady state processes the PEI balance of a chemical process is as follows [
45],
Implying that the rate of PEI generation (
gen) is as follows,
where according to Young and Cabezas [
45],
In these Eqns.,
and
are the rate of PEIs into and out of a system due to chemical interactions within the system respectively;
and
are the rate of PEI out and into a system due to energy generation processes within the system respectively;
and
are the PEI out of a system as a result of the release of waste energy due to energy generation and chemical processes within a system respectively;
and
are the rate of waste energy emission from the chemical process and energy generation process. Typically, in the context of the WAR algorithm, the
and
terms are considered as negligible since chemical process plants do not produce significant ‘unwanted’ energy into the environment [
45]. Furthermore, for chemical plants the potential environmental impact associated with the emission of mass such as the spills of large mass of toxic chemicals or large masses of gaseous pollutants is usually much greater than that associated with the emission of energy such as rise in the temperature of the surrounding environment, Young and Cabezas [
45] also emphasize that the rate of PEI out of a system as a result of energy generation processes (
) are largely due to the generation of masses of gaseous pollutants such as SO
2, NO
x, CO and CO
2 which are a consequence of the combustion of solid, liquid and gaseous fuels. is the mass flow rate of stream
i output stream,
is the mass flow rate of stream
i input stream,
xki is the mass fraction of component
k in stream
i and
øk is the potential environmental impact due to component
k. The parameter,
øk, was determined by summing the specific potential environmental impact of component
k over all of the possible impact categories
l as follows [
45],
In Equation (25),
αi represents the relative weighting factor of impact categories,
l and is indicative of the importance of each impact category to the analysis. In this study, however all impact categories were considered equally significant. The approach of considering the impact categories equally significant constitutes the accepted approach when the case studies being assessed have no specific site in mind [
41,
42,
45]. The impact categories were classified largely based on extensive studies undertaken by Heijungs et al. [
47]. The potential environmental impacts of each chemical utilised in the chemical process were obtained from a chemical database similar to the database present in the ChemCad 4.0 software (ChemCad, Houston, TX, USA). The impact categories
l embeded in the WAR algorithm and the measure of impact category are presented in
Table 2. Normalisation of the impact scores of the potential environmental impact due to component
k, within each impact category were undertaken by the WAR algorithm to ensure that values of different pollution categories contain the same units to allow for their combination [
40]. The normalisation of the impact scores also ensures that values from different categories will have on the average equivalent scores thus eliminating unintentional bias in the calculation of the PEI indexes [
40]. The normalisation of the impact scores of the potential impact, due to component
k is achieved by the WAR algorithm using the following normalisation scheme [
40],
where (
Score)
kl represents the impact score or value of chemical
k on some arbitrary scale for impact category
l and [(
Score)
k]
l represents the average impact score or value of all chemicals in category
l.
In the present study only the energy inputs (heating utilities) required for the ISHSE and the DST biodiesel production processes were assumed to be provided by natural gas as a fossil-based, non-renewable energy source. The WAR software is freely available on the EPA website. Using information on participating chemical components, inlet, product and waste streams, and total energy consumption, the PEI per kg of products of both processes have been determined.