High Reynolds Number Conditions Research Articles

Macroscopic turbulent flow via hard sphere potential

In recent works, we proposed a hypothesis that the turbulence in gases could be produced by particles interacting via a potential and examined the proposed mechanics of turbulence formation in a simple model of two particles for a variety of different potentials. In this work, we use the same hypothesis to develop new fluid mechanics equations that model turbulent gas flow on a macroscopic scale. The main difference between our approach and the conventional formalism is that we avoid replacing the potential interaction between particles with the Boltzmann collision integral. Due to this difference, the velocity moment closure, which we implement for the shear stress and heat flux, relies upon the high Reynolds number condition rather than the Newton law of viscosity and the Fourier law of heat conduction. The resulting system of equations of fluid mechanics differs considerably from the standard Euler and Navier–Stokes equations. A numerical simulation of our system shows that a laminar Bernoulli jet of an argon-like hard sphere gas in a straight pipe rapidly becomes a turbulent flow. The time-averaged Fourier spectra of the kinetic energy of this flow exhibit Kolmogorov’s negative five-thirds power decay rate.

Using a Modified Turian-Yuan Model to Enhance Heterogeneous Resistance in Municipal Sludge Transportation Pipeline.

Based on the Turian–Yuan heterogeneous resistance model, the simulation results of three urban sludge pipelines with a volumetric concentration of 2.38, 3.94, and 5.39% were analyzed. The reasons for the large deviation of the simulation results under high Reynolds number conditions were also analyzed. The results showed that the deviation of the simulation was mainly caused by the difference between the sludge volumetric concentration (CV), the settlement resistance coefficient (CD), and the values of the two parameters in the Turian–Yuan heterogeneous resistance model. Consequently, it was necessary to optimize the index m1 of CV and the index m2 of CD. Taking mean square deviation as the objective function, using Matlab programming, the abovementioned two indexes were optimized by the simulated annealing algorithm. The optimized index m1 of CV was 0.887, and the index m2 of CD was −0.162. Hence, a modified Turian–Yuan heterogeneous resistance model was obtained. The model verified that the minimum value of the regression coefficient, R2, of the simulated value reached 0.9701, proving that, the model can be used to simulate the heterogeneous resistance of urban sludge pipeline transportation.

An adaptive plate at flapping wing’s trailing edge in promoting energy extraction performance

As is known, the maximum thickness position plays a key role in determining the energy extraction properties for flapping wing airfoils. Previous work demonstrated that the energy extraction efficiency might increase when we move the maximum thickness position towards the leading-edge direction. Inspired from this phenomenon, we propose in this work the addition of a flat plate at the tail of the flapping wings. Under two different Reynolds numbers (Re = 1100 and Re = 5×105), the effect of adding a rigid flat plate or a flexible flat plate to the trailing edge of the flapping wing on the flapping wing’s energy efficiency was investigated. The corresponding results were also compared with that of the original flapping wing. Results show that the addition of a rigid flat plate is sufficient for the amelioration of the flapping wing’s energy extraction efficiency in low Reynolds number cases. At high Reynolds number cases, the flapping wing with flexible flat plate trailing edge can maintain a relatively high energy extraction efficiency over a wide range of frequency, while the flapping wing with rigid flat plate trailing edge can only improve the energy extraction efficiency in the low frequency region. Therefore, we propose the application of flapping wings with flexible flat plate trailing edge for increasing the energy extraction efficiency under high Reynolds number conditions.

Low-Speed DSMC Simulations of Hotwire Anemometers at High-Altitude Conditions

Numerical simulations of hotwire anemometers in low-speed, high-altitude conditions have been carried out using the direct simulation Monte Carlo (DSMC) method. Hotwire instruments are commonly used for in-situ turbulence measurements because of their ability to obtain high spatial and temporal resolution data. Fast time responses are achieved by the wires having small diameters (1–5 μm). Hotwire instruments are currently being used to make in-situ measurements of high-altitude turbulence (20–40 km). At these altitudes, hotwires experience Knudsen number values that lie in the transition-regime between slip-flow and free-molecular flow. This article expands the current knowledge of hotwire anemometers by investigating their behavior in the transition-regime. Challenges involved with simulating hotwires at high Knudsen number and low Reynolds number conditions are discussed. The ability of the DSMC method to simulate hotwires from the free-molecular to slip-flow regimes is demonstrated. Dependence of heat transfer on surface accommodation coefficient is explored and discussed. Simulation results of Nusselt number dependence on Reynolds number show good agreement with experimental data. Magnitude discrepancies are attributed to differences between simulation and experimental conditions, while discrepancies in trend are attributed to finite simulation domain size.

