Investigation of species-mass diffusion in binary-species boundary layers at high pressure using direct numerical simulations

2021 ◽  
Vol 928 ◽  
Author(s):  
Takahiko Toki ◽  
Josette Bellan

Direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed with special attention to species-mass diffusion. The working fluids are nitrogen or a mixture of nitrogen and methane. Mean profiles and turbulent fluctuations of mass fraction show that their qualitative characteristics are different from those of streamwise velocity and temperature, due to the different boundary conditions. In a wall-parallel plane near the wall, the streamwise velocity and temperature have streaky patterns and the fields are similar. However, the mass fraction field at the same location is different from the streamwise velocity and temperature fields indicating that species-mass diffusion is not similar to the momentum and thermal diffusion. In contrast, at the centre and near the edge of the boundary layer, the mass fraction and temperature fields have almost the same pattern, indicating that the similarity between thermal and species-mass diffusion holds away from the wall. The lack of similarity near the wall is traced to the Soret effect that induces a temperature-gradient-dependent species-mass flux. As a result, a new phenomenon has been identified for a non-isothermal binary-species system – uphill diffusion, which in its classical isothermal definition can only occur for three or more species. A quadrant analysis for the turbulent mass flux reveals that near the wall the Soret effect enhances the negative contributions of the quadrants. Due to the enhancement of the negative contributions, small species-concentration fluid tends to be trapped near the wall.

1997 ◽  
Vol 336 ◽  
pp. 379-409 ◽  
Author(s):  
PEDRO L. GARCÍA-YBARRA ◽  
JOSE L. CASTILLO

The concentration distribution of massive dilute species (e.g. aerosols, heavy vapours, etc.) carried in a gas stream in non-isothermal boundary layers is studied in the large-Schmidt-number limit, Sc[Gt ]1, including the cross-mass-transport by thermal diffusion (Ludwig–Soret effect). In self-similar laminar boundary layers, the mass fraction distribution of the dilute species is governed by a second-order ordinary differential equation whose solution becomes a singular perturbation problem when Sc[Gt ]1. Depending on the sign of the temperature gradient, the solutions exhibit different qualitative behaviour. First, when the thermal diffusion transport is directed toward the wall, the boundary layer can be divided into two separated regions: an outer region characterized by the cooperation of advection and thermal diffusion and an inner region in the vicinity of the wall, where Brownian diffusion accommodates the mass fraction to the value required by the boundary condition at the wall. Secondly, when the thermal diffusion transport is directed away from the wall, thus competing with the advective transport, both effects balance each other at some intermediate value of the similarity variable and a thin intermediate diffusive layer separating two outer regions should be considered around this location. The character of the outer solutions changes sharply across this thin layer, which corresponds to a second-order regular turning point of the differential mass transport equation. In the outer zone from the inner layer down to the wall, exponentially small terms must be considered to account for the diffusive leakage of the massive species. In the inner zone, the equation is solved in terms of the Whittaker function and the whole mass fraction distribution is determined by matching with the outer solutions. The distinguished limit of Brownian diffusion with a weak thermal diffusion is also analysed and shown to match the two cases mentioned above.


2008 ◽  
Vol 600 ◽  
pp. 403-426 ◽  
Author(s):  
P. BURATTINI ◽  
S. LEONARDI ◽  
P. ORLANDI ◽  
R. A. ANTONIA

