scholarly journals Triple-deck and direct numerical simulation analyses of high-speed subsonic flows past a roughness element

2015 ◽  
Vol 774 ◽  
pp. 311-323 ◽  
Author(s):  
G. Mengaldo ◽  
M. Kravtsova ◽  
A. I. Ruban ◽  
S. J. Sherwin

This paper is concerned with the boundary-layer separation in subsonic and transonic flows caused by a two-dimensional isolated wall roughness. The process of the separation is analysed by means of two approaches: the direct numerical simulation (DNS) of the flow using the Navier–Stokes equations, and the numerical solution of the triple-deck equations. Since the triple-deck theory relies on the assumption that the Reynolds number ($\mathit{Re}$) is large, we performed the Navier–Stokes calculations at $\mathit{Re}=4\times 10^{5}$ based on the distance of the roughness element from the leading edge of the flat plate. This $\mathit{Re}$ is also relevant for aeronautical applications. Two sets of calculation were conducted with the free-stream Mach number $\mathit{Ma}_{\infty }=0.5$ and $\mathit{Ma}_{\infty }=0.87$. We used different roughness element heights, some of which were large enough to cause a well-developed separation region behind the roughness. We found that the two approaches generally compare well with one another in terms of wall shear stress, longitudinal pressure gradient and detachment/reattachment points of the separation bubbles (when present). The main differences were found in proximity to the centre of the roughness element, where the wall shear stress and longitudinal pressure gradient predicted by the triple-deck theory are noticeably different from those predicted by DNS. In addition, DNS predicts slightly longer separation regions.

2018 ◽  
Vol 2018 ◽  
pp. 1-16 ◽  
Author(s):  
Mongkol Kaewbumrung ◽  
Somsak Orankitjaroen ◽  
Pichit Boonkrong ◽  
Buraskorn Nuntadilok ◽  
Benchawan Wiwatanapataphee

A mathematical model of dispersed bioparticle-blood flow through the stenosed coronary artery under the pulsatile boundary conditions is proposed. Blood is assumed to be an incompressible non-Newtonian fluid and its flow is considered as turbulence described by the Reynolds-averaged Navier-Stokes equations. Bioparticles are assumed to be spherical shape with the same density as blood, and their translation and rotational motions are governed by Newtonian equations. Impact of particle movement on the blood velocity, the pressure distribution, and the wall shear stress distribution in three different severity degrees of stenosis including 25%, 50%, and 75% are investigated through the numerical simulation using ANSYS 18.2. Increasing degree of stenosis severity results in higher values of the pressure drop and wall shear stresses. The higher level of bioparticle motion directly varies with the pressure drop and wall shear stress. The area of coronary artery with higher density of bioparticles also presents the higher wall shear stress.


2011 ◽  
Vol 8 (64) ◽  
pp. 1594-1603 ◽  
Author(s):  
A. Kazakidi ◽  
A. M. Plata ◽  
S. J. Sherwin ◽  
P. D. Weinberg

Atherosclerotic lesions have a patchy distribution within arteries that suggests a controlling influence of haemodynamic stresses on their development. The distribution near aortic branches varies with age and species, perhaps reflecting differences in these stresses. Our previous work, which assumed steady flow, revealed a dependence of wall shear stress (WSS) patterns on Reynolds number and side-branch flow rate. Here, we examine effects of pulsatile flow. Flow and WSS patterns were computed by applying high-order unstructured spectral/hp element methods to the Newtonian incompressible Navier–Stokes equations in a geometrically simplified model of an aorto-intercostal junction. The effect of pulsatile but non-reversing side-branch flow was small; the aortic WSS pattern resembled that obtained under steady flow conditions, with high WSS upstream and downstream of the branch. When flow in the side branch or in the aortic near-wall region reversed during part of the cycle, significantly different instantaneous patterns were generated, with low WSS appearing upstream and downstream. Time-averaged WSS was similar to the steady flow case, reflecting the short duration of these events, but patterns of the oscillatory shear index for reversing aortic near-wall flow were profoundly altered. Effects of reverse flow may help explain the different distributions of lesions.


2021 ◽  
Vol 7 (9) ◽  
pp. 1546-1556
Author(s):  
Hatim El Assad ◽  
Benaissa Kissi ◽  
Rhanim Hassan ◽  
Parron Vera Miguel Angel ◽  
Rubio Cintas Maria Dolores ◽  
...  

In the area of civil engineering and especially hydraulic structures, we find multiple anomalies that weakens mechanical characteristics of dikes, one of the most common anomalies is erosion phenomenon specifically pipe flow erosion which causes major damage to dam structures. This phenomenon is caused by a hole which is the result of the high pressure of water that facilitate the soil migration between the two sides of the dam. It becomes only a question of time until the diameter of the hole expands and causes destruction of the dam structure. This problem pushed physicist to perform many tests to quantify erosion kinetics, one of the most used tests to have logical and trusted results is the HET (hole erosion test). Meanwhile there is not much research regarding the models that govern these types of tests. Objectives: In this paper we modeled the HET using modeling software based on the Navier Stokes equations, this model tackles also the singularity of the interface structure/water using wall laws for a flow turbulence. Methods/Analysis: The studied soil in this paper is a clay soil, clay soil has the property of containing water more than most other soils. Three wall laws were applied on the soil / water interface to calculate the erosion rate in order to avoid the rupture of such a structure. The modlisitation was made on the ANSYS software. Findings: In this work, two-dimensional modeling was carried of the soil.in contrast of the early models which is one-dimensional model, the first one had shown that the wall-shear stress which is not uniform along the whole wall. Then using the linear erosion law to predict the non-uniform erosion along the whole length. The previous study found that the wall laws have a significant impact on the wall-shear stress, which affects the erosion interface in the fluid/soil, particularly at the hole's extremes. Our experiment revealed that the degraded profile is not uniform. Doi: 10.28991/cej-2021-03091742 Full Text: PDF