Numerical Investigation on Proppant–Water Mixture Transport in Slot under High Reynolds Number Conditions

Water hydraulic fracturing involves pumping low viscosity fluid and proppant mixture into the artificial fracture under a high pumping rate. In that high Reynolds number conditions (HRNCs, Re > 2000), the turbulence effect is one of the key factors affecting proppant transportation and placement. In this paper, a Eulerian multiphase model was used to simulate the proppant particle transport in a parallel slot under HRNCs. Turbulence effects in high pumping rates and frictional stress among the proppant particles were taken into consideration, and the Johnson-Jackson wall boundary conditions were used to describe the particle-wall interaction. The numerical simulation result was validated with laboratory-scale slot experiment results. The simulation results demonstrate that the pattern of the proppant bank is significantly affected by the vortex near the wellbore, and the whole proppant transport process can be divided into four stages under HRNCs. Furthermore, the proppant placement structure and the equilibrium height of proppant dune under HRNCs are comprehensively discussed by a parametrical study, including injection position, velocity, proppant density, concentration, and diameter. As the injection position changes from the lower one to the top one, the unpropped area near the entrance decrease by 7.1 times, and the equilibrium height for the primary dune increase by 5.3%. As the velocity of the slurry jet increases from 2 m/s to 5 m/s (Re = 2000–5000), the vortex becomes stronger, so the non-propped area near the inlet increase by 5.3 times, and the equilibrium height decrease by 5.2%. The change of proppant properties does not significantly change the vortex; however, the equilibrium height is affected by the high-speed flush. Thus, the conventional equilibrium height prediction correlation is not suitable for the HRNCs. Therefore, a modified bi-power law prediction correlation was proposed based on the simulation data, which can be used to accurately predict the equilibrium height of the proppant bank under HRNCs (mean deviation = 3.8%).

A new approach for handling body motion by combining a grid deformation method and an overset grids technique

A new approach for handling body motion by combining the grid deformation and the overset grids methods is proposed. The present approach can completely avoid the computational load of the dynamic overset grids method. The regions where the overset relations are composed move with the body motion, maintaining the geometrical relation to determine for the weight of the overset interpolation. The present approach is applied to the cases with 2D single and multiple cylinders with motion. The scaled-Jacobian before and after the grid deformation is examined. The values of the scaled-Jacobian maintain at least 0.35 after the grid deformation, and the velocity and vorticity contour lines are smoothly connected on the computational grids, thus, the effectiveness of the present approach is revealed. Additionally, the verification using the computational grids with three different resolutions is performed on the case with the single cylinder. The present method which combines the grid deformation and overset grids methods succeeds to reduce the CPU time approximately from 40% to 65% of the time using the dynamic overset grids method. Next, the present approach is applied to the computation of the flow around a ship with an appendage, including ship motion in regular head waves with a high Reynolds number condition. The present method reproduces the impact pressure when points on the hull surface enter the water and the difference in the peak values of the pressure at different positions. The present method succeeds in reducing the CPU time to 40% of the case using the dynamic overset grids method at the case with the ship motions in the regular waves.

Contribution of viscosity to the circulation deposition in the Richtmyer–Meshkov instability

This study focuses on the process of the circulation deposition in the Richtmyer–Meshkov instability (RMI). The growth rate of circulation and its sources are theoretically and numerically studied to reveal the physical mechanism of the viscosity in the circulation deposition process. We derive a predicting model of the circulation rate for RMI. More importantly, all the contributing sources are separately predicted. Particularly, the viscous source, which previously lacked theoretical or numerical investigations, is efficiently predicted. The RMI problems in a large range of initial conditions are simulated with the direct simulation Monte Carlo (DSMC) method to verify our predicting model and further reveal the circulation deposition mechanism. The DSMC simulations provide reliable quantification of the circulation deposition (especially viscous contribution) for RMI due to its molecular nature. Our model predicts the circulation rate, baroclinic and viscous sources accurately for all the cases in comparison with the simulations. A new physical insight into the mechanism of viscosity in RMI is provided. Unlike the previous understandings that nearly all circulation deposition in RMI comes from the baroclinic source, this study reveals the hidden positive contribution of the viscous source, especially for high Mach number conditions (up to 11 % of total circulation rate). For RMI, the large viscosity gradient inside the shock waves plays a crucial role in the circulation deposition even under high Reynolds number conditions. Our study also provides exciting opportunities to further understand the viscous contribution to the vorticity dynamics in the reshocked RMI and shock wave–turbulence interactions.