The turbulent flow in a two-dimensional channel with roughness on one wall is investigated using experiments and direct numerical simulations (DNS). The elements have a square cross-section with height k=0.1H (H is the channel half-width) and a streamwise spacing of 4k. The Reynolds number Reτr, based on the friction velocity at the rough wall and H, is in the range 300–1100. Particular attention is given to the rough-wall side. Measured turbulence intensities, length scales, leading terms in the turbulent kinetic energy budget, and velocity spectra are compared with those obtained from the DNS. Close agreement is found, yielding support for the simplifying assumptions in the experiment (notably local isotropy and Taylor's hypothesis) and the adequacy of the spatial resolution in the simulation. Overall, the profiles of the Reynolds normal stresses on the roughness side are almost independent of Reτr, when normalized by outer variables. Energy spectra at different locations above the rough wall collapse well at high wavenumbers, when normalized by Kolmogorov scales. In contrast to previous studies, a region of negative energy production near the location of the maximum streamwise velocity is not observed. Comparison with a smooth-wall channel, at similar values of the friction-velocity Reynolds number, highlights differences only in the streamwise velocity component near the wall.


2014 ◽  
Vol 753 ◽  
pp. 360-401 ◽  
Author(s):  
R. Vicquelin ◽  
Y. F. Zhang ◽  
O. Gicquel ◽  
J. Taine

AbstractThe role of radiative energy transfer in turbulent boundary layers is carefully analysed, focusing on the effect on temperature fluctuations and turbulent heat flux. The study is based on direct numerical simulations (DNS) of channel flows with hot and cold walls coupled to a Monte-Carlo method to compute the field of radiative power. In the conditions studied, the structure of the boundary layers is strongly modified by radiation. Temperature fluctuations and turbulent heat flux are reduced, and new radiative terms appear in their respective balance equations. It is shown that they counteract turbulence production terms. These effects are analysed under different conditions of Reynolds number and wall temperature. It is shown that collapsing of wall-scaled profiles is not efficient when radiation is considered. This drawback is corrected by the introduction of a radiation-based scaling. Finally, the significant impact of radiation on turbulent heat transfer is studied in terms of the turbulent Prandtl number. A model for this quantity, based on the new proposed scaling, is developed and validated.


2014 ◽  
Vol 760 ◽  
pp. 63-94 ◽  
Author(s):  
M. J. Philipp Hack ◽  
Tamer A. Zaki

AbstractThe influence of harmonic spanwise wall motion on bypass transition in boundary layers is investigated using direct numerical simulations. It is shown that the appropriate choice of the forcing parameters can achieve a substantial stabilization of the laminar flow regime. However, an increase of the forcing amplitude or period beyond their optimal values diminishes the stabilizing effect, and leads to breakdown upstream of the unforced case. For the optimal wall-oscillation parameters, the reduction in propulsion power substantially outweighs the power requirement of the forcing. The mechanism of transition delay is examined in detail. Analysis of the pre-transitional streaks shows that the wall oscillation substantially reduces their average amplitude, and eliminates the most energetic streaks. As a result, the secondary instabilities that precede breakdown to turbulence are substantially weakened – an effect demonstrated by linear stability analyses of flow fields from direct numerical simulations. The outcome is transition delay owing to a significant reduction in the frequency of occurrence of turbulent spots and a downstream shift in their average inception location. Finally, it is shown that the efficiency of the forcing can be further improved by replacing the sinusoidal time dependence of the wall oscillation with a square wave.


2012 ◽  
Vol 693 ◽  
pp. 28-56 ◽  
Author(s):  
Suman Muppidi ◽  
Krishnan Mahesh

AbstractDirect numerical simulations are used to study the laminar to turbulent transition of a Mach 2.9 supersonic flat plate boundary layer flow due to distributed surface roughness. Roughness causes the near-wall fluid to slow down and generates a strong shear layer over the roughness elements. Examination of the mean wall pressure indicates that the roughness surface exerts an upward impulse on the fluid, generating counter-rotating pairs of streamwise vortices underneath the shear layer. These vortices transport near-wall low-momentum fluid away from the wall. Along the roughness region, the vortices grow stronger, longer and closer to each other, and result in periodic shedding. The vortices rise towards the shear layer as they advect downstream, and the resulting interaction causes the shear layer to break up, followed quickly by a transition to turbulence. The mean flow in the turbulent region shows a good agreement with available data for fully developed turbulent boundary layers. Simulations under varying conditions show that, where the shear is not as strong and the streamwise vortices are not as coherent, the flow remains laminar.


Sign in / Sign up

Export Citation Format

Share Document