ccitonline.com

Assalamualaikum Wr. Wb.

Nama saya Muhammad Sabiq Nur Rachman (NPM 2306228724), mahasiswa dari kelas Metode Numerik 01. Laporan ini disusun dengan tujuan menerapkan dan mengintegrasikan berbagai metode numerik dalam kajian fenomena aliran fluida yang dipengaruhi oleh keberadaan vortex generator (VG). Meningkatnya tuntutan terhadap efisiensi energi dan optimalisasi performa sistem aerodinamika mendorong perlunya pemahaman yang lebih mendalam terhadap dinamika aliran. Dalam konteks ini, pemodelan berbasis Computational Fluid Dynamics (CFD) menawarkan pendekatan komputasional yang efektif untuk mengevaluasi fenomena fluida kompleks yang sulit diakses melalui metode eksperimental konvensional.

Dalam penyusunan laporan ini, saya mengadopsi kerangka berpikir Deep Awareness of I 5 (DAI5) sebagai pendekatan konseptual yang mendasari seluruh tahapan analisis. Pendekatan ini tidak hanya membantu dalam memahami aspek teknis dari simulasi, tetapi juga memberikan arah berpikir yang sistematis, reflektif, serta berorientasi pada pemecahan masalah. DAI5 menjadi kerangka kerja yang memfasilitasi proses mulai dari kesadaran individu, perumusan tujuan yang jelas, penetapan asumsi dasar, penyederhanaan fenomena kompleks menjadi model ideal, hingga perencanaan langkah-langkah aplikatif dalam proses simulasi.

1. Deep Awareness of I (Kesadaran Diri Mendalam)

Sebagai mahasiswa Teknik Mesin dengan minat khusus pada dinamika fluida komputasional, saya memiliki ketertarikan terhadap interaksi kompleks antara fluida dan struktur. Kesadaran ini mendorong saya untuk memahami bagaimana manipulasi aliran udara—khususnya melalui penerapan vortex generator—dapat meningkatkan efisiensi aerodinamis. Simulasi ini saya pandang sebagai medium penghubung antara teori yang dipelajari secara akademik dengan tantangan rekayasa yang bersifat praktis. Oleh karena itu, studi ini tidak hanya merupakan latihan akademik, tetapi juga sarana untuk melatih intuisi teknis dalam memvisualisasikan fenomena fluida yang tidak dapat diamati secara langsung serta untuk menyumbangkan solusi desain yang lebih efektif dan inovatif.

2. Intention (Tujuan dan Niat Awal)

Tujuan utama dalam studi ini adalah mengeksplorasi penerapan Finite Element Method (FEM) dalam simulasi aliran fluida di sekitar vortex generator menggunakan perangkat lunak Siemens STAR-CCM+. Tujuan khusus yang ingin dicapai meliputi:

• Analisis Pola Aliran: Mengkaji perubahan karakteristik aliran udara akibat kehadiran vortex generator, khususnya dalam hal distribusi kecepatan dan tekanan di sekitar permukaan objek aerodinamis.
• Evaluasi Efektivitas VG: Menilai kontribusi vortex generator dalam menunda atau mengurangi separasi lapisan batas, serta implikasinya terhadap peningkatan kinerja aerodinamis secara keseluruhan.
• Penguatan Pemahaman Numerik: Mengembangkan kemampuan analitis dan komputasional melalui pemanfaatan STAR-CCM+ sebagai alat bantu untuk memodelkan fenomena fluida kompleks yang sulit dijangkau melalui pengamatan eksperimental konvensional.

3. Initial Thinking (Pemikiran Awal dan Perumusan Masalah)

Langkah awal dalam memformulasikan permasalahan aliran fluida ini mengacu pada penerapan prinsip dasar fisika fluida ke dalam bentuk formulasi matematis. Persamaan-persamaan fundamental dalam Computational Fluid Dynamics (CFD) menjadi landasan bagi representasi numerik fenomena yang dikaji.

Salah satu persamaan utama yang digunakan adalah:

• Persamaan Kontinuitas, yang merepresentasikan hukum kekekalan massa. Untuk kasus aliran tunak (steady-state) dan inkompresibel (di mana densitas dianggap konstan), persamaan ini menyatakan bahwa laju massa yang masuk ke suatu volume kontrol harus sebanding dengan laju massa yang keluar. Dalam bentuk diferensial tiga dimensi, persamaan ini dinyatakan sebagai:

di mana u,v,w adalah komponen kecepatan dalam arah x,y,z.

Persamaan ini menjadi syarat penting untuk menjamin konservasi massa dalam simulasi aliran fluida yang dilakukan.