Multipoint design optimization for a controlled diffusion airfoil compressor cascade

Under high altitude and low Reynolds number conditions, the aerodynamic performance of compressor cascades deteriorates drastically. In this paper, an optimally designed system combining class-shape-transformation method, S1 surface flow solver and whale optimization algorithm was established to achieve for a controlled diffusion airfoil, called MANGHH. The aim of this work is to improve our understanding of the loss mechanism for the original cascade and optimal cascade under different inflow conditions. The study shows that the total pressure loss of the optimal cascade at an angle of attack of −4°, 0°, and 6° decreases by 55.9%, 16.1%, and 16.3%, respectively, compared with the original controlled diffusion airfoil. The range of the available low loss incidence improves significantly. At different incidences, the optimal cascade moves the blade loading forward compared with that of the original controlled diffusion airfoil while reducing the growth rate of the boundary layer thickness, eliminating a wide range of flow separations. The optimal cascade reduces the total pressure loss mainly by reducing trailing edge mixing loss compared with that of the original controlled diffusion airfoil. Under different inlet Mach number conditions, a laminar separation bubble appears on the suction surface of the original controlled diffusion airfoil. As the inlet Mach number increases, the position of the laminar separation bubble moves slightly upstream, while the length and depth of the laminar separation bubble increase. Fortunately, the total pressure loss of the optimal cascade decreases significantly compared with that of the original controlled diffusion airfoil. Under different incoming turbulence intensity conditions, the total pressure loss of the optimal cascade is always lower than that of the original controlled diffusion airfoil. As the incoming turbulence intensity increases, the total pressure loss of the original controlled diffusion airfoil decreases first and then increases. However, the total pressure loss of the optimal cascade increases with increasing incoming turbulence intensity due to the improvement of the turbulence dissipation capacity.

Analysis of heat transfer in grooved plain carbon steel tube for solar applications

A computational fluid dynamics study is conducted in pressure drop and thermal performance of a plain material tube and plain material tube enriched with square, circular and trapezoidal internal grooves. The computational fluid dynamics (CFD) is an authoritative method to evaluate the thermal properties for high Reynolds number condition in water heating using solar energy. A transition in flow from laminar to turbulent condition there prevails an enhancement in Nusselt number (Nu), thermal performance and pressure drop. However, for turbulent flow condition an increase in flow rate and Reynolds number (Re) the value of friction factor decreases for plain material tube and internal grooved tubes. The computational fluid dynamics results indicates that trapezoidal shaped grooving tubes attains the maximum value of drop in pressure, thermal performance and Nusselt number as compared to plain material tube and other grooved tubes pertaining to turbulent flow condition of Re from 5500 to 11500. The output signifies that the highest enhanced value of reduction in pressure attained in trapezoidal, circular and square grooving tubes are 68%, 65% and 63%. The CFD results also indicate that the highest value of increase of heat transfer rate (Q) is attained in trapezoidal grooved tube as 27%. The enhancement in Q value pertaining to circular and square grooved tubes are 21% and 14%.

Heat Transfer in a Rotating, Blade-Shaped Serpentine Cooling Passage With Discrete Ribbed Walls at High Reynolds Numbers

Abstract This paper experimentally investigates the effect of rotation on heat transfer in a typical turbine blade, three-pass, serpentine coolant channel with discrete ribbed walls at high Reynolds numbers. To achieve the high Reynolds number (Re → 190,000) and low rotation number conditions, pressurized Freon R-134a vapor is utilized as the working fluid. Cooling flow in the first passage is radial outward; after the 180 deg tip turn, the flow is radial inward through the second passage; and after the 180 deg hub turn, the flow is radial outward in the third passage. The effects of rotation on the heat transfer coefficients were investigated at rotation numbers as low as 0.07 and Reynolds numbers from 85,000 to 187,000 (based on the first passage geometry and flow conditions). Heat transfer coefficients were measured using thermocouples embedded in copper plates to provide regionally averaged heat transfer coefficients. Heat transfer enhancement due to rotation is observed on the first passage, pressure-side with radially outward flow and the second passage, suction-side with radially inward flow, but a reduction in heat transfer is observed on the third passage pressure-side with radially outward flow. In addition, results from the discrete, broken ribs are compared with those from the same serpentine coolant passage with conventional, angled ribbed walls. A significant increase in the heat transfer due to the discrete ribs is observed in the first passage. These results can be useful for understanding real rotor blade coolant passage heat transfer under high Reynolds number and low rotation number conditions.

