Upwind Discretization Research Articles

Heat transfer analysis of CuO-ZnO/water hybrid nanofluid in a Shell and tube heat exchanger with various tube shapes

Heat exchangers are widely used in various applications, including heating and cooling systems. Adding nanoparticles to the base working fluid can enhance the fluid's thermophysical properties, thereby increasing the thermal performance of heat exchangers. Further heat transfer enhancement can be done by using a hybrid nanofluid, which contains two or more different nanoparticles. The findings of the paper can be utilized to enhance the thermal efficiency of a shell and tube heat exchanger with different tube shapes and by using CuO-ZnO/ water hybrid nanofluid as a working fluid. Numerical simulations were conducted to model thermal heat exchangers with different configurations. Copper tubes of circular, hexagonal, and elliptical geometry were considered for the experiment. The model is solved using Ansys software with a finite volume approach. The solver employs an upwind-based multidimensional linear approach and upwind discretization schemes. The water inlet velocity in the Shell ranges from 0.5m/sec to 3.2 m/sec at a Reynold number of 10,000 to 15,000 whereas the velocity of cold fluid in the tube varies from 1.4 to 2 m/sec. All experiments were conducted with a 0.01–0.1 %vol % CuO-ZnO (80:20)/ water hybrid nanofluid. The hybrid nanofluid on the tube side is used to cool hot fluid (water) on the shell side. Hybrid nanofluid inlet temperature was maintained at 30°C to cool hot water at 60°C. Hexagonal tube geometry increases the contact surface area, resulting in a 22.11 % increase in heat transfer rate compared to round tubes and an 18.42 % increase compared to round tubes. Experimental and numerical results for the convective heat transfer coefficient, Nusselt number, and pressure drop were also reported.

Impact of mixed boundary conditions and nonsmooth data on layer‐originated nonpremixed combustion problems: Higher‐order convergence analysis

AbstractThis work explores the theoretical and computational impacts of mixed‐type flux conditions and nonsmooth data on boundary/interior layer‐originated singularly perturbed semilinear reaction–diffusion problems. Such problems are prevalent in nonpremixed combustion models and catalytic reaction models. The inclusion of arbitrarily small diffusion terms results in boundary layers influenced by specific flux conditions normalized by perturbation parameters. We have demonstrated theoretically that the sharpness of the boundary layer is significantly reduced when this normalization is independent of the diffusion parameter. In addition, the presence of a nonsmooth source function gives rise to an interior layer in the current problem. We show that using upwind discretizations for mixed boundary fluxes achieves nearly second‐order accuracy if the first derivatives are not normalized concerning perturbation parameters. This outcome arises from the bounds of a prior derivative of the continuous solution. Furthermore, it is theoretically shown that nearly second‐order uniform accuracy can be attained for reaction‐dominated semilinear problems using a hybrid scheme at the discontinuity point. To ensure the uniform stability of the discrete solution, a transformation is necessary for the corresponding discrete problem. Theoretical results are supported by various experiments on nonlinear problems, illustrating the pointwise rates and highlighting both linear and higher‐order accuracy at specific points.

Comparison of heat transfer characteristics of a heat exchanger with straight helical tube and a heat exchanger with coiled flow reverser

Vertical shell and coiled tube heat exchangers are an important aspect of numerous industrial processes and are renowned for their high heat transfer efficiency, versatility, compact design, easy manufacturability and maintenance. As a novel technique, a coiled flow reverser was used to enhance the thermal performance of shell and coiled tube heat exchangers. Therefore, the impact of employing a coiled flow reverser on the heat transfer characteristics and pressure drop in a counter-current shell and tube heat exchanger is investigated using both numerical and experimental approaches. The results are compared to a straight helical coiled tube heat exchanger. The numerical simulations were carried out by employing computational fluid dynamics (CFD) techniques, and the SIMPLE scheme was implemented through the second-order upwind discretization method. The laboratory device used for experimentation features a shell with an inner diameter of 15 cm, and the coil pitch and coil diameter of both coiled tubes are set at 3 cm and 12.1 cm, respectively. Results reveal that the overall heat transfer coefficient of the heat exchanger with a coiled flow reverser is 39.5 % to 49.9 % higher than that of the heat exchanger with a straight helical tube. The effectiveness ratio of the heat exchanger with coiled flow reverser to the straight helical tube heat exchanger ranges from 1.254 to 1.379. Optimal results are observed at a cold flow rate of 6 LPM and a hot flow rate of 1 LPM for both heat exchangers. Additionally, the coiled flow reverser enhances heat transfer performance at the cost of increased pressure drop. Numerical results are in good agreement with experimental findings, establishing the reliability of the study.