4. Idealization (Pemodelan dan Penyederhanaan Awal)

Dalam rangka mencapai keseimbangan antara akurasi fisik dan efisiensi komputasional, beberapa asumsi dan penyederhanaan diterapkan dalam proses pemodelan simulasi ini. Asumsi-asumsi tersebut dirancang untuk menyederhanakan kompleksitas masalah fisik tanpa menghilangkan karakteristik penting dari fenomena yang dikaji, yaitu:

• Aliran Tunak (Steady-State): Aliran diasumsikan bersifat tunak, yakni tidak mengalami perubahan terhadap waktu. Asumsi ini memfokuskan analisis pada kondisi aliran yang telah mencapai keseimbangan, sehingga memudahkan proses konvergensi numerik.
• Geometri Simetris (jika berlaku): Untuk kasus geometri yang memiliki bidang simetri, hanya setengah domain fisik yang dimodelkan. Hal ini memungkinkan pengurangan signifikan pada beban komputasi, dengan menerapkan kondisi batas simetri pada bidang tersebut.
• Properti Fluida Konstan: Karakteristik fluida (dalam hal ini udara) seperti densitas dan viskositas dianggap konstan, dengan asumsi bahwa tidak terdapat variasi suhu atau tekanan yang cukup signifikan untuk mempengaruhi properti termofisikanya.
• Pemodelan Turbulensi: Fenomena turbulensi disimulasikan menggunakan model turbulensi k-ε, yang merepresentasikan efek rata-rata turbulensi tanpa menyelesaikan setiap skala turbulen secara langsung. Pendekatan ini menawarkan efisiensi komputasional yang lebih tinggi sembari mempertahankan akurasi makroskopik dari perilaku aliran.

Penyederhanaan ini merupakan strategi yang umum digunakan dalam simulasi numerik untuk memungkinkan proses penyelesaian dalam waktu yang wajar serta tetap relevan terhadap tujuan analitis studi ini.

5. Instruction Set (Langkah-Langkah Simulasi)

Proses simulasi dilakukan secara bertahap melalui perangkat lunak Siemens STAR-CCM+, dengan alur kerja sebagai berikut:

• Pembuatan Geometri: Objek aerodinamis berupa profil menyerupai airfoil dengan penambahan vortex generator dimodelkan secara tiga dimensi. Geometri ini dibangun secara presisi dengan mempertimbangkan dimensi dan konfigurasi penempatan vortex generator, kemudian diimpor ke dalam STAR-CCM+ untuk keperluan simulasi.
• Pembuatan Mesh: Domain fluida di sekitar objek didiskretisasi menjadi elemen-elemen kecil menggunakan polyhedral mesh. Kerapatan mesh diperhalus secara lokal di area kritis seperti di sekitar vortex generator dan lapisan batas untuk menangkap gradien kecepatan dan tekanan secara lebih akurat.
• Penetapan Kondisi Batas (Boundary Conditions):
  • Inlet: Aliran udara masuk ditetapkan dengan kecepatan 70 m/s dan tingkat turbulensi sesuai parameter standar.
  • Outlet: Ditetapkan dengan kondisi tekanan gauge nol (0 Pa) sebagai batas keluaran.
  • Wall: Kondisi no-slip diterapkan pada seluruh permukaan padat, mengasumsikan kecepatan fluida relatif terhadap permukaan adalah nol.
  • Symmetry Plane: Pada geometri yang simetris, bidang tengah diberi kondisi batas simetri untuk mengurangi kompleksitas simulasi.
• Konfigurasi Model Fisika: Simulasi dilakukan dalam rezim steady-state dengan mengaktifkan model turbulensi k-ε, yang sesuai untuk aliran dengan karakteristik separasi dan interaksi dinding.
• Solusi dan Iterasi: Simulasi dijalankan dengan kriteria konvergensi berdasarkan nilai residual yang ditargetkan di bawah 10^{-5}hingga 10^{-6}, serta pemantauan kestabilan parameter penting seperti gaya aerodinamis. Jumlah iterasi yang dibutuhkan untuk mencapai konvergensi pada kasus ini berkisar 500 iterasi.
• Post-Processing: Hasil simulasi dianalisis secara visual dan kuantitatif. Data yang dikaji meliputi distribusi kecepatan dan tekanan, pola aliran (streamline), serta perhitungan gaya drag dan koefisien drag untuk menilai performa aerodinamis dari konfigurasi yang disimulasikan.

Analisis Hasil Simulasi

