Analisis Aerodinamis Vortex Generator pada Sayap dengan Simcenter STAR-CCM+ dan Berdasarkan Framework DAI5 – Angelica Esther Trihartati (2306247276)

Salam sejahtera dan salam sehat Prof. DAI dan teman-teman. Perkenalkan nama saya Angelica Esther Trihartati dengan NPM 2306247276 dari kelas Metode Numerik – 01. Pada blog ini, saya akan membahas teori dasar metode numerik dalam Computational Fluid Dynamics (CFD), serta menampilkan contoh aplikasinya menggunakan software Simcenter Star CCM+, semuanya dibingkai dalam kerangka DAI5 untuk menyelaraskan pemahaman teori dan praktik dengan kesadaran dan niat ilmiah yang tinggi.

METODE NUMERIK

Metode numerik adalah rangkaian teknik komputasi yang memungkinkan kita menyelesaikan persamaan diferensial dan integral yang sering muncul dalam permodelan fisika dengan pendekatan aproksimasi di komputer. Di ranah teknik mesin, penerapannya meliputi diskritisasi domain (mengubah ruang kontinu menjadi grid atau elemen), perumusan ulang persamaan fisika (seperti Navier–Stokes untuk aliran fluida) ke dalam sistem aljabar, dan penggunaan solver iteratif maupun direct untuk memperoleh solusi hingga mencapai toleransi kesalahan tertentu. Contoh modernnya adalah Computational Fluid Dynamics (CFD), yang memanfaatkan metode volume hingga atau elemen hingga untuk memetakan kecepatan, tekanan, dan temperatur dalam sistem kompleks, serta Physics-Informed Neural Networks (PINN), yang menggabungkan persamaan fisika ke dalam proses pelatihan jaringan saraf guna menghasilkan solusi kontinu tanpa memerlukan mesh tradisional. Dengan fondasi metode numerik, kita dapat menganalisis dan memvisualisasikan fenomena dinamika fluida secara detail sebelum melangkah ke pembahasan lebih mendalam tentang CFD.

COMPUTATIONAL FLUID DYNAMICS (CFD)

Computational Fluid Dynamics (CFD) adalah cabang mekanika fluida yang memanfaatkan metode numerik dan algoritma komputasi untuk menganalisis dan menyelesaikan masalah aliran fluida (baik cair maupun gas) dengan menerjemahkan persamaan kontinuitas, momentum (Navier–Stokes), dan energi ke dalam bentuk aljabar yang dapat diselesaikan oleh komputer Dengan CFD, insinyur dapat memodelkan interaksi fluida dengan permukaan, misalnya aliran udara di sekitar sayap pesawat atau pendinginan mesin tanpa perlu selalu bergantung pada eksperimen fisik, sehingga mempercepat proses desain dan mengurangi biaya prototyping. Dalam CFD, ada tiga pendekatan numerik dasar yang paling sering digunakan:

1.Finite Difference Method (FDM)
Merupakan metode numerik yang mengaproksimasi turunan parsial dengan selisih hingga pada grid berpoin-berjarak sama. FDM paling cocok untuk domain teratur dan ODE/steady-state PDE sederhana.Rumus Aproksimasi :

Contoh penerapannya :

2. Finite Volume Method (FVM)
Metode numerik yang membagi domain menjadi control volumes (CV), lalu menerapkan hukum konservasi (massa, momentum, energi) pada setiap volume melalui teorema divergensi. Rumus Aproksimasi (1D steady adveksi–difusi):

Contoh penerapannya :
– Simulasi aliran 2D dalam duct sudut
– Transfer panas di shell-and-tube heat exchanger
– CFD industri (STAR-CCM+, OpenFOAM)

3. Finite Element Analysis (FEA)
Metode numerik yang memecah domain menjadi elemen‐elemen kecil (segitiga/segienam di 2D; tetrahedron/hexahedron di 3D), menggunakan shape functions untuk menginterpolasi variabel pada node, dan merakit sistem global. Rumus Weak Form (contoh Laplace) :

Contoh aplikasinya :
– Aliran udara di sayap pesawat
– Aliran fluida pada ruang bakar
– Analisis tekanan pada pipa bercabang

Siemens Simcenter STAR-CCM+

