ฟองแคเวิเทชัน: กำเนิด การเติบโต และการระเบิด

การเดินทางของฟองแคเวิเทชันเริ่มต้นที่จุดที่เรียกว่าการนิวคลีเอชัน ในของเหลวทั่วไปมักมีฟองแก๊สขนาดเล็กละลายอยู่แล้ว รวมถึงอนุภาคของแข็งขนาดเล็กและรอยขรุขระบนพื้นผิวภาชนะ สิ่งเหล่านี้ทำหน้าที่เป็นแหล่งกำเนิดเริ่มต้นให้ฟองขยายตัวขึ้นมา

เมื่อฟองเกิดขึ้นแล้ว วัฏจักรการอัดและขยายของคลื่นเสียงจะกระทำต่อฟองอย่างต่อเนื่อง ในระหว่างการขยายตัว แรงดันภายในฟองต่ำลงทำให้แก๊สจากของเหลวรอบข้างแพร่เข้ามา และในระหว่างการหดตัว แรงดันที่สูงขึ้นจะไล่แก๊สออกไป แต่เนื่องจากพื้นที่ผิวของฟองขณะขยายตัวมากกว่าขณะหดตัว จึงมีแก๊สไหลเข้ามาสุทธิมากกว่าที่ออกไป ปรากฏการณ์นี้เรียกว่า rectified diffusion และทำให้ฟองโตขึ้นเรื่อยๆ ในแต่ละรอบ

เมื่อฟองโตเกินขนาดวิกฤต มันจะไม่สามารถเติบโตต่อไปได้อย่างสมดุลอีกต่อไป และจะพังทลายลงอย่างรุนแรงในระหว่างวัฏจักรการอัด เราเรียกการพังทลายแบบรวดเร็วและรุนแรงนี้ว่า transient cavitation หรือ inertial cavitation ซึ่งเป็นต้นเหตุของผลทางเคมีและกายภาพที่น่าสนใจทั้งหมด ในทางตรงข้าม หากฟองสามารถแกว่งอยู่ที่รัศมีเฉลี่ยคงที่ได้หลายรอบโดยไม่พัง เราเรียกว่า stable cavitation หรือ non-inertial cavitation ซึ่งสังเกตพบในระบบที่มีแรงดันเสียงต่ำกว่า

ผลกระทบทางกายภาพที่เกิดจากการพังทลายของฟอง

การพังทลายของฟองแคเวิเทชันไม่ได้เป็นเหตุการณ์ที่เกิดขึ้นเงียบๆ หากฟองพังทลายใกล้กับพื้นผิวของแข็ง การพังทลายจะไม่สมมาตร ทำให้เกิดไมโครเจ็ต ซึ่งเป็นกระแสของเหลวความเร็วสูงพุ่งทะลุเข้าสู่พื้นผิว ด้วยความเร็วที่สามารถสูงถึง 4,000 เมตรต่อวินาที แรงกระแทกนี้คือสาเหตุของปรากฏการณ์ที่รู้จักกันดีในอุตสาหกรรมการต่อเรือมาตั้งแต่ศตวรรษที่ 19 นั่นคือการกร่อนเป็นรูปแอ่งบนใบพัดเรือ

นอกจากไมโครเจ็ต การพังทลายยังปล่อยคลื่นกระแทก (shockwave) ออกมาด้วย คลื่นกระแทกเหล่านี้ช่วยเพิ่มการผสมในระดับจุลภาค (micromixing) และเร่งการถ่ายเทมวลสารได้อย่างมีนัยสำคัญ ในการประยุกต์บางอย่าง ยังสามารถใช้คลื่นกระแทกเหล่านี้เพื่อทำให้ผลึกของแข็งแตกเป็นชิ้นเล็กๆ ในกระบวนการที่เรียกว่า sonofragmentation หรือนำไปใช้ในการทำความสะอาดอุปกรณ์ด้วยคลื่นเสียง

ปัจจัยที่กำหนดความรุนแรงของ Acoustic Cavitation

ความเข้าใจว่าอะไรทำให้แคเวิเทชันรุนแรงขึ้นหรืออ่อนลงมีความสำคัญอย่างยิ่งต่อการออกแบบระบบ ปัจจัยหลักที่ส่งผลได้แก่

แก๊สที่ละลายในของเหลว

แก๊สที่ละลายอยู่ในของเหลวมีบทบาทคู่ คือทั้งช่วยสร้างแหล่งนิวคลีเอชันและส่งผลต่อความเข้มของการพังทลาย กฎทั่วไปที่ใช้กันคือ แก๊สที่มีค่า specific heat ratio สูงกว่าจะให้ผลของแคเวิเทชันที่รุนแรงกว่า ซึ่งหมายความว่าแก๊สเชิงอะตอมเดี่ยว เช่น อาร์กอน จะให้การปลดปล่อยพลังงานที่มีประสิทธิภาพกว่าแก๊สแบบอะตอมคู่ เช่น ไนโตรเจน