Author(s):  
Katsuaki Shirai ◽  
Keisuke Tsuru ◽  
Shinnosuke Obi

We conducted a performance prediction for an optical wall shear stress sensor with using the velocity data of a direct numerical simulation. The Doppler signals were generated with respect to the path of tracer particles passing through the measurement volume. A signal processing technique was proposed to estimate the magnitude and yaw angle of local wall shear stress simultaneously from each Doppler signal. The simulated Doppler signals were processed with the technique, however the accuracy of estimating the yaw angle is not sufficient. In contrast, the estimated magnitude of wall shear stress showed a good agreement with the direct estimate from the DNS data if the yaw angle was accurately estimated. The measurement accuracy of the sensor mainly depends on estimating the yaw angle of each tracer particle. Another technique for detecting the yaw angle is needed for the accurate measurement of both the yaw angle and magnitude of local wall shear stress.


1993 ◽  
Vol 115 (4A) ◽  
pp. 412-417 ◽  
Author(s):  
Masahide Nakamura ◽  
Wataru Sugiyama ◽  
Manabu Haruna

An experiment on the fully developed sinusoidal pulsatile flow at transitional Reynolds numbers was performed to evaluate the basic characteristics of the wall shear stress. In this experiment, the wall shear stress was calculated from the measured section averaged axial velocity and the pressure gradient by using the section averaged Navier-Stokes equation. The experimental results showed that the ratio of the amplitude of the wall shear stress to the amplitude of the pressure gradient had the maximum value when the time averaged Reynolds number was about 4000 and the Womersley number was about 10. As this condition is close to the blood flow condition in the human aorta, it is suggested that the parameter of the aorta has an effect to increase the amplitude of the wall shear stress acting on the arterial wall.


1992 ◽  
Vol 114 (4) ◽  
pp. 515-520 ◽  
Author(s):  
L. H. Back ◽  
D. W. Crawford

Wall shear stress estimates from laminar boundary layer theory were found to agree fairly well with the magnitude of shear stress levels along coronary artery constrictions obtained from solutions of the Navier Stokes equations for both steady and pulsatile flow. The relatively simple method can be used for in vivo estimates of wall shear stress in constrictions by using a vessel shape function determined from a coronary angiogram, along with a knowledge of the flow rate.


2011 ◽  
Vol 64 (2) ◽  
Author(s):  
Giancarlo Alfonsi

The direct numerical simulation of turbulence (DNS) has become a method of outmost importance for the investigation of turbulence physics, and its relevance is constantly growing due to the increasing popularity of high-performance-computing techniques. In the present work, the DNS approach is discussed mainly with regard to turbulent shear flows of incompressible fluids with constant properties. A body of literature is reviewed, dealing with the numerical integration of the Navier-Stokes equations, results obtained from the simulations, and appropriate use of the numerical databases for a better understanding of turbulence physics. Overall, it appears that high-performance computing is the only way to advance in turbulence research through the front of the direct numerical simulation.


1951 ◽  
Vol 18 (1) ◽  
pp. 95-100
Author(s):  
Donald Ross ◽  
J. M. Robertson

Abstract As an interim solution to the problem of the turbulent boundary layer in an adverse pressure gradient, a super-position method of analysis has been developed. In this method, the velocity profile is considered to be the result of two effects: the wall shear stress and the pressure recovery. These are superimposed, yielding an expression for the velocity profiles which approximate measured distributions. The theory also leads to a more reasonable expression for the wall shear-stress coefficient.


2019 ◽  
Vol 878 ◽  
pp. 5-36 ◽  
Author(s):  
Yuji Hattori ◽  
Francisco J. Blanco-Rodríguez ◽  
Stéphane Le Dizès

The linear instability of a vortex ring with swirl with Gaussian distributions of azimuthal vorticity and velocity in its core is studied by direct numerical simulation. The numerical study is carried out in two steps: first, an axisymmetric simulation of the Navier–Stokes equations is performed to obtain the quasi-steady state that forms a base flow; then, the equations are linearized around this base flow and integrated for a sufficiently long time to obtain the characteristics of the most unstable mode. It is shown that the vortex rings are subjected to curvature instability as predicted analytically by Blanco-Rodríguez & Le Dizès (J. Fluid Mech., vol. 814, 2017, pp. 397–415). Both the structure and the growth rate of the unstable modes obtained numerically are in good agreement with the analytical results. However, a small overestimation (e.g. 22 % for a curvature instability mode) by the theory of the numerical growth rate is found for some instability modes. This is most likely due to evaluation of the critical layer damping which is performed for the waves on axisymmetric line vortices in the analysis. The actual position of the critical layer is affected by deformation of the core due to the curvature effect; as a result, the damping rate changes since it is sensitive to the position of the critical layer. Competition between the curvature and elliptic instabilities is also investigated. Without swirl, only the elliptic instability is observed in agreement with previous numerical and experimental results. In the presence of swirl, sharp bands of both curvature and elliptic instabilities are obtained for $\unicode[STIX]{x1D700}=a/R=0.1$, where $a$ is the vortex core radius and $R$ the ring radius, while the elliptic instability dominates for $\unicode[STIX]{x1D700}=0.18$. New types of instability mode are also obtained: a special curvature mode composed of three waves is observed and spiral modes that do not seem to be related to any wave resonance. The curvature instability is also confirmed by direct numerical simulation of the full Navier–Stokes equations. Weakly nonlinear saturation and subsequent decay of the curvature instability are also observed.


Sign in / Sign up

Export Citation Format

Share Document