Turbulent Zone Research Articles

Experimental investigation on flow and heat transfer characteristics of helium in rectangular narrow slit channel

Helium is utilized as coolant in high temperature gas cooled reactors. Rectangular narrow slit channels are commonly present in the core. The flow and heat transfer properties of helium significantly affect the core temperature and flow distribution. This paper experimentally investigates the flow and heat transfer properties of high temperature helium in rectangular narrow slit channel. The Reynolds numbers ranged from 468 to 9357, the temperature ratio of wall to bulk from 0.91 to 1.22, the helium outlet temperature up to 945 K, and the maximum heat flux density up to 0.101 MW/m2. In the experiment, the total and local convective heat transfer coefficients, along with the friction factors, are measured. Correlation equations for the Nusselt number and friction factor with respect to Reynolds number are derived. The primary findings are as below: the friction factors of helium in rectangular narrow slit channels are significantly larger than those calculated by current empirical correlations. In the turbulent zone, the measured local Nusselt numbers are in excellent accordance with the Gnielinski correlation; however, a significant discrepancy exists in the laminar zone. New flow heat transfer correlation equations are proposed on the bases of experimental data. Upon comparison, the total Nusselt numbers fall in the error of 10 %, and almost all local Nusselt numbers and the friction factors fall in the error of 20 %.

Evolution of amplitude and longitude phase of tachocline Rossby waves diffusing to the photosphere

ABSTRACT Physics of magnetohydrodynamic (MHD) Rossby waves in the tachocline-layer were studied by Dikpati et al., using a fluid-particle-trajectory approach along with solving vorticity and induction equations. By extending that model to include a hydrodynamic turbulent convection zone (CZ), we examine how MHD Rossby waves generated in the tachocline might diffuse upward through the CZ to solar surface. We find that pure hydrodynamic Rossby wave amplitudes decline with height due to viscous diffusion at a rate that is independent of viscosity and increases with longitude wavenumber. Fast MHD Rossby waves amplitude declines faster with height for increasing toroidal field, due to their longitude-phase shifting with height, which increases dissipation of kinetic energy in the wave velocities. Slow MHD Rossby waves decline even faster with height because their longitude-phase shifts more rapidly with height, due to their slow phase speed. We conclude that low wavenumber HD and fast MHD Rossby waves, originating in the tachocline, might be detected at the photosphere, but slow MHD Rossby waves should be virtually impossible to detect. We infer from fluid particle trajectories that wave amplitudes declining with height and longitude phase shifting with height associated with decline, implies a powerful mechanism for tangling of magnetic fields, distinct from convective turbulence effects. This could cause a sustained or dissipative local dynamo action triggered by Rossby waves.

New generalized expressions for forced convective heat transfer coefficients across the flat plate

This study presents a novel generalized model for forced convective heat transfer coefficients (CHTC) over flat plates, addressing key limitations in existing models. Current expressions often provide averaged values and fail to capture the local variations in heat transfer, leading to inaccuracies, particularly in transitional zones. Moreover, there is no consensus on boundary layer classifications, which further complicates the development of reliable CHTC models. To address these issues, we conduct extensive CFD simulations to systematically analyze the combined effects of plate length (L), wind velocity (v), position along the plate (x), and temperature difference (Δt) on local CHTC. The simulations span a wide range of conditions (0.25 m ≤ L ≤ 8 m, 3 m/s ≤ v ≤ 25 m/s, and 5 °C ≤ Δt ≤ 80 °C), creating a comprehensive database of local CHTC values. Based on this data, a new generalized CHTC expression is derived as a function of L, v, x, and Δt, accurately predicting CHTC across laminar, transitional, and turbulent zones of flat plats. The accuracy of proposed CHTC expression is verified through thorough in-sample and out-of-sample evaluations, achieving R2 values between 0.973 and 0.999, with maximum deviations ranging from 2.6 % to 8.5 %. These research methods and conclusions can provide valuable references for optimizing thermal processes in applications such as industrial heat exchangers, electronics cooling, and building heating systems, where precise heat transfer modeling is critical for enhancing efficiency and system design.