Bagian ini menyajikan interpretasi terhadap hasil-hasil utama dari simulasi CFD yang telah dilakukan, dengan fokus pada visualisasi pola aliran dan karakteristik distribusi parameter aliran di sekitar objek aerodinamis.

1. Visualisasi Pola Aliran dan Distribusi Kecepatan

Gambar berikut menyajikan representasi streamline yang menggambarkan pola aliran udara di sekitar permukaan objek aerodinamis yang dilengkapi dengan vortex generator (VG). Distribusi magnitudo kecepatan turut divisualisasikan untuk mengidentifikasi zona-zona percepatan dan perlambatan aliran.

Keberadaan vortex generator menghasilkan gangguan terkontrol pada aliran utama, membentuk struktur pusaran (vortex) yang membantu mencampurkan aliran lapisan batas dengan aliran utama. Hal ini terlihat dari perubahan arah dan intensitas streamline yang menunjukkan peningkatan energi kinetik di dekat permukaan, serta penundaan terjadinya separasi aliran. Daerah dengan gradien kecepatan tinggi teramati di sekitar trailing edge vortex generator, mengindikasikan keberhasilan mekanisme vortex dalam meningkatkan reattachment aliran.

Dari hasil visualisasi, pola streamline yang ditampilkan dengan warna hijau dan biru muda memperlihatkan karakteristik aliran yang kompleks di sekitar permukaan objek aerodinamis. Teridentifikasi zona percepatan aliran di bagian atas objek—ditunjukkan oleh warna hijau terang—dengan kecepatan mencapai sekitar ±20,8 m/s sesuai dengan skala legenda. Sebaliknya, di area yang terhalang geometri, terjadi perlambatan aliran yang ditandai dengan perubahan densitas streamline dan penurunan magnitudo kecepatan.

Salah satu temuan utama adalah terbentuknya struktur vorteks secara konsisten di daerah hilir masing-masing vortex generator (VG). Pola aliran yang berputar—yang mengindikasikan keberadaan vorteks—menunjukkan adanya fenomena pencampuran fluida secara intens di wilayah tersebut. Vorteks ini berfungsi sebagai media transfer momentum, dengan cara mengalirkan energi kinetik dari aliran bebas yang berkecepatan tinggi ke dalam lapisan batas yang kecepatan alirannya lebih rendah. Mekanisme ini berperan penting dalam menghambat atau menunda terjadinya separasi aliran dari permukaan padat, sehingga meningkatkan efisiensi aerodinamis dan mengurangi risiko fenomena stall pada profil aerodinamis.

2. Mekanisme Pembentukan Vorteks Lokal oleh Vortex Generator

Untuk mengkaji lebih dalam proses pembangkitan vorteks secara lokal, dilakukan visualisasi terfokus pada area di sekitar satu atau beberapa vortex generator. Representasi ini bertujuan untuk mengidentifikasi distribusi kecepatan dan arah aliran secara lebih rinci, serta bagaimana struktur vorteks terbentuk akibat gradien tekanan dan interaksi antara aliran utama dengan geometri VG.

Hasil menunjukkan bahwa VG memicu gradien kecepatan lateral yang signifikan, menghasilkan pusaran aliran (streamwise vortex) yang stabil dan terarah. Vorteks ini muncul segera setelah trailing edge VG dan menunjukkan kecenderungan bertahan dalam jarak tertentu di sepanjang permukaan, tergantung pada bentuk dan orientasi VG yang digunakan. Pola ini menegaskan bahwa vortex generator secara efektif memodifikasi struktur aliran di wilayah dekat dinding (near-wall region), yang secara langsung berdampak pada distribusi tegangan gesek dan tekanan di permukaan.

3. Distribusi Tekanan

Meskipun tidak disajikan secara eksplisit dalam bentuk kontur tekanan permukaan, legenda tekanan pada hasil visualisasi kecepatan (dengan rentang −527 Pa hingga +228 Pa) memberikan indikasi penting mengenai karakteristik distribusi tekanan di dalam domain simulasi.

Nilai tekanan maksimum, yang mencapai sekitar +228 Pa, umumnya terkonsentrasi di wilayah stagnasi—yaitu pada bagian depan leading edge objek aerodinamis—serta di sisi depan vortex generator. Hal ini konsisten dengan fenomena perlambatan aliran yang menghasilkan peningkatan tekanan dinamis di area tersebut. Sebaliknya, tekanan minimum hingga sekitar −527 Pa teridentifikasi pada area permukaan atas yang melengkung, di mana percepatan aliran terjadi, serta pada bagian hilir vortex generator sebagai akibat dari terbentuknya zona hisap oleh struktur vorteks.

