Assalamualaikum Wr. Wb.
Puji syukur penulis panjatkan ke hadirat Tuhan Yang Maha Esa atas limpahan rahmat dan karunia-Nya sehingga penulis dapat menyelesaikan laporan tugas bertajuk “Numerical Aerodynamics of 3D Airfoil: A Finite Element-Based CFD Simulation Using STAR-CCM+” dengan baik dan lancar. Tugas ini disusun sebagai bagian dari pemenuhan mata kuliah Metode Numerik, yang secara khusus menekankan pada pemodelan numerik dan pendekatan metode elemen hingga (Finite Element Method/FEM) dalam konteks rekayasa teknik mesin.
Penulis juga menyampaikan apresiasi yang sebesar-besarnya kepada Prof. DAI, yang telah memberikan kesempatan berharga untuk mendalami kembali topik Computational Fluid Dynamics (CFD) secara praktikal menggunakan STAR-CCM+. Melalui tugas ini, penulis tidak hanya mampu mengaplikasikan prinsip-prinsip numerik dalam pemecahan masalah aliran fluida, tetapi juga memperoleh pemahaman menyeluruh tentang bagaimana pendekatan FEM diimplementasikan secara nyata dalam simulasi aerodinamika berbasis perangkat lunak industri.
1. Pendahuluan
Perkembangan teknologi simulasi numerik memungkinkan analisis performa aerodinamika dilakukan secara digital menggunakan pendekatan Computational Fluid Dynamics (CFD). Salah satu kasus penting dalam CFD adalah simulasi aliran udara di sekitar airfoil tiga dimensi (3D), yang secara langsung memengaruhi performa sayap pesawat terbang.
Melalui perangkat lunak Simcenter STAR-CCM+, studi ini menyimulasikan distribusi tekanan dan kecepatan di sekitar airfoil, menghitung gaya angkat dan hambat, serta menganalisis konvergensi solusi numerik menggunakan pendekatan Finite Element Method (FEM). Proses ini juga dievaluasi melalui lensa 5i dari DAI5 (Deep Awareness, Intention, Initial Thinking, Idealization, Instruction Set) untuk merefleksikan kesiapan digitalisasi industri.
2. Tujuan
- Melakukan simulasi aliran udara pada sayap 3D menggunakan STAR-CCM+.
- Menganalisis distribusi tekanan dan kecepatan.
- Menghitung koefisien tekanan serta mendeteksi fenomena aerodinamika penting seperti lift dan drag.
- Menunjukkan penerapan Finite Element Method dalam proses numerik CFD.
- Mengkaitkan proses pengerjaan dengan analisis DAI5 sebagai cermin penerapan digitalisasi rekayasa.
3. Landasan Teori
Simulasi Computational Fluid Dynamics (CFD) merupakan pendekatan numerik untuk menyelesaikan dan menganalisis permasalahan aliran fluida dengan bantuan komputasi. Di balik visualisasi aliran dan kontur tekanan yang dihasilkan, terdapat struktur matematis dan numerik yang kompleks, yaitu sistem persamaan diferensial parsial (PDE) yang merepresentasikan hukum fisika fundamental: konservasi massa, momentum, dan energi.
Dalam studi aerodinamika airfoil ini, software Simcenter STAR-CCM+ digunakan untuk mensimulasikan perilaku aliran udara di sekitar airfoil 3D dengan menyelesaikan persamaan Navier-Stokes sebagai model utama. Persamaan ini tidak memiliki solusi eksak untuk kebanyakan kasus praktikal, sehingga pendekatan numerik menjadi satu-satunya jalan.
1. Persamaan Governing
a. Konservasi Massa (Continuity Equation):
b. Konservasi Momentum (Navier-Stokes):
c. Konservasi Energi:
Di mana:
- v : vektor kecepatan fluida
- p : tekanan
- ฯ : densitas
- ฮผ : viskositas dinamik
- h : entalpi
- T : suhu
- ฮฆ : dissipasi energi mekanik
2. Diskretisasi Domain: Finite Element Method (FEM)
Dalam STAR-CCM+, pendekatan Finite Element diterapkan secara implisit melalui meshing domain dan penyelesaian variabel aliran secara lokal pada setiap elemen.
a. Pembentukan Mesh
Domain kontinu (ฮฉ) dipecah menjadi elemen diskrit (ฮฉeโ):
Fungsi-fungsi basis lokal (Ni(x)) digunakan untuk mendekati distribusi kecepatan, tekanan, dan energi dalam setiap elemen:
b. Weak Form dan Galerkin Method
Untuk memastikan kestabilan numerik, persamaan diferensial diubah ke bentuk integral (weak form), lalu diaproksimasi dengan metode Galerkin:
Dengan R(u) sebagai residual persamaan dan www sebagai fungsi uji.
3. Penyelesaian Sistem Persamaan
Sistem non-linear besar hasil diskretisasi diselesaikan dalam bentuk:
[K]{u}={f}
- [K]: matriks kekakuan
- {u}: vektor solusi
- {f}: gaya luar, termasuk efek batas
Solver iteratif (segregated atau coupled) menyelesaikan sistem ini hingga konvergen berdasarkan nilai residual.
4. Model Turbulensi dan Transfer Energi
Dalam kasus airfoil, model turbulensi seperti k-ฮต atau k-ฯ SST digunakan, sehingga perlu menyelesaikan dua persamaan tambahan:
- Turbulent kinetic energy (k):
- Turbulent dissipation rate (ฮต atau ฯ) dengan bentuk serupa.
Hasil dari model ini berdampak langsung pada profil tekanan dan kecepatan, serta prediksi gaya angkat dan hambatan.
4. Metodologi Simulasi
3.1 Tahap Pre-processing (Pembuatan Domain dan Meshing)
Gambar 1: Diskretisasi Mesh
- Geometri 3D airfoil dimodelkan dalam domain fluida.
- Boundary condition: Inlet (velocity inlet), outlet (pressure outlet), dinding no-slip (airfoil), dan simetri.
- Mesh polyhedral digunakan dengan refinement lokal di sekitar permukaan airfoil.
- Jenis elemen mesh mendekati elemen hexahedral terstruktur untuk efisiensi kalkulasi.
FEM Relevansi:
Setiap sel dalam mesh mewakili elemen volume hingga FEM menyelesaikan sistem persamaan diferensial parsial berbasis Galerkin Weak Form:
Dengan ฯi sebagai fungsi basis lokal di setiap elemen.
3.2 Solving (Proses Iteratif)
Gambar 2: Plot Residuals
Solver menyelesaikan persamaan konservasi:
- Massa:
- Momentum:
- Energi: (jika compressible)
- Model turbulensi: kโฮต (Tke & Tdr).
- Konvergensi dicapai saat residual <10^-5.
FEM Relevansi:
Solusi dihitung melalui sistem linier diskret:
[K]{u}={f}
dengan [K] adalah stiffness matrix yang dibentuk dari integral over elemen domain.
3.3 Post-processing
Gambar 3: Grafik Koefisien Tekanan (Cp)
- Distribusi tekanan pada permukaan airfoil menunjukkan efek Bernoulli: tekanan rendah di permukaan atas menghasilkan lift.
- Rumus:
Gambar 4: Kontur Tekanan
- Perbedaan tekanan jelas terlihat: warna biru di atas (low pressure), merah di bawah (high pressure).
- Leading edge menunjukkan titik stagnasi (tekanan maksimum).
Gambar 5: Kontur Kecepatan
- Kecepatan tinggi di atas airfoil (area merah), kecepatan rendah di bawah dan wake region di belakang.
FEM Relevansi (Output Pasca-Proses):
Distribusi nilai-nilai hasil simulasi (velocity, pressure) adalah hasil interpolasi solusi {u} dari sistem FEM di seluruh domain.
5. Analisis DAI5 dalam Proyek
4.1 Deep Awareness
Pemahaman mendalam tentang fenomena fisik (aliran udara, tekanan, turbulensi) menjadi landasan pemodelan. Kesadaran bahwa simulasi CFD mewakili sistem fisik nyata mendorong kehati-hatian dalam diskretisasi dan pemilihan model turbulensi.
4.2 Intention
Tujuan dari tugas ini bukan hanya menyelesaikan simulasi, tetapi memahami proses matematis dan numerik (FEM) di baliknya. Intensi utama adalah menganalisis perilaku aerodinamika airfoil melalui pendekatan numerik.
4.3 Initial Thinking
Perencanaan awal mencakup:
- Pemilihan geometri sederhana namun representatif.
- Asumsi steady flow, inkompresibel, dan turbulent.
- Strategi pemilihan mesh refinement di zona gradien tinggi.
4.4 Idealization
Model disederhanakan dari geometri nyata pesawat menjadi airfoil 3D terisolasi dengan asumsi simetri dan boundary yang ideal. Ini bertujuan menyeimbangkan antara kompleksitas fisik dan keterbatasan komputasi.
4.5 Instruction Set
Software STAR-CCM+ digunakan berdasarkan SOP numerik: setup physics model, boundary conditions, meshing, solving, hingga post-processing. Proses ini sepenuhnya terdigitalisasi dan menjadi refleksi dari pengotomatisasian proses rekayasa digital.
6. Kesimpulan
Simulasi CFD pada airfoil 3D menggunakan STAR-CCM+ berhasil menunjukkan fenomena aerodinamika utama: lift akibat tekanan rendah di atas airfoil dan wake region akibat viskositas fluida.
Dari perspektif metode numerik, seluruh proses mensyaratkan pemahaman mendalam tentang Finite Element Method sebagai dasar penyelesaian sistem PDE dalam domain terdiskretisasi.
Penerapan DAI5 menunjukkan bahwa kegiatan ini tidak sekadar tugas teknis, melainkan bagian dari transformasi rekayasa menuju digitalisasi cerdas dan sadar dalam era Industri 5.0.
Terima kasih atas perhatiannya,
Wassalamualaikum Wr. Wb.