Steady secondary flow in a turbulent boundary layer past a slender axisymmetric body

The turbulent boundary layer in a viscous incompressible fluid developing longitudinally past the surface of a thin cone or cylinder at a finite distance from the laminar–turbulent transition zone is studied. The characteristic Reynolds number, determined from the external flow velocity and the length of the body, is assumed to be large, and the thickness of the boundary layer is small and comparable to the radius of the body. The asymptotic method of multiple scales is used to find solutions to the Navier–Stokes equations. Instead of the traditional decomposition of the solution into time-averaged values and their fluctuations, the velocities and pressure are expressed as an asymptotic series consisting of steady and perturbed terms. As a result, the viscous steady flow (‘secondary’) that arises in the boundary layer as a mandatory component of fast turbulent fluctuations was described. Analytical and numerical solutions for the radial steady velocity are presented, describing the self-induced suction of fluid from the external flow into the boundary layer. Further analytical solutions are obtained for the longitudinal and circumferential velocities, which differ markedly from the laminar regime. The solutions found are somewhat similar to the degenerate (one-dimensional) case of self-sustaining longitudinal thin structures in turbulent shear flows. A qualitative comparison with direct numerical simulations is presented.

Local-Variable-Based Transition Model for Hypersonic Flows Considering Wall Mass Injection

Pyrolysis gas injection during the ablation process will affect the hypersonic boundary-layer transition. This paper proposes a three-equation γ-Reθ-fb transition model that considers wall mass injection. First, the boundary conditions of wall mass injection are established and verified. Second, based on analysis of the compressible boundary-layer solutions under wall mass injection, new transition correlations are developed. This allows an injection factor transport equation to be constructed using strictly local variables, and this equation is coupled with the γ-Reθ transition model. Additionally, to model the wall mass injection effects in the fully turbulent zone, the wall boundary conditions for the specific turbulence dissipation rate are amended. The results of numerical simulations show that the transition onset locations and the changes in the heat transfer rate are reasonably consistent with experimental data.

Impact Mill

An impact mill has been developed to produce powders from shavings of refractory metals using impact grinding method for reuse in electrometallurgy in devices with screw feed, for example, in 3D printers. The proposed device provides high uniformity of grinding with a minimum content of dust fraction and impurity content at low technical and economic costs. The result is achieved by using a Laval nozzle, which operates in the supersonic jet formation mode. In the area of the first Mach disk there are rod fenders arranged in a cascade, and the impact plate is located in the turbulence zone and is equipped with winglets with holes for separating crushed metal.

Near-field propagation of a flat-topped gaussian beam: analysis in weakly ‎turbulent atmospheres

We investigate the analysis of optical communication employing the shape beams, definitely concentrating on the effect of factors on the physical characteristics of a Flat Top Gaussian (FTG) beam. The propagation of a beam of laser light in the atmosphere layer is disposed to being affected by many optical phenomena such as scattering, absorption, and turbulence. These phenomena arise from differences in the scintillation index and the intensity modes realized in the source and receiver planes. This work is a numerical investigation of the FTG laser beam as it beams trekking through a zone of low turbulence using open-source software. This simulation will be conducted using a mathematical model derived from the split-step beam propagation technique. The intensity distributions of the source plane and the average intensity is received in air turbulence are computed, with an extra contour at the transducer plane. The Rytov approach to develops the scintillation index, structural constant, source size, and supplementary parameters to measure the weak turbulent model. Furthermore, these factors are studied in the context of near-field propagation. Also, the impression of the scintillation and wander of the beam is assessed. All simulated results are analyses and contrasted with the TEM00 Gaussian beam. At last, these discoveries are compared with the data obtained in the experimental piece of the study, additionally by applied in laser applications.

Effect of cloud chemistry on seasonal variations of sulfate and its precursors in China

