Direct Numerical Simulation Research Articles

Steric effect induces enhancement of electroconvective flow near electrochemical interfaces.

Electroconvection, occurring near electrochemical interfaces, propels the movement of ions and water, leading to intricate phenomena rooted in the fine interplay between fluid, voltage, and ion. Here, neglecting ionic interactions, by incorporating the steric term into the Poisson-Nernst-Planck-Stokes coupling equation, direct numerical simulations of electroconvective vortex near nanoslot-bulk interfaces are conducted. For the steric effect, the steric number is introduced to discuss the factors and laws affecting the vortex. We illustrate the substantial enhancement of electroconvective vortex due to the steric effect of ions within the nanoslot. Upon increasing the steric number, the cation concentration in the nanoslot is enhanced, resulting in the expansion of the electric double layer (EDL). The EDLs on the walls inside the nanoslot come into contact with each other, causing the EDLs to overlap, consequently increasing the total charge within the EDLs inside the nanoslot. This EDL overlap enhances the charge density of the extended space charge layer, leading to the enhancement of the electroconvective vortex. Further, our scaling analysis, corroborated by direct numerical simulation and existing data, establishes the scaling of slip velocity, jointly regulated by the steric number and voltage difference. By modulating the membrane transport characteristics, the steric effect reduces flow structure size and flux fluctuations, which offers new perspectives for manipulating ion transport and flow instability.

Magnus force reduction in a shear-thinning fluid

This study aims to investigate the impact of fluid shear-thinning on the Magnus forces acting on a rotating cylinder or a sphere immersed in an unbounded flow using direct numerical simulation. The Carreau model is employed to represent the shear-thinning fluid, with the considered Reynolds number (Re) ranging from 0.01 to 100, Carreau number (Cu) from 0 to 100, power-law index (n) from 0.1 to 1, and viscosity ratio (β) from 0.001 to 0.5. The rotation rate (α) is fixed at 6. A characteristic Reynolds number, Rec, based on a viscosity evaluated at the characteristic shear rate, γ˙α=αU∞/2a, is introduced. It is found that, at a constant Rec, compared to that in a Newtonian fluid, the Magnus force exerted on the rotating cylinder or sphere in the shear-thinning fluid is reduced. This reduction results from the difference in viscosity distribution between the upper and lower sides of the cylinder or sphere. Furthermore, our analysis demonstrates that the logarithmic reduction in the Magnus force coefficient can be expressed as a linear combination of the logarithm of the strain rate difference and the logarithm of the shear strain rate sensitive function at two limit states, Cu→0 or Cu→∞. This work may be helpful to deepen the understanding of complex rheological behavior encountered in swirling flow hydrodynamics.

Modulating the electrode pore structure using the magnetic field for reduced local-oxygen transport resistance in polymer electrolyte membrane fuel cell

We utilize a magnetic field to align the pore structure of the PEMFC cathode catalyst layer in the direction of O2 infiltration. This alignment reduces the oxygen diffusion resistance within the catalyst layer, thereby augmenting the performance of the PEMFC. Membrane Electrode Assemblies (MEAs) with magnetically aligned cathodes demonstrate a 50% reduction in O2 transfer resistance compared to conventional MEAs, as verified through electrochemical experiments and it makes a performance enhancement in the I-V curve. Beyond electrochemical experimentation, electrode pore analysis utilizing 3D reconstruction and Lattice Boltzmann Direct Numerical simulation (LB-DNS) for internal flow field analysis within pores, reveals that alignment of the electrode structure in the through-plane direction enhanced pore continuity. This continuity and the resultant minimized proportion of dead pores are identified as the causative factors for the observed performance improvements.

Convection induced by centrifugal and Coriolis buoyancy in a rotating Hele-Shaw reactor

The study of heat and mass transfer in a Hele-Shaw cell rotating around a perpendicular axis has various advanced technological applications. These include the design of microfluidic devices and continuous-flow chemical microreactors, to name a couple. In this setup configuration, the quasi-two-dimensional design allows for recording the density field using optical methods, and the rotation enables control of this field through spatially distributed inertial forces. As is known, in the limit of an infinitely thin layer, the Coriolis force vanishes within a standard mathematical model. However, experimental observations of fluid flow in a rotating Hele-Shaw cell indicate the opposite. In this paper, we show that the correct derivation of the equation of motion under the Hele-Shaw approximation leads to the appearance of a Boussinesq-type term for the Coriolis force. To study the effect of the Coriolis buoyancy, we consider the problem of fluid stability during the internal generation of a transfer component, which can be either the concentration of the dissolved substance or the temperature of the medium. The careful study of system dynamics involves linear stability analysis, weakly nonlinear analysis, and direct numerical simulation. The general properties of the disturbance spectrum are analyzed. The branching of solutions near the first bifurcation is studied using the technique of multiple time scales. A stationary convection is replaced by an oscillatory one under the action of the Coriolis force, as demonstrated by weakly nonlinear analysis. Finally, we investigate the nonlinear dynamics using direct numerical simulation.

