Velocity Fluctuations Research Articles

Effects of Thermal Expansion and Pressure on Far-Field Fluctuations of Heterogeneous Propellants

In this work, we carry out three-dimensional mesoscale simulations of heterogeneous solid propellant combustion. We solve the reactive low-Mach-number equations in the gas phase with complete coupling to the solid phase. The model takes into account thermal expansion and deformation in the solid phase by using a hypoelastic law in the quasi-static limit. To account for morphology, we select two different propellant formulations with different particle size distributions and present the results as a function of pressure. The thermomechanical behavior is accessed by examining quantities such as strain and stress in the propellant as a function of pressure and propellant morphology. We also show that temperature and velocity fluctuations exist in the far field above the propellant surface and that these fluctuations can be significant. To better understand the nature of these fluctuations, we vary the pressure and make relevant plots of normal velocity and temperature probability density functions, as well as time autocorrelations. Such descriptions are necessary to account for the coupling between the mesoscale and the macroscale, where the fluctuations at the mesoscale can affect quantities at the macroscale, such as head-end pressure, trigger parietal vortex shedding, and aeroacoustics.

Fluctuation-response relation as a probe of long-range correlations in nonequilibrium quantum and classical fluids.

The absence of a simple fluctuation-dissipation theorem is a major obstacle for studying systems that are not in thermodynamic equilibrium. We show that for a fluid in a nonequilibrium steady state characterized by a constant temperature gradient the commutator correlation functions are still related to response functions; however, the relation is to the bilinear response of products of two observables, rather than to a single linear response function as is the case in equilibrium. This modified fluctuation-response relation holds for both quantum and classical systems. It is both motivated and informed by the long-range correlations that exist in such a steady state and allows for probing them via response experiments. This is of particular interest in quantum fluids, where the direct observation of fluctuations by light scattering would be difficult. In classical fluids it is known that the coupling of the temperature gradient to the diffusive shear velocity leads to correlations of various observables, in particular temperature fluctuations, that do not decay as a function of distance, but rather extend over the entire system. We investigate the nature of these correlations in a fermionic quantum fluid and show that the crucial coupling between the temperature gradient and velocity fluctuations is the same as in the classical case. Accordingly, the nature of the long-ranged correlations in the hydrodynamic regime also is the same. However, as one enters the collisionless regime in the low-temperature limit the nature of the velocity fluctuations changes: they become ballistic rather than diffusive. As a result, correlations of the temperature and other observables are still singular in the long-wavelength limit, but the singularity is weaker than in the hydrodynamic regime.

An LES investigation on turbulent flow in a centrifugal impeller affected by normal-distributed inflow perturbations

The inflow of a rotating centrifugal impeller is normally perturbed by an upstream stationary component; therefore, the development of turbulent flow is different from the case with steady and uniform inflow. In this work, we performed a large-eddy simulation on turbulent flow in a centrifugal impeller, considering perturbation from the inflow and emphasizing the development of perturbation and its influence on flow in the impeller. The inflow perturbation is applied for the streamwise (w-) velocity and is time-varying as generated by a random number generator. A normal-distributed pattern of perturbation is always assumed with the intensity of perturbation, defined as the ratio between the perturbation amplitude and the mean velocity, set as fv = 0%, 5%, 10%, and 20%, where fv denotes the perturbation intensity. The inflow perturbation notably affects the passage flow. The velocity fluctuation and secondary flow increase in intensity as the perturbation intensity increases from fv = 0% to 10%, while a further increase to fv = 20% slightly weakens the velocity fluctuation. Although this phenomenon is less obvious in terms of the time-averaged characteristics of velocity, the Reynolds stress terms CtCa and CrCa under time-averaging still reflect a clear variation trend, and the Reynolds stresses are observed significantly on the blade suction surface.

Parametric study of large-eddy simulation to capture scaling laws of velocity fluctuations in neutral atmospheric boundary layers