Cloud chemistry is of paramount importance in the secondary production of atmospheric aerosols, influencing the spatial-temporal distribution of gases and aerosols in the atmosphere. Using WRF/CUACE (China Meteorological Administration Unified Atmospheric Chemistry Environment), this study assesses the seasonal impacts of cloud chemistry on the concentrations of SO2, sulfate, as well as two oxidizers, H2O2 and O3, in the most east-central areas of China, including four key pollution zones (the North China Plain (NCP), the Yangtze River Delta (YRD), the Pearl River Delta (PRD), and the Sichuan Basin (SCB)). Near the surface, H2O2-oxidation was the dominant pathway for cloud chemistry in four key pollution zones in four seasons. H2O2 consumption is most pronounced in summer, especially in the SCB and NCP, while O3 consumption peaks in autumn, particularly in the PRD and southeastern coastal areas. While at higher altitudes, oxidation by O3 and H2O2 is compatible with the cloud chemistry process. Near the surface, cloud chemistry consumes SO2 ranging from approximately 0.1 ppb–5.0 ppb, resulting in the generation of about 6.0–25.0 μg m−3 of sulfate. Higher SO2 reduction and sulfate increase are in both summer and winter, especially for the SCB and NCP in summer, and the SCB in winter. Vertically, the cloud chemistry process primarily concentrates its influence on SO2 and sulfate concentrations below 5 km, particularly within the turbulent zone of the troposphere below 2 km in all the four pollution zones and four seasons. The most notable seasonal variation occurs in the NCP compared to other zones. This study also shows that cloud chemistry effectively improves the seasonal simulation accuracy of SO2 and sulfate, resulting in improved correlation and a notable reduction in RMSE.

СОВРЕМЕННОЕ СОСТОЯНИЕ ИССЛЕДОВАНИЙ И ПРОГНОЗИРОВАНИЯ ВЛИЯЮЩЕЙ НА АВИАЦИЮ ТУРБУЛЕНТНОСТИ В СВОБОДНОЙ АТМОСФЕРЕ

Three main physical processes are responsible for generation of turbulence affecting the aircraft: (1) hydrodynamic instability of the airflow stratified in wind and temperature (Kelvin-Helmholtz instability), which manifests itself in origination of turbulent patches called clear-air turbulence (CAT) in the free atmosphere; (2) breaking of large-amplitude mountain waves over the obstacles and at their leeward slopes, which causes so called orographic turbulence or mountain wave turbulence (MWT); (3) outflow of highly turbulent air from convective storms, which leads to the formation in their vicinity of turbulent zones of deep convection, that is, convectively induced turbulence (CIT). In the paper, the state of the art in studying and forecasting the turbulence generated by the mentioned physical processes is considered. The present short-range forecasting and nowcasting techniques are characterized by transition to the universal standard for intensity of aviation turbulence, that is, to the energy dissipation rate (EDR) defined in terms of the theory of homogeneous and isotropic turbulence, and by taking into account the informativity of empirical turbulence diagnostics. For the CIT nowcasting, special techniques for convection diagnosis are developed based on the satellites, radars, and lightning registers data.

Numerical investigation on the enhanced heat transfer characteristics of non-Newtonian viscoelastic fluids in microchannel with chip arrays

Elastic turbulence is an efficient heat transfer method at microscale. In this work, we conduct a systematic investigation on the flow and heat transfer performances of viscoelastic polymer solutions in microchannels with chip arrays by means of high-fidelity numerical simulations. Focusing on the effects of fluid elasticity [i.e. polymer relaxation time (λ)] and flow Reynolds number (Re), we find that the synergistic interplay between inertial and elastic forces is strategically exploited to amplify convection, thereby facilitating an enhanced heat transfer performance. For highly elastic fluid, the Nusselt number (Nu) can be improved by 76.13 % compared to that of less elastic fluid, while the pressure drop is reduced appropriately. In addition, as the fluid elasticity increases, the pressure drop decreases and eventually converges to a constant. In the zone of elastic turbulence, the fluid reinforces the convective heat transfer, which can be ascribe to the deformation of long-chain polymer molecules. Particularly, we consider the effect of temperature-dependent λ on the heat transfer performance. Our results demonstrate that the increased temperature-dependent λ decreases the value of Nu, and the variation between values of variable and constant λ is reduced with increasing Re or λ. Hence, in some cases of large fluid elasticity or high Re, it is possible to safely ignore the impact of temperature-dependent λ with an acceptable accuracy. The present work provides important insights into the enhanced heat transfer via elastic turbulence, which could guide the applications in the field of chip cooling.

