REVIEW 3 major objections 5 minor 12 references
Numerical Investigation of the Effect of an Oblique Flow Entry on the Pressure Losses in Square Channels
T0 review · 3 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read Oblique entry trips turbulence in square channels at Re≈2000
desk verdict Solid engineering CFD paper with a real new finding (length-dependent KObl) and an under-supported causal claim about shear-layer-driven transition. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the shear layer that forms between the recirculation zone (created by flow separation at the channel entrance) and the freestream jet. Its instability produces coherent vortical structures that shed and persist downstream, acting as a self-sustained turbulence source. The argument is carried by the dimensionless oblique pressure-loss coefficient KObl = ΔPObl / (ρU1²/2), defined by subtracting the axial-entry pressure drop; the paper shows that this subtraction is only valid while both flows are in the same regime. An auxiliary but decisive piece of machinery is the leading-edge rounding radius r/d, which the authors measure from a real cordierite catalytic converter (about 2.9% of channel width) and vary in simulation, showing that it controls the shear-layer instability and therefore KObl.
What would settle it
Measure the pressure drop and velocity fluctuations across a square channel with a 40° oblique entry at Rec ≈ 2000 while systematically increasing the turbulence intensity of the inlet flow from near zero to several percent; if the sharp rise in KObl shifts to lower Reynolds numbers or changes magnitude with inlet turbulence, the claim that the shear layer alone triggers transition is falsified.
Extended reading notes
Core claim
On its own terms, the paper's central claim is that an obliquely entering flow into a square channel is not a laminar flow at moderate channel Reynolds numbers. At entry angles above roughly 20 degrees and Rec around 2000, the recirculation zone at the channel entrance creates a shear layer that becomes unstable and continuously sheds coherent eddies; these persist along the channel and drive a transition to turbulence. This transition is triggered by the shear layer itself, with no turbulence prescribed at the inlet. Because the axial-entry flow at the same Reynolds number remains laminar and developing, subtracting the axial pressure drop (as the definition of KObl does) no longer isolates the entrance loss: KObl then contains the difference between laminar and turbulent friction over the whole channel length, making it length-dependent. The paper further shows that rounding the leading edge, as found in real cordierite monoliths, stabilises the shear layer, delays the transition, and lowers KObl, with a radius of 0.01d bringing LES into agreement with experiments. Laminar flow models, steady or unsteady, fail in the transitional cases, while the k–ω SST RANS model predicts pressure losses within 5% of experiments for most conditions, except where transition occurs (differences up to 40%).
Load-bearing premise
The inlet is a flat, uniform, turbulence-free laminar velocity profile; if real upstream flow contains residual turbulence or pulsation, the shear-layer transition could occur earlier or by a different mechanism, changing both the predicted losses and the comparison with experiments.
Editorial extensions
If this is right
- Oblique entry into square channels cannot be treated as laminar at moderate Reynolds numbers: for entry angles above about 20°, the shear layer generates turbulence even at Rec ≈ 2000, so design correlations must include a transition term.
- The oblique pressure-loss coefficient KObl depends on channel length when oblique flow is turbulent and axial flow is laminar; a single KObl per angle is not valid.
- Sharp-edged simulations overestimate pressure losses at higher Reynolds numbers and angles; rounding the leading edge to realistic values (r ≈ 0.01d–0.05d) reduces losses and brings LES into agreement with experiments.
- For engineering pressure-drop predictions, RANS with the k–ω SST model is accurate to within 5% for most conditions; steady or unsteady laminar models are unsuitable for oblique entry even at Rec = 2000.
- The turbulence generated by oblique entry could be exploited for enhanced mixing in monolith channels, at the cost of increased pressure drop.
Reading between the lines
- The same shear-layer mechanism should appear in channels of other cross-sections (e.g., circular or hexagonal) whenever the inlet flow has a transverse component; the transition threshold likely scales with the ratio of transverse to axial momentum rather than with Reynolds number alone.
- The flat, turbulence-free inlet used here probably makes the predicted transition a conservative lower bound: residual upstream turbulence or pulsation in real engines would likely trip the shear layer earlier, so the length-dependence of KObl may be even more pronounced in service.
- If KObl is length-dependent, device-level pressure-drop models for monoliths should either resolve the entrance region explicitly or use an entry-loss term that accounts for the distance over which turbulent friction acts, rather than a fixed entry-loss coefficient.
- A testable design extension follows directly: vary the inlet chamfer or rounding radius as a design parameter to trade pressure drop against mixing; the paper's data suggest a radius around 1% of channel width already suppresses most of the shear-layer loss.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents CFD simulations of a single square monolith channel with flow entering obliquely. LES (dynamic Smagorinsky), RANS k-omega SST with a transition model, and steady/unsteady laminar models are compared over Rec=500-3000 and alpha=0-70 degrees. The authors report that oblique entry creates a recirculation region and shear layer at the channel entrance, that the shear layer sheds coherent structures and promotes transition to turbulence even at Rec~2000, that this makes the oblique-entry loss coefficient KObl length-dependent, and that sharp-edged simulations overpredict experimental losses because the real leading edges are rounded. The simulations are compared with the Shah-London developing-flow correlation, the K-W, Persoons, and Quadri correlations, and experimental data from Quadri et al. and Mat Yamin.
Significance. If the central transition claim is correct, the paper overturns a common laminar-flow assumption for monolith channels at moderate Reynolds numbers and identifies a physical mechanism (shear-layer shedding from oblique entry) that explains both the Reynolds-number and length dependence of KObl. The study has several genuine strengths: the axial-entry cases are validated against the Shah-London developing-flow correlation; the outlet-length effect is checked; the LES mesh study reports both global pressure-drop convergence and the Celik LESIQnu index; and no constants are fitted to the target pressure-loss data, with the edge-rounding radius being varied and the best-agreeing value identified post hoc. The experimental comparison over four channel lengths is a valuable addition to the literature. However, the central claim rests on LES with a disturbance-free inlet and only two residence times of averaging, and the mesh study does not assess the convergence of the unsteady quantities used to infer shedding. These issues are fixable but currently leave the main physical claim underdetermined.
major comments (3)
- [III.A.3 / II.E / II.F] The central claim that the oblique-entry shear layer produces continuous shedding of eddies that persist downstream at Rec~2000 rests on LES in which the inlet is a flat, turbulence-free profile (Sec. II.E). Since no physical disturbances are seeded, the perturbations that initiate the shear-layer instability are numerical in origin. The mesh assessment of Sec. II.F establishes convergence of the s1-s3 pressure drop and LESIQnu, but it does not test convergence of the fluctuation statistics used to support the shedding claim: shedding frequency, spectral content, transition-onset location, or the streamwise evolution of U'_RMS. The observed transition could therefore depend on mesh-dependent numerical noise, particularly in a subcritical regime, and the causal attribution to the oblique-entry shear layer is underdetermined. I recommend either a mesh-4 replicate showing the same shedding frequency and U'_RMS streamwise evolution, or a test with a small controlled inlet disturbance to confirm that the instability is physical rather than triggered by numerical noise.
- [II.C / III.A.3 / III.B] All reported time-averaged pressures and fluctuation profiles are averaged over only two residence times of the full domain (Sec. II.C). No convergence check of the running time average is reported for the pressure drop or U'_RMS, and no shedding frequency or Strouhal number is examined. If low-frequency or intermittent transitional behavior is present, two residence times may alias the statistics and materially change KObl. Please report the time-history of the running mean over a longer interval (for example, ten residence times) for a representative transitional case, and quantify the uncertainty of the time-averaged pressure drop used in Figs. 20-26.
- [III.B.5 / IV] The conclusion that leading-edge rounding explains the discrepancy with experiments is based on rounded-edge simulations at only alpha=0 and 40 degrees, with the radius varied, and on a single comparison at Rec=1500 and alpha=45 degrees for L=27d. The statement that the r=0.01d case agrees within experimental error is thus not demonstrated across the parameter range over which the discrepancy appears (for example, alpha=30-60 degrees, Rec=2000-3000, and all monolith lengths). Either extend the rounded-edge matrix or temper the conclusion to a demonstration of sensitivity rather than a full quantitative explanation of the experimental differences.
minor comments (5)
- [Eq. (22)] In Eq. (22), setting epsilon = 1 for a smooth channel is dimensionally unclear; since epsilon is a roughness height, smooth walls should correspond to epsilon = 0, and the plotted comparison should state the relative roughness used.
- [III.A.2] There is a typo in the text near Fig. 13: 'seperation' should be 'separation'.
- [III.A.3] Figure 15 defines U'_RMS but does not state the time window over which the root-mean-square is computed; this is important given the two-residence-time averaging stated in Sec. II.C.
- [III.B.3 / Figs. 23-26] The experimental data points and error bars in Figs. 23-26 are not described in terms of the number of repeated measurements or the source of uncertainty; a sentence explaining the experimental uncertainty would help the reader judge the agreement.
- [II.D.1 / II.D.2] The numerical schemes in Star CCM+ are described only by order of accuracy; specifying the actual discretization schemes (e.g., bounded central differencing for LES, upwind for RANS) would improve reproducibility.
Circularity Check
One post hoc edge-radius selection; central transition and KObl-length claims remain independent.
-
fitted input called prediction
[Section III.B.5, 'Effect of Rounded Walls' (Figs. 31-33)]
"To study the effect of the rounding of the leading edge on the channel pressure losses, three new geometries were produced with the leading edge rounded with a radii of r = 0.01d, 0.025d, and 0.05d ... For the case with r = 0.01d, the agreement with experiments is within the experimental error. This confirms that the presence of a rounded edge in the experimental setup could be the reason for the differences between measured and predicted pressure losses for higher Reynolds numbers and oblique entry angles."
The rounding radius that gives agreement is selected by scanning three candidate values (0.01d, 0.025d, 0.05d) and then the matching case, r = 0.01d, is used to 'confirm' that edge rounding explains the LES/experiment discrepancy. This is a post hoc fit: the agreement is forced by selection, not an independent prediction. Notably, the microscope measurement reported in the same section gives a radius of about 2.92% of channel width (~0.03d), not the 0.01d value that matched the data, reinforcing that the matching radius was chosen by fit rather than by the independent geometric measurement. The paper's central claims about shear-layer shedding and the length dependence of KObl do not depend on this radius choice.
full rationale
The core derivation chain is self-contained. No constants are fitted to the target pressure-loss data: axial-entry cases are validated against Shah-London developing-flow losses and Colebrook turbulent friction, and KObl comparisons use external correlations (K-W, Persoons, Quadri) and experimental datasets (Quadri et al., Mat Yamin) rather than values tuned inside the simulations. The persistent-eddies and transition claims are direct LES observations made with a flat, turbulence-free inlet; they are not defined into existence by the pressure-loss formula. The paper also explicitly acknowledges that Eq. 8 includes downstream laminar/turbulent friction differences, so the length dependence of KObl is a consequence the authors identify and explain rather than a hidden input. Same-group references (Quadri et al., Mat Yamin) are experimental evidence, not load-bearing theoretical authority, so they do not constitute circularity. The only mild circular element is the post hoc selection of the rounded-edge radius that matches experiment, which is secondary to the central claims and therefore warrants only a low score.
Assumptions & free parameters
free parameters (2)
- Leading-edge rounding radius, r =
0.01d (best agreement with experiments); measured representative radius about 0.029d
- Time-averaging window =
2 residence times
assumptions (6)
- domain assumption Air is modeled as incompressible, isothermal, Newtonian fluid.
- domain assumption Flat, uniform, turbulence-free velocity profile at the inlet boundary.
- domain assumption Side walls are translational periodic boundaries, representing an infinite array of identical channels.
- domain assumption Channel walls are smooth with no-slip and zero roughness.
- domain assumption The measured edge profile from one cordierite monolith is representative of the experimental monoliths used by Quadri et al. and Mat Yamin.
- standard math Navier-Stokes, LES-filtered equations, and the Dynamic Smagorinsky SGS model as implemented in Star CCM+ correctly represent the flow physics.
Cite this review
Pith. "Pith review of Numerical Investigation of the Effect of an Oblique Flow Entry on the Pressure Losses in Square Channels." pith.science (2026). https://pith.science/paper/CX3OYSQR
@misc{pith2026250207678,
author = {Pith},
title = {Pith review of: Numerical Investigation of the Effect of an Oblique Flow Entry on the Pressure Losses in Square Channels},
year = {2026},
howpublished = {\url{https://pith.science/paper/CX3OYSQR}},
note = {Machine review of arXiv:2502.07678}
}
read the original abstract
Flows in square channels are common in applications, such as automotive after-treatment systems and heat exchangers. Flows with axial flow entry are well understood, but for oblique flow entry, there is no clarity on the additional pressure loss magnitude or the flow regime. Laminar flow is often assumed, even though flow separation at the channel entrance can cause a transition to turbulence. Here, the impact of oblique flow entry on the flow is investigated using LES (Large Eddy Simulation) and RANS (Reynolds Averaged Navier Stokes) models, and their advantages and limitations are identified. The LES simulations show that the shear layer at the channel entrance produces continuous shedding of eddies that persist downstream even at moderate channel Reynolds numbers (~2000). The predicted pressure losses mostly agree with experimental data. The differences observed for some parameters are attributed to the difficulty of accurately replicating the experimental geometry. It is shown that LES results are susceptible to the rounding of the leading edge (present in experiments). Including edge rounding improves the pressure predictions. RANS simulations predicted pressure losses within 5% of experimental values for most cases, apart from where transitional flow was observed (resulting in differences up to 40%). This study provides insights into the flow structure and sources of pressure losses in square channels and highlights the importance of understanding key flow and geometric features when using LES to predict complex flows involving flow separation and shear layers.
Figures
Figures from the paper (21 more)
Reference graph
Works this paper leans on
-
[1]
29 The mesh cell size acts as a filter sorting the eddies by size
LES Setup The LES approach segregates the turbulent flow structures (eddies) into those that are sufficiently large to be resolved on the mesh, and the remaining, smaller, eddies which are mod- elled by a sub-grid scale turbulence model. 29 The mesh cell size acts as a filter sorting the eddies by size. Consequently, a more refined mesh results in a small...
-
[2]
RANS Setup LES is considered a superior modelling tool for capturing flows with shear; however, RANS models are frequently used because of their computational efficiency. To compare the ac- curacy of different models for the flow in the considered range of Reynolds numbers, a RANS model ( k − ω SST) has been used for several configurations. The RANS gover...
-
[3]
Laminar Setup Both the LES and RANS models account for turbulence and therefore perform favourably when turbulence is present. Since some simulations are conducted at Reynolds numbers associated with laminar flow transitioning to turbulence, these models may artificially induce turbulence under conditions where it would not develop in reality. To assess w...
work page 2000
-
[4]
This causes the formation of a vena contracta immediately post-entrance (see Fig
Recirculation Region and Shear Layer In the case of an axial entry, the flow path contracts as it en- ters the channel. This causes the formation of a vena contracta immediately post-entrance (see Fig. 8). For contraction ratios typical to monoliths and catalytic converters ( ≈ 0.75 - 0.9), the effect is not very pronounced, and the separation zones forme...
-
[5]
Secondary Flow Formation of secondary vortices in non-circular ducts has been observed in many experimental and numerical studies,36–38 and is attributed to the transverse turbulent stress gradients. These vortices play a significant role in shaping the overall flow structure and enhancing mixing within the channel.19,21 At lower Reynolds numbers, the sec...
work page 2000
-
[6]
Turbulence Development The discussed flow features (flow acceleration due to sepa- ration, shear layer, and secondary flow formation) contribute to the transition of the flow to turbulence with increasing Reynolds number and/or oblique entry angle. Fig. 15 shows the average cross-sectional distribution of resolved velocity Numerical Investigation of the E...
work page 2000
-
[7]
Pressure Drop for Axial Flow Entry Simulations Fig. 16 shows the pressure distribution in the domain for axial entry ( α = 0◦) simulations with Rec = 1000, 2000, and
work page 2000
-
[9]
Pressure Drop for Oblique Flow Entry Simulations Figs. 18 and 19 show the pressure distribution along the channel for cases with an oblique flow entry for Rec = 1000 and 3000. For lower angles, the cross-sectional average pres- sure drops at the entrance to the channel, but recovers from the contraction and separation of the flow back to the developing fl...
work page 2000
Show all 12 references
-
[10]
Four common correlations used are detailed in Table IV
Comparison with Experimental Data and Existing Correlations The results can be compared with correlations and exper- imental data that exist within the literature. Four common correlations used are detailed in Table IV. Values of A and n for the correlation by Quadri et al.13 ...
2000
-
[11]
III A that for higher values of Reynolds number and α, large vortical structures form within the shear layer
Effect of Turbulence Model It has been shown in Sec. III A that for higher values of Reynolds number and α, large vortical structures form within the shear layer. This suggests that the shear layer makes a sig- nificant contribution to the total pressure loss, in addition to t...
2000
-
[12]
Monolith structures, materials, properties and uses,
Effect of Rounded Walls The geometry of the leading edge at the channel entrance plays a crucial role in determining the behaviour of fluid flow around it.4 A rounded edge allows for a gradual change in the flow direction, reducing the adverse pressure gradient and en- abling ...
2001 arXiv
-
[3000]
To assess the importance of the boundary layer develop- ment, and the effects of the oblique entry, the results are com- pared with the established correlations for frictional losses in a square channel with laminar and turbulent flow. For laminar flow, the fully developed, fr...
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.