Abstract:
การเปลี่ยนแปลงทางเคมีในปฏิกรณ์ฟลูอิดไดซ์เบดที่สภาวะไม่คงตัวและอุณหภูมิที่จะติดตามได้โดยอาศัยแบบจำลองคณิตศาสตร์ที่ประกอบด้วยสมการอนุรักษ์มวลร่วมกับสมการการเคลื่อนที่ของก๊าซและสมการการเคลื่อนที่ของอนุภาค ในสมการการเคลื่อนที่ของอนุภาคจะอาศัยแบบจำลองชนิดรายเม็ด (Distinct Element Method) เพื่อคำนวณแรงกระทำต่ออนุภาค ตัวอย่างปฏิกิริยาที่ใช้ในการจำลองแบบเป็นปฏิกิริยาการสลายก๊าซโอโซนบนตัวเร่งปฏิกิริยาเหล็กออกไซด์ ผลการจำลองแบบพบว่าปฏิกิริยาส่วนใหญ่จะเกิดขึ้นในบริเวณเบด อนุภาคตัวเร่งปฏิกิริยาหนาแน่น โดยการเปลี่ยนแปลงทางเคมีจะมีเพิ่มขี้นตามแนวความสูงและรัศมีปฏิกรณ์ที่เพิ่มขึ้น ข้อมูลที่ได้จากการจำลองแบบสอดคล้องกับผลการทดลองของ Piccinini และคณะในปี 1983 และให้ผลที่ใกล้เคียงกับผลการทดลองมากกว่าแบบจำลองชนิดสองส่วนของ Mathur และ Lim ในปี 1974 นอกจากนี้ได้ศึกษาอิทธิพลของความเร็วก๊าซป้อน ขนาดอนุภาคตัวเร่งปฏิกิริยาและขนาดปฏิกรณ์ที่มีผลต่อการเปลี่ยนแปลงทางเคมี เมื่อความเร็วของก๊าซป้อนเพิ่มขึ้นจาก 2.07 เป็น 4.13 เมตรต่อวินาที พบว่าจะเกิดการเปลี่ยนแปลงทางเคมีลดลงจาก 0.33 เป็น 0.20 ในการศึกษาระบบที่มีตัวเร่งปฏิกิริยาต่างขนาดกันคือ 3.1 4.4 และ 6.0 มิลลิเมตร การเปลี่ยนแปลงทางเคมีจะเกิดขึ้นมากในปฏิกรณ์ที่ใช้ตัวเร่งปฏิกิริยาขนาดเล็กกว่า การเปลี่ยนแปลงทางเคมีจะเพิ่มขึ้นจาก 0.18 เป็น 0.30 เมื่อปฏิกรณ์มีขนาดเส้นผ่านศูนย์กลางเพิ่มขึ้นจาก 0.10 เป็น 0.22 เมตรNumerical simulation models were developed for investigating of unsteady state mass transfer throughout an isothermal catalytic gas-solid fluidized bed reactor. The mathematical models were based on the equations of mass conservation coupling with equations of gas motion and particle motion. The realistic model of Distinct Element Method (DEM), which calculates the contact force on the individual particles, was used for prediction of particle motion. The decomposition of ozone on iron oxide catalyst in a fluidized bed reactor has been simulated. The unconverted ozone fraction in gas and particle phase as a function of time was illustrated. The simulation results obtained from this model and the two-region model developed by Mathur and Lim, 1974 were compared with experimental results of Piccinini et al. in 1983. It was found that the predictions of this model are in better agreement with the experimental data than that of the two-region model. In addition, the effects of superficial gas velocities, catalyst sizes and column sizes were studied. Increasing the gas velocities from 2.07 to 4.13 m/s, the averaged conversion decreases from 0.33 to 0.20. In the simulation of the beds containing the catalyst particles with the diameters of 3.1, 4.4 and 6.0 mn, the highest conversion was found in the bed containing smallest catalyst particles. The averaged conversion increases from 0.18 or 0.30 as the bed diameter increases from 0.10 m to 0.22m.