The development of digital twins for wind farms often involves the use of large-eddy simulation (LES) to model atmospheric boundary layers. Existing LES solvers primarily focus on accurately capturing streamwise fluctuations. They, however, overlook the less energetic cross-stream fluctuations, which play a crucial role in wind turbine wake evolution. In this study, we conduct a systematic parametric study and incorporate changes in an open-source LES solver. The improved solver is able to predict all three components of velocity fluctuations in alignment with the scaling laws derived from the attached-eddy hypothesis. In particular, we examine the impact of (i) the subgrid-scale model, (ii) the wall model, (iii) the von Kármán constant, and (iv) the grid-cell aspect ratio. We find that although all these factors influence the prediction of velocity fluctuations, the grid-cell aspect ratio has the greatest effect on the spanwise and vertical velocity components. Notably, utilizing nearly isotropic grid cells leads to the best alignment of all three velocity component fluctuations with the scaling laws. Spectral analysis further demonstrates that the present LES solver accurately predicts the characteristic length scales for all velocity fluctuation components, making it a reliable tool for obtaining turbulent inflow conditions for wind farm modeling.

An experimental study of the flow structure in arrays of cold bidirectional swirling jets

Jets arrays have become a promising technology for industrial applications, including heat transfer and combustion processes. This paper presents the first experimental results of arrays of bidirectional swirling jets, including their flow structures and specific turbulent interactions. Particle image velocimetry of non-reacting linear arrays showed the formation of recirculation zones within each jet of the array. These regions are the most stable, with fluctuations less than 0.002 of the bulk inlet velocity Vin. In addition, jets merging locations also show a significant decrease in velocity fluctuations, with values V′/Vin < 0.005, which is most clearly seen at a smaller pitch distance. However, the highest turbulent fluctuations occur in shear layers, reaching values of V'/Vin ≈ 0.02 for linear arrays and V'/Vin ≈ 0.03 for planar configurations. The flow structure of planar arrays is more complex, leading to the formation of secondary vortex structures at locations of jet interaction. This results in a significant decrease in the flow swirl number, from an initial value of Sin = 2.1 in the vortex chamber, to S1 = 0.36–0.49 at z/dout2 = 1 and S2 = 0.19–0.27 at z/dout2 = 2 beyond the outlet nozzle. Strong cross-sectional motions are observed in planar arrangements, leading to the development of additional regions of negative axial velocity between the outlet nozzles. The boundaries of these regions have nearly zero axial velocity values and can provide reliable ignition and flame stabilization in the case of reactive bidirectional jets. Therefore, the flow structure in planar arrays forms a large stabilization zone, and the obtained results can be used to develop a new multipoint combustion technique.

Tracking the Nearfield Evolution of an Initially Shallow, Neutrally Buoyant Plane Jet Over a Sloping Bottom Boundary

AbstractUnderstanding coastal plane jets which occur when a body of water discharges into an ocean or a lake through a channel or outlet is important, since they play a significant role in sediment, nutrient, and pollutant exchange. This study investigates the nearfield of initially shallow, neutrally buoyant plane jets, bounded by a free surface and a sloping bottom (Sloping Bottom Jet; SBJ) that issue into a laterally unconfined quiescent ambient both experimentally and numerically, and compares them with a plane jet flowing over a horizontal bottom (Horizontal Bottom Jet; HBJ). Results revealed that, different from the HBJ, the width and centerline velocity of SBJ decrease near the mouth. The SBJ width gradually increases after that as the transverse longitudinal velocity profile progressively transforms from a “top‐hat” into a Gaussian distribution. Once the Gaussian distribution is established, both jets diverge and centerline velocity decreases. Shear layers are generated on the sides of both jets with Kelvin Helmholtz‐type Coherent Structures (KHCS) developing inside. KHCS produce periodic velocity fluctuations with a Strouhal number of ∼0.079 and contribute significantly to momentum exchange and turbulent kinetic energy production. Since the thickness of the SBJ increases longitudinally, the vertical extent of KHCS also increases. When the two shear layers meet and merge at the centerline, they cause a flapping motion of the jet. This location is closer to the jet mouth for SBJs than for the HBJ. These findings demonstrate that a sloping bottom modifies the flow field from quasi‐2D for the HBJ to strongly 3D for SBJs.

