2. The mathematical model for biodi

2. The mathematical model for biodiesel production using
chemical catalyst
To describe a simple mathematical model for the alkali catalytic
transesterification reaction of Jatropha Curcas oil, the following
assumptions have been made.
(A1) Transesterification of vegetable oil with methanol consists of
three stepwise and reversible reactions [13,40].
(A2) Since there is only little water (0.2%, w/w) in the reaction
mixtures and negligible free fatty acid is detected in the system,
the hydrolysis reaction could be negligible [41].
(A3) Mixture of Jatropha Curcas waste food oil that contains 1 wt%
FFAs. TG, BD, FFA might take part in saponification with the
catalyst KOH [42].
(A4) Base catalyzed transesterification reaction for biodiesel production
form Jatropha Curcas oil with methanol (AL) involves
three reversible steps [40]. Here, in the first reversible step,
triglycerides (TG) reacts with alcohol (AL) to produce diglycerides
(DG) which further reacts with AL to produce monoglycerides
(MG) in the second reversible step. Finally, MG
reacts with AL to give glycerol (GL) as by product. At each
reaction step, one molecule of biodiesel (BD) is produced
for each molecule of AL consumed. The kinetic steps are
shown by the following sequential reaction
TG þ AL
k1
k2
DG þ BD;
DG þ AL
k3
k4
MG þ BD;
MG þ AL
k5
k6
GL þ BD:
ð1Þ
It is shown above that the alcoholysis kinetic reaction scheme
consists of three reversible reactions. Beside, there are also
TG, FFA and BD saponification reactions. It is shown below
that the saponification kinetic reaction scheme consists of
mainly three reactions [4].
TG þ K !
k7
S þ GL
BD þ K !
k8
S þ AL
F þ K !
k9
S þ W:
ð2Þ
504 F.A. Basir et al. / Fuel 158 (2015) 503–511
(A5) Here k1; k3; k5 are forward reaction rate constants and k2; k4
and k6 are backward reaction rate constants. The values of
k1 to k6 are shown in Table 1. We denote the concentration
of TG, DG, MG, BD, AL, FFA, K, S, W and GL by
xT ; xD; xM; xB; xA; xF ; xK ; xS; xW and xG respectively. The dependency
of reaction rate constants on the temperature is
expressed by the Arrhenius equation as below [43]
ki ¼ aie
bi
T : ð3Þ
T is the reaction temperature, ai is the frequency factor, and
bi ¼ Eai
R ;
in which Eai is the activation energy for each component and
R is the universal gas constant. Using the values of ai and bi
[43], we obtain the values of ki.
(A6) As mixing intensity in the reaction system directs the mass
transfer limitations between phases, thus mechanical stirring
is one of the most effective factors in transesterification
reaction. Here, we use ks as the mass transfer rate constant
and its unit is min1 and the term has been defined as below
[33]:
ks ¼ a
1 þ expðbðN cÞÞ; ð4Þ
where N is the speed of stirrer and a, b and c are constants.
The term is used in our model by the expression
ksxBð1 xB
Bmax Þ. Here xB denotes the concentration of biodiesel
and Bmax represents maximum production of biodiesel. Here
xB denotes the concentration of biodiesel and Bmax represents
maximum biodiesel production in an ideal situation which is
defined as system having no mass transfer resistance. The
term has been used logistically because with the increase
of stirrer speed, the mass transfer resistance decreases and
beyond a certain stirrer speed the mass transfer resistance
is negligible. This is also evident from experimental observations
of other workers [44].
Based on the above assumptions, the differential equations
characterizing the stepwise base catalyzed reactions may be enunciated
follows:
dxB
dt ¼ k1xT xA k2xDxB þ k3xDxA k4xMxB þ k5xMxA k6xGxB
k8xBxK þ kSxB 1 xB
Bmax
dxT
dt ¼ k1xT xA þ k2xDxB k7xT xK ;
dxD
dt ¼ k1xT xA k2xDxB k3xDxA þ k4xMxB;
dxM
dt ¼ k3xDxA þ k4xMxB k5xMxA þ k6xGxB;
dxA
dt ¼ k1xT xA þ k2xDxB k3xDxA þ k4xMxB k5xMxA þ k6xGxB;
dxG
dt ¼ k5xMxA þ k6xGxB þ k7xT xK :
dxK
dt ¼ dxS
dt ¼ ðk7xT xK þ k8xBxK þ k9xK xF Þ
dxF
dt ¼ dxW
dt ¼ k9xF xK ; ð5Þ
with initial conditions,
xBð0Þ ¼ 0; xT ð0Þ ¼ xT0 ; xDð0Þ ¼ 0; xMð0Þ ¼ 0; xAð0Þ ¼ xA0 ; xF ð0Þ ¼ 0;
xK ¼ xK0 , xS ¼ 0; xW ¼ 0 and xGð0Þ ¼ 0.
3. The mathematical model for biodiesel production using
enzymatic reaction
To describe a simple mathematical model for enzyme catalytic
transesterification reaction of Jatropha Curcas oil, the following
assumptions have been taken
(A) Biocatlytic catalyzed transesterification (using Lipase) of
Jatropha Curcas oil with an alcohol (AL) may be proposed
to involve two-step mechanisms. The first step consists of
hydrolysis of Jatropha Curcas oil or TG to produce acylated
enzyme (AcE) and release of glycerol as only by product
through a complex C1 ðE:TGÞ. Here E stands for enzyme.
The second step is the esterification of methanol (AL)
with AcE to form the desired product i.e. BD with the
release of free enzyme (E) through a second complex C2
ðAcE:ALÞ [45].
(A2) To account for the inhibition of reaction by methanol in the
system kinetics, here, a competitive inhibition is assumed
when a AL molecule combines with the E to produce a
dead-end complex, complex C3 (E.AL) [18,22,26,46]. All the
mechanistic steps for the biodiesel production can be represented
by the following sequence of reactions:
E þ TG
r1
r1
½E:TG
r2
r2
AcE þ GL;
AcE þ AL
r3
r3
½AcE:AL
r4
r4
E þ BD;
E þ AL
r5
r5
½E:AL: ð6Þ
Here, r1 and r1;r2 and r2 are the rate constants for the
reversible formation of complex C1, acylated enzyme and
by product glycerol respectively in the first step of biodiesel
formation. r3 and r3;r4 and r4 are the rate constants for
the reversible formation of complex C2 and biodiesel formation
respectively in the final step. For the inhibition reaction,
r5 and r5 are the rate constants for the reversible formation
of dead-end complex C3.
(A3) As without stirring, the rate of transesterification reaction is
very slow and incomplete due to mass transfer resistance
[26,47], so the effect of stirring is included in the reaction
kinetics in a logistic fashion as earlier model by the two
expressions, ksxB 1 xB
Bmax   and ksxC3 1 xC3
C3max  , where C3max
is the maximum concentration of dead-end complex.
We denote the concentration of TG; E; AcE; C1; C2; C3; AL; BD and
GL as xT ; xE; xAcE; xC1; xC2; xC3; xA; xB and xG respectively. Now from the
above assumptions with the above reaction mechanism followed
by law of mass action we can formulate the set of differential equations
given below:
dxE
dt ¼ r1xT xE þ r1xC1 þ r4xC2 r4xExB r5xExA þ r5xC3;
dxT
dt ¼ r1xT xE þ r1xC1 r5xAcE þ r5xC3;
dxAcE
dt ¼ r2xC1 r2xAcExG r3xF xA þ r3xC2 r5xAcE þ r5xC3;
dxB
dt ¼ r4xC2 r4xExB þ ksxB 1 xB
Bmax ;
F.A. Basir et al. / Fuel 158 (2015) 503–511 505
dxA
dt ¼ r3xAcExA þ r3xC2 r5xExA þ r5xC3;
dxC1
dt ¼ r1xT xE r1xC1 r2xC1 þ r2xAcExG;
dxC2
dt ¼ r3xAcExA r3xC2 r4xC2 þ r4xExB;
dxC3
dt ¼ r5xExA r5xC3 þ ksxC3 1 xC3
C3max ;
dxG
dt ¼ r2xC1 r2xAcExG; ð7Þ
with the initial conditions
xEð0Þ ¼ xE0 ; xC1ð0Þ ¼ 0; xAcEð0Þ ¼ 0;
xAð0Þ ¼ xA0 ; xT ð0Þ ¼ xT0 ; xC2ð0Þ ¼ 0;
xBð0Þ ¼ 0; xGð0Þ ¼ 0 and xC3ð0Þ ¼ 0: ð8Þ
The last term, in expression dxB
dt and dxC3
dt , represents only the mass
transfer rate or indirectly the contribution of the mass transfer
resistance to the overall resistance. Here C3max is the maximum concentration
of the dead-end complex xC3.
4. The optimal control problem for enzymatic reaction system
Our main objective is to control the rate of stirrer (i.e. mass
transfer rate), so that we can get maximum production of biodiesel.
Also, it is our object to keep the cost function as low as possible.
We use a control variable uðtÞ, which represents the stirring activator
input at time t satisfying 0 6 uðtÞ 6 1. Here uðtÞ represents control
input with values normalized to be between 0 and 1 [48]. Also
uðtÞ ¼ 1 represents the maximal use of stirrer and uðtÞ  0, which
signifies no change in stirrer rotation.
Based on the above assumptions, the control induced system
corresponding to the system (7) would be:
dxE
dt ¼ r1xT xE þ r1xC1 þ r4xC2 r4xExB;
dxT
dt ¼ r1xT xE þ r1xC1 r5xAcE þ r5xC3;
dxAcE
dt ¼ r2xC1 r2xAcExG r3xF xA þ r3xC2 r5xAcE þ r5xC3;
dxB
dt ¼ k4xC2 r4xExB þ uðtÞksxB 1 xB
Bmax ;
dxA
dt ¼ r3xAcExA þ r3xC2;
dxC1
dt ¼ r1xT xE r1xC1 r2xC1 þ r2xAcExG;
dxC2
dt ¼ r3xAcExA r3xC2 r4xC2 þ r4xExB;
dxC3
dt ¼ r5xAcE r5xC3 þ uðtÞksxC3 1 xC3
C3max ;
dxG
dt ¼ r2xC1 r2xAcExG; ð9Þ
with initial conditions as given by Eq. (8).
The cost function is thus formulated as
J½uðtÞ ¼ Z tf
t0
Pu2
ðtÞ Qx2
BðtÞ þ Rx2
C3
h idt ð10Þ
Here, we want to maximize biodiesel production and reduce the
inhibition of enzyme i.e. minimize the concentration of dead-end
complex xC3. The parameter Pð> 0Þ is the weight constant on the
benefit of the cost of production and Q and R are the penalty
multipliers.
Now, we want to find out the optimal control u such that
Jðu
Þ ¼ min JðuÞ : u 2 U
where U is the admissible control set defined by
U ¼fuðtÞ : uðtÞ is measurable; 0 6 uðtÞ 6 1;t 2 ½ti;tf g
ð11Þ
Here we use Pontryagin Minimum Principle [49,50] to find uðtÞ.
The Hamiltonian is formulated as,
H ¼ Pu2
ðtÞ Qx2
BðtÞ þ Rx2
C3 þX9
i¼1
nifi; ð12Þ
where n1 n9 are adjoint variables and fi, i = 1, 2, ..., 9 are the right
side of Eq. (9), i.e. f 1 ¼ r1xT xE þ r1xC1 þ r4xC2 r4xExB, etc.
Theorem. If the given optimal control uðtÞ and the solution
x
E; x
T ; x
AcE; x
B; x
A; x
C1; x
C2; x
C3
  of the corresponding system (Eq. 7)