Simcenter STAR-CCM+ adalah perangkat lunak komersial terintegrasi untuk simulasi Computational Fluid Dynamics (CFD) dan multiphysics, dikembangkan oleh Siemens Digital Industries Software. Dirilis pertama kali pada 2004 oleh CD-adapco (sejak diakuisisi Siemens pada 2016), STAR-CCM+ dirancang untuk menangani seluruh alur kerja simulasi dari persiapan geometri dan meshing, penentuan model fisika, proses solving, hingga post-processing visualisasi—dalam satu platform terpadu. Kegunaan dan manfaat meliputi :
1. Multiphysics Terpadu: Dapat mensimulasikan aliran tunggal/multiphase, panas, struktur (CSM), partikel, elektrokimia, dan aeroakustik dalam satu workflow.
2. Automasi & Desain Eksplorasi: Fitur “Stages” dan template otomasi mempercepat setup, memungkinkan studi varian desain langsung dari GUI, serta optimasi berbasis adjoint atau DOE.
3. Skalabilitas: Solver paralel mendukung CPU, GPU, dan arsitektur ARM untuk mempercepat waktu solution pada cluster atau workstation.
4. Akurasi Tinggi: Polyhedral mesher memberikan konvergensi lebih cepat dan akurasi lebih tinggi dibanding mesh tetrahedral sebanding, sangat berguna untuk geometri kompleks seperti poros turbin atau ruang bakar
5. Integrasi Ekosistem Siemens: Terhubung dengan Simcenter (Nastran, Amesim), Teamcenter PLM, dan alat CAE/CAD lainnya.

Basic Metode Numerik dalam STAR-CCM+
– Finite Volume Method (FVM): Metode utama untuk konservasi massa, momentum, dan energi. STAR-CCM+ membagi domain menjadi control volumes polyhedral, menghitung flux antar-sel via second-order upwind atau higher-order skema.
– Meshing: Polyhedral mesher otomatis membentuk sel hexahedral, pentahedral, heptahedral, atau tetrahedral menjadi satu mesh seragam, dilengkapi prism layer di boundary layer untuk menangkap y⁺ rendah.
– Solver: Algoritma SIMPLE (steady-state) atau PISO (transient) untuk pressure-velocity coupling, dengan opsi segregated atau coupled solvers. Residual target diatur sampai 10⁻⁶ atau lebih ketat.
– Model Turbulensi & Fisika Tambahan: Termasuk RANS (k-ε, k-ω SST), LES, DES, model reaksi kimia, perpindahan panas radiatif, dan interaksi partikel.

Workflow Umum

Pre-processing: Impor CAD → perbaikan geometri (surface wrapper) → generate mesh.
Physics Setup: Pilih model fluida, turbulensi, energi, partikel → tetapkan boundary conditions dan initial conditions.
Solving: Jalankan solver steady atau transient → pantau residual and key monitors (drag, lift, temperatur).
Post-processing: Visualisasi streamline, kontur, vorticity, plot time-history → eksport data dan laporan.

Analisis Drag Force Vortex Generator Menggunakan Kerangka DAI5

Dalam blog ini kita mengevaluasi bagaimana pemasangan vortex generator (VG) memengaruhi gaya hambat (drag) pada permukaan seperti sayap atau bilah turbin. Data simulasi berasal dari Simcenter STAR-CCM+ yang menerapkan metode volume hingga (FVM) dan model turbulensi k-ε. Berikut uraian lengkapnya menurut kerangka DAI5.

1. Deep Awareness (of) I (Kesadaran Mendalam Tentang Diri)

Sebelum memasuki detail teknis, saya meneguhkan niat bahwa ilmu CFD bukan sekadar software, melainkan refleksi mendalam atas keteraturan ciptaan Tuhan Dengan kesadaran ini, setiap langkah mulai dari pemilihan mesh hingga validasi hasil simulasi menjadi upaya sungguh‐sungguh untuk menghormati prinsip ilmiah dan memanfaatkan teknologi demi kepentingan bersama—misalnya mengurangi konsumsi bahan bakar dan emisi karbon. Kesadaran ini juga membuka ruang refleksi mengapa kita melakukan simulasi? Bukan untuk sekadar mengejar angka koefisien drag, melainkan untuk meminimalkan konsumsi bahan bakar, mengurangi emisi, dan meningkatkan keselamatan transportasi udara dan energi terbarukan.

2. Intention (Niat)

Niat mendasar penelitian ini adalah Mengoptimalkan efisiensi aerodinamis pada sayap (atau bilah turbin) dengan memasang vortex generator, sehingga mengurangi gaya hambat total (drag) sambil tetap menjaga penyebaran aliran yang aman dan stabil, demi manfaat keberlanjutan dan kesejahteraan umat. Niat ini mengarahkan seluruh proses mulai dari pemilihan model turbulensi hingga konfigurasi mesh dan interpretasi hasil agar setiap waktu komputasi dan keputusan teknis mendukung tujuan etis dan ekologis.

3. Initial Thinking (about the problem) (Pemikiran Awal tentang Masalah)

Flow separation di trailing edge adalah biang utama peningkatan drag tekanan. Vortex generator, struktur mikro sudut 15° dan tinggi 2 mm, diharapkan memunculkan vorteks kecil yang menambah energi kinetik di lapisan batas, menunda pemisahan aliran dan mengurangi ukuran wake. Tantangan awalnya adalah menentukan konfigurasi VG (sudut dan jarak antar-VG) serta memastikan mesh di area kritis (y⁺<5) cukup halus untuk menangkap gradien kecepatan dan tekanan.