CFD 기법을 활용한 공기층에 의한 마찰항력 감소 현상 연구

The flow pattern of air layers and skin-friction drag reduction by air injection are investigated to find the suitable multiphase flow model using unstructured finite-volume CFD solver for the Reynolds-averaged Navier-Stokes equations. In the present computations, two different multiphase flow modeling approaches, such as the Volume of Fluid (VOF) and the Eulerian Multi-Phase (EMP), are adopted to investigate their performances in resolving the two-phase flow pattern and in estimating the frictional drag reduction. First of all, the formation pattern of air layers generated by air injection through a circular opening on the bottom of a flat plate are investigated. These results are then compared with those of M kiharju’s experimental results. Subsequently, the quantitative ratios of skin-friction drag reduction including the behavior of air layers, within turbulent boundary layers in large scale and at high Reynolds number conditions, are investigated under the same conditions as the model test that has been conducted in the US Navy

Fluid-Structure-Interaction Analysis of an Aero Hydraulic Pipe Considering Friction Coupling

A hydraulic pipe is a basic component of the aero hydraulic system, while high pressure and high speed lead to the enhancement of multiple field coupling effects, such as the structure, fluid, heat, and sound of the aero hydraulic pipe. A larger elastic deformation of the pipe and the more common unsteady flow of fluid can aggravate the vibration and noise of the hydraulic pipe system. Based on the 14-equation model of fluid-structure-interaction (FSI) vibration of aero hydraulic pipe, a model suitable for a wide Reynolds number range is proposed by comparing and analyzing various kinds of friction coupling models. This model is used for the FSI dynamic behavior analysis of aero hydraulic pipes, which improves the accuracy of FSI frequency-domain solutions under high Reynolds number conditions. Second, based on the 14-equation model of FSI vibration of aero hydraulic pipe and the transfer matrix method, the mathematical model of airframe deformation load of an aero hydraulic pipe is established, and the evolution law of the FSI dynamic behavior of the aero hydraulic pipe under different airframe deformation loads is analyzed. All the above analysis results were verified by the experimental analysis. The error of the analysis results of the FSI vibration frequency is smaller than 5% compared with the experimental results. The results show that the proposed method is feasible and effective for predicting the hydrodynamic behavior of the aero hydraulic pipes under load conditions.

Settling behavior of non-spherical particles in power-law fluids: Experimental study and model development

Solid-particle settling occurs in many natural and industrial processes, such as in the transportation of drilling cuttings and fracturing proppant. Knowledge of the drag coefficient and settling velocity of cuttings and proppant is of significance to hydraulics design, wellbore cleanout, and fracture optimization. We conducted 553 tests to investigate the settling characteristics of spherical and non-spherical particles in power-law fluids. Three major particle shapes (spherical, cubic, and cylindrical) and eight different particle sphericities were used to simulate cuttings and proppant, and power-law fluids were applied to simulate drilling and fracturing fluids. Based on the data analysis, a new drag coefficient–particle Reynolds number correlation was developed to determine the drag coefficient in a power-law fluid for spherical and non-spherical particles. The drag coefficient increases as the sphericity decreases for the same particle Reynolds number. For a specific particle shape, the drag coefficient decreases as the particle Reynolds number increases, but the decreasing trend is reduced at high particle Reynolds number conditions. An explicit settling-velocity equation was proposed to calculate the settling velocity of spherical and non-spherical particles in power-law fluids by considering the effect of sphericity. A suitable range for the proposed model is 0.0001<Re<200, 0.471<ϕ<1, and 0.505<n<1. An illustrative example is presented to show how to calculate the drag coefficient and settling velocity in power-law fluids with given particle and fluid properties.

Numerical Study of Heat Transfer Enhancement of Internal Flow Using Double-Sided Delta-Winglet Tape Insert