minimize JðuÞ over U, then there exists adjoint variables n1 n9 which
satisfying the following equations,
dn1
dt ¼ r1xT ðn1 n6Þ þ r4xBðn4 n7Þ þ r5xAðn5 n8Þ;
dn2
dt ¼ r1xEðn1 þ n2 n6Þ;
dn3
dt ¼ r3xAðn5 n7Þ þ r5xT ðn8Þ þ r2xGðn3 þ n9 n6Þ;
dn4
dt ¼ 2QxB þ n4 r4xE þ
ksxBuðtÞ
Bmax
ksuðtÞ 1 xB
Bmax  n7r4xE;
dn5
dt ¼ r3xAcEðn5 n7Þ þ r5xEðn1 n8Þ;
dn6
dt ¼ r1ðn6 n2Þ þ r2ðn6 n3 n9Þ;
dn7
dt ¼ r4n4 r3 þ r3n7;
dn8
dt ¼ r5ðn8 n2 n3Þ þ n8uðtÞks 1 xC3
C3max
þ n8uðtÞksxC3 1 1
C3max ;
dn9
dt ¼ r2xAcEðn3 n6 þ n9Þ: ð13Þ
along with the transversality condition niðtfÞ ¼ 0 for i ¼ 1 to 9.
Further, uðtÞ can be written as,
u
ðtÞ ¼ ma

ผลลัพธ์ (ไทย) 1:

