1001Ferramentas
🔥 Calculators

Karlovitz Number

Computes the Karlovitz number of a turbulent premixed flame, Ka = (u'/S_L)^1.5 · (ℓ_t/δ_L)^−0.5, comparing the chemical time of the flame front with the turnover time of the smallest eddies (Kolmogorov scale). Ka below 1 means the laminar flame structure survives the turbulence (wrinkled and corrugated regimes); between 1 and 100 eddies penetrate the preheat zone and thicken the flame; above 100 the reaction zone itself is broken. Together with the Damköhler number it forms the axes of the Borghi diagram, used to pick combustion models in CFD. Peters' form is adopted; integral scale and flame thickness must share the same unit. Enter the velocity fluctuation, the laminar flame speed, the integral length scale and the laminar flame thickness.

Result

Karlovitz Number: locating a flame on the Borghi diagram

Before picking a combustion model for a simulation, someone has to decide whether the flame still behaves as a thin sheet wrinkled by the flow or whether turbulence has already worked its way inside it. Getting that wrong is expensive: a flamelet model used in a regime where eddies have torn into the reaction zone predicts a stability the burner does not have, under-predicts CO and misses blow-off. The Karlovitz number is the vertical axis of the Borghi diagram, and it settles the question with a single figure.

The tool uses Peters' form, Ka = (u'/S_L)^1.5 · (ℓ_t/δ_L)^−0.5, picked because it needs only four quantities that any measurement or LES already provides. Here u' is the turbulent velocity fluctuation, roughly 0.5 to 10 m/s in industrial burners; S_L is the laminar flame speed, 0.4 m/s for stoichiometric methane-air and above 2 m/s for hydrogen; ℓ_t is the integral length scale, typically the order of the nozzle diameter; δ_L is the laminar flame thickness, near 0.5 mm at 1 atm. Below Ka = 1 the laminar structure survives; from 1 to 100 the smallest eddies thicken the preheat zone; above 100 the reaction layer itself breaks up.

The expression assumes homogeneous isotropic turbulence with a fully developed Kolmogorov cascade. At low turbulent Reynolds number, close to a wall or inside a recirculation bubble, the value it returns is a rough estimate at best. It carries no information about Lewis number, curvature-induced stretch or heat loss, which happen to be exactly the effects that govern extinction around Ka = 1. Units deserve the same attention: ℓ_t and δ_L must share one unit, since mixing millimetres with metres divides Ka by 31.6, the square root of a thousand.

Frequently asked questions

What does the Ka = 3.536 shown on page load mean?
With u' = 2 m/s, S_L = 0.4 m/s, a 5 mm integral scale and a 0.5 mm flame thickness, the value 3.536 sits between 1 and 100, in the thin reaction zone regime. Kolmogorov eddies are already smaller than the flame thickness and penetrate the preheat layer, thickening it, while the reaction layer itself — about a tenth of δ_L — stays intact. Flamelet models still work here if a thickening correction is applied; below Ka = 1 they would need no such caveat.
Can I enter the integral scale in millimetres and the thickness in metres?
No. Nothing gets converted: the two numbers are simply divided, and the ratio has to come out dimensionless. Typing 5 for the scale while thinking in millimetres and 0.0005 for the thickness while thinking in metres turns the ratio into 10,000 instead of 10, dropping Ka from 3.536 to 0.112 — enough to file a thin reaction zone flame as a corrugated one. The same rule applies to u' and S_L, which both belong in metres per second.
Which entries does the calculator reject, and what shows up then?
Flame speed, integral scale and flame thickness must all be strictly positive; a zero or negative value in any of them raises the message 'Check the values you entered.' and the result reverts to a dash. A fluctuation of u' = 0 is accepted on purpose and returns Ka = 0, the laminar limit. Clearing a field counts as incomplete rather than wrong: the dash comes back until the box is filled again, and the number refreshes on every keystroke, with no need to press Calculate.

Related Tools

🔥

Zeldovich Number

