Numerical Modelling of Combustion: By: Djatmiko Bagus Maulana. E
Numerical Modelling of Combustion: By: Djatmiko Bagus Maulana. E
Numerical Modelling of Combustion: By: Djatmiko Bagus Maulana. E
COMBUSTION
By : Djatmiko Bagus Maulana. E
Reaksi Pembakaran
Pembakaran Ideal 𝐶𝑥 𝐻𝑦 + 𝑂2 → 𝐶𝑂2 + 𝐻2 𝑂
GEOMETRY RESULT
INPUT DATA
CFD RESEARCH METHOD
1. PROBLEM 2. FLOW REGIMES 3. VARIABLE
IDENTIFICATION IDENTIFICATIONS LISTING
5. PRE PROCESSING
-Geometry
6. CALCULATION 4. GOVERNING
-Meshing
METHOD EQUATIONS
- Boundary Input
-Material Properties
8. POST 9. NUMERICAL
7. CALCULATION
PROCESSING VALIDATION
1. PROBLEM IDENTIFICATION
“NON-PREMIXED
COMBUSTION”
2. FLOW REGIMES IDENTIFICATION
Flow
Chemical Conduction
Equilibrium Convection
Chemical Heat
Reaction Combustion Transfer
Kinematics Radiation
Turbulent
Mixing
3. VARIABLE LISTING
HEAT TURBULENCE 2 CHEMICAL
FLOW 4 TRANSFER 1 REACTION 2
• Tekanan (P) • Suhu (T) • Turbulence Kinetic • Mixture Fraction
• Kecepatan x (Vx) Energy (K) (f)
• Kecepatan y (Vy) • Turbulence • Mixture Fraction
• Kecepatan z (Vz) Dissipasion Energy Variance (f2)
(Ɛ)
In Out
• Transient
• Convection
• Diffusion
Gen • Source
4. GOVERNING EQUATIONS
”Governing equation merupakan persamaan yang mengatur hukum-hukum konservasi di alam semesta”
KekekalanMassa
Kekekalan Momentum(x,y,z)
KekekalanEnergi
KekekalanTKE
KekekalanTDE
KekekalanMixtureFraction
Kekekalan Mixture
FractionVariance
4. GOVERNING EQUATIONS
How toconfigure governing equation in fluent interface? (1)
General GoverningEquation
Steady =0
Transient
Absolute =v
Relative =v +v0
Pressure-Based 𝝆 = 𝝆 (𝒇)
Densitybased-Based 𝝆 = 𝝆 (𝑴𝒂)
4. GOVERNING EQUATIONS
How toconfigure governing equation in fluent interface? (2)
Kekekalan
TKEdanTDE
4. GOVERNING EQUATIONS
TURBULEN KINETIK ENERGI VS TURBULEN DISIPASI ENERGI
TDE
TKE TKE
TKE
Kekekalan
Mixture
Fraction &
Variance
NON-PREMIXED COMBUSTION
Species
• Menghitung seluruh
Transport mass fraction species
• Membutuhkan waktu
perhitungan yangcukup
lama
Non-
Premixed
Combustion
• Hanya menghitung
mixture fraction(f)
Mixture
Fraction • Waktuperhitungannya
singkat
MENGAPA MIXTURE FRACTION ?
Mixture Fraction memiliki relasi linear yang
Temperature
(T)
unik terhadap berbagai variable pada
spesies.(ANSYS, 2017)
Mixture
Air Mole
Fraction(Xa) Fraction Mean
Density (𝜌)
(f)
Fuel Mole
Fraction
(Xf) LeiWang, 2007
SKEMA INLET MIXTURE FRACTION
5. PRE PROCESSING
gasFlowRate (kg/s) 0.505 [Wall]
airFlowRate (kg/s) 27.0485
7 OutletInner
gasInletTemperature (K) 293 [Pressure
airInletTemperature (K) 635
Outlet]
[Periodic]
Inlet-Fuel
Inlet-Air [Massflow inlet]
[Massflow inlet] f =1 CH4 90.45%
C2H6 2.12%
f =0 C3H8 2.66%
N2 78.99% C4H10 1.22%
O2 21.01% C5H12 0.39%
C6H14 0.26%
CO2 1.96%
N2 0.95%
6. CALCULATION
METHOD
Semi-Implicit
Metode perhitungan
Method for SIMPLE Algorithm
Pressure menggunakannilai P
(Pressure) sebagai
Linked factor koreksi untuk
Equations nilai-nilai lain.
6. CALCULATION METHOD
Under
Relaxation
URFtidakmemengaruhi
Factor tingkatakurasi hasil.
Pemilihan URF
berdasarkan TRIAL &
ERROR, nilainya selalu
berbedauntukberbagai
kasus.
Under Relaxation Factor (URF) merupakan konstanta yang berfungsi untuk menentukan nilai
prediksi pada iterasi berikutnya. Nilai URF diubah-ubah agar mencapai kestabilan nilai residual
dan konsistensinilai.
8. POST PROCESSING
Polyhedral Mesh
8. POST PROCESSING
DistribusiTemperatursepanjangsumbuZ DistribusiTemperatursepanjangsumbuY
8. POST PROCESSING
Temperaturkeluar CC:
T (Max) =2411.6062 K
Distribusi T (Mean)=1781.9636 K
Temperatur
T(Min)=1046.7605K
Streamline
8. POST PROCESSING
DistribusiCH4
DistribusiCO
8. POST PROCESSING
DistribusiCO2
DistribusiO2
8. POST PROCESSING
EmisiOutlet
DistribusiH2O
9. NUMERICAL VALIDATION
Walaupun residual
sudah di bawah kriteria,
perlu dipastikan juga
salah satu nilai
tercapture (misal suhu
outlet) telah konsisten.
Jika nilai telah konsisten,
maka kalkulasi dapat
diberhentikan.
Massflow Difference
(A+B)-(C+D)=0