EAF–LF Off-Gas Duct Elbow

A DN1600 elbow in a steel plant’s off-gas duct had to be checked for pressure loss and strength: it loses about half the pressure the handbook predicts, and its safety factor against buckling is about 35, against a target of 4. For Dejpod Sanaat Sazeh: CAD, mesh-verified CFD, and structural FEA loaded by the CFD’s own pressure field.

Project
Mechanical design & simulation internship
Company
Dejpod Sanaat Sazeh (DEJPOD), Isfahan
My role
Design basis, CAD, CFD and structural FEA
Tools
SolidWorks, ANSYS Fluent, ANSYS Mechanical
Method
SST k-ω RANS CFD + shell FEA
3D render of the DN1600 duct elbow: a straight leg, a 90-degree long-radius bend and a second straight leg, with a flange at each end
Duct elbow, my CAD model Interactive 3D view
Interactive render of the duct elbow solid geometry, S235JR carbon structural steel (EN 10025-2).
  • 35.4PaCFD pressure drop, refined mesh: 85% of the corrected benchmark
  • 0.153corrected elbow loss coefficient K, about half the generic 0.3
  • ≈34.7buckling safety factor, the governing check (target 4)
  • ≈540yield safety factor, both load cases

Skills shown

  • CAD Modeling & Shell Idealization (SolidWorks)
  • CFD Meshing & Solving (ANSYS Fluent)
  • Structural FEA (ANSYS Mechanical)
  • Turbulence-Model Selection (Decision Matrix)
  • Mesh Convergence & Verification
  • CFD-to-FEA Load Transfer
  • Pressure-Vessel / Buckling Design (EN 13445-3)
  • Engineering Judgment & Risk Flagging
  • Technical Reporting & Traceability

How I did it

Dejpod Sanaat Sazeh designs and delivers dust-collection and fume-treatment systems for the steel, iron-making and mining sectors. This internship project set out to build and prove a combined CFD/FEA workflow for the company’s duct systems, using one representative EAF-LF off-gas duct section as the case.

  1. Define the problem and choose the analysis approach explicitly, rather than by default assumption.
  2. Fix a complete, traceable design basis for the duct section, anchored to a real reference plant.
  3. Run a mesh-independence-checked CFD analysis of the section’s aerodynamics.
  4. Run a structural (FEA) analysis of the elbow under both an idealized uniform pressure and the CFD’s own resolved wall-pressure field.
  5. Check every stage against independent hand calculation, and close the study out against its own starting specification.
  1. The duct section

    Design basis traced to a reference plant

    The section models one branch of a primary EAF-LF dedusting line: a straight inlet run, a 90-degree long-radius elbow, and a straight outlet run. Flow rate and duct size were derived from a real reference plant rather than assumed, and every value in the design-basis table (technical details) is traced back to that reference or to a stated engineering assumption.

    Internal diameter (DN1600)
    1.60 m
    Design velocity
    21.1 m/s, Re ≈ 9.7×105
    Design gas temperature
    200 °C
    System design (suction) pressure
    3.0 kPa
  2. Method: turbulence-model selection

    Weighted decision matrix + sensitivity check

    The elbow introduces streamline curvature and a risk of local flow separation on the inner wall. Rather than default to one closure, two candidate RANS turbulence models were compared with a weighted decision matrix, then stress-tested with a sensitivity check on the weights themselves.

    Weighted score = Σ (criterion weight × score)each criterion scored 1–3

    • SST k-ω selected: better suited to curved, separation-prone flow, and its near-wall treatment adapts across a wide range of mesh resolution instead of needing a narrow y+ window.
    • A sensitivity check on the weights confirmed the result is stable. SST k-ω only loses out if raw computational economy is made the dominant priority, which it is not for this study.
    • The decision was carried forward without revisiting through the CFD build, run, and verification stages.
    • A risk register scored seven technical uncertainties by probability × severity before any modeling began, each with a stated mitigation carried into the analysis.
    Turbulence-model decision matrix
    CriterionWeightRealizable k-εSST k-ω
    Accuracy for curvature / separation0.3023
    Computational / mesh cost0.2032
    Near-wall mesh-targeting robustness0.3023
    General solver robustness0.2033
    Weighted total1.002.402.80
  3. CFD model: mesh and setup

    ANSYS Fluent · hexahedral O-grid, y+ ≈ 1

    The fluid volume, the space the gas occupies from inlet to outlet, was meshed as a hexahedral O-grid with a boundary-layer inflation targeting a wall y+ near 1, so the near-wall region is resolved directly rather than through a wall function.

    Base and refined mesh statistics
    ControlBase meshRefined mesh
    Elements≈ 1.04 million≈ 4.64 million
    Wall y+ (average)0.941.05
    Orthogonal quality (minimum)0.570.83
    Skewness (maximum)0.440.44 (unchanged)
    Base hexahedral surface mesh of the duct fluid volume, full domain
    Base hexahedral surface mesh, full domain. Select to enlarge.
    Close-up of the hexahedral surface mesh on the elbow
    Close-up of the surface mesh on the elbow.
    Cross-sectional view of the mesh through the elbow
    Cross-sectional view of the mesh through the elbow.
  4. Flow-field insights

    CFD post-processing

    Beyond the pressure-drop number, detailed post-processing shows what is physically happening through the bend, and where that matters for the duct’s structural and wear design.

    Six independent results point to the same finding. Wall and interior pressure, wall shear stress, turbulent kinetic energy, streamline topology, and helicity at two downstream stations all converge on one conclusion: the elbow generates a persistent secondary (Dean-vortex) flow that does not fully die out within the 8 m downstream leg. Peak wall pressure and peak wall shear also coincide at the same location, the outer wall at the bend entrance, which flags that specific point as the one to consider for a wear-resistant liner in abrasive, dust-laden service.

    Streamlines through the elbow, main view with two inset angles showing braiding downstream of the bend
    Streamlines released from a uniform grid at the inlet, traced to the outlet. The braided pattern downstream of the bend is a direct visual signature of secondary (rotational) flow. Select to enlarge.
    Velocity magnitude contour on the bend symmetry plane, a high-speed core hugging the inner wall
    Velocity magnitude, symmetry-plane cut. Peak: 28.4 m/s at the inner wall. Select to enlarge.
    Static wall pressure contour on the duct, high pressure on the outer bend wall and low pressure on the inner bend wall
    Static (gauge) wall pressure, full duct length. Range: −149.7 Pa to +93.2 Pa. Select to enlarge.
  5. Structural analysis (FEA)

    ANSYS Mechanical · shell model from the SolidWorks solid

    The elbow and its adjoining straight stubs, idealized as a shell from the CAD solid, were checked under two load cases: the uniform system design pressure, and the CFD’s own resolved, non-uniform wall-pressure field.

    A. Uniform design pressure
    0.436 MPa peak stress, 8.9 µm peak deformation
    B. CFD-resolved pressure
    0.429 MPa peak stress, 10.5 µm peak deformation
    C. Buckling (governing)
    Mode 1 load multiplier 34.661 at 3.0 kPa
    ≈540yield safety factor, both load cases≈34.7buckling safety factor, governing check4target buckling SF (uplifted from a baseline of 3)