อุณหภูมิและความดันของระบบ

อุณหภูมิที่สูงขึ้นทำให้ค่าแรงดันไอของของเหลวสูงขึ้น ส่งผลให้แคเวิเทชันเกิดขึ้นได้ง่ายขึ้น แต่ขณะเดียวกันก็ลดความรุนแรงของการพังทลายลงเพราะไอน้ำภายในฟองต้านทานการหดตัว ในทางกลับกัน การเพิ่มความดันรวมของระบบทำให้เริ่มแคเวิเทชันยากขึ้น แต่เมื่อเกิดขึ้นแล้วจะรุนแรงกว่า การปรับสมดุลระหว่างสองปัจจัยนี้จึงเป็นโจทย์สำคัญในการออกแบบ sonochemical reactor

สมบัติทางกายภาพของของเหลว

ของเหลวที่มีความหนืดสูงมักให้ความเข้มของแคเวิเทชันที่สูงกว่า ขณะที่ของเหลวที่ระเหยง่ายมากอาจไม่เกิดแคเวิเทชันที่มีประสิทธิภาพเลย แรงตึงผิวก็มีผลด้วยเช่นกัน เนื่องจากเกี่ยวข้องโดยตรงกับแรง Laplace ที่ต้านการขยายตัวของฟอง

ความถี่และกำลังของอัลตราซาวนด์

ความถี่สูงขึ้นทำให้วัฏจักรการอัดและขยายสั้นลง ฟองจึงมีเวลาเติบโตน้อยลง ส่งผลให้แคเวิเทชันแบบ transient ลดลง แต่ก็เพิ่มอัตราการสร้างฟองด้วย งานวิจัยหลายชิ้นพบว่าการเพิ่มความถี่มักให้ฟองแบบ stable มากขึ้น ซึ่งเหมาะสำหรับบางการประยุกต์ใช้ ส่วนกำลังเสียงก็มีค่าเหมาะสมที่ให้ประสิทธิภาพสูงสุด เพราะหากกำลังสูงเกินไป ฟองที่เกิดขึ้นหนาแน่นรอบปลาย transducer อาจกลับมาบดบังการส่งพลังงานไปยังส่วนอื่นของระบบ

แบบจำลองเชิงตัวเลขสำหรับระบบอัลตราซาวนด์

การจำลองระบบอัลตราซาวนด์ด้วยคอมพิวเตอร์ช่วยให้นักวิจัยเข้าใจการกระจายความเข้มของเสียงในระบบ ซึ่งเป็นเรื่องที่ยากมากหากจะวัดด้วยวิธีทดลองอย่างเดียว ไม่ว่าจะเป็น calorimetry การทดสอบการกัดกร่อน หรือการวิเคราะห์อัตราการเกิดปฏิกิริยา การจำลองช่วยลดทรัพยากรที่ใช้ในการพัฒนาและออกแบบเทคโนโลยีใหม่ได้มาก

แบบจำลองเชิงเส้น: สมการ Helmholtz

จุดเริ่มต้นที่นิยมใช้มากที่สุดคือสมการ Helmholtz ซึ่งอธิบายการแพร่กระจายของคลื่นเสียงในของเหลวที่เป็นเนื้อเดียวกัน ในกรณีที่ระบบทำงานที่ความถี่เดียวคงที่ สมการนี้จะลดรูปลงเหลือ

∇²p + k²p = 0

โดยที่ p คือแอมพลิจูดของแรงดันเสียงที่ไม่ขึ้นกับเวลา และ k คือเลขคลื่น (wave number) ซึ่งขึ้นกับความถี่เชิงมุมและความเร็วเสียงในตัวกลาง

สมการ Helmholtz มีข้อดีสำคัญคือแก้ได้ไม่ยาก มีให้ใช้ใน FEM software เชิงพาณิชย์หลายตัว และให้ข้อมูลที่มีประโยชน์ว่าแรงดันเสียงสูงสุดอยู่ที่ตำแหน่งใด ซึ่งช่วยประเมินว่าแคเวิเทชันจะเกิดขึ้นบริเวณไหนได้ อย่างไรก็ตาม ข้อจำกัดสำคัญคือสมการนี้สมมติให้การแพร่กระจายเป็นเชิงเส้น ซึ่งไม่สอดคล้องกับความเป็นจริงในระบบที่มีฟองแคเวิเทชันอยู่มาก

เงื่อนไขขอบเขต: กุญแจสู่ความแม่นยำ