High-fidelity simulations of boundary layer transition on an underwater axisymmetric body with stability theory

During the cruise status of an underwater high-speed vehicle, the flow around the vehicle's head typically transitions from a laminar to a turbulent state, triggering flow noise that can interfere with the normal operation of sonar. In order to accurately investigate this flow noise through numerical simulation, a high-fidelity turbulent flow field solution is essential. Common traditional turbulence numerical simulation methods, such as unsteady Reynolds-averaged simulation and large eddy simulation (LES), struggle to capture high-frequency turbulent fluctuations accurately due to their inability to directly resolve small-scale eddy structures, which results in compromising the simulation accuracy of high-frequency flow noise. To address this issue, this paper employs direct numerical simulation (DNS) to achieve high-fidelity resolution of the turbulent flow field, thereby enabling a more accurate assessment of flow noise distribution on the vehicle's surface. Meanwhile, considering significant computational resources required to solve the entire flow field in an underwater high Reynolds number environment, this study also incorporates the fixed transition modeling method and stability theory to confine the DNS computational domain to the vicinity of the transitional zone to improve simulation efficiency. Comparative analysis of flow noise monitoring results in the laminar, transitional, and turbulent zones revealed that the flow noise source in the laminar zone exhibits the lowest amplitude across all frequencies, while the flow noise source in the transitional zone features the highest amplitude, approximately 10 dB higher than that in the turbulent zone. Moreover, significant amplitudes in high-frequency components (above 30 kHz) are detected in both the transition and turbulent zones. Additionally, this study employs LES with the Smagorinsky model to simulate the flow field within the same computational domain as DNS, demonstrating the limitations of the Smagorinsky model-based LES in capturing high-frequency flow noise.

Numerical study on the effect of high efficient cooling nozzles and varying cooling intensity on metallurgical transport behaviors during the slab continuous casting

In this study, the spray characteristics of the cooling water flux of the traditional single nozzle and novel dual nozzle were innovatively and effectively incorporated into a 3D/2D flow-temperature-concentration segmented model. The model was used to investigate the effects of spray characteristics on the flow, heat distribution, solute transport, solidified shell, and mushy zone of steel in the continuous casting. The results showed that various flow patterns of liquid steel in the turbulent zone significantly affected the temperature and carbon concentration distribution. Until Zone 7, the cooling water fluxes in the six cases remained unchanged. The peak temperatures of cases 1 and 4 in Zone 7 were 1253.11 and 1273.51 K, respectively, indicating that the spray characteristic was the primary cause of the variations in slab surface temperatures. The cooling water fluxes in the six cases change from Zone 8 onward. The six cases had differences of −21.39, 17.87, 46.95, −22.08, 16.86, and 45.68 K between the initial and final temperatures in Zone 8, respectively, meeting the requirement of keeping the maximum temperature recovery rate of slab surface under 100 °C/m. However, the final temperatures for cases 2, 3, and 6 were 1225.62, 1196.50, and 1218.91 K, respectively. These temperatures fall within a realistic plastic temperature range, which must be higher than 1220 K. As a result, the plastic cracking in the slab was possible. The maximum temperature gradient differences of the six cases between feature lines at the end of Zone 8 were 0.793, 0.814, 0.829, 0.185, 0.179, and 0.179 K/mm. These results showed that the optimization effect on the temperature gradient difference was insignificant as the cooling intensity rose. However, the carbon concentration uniformity was marginally improved by increasing the cooling water flux. Finally, a state-owned steel company in China chose case 5 (dual nozzle with moderate cooling intensity at the slab solidified end) as the optimum option for the slab continuous casting caster. The dual nozzle enhanced and promoted the efficient and homogeneous production of the metallurgical process.