How I checked the results

A base mesh and a refined mesh were both run to measure how much the pressure-drop result still depends on the mesh. Where the CFD and the hand-calculated benchmark disagreed by more than mesh refinement alone could explain, the benchmark itself was audited and found to need a correction, not the CFD.

35.4Pa

final CFD-predicted pressure drop, refined mesh, 85% of the corrected theoretical benchmark

The base and refined meshes were also rerun in a controlled comparison that varied only one solver setting at a time, which confirmed mesh density, not solver configuration, was driving the shift between runs. A two-point Richardson extrapolation from the two mesh states places the mesh-independent result at ≈ 38–40 Pa, or 91–96% of the corrected benchmark.

Pressure drop by mesh state, against the theoretical benchmark

Total pressure drop, inlet to outlet, in Pa. Hover a bar for details.

Total pressure drop: base mesh 27.98 Pa, refined mesh 35.43 Pa, Richardson estimate 39.5 Pa, against a corrected benchmark of 41.7 Pa and the superseded original benchmark of 66.1 Pa.

The hand-calc benchmark needed a correction, not the CFD

The first comparison suggested the CFD was under-predicting pressure drop by a wide margin. Rather than accept that gap, the benchmark itself was audited: its elbow loss coefficient (K = 0.3) had been taken from a generic reference table built around a small, roughly 50 mm reference pipe, about thirty-two times smaller than this 1.6 m duct, and the same reference (Crane TP-410) states its coefficients are diameter-dependent. Recomputed with Crane’s own diameter-dependent method for this duct’s actual size, the corrected coefficient is K ≈ 0.153, roughly half the original value.

66.1 Pageneric K, superseded41.7 Padiameter-corrected K
Independent theory against the CFD result, with status
QuantityTheory (independent)CFD resultStatus
Near-wall y+Target ≈ 11.045Excellent agreement
Straight-duct friction5.68 Pa8.13 PaExplainable, developing flow
Peak velocity25.8 m/s28.4 m/sConsistent, curvature effect
Wall shear stress0.504 Pa1.0–1.3 PaConsistent, curvature effect
Total pressure drop41.7 Pa (corrected)≈39.5 Pa extrapolated (94.7%); 35.4 Pa on the refined meshGood, near mesh-independent

The finished study was then checked back against its own starting specification, item by item: see the integration and specification check in the technical details.

Engineering judgment and recommendations