ส่วนที่ซับซ้อนในการจำลองระบบอัลตราซาวนด์คือการกำหนดเงื่อนไขขอบเขต (boundary conditions) ให้สมจริง งานวิจัยในช่วงแรกมักสมมติให้ผนังภาชนะเป็น infinitely rigid (ไม่เปลี่ยนรูปเลย) หรือ infinitely soft (เปลี่ยนรูปได้อย่างอิสระ) ทั้งสองอย่างเป็นการลดรูปปัญหาที่ผิดความเป็นจริง

ผลวิจัยพบว่าสำหรับผนังกระจกบาง (2 มม.) สมบัติใกล้เคียงกับขอบเขตแบบ soft ส่วนผนังหนา (7 มม.) ใกล้เคียงกับขอบเขตแบบ rigid มากกว่า และการใช้แบบจำลอง fluid-structure interaction ที่จำลองการเสียรูปของผนังด้วยจะให้ผลที่แตกต่างออกไปอย่างมีนัยสำคัญในด้านตำแหน่งของ resonance frequencies และ antinode สำหรับ transducer เองก็มักกำหนดเป็น pressure source หรือ displacement amplitude boundary แทนที่จะจำลองโครงสร้างทางกายภาพทั้งหมดของ piezoelectric element

แบบจำลองเชิงไม่เส้น: เมื่อ Helmholtz ไม่เพียงพอ

เมื่อแคเวิเทชันเริ่มเกิดขึ้น ฟองจะเปลี่ยนสมบัติเชิงเสียงของของเหลวอย่างมีนัยสำคัญ ตัวอย่างเช่น ความเร็วเสียงในน้ำปกติอยู่ที่ประมาณ 1,500 เมตรต่อวินาที แต่ในของเหลวที่มีฟองปะปนอยู่จำนวนมาก ความเร็วเสียงอาจลดลงเหลือเพียง 20 เมตรต่อวินาที ความแตกต่างระดับนี้ทำให้การใช้แบบจำลองเชิงเส้นให้ผลที่คลาดเคลื่อนมาก

แบบจำลองที่ก้าวหน้ากว่า เช่น สมการ Caflisch ซึ่งอธิบายการอนุรักษ์มวลและโมเมนตัมของของเหลวที่มีฟอง นำเสนอแนวทางที่ครอบคลุมกว่า แต่ก็มีความซับซ้อนสูงและต้องการทรัพยากรการคำนวณมาก โดยเฉพาะเมื่อต้องจำลองระบบที่มี inertial cavitation รุนแรง อีกแนวทางหนึ่งที่ได้รับความนิยมคือ nonlinear Helmholtz equation ของ Louisnard ซึ่งคำนวณ attenuation โดยตรงจากพลังงานที่สูญเสียในฟองแต่ละฟอง แนวทางนี้สามารถทำนายโครงสร้างฟองแบบกรวย (conical bubble structures) ใต้ ultrasonic horn ได้ ในขณะที่แบบจำลองเชิงเส้นทำไม่ได้

ตาราง: เปรียบเทียบแบบจำลองหลักในการจำลองระบบอัลตราซาวนด์

ตารางสรุปจาก Tiong et al. (2025): เปรียบเทียบแบบจำลองหลักสำหรับระบบอัลตราซาวนด์

ซอฟต์แวร์และเครื่องมือสำหรับการจำลองระบบอัลตราซาวนด์

การพัฒนา FEM solver เชิงพาณิชย์ที่มี acoustic module ในตัวได้เปิดประตูให้นักวิจัยหลายกลุ่มสามารถจำลองระบบเสียงได้สะดวกขึ้นมาก

ซอฟต์แวร์ที่ถูกใช้มากที่สุดในงานวิจัยปัจจุบันคือ COMSOL Multiphysics (เดิมชื่อ FEMLAB) เพราะมี acoustic module พร้อมใช้งาน และยังสามารถเชื่อมโยงฟิสิกส์หลายระบบไว้ด้วยกัน เช่น เสียง การไหล ความร้อน และปฏิกิริยาเคมี ในการศึกษาระบบ sonochemical reactor ยุคใหม่ มีการใช้ COMSOL จำลองทั้งสนามแรงดันเสียง โปรไฟล์การไหล และจลนพลศาสตร์ปฏิกิริยาไว้ในการจำลองเดียวกัน

นอกจาก COMSOL ยังมีซอฟต์แวร์อื่นที่ถูกนำมาใช้ เช่น ATILA code สำหรับระบบ textile treatment, PZflex code สำหรับ cylindrical reactor, Sysnoise สำหรับการจำลองสนามเสียงก่อนนำผลไปใช้ใน Fluent เพื่อวิเคราะห์การไหล และ Aquilon สำหรับการจำลองระบบ 2 เฟสที่มีฟองอยู่ด้วย

ความก้าวหน้าล่าสุดและทิศทางในอนาคต