MHD Inertial and Energy-containing Range Turbulence Anisotropy in the Young Solar Wind

We study solar wind turbulence anisotropy in the inertial and energy-containing ranges in the inbound and outbound directions during encounters 1–9 by the Parker Solar Probe (PSP) for distances between ∼21 and 65 R ⊙. Using the Adhikari et al. approach, we derive theoretical equations to calculate the ratio between the 2D and slab fluctuating magnetic energy, fluctuating kinetic energy, and the outward/inward Elsässer energy in the inertial range. For this, in the energy-containing range, we assume a wavenumber k −1 power law. In the inertial range, for the magnetic field fluctuations and the outward/inward Elsässer energy, we consider that (i) both 2D and slab fluctuations follow a power law of k −5/3, and (ii) the 2D and slab fluctuations follow the power laws with k −5/3 and k −3/2, respectively. For the velocity fluctuations, we assume that both the 2D and slab components follow a k −3/2 power law. We compare the theoretical results of the variance anisotropy in the inertial range with the derived observational values measured by PSP, and find that the energy density of 2D fluctuations is larger than that of the slab fluctuations. The theoretical variance anisotropy in the inertial range relating to the k −5/3 and k −3/2 power laws between 2D and slab turbulence exhibits a smaller value in comparison to assuming the same power law k −5/3 between 2D and slab turbulence. Finally, the observed turbulence energy measured by PSP in the energy-containing range is found to be similar to the theoretical result of a nearly incompressible/slab turbulence description.

Investigating the use of 3-component-2-dimensional particle image velocimetry fields as inflow boundary condition for the direct numerical simulation of turbulent channel flow

Direct numerical simulation (DNS) of turbulent wall-bounded flows requires long streamwise computational domains to establish the correct spatial evolution of large-scale structures with high fidelity. In contrast, experimental measurements can relatively easily capture large-scale structures but struggle to resolve the dissipative flow scales with high fidelity. One methodology to overcome the shortcomings of each approach is by incorporating experimental velocity field measurements into DNS as an inflow boundary condition. This hybrid approach combines the strengths of DNS and experimental measurements, allowing for a reduction in the streamwise computational domain and accelerated development of large-scale structures in turbulent wall-bounded flows. To this end, this paper reports the results of an investigation to establish the impact of limited spatial resolution and limited near-wall experimental inflow data on the DNS of a wall-bounded turbulent shear flow. Specifically, this study investigates the spatial extent required for the DNS of a turbulent channel flow to recover the turbulent velocity fluctuations and energy when experimental inflow data is typically unable to capture fluctuations down to the viscous sub-layer or the smallest viscous scales (i.e. the Kolmogorov scale or their surrogate viscous scale in wall-bounded turbulent shear slows) is used as the inflow to a DNS. A time-resolved numerically generated experimental field is constructed from a periodic channel flow DNS (PCH-DNS) at Reτ=\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re_{\ au } =$$\\end{document} 550 and 2300, which is subsequently used as the inflow velocity field for an inflow–outflow boundary conditions DNS. The time-resolved experimental inflow field is generated by appropriately filtering the small scales from the PCH-DNS velocity by integrating over a spatial domain that is representative of a particle image velocimetry interrogation window. This study shows that the recovery of small scales requires a longer domain as the spatial resolution at the inflow decreases with all flow scales recovered and their correct scale-dependent energy is re-established once the flow has developed for 3 channel heights.

Direct numerical simulations of suspension of disk-shaped particles