Computes the Zeldovich number of a flame, β = E_a·(T_b − T_u) ÷ (R·T_b²), the activation energy made dimensionless by the temperature rise across the flame front. It measures how sensitive the reaction rate is to a small temperature change: a high β (typically 8 to 12 for hydrocarbons) means the reaction is concentrated in a very thin layer near the flame temperature, which justifies the large-activation-energy assumption of asymptotic flame theory and the extinction and cellular-instability criteria. The universal gas constant R = 8.314 J/(mol·K) is adopted, with activation energy in J/mol and temperatures in kelvin. Enter the activation energy, the burned gas temperature and the unburned gas temperature.

🌡️

Adiabatic Flame Temperature

Estimate the adiabatic flame temperature, T_ad = T_initial + LHV ÷ (m × cp), from the lower heating value LHV (kJ/kg fuel), the mass of combustion products per kg fuel m, the average specific heat of the gases cp (kJ/kg·K) and the initial temperature. The result, in °C, is the maximum theoretical temperature the gases would reach if all the combustion energy heated the products, with no heat loss. It is an upper bound: real flames are cooler (radiation losses, dissociation, excess air). It sets the thermal severity on materials and NOx formation. Enter the LHV, the gas mass, the cp and the initial temperature.

🌡️

Elenbaas Number

Computes the Elenbaas number of a channel formed by two heated vertical parallel plates, El = g·β·ΔT·b⁴ ÷ (ν·α·H), which is the Rayleigh number based on the plate spacing b multiplied by the ratio b/H. It governs natural convection in finned heat sinks, electronics enclosures and solar collectors: a very small El means a narrow, tall channel where the boundary layers merge and choke the flow, while a large El means the plates behave as isolated. The spacing that maximises a heat sink's total dissipation falls near El ≈ 46 for isothermal plates, a classic natural-convection fin design criterion. Gravity is taken as 9.80665 m/s², β is the fluid thermal expansion coefficient and H the plate height. Enter the expansion coefficient, the temperature difference, the spacing, the kinematic viscosity, the thermal diffusivity and the plate height.

🔇

Acoustic Barrier Attenuation (Maekawa)

Computes how much an acoustic barrier cuts the noise reaching a receiver, using Maekawa's empirical formula, Attenuation = 10 × log₁₀(3 + 20N), where N is the Fresnel number, equal to twice the path difference divided by the wavelength — that is, N = 2 × path difference × frequency ÷ speed of sound. The path difference is the extra distance sound must travel to go over the top of the barrier instead of straight from source to receiver, and it is the only geometric input the formula needs. The result, in decibels, is what the barrier subtracts from the level that would arrive without it: with N equal to zero, meaning the receiver sits exactly on the line of sight to the top edge, attenuation is already 4.8 dB, and in practice the gain saturates between 20 and 24 dB because sound eventually flanks around the sides and passes through the panel. Since N grows with frequency, the same barrier is far more effective at high frequencies than at low ones: once N is large, every octave adds about 3 dB, which is why enclosing a compressor kills the hiss and barely touches the rumble. Enter the path difference, the frequency and the speed of sound.

Flash Guide Number by Distance

Returns ideal aperture from flash Guide Number and subject distance.

🌀

Iribarren Number

Calculate the Iribarren number (surf similarity parameter), ξ = tan(β) ÷ √(H ÷ L), from the beach slope tan(β), the wave height H and the wavelength L. The dimensionless result classifies the wave breaking type and runup: ξ < 0.5 indicates spilling breakers (flat beaches); 0.5 < ξ < 3.3, plunging (the wave 'barrels'); ξ > 3.3, surging/collapsing (steep beaches). It is widely used in coastal engineering to predict runup, structure overtopping and riprap stability. Enter the beach slope, wave height and wavelength.

The results provided by this tool are for general informational and educational purposes only and do not constitute professional, financial, medical, legal, tax or accounting advice. Always confirm important decisions with a qualified professional and official sources.