A numerical study was performed to investigate the thermal performance characteristics of an enhanced tube heat exchanger fitted with punched delta-winglet vortex generators. Computational fluid dynamics modeling was applied using the k–ε renormalized group turbulence model. Benchmarking was performed using the results of the experimental study for a similar geometry. Attack angles of 30°, 50°, and 70° were used to investigate the heat transfer and pressure drop characteristics of the enhanced tube. Flow conditions were considered in the turbulent region in the Reynolds number range of 9100 to 17,400. A smooth tube was employed for evaluating the increment in the Nusselt number and the friction factor characteristics of the enhanced tube. The results show that the Nusselt number, friction factor, and thermal performance factor have a similar tendency. The presence of this insert offers a higher thermal performance factor as compared to that obtained with a plain tube. Vortex development in the flow structure aids in generating a vortex flow, which increases convective heat transfer. In addition, as the angle is varied, it is observed that the largest attack angle provides the highest thermal performance factor. The greatest increase in the Nusselt number and friction factor, respectively, was found to be approximately 3.7 and 10 times greater than those of a smooth tube. Through numerical simulations with the present simulation condition, it is revealed that the thermal performance factor approaches the value of 1.1. Moreover, the numerical and experimental values agree well although they tend to be different at high Reynolds number conditions. The numerical and experimental values both show similar trends in the Nusselt number, friction factor, and thermal performance factor.

Effects of heat release and fuel type on highly turbulent premixed jet flames

An experimental study on the fuel type and attendant heat release effects on piloted turbulent premixed jet flames was carried out. The investigation focused on four fuels, namely methane, ethylene, n-heptane, and toluene, five lean to stoichiometric equivalence ratios, and two jet Reynolds numbers. The average flame height and the global turbulent consumption speed were scaled against the laminar flame speed, the maximum heat release rate, and the laminar flame thickness for all fuels and conditions. Results showed that for the different fuel types and for the lower equivalence ratios, the average flame height does not scale well with any of the aforementioned parameters, while the global consumption speed was determined to scale well with laminar flame speed and maximum heat release rate. The thickness of the shear layer was also characterized through detailed particle image velocimetry measurements and was found to decrease with increasing heat release when the jet diameter is used to scale the axial distance. However, when a density-based momentum diameter is used as the scaling parameter of the axial distance, the effect of heat release is suppressed. Additionally, the growth rates of the shear layer thickness could not be scaled between the different fuels using either the laminar flame speed or the maximum heat release rate. Finally, the turbulent kinetic energy and turbulent shear stress development in the shear layer was determined also to be different between the fuels despite keeping either the laminar flame speed or maximum heat release rate constant, with methane flames most notably showing a sharper decay in intensity with axial distance. The fuel effects on the turbulent kinetic energy and the turbulent shear stress were more pronounced at the higher Reynolds number conditions.

Wavepacket models for subsonic twin jets using 3D parabolized stability equations

An extension of the classical parabolized stability equations to flows strongly dependent on the two cross-stream spatial directions and weakly dependent on the streamwise one is applied to model the large-scale structures present in twin-jet configurations. The existence of these unsteady flow structures, usually referred to as wavepackets, has been demonstrated in the literature for both subsonic and supersonic round jets, along with their relation to the generation of highly directional noise emitted in the aft direction. The present study considers twin-jet configurations with different separations at high Reynolds number and subsonic conditions. The existing instability modes for the twin-jet mean flow, their dependence on the separation of the two jets, and the interaction between the wavepackets originating from the two jets is investigated here. Arising from the axisymmetric mode for single round jets, two dominant modes are found for twin jets: a varicose one, relatively insensitive to jets' proximity, but likely to be efficient in radiating noise; a sinuous one, whose amplification is strongly dependent on the jets' distance, and which can be expected to produce weaker acoustic signatures.

Flow turbulence topology in regular porous media: From macroscopic to microscopic scale with direct numerical simulation

Microscopic analysis of turbulence topology in a regular porous medium is presented with a series of direct numerical simulation. The regular porous media are comprised of square cylinders in a staggered array. Triply periodic boundary conditions enable efficient investigations in a representative elementary volume. Three flow patterns—channel with sudden contraction, impinging surface, and wake—are observed and studied quantitatively in contrast to the qualitative experimental studies reported in the literature. Among these, shear layers in the channel show the highest turbulence intensity due to a favorable pressure gradient and shed due to an adverse pressure gradient downstream. The turbulent energy budget indicates a strong production rate after the flow contraction and a strong dissipation on both shear and impinging walls. Energy spectra and pre-multiplied spectra detect large scale energetic structures in the shear layer and a breakup of scales in the impinging layer. However, these large scale structures break into less energetic small structures at high Reynolds number conditions. This suggests an absence of coherent structures in densely packed porous media at high Reynolds numbers. Anisotropy analysis with a barycentric map shows that the turbulence in porous media is highly isotropic in the macro-scale, which is not the case in the micro-scale. In the end, proper orthogonal decomposition is employed to distinguish the energy-conserving structures. The results support the pore scale prevalence hypothesis. However, energetic coherent structures are observed in the case with sparsely packed porous media.