This study investigates the dynamics of disk-shaped particles using direct numerical simulations with the smoothed profile method for rigid particles. These disk-shaped particles are formed by joining the spherical beads and are allowed to settle/sediment in a Newtonian fluid. The concentration effects of the mono-dispersed particles are studied in the Stokes regime, varying the volume fraction (ϕ) from 0.0003 to 0.1. Strong inhomogeneities in the system were noticed, producing multiple peaks in the radial distribution function caused by the orientation preference of particles, while settling. A histogram analysis of the particles' orientation angle suggests that particles prefer horizontal orientation at very low volume fractions and then start orienting vertically with subsequent increase in the volume fraction. Average settling velocity increases initially till volume fraction 0.001, creating a local maxima, and then decreases monotonically following the Richardson–Zaki law. It was also found that velocity fluctuations increased with increasing volume fraction, following the ϕ1/3 trend. These fluctuations are smaller than those of rod-like particles and larger than spherical particles, though the qualitative trend is quite similar.

Implementation of a PAFV turbulence model for airfoil flow simulation on OpenFOAM

To further develop a more effective turbulence model and to improve the calculation accuracy of the flow around airfoil, a new PAFV turbulence model has here been constructed by using a deformation rate tensor and the grouping of an average fluctuation velocity. To evaluate the applicability of the PAFV turbulence model, the numerical calculations of flow around the airfoil have here been implemented, which was based on the OpenFOAM calculation platform. On the basis of grid independence research, the model was used to calculate the low-speed flow-around problem for the plano-convex airfoil NACA4412 and the symmetric airfoil NACA0012. It was also compared with the S-A (Spalart-Allmaras) and SST (Shear Stress Transport) k-ω turbulence models. Firstly, the maximum lift angle-of-attack case of the NACA4412 airfoil was calculated. Thereafter, numerical calculations were performed for the flow around the airfoil in the multi-angle-of-attack case of the NACA0012 airfoil. The results showed that the NACA4412 airfoil had an obviously separated vortex at the trailing edge of the airfoil at the maximum lift angle of attack. Also, there was a certain velocity loss downstream of the trailing edge, which was calculated by all three models. However, the results of the PAFV turbulence model were found to be better than those of the S-A and SST turbulence models. The three turbulence models showed comparable accuracies for the calculations of the surface pressure coefficients of the NACA0012 airfoil. However, the S-A and SST k-ω turbulence models were slightly better for the calculations of the mean velocity profiles of the NACA0012 airfoil. Also, the PAFV turbulence model was more accurate for the calculations of the lift and drag coefficients. In conclusion, the PAFV model can make effective predictions of the airfoil low-speed flow around the problem at hand, which in turn preliminarily verifies the applicability of this turbulence model for the low-speed flow around the airfoil problems.

Numerical evaluation of the bubble dynamic influence on the characteristics of multiscale cavitating flow in the bluff body wake

Unsteady cavitating flow around the bluff body widely existed in practical engineering fields. This paper aims to investigate the physics involved in the multi-scale cavitating flow in the wake of a wedge-shaped bluff body, with special emphasis on the bubble effects. Numerical simulation is conducted using the traditional Eulerian-Eulerian (E-E) and newly developed Eulerian-Lagrangian (E-L) methods in OpenFOAM. It is found that, compared with the experimental results, the E-L method can predict more accurate cavitation characteristics than the E-E method due to the contribution of discrete bubbles. Bubbles can significantly influence the behaviors of the multi-scale cavitating flow. Specially, analysis of pressure fluctuations indicates that the bubbles induce a stronger power spectral density in the middle- and high-frequency regions. With the bubble effects considered, the vortex structure is stretched and more extensive in the near wake regions. The strength of time-averaged vorticity distribution decreases in the far wake region. Further analysis of the vortex transport equation shows that bubble dynamics significantly increase the intensity of the stretching term by enhancing the spanwise velocity fluctuations. The intensity of the baroclinic torque term predicted by the E-L method is higher than that by the E-E method near the center line due to the influence of bubbles. On the other hand, the analysis of turbulent kinetic energy distribution indicates that bubbles can induce turbulence by increasing cross-stream velocity fluctuations. The increment of the shear Reynold stress, u′xu′y‾, suggests that the bubbles intensify the coupling between the streamwise and cross-stream velocity fluctuations.