2. แบบจำลองทางคณิตศาสตร์สำหรับการใช้ผลิตไบโอดีเซลสารเคมีเศษเพื่ออธิบายแบบจำลองทางคณิตศาสตร์อย่างง่ายสำหรับด่างตัวเร่งปฏิกิริยาเพิ่มปฏิกิริยาของน้ำมันสบู่ดำ Curcas ต่อไปนี้ได้ทำการสมมติฐาน(A1) ของน้ำมันพืชกับเมทานอลประกอบด้วยสาม stepwise และผันกลับได้ปฏิกิริยา [13,40](A2) เนื่องจากมีเพียงเล็กน้อยน้ำ (0.2%, w/w) ในปฏิกิริยาส่วนผสมและกรดไขมันอิสระระยะตรวจพบในระบบไฮโตรไลซ์ปฏิกิริยาอาจเป็นระยะ [41](A3) ส่วนผสมของสบู่ดำ Curcas เสียน้ำมันอาหารที่ประกอบด้วย 1 wt %FFAs TG, BD, FFA อาจมีส่วนร่วมในการสะพอนิฟิกับการเศษเกาะ [42](A4) พื้นฐานของปฏิกิริยาในกระบวนการผลิตไบโอดีเซลเกี่ยวข้องกับแบบฟอร์มน้ำมันสบู่ดำ Curcas ด้วยเมทานอล (AL)สามย้อนกลับขั้นตอน [40] ที่นี่ ในขั้นตอนแรกใหม่ทั้งหมดระดับไตรกลีเซอไรด์ (TG) ทำปฏิกิริยากับแอลกอฮอล์ (AL) ผลิต diglycerides(กิจ) ที่ไปทำปฏิกิริยากับอัลผลิต monoglycerides(MG) ในขั้นตอนสองที่สามารถย้อนกลับ สุดท้าย มิลลิกรัมทำปฏิกิริยากับอัลให้กลีเซอร (GL) โดยผลิตภัณฑ์ ในแต่ละผลิตขั้นตอนปฏิกิริยา โมเลกุลหนึ่งของไบโอดีเซล (BD)สำหรับแต่ละโมเลกุลของอัลการบริโภค ขั้นตอนเดิม ๆแสดงปฏิกิริยาต่อไปนี้ตามลำดับÞ TG อัลk1k2Þกิจ BDÞกิจอัลk3k4Þ MG BDMG þ ALk5k6GL þ BD:ð1Þมันจะแสดงอยู่เหนือโครงร่างเคลื่อนไหวปฏิกิริยา alcoholysisประกอบด้วยสามปฏิกิริยาผันกลับได้ ข้าง มีปฏิกิริยาสะพอนิฟิ TG, FFA และ BD ที่แสดงด้านล่างว่า โครงร่างเคลื่อนไหวปฏิกิริยาสะพอนิฟิประกอบด้วยสามส่วนใหญ่ปฏิกิริยา [4]TG þ Kk7S þ GLBD þ Kk8Þ S อัลF þ Kk9Þ S ไร w:ð2ÞAl. et Basir หลังการเก็บเกี่ยวที่ 504 / เชื้อเพลิง 158 (2015) 503-511(A5) K1 ที่นี่ k3 k5 เป็นค่าคงที่อัตราปฏิกิริยาไปข้างหน้าและ k2 k4และ k6 มีค่าคงที่อัตราปฏิกิริยาย้อนหลัง ค่าของk1 กับ k6 จะแสดงในตารางที่ 1 เราแสดงความเข้มข้นของ TG กิจ MG, BD, AL, FFA, K, S, W และ GL โดยxT xD xM xB ซา คุณ xK xS xW และ xG ตามลำดับ การขึ้นต่อกันค่าคงที่อัตราปฏิกิริยากับอุณหภูมิคือแสดง โดยสมการของอาร์เรเนียสเป็นด้านล่าง [43]กี่¼ aiebiT: ð3ÞT คือ อุณหภูมิปฏิกิริยา ไอมีตัวคูณความถี่ และbi ¼ EaiRซึ่ง Eai ถูก พลังงานกระตุ้นสำหรับแต่ละส่วน และR คือ ค่าคงสากลของก๊าซ ใช้ค่า ai และ bi[43], เรารับค่า ki(A6) เป็นการผสมความเข้มในปฏิกิริยา ระบบชักจูงมวลชนข้อจำกัดการโอนย้ายระหว่างระยะ เครื่องจักรกลจึงกวนเป็นหนึ่งในปัจจัยที่มีประสิทธิภาพสูงสุดในการเพิ่มปฏิกิริยาการ ที่นี่ เราใช้ ks เป็นค่าคงอัตราการถ่ายโอนมวลและหน่วยของ min1 และระยะเวลาที่กำหนดไว้ด้านล่าง[33]:ks ¼เป็น1 þ expðbðN cÞÞ ð4Þโดยที่ N คือ ความเร็วของช้อนคน และ a, b และ c เป็นค่าคงที่คำที่ใช้ในรุ่นของเรา โดยนิพจน์ksxBð1 xBBmax Þ นี่ xB แสดงความเข้มข้นของไบโอดีเซลและ Bmax แทนผลิตไบโอดีเซลสูงสุด ที่นี่xB แสดงความเข้มข้นของไบโอดีเซล และ Bmax แทนผลิตไบโอดีเซลสูงสุดในสถานการณ์ที่เหมาะซึ่งกำหนดเป็นระบบที่มีความต้านทานการถ่ายโอนมวลไม่ ที่ระยะใช้วันนี้เนื่องจาก มีการเพิ่มขึ้นช้อนคนความเร็ว มวลโอนลดลงความต้านทาน และนอกเหนือจากความบางช้อนคนเร็วต้านทานการถ่ายโอนมวลเป็นระยะ นี้ก็เห็นได้จากการสังเกตทดลองของคนที่อื่น ๆ [44]ตามสมมติฐานข้างต้น สมการเชิงอนุพันธ์กำหนดลักษณะของปฏิกิริยากระบวนพื้นฐาน stepwise อาจจะ enunciatedดังนี้:dxBdt ¼ k1xT ซา k2xDxB þ k3xDxA k4xMxB þ k5xMxA k6xGxB k8xBxK þ kSxB 1 xBBmax dxTdt ¼ k1xT xA þ k2xDxB k7xT xKdxDdt ¼ k1xT ซา k2xDxB k3xDxA þ k4xMxBdxMdt ¼ k3xDxA þþ k4xMxB k5xMxA k6xGxBdxAdt ¼ k1xT xA þ k2xDxB k3xDxA þ k4xMxB k5xMxA þ k6xGxBdxGdt ¼ k5xMxA þ k6xGxB þ k7xT xK:dxKdxS dt ¼dt ¼ ðk7xT xK þ k8xBxK þ k9xK xF ÞdxFdxW dt ¼dt ¼ k9xF xK ð5Þกับเงื่อนไขเริ่มต้นxBð0Þ ¼ 0 xT ð0Þ ¼ xT0 xDð0Þ ¼ 0 xMð0Þ ¼ 0 xAð0Þ ¼ xA0 คุณ ð0Þ ¼ 0xK ¼ xK0, xS ¼ 0 xW ¼ 0 และ xGð0Þ ¼ 03.แบบจำลองทางคณิตศาสตร์สำหรับการใช้ผลิตไบโอดีเซลปฏิกิริยาที่เอนไซม์ในระบบเพื่ออธิบายแบบจำลองทางคณิตศาสตร์อย่างง่ายสำหรับตัวเร่งปฏิกิริยาของเอนไซม์เพิ่มปฏิกิริยาของน้ำมันสบู่ดำ Curcas ต่อไปนี้มีการใช้สมมติฐาน(A) Biocatlytic กระบวนเพิ่ม (ใช้เอนไซม์ไลเปส) ของอาจเสนอน้ำมัน Curcas สบู่ดำกับเครื่องดื่มแอลกอฮอล์ (AL)การเกี่ยวข้องกับกลไกสองขั้นตอน ขั้นตอนแรกประกอบด้วยไฮโตรไลซ์น้ำมันสบู่ดำ Curcas หรือ TG ผลิต acylatedเอนไซม์ (เอซ) และรุ่นของกลีเซอรเพียงโดยผลิตภัณฑ์ผ่าน C1 ที่ซับซ้อน ðE:TGÞ นี่อียืนสำหรับเอนไซม์ขั้นตอนสองเป็น esterification ของเมทานอล (AL)มี AcE แบบสินค้าที่ระบุเช่น BD มีการของฟรีเอนไซม์ (E) ผ่าน C2 ซับซ้อนเป็นที่สองðAcE:ALÞ [45](A2) การยับยั้งปฏิกิริยา ด้วยเมทานอลในการจลนพลศาสตร์ของระบบ ที่นี่ ยับยั้งแข่งขันสันนิษฐานเมื่อโมเลกุลอัลเด่น E เพื่อผลิตเป็นdead-end ซับซ้อน ซับซ้อน C3 (E.AL) [18,22,26,46] ทั้งหมดนี้สามารถแสดงกลไกการทำขั้นตอนการผลิตไบโอดีเซลโดยลำดับของปฏิกิริยาต่อไปนี้:E þ TGr1r1½E:TGr2r2Þเอ GLÞเออัลr3r3½AcE:ALr4r4Þ E BDE þ ALr5r5½E:AL: ð6Þที่นี่ r1 และ r1, r2 และ r2 มีค่าคงที่อัตราในการกลับก่อตัวของซับซ้อน C1 เอนไซม์ acylated และโดยกลีเซอรผลิตภัณฑ์ในขั้นตอนแรกของไบโอดีเซลตามลำดับการก่อ r3 และ r4; r3 และ r4 มีค่าคงที่อัตราการการก่อตัวกลับของผู้แต่ง C2 และไบโอดีเซลที่ซับซ้อนในขั้นตอนสุดท้ายตามลำดับ สำหรับปฏิกิริยายับยั้งr5 และ r5 จะคงอัตราการก่อตัวกลับของ dead-end C3 ซับซ้อน(A3) โดยไม่ต้องกวน อัตราของปฏิกิริยาเพิ่มเป็นช้า และไม่สมบูรณ์เนื่องจากความต้านทานการถ่ายโอนมวล[26,47], เพื่อผลของการกวนอยู่ในปฏิกิริยาจลนพลศาสตร์ในโลจิสติกเป็นรุ่นก่อนหน้านี้โดยทั้งสองนิพจน์ xB ksxB 1Bmax ksxC3 และ 1 xC3C3max ที่ C3maxมีความเข้มข้นสูงสุดของคอมเพล็กซ์ dead-endเราแสดงความเข้มข้นของ TG E เอซ C1 C2 C3 อัล BD และGL เป็น xT xE xAcE xC1 xC2 xC3 ซา xB และ xG ตามลำดับ ตอนนี้จากตามข้างบนสมมติฐานด้วยกลไกปฏิกิริยาข้างต้นตามกฎหมายของการกระทำโดยรวม เราสามารถกำหนดชุดของสมการเชิงอนุพันธ์กำหนดให้ด้านล่าง:dxEdt ¼ r1xT xE þ r1xC1 þ r4xC2 r4xExB r5xExA þ r5xC3dxTdt ¼ r1xT xE þ r1xC1 r5xAcE þ r5xC3dxAcEdt ¼ r2xC1 r2xAcExG r3xF ซาþ r3xC2 r5xAcE þ r5xC3dxBdt ¼ r4xC2 r4xExB þ ksxB 1 xBBmaxหลังการเก็บเกี่ยว Basir et al. / เชื้อเพลิง 158 (2015) 503-511 505dxAdt ¼ r3xAcExA þþ r3xC2 r5xExA r5xC3dxC1dt ¼ r1xT xE r1xC1 r2xC1 þ r2xAcExGdxC2dt ¼ r3xAcExA r3xC2 r4xC2 r4xExB þdxC3ksxC3 dt ¼ r5xExA r5xC3 þ 1 xC3C3maxdxGdt ¼ r2xC1 r2xAcExG ð7Þกับเงื่อนไขเริ่มต้นxEð0Þ ¼ xE0 xC1ð0Þ ¼ 0 xAcEð0Þ ¼ 0xAð0Þ ¼ xA0 xT ð0Þ ¼ xT0 xC2ð0Þ ¼ 0xBð0Þ ¼ 0 xGð0Þ ¼ 0 และ xC3ð0Þ ¼ 0: ð8Þระยะสุดท้าย ในนิพจน์ dxBdt และ dxC3dt แสดงเฉพาะมวลอัตราการถ่ายโอนหรือโดยทางอ้อมของการถ่ายโอนมวลความต้านทานต่อความต้านทานรวม C3max เป็นความเข้มข้นสูงสุดของ xC3 ซับซ้อน dead-end4.ปัญหาระบบปฏิกิริยาที่เอนไซม์ในระบบควบคุมที่เหมาะสมวัตถุประสงค์หลักของเราคือการ ควบคุมอัตราการช้อนคน (เช่นมวลอัตราการถ่ายโอน), เพื่อให้เราได้สูงสุดผลิตไบโอดีเซลยัง เป็นวัตถุของเราเพื่อให้ฟังก์ชันต้นทุนต่ำที่สุดเราใช้การควบคุมตัวแปร uðtÞ ซึ่งแสดงถึง stirring activatorป้อนข้อมูลที่เวลา t ความพึงพอใจ 0 6 6 1 uðtÞ ที่นี่ uðtÞ แทนการควบคุมป้อนค่าที่ตามปกติจะอยู่ระหว่าง 0 และ 1 [48] นอกจากนี้uðtÞ ¼ 1 แทนใช้ช้อนคนและ uðtÞ 0 สูงสุดที่หมายถึงการเปลี่ยนแปลงหมุนช้อนคนตามสมมติฐานข้างต้น การควบคุมทำให้เกิดระบบที่สอดคล้องกับระบบ (7) จะเป็น:dxEdt ¼ r1xT xE þ r1xC1 þ r4xC2 r4xExBdxTdt ¼ r1xT xE þ r1xC1 r5xAcE þ r5xC3dxAcEdt ¼ r2xC1 r2xAcExG r3xF ซาþ r3xC2 r5xAcE þ r5xC3dxBdt ¼ k4xC2 r4xExB þ uðtÞksxB 1 xBBmaxdxAdt ¼ r3xAcExA þ r3xC2dxC1dt ¼ r1xT xE r1xC1 r2xC1 þ r2xAcExGdxC2dt ¼ r3xAcExA r3xC2 r4xC2 r4xExB þdxC3uðtÞksxC3 dt ¼ r5xAcE r5xC3 þ 1 xC3C3maxdxGdt ¼ r2xC1 r2xAcExG ð9Þมีเงื่อนไขเบื้องต้นที่กำหนดโดย Eq. (8)ฟังก์ชันต้นทุนจึงมีสูตรเป็นJ½uðtÞ ¼ Z tft0Pu2ðtÞ Qx2BðtÞ þ Rx2C3h idt ð10Þที่นี่ เราต้องการขยายการผลิตไบโอดีเซล และลดการยับยั้งเอนไซม์เช่นลดความเข้มข้นของ dead-endxC3 ซับซ้อน พารามิเตอร์ Pð > 0Þ เป็นค่าคงของน้ำหนักประโยชน์ของต้นทุนผลิต Q และ R มีเมลล์multipliersตอนนี้ เราต้องการค้นหายูควบคุมเหมาะสมที่JðuÞนาที¼ JðuÞ: U u 2โดยที่ U คือ ชุดควบคุม admissible กำหนดโดยU ¼fuðtÞ: uðtÞ เป็นวัด 0 6 uðtÞ 6 1; t 2 ½ti; tf gð11Þที่นี่เราใช้หลักน้อย Pontryagin [49,50] หา u ðtÞHamiltonian ที่เป็นสูตรเป็นH ¼ Pu2ðtÞ Qx2BðtÞ þ Rx2C3 þX9i¼1nifi ð12Þn1 n9 จะแปร adjoint และเน็ตไร้สาย ฉัน = 1, 2,..., 9 อยู่ด้านขวาด้านของ Eq. (9), เช่น f 1 ¼ r1xT xE þ r1xC1 þ r4xC2 r4xExB ฯลฯทฤษฎีบท ถ้า ðtÞ u ควบคุมให้เหมาะสมและการแก้ไขxE xT xเอซ xB xA xC1 xC2 xC3ระบบที่เกี่ยวข้อง (Eq. 7)ลด JðuÞ ผ่าน U แล้วมี adjoint แปร n1 n9 ซึ่งความพึงพอใจของสมการต่อไปนี้dn1dt ¼ r1xT ðn1 n6Þ þ r4xBðn4 n7Þ þ r5xAðn5 n8Þdn2dt ¼ r1xEðn1 þ n2 n6Þdn3dt ¼ r3xAðn5 n7Þ þ r5xT ðn8Þ þ r2xGðn3 þ n9 n6Þdn4dt ¼ 2QxB þ n4 r4xE þksxBuðtÞBmax xB ksuðtÞ 1Bmax n7r4xEdn5dt ¼ r3xAcEðn5 n7Þ þ r5xEðn1 n8Þdn6dt ¼ r1ðn6 n2Þ þ r2ðn6 n3 n9Þdn7dt ¼ r4n4 r3 þ r3n7dn8dt ¼ r5ðn8 n2 n3Þ þ n8uðtÞks 1 xC3C3max þ n8uðtÞksxC3 1 1C3maxdn9dt ¼ r2xAcEðn3 n6 þ n9Þ: ð13Þพร้อมกับ transversality เงื่อนไข niðtfÞ ¼ 0 หา¼ 1-9สามารถเขียน u ðtÞ เป็น เพิ่มเติมuma ðtÞ ¼