Enhanced thermal effectiveness of square duct with V-type double-baffles: Numerical study

The article puts forward three-dimensional computational research on heat transmission augmentation within a square channel containing 45o V-type double-baffles positioned on the lower and top parts at regular intervals in the turbulence zone for Reynolds numbers (Re) that vary from 3000 to 20,000. The primary goal of this research is to increase the thermal effectiveness and relative Nusselt number (Nu/Nu0), in order to conserve energy and reduce the size of the heating or cooling system. The simulations utilize a finite volume approach in common with the SIMPLE algorithm, whereas the turbulent model used is the realizable k–ε. The baffles are designed to be separated vertically for reducing pressure loss. Both single V-baffles and double V-baffles have four relative pitches (PR = 0.4, 0.5, 0.6, and 1.0) and height/blockage ratios (BR = 0.05, 0.1, 0.15, and 0.2), with a fixed attack angle (α) of 45o. The computational findings show that both V-baffles are capable of producing the primary vortices, but only the double V-baffles have the ability to provide the impinging streams onto the wall, cooling the region behind the baffles. This suggests that the double V-baffles not only boost heat transmission but also reduce frictional loss. When compared to a single V-baffle, the double ones enhance heat transfer by an average of 1.04–9.94% while decreasing frictional loss by an average of 9.88–31.73%. The thermal effectiveness factor (TEF) of the double V-baffles ranges from 1.03 to 3.21, and its peak value of around 3.21 is for PR = 0.4, BR = 0.05, at lower Re.

Turbulence generation in transonic shock diffraction

In this Letter, we have quantified the zones of turbulence in transonic shock diffraction across a 90° step corner using Lamb vector analysis. We have analyzed the nature of the three-dimensional turbulence structures using invariants of the velocity gradient tensor. From our results, we conclude that the vortex-induced shock, the lambda shock on the shear layer, and the Kelvin–Helmholtz instability forming on the shear layer are the main contributors in turbulence generation. In our study, we have shown that two-dimensional simulations can resolve the gasdynamic behavior accurately, though three-dimensional simulations are required to understand turbulent structures.

Investigation on the microstructure and mechanical properties of the turbulent zone in friction stir welded 7075 aluminum alloy medium thickness plate joints

Investigation on the microstructure and mechanical properties of the turbulent zone in friction stir welded 7075 aluminum alloy medium thickness plate joints

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10−3 to 10−1, which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

Hybrid aeroacoustic investigation of turbulent 90° pipe bend flow with source terms from Large-Eddy Simulation

The internal, low Mach number, turbulent flow-induced acoustic field of a cylindrical 90° pipe bend with a curvature radius Rm/D=1.02 is investigated using hybrid Computational Aero-Acoustics (CAA). Incompressible Large-Eddy Simulation (LES) is used to computationally describe the turbulent flow field at bulk Reynolds numbers ReB=5000,12000,26000 and 66000. The incompressible LES predicts all salient flow features, like the separation from the inner wall of the bend followed by a turbulent recirculation zone, in good agreement with literature. The internal flow-induced noise is predicted by Acoustic Propagation Simulations (APS) applying Lighthill’s Equation (LE), Ribner’s Dilatation Equation (RDE) and the Perturbed Convective Wave Equation (PCWE), provided with acoustic source terms obtained from the incompressible LES. The turbulent recirculation zone is the most active region for generating acoustic sources. Among the two constituent components of the PCWE source term, the local pressure time derivative term turned out as clearly dominant in this region, so that PCWE and RDE predicted similar acoustic fields. The sub-grid scale (SGS) model applied in LES notably contributes only at high frequencies to the total LE source term, starting from f=2000 Hz for the lower Reynolds numbers and beyond for the higher. The mapping of source terms from the boundary layer resolving LES grid onto the typically coarser preferably uniform APS grid leads to lower high-frequency signal amplitudes for coarser APS mesh size in regions with small turbulent structures. The impact of this amplitude reduction on the predicted acoustics sound power levels was still very small, exhibiting relative differences less than 0.7% for the different APS meshes. Despite the LES-specific limitations expected from the decrease of directly resolved content for higher Reynolds numbers, the use of increasingly underresolved source terms did not significantly lower the accuracy of finally predicted acoustics. The comparison of the acoustic results against data from dedicated measurements generally showed very good agreement for all Reynolds numbers.