Effect of slug characteristics on the nonlinear dynamic response of a long flexible fluid-conveying cylinder

Hydrocarbon flows in a marine riser may appear in multiple phases with varying flow patterns, among which, slug flows are known to exhibit complex flow characteristics due to fluctuations in multiphase mass, velocities, and pressure changes. Bulk of the literature concerning internal single and two-phase flow induced vibrations have largely focused on the use of a linearized tension beam model. In this study, the fluid-structure interaction phenomena of a submerged long flexible cylinder conveying two-phase slug flows is investigated taking into account the geometric and hydrodynamic nonlinearities. A semi-empirical theoretical model which consists of nonlinear structural equations of coupled out-of-plane, in-plane and axial structural motions is presented for the analysis and prediction of two-phase slug flow induced vibrations (SIV). The model includes equations involving centrifugal and Coriolis forces to capture mass variations in the slug flow regime. A numerical space-time finite difference approach combined with a frequency domain analysis is used for analyzing the highly nonlinear three-dimensional responses. The model assumes constant geometric properties across the span of the cylinder and utilizes an idealized slug unit concept, wherein slug properties are considered to be fully developed and undisturbed by pipe oscillations. Model validations are performed through comparisons with published experimental internal flow-induced vibration results. Parametric study captures the effects of several key slug characteristics on the vibration responses of a fluid-conveying cylinder. The results demonstrate amplitude-modulation response, mean displacements, three-dimensional cylinder displacements due to geometric nonlinear couplings and the influence of internal flow velocities. Higher dominant modes and oscillation frequencies were observed when the cylinder experiences large amplitude motions at specific slug formations. A new dimensionless parameter is introduced to predict scenarios of high amplitude oscillations for long flexible cylinders carrying two-phase slug flows.

Modeling Internal Flow Patterns of Sessile Droplets on Horizontally Vibrating Substrates

A three-dimensional Navier–Stokes and continuity equation model is employed to numerically predict the resonant modes of sessile droplets on horizontally vibrating substrates. A dynamic contact angle model is implemented to simulate the contact angle variations during vibrations. The four resonant modes (n = 1, 2, 3 and 4) of a droplet under horizontal vibrations are investigated. Simulations are compared to experimental results for validation. Excellent agreement is observed between predicted results and experiments. The model is used to simulate the internal flow patterns within the droplet under resonant modes. It is found that the flow in all four resonant modes can be divided into the Stokes region, the gas–liquid interface region, and the transition region located in between. Numerical simulations show that the average velocity within the droplet increases with the increase in frequency, while the fluctuations in average velocity after reaching the steady state show different trends with the increase in frequency. It is also found that with an increase in the order of resonant modes, the contact angle difference between the two sides of the droplet increases, and the contact angle difference of the droplet is maximized when the applied frequency is the resonant frequency of the specified mode.

Research on the Mechanism of Flow-Induced Vibration in the Cooling System of a Double Crystal Monochromator

To investigate the mechanism of flow-induced vibrations in the cooling system of a double crystal monochromator (DCM), this paper utilizes a multi-physics numerical simulation approach, employing ANSYS and FLUENT platforms to simulate the flow state of liquid nitrogen in the cooling system and explore the amplitude response of the DCM. Initially, simulations were conducted to examine the flow state of liquid nitrogen with varying frequency and amplitude pulsations. Subsequently, modal analysis was employed to investigate the amplitude response of the DCM in the pitch direction vibrations under pulsating excitation. Finally, this research investigated the influence of high heat load-induced liquid nitrogen boiling on a DCM. The results indicate that pipe resistance is the fundamental cause of vibration induced by pulsating excitation. Low-frequency excitation enhances the amplification factor of DCM vibration. In contrast, due to the rapid conversion of fluid kinetic energy to pressure potential energy, high-frequency excitation increases the pulsation amplitude in the pipe. Additionally, there is a linear relationship between the amplitude of liquid nitrogen velocity fluctuations and the response amplitude of a DCM. The slug flow formed after liquid nitrogen boiling generates low-frequency pulse signals, and intermittent fluid impacts cause significant vibrations in the DCM. These research findings provide a reference for the analysis and design of ultra-high-stability DCM cooling systems.