ในช่วงทศวรรษที่ผ่านมา มีพัฒนาการสำคัญหลายประการในการจำลองระบบอัลตราซาวนด์ กลุ่มวิจัยของ Louisnard ได้พัฒนาแบบจำลองที่ทำนาย acoustic attenuation จากการพังทลายของฟอง โดยสรุปว่าแนวทางที่ใช้หลักการ primary Bjerknes force ให้การแทนค่าสภาพแวดล้อมเชิงเสียงในชีวิตจริงได้ดีกว่าวิธีเดิม

กลุ่มวิจัยของ Vanhille ได้ทำการศึกษาระบบ nonlinear acoustic waves อย่างต่อเนื่อง ทั้งใน stationary และ transient phase รวมถึงงานที่พิจารณาผลของความร้อนต่อระบบเสียงด้วย งานวิจัยกลุ่มนี้ช่วยเสริมความเข้าใจเกี่ยวกับพฤติกรรมคลื่นที่ซับซ้อนในของเหลวที่มีฟอง

ในส่วนของทีมวิจัยที่มีบทบาทในการเขียนบทความทบทวนนี้เอง ได้แสดงให้เห็นว่าสามารถประยุกต์ใช้แบบจำลอง bubbly liquid แบบ semi-empirical ในการจำลองระบบ multi-frequency ได้เป็นครั้งแรก ซึ่งเปิดโอกาสให้นักออกแบบ sonoreactor ใช้ความถี่หลายตัวพร้อมกันเพื่อเพิ่มประสิทธิภาพของระบบ

ความท้าทายที่ยังคงอยู่คือการเปรียบเทียบผลการจำลองกับข้อมูลทดลองจริงในระดับ pilot scale รวมถึงการจำลองระบบ multi-phase ที่ซับซ้อน และการจับภาพผลของ nonlinear effects ในระบบขนาดใหญ่ที่ต้องการ scale-up สำหรับการใช้งานอุตสาหกรรม

บทสรุป: จากฟองอากาศจิ๋วสู่นวัตกรรมระดับอุตสาหกรรม

บทความทบทวนวรรณกรรมฉบับนี้พาเราเห็นเส้นทางพัฒนาการของการจำลองระบบ acoustic cavitation ตั้งแต่สมการ Rayleigh-Plesset ที่เป็นรากฐานในการอธิบายพลวัตของฟองเดี่ยว มาจนถึงแบบจำลองที่ซับซ้อนอย่าง Commander-Prosperetti และ nonlinear Helmholtz equation ที่สามารถรับมือกับระบบหลายเฟสและผลของฟองจำนวนมากได้

สิ่งที่โดดเด่นตลอดวิวัฒนาการนี้คือความตึงเครียดระหว่างความแม่นยำและความซับซ้อนในการคำนวณ แบบจำลองที่ง่ายกว่าช่วยให้ออกแบบและ optimize ระบบได้รวดเร็ว แต่อาจพลาดปรากฏการณ์สำคัญที่เกิดขึ้นในระบบจริง ขณะที่แบบจำลองที่ครอบคลุมกว่าต้องการทรัพยากรคำนวณสูงและมีความยุ่งยากในการตั้งค่า

สำหรับอุตสาหกรรมและนักวิจัยที่กำลังพัฒนา sonochemical reactor รุ่นต่อไป ความเข้าใจในข้อจำกัดและจุดแข็งของแต่ละแบบจำลองจึงเป็นสิ่งจำเป็น เพื่อเลือกวิธีการที่เหมาะสมที่สุดกับโจทย์ที่ต้องการตอบ ทั้งนี้ความก้าวหน้าในการจำลองเชิงตัวเลขยังคงเป็นกุญแจที่จะปลดล็อกการประยุกต์ใช้ acoustic cavitation ในระดับอุตสาหกรรมที่กว้างขวางขึ้นในอนาคต

เอกสารอ้างอิง (References)

Tiong, T. J., Chu, J. K., & Tan, K. W. (2025). Advancements in Acoustic Cavitation Modelling: Progress, Challenges, and Future Directions in Sonochemical Reactor Design. Ultrasonics Sonochemistry, 112, 107163.

Michlovitz, S. L. (2011). Modalities for Therapeutic Intervention (5th ed.). F.A. Davis Company.

Commander, K. W., & Prosperetti, A. (1989). Linear pressure waves in bubbly liquids: Comparison between theory and experiments. Journal of the Acoustical Society of America, 85(2), 732-746.

Louisnard, O. (2012). A simple model of ultrasound propagation in a cavitating liquid. Part I: Theory, nonlinear attenuation and traveling wave generation. Ultrasonics Sonochemistry, 19(1), 66-76.

Thompson, L. H., & Doraiswamy, L. K. (1999). Sonochemistry: Science and Engineering. Industrial & Engineering Chemistry Research, 38(4), 1215-1249.