Distribusi tekanan ini berimplikasi langsung terhadap gaya-gaya aerodinamis yang bekerja pada objek, khususnya gaya angkat (lift) dan hambat (drag). Variasi tekanan yang tajam antara sisi hulu dan hilir vortex generator menunjukkan kontribusi signifikan VG dalam memodifikasi medan tekanan di sekitarnya.

4. Analisis Konvergensi Simulasi

Stabilitas dan validitas hasil simulasi CFD sangat ditentukan oleh tingkat konvergensi yang dicapai selama proses iterasi numerik. Data numerik yang diperoleh dari Siemens Star CCM+ memperlihatkan perilaku residual untuk setiap persamaan yang diselesaikan—termasuk kontinuitas, momentum, dan parameter turbulensi (k dan ε).

Grafik residual menunjukkan bahwa seluruh parameter mengalami penurunan nilai secara progresif dan mencapai ambang batas konvergensi yang telah ditentukan (biasanya <10⁻⁵), menunjukkan bahwa solusi telah stabil dan bebas dari fluktuasi numerik yang signifikan. Selain residual, parameter fisik seperti gaya aerodinamis dan distribusi tekanan juga menunjukkan kestabilan nilai pada iterasi akhir, memperkuat validitas konvergensi.

Dengan demikian, dapat disimpulkan bahwa proses simulasi telah mencapai kondisi tunak yang layak untuk dianalisis lebih lanjut dan digunakan sebagai dasar interpretasi teknis.

Analisis Konvergensi Simulasi (Lanjutan)

Data yang ditampilkan pada Gambar 4 merepresentasikan nilai residual untuk sejumlah persamaan konservasi, termasuk kontinuitas, momentum (sumbu X, Y, dan Z), energi kinetik turbulen (TKE), dan laju disipasi turbulen (ε atau TDR), serta parameter hasil seperti Drag Force dan Drag Coefficient. Meskipun hanya menampilkan sebagian kecil dari rentang iterasi (iterasi ke-491 hingga 500), semua residual menunjukkan penurunan hingga ke orde 10⁻⁶ hingga 10⁻⁷, yang menandakan bahwa solusi numerik telah mencapai tingkat kestabilan yang sangat tinggi.

Stabilitas nilai drag coefficient dan gaya hambat (drag force) pada iterasi-iterasi akhir mengindikasikan tidak adanya fluktuasi signifikan, yang merupakan tanda bahwa sistem telah mencapai kondisi tunak (steady state) dan hasil simulasi dapat dianggap konvergen secara numerik. Dengan tingkat konvergensi yang tinggi ini, dapat disimpulkan bahwa simulasi menghasilkan solusi yang akurat dan dapat diandalkan dalam merepresentasikan fenomena fisik yang dianalisis.

Kesimpulan

Berdasarkan seluruh tahapan simulasi dan analisis yang dilakukan, dapat disimpulkan bahwa pendekatan numerik berbasis Finite Element Method (FEM) dengan perangkat lunak Siemens Star CCM+ telah berhasil digunakan untuk memodelkan interaksi aliran udara dengan vortex generator (VG) secara efektif.

Hasil simulasi menunjukkan bahwa penambahan VG secara signifikan memengaruhi struktur aliran, terutama melalui pembangkitan vorteks kuat di bagian hilir. Vorteks ini berperan penting dalam menciptakan pencampuran antara aliran berenergi tinggi di luar lapisan batas dengan aliran yang berenergi rendah di dekat permukaan padat. Mekanisme ini terbukti dapat menunda atau menghambat terjadinya separasi aliran, yang berdampak langsung pada peningkatan performa aerodinamis suatu sistem.

Distribusi kecepatan dan tekanan yang diperoleh dari simulasi konsisten dengan prinsip-prinsip dasar aerodinamika. Hal ini diperkuat oleh indikator konvergensi numerik yang sangat baik, ditunjukkan oleh residual yang rendah dan parameter aerodinamis yang stabil, yang secara keseluruhan mendukung validitas hasil simulasi.

Dengan demikian, studi ini menegaskan bahwa vortex generator merupakan elemen pasif yang efektif dalam pengendalian aliran, serta bahwa pendekatan CFD merupakan metode yang andal untuk mendukung desain dan optimasi sistem aerodinamika secara komputasional.

Daftar Pustaka:

• Schlichting, H., & Gersten, K. (2016). Boundary-Layer Theory (9th ed.). Springer.
• White, F. M. (2011). Fluid Mechanics (7th ed.). McGraw-Hill.