From clusters of moving molecules to continua: Material elements as open systems

We discuss the discrete-to-continuum transition in the description of matter, starting from a cluster of flowing molecules with equal mass and ending up with a non-simple fluid. We account for the fluctuations beyond a local affine approximation of the velocity distribution of molecules within a space window adopted to compute some prominent statistics. The resulting continuum picture accounts for local mass variation in each space window corresponding to a point in the continuum scale. From a statistical viewpoint, every material element is thus considered as a grand-canonical ensemble. So-called C-derivatives account for macroscopic-to-mesoscopic relative motion. When considered for second-rank tensors, they extend Truesdell’s derivative and reduce to Oldroyd’s one when the macroscopic-to-mesoscopic relative motion vanishes. Fluctuations beyond an affine component best approximating the kinetic energy are summarized into a second-rank symmetric tensor whose time variation enters a balance equation governing transfer from velocity fluctuations to heat. Eventually, we discuss essential elements of thermodynamics in the present setting. What emerges is the possibility of a non-Fourier type heat transfer. The results address computational schemes for field representations of sparse phase dynamics, such as granular materials, and the one of bodies with transport of scattered molecules, such as pollutants in fluids or proteins in biological tissues.

Effect of Angular Speed Variations on the Nonlinear Vibrations of a Rotational Spring-Mass System

A rotating spring-mass system is considered using polar coordinates. The system contains a cubic nonlinear spring with damping. The angular velocity harmonically fluctuates about a mean velocity. The dimensionless equations of motion are derived first. The velocity fluctuations resulted in a direct and parametric forcing terms. Approximate analytical solutions are sought using the Method of Multiple Scales, a perturbation technique. The primary resonance and the principal parametric resonance cases are investigated. The amplitude and frequency modulation equations are derived for each case. By considering the steady state solutions, the frequency response relations are derived. The bifurcation points are discussed for the problems. It is found that speed fluctuations may have substantial effects on the dynamics of the problem and the fluctuation frequency and amplitude can be adjusted as passive control parameters to maintain the desired responses.

PIV Experimental Research and Numerical Simulation of the Pigging Process

The initial running speed of the pig during gas–liquid two-phase pipeline pigging can significantly influence the velocities of both gas and liquid phases within the pipeline. However, due to the complexity and limited understanding of these velocity variations, developing an effective operational plan for pigging becomes challenging. To enhance pigging efficiency and effectively seal the pig, it is crucial to monitor the velocity variations in the gas–liquid phase within the pipeline. In this study, an experimental platform was established to facilitate precise observation of these variations. Particle image velocimetry (PIV) technology was employed for a comprehensive understanding of gas–liquid two-phase velocities during pig operation in the pipeline. The experimental results demonstrated that increasing both the velocity and the initial liquid level height of the pig resulted in a corresponding augmentation of velocity fluctuation range. Specifically, at a holdup rate of 30%, there was a 10% reduction in the maximum liquid-phase velocity, while at a holdup rate of 25%, this reduction amounted to 16% compared to the pigging velocity.

Operational mechanism of valved-pulsejet engines

The study investigates twelve configurations of the self-aspirated valved pulsejet, focusing on its operational mechanism. It begins by characterizing the engine's acoustic field with varied boundary conditions, revealing an extended effective acoustic length beyond the combustion chamber and tailpipe sections when the valve is open. Radial and lateral velocity fluctuations are hypothesized as causes for vortex structure generation. Analysis of reacting pulsejets using gasoline and ethanol fuels shows that the operating frequencies are unaffected by fuel type, with geometric arrangement as the primary frequency factor. Fluctuations in dynamic pressure, microphone readings, and thrust characterize engine stability, impacted by low-frequency modes and higher harmonics. Tailpipe length emerges as the key geometric factor for performance enhancement with the medium size being the optimal option. Pressure field analysis demonstrates shock wave generation during ignition, serving as an excitation mechanism. Notably, different configurations with varied operating frequencies can yield equivalent thrust, indicating that thrust production is not solely dependent on pressure rise during combustion, despite similar pressure rise observed across most cases.