Large/small eddy simulations: A high-fidelity method for studying high-Reynolds number turbulent flows

Direct numerical simulations (DNS) are one of the main ab initio tools to study turbulent flows. However, due to their considerable computational cost, DNS are primarily restricted to canonical flows at moderate Reynolds numbers, in which turbulence is isolated from the realistic, large-scale flow dynamics. In contrast, lower fidelity techniques, such as large eddy simulations (LES), are employed for modeling real-life systems. Such approaches rely on closure models that make multiple assumptions, including turbulent equilibrium, small-scale universality, etc., which require prior knowledge of the flow and can be violated. We propose a method, which couples a lower-fidelity, unresolved, time-dependent calculation of an entire system (LES) with an embedded small eddy simulation (SES) that provides a high-fidelity, fully resolved solution in a sub-region of interest of the LES. Such coupling is achieved by continuous replacement of the large SES scales with a low-pass filtered LES velocity field. The method is formulated in physical space, with no assumptions of equilibrium, small-scale structure, and boundary conditions. A priori tests of both steady and unsteady homogeneous, isotropic turbulences are used to demonstrate the method's accuracy in recovering turbulence properties, including spectra, probability density functions of the intermittent quantities, and sub-grid dissipation. Finally, SES is compared with two alternative approaches: one embedding a high-resolution region through static mesh refinement and a generalization of the traditional volumetric spectral forcing. Unlike these methods, SES is shown to achieve DNS-level accuracy at a fraction of the cost of the full DNS, thus opening the possibility to study high-Re flows.

Three-dimensional wake transition of rectangular cylinders and temporal prediction of flow patterns based on a machine learning algorithm

This study investigates the three-dimensional (3D) wake transition in unconfined flows over rectangular cylinders using direct numerical simulation (DNS). Two different cross-sectional aspect ratios (AR) and Reynolds numbers (Re) are scrutinized: AR = 0.5 at Re = 200 and AR = 3 at Re = 600. The investigation focuses on characterizing the flow patterns and forecasting their temporal evolution utilizing the proper orthogonal decomposition (POD) technique coupled with a long short-term memory (LSTM) network. The DNS results reveal the emergence of an ordered mode A for AR = 3, attributed to the stabilizing effect of the elongated AR. On the other hand, the case with a smaller AR (= 0.5) exhibits a mode-swapping regime characterized by modes A and B's distinct and simultaneous manifestation. The spanwise wavelengths of mode A and mode B are approximately 4.7 and 1.2 D for AR = 0.5, while the spanwise wavelength of mode A is 3.5 D for AR = 3. The POD serves as a dimensionality reduction technique, and LSTM facilitates temporal prediction. This algorithm demonstrates satisfactory performance in predicting the flow patterns, including the instabilities of modes A and B, across both transverse and spanwise directions. The employed algorithm adeptly predicts the pressure time series surrounding the cylinders. The duration for training the algorithm is only about 0.5% of the time required for DNS computations. This research, for the first time, demonstrates the effectiveness of the POD–LSTM algorithm in predicting complex 3D instantaneous wake transition patterns for flow past rectangular cylinders.

Investigation of porous surface on second mode in hypersonic boundary layer

The second mode is the main instability factor in hypersonic flow and shows the characteristics of acoustic disturbance. The stabilization effect of a porous surface on the second mode in a hypersonic boundary layer is studied according to the acoustic disturbance characteristics—including the flow characteristics of the porous surface, the effects of porous parameters, and the effect of mechanics—by direct numerical simulation. A rectangle is the basic structure of pores (slots) with a two-dimensional porous surface. The results show that the porosity (φ) and the number of pores (n) are the primary and secondary primary factors affecting the inhibition effect on the second mode by the porous surface. By influencing the amplitude and phase distribution of the pressure disturbance, the porous surface affects the wall normal transport between the mainstream and the porous surface in the near wall region, and it finally inhibits the growth of the second mode.