ผลลัพธ์ (ไทย) 2:

2. แบบจำลองทางคณิตศาสตร์สำหรับการผลิตไบโอดีเซลโดยใช้ตัวเร่งปฏิกิริยาเคมีเพื่ออธิบายแบบจำลองทางคณิตศาสตร์ที่เรียบง่ายสำหรับตัวเร่งปฏิกิริยาด่างปฏิกิริยาtransesterification ของน้ำมันสบู่ดำตามสมมติฐานที่ได้รับการทำ. (A1) Transesterification ของน้ำมันพืชที่มีเมทานอลประกอบด้วยสามขั้นตอนและย้อนกลับปฏิกิริยา [13,40]. (A2) เนื่องจากมีน้ำเพียงเล็กน้อย (0.2% w / w) ในการทำปฏิกิริยาผสมและกรดไขมันที่สำคัญฟรีมีการตรวจพบในระบบปฏิกิริยาการย่อยสลายที่อาจจะเล็กน้อย[41]. ( A3) ส่วนผสมของน้ำมันอาหารขยะสบู่ดำที่มีน้ำหนัก 1% FFAs TG, BD, FFA อาจจะมีส่วนร่วมในสะพอกับตัวเร่งปฏิกิริยาเกาะ[42]. (A4) ฐานเร่งปฏิกิริยา transesterification สำหรับผลิตไบโอดีเซลน้ำมันสบู่ดำรูปแบบกับเมทานอล(AL) ที่เกี่ยวข้องกับสามขั้นตอนย้อนกลับได้[40] ที่นี่ในขั้นตอนแรกย้อนกลับ, ไตรกลีเซอไรด์ (TG) ทำปฏิกิริยากับเครื่องดื่มแอลกอฮอล์ (AL) ในการผลิต Diglycerides (DG) ซึ่งต่อไปจะทำปฏิกิริยากับอัลในการผลิต monoglycerides (MG) ​​ในขั้นตอนที่ย้อนกลับที่สอง สุดท้าย MG ทำปฏิกิริยากับอัลจะให้กลีเซอรอล (GL) โดยผลิตภัณฑ์ ในแต่ละขั้นตอนของปฏิกิริยาหนึ่งโมเลกุลของน้ำมันไบโอดีเซล (BD) คือการผลิตโมเลกุลของอัลแต่ละบริโภค ขั้นตอนการเคลื่อนไหวจะแสดงโดยปฏิกิริยาตามลำดับต่อไปนี้TG þ AL k1 k2 DG þ BD; DG þ AL k3 k4 MG þ BD; MG þ AL K5 K6 GL þ BD: ð1Þมันแสดงให้เห็นข้างต้นว่าโครงการปฏิกิริยาการเคลื่อนไหว alcoholysis ประกอบด้วย สามปฏิกิริยาย้อนกลับ นอกจากนี้ยังมีTG, FFA และ BD ปฏิกิริยาสะพอ มันแสดงให้เห็นด้านล่างว่าโครงการปฏิกิริยาสะพอพลังงานจลน์ประกอบด้วยสามส่วนใหญ่เกิดปฏิกิริยา[4]. TG þ K! k7 S þ GL BD þ K! K8 S þอัลF þ K! k9 S þ W: ð2Þ 504 เอฟเอ Basir et al, . เชื้อเพลิง / 158 (2015) 503-511 (A5) นี่ k1; k3; K5 ไปข้างหน้ามีค่าคงที่อัตราการเกิดปฏิกิริยาและ k2; k4 และ K6 มีค่าคงที่อัตราการเกิดปฏิกิริยาย้อนกลับ ค่าของk1 K6 จะแสดงในตารางที่ 1 เราแสดงถึงความเข้มข้นของTG, DG, MG, BD, AL, FFA, K, S, W และ GL โดยXT; xD; xM; xB; xA; XF; XK; xS; และ XW xG ตามลำดับ พึ่งพาของค่าคงที่อัตราการเกิดปฏิกิริยากับอุณหภูมิจะแสดงโดยสมArrhenius ดังต่อไปนี้ [43] ki ¼ AIE สองT: ð3Þ T คืออุณหภูมิปฏิกิริยา ai เป็นปัจจัยความถี่และสอง¼ Eai R; ที่ Eai เป็น พลังงานกระตุ้นสำหรับแต่ละองค์ประกอบและR คือค่าคงที่ของก๊าซสากล ใช้ค่าไอและสอง[43] เราได้รับค่า ki. the (A6) ในฐานะที่ผสมความเข้มในระบบของปฏิกิริยาชี้นำมวลข้อจำกัด การโอนระหว่างขั้นตอนจึงกวนกลเป็นหนึ่งในปัจจัยที่มีประสิทธิภาพมากที่สุดในtransesterification ปฏิกิริยา ที่นี่เราใช้แคนซัสเป็นการถ่ายโอนมวลอัตราดอกเบี้ยคงที่และหน่วยของมันคือ min1 และระยะได้ถูกกำหนดให้ดังต่อไปนี้ [33]: KS ¼ 1 þexpðbðNcÞÞ; ð4Þที่ N คือความเร็วในการกวนและ A, B และ C เป็นค่าคงที่. เป็นคำที่นำมาใช้ในรูปแบบของเราโดยการแสดงออกksxBð1 xB Bmax Þ นี่ xB หมายถึงความเข้มข้นของไบโอดีเซลและBmax แสดงให้เห็นถึงการผลิตสูงสุดของไบโอดีเซล นี่xB หมายถึงความเข้มข้นของไบโอดีเซลและ Bmax แสดงให้เห็นถึงการผลิตไบโอดีเซลสูงสุดสถานการณ์ที่เหมาะซึ่งถูกกำหนดให้เป็นระบบที่ไม่มีความต้านทานการถ่ายเทมวล ระยะได้รับการใช้ logistically เพราะมีการเพิ่มขึ้นของความเร็วในการกวน, ความต้านทานการถ่ายเทมวลลดลงและเกินความเร็วกวนบางต้านทานการถ่ายเทมวลเป็นเล็กน้อย นอกจากนี้ยังเห็นได้ชัดจากการสังเกตการทดลองของคนงานอื่น ๆ [44]. จากสมมติฐานข้างต้นสมการเชิงอนุพันธ์พัฒนาการฐานแบบขั้นตอนเร่งปฏิกิริยาอาจจะกำหนดดังนี้DXB dt ¼ k1xT xA k2xDxB þ k3xDxA k4xMxB þ k5xMxA k6xGxB k8xBxK þ kSxB 1 xB Bmax DXT dt ¼ k1xT xA þ k2xDxB k7xT XK; DXD dt ¼ k1xT xA k2xDxB k3xDxA þ k4xMxB; DXM dt ¼ k3xDxA þ k4xMxB k5xMxA þ k6xGxB; DXA dt ¼ k1xT xA þ k2xDxB k3xDxA þ k4xMxB k5xMxA þ k6xGxB; DXG dt ¼ k5xMxA þ k6xGxB þ k7xT XK: dxK dt ¼ dxS dt ¼ðk7xT XK þ k8xBxK þ k9xK XF Þ DXF dt ¼ dxW dt ¼ k9xF XK; ð5Þมีเงื่อนไขเริ่มต้นxBð0Þ¼ 0; XT ð0Þ¼ xT0; xDð0Þ¼ 0; xMð0Þ¼ 0; xAð0Þ¼ xA0; XF ð0Þ¼ 0; XK ¼ xK0, xS ¼ 0; XW ¼ 0 xGð0Þ¼ 0. 3 แบบจำลองทางคณิตศาสตร์สำหรับการผลิตไบโอดีเซลโดยใช้ปฏิกิริยาเอนไซม์เพื่ออธิบายแบบจำลองทางคณิตศาสตร์ที่ง่ายสำหรับตัวเร่งปฏิกิริยาเอนไซม์ปฏิกิริยาtransesterification ของน้ำมันสบู่ดำตามสมมติฐานที่ได้รับ(A) Biocatlytic เร่ง transesterification (โดยใช้เอนไซม์ไลเปส) ของน้ำมันสบู่ดำที่มีเครื่องดื่มแอลกอฮอล์(AL ) อาจมีการเสนอที่จะเกี่ยวข้องกับกลไกสองขั้นตอน ขั้นตอนแรกประกอบด้วยการย่อยสลายของน้ำมันสบู่ดำหรือ TG acylated การผลิตเอนไซม์(ACE) และการเปิดตัวของกลีเซอรอลเป็นโดยเฉพาะสินค้าที่ผ่านที่ซับซ้อนC1 DE: TGÞ ที่นี่ E ย่อมาจากเอนไซม์. ขั้นที่สองคือ esterification ของเมทานอล (AL) กับเอซในรูปแบบสินค้าที่ต้องการเช่น BD กับการเปิดตัวของเอนไซม์ฟรี(E) ผ่าน C2 ซับซ้อนสองDace. alth [45] (A2) บัญชีสำหรับการยับยั้งการเกิดปฏิกิริยาโดยเมทานอลในที่จลนศาสตร์ระบบที่นี่ยับยั้งการแข่งขันจะถือว่าเมื่อโมเลกุลAL รวมกับ E ในการผลิตที่ซับซ้อนสิ้นตายซับซ้อนC3 (E.AL) [18,22,26 46] ทุกขั้นตอนกลไกสำหรับการผลิตไบโอดีเซลสามารถแสดงโดยลำดับต่อไปของปฏิกิริยา: E þ TG r1 r1 ½E: TG? r2 r2 ACE þ GL; ACE þ AL r3 r3 ½AcE: AL r4 r4 E þ BD; E þ AL R5 R5 ½E: AL ?: ð6Þนี่r1 และ r1; r2 และ r2 เป็นค่าคงที่อัตราการก่อย้อนกลับของC1 ซับซ้อนเอนไซม์ acylated และจากกลีเซอรอลผลิตภัณฑ์ตามลำดับในขั้นตอนแรกของไบโอดีเซลก่อ r3 และ r3; r4 r4 และเป็นค่าคงที่อัตราการก่อตัวของC2 ย้อนกลับที่มีความซับซ้อนและการพัฒนาไบโอดีเซลตามลำดับในขั้นตอนสุดท้าย สำหรับปฏิกิริยายับยั้ง, R5 R5 และเป็นค่าคงที่อัตราการก่อย้อนกลับของตายสิ้นซับซ้อนC3. (A3) ในฐานะที่เป็นโดยไม่ต้องกวนอัตราการเกิดปฏิกิริยา transesterification เป็นช้ามากและไม่สมบูรณ์เนื่องจากการต้านทานการถ่ายเทมวล[26,47] ดังนั้นผลของการกวนจะรวมอยู่ในปฏิกิริยาจลนศาสตร์ในแฟชั่นโลจิสติกเป็นรูปแบบก่อนหน้านี้โดยทั้งสองสำนวนksxB 1 xB Bmax? ? และ ksxC3 1 xC3 C3max? ? ที่ C3max. เป็นความเข้มข้นสูงสุดของความซับซ้อนสิ้นตายเราแสดงถึงความเข้มข้นของ TG นั้น E; เอซ; C1; C2; C3; อัล; BD และGL เป็น XT; XE; xAcE; xC1; xC2; xC3; xA; xB และ xG ตามลำดับ ตอนนี้จากสมมติฐานข้างต้นด้วยกลไกการเกิดปฏิกิริยาดังกล่าวข้างต้นตามตามกฎหมายของการกระทำมวลเราสามารถกำหนดชุดของสมการเชิงอนุพันธ์ได้รับดังนี้DXE dt ¼ r1xT XE þ r1xC1 þ r4xC2 r4xExB r5xExA þ r5xC3; DXT dt ¼ r1xT XE þ r1xC1 r5xAcE þ r5xC3; dxAcE dt ¼ r2xC1 r2xAcExG r3xF xA þ r3xC2 r5xAcE þ r5xC3; DXB dt ¼ r4xC2 r4xExB þ ksxB 1 xB Bmax; เอฟเอ Basir et al, เชื้อเพลิง / 158 (2015) 503-511 505 DXA dt ¼ r3xAcExA þ r3xC2 r5xExA þ r5xC3; dxC1 dt ¼ r1xT XE r1xC1 r2xC1 þ r2xAcExG; dxC2 dt ¼ r3xAcExA r3xC2 r4xC2 þ r4xExB; dxC3 dt ¼ r5xExA r5xC3 þ ksxC3 1 xC3 C3max; DXG dt ¼ r2xC1 r2xAcExG; ð7Þกับเงื่อนไขเริ่มต้นxEð0Þ¼ xE0; xC1ð0Þ¼ 0; xAcEð0Þ¼ 0; xAð0Þ¼ xA0; XT ð0Þ¼ xT0; xC2ð0Þ¼ 0; xBð0Þ¼ 0; xGð0Þ¼ 0 xC3ð0Þ¼ 0: ð8Þระยะที่ผ่านมาในการแสดงออกของDXB dt และ dxC3 dt, เป็นเพียงมวลอัตราการถ่ายโอนหรือทางอ้อมผลงานของการถ่ายโอนมวลความต้านทานต่อความต้านทานโดยรวม นี่ C3max เป็นความเข้มข้นสูงสุดของผู้ตายสิ้นxC3 ซับซ้อน. 4 ปัญหาการควบคุมที่เหมาะสมสำหรับระบบปฏิกิริยาของเอนไซม์วัตถุประสงค์หลักของเราคือการควบคุมอัตราการกวน(เช่นมวลอัตราการถ่ายโอน) เพื่อที่เราจะได้รับการผลิตสูงสุดของไบโอดีเซล. นอกจากนี้ก็เป็นวัตถุของเราเพื่อให้การทำงานของค่าใช้จ่ายที่ต่ำที่สุดเท่าที่เป็นไปได้ . เราใช้การควบคุมตัวแปรuðtÞซึ่งแสดงให้เห็นถึงกระตุ้นกวนอินพุทที่เวลา t ความพึงพอใจ 0 6 uðtÞ 6 1. uðtÞนี่แสดงให้เห็นถึงการควบคุมการป้อนข้อมูลที่มีค่าปกติจะอยู่ระหว่าง0 และ 1 [48] นอกจากนี้uðtÞ¼ 1 แสดงให้เห็นถึงการใช้งานสูงสุดของกวนและuðtÞ? 0 ซึ่ง. หมายถึงไม่มีการเปลี่ยนแปลงในการหมุนกวนจากสมมติฐานการควบคุมระบบการเหนี่ยวนำให้สอดคล้องกับระบบ (7) จะเป็น: DXE dt ¼ r1xT XE þ r1xC1 þ r4xC2 r4xExB; DXT dt ¼ r1xT XE þ r1xC1 r5xAcE þ r5xC3; dxAcE dt ¼ r2xC1 r2xAcExG r3xF xA þ r3xC2 r5xAcE þ r5xC3; DXB dt ¼ k4xC2 r4xExB þuðtÞksxB 1 xB Bmax; DXA dt ¼ r3xAcExA þ r3xC2; dxC1 dt ¼ r1xT XE r1xC1 r2xC1 þ r2xAcExG; dxC2 dt ¼ r3xAcExA r3xC2 r4xC2 þ r4xExB; dxC3 dt ¼ r5xAcE r5xC3 þuðtÞksxC3 1 xC3 C3max; DXG dt ¼ r2xC1 r2xAcExG; ð9Þกับเงื่อนไขเริ่มต้นตามที่กำหนดโดยสมการ (8). ฟังก์ชั่นค่าใช้จ่ายเป็นสูตรจึงเป็นJ½uðtÞ? ¼ Z TF t0 Pu2 ðtÞ Qx2 BðtÞþ RX2 C3 ชั่วโมง idt ð10Þที่นี่เราต้องการที่จะเพิ่มการผลิตไบโอดีเซลและลดการยับยั้งการทำงานของเอนไซม์คือลดความเข้มข้นของตายสิ้นxC3 ซับซ้อน พารามิเตอร์ PD> 0 คือน้ำหนักคงที่ที่ได้รับประโยชน์จากต้นทุนการผลิตและQ และ R ที่มีการปรับตัวคูณ. ตอนนี้เราต้องการที่จะหาการควบคุมที่ดีที่สุด u? เช่นที่Jðu? Þ¼นาทีJðuÞ: ยู 2 ท่านที่U คือชุดควบคุมที่ยอมรับที่กำหนดโดยU ¼fuðtÞ: uðtÞเป็นที่วัด; 0 6 uðtÞ 6 1 2 t ½ti; TF กรัม? ð11Þที่นี่เราใช้หลักการขั้นต่ำ Pontryagin [49,50] เพื่อหายูðtÞ?. มิลโตเนียนเป็นสูตรที่เป็นH ¼ Pu2 ðtÞ Qx2 BðtÞþ RX2 C3 þX9i¼1 nifi; ð12Þที่ n1 N9 เป็นตัวแปร adjoint และ fi, i = 1, 2, ... , 9 ขวาด้านข้างของสมการ (9) คือฉ 1 ¼ r1xT XE þ r1xC1 þ r4xC2 r4xExB ฯลฯทฤษฎีบท หากได้รับการควบคุมที่ดีที่สุดยูðtÞและการแก้ปัญหา? x? E; ? x T; ? x ACE; x? B; x?; x? C1; x? C2; x? C3? ? ระบบที่สอดคล้องกัน (. สม 7) ลดJðuÞมากกว่า U นั้นมีอยู่ตัวแปร adjoint n1 N9 ซึ่งความพึงพอใจของสมการต่อไปนี้DN1 dt ¼ r1xT DN1 n6Þþr4xBðn4n7Þþr5xAðn5n8Þ; dn2 dt ¼r1xEðn1þ n2 n6Þ; dn3 dt ¼r3xAðn5n7Þþ r5xT ðn8Þþr2xGðn3þ N9 n6Þ; dn4 dt ¼ 2QxB þ N4 r4xE þksxBuðtÞ Bmax ksuðtÞ 1 xB Bmax? n7r4xE; DN5 dt ¼r3xAcEðn5n7Þþr5xEðn1n8Þ; DN6 dt ¼r1ðn6n2Þþr2ðn6 n3 n9Þ; dn7 dt ¼ r4n4 r3 þ r3n7; DN8 dt ¼r5ðn8 n2 n3Þþn8uðtÞks 1 xC3 C3max þn8uðtÞksxC3 1 1 C3max; dn9 dt ¼r2xAcEðn3 n6 þn9Þ: ð13Þพร้อมกับสภาพtransversality niðtfÞ¼ 0 สำหรับฉัน¼ 1 ถึง 9? นอกจากนี้ยูðtÞสามารถเขียนเป็น, u? ðtÞ¼แม่