The impact of wedge and wavy-wall combustor side walls on shock-induced mixing and combustion in a hydrogen-fuelled model scramjet combustor

The mixing and combustion processes are linked and challenging for a supersonic combustion engine. Supersonic air doesn't spend much time in the combustor, making it challenging to blend fuel with supersonic air. The expansion of shear-mixing layer with shock waves and their interactions are the primary characteristics of mixing augmentation. This paper investigates the same methodology for a model scramjet combustion chamber with wedge and wavy sides. The novel scramjet combustion models include wedge and way-wall sidewalls based on the basic scramjet model. By combining the finite volume approach with a second-order upwind discretization method, the Reynolds-averaged Navier-Stokes equations are numerically analyzed. The interplay between wedge and wavy-wall effects and free stream flows is taken into consideration by using the shear stress transport (SST) k−ω turbulence model. The effect of the wedge and wavy wall combustor has been evaluated by analysing the shock production and their interactions with the mixing layer and internal flow structure. The turbulence characteristics are also examined to analyse the degree of flow disturbance for wedge and wavy-wall combustion. The wavy wall combustor model results in a 14.5% increase in turbulence intensity, which raises air-fuel mixing and there by improves mixing and combustion efficiency by 12.59% and 11.23%, compared to the base wedge type combustor.

Quasi-Newton positivity-preserving scheme for the Spalart-Allmaras turbulence model using unstructured grids

A new implicit method for solving the Spalart-Allmaras (SA) turbulence model is proposed using a formal second-order upwind scheme on unstructured grids for steady-state flows. The mean-flow set of equations and the SA model equation are solved in a loosely coupled manner. While the Jacobian of the mean-flow equations is derived from a low-order approximation of the residual (i.e. the defect-correction approach), the SA Jacobian is obtained by modifying a direct linearization of the second-order residual. The modified Jacobian is designed to form an M-matrix without the artificial time derivative commonly used for steady-state problems. Namely, no time step is involved in solving the SA model equation. Thanks to a special design of the M-matrix, the positivity of the SA model working variable and the convergence of the SA model linearized problem (fixed-point iterative problem) is guaranteed. To further enhance the overall stability of the flow solver, the Runge-Kutta implicit algorithm is employed, allowing CFL numbers as high as one to two thousand for the mean-flow equations. A second-order upwind scheme is achieved with a least-square method using a limiter to suppress oscillation in the solution. However, second-order upwind scheme discretization of turbulence models may lack robustness, especially when using unstructured grids. The limiter can hinder the convergence of the non-linear iterations, especially of the turbulence model. A square root transformation is applied to the dependent variable of a baseline SA model to improve the non-linear convergence characteristics. Three variants of the SA model are implemented and tested together with two limiters. Numerical simulations of a broad range of aerodynamic applications are conducted. The transformed variant exhibits the least sensitivity to the limiter type.

The effect of alphabet-shaped twisted tape on thermo-fluidic properties in a shell-tube type heat exchanger: A comparative study utilizing a two-phase mixture model

In the present Computational Fluid Dynamic (CFD) simulation, the impact of various shapes of alphabet twisted tape in a shell –tube type of heat exchanger on the pressure drop, temperature, velocity distribution, and heat transfer coefficient was studied numerically. The best shape of the twisted tape was selected by thermo-hydraulic analysis. SIMPLE algorithm and Finite Volume Method (FVM) were then implemented to solve the conversion formulas. Conversion equations were discretized using second-order upwind discretization. Comparing the results of different types of twisted tape, Case I (W-shaped tape) exhibited the best performance and the utmost pressure differences. Furthermore, the utmost heat transfer factor was obtained by W shape twisted tape. Also, the helical baffles illustrated the best performance of the shell-part baffles in the shell-tube type of heat exchanger. Comparison between single and two-phase flow by utilizing mixture model for heat transfer factor and pressure difference reveal the negligible changes. On the quant part, FOM amount have been shown a great significance (1.3439) in the shell part for helical baffle type.

Parameter uniform higher order numerical treatment for singularly perturbed Robin type parabolic reaction diffusion multiple scale problems with large delay in time