Slip effects on the receptivity of supersonic flat-plate boundary layer to freestream acoustic waves

The receptivity phase located upstream from the neutral point might be significantly affected by local rarefaction effects (especially surface slip effects) in terms of the boundary-layer transition of near-space hypersonic vehicles. In this paper, the receptivity of a supersonic flat-plate boundary layer to freestream acoustic waves in no-slip and slip flows is analyzed using direct numerical simulations and linear stability theory. The Maxwell–Smoluchowski velocity-slip and temperature-jump boundary conditions are adopted at the wall to account for surface slip effects. A Mach 4.5 flow at different wall-cooling degrees is mainly analyzed, and another Mach-6 case is presented, both with freestream unit Reynolds number on the order of 1×106/m. The main goal is to clarify the qualitative and quantitative influence of surface slip effects on the receptivity phase under different conditions. The results show that the receptivity mechanism in the slip flow is similar to that in the no-slip flow. That is, the mode S or F is excited near the leading edge due to synchronization with slow or fast acoustic waves, and the Mack second mode is excited further downstream after synchronization between modes S and F. However, the slip effects lead to distinctly quantitative differences in receptivity. The slip effects have little influence on the excitation of mode S or F near the leading edge but largely affect the evolution (intermodal exchange) of modes S and F as propagating downstream. Consequently, as for the receptivity to slow acoustic waves, the slip effects play a stabilizing role in receptivity when mode S is stable while a destabilizing role when mode S converts to the first mode in the upstream. As for receptivity to fast acoustic waves, as slip degree increases, the slip effects initially stabilize and then destabilize the receptivity, where the receptivity coefficient of the tested slip case can increase by 25% compared with the no-slip case.

Direct numerical simulation of fully-developed supersonic turbulent channel flows with dense vapors

This work aims to investigate the impact of the molecule-complexity effect and the non-ideal effect on wall-bounded turbulent flows by applying direct numerical simulation (DNS) to fully-developed channel flows of two typical organic vapors: R1233zd(E) and octamethyltrisiloxane (MDM). For each vapor, three thermodynamic states are analyzed: one in the dilute-gas region, one near the saturation line, and one in the supercritical region. For mean flow fields, it is found that, due to smaller Prandtl and Eckert numbers, both the molecule-complexity effect and the non-ideal effect reduce the mean temperature rise from the cold wall to the channel center. Meanwhile, the molecule-complexity effect weakens the mean density drop, while the non-ideal effect strengthens the drop. Furthermore, once the density and viscosity variations are considered, the mean streamwise velocity profiles of dense vapors are practically the same as the ideal gas. For turbulent fluctuations, it is found that the correlations between T′, p′, and ρ′ in dense vapors are more complicated than the ideal gas: for the ideal gas, fluctuations are dominated by “vorticity mode”; hence, ρ′ and T′ are strongly related to u′ but independent of p′; however, for dense vapors, “acoustic mode” can also play an important role. A newly derived equation illustrates that, through the “acoustic mode,” the molecule-complexity effect obviously enhances the positive correlation between ρ′ and p′, while the non-ideal effect can enhance the positive correlation between T′ and p′. Further analysis of instantaneous flow fields shows that p′ is isotropic. The isotropic character affects fluctuation magnitudes but has limited effect on the specified wall-direction turbulent transport. Consequently, Walz's equation and Reynolds analogy in terms of enthalpy are still valid. Finally, a comparison between the DNS energy budget and k equation of Reynolds-averaged Navier–Stokes (RANS) model has been carried out. Results show that obvious deviation happens on the production term in spite of the careful selection of eddy viscosity model.

Direct numerical simulation of skin friction drag reduction on supersonic turbulent boundary layers with micro-blowing

The flow control mechanism and skin friction drag reduction characteristics of micro-blowing on a Ma2.25 supersonic turbulent boundary layers are investigated through direct numerical simulations, and the effects of blowing intensity and micro-hole arrangement on turbulent structure and skin friction drag in the local control region and downstream region are considered. The results show that the skin friction drag decreases remarkably in the control region under the influence of micro-blowing, and a certain drag reduction can still be maintained in the downstream region. The drag reduction performance in the control region is jointly determined by blowing intensity and micro-hole arrangement. The drag reduction performance of the staggered arrangement is 5.7% and 11.1% higher than that of the inline arrangement at blowing intensities of 0.2% and 0.5%, respectively. However, it is found that the drag reduction in the downstream region is only determined by the blowing intensity and almost independent of the micro-hole arrangement. The effect of micro-blowing on turbulent structures is quantitatively characterized by energy spectrum analysis, and it shows that the streamwise scales of the near-wall streaks are significantly reduced under the influence of micro-blowing. In addition, the compressibility of fluids and the local reverse transfer in the strong expansion region are significantly improved under the influence of micro-blowing. These effects should be considered when performing Large Eddy Simulation modeling of supersonic turbulent boundary layers with micro-blowing.