ผลลัพธ์ (ไทย) 3:

2 . แบบจำลองทางคณิตศาสตร์สำหรับการผลิตไบโอดีเซลโดยใช้ตัวเร่งปฏิกิริยาเคมี

อธิบายแบบจำลองทางคณิตศาสตร์อย่างง่ายสำหรับการเกิดปฏิกิริยาด่าง
กระบวนการทรานส์เอสเทอริฟิเคชั่นของสบู่ดำน้ำมัน สมมติฐานที่ได้รับการทำดังต่อไปนี้
.
( A1 ) ตัวเร่งปฏิกิริยาทรานส์เอสเทอริฟิเคชันของน้ำมันพืชกับเมทานอลประกอบด้วย
3 = [ 13,40 และปฏิกิริยาย้อนกลับ ] .
( A2 ) ตั้งแต่ มีเพียงน้ำเล็กน้อย ( 0.2 %w / w ) ในปฏิกิริยา
ผสมและไม่มีกรดไขมันอิสระที่ตรวจพบในระบบการย่อย อาจจะกระจอก
ปฏิกิริยา [ 41 ] .
( A3 ) ส่วนผสมของสบู่ดำเศษอาหาร น้ำมันที่ประกอบด้วย 1 เปอร์เซ็นต์โดยน้ำหนัก
ffas . TG , BD , FFA อาจมีส่วนร่วมในการสปอนนิฟิเคชั่นด้วย
ตัวเกาะ [ 42 ] .
( A4 ) ฐานเร่งปฏิกิริยาสำหรับกระบวนการทรานส์เอสเทอริฟิเคชั่น
การผลิตไบโอดีเซลน้ำมันจากสบู่ดำกับเมทานอล ( AL ) เกี่ยวข้องกับ
3 , ขั้นตอน [ 40 ] อยู่ในขั้นตอนที่กลับก่อน
ไตรกลีเซอไรด์ ( TG ) ทำปฏิกิริยากับแอลกอฮอล์ ( AL ) เพื่อผลิตไดกลีเซอไรด์
( DG ) ซึ่งต่อไปจะทำปฏิกิริยากับ อัล ผลิต monoglycerides
( มิลลิกรัม ) ในขั้นตอนการ 2 ในที่สุด มก
ปฏิกิริยากับอัลให้กลีเซอรอล ( GL ) โดยผลิตภัณฑ์ ในแต่ละ
ปฏิกิริยาขั้นตอนหนึ่งโมเลกุลของไบโอดีเซล ( BD ) ผลิต
สำหรับแต่ละโมเลกุลของ อัล ผลาญ ขั้นตอนที่ 1 จะแสดงปฏิกิริยา