PIV investigation on the effect of gas distributor design for fluidization of Geldart A spherical particles

This is a hydrodynamics study on gas distribution during fluidization of the industrially relevant Geldart A type particles. The impact of 10 different distributor designs comprising perforated distributors (pitch and orifice arrangement pattern varied) and nozzle distributors (nozzle size varied) were studied experimentally on a bench-scale unit. Particle flow homogeneity in air was determined by a novel Particle Image Velocimetry (PIV) analysis of the freeboard velocity field. Over 50 PIV particle phase vector field images are presented and analysed. Bed expansion, Particle attrition with and without fines, velocity profiles with and without fines and effect of fines on bed expansion were ascertained. Best distributor design for a unit such as Fluid Catalytic Cracking (FCC) regenerator and FCC riser were found to be 0.8 cm and 0.6 cm nozzle plates distributors respectively. Confirmation with previous study revealed a correlation between homogeneous distribution of particulate flow and lower number of turbulence zones.

Experimental study and turbulence dissipative scale modelling of the rapid micromixing of impinging thin sheets of liquids

Previous studies have shown that the impingement of thin liquid sheets produces high energy dissipation rates due to the release of kinetic energy in a very small mass of liquid (0.0001–0.01 g), even though flowrates are on the order of L/min. Rapid micromixing occurs because the dissipated energy leads to a substantial reduction in the initial segregation size scale of the liquids, which is the single-sheet thickness at impingement (~ 100 μm). In the present study, the micromixing was investigated by following the progress of an acid-base neutralization accompanied by a change in fluorescence intensity of a fluorophore. Micromixing was modeled using a framework that assumes diffusion and reaction of species occur within slabs of fixed thickness (2L). The slabs are fixed in size because the released kinetic energy is dissipated within one turnover time of the large energy-containing eddies produced in the turbulent impingement zone. A simulation, which included a module for calculating the fluorescence intensity, determined 2L for experimental energy dissipation rates (ϵ) ranging from 4×104 to 7.7×106 W/kg. 2L was found to lie in the range of turbulent dissipative scales less than the Taylor microscale but greater than the Kolmogorov microscale. 2L/v′ was found to be a function of ReΛ−1/4 and ν/ϵ1/2, where v′ is the velocity associated with the turbulent kinetic energy, ReΛ is the large-eddy Reynolds number and ν is the kinematic viscosity. Similar to η/Λ, 2L/Λ is a function of ReΛ−3/4, where ηis the Kolmogorov microscale and Λ is the size of the large energy-containing eddies.

Experimental investigation and performance prediction of SAH using different arc rib roughness geometries – A comparative study

This study aims to investigate and compare the thermal performance of a solar air heater using a passive technique to enhance heat transfer between the absorber plate and the flowing fluid. The technique involves generating turbulence near the heat transferring surface through the use of artificial rib roughness. The study focuses on two different novel roughness geome-tries: full symmetrical arc rib roughness and half symmetrical arc rib roughness. By introducing additional gaps and varying the number of gaps in the roughness geometries, the study examines their effects on the solar air heaters thermal performance. The artificially roughened surface creates different turbulent zones, which are essential to the development of different types of turbulence in the vicinity of the heat transferring surface. The study finds that an optimal escalation in Nusselt number and friction factor by 2.36 and 3.45 times, respectively, occurs at certain gap numbers as 6 and ng as 5 for full symmetrical arc rib roughness. The maximum thermal-hydraulic performance parameter of 1.66 is attained at a Reynolds number of 6 000. The study also conducts correlation, mathematical modeling, and performance prediction under different operating circumstances.