An investigation of anisotropy in the bubbly turbulent flow via direct numerical simulations

We investigated the effects of bubble count, flow direction, and Eötvös number on deformable bubbles in turbulent channel flow. For a given shear Reynolds number Re = 180 and fixed bubble volume fractions (1.263% and 2.525%), we conducted a series of direct numerical simulations using a coupled level-set and volume-of-fluid solver to evaluate their impact on bubble volume fraction distribution, velocity fields, and turbulence characteristics. Each aspect was studied based on the microscopic equations of two-phase flow, and the accuracy of the modeling terms used in current Reynolds-averaged Navier–Stokes equation (RANS) models was assessed. The influence on the anisotropic state was analyzed using the Lumley triangle, and the anisotropy of Reynolds stresses was captured through the exact balance equations. The results indicate that in upward flow, bubbles tend to accumulate near the wall, with smaller Eötvös numbers leading to closer proximity to the wall and greater attenuation of the liquid-phase velocity. This distribution enhances energy dissipation and turbulence isotropy. In downward flow, bubbles cluster in the channel center, generating additional pseudo-turbulence and attenuating energy in the buffer layer. Moreover, the interfacial transfer of turbulent energy, as currently modeled in RANS, is found to be inadequate for upward flows.

An improved viscoelastic k- ε -v2-f turbulence model for intermediate and high drag reduction regimes

In this work, an improved anisotropic k-ε-v2-f model based on the finite extensible nonlinear elastic model with the Peterlin approximation for viscoelastic channel flows is proposed. This model is tested using direct numerical simulation (DNS) data for friction Reynolds numbers (Reτ) in the range of 120–1000, friction Wiesenberg numbers (Wiτ) in the range of 25–116, viscosity ratios (β) in the range of 0.6–0.9, and maximum polymer extensibility values (L2) in the range of 900–14 400. The flow characteristics of viscoelastic fluids with various parameters obtained from the new model agree well with existing DNS results. By adding closures for the flow, shear, and transverse components, the incomplete prediction of nonlinear terms from interactions between the fluctuating components of the conformation and velocity gradient tensors is improved. Compared with DNS results, these closures can fully obtain each component and significantly improve the accuracy of the flow direction component in the intermediate and high drag reduction regimes. Furthermore, the model in this paper retains the advantages of an anisotropic model, does not require a damping function, is simple to construct, and is easily extended to a variety of bounded flows.

Direct numerical simulation of dispersion and mixing in gas–liquid Dean-Taylor flow with influence of a 90° bend

Gas-liquid capillary flow finds widespread applications in reaction engineering, owing to its ability to facilitate precise control and efficient mixing. Incorporating compact and regular design with Coiled Flow Inverter (CFI) enhances process efficiency due to improved mixing as well as heat and mass transfer leading to a narrow residence time distribution. The impact of Dean and Taylor flow phenomena on mixing and dispersion within these systems underscores their significance, but is still not yet fully understood. Direct numerical simulation based on finite element method enables full 3D resolution of the flow field and detailed examination of laminar flow profiles, providing valuable insights into flow dynamics. Notably, the deflection of flow velocity from the center axis contributes is followed by tracking of particle with defined starting positions, aiding in flow visualization and dispersion characterization. In this CFD study, the helical flow with the influence of the centrifugal force and pitch (Dean flow) as well as the capillary two-phase flow (Taylor bubble) is described and characterized by particle dispersion and related histograms. Future prospects in this field include advancements in imaging techniques to capture intricate flow patterns, as well as refined particle tracking methods to better understand complex flow behavior.

Wake flow behind a transversely rotating sphere confined in a pipe at low Reynolds numbers