( TG þอัล K1 K2

DG DG þ BD ต่อไปนี้ ; K3


þลแห่งมก. þ BD ;

เสียใจต่อþอัล K6

GL þ BD :
ð 1 Þ
มันแสดงข้างบนที่แกล้วกล้าจลนศาสตร์ปฏิกิริยา โครงการ
ประกอบด้วยสามปฏิกิริยาย้อนกลับได้ นอกจากนั้นยังมี
TG , FFA BD สปอนนิฟิเคชั่นและปฏิกิริยามันเป็นที่แสดงด้านล่างว่าจลนศาสตร์ปฏิกิริยาสปอนนิฟิเคชั่น

โครงการประกอบด้วยสามส่วนใหญ่ปฏิกิริยา [ 4 ] .
TG þ k
K7
s
þ GL BD þ k

ของกานต์แอร์þอัล
F þ k

þ K9 s W :

แต่Þð 2 บ. บาซีร์ et al . 158 / เชื้อเพลิง ( 2015 ) 503 – 511
( A5 ) K1 K3 ; ; เสียใจเป็นค่าคงที่อัตราการเกิดปฏิกิริยาไปข้างหน้าและ K2 ; K6
แห่งและมีย้อนกลับอัตราการเกิดปฏิกิริยาค่าคงที่ ค่า K1 K6
เพื่อจะแสดงในตารางที่ 1