In this paper, we address a class of boundary layer originated singularly perturbed parabolic reaction-diffusion problems with Robin boundary conditions having large time delay; for the higher order numerical analysis. The boundary layer originated nonuniform mesh, is generated by adaptive moving mesh approach based on equidistribution principle. The mesh is obtained by solving a nonlinear system; involving the discrete problem and the discrete nonlinear mesh generating functions. We provide the theoretical mesh structure. The difference scheme for Robin type boundary conditions having normalized flux; uses a combination of space and time discretizations so that the order of convergence in space is higher than the usual linear accurate upwind discretization in space. The present theoretical convergence analysis requires less regularity of the solution and proves to be exactly second order parameter uniform accurate in space. Numerical experiments highly validate the theoretical findings and also show that the present approach is better than the well-known upwind discretizations.

Investigation of the effects of volume change on flow structure and acoustic in a silencer

This study numerically examined the propellant flow from gunfire using the kω-SST turbulence model and their sound attenuation using the Ffowcs-Williams and Hawkins equations (FW-H). For simulation, a pressure-based solver and 3D axisymmetric geometry were used. The second-order implicit time approach and the second-order upwind scheme spatial discretization were used in the simulation. The maximum exit pressure was 3.748 MPa for the suppressor with a length of 70 mm and diameter of 20 mm. However, when the diameter suppressor increased by 1/6, the maximum exit pressure was reduced to 3.4961 MPa. When the length increased by 1/6, the maximum pressure became 3.3636 MPa. Lastly, when the diameter and length were increased by 1/6, the maximum exit pressure became 3.177 MPa. For this suppressor, 20.835 dB (12.29%) sound pressure level attenuation was achieved with 16.823 MPa (84.115%) overpressure reduction and 484.86 K or 32.32% temperature reduction. Generally, the attenuation increased with the increase in the suppressor’s internal volume.

Effect of jet angle and jet pitch on the Thermo-Hydraulic performance of solar air heater having absorber plate with jet impingement

This numerical investigation focuses on enhancing the heat transfer characteristics and thermo-hydraulic performance of a solar air heater by employing jet impingement over the absorber plate. The study explores the influence of two geometrical parameters, namely jet pitch ratio (P/Dh) and jet angle ratio (α/90), on the heat transfer, frictional losses, and overall performance of the solar air heater with jet impingement (SAHJI) under various flow conditions, represented by Reynolds numbers (Re) ranging from 3500 to 17500. The simulations are conducted using the RNG k-ɛ model with second-order upwind discretization. A comparative analysis is performed between the SAHJI system and a conventional smooth solar air heater (SPSAH) configuration, revealing a significant enhancement in thermal–hydraulic performance with the use of jet impingement. The investigation identifies the maximum thermo-hydraulic performance parameter, reaching a value of 4.12. This peak performance is achieved at specific values of jet pitch ratio (P/Dh) of 1.30 and jet angle ratio (α/90) of 0.094, respectively, at a Reynolds number of 7500. These findings demonstrate the potential of jet impingement as a heat transfer enhancement technique in solar air heaters, contributing to improved overall system efficiency.

A Discontinuous Galerkin and Semismooth Newton Approach for the Numerical Solution of Bingham Flow with Variable Density

Abstract This paper is devoted to the study of Bingham flow with variable density. We propose a local bi-viscosity regularization of the stress tensor based on a Huber smoothing step. Next, our computational approach is based on a second-order, divergence-conforming discretization of the Huber regularized Bingham constitutive equations, coupled with a discontinuous Galerkin scheme for the mass density. We take advantage of the properties of divergence-conforming and discontinuous Galerkin formulations to effectively incorporate upwind discretizations, thereby ensuring the stability of the formulation. The stability of the continuous problem and the fully discrete scheme are analyzed. Further, a semismooth Newton method is proposed for solving the obtained fully discretized system of equations at each time step. Finally, several numerical examples that illustrate the main features of the problem and the properties of the numerical scheme are presented.

Exact Nonlinear Spectral Analysis of Nonsmooth WENO-JS Solutions

The fifth-order accurate Weighted Essentially Nonoscillatory space discretization developed by Jiang and Shu (WENO-JS) is studied theoretically. An exact Nonlinear Spectral Method (NSM) is developed based on an innovative yet simple methodology. The NSM explains the behaviour of nonsmooth solutions because it is valid for arbitrary modified wave numbers (MWN). The NSM clarifies the effects of the time integration methods and the Courant number. The mode isolation assumption, extensively used to analyse WENO-JS, is elucidated, and several novel findings are presented. The improved performance of the combination of WENO-JS with the forward Euler time integration method, compared to the Linear Fifth-Order Upwind discretization, is illustrated. The overdamping of the combination of WENO-JS with the popular third-order total variation diminishing Runge-Kutta method is discovered. Thus, the NSM covers several gaps in the current literature.