Direct numerical simulations are performed to investigate the flow past a transversely rotating sphere confined in a pipe. Three sphere diameters of d=0.2D, 0.5D, and 0.8D (D is the pipe diameter) and sphere Reynolds numbers (based on the sphere cross-sectional mean velocity and d) of Res=250–500 are considered. The simulation results highlight the combined effects of the blockage ratio (BR=d/D) and the nondimensionalized rotation rate Ω in the range of 0.2–1.2 on the wake flow. For each BR, with increasing Res and Ω, the wake flow undergoes a transition of the flow pattern from steady symmetry and unsteady symmetry to unsteady asymmetry. For BR=0.2, the variation in flow behaviors with Ω is similar to those reported in the previous studies on a transversely rotating sphere subjected to a uniform flow. With increasing BR from 0.2 to 0.5, the unsteady asymmetric wake flow is characterized by distorted near-wake vortex rings and a series of staggered inward rolled-up vortex arcs on the double streamwise vortex threads, which results in a low-frequency modulation of the unsteady wake flow. These rolled-up vortex arcs are connected with downstream quasi-streamwise vortices and form hairpin-like vortices. Dynamic Mode Decomposition (DMD) analysis shows that the distorted near-wake vortex rings are associated with the antisymmetric DMD mode at the adjacent frequency around the symmetric primary DMD mode. With BR=0.8, there is a drastic increase in the sphere's force coefficients. A more complex wake flow behavior is found due to a strong interaction between the wake flow and the near-pipe-wall flow. At Ω=1.2, dominant asymmetric flow structures are characterized by impingement of the vortex filaments. The low-frequency modulation of these wake structures is also analyzed using the DMD analysis.

An explicit multilevel power turbulent wall function based on the de-thresholding Douglas–Peucker algorithm

Wall functions are extensively applied in engineering simulations with turbulence. They facilitate a significant increase in the scale of the grids next to the wall, which in turn reduces the total number of grids needed. This optimization enhances computational efficiency, making the simulation process more effective and streamlined. However, the current commonly used wall functions, such as the Spalding wall function, are an implicit expression that needs to be solved iteratively, which affects the computational efficiency, and the multilayer segmented wall function is not smoothly articulated, which affects the fidelity. In this study, based on flat plate direct numerical simulation (DNS) data, combined with structural ensemble dynamics theory, the de-thresholding Douglas–Peucker algorithm is introduced to construct an explicit wall function expression in the form of multilevel power exponential concatenated multiplication. The comparison of the new wall function against DNS data reveals that it demonstrates superior fitting accuracy in contrast to the traditional ones, and eliminates the need for manual calibration, reduces subjective influence, and enhances reliability. The numerical simulation outcomes for the flat plate boundary layer and a series of airfoils showcase the new wall function's exceptional accuracy, which not only meets but also surpasses the demanding standards of engineering practice.

Hypersonic Stability and Breakdown Measurements on the AFRL/AFOSR BOLT II Flight Geometry

The Air Force Research Laboratory/Air Force Office of Scientific Research (AFRL/AFOSR) introduced the Boundary Layer Turbulence (BOLT II) flight experiment to further the understanding and modeling of a hypersonic turbulent boundary layer on a geometry that features concave surfaces with swept leading edges. In support of the flight experiment, this paper identifies and documents the transition instabilities and quantifies the breakdown to turbulence on a 25%-scale BOLT II model in hypersonic flow. Experiments were conducted in the conventional Actively Controlled Expansion wind tunnel facility and the Mach 6 Quiet Tunnel located at the Texas A&M University National Aerothermochemistry and Hypersonics Laboratory. Global surface heating was viewed using infrared thermography with on-surface measurements made using high-frequency Kulite and PCB Piezotronics surface pressure transducers. Modal growth was observed among the sensors both upstream and downstream of the model. Off-surface measurements were made in the mixed-mode region utilizing constant-temperature hot-film anemometry in the Mach 6 Quiet Tunnel. Comparisons were made to the Direct Numerical Simulation (DNS) results from the University of Minnesota. Amplified content was observed between f=53 and 75 kHz in the root-mean-square voltage fluctuations of the hot-film measurements.

Effects of Mach number on space-time characteristics of wall pressure fluctuations beneath turbulent boundary layers