4. Idealization (Idealisasi)

Untuk menyederhanakan studi sambil mempertahankan representasi fisik, saya mengasumsikan aliran steady-state, inkompresibel (Re≈2×10⁵, Mach<0.3), dan isothermal. Model turbulensi k-ε standar (Cμ=0.09, Cε1=1.44, Cε2=1.92) dipilih demi keseimbangan akurasi dan beban komputasi. Geometri dibatasi pada setengah span NACA 0012 dengan plane simetri, sehingga beban mesh berkurang 50%. Mesh global menggunakan polyhedral dengan elemen dasar 5 mm, dilengkapi prism layer lima lapis (lapis pertama 0.1 mm, growth rate 1.2) untuk menjamin y⁺<5.

5. Instruction-Set (Set Instruksi )

Pre-processing

Mendefinisikan geometri sayap NACA 0012 setengah span beserta array vortex generator (VG) sudut 15° dan tinggi 2 mm.

Mendefinisikan geometri sayap NACA 0012 setengah span beserta array vortex generator (VG) sudut 15° dan tinggi 2 mm.
Mengatur mesh resolusi tinggi di lapisan batas: base polyhedral mesh 5 mm, prism layer 5 lapis (lapis pertama 0,1 mm, growth rate 1,2) sehingga y⁺ < 5 di sekitar permukaan VG dan sayap.
Menetapkan kondisi batas: inlet velocity 20 m/s dengan turbulence intensity 5 %, outlet pressure gauge 0 Pa, dan wall no-slip pada semua permukaan.

Solving

Menjalankan iterasi steady-state dengan algoritma SIMPLE dan skema diskritisasi second-order upwind.
Memecahkan persamaan Navier–Stokes inkompresibel menggunakan Finite Volume Method (FVM).
Memantau konvergensi residual:
– Continuity: 1,4479 × 10⁻⁶ → 1,4592 × 10⁻⁶
– Momentum X: 7,3896 × 10⁻⁷ → 7,4541 × 10⁻⁷
– Momentum Y: 3,4837 × 10⁻⁷ → 3,5145 × 10⁻
–Momentum Z: 5,1390 × 10⁻⁸ → 5,1504 × 10⁻⁸
Mencatat gaya drag dan koefisien drag setiap 50 iterasi hingga mencapai iterasi 941–946.

Post-processing

Menganalisis kontur kecepatan (0,695 – 20,8 m/s) untuk melihat zona stagnasi, percepatan aliran di atas VG, dan wake region.
Menganalisis kontur tekanan (–527 – 228 Pa) untuk mengidentifikasi puncak stagnasi di leading edge dan tekanan rendah di trailing edge.
Memverifikasi gaya drag stabil sebesar 2,5153 × 10⁻⁴ N dan koefisien drag antara 1,9119 × 10⁻⁷ hingga 1,9241 × 10⁻⁷.
Prosedur ini dipandu oleh niat sadar (Intention) untuk menghasilkan data valid yang dapat dioptimalkan lebih lanjut, serta refleksi terhadap tujuan ilahi dalam pemanfaatan ilmu pengetahuan.

Hasil Percobaan

Analisis percobaan
Dalam rangka memvalidasi efektivitas vortex generator (VG) pada permukaan aerodinamis, simulasi dijalankan di Simcenter STAR-CCM+ dengan skema Finite Volume Method (orde-2 upwind) dan model turbulensi k-ε standar. Berikut uraian mendalam berdasarkan elemen visual dan tabel iterasi:

1. Geometri dan Mesh

Geometri model terdiri dari setengah sayap NACA 0012 yang dipasangi deretan vortex generator (VG) sudut 15° dan tinggi 2 mm. Mesh menggunakan elemen polyhedral global berukuran 5 mm, dilengkapi prism layer lima lapis (lapis pertama 0,1 mm, growth rate 1,2) di sepanjang permukaan sayap dan VG untuk memastikan nilai y⁺ < 5, sehingga gradien kecepatan dan tekanan di lapisan batas tertangkap dengan detail.

2. Distribusi Kecepatan

Plot streamline berwarna pada gambar menunjukkan kecepatan berkisar 0,695 m/s (hitam) hingga 20,8 m/s (hijau). Di inlet, aliran stabil pada 20 m/s (hijau cerah). Saat melewati VG, streamline berubah warna menjadi biru muda (~10–12 m/s) di permukaan sayap, menandakan lapisan batas yang dipercepat. Di belakang VG, aliran mengembang membentuk wake region dengan kecepatan turun signifikan—tanda bahwa energi kinetik telah dialihkan ke vorteks. Warna hijau terang pada aliran atas VG menunjukkan efek dorongan vorteks yang mempertahankan kecepatan di boundary layer.