Large eddy simulations of a turbulent flow with hybrid nanofluid subjected to symmetric and asymmetric heating

The present study includes a comprehensive numerical analysis of the usage of nanofluids inside a thermal receiver to assess the possibility of enhancing the efficiency of this device. All simulations are performed on a small portion of a thermal receiver, similar to a bi-periodic channel subjected to symmetric and asymmetric heating. To generate conventional and hybrid nanofluids, nanoparticles of Al2O3 and Al2O3-Cu are dispersed in the base fluid H2O. The values of nanoparticles' diameter is 20nm and volume fraction of the nanoparticles is 2%. The Reynolds numbers used in this study are 12 000 and 18 000, which respectively correspond to mean friction Reynolds numbers of 140 and 260. The thermal Large Eddy Simulation (T-LES) approach is used to simulate the fluid flow and heat transfer inside the computational domain. The results are compared with those of prior numerical and experimental studies to validate the method and the numerical setup. A high level of agreement is achieved. In the case of symmetric and asymmetric heating, the dimensionless findings reveal that adding nanoparticles to the base fluid has a negligible impact on the velocity fluctuations. Nanoparticles make the base fluid denser and more viscous but also more conductive, allowing it to minimize kinetic energy while boosting heat transfer rate. Results also show that using hybrid nanofluids contributes better to improving the heat transfer rate compared to the base fluid and the conventional nanofluid. Introducing nanoparticles into the base fluid decreases its temperatures slightly, especially when Al2O3-Cu nanoparticles are used instead of Al2O3 nanoparticles. In the asymmetric case, water's turbulent heat flux values are greater than nanofluids and hybrid nanofluid, implying that more conductive and viscous liquids have the weakest turbulent heat flux. In the circumstances of symmetric and asymmetric heating, the variation of the velocity fluctuations is almost the same, but the behavior of the temperature–velocity correlation is entirely different. Finally, results show that symmetric heating leads to a higher heat transfer rate than asymmetric heating. The suggested research's findings may be utilized to promote the usage of nanofluids in Concentrated Solar Power (CSP) systems.

On the streamwise velocity, temperature and passive scalar fields in compressible turbulent channel flows: a viewpoint from multiphysics couplings

It is generally believed that the velocity and passive scalar fields share many similarities and differences in wall-bounded turbulence. In the present study, we conduct a series of direct numerical simulations of compressible channel flows with passive scalars and employ the two-dimensional spectral linear stochastic estimation and the correlation function as diagnostic tools to shed light on these aspects. Particular attention is paid to the relevant multiphysics couplings in the spectral domain, i.e. the velocity–temperature ( $u-T$ ), scalar–temperature ( $g-T$ ) and velocity–scalar ( $u-g$ ) couplings. These couplings are found to be utterly different at a given wall-normal position in the logarithmic and outer regions. Specifically, in the logarithmic region, the $u-T$ and $u-g$ couplings are tight at the scales that correspond to the attached eddies and the very large-scale motions (VLSMs), whereas the $g-T$ coupling is robust in the whole spectral domain. In the outer region, $u-T$ and $u-g$ couplings are only active at the scales corresponding to the VLSMs, whereas the $g-T$ coupling is diminished but still strong at all scales. Further analysis indicates that although the temperature field in the vast majority of zones in a channel can be roughly treated as a passive scalar, its physical properties gradually deviate from those of a pure passive scalar as the wall-normal height increases due to the enhancement of the acoustic mode. Furthermore, the deep involvement of the pressure field in the self-sustaining process of energy-containing motions also drives the streamwise velocity fluctuation away from a passive scalar. The current work is an extension of our previous study (Cheng & Fu, J. Fluid Mech., vol. 964, 2023, A15), and further uncovers the details of the multiphysics couplings in compressible wall turbulence.