Wall pressure fluctuations beneath turbulent boundary layers are a fundamental source of aerodynamic noise by exciting the wall structure, with their space-time characteristics serving as the basic ingredient for predicting the wall structural response. To this end, direct numerical simulations of fully developed compressible turbulent boundary layers at Mach numbers of 0.5, 1.2, and 2.0 are conducted to investigate wall pressure fluctuations comprehensively. The effects of Mach number on the single-point statistics of wall pressure fluctuations, such as the root mean square, skewness and flatness factors, probability density function, and frequency spectrum, are assessed to be very weak. Regarding the space-time characteristics, the convection velocity Uc determined by the space-time correlation of wall pressure fluctuations increases slightly with the Mach number, which only reflects the convective behavior of turbulent vortices. On the wavenumber–frequency spectrum, characteristic peaks of both the acoustic wave and convective vortices are identified. At Mach 0.5, the peaks of the fast (Uc+c) and slow (Uc−c) acoustic waves are unattached to others with c denoting acoustic speed, while only the peak of the fast acoustic wave is distinguishable from the convective peak at Mach 1.2 and 2.0. Due to the aerodynamic heating at supersonic conditions, the thermal effect on acoustic speed should be taken into account in determining the acoustic wavenumber. By introducing a convective Prandtl–Glauert parameter, a refined relation is proposed to provide a more accurate depiction of the acoustic domain in the wavenumber–frequency spectrum.

Vector rogue waves in spin-1 Bose–Einstein condensates with spin–orbit coupling

Abstract We analytically and numerically study three-component rogue waves (RWs) in spin-1 Bose–Einstein condensates with Raman-induced spin–orbit coupling (SOC). Using the multiscale perturbative method, we obtain approximate analytical solutions for RWs with positive and negative effective masses, determined by the effective dispersion of the system. The solutions include RWs with smooth and striped shapes, as well as higher-order RWs. The analytical solutions demonstrate that the RWs in the three components of the system exhibit different velocities and their maximum peaks appear at the same spatiotemporal position, which is caused by SOC and interactions. The accuracy of the approximate analytical solutions is corroborated by comparison with direct numerical simulations of the underlying system. Additionally, we systematically explore existence domains for the RWs determined by the baseband modulational instability (BMI). Numerical simulations corroborate that, under the action of BMI, plane waves with random initial perturbations excite RWs, as predicted by the approximate analytical solutions.

Large eddy simulation analysis of a model reactive tracer through spatial filtering

Large eddy simulations (LES) provide a methodology for both analyzing and simulating multi-scale flows when the smallest scales of motion cannot be resolved. Within environmental flows there exist numerous biogeochemical processes involving tracers undergoing reactions. In this study, we perform an a posteriori LES analysis on a direct numerical simulation of an idealized model reactive tracer subjected to three-dimensional turbulence induced by a Rayleigh–Taylor instability. The governing equations, including an advection–diffusion–reaction equation for the reactive tracer, are filtered, and the resulting sub-filter-scale terms are expressed in terms of interactions between scales. The procedure is demonstrated for a generalized degree N polynomial reaction function. Various spectral filters are applied to the data and compared. The preferential choice is to use the widest filter possible with a smoothed cutoff. The sub-filter-scale reaction term that results from filtering the reaction function is considered for each of the filter choices. When using a particularly harsh filter, local balances are found for the resolved scale and cross-scale components of the sub-filter-scale reaction term. The same result is shown for the vertical sub-filter-scale flux for both a reactive and a passive tracer. The components of the sub-filter-scale reaction and vertical flux terms involving interactions at the sub-filter-scale do not show any evidence of local balances and are distributed around the fine turbulent structures in the flow. This suggests that parameterizations for the sub-filter-scale terms would benefit from considering event specific dynamics.

Scaling of coherent structures in compressible wall-bounded turbulence

Semi-local scales have been widely used in compressible wall-bounded turbulence, but it is still unclear whether they are applicable to the scaling of coherent structures, especially under conditions of high Mach number and cold wall temperature. By scrutinizing the direct numerical simulation dataset at different Mach numbers and wall temperatures, this paper demonstrates that the coherent structures normalized by semi-local scales are universal in size. In addition to this, we find that the ratios of Kolmogorov scales to semi-local scales are independent of Mach number and wall temperature. Thus, Kolmogorov scales can achieve the same scaling effect as the semi-local scales. The velocity spectra are also compared to verify the current scaling method quantitatively. A method to determine the threshold for the vortex identification criterion is proposed, allowing the same threshold for different cases to obtain vortices of similar size. The scaling of other statistics including turbulent kinetic energy, streamwise Reynolds normal stress, and root mean square of fluctuating vorticity is also investigated. A new velocity scale is proposed based on the total-stress-based transformation for mean streamwise velocity, which can collapse the profiles of these statistics more accurately than the semi-local velocity scale. The present paper demonstrates that through appropriate normalization, the structures and statistics of compressible turbulence become universal, reaffirming the validity of Morkovin's hypothesis even for the present high Mach number and cold wall cases.