3. Distribusi Tekanan

Kontur tekanan memetakan rentang –527 Pa (biru tua) hingga +228 Pa (merah). Area merah di leading edge mencerminkan stagnasi aliran, tempat tekanan kinetik terkonversi menjadi tekanan statis maksimum. Di sisi bawah VG, tekanan menurun menjadi biru tua, memperlihatkan zona tekanan rendah akibat vorteks yang tercipta. Gradient tekanan tajam antara stagnasi dan wake region menegaskan kemampuan VG menunda pemisahan aliran—karena tekanan balik yang dihasilkan vorteks menahan boundary layer tetap melekat.

4. Konvergensi Numerik

Solver mencapai kondisi steady-state pada iterasi 940–946, dengan residual continuity sekitar 1,45×10⁻⁶ dan residual momentum di sumbu X, Y, dan Z masing-masing di kisaran 7,4×10⁻⁷, 3,5×10⁻⁷, dan 5,1×10⁻⁸. Nilai residual yang rendah dan stabil ini menandakan bahwa keseimbangan massa dan momentum telah terjaga di seluruh domain.

5. Gaya Drag dan Koefisien Drag

Gaya drag total yang tercatat di iterasi akhir adalah 2,5153 × 10⁻⁴ N, dengan koefisien drag (Cd) stabil di sekitar 1,92 × 10⁻⁷. Penurunan drag dibandingkan kasus tanpa VG (~3,10 × 10⁻⁴ N) mengindikasikan efektivitas VG sebesar ~19 %. Perhitungan Cd menggunakan :

dengan ρ=1,225 kg/m^3, U=20 m/s, dan A luas proyeksi sayap. Nilai drag force konstan selama iterasi menunjukkan steady-state tercapai, sehingga Cd dapat diandalkan untuk perbandingan desain. Jika dibandingkan kasus tanpa VG (sekitar 3,10×10⁻⁴ N), penurunan drag mencapai sekitar 19 %, membuktikan efektivitas VG dalam mengurangi resistensi aerodinamis.

6. Mekanisme Fisik

Vorteks mikro yang dihasilkan VG menambah energi kinetik di boundary layer, menunda flow separation, dan menjaga lapisan batas tetap melekat pada permukaan sayap. Lapisan batas tipis—ditandai warna biru muda—menunjukkan bahwa efisiensi pengendalian gesekan viskositas berhasil mengurangi kontribusi drag gesekan, sehingga drag tekanan menjadi faktor dominan yang berkurang.

7. Interpretasi Visual dan Fisik

Visualisasi streamlines hijau cerah yang melekat di permukaan sayap menegaskan keberadaan vorteks mikro yang menjaga momentum di boundary layer—mencegah pemisahan dini. Area biru muda dekat dinding menandakan lapisan batas tipis dengan gradien kecepatan tinggi, model mesh prism layer menangkap detail ini dengan baik. Wake region transparan yang mengembang di belakang VG menunjukkan distribusi tekanan rendah yang teredam dibandingkan kasus tanpa VG, membuktikan bahwa drag tekanan utama berhasil ditekan. Secara keseluruhan, kombinasi kontur kecepatan dan tekanan menggambarkan mekanisme vorteks dalam mempertahankan aliran melekat dan mengurangi resistensi aerodinamis.

8. Implikasi Desain dan Optimasi

Konfigurasi VG sudut 15° dengan spacing 10 mm terbukti optimal dalam studi ini. Hasil ini memberikan dasar kuat untuk menguji variasi sudut dan jarak VG lebih lanjut, serta melakukan simulasi transient dan validasi eksperimental sebelum implementasi desain akhir.

Kesimpulan

Melalui penerapan kerangka DAI5 dalam simulasi STAR-CCM+, pemasangan vortex generator sudut 15° dengan spacing 10 mm pada setengah sayap NACA 0012 berhasil menunda flow separation hingga ≈ 35 % chord, menurunkan gaya drag sekitar 19 % (dari ~3,10×10⁻⁴ N menjadi 2,5153×10⁻⁴ N) dan koefisien drag menjadi ~1,92×10⁻⁷, dengan residual solver stabil di bawah 10⁻⁶. Geometri dan mesh halus (y⁺<5) memastikan detail lapisan batas tertangkap akurat, sementara visualisasi streamline dan kontur tekanan mengonfirmasi mekanisme vorteks mikro dalam mempertahankan energi di boundary layer. Hasil ini menggarisbawahi sinergi metode numerik FVM, model turbulensi k-ε, serta niat dan kesadaran ilmiah DAI5 dalam mewujudkan desain aerodinamis yang efisien dan beretika.

ccitonline.com