Design and Optimisation of a Low Reynolds Number Airfoil for Small Horizontal Axis Wind Turbines

Horizontal axis wind turbines (HAWTs) are scaled down to incorporate rotor blades that usually have a diameter ranging from two to four meters in length. A common misconception with regard to selection of airfoils followed by subsequent designing of rotor blades involves the use of NACA airfoils and other conventional high Reynolds number airfoils. Micro horizontal axis wind turbines usually operate at low Reynolds number conditions along the blade length. Conventional aerodynamic schemes cannot be applied to rotor blades operating under low Reynolds number conditions as compared to those occurring under high Reynolds number conditions since certain unusual aerodynamic phenomena predominate in the case of the former. The efficiency of a wind turbine is largely dependent on blade optimization, which is why airfoil selection of the rotor blade is of considerable importance. The difference between high and low Reynolds number operation is the onset of boundary layer transition. In the case of high Reynolds number operation, as is the case of aircraft propellers and other high speed turbines, boundary layer transition takes place before laminar separation, which is in direct contrast to low Reynolds number boundary layer phenomena, wherein laminar separation takes place before boundary layer transition.

A geometric study on shell side heat transfer and flow resistance of a six-start spirally corrugated tube

Heat transfer enhancement is of great importance for energy efficiency improvement. The utilization of spirally corrugated tubes is one of the efficient ways to strengthen heat transfer. In this article, based on a validated numerical model, the effects of geometric parameters of a six-start spirally corrugated tube, including the pitch p and the corrugation depth e, on the shell side heat transfer and flow resistance performance are numerically investigated, in high Reynolds number conditions ranging from 10,000 to 60,000. The shell side secondary flow velocity distribution, longitudinal vortex distribution, and temperature distribution of a six-start spirally corrugated tube are presented, respectively. In addition, the heat transfer and flow resistance characteristics are evaluated by comparing the Nusselt number and the flow resistance coefficient with these of smooth tubes. Results show that the utilization of six-start spirally corrugated tubes can enhance the heat transfer performance at the expense of an increase of the flow resistance. However, with the same geometric parameters, the Nusselt number increases and the flow resistance coefficient decreases as Reynolds number increases. With the pitch increasing, the Nusselt number and the flow resistance coefficient decrease at a fixed Reynolds number. In contrast, as the corrugation depth increasing, the Nusselt number changes irregularly, and the flow resistance coefficient increases. Finally, correlations for the shell side Nusselt number and flow resistance coefficient of the six-start spirally corrugated tube are established. This work is of significance for engineers and scientists focusing on the heat transfer and the flow resistance characteristics of spirally corrugated tubes and their applications.

Calculation of proppant-carrying flow in supercritical carbon dioxide fracturing fluid

As a new clear and efficient fracturing fluid, supercritical carbon dioxide (SC-CO2) has many advantages in the development of unconventional oil and gas resources, such as no harm to the reservoir, improving the reservoir permeability and fluid flow resistance, quickly completely flowing back after fracturing. Based on the principle of suspension energy balance, an equilibrium flow rate model of proppant-carrying flow of SC-CO2 is derived, and calculation errors are 3.32% and 4.54%, respectively, compared with equilibrium flow rates obtained from experiments under different density and particle sizes. By considering the change of fluid internal energy and flow work, filtration characteristics, physical parameters, and SC-CO2–rock adsorption in the SC-CO2 fracturing process, a proppant-carrying flow calculation model of the SC-CO2 fracturing fluid is established based on temperature–pressure field coupling, fracture and filtration zone coupling, and SC-CO2 fluid and proppant coupling. In the experiment, the heights of nine groups of proppant embankments were measured, and the average model prediction error value is 3.43%. SC-CO2 has a similar proppant-carrying capacity to that of clear water under the high Reynolds number condition. Based on the proppant-embankment prediction model of SC-CO2 fracturing, through analyzing proppant-carrying characteristics of the SC-CO2 fracturing fluid, the proppant concentration distribution and proppant embankment section are obtained at different fracturing times. According to the calculation, as the fracturing fluid pumping rate increases, proppant density and proppant particle size will decrease, the viscosity of the SC-CO2 fracturing fluid system will increase, and proppant accumulation height can be reduced over a small range, while proppant accumulation length will experience a greater increase along the fracture length direction. During field operation of SC-CO2 fracturing, the proppant pumping process should be optimized according to the proppant embankment section and proppant concentration should be choose increased in small steps to meet fracture propping requirements.