Beyond the pass/fail result, the study surfaces where a shortcut would have gone wrong, and what to check next.

  • Flagged

    Check diameter-dependence before trusting a generic K-factor

    A textbook elbow loss coefficient looked like a CFD discrepancy until it was traced to the table’s own small-pipe basis. Worth carrying into any future large-diameter ductwork: check a generic coefficient’s stated range before applying it, rather than after.

  • Erosion risk now has a specific location

    Peak wall pressure and peak wall shear coincide at the outer wall, bend entrance. For hot, dust-laden service, that is the natural place to specify a wear-resistant liner or added wall allowance in a production design.

  • Buckling governs by a wide margin, yielding does not

    With a buckling FS of roughly 35 (still 15–17 after a conservative imperfection derating) against a target of 4, and a yield FS over 500, future duct sizing on this line can treat external-pressure buckling as the controlling check from the outset.

  • Two checks recommended before production

    A dedicated FEA mesh-convergence study, mirroring the CFD phase’s base-versus-refined comparison, and a thermal-loading pass at the 200 °C operating temperature, deferred here on the internship’s time budget, are the two clearest next steps.

Technical details

The full design basis, the three structural load cases and the specification check, for engineers who want to go deeper.

Design basis Geometry, flow, thermal and structural inputs, material and flange data
Design basis of the duct section
ParameterValue
Geometry
Nominal internal diameter1.60 m
Shell wall thickness6 mm
Elbow angle / bend radius90° / 2.40 m (R = 1.5D)
Upstream / downstream straight length5.0 m / 8.0 m
Flow
Design flow rate42.50 m³/s
Design velocity21.1 m/s
Reynolds number≈ 9.7 × 105
Thermal / structural
Design gas temperature200 °C
Design (suction) pressure−3.0 kPa
Material
Shell materialS235JR (EN 10025-2)
Flange
Outer diameter / bolt-circle diameter1.80 m / 1.71 m
Bolts32 × M20
Structural load cases A. Uniform design pressure, B. CFD-resolved pressure, C. Buckling check (governing)

A. Uniform design pressure

Load case

The elbow and adjoining straight stubs, idealized as a shell from the SolidWorks solid, are loaded with the uniform −3.0 kPa system design pressure used to size the duct wall in the design basis.

0.436 MPapeak stress, clean elbow reading8.9 µmpeak deformation

Stress stays well below yield everywhere. A separate local maximum at the model’s cut edges is a known artifact of the simplified boundary condition there, confirmed by probing the elbow directly rather than reading the raw global maximum.

Stress contour on the elbow shell, uniform design pressure case
Stress contour, uniform-pressure case. Select to enlarge.

B. CFD-resolved pressure

Load case

The same elbow is reloaded with the CFD’s own converged, spatially non-uniform wall-pressure field (shown in step 4) in place of the uniform assumption, on top of the system design pressure. On import, the field’s minimum matched the CFD’s own minimum almost exactly (−149.74 Pa vs. −149.7 Pa).

0.429 MPapeak stress, clean elbow reading10.5 µmpeak deformation

The realistic load redistributes stress rather than raising it: elevated at the outer wall, relieved at the inner wall, tracing directly to the bend’s flow asymmetry from step 4. The verdict does not change, but the result is now a physically explained one rather than an assumed one.

CFD wall pressure field mapped onto the FEA shell model of the elbow
The CFD wall-pressure field, imported directly onto the shell. Select to enlarge.

C. Buckling check (governing failure mode)

Governing check

Yielding carries an overwhelming margin in both load cases, so buckling under external (suction) pressure, identified as the governing mode from the very first hand calculation in the design basis, is the check that actually sizes this duct.

34.661Mode 1 load multiplier at 3.0 kPa

This is an idealized linear eigenvalue result, not a pressure the real shell can be assumed to reach. Real thin shells buckle at a lower pressure once fabrication imperfections are present, which is exactly why the design basis itself uplifted its target safety factor from a baseline of 3 to 4 for this duct’s field-rolled, site-welded construction. Even a substantial, conservative derating for imperfection sensitivity leaves a margin of roughly 15–17, still comfortably above that target of 4.

Linear eigenvalue buckling Mode 1 shape, multi-lobe ovalization of the outlet leg
Buckling Mode 1: multi-lobe ovalization, concentrated on the longer (8 m) outlet leg. Select to enlarge.
Integration and specification check The finished study checked back against its own starting specification: five items

The two analysis phases connect directly: the realistic load case in step 5 is the CFD’s own converged wall pressure, carried into the structural model. The finished study was then checked back against its own starting specification, item by item.

Specification check: item, how it was met, and status
Specification itemHow it was metStatus
Standards basisISO/EN used throughout wherever an equivalent exists (ISO 7005-1, EN 10025-2, ISO 898-1, EN 13445-3).Satisfied
TraceabilityEvery numeric input tagged back to a cited reference or a stated assumption.Satisfied
VerificationEvery CFD and FEA result checked against independent hand calculation.Satisfied
Project disciplineSingle representative geometry, isothermal, single-phase, as scoped at the outset.Satisfied
ReusabilityBuilt on general, parametric formulas rather than hard-coded numbers, so the same workflow rescales to a different duct.Satisfied

Questions about this project?

I’m happy to walk through the setup and results in more detail.