Optimisation of the Aerodynamic Shape of a Monorail Suspended Unibus

The aerodynamic optimisation of the shape of a monorail suspended unibus of string transport has been described within the study on the influence of geometric and structural elements on aerodynamic characteristics. The estimate was carried out through a comparative analysis of the features of two shapes of a model with a further change and recalculation of the model being finalised. The comparison was focused on the drag force and the shape drag coefficient. The calculations used a gas dynamics model based on the Reynolds equations using the Menter SST-k-ω (shear stress transport) turbulence model. To solve the equations to find the design functions, an upwind second-order discretisation scheme was used applying the «pressure–velocity» refinement procedure in the framework of SIMPLE algorithm of Patankar–Spalding with the ANSYS Computational Fluid Dynamics Software. The dimensions of the computational domain were chosen considering the geometric dimensions of the 3D model of the shape. Boundary conditions were identified in the solver. The simulation was carried out for the case of motion of a vehicle at a constant speed. The calculations have shown the importance and influence of the geometry of the transition sections of the vehicle body, the mandatory use of wheel fairings and the advantages of the S-shaped tail. The proposed design optimisation made it possible to reduce the drag force and coefficient by 16,9 %. The studies have resulted in selection of the optimal vehicle model which has the lowest aerodynamic drag coefficient, which made it possible to improve the energy efficiency of the system and its environmental friendliness, and consequently, the profit potential of the transportation process.

Computational analysis and evaluation of thermal comfort environment of a space by reducing the inside air temperature

The tents are used as temporary shelters and evenly they sometimes become a home for a number of years when the returns to permanent shelter is very slow or unaffordable for any person so that the comfort of these shelters are very much essential like home. This paper is primarily concerned with identifying the objectives of reducing the inside air temperature of shelter using different retro-reflective materials. In the present work computational fluid dynamics analysis have been performed for a shelter under the summer climatic conditions to reduce the inside air temperature of shelter using different retro-reflective materials like polyvinyl chloride, wood-fibre, polystyrene, aerogel and phenolic foam. For CFD analysis free stream temperature of air is set 311 K taken with constant heat flux of 900 W/m2 solar radiations is used on roof, second order upwind discretization scheme is used to approximate the solution for momentum and pressure. Result show that that minimum temperature variation of 4.8 degreesK by using retro-reflective material as phenolic-foam and 18.74 % temperature percentage variation compared with other reflective materials.

Free surface dynamics of MHD third‐grade fluid model over a heated stretching sheet with variable fluid properties

AbstractA thorough investigation of MHD third‐grade differential‐type fluid flow over a heated stretching sheet is performed in this work. In particular, we analyze the film thinning process, when the thermal sensitive fluid parameters vary due to the effect of heat supplied to the stretching sheet. Starting with a two‐dimensional (2D) free surface boundary value problem of non‐Newtonian third‐grade fluid, we present a systematic derivation of a 1D transient thin‐film height equation using longwave analysis with respect to the small aspect ratio of the fluid domain. The derived model is used to study the impact of Newtonian and non‐Newtonian parameters with variable fluid properties on the thin film height. The model is discretized using an upwind discretization in space and implicit time integration to guarantee first‐order convergence. The model is analyzed thoroughly with the help of numeric computing software MATLAB. The existing findings for a Newtonian fluid are in excellent agreement with derived evidence. In comparison to Newtonian fluid, the study finds that the third‐grade parameter causes thinning under different parametric restrictions. Simulations on the coupling effect explain that, the film thickness can be reduced with a high Marangoni number for highly viscous fluids. Also, since the effect of the conductivity parameter can be reduced at a low Prandtl number, the fluid shows a thinning effect. The film thinning rate, on the other hand, is reduced by the magnetic field.

Longwave modeling of thin film flow of a generalized second-grade fluid down a slanted plate

The free surface dynamics of a thin film of a generalized second-grade fluid flowing down a slanted plate subjected to the action of gravity have been studied. A nonlinear evolution equation for the dynamics of the thin liquid film is derived using the long-wavelength approximation. The model allows to investigate the impact of fluid rheology and geometrical parameters on the thin film’s height. The nonlinear dynamics equation is implicitly approximated on a uniform grid by applying the finite volume discretization method, which uses upwind discretization of the flux function and first-order discretization of the time differential term. The investigation reveals that the free surface deformation trends are similar for the Newtonian and non-Newtonian fluids but the non-Newtonian fluid properties substantially influenced the size and shape of the deformation. It is also found that the shear-thinning fluid moves faster down to the plate compared to the Newtonian and shear-thickening fluid.

Entropy generation analysis of a dual-height plate-fin heat sink with Nanofluid coolant applying two-phase mixture model

In this work, the two-phase mixture model is applied to solve the velocity field, the pressure field and obtain the entropy generation rates in a dual-height plate-fin heatsink. The governing equations were solved using the second-order upwind discretization scheme as well as SIMPLE scheme used in coupling the pressure and velocity fields. The Re, nanoparticle concentration (ψ), secondary fin height (hsf), and fin spacing (dsf) were considered within the ranges of 500-2000, 0%-1%, 5 mm-17.5 mm, and 0.25 mm-0.5 mm, respectively. The results demonstrated that S˙th almost decrease by 30% as hsf increases from 5 mm to 17.5 mm. The increase in hsf escalates S˙fr by 308% and 238% at Re numbers of 500 and 2000, respectively. In addition, the escalation of Re from 500 to 2000 increases S˙fr by 2125% and 1748% for hsf 5 mm and 17.5 mm, respectively. The fin spacing analysis revealed that the lowest S˙th and highest S˙fr belong to dsf=0.5 mm. The lower the dsf the higher the backward flow and mixing flow at the heatsink outlet, thereby delay in flow passage and intensification of S˙th and S˙fr.

Numerically Detection Fluid Characteristic Effects in Porous Media for Plastic Manufacturing Process Reconstruction

Research work aimed to prove the cooling effectiveness and efficiency of heating numerically. It was carried out by reconstructing the mold to reduce condensation while keeping the production cycle time. The simulation was done using CFD Fluent, Finite Volume Method (FVM) with a solution method of pressure-velocity, SIMPLE coupling, and second-order upwind discretization scheme. The model's boundary condition represents plastic packaging's production process. The simulation was done to optimize performance during production. Reconstruction and simulation were undertaken. The initial investigation was to produce data for cooling time versus the whole production time. The research result indicated that the average cooling time was 3/5 of the total production time. There were 15.5706 of 26.4206 sec and 15.641 of 26.491 sec, respectively. Further investigation was focused on investigating heat absorption and transfer through cells in the matrix. The porous effect showed temperature gradient reduction between cooled and cooler temperatures by simulation in color. The phenomenon predicted that the Ra contributed to the rate of temperature decrease revealed. There is also an apparent increase in temperature interfaces is evidence. The porous media application with The Rayleigh Number ≥ 105 retained the potential cooling effectivity and heating efficiency. By simulation, a series of discussions of the result of research, the potential may apply in the plastic mold to reconstruct the space and its methods using filling-porous. It is relevant because the heating ability in the phase change of fluid that creates particular characteristics has the potential to make effective cooling.

Effect of Magnetic Baffles and Magnetic Nanofluid on Thermo-Hydraulic Characteristics of Dimple Mini Channel for Thermal Energy Applications

The combined effect of a magnetic baffle and a dimple turbulator on the heat transfer and pressure drop is investigated computationally in a mini channel. Fe3O4 magnetic nanofluid is used as a working fluid. The Reynolds number (Re) is varied from 150 to 210 and the magnetic field intensities range from 1200 G to 2000 G. Finite-volume based commercial computational fluid dynamics (CFD) solver ANSYS-Fluent 18.1 was used for the numerical simulations. A laminar viscous model is used with pressure-velocity coupling along with second-order upwind discretization and QUICK scheme for discretizing the momentum and energy equations. The results show that there is an increase of 3.53%, 10.77%, and 25.39% in the Nusselt numbers when the magnetic fields of 1200 G, 1500 G and 2000 G, respectively, are applied at x = 15 mm, as compared to the flow without a magnetic field when the pitch = 10 mm. These values change to 1.51%, 6.14% and 18.47% for a pitch = 5 mm and 0.85%, 4.33%, and 15.25% for a pitch = 2.5 mm, when compared to the flow without a magnetic field in the respective geometries. When the two sources are placed at x = 7.5 mm and 15 mm, there is an increase of 4.52%, 13.93%, and 33.08% in the Nusselt numbers when magnetic fields of 1200 G, 1500 G, and 2000 G are applied when the pitch = 10 mm. The increment changed to 1.82%, 8.16%, and 22.31% for a pitch = 5 mm and 1.01%, 5.96%, and 21.38% for a pitch = 2.5 mm. This clearly shows that the two sources at the front have a higher increment in the Nusselt numbers compared to one source, due to higher turbulence. In addition, there is a decrease in the pressure drop of 10.82%, 16.778%, and 26.75% when magnetic fields of 1200 G, 1500 G, and 2000 G, respectively, are applied at x = 15 mm, as compared to flow without magnetic field when the pitch = 10 mm. These values change to 2.46%, 4.98%, and 8.54% for a pitch = 5 mm and 1.62%, 3.52%, and 4.78% for a pitch = 2.5 mm, when compared to flow without magnetic field in the respective geometries. When two sources are placed at x = 7.5 mm and 15 mm, there is an decrease of 19.02%, 31.3%, and 50.34% in the pressure drop when the magnetic fields of 1200 G, 1500 G and 2000 G are applied when the pitch = 10 mm. These values change to 4.18%, 9.52%, and 16.52% for a pitch = 5 mm and 3.08%, 6.88%, and 14.88% for a pitch = 2.5 mm. Hence, with the increase in the magnetic field, there is a decrease in pressure drop for both the cases and the pitches. This trend is valid only at lower magnetic field strength, because the decrease in the pressure drop dominates over the increase in pressure drop due to turbulence.

Novel discretization strategies for the 2D non-Newtonian resistance term in geophysical shallow flows

In the context of two-dimensional models for complex geophysical surface flows such as debris flows, muddy slurries, oil spills over land, hyperconcentrated floods, lava flows, etc, depth-averaged rheological models relate the shear stress state within the fluid column to the depth-averaged local flow features. Despite it is the most influencing term on the mobility of complex shallow flows, the numerical treatment of the resistance contribution to the flow momentum is still a challenging topic, especially when dealing with 2D large-scale applications. In this work, two novel strategies for the explicit upwind discretization of generalized non-Newtonian resistance terms in two-dimensional numerical models are proposed, called integral and differential approaches. These new strategies are applicable to generalized rheological formulations in any type of mesh topology. Results from benchmark tests running in orthogonal, triangle structured and triangle unstructured meshes demonstrate that both approaches represent an improvement for the explicit upwind integration of the 2D resistance force compared with previous procedures. It is shown that the alignment of the flow with the mesh main-axis, which has been previously attributed to faults of 2D FV numerical methods and insufficient mesh refinements, is directly related to the loss of the rotational invariance of the integrated resistance force. This is caused by the erroneous procedure for including the 2D resistance term into the local flux balance at the cell edges. Furthermore, a novel implicit centered method for the integration of the 2D resistance force has also been derived for the quadratic frictional non-linear resistance formulation. Despite the implicit procedure fails to converge to steady uniform flow states, the differential explicit upwind and the implicit centered methods show similar level of accuracy, robustness and computational efficiency for transient 2D frictional visco-plastic flows.

Hydrogen fueled scramjet combustor with a wavy-wall double strut fuel injector

The demand for the enhancement of scramjet performance is increasing worldwide due to its tremendous and upcoming applications. The major challenge in scramjet engines is the mixing of fuel and air and succeeding combustion at supersonic speed. In a scramjet combustor flow field, the development of high turbulence (flow disturbance), shear-mixing layer, and oblique shock waves are the key factors to magnify the mixing efficiency of fuel and air at the significantly less resident time of air in the combustion chamber. Based on these techniques, an innovative strut fuel injector with a wavy wall surface and double strut has been investigated for possible mixing enhancement of fuel and air. These results are compared with the basic single wedge strut available in the open literature to identify the performance improvement. The computational geometries are modeled with Reynolds-averaged Navier Stokes equations, SST k-ω model, and global one-step reaction mechanism. All the flow modeling has been solved with finite volume and a second ordered upwind discretization scheme. The evaluation of scramjet internal flow structure and numerical results found that the double strut with a wavy-wall shape improved the mixing efficiency and combustion phenomenon significantly. The mixing and combustion profiles reveal that wavy-wall double strut generates an early complete mixing within a distance of 0.17 m with increased mixing and combustion efficiency of 18% and 20%, respectively.