scholarly journals Linear and nonlinear dynamics of pulsatile channel flow

2017 ◽  
Vol 815 ◽  
pp. 435-480 ◽  
Author(s):  
Benoît Pier ◽  
Peter J. Schmid
Keyword(s):  
Nonlinear Dynamics ◽  
Numerical Simulations ◽  
Flow Rate ◽  
Stokes Equations ◽  
Base Flow ◽  
Direct Numerical Simulations ◽  
Mean Flow ◽  
Stabilizing Effect ◽  
Linear And Nonlinear ◽  
Pulsation Cycle

The dynamics of small-amplitude perturbations, as well as the regime of fully developed nonlinear propagating waves, is investigated for pulsatile channel flows. The time-periodic base flows are known analytically and completely determined by the Reynolds number $Re$ (based on the mean flow rate), the Womersley number $Wo$ (a dimensionless expression of the frequency) and the flow-rate waveform. This paper considers pulsatile flows with a single oscillating component and hence only three non-dimensional control parameters are present. Linear stability characteristics are obtained both by Floquet analyses and by linearized direct numerical simulations. In particular, the long-term growth or decay rates and the intracyclic modulation amplitudes are systematically computed. At large frequencies (mainly $Wo\geqslant 14$), increasing the amplitude of the oscillating component is found to have a stabilizing effect, while it is destabilizing at lower frequencies; strongest destabilization is found for $Wo\simeq 7$. Whether stable or unstable, perturbations may undergo large-amplitude intracyclic modulations; these intracyclic modulation amplitudes reach huge values at low pulsation frequencies. For linearly unstable configurations, the resulting saturated fully developed finite-amplitude solutions are computed by direct numerical simulations of the complete Navier–Stokes equations. Essentially two types of nonlinear dynamics have been identified: ‘cruising’ regimes for which nonlinearities are sustained throughout the entire pulsation cycle and which may be interpreted as modulated Tollmien–Schlichting waves, and ‘ballistic’ regimes that are propelled into a nonlinear phase before subsiding again to small amplitudes within every pulsation cycle. Cruising regimes are found to prevail for weak base-flow pulsation amplitudes, while ballistic regimes are selected at larger pulsation amplitudes; at larger pulsation frequencies, however, the ballistic regime may be bypassed due to the stabilizing effect of the base-flow pulsating component. By investigating extended regions of a multi-dimensional parameter space and considering both two-dimensional and three-dimensional perturbations, the linear and nonlinear dynamics are systematically explored and characterized.

Evolution of enstrophy in shock/homogeneous turbulence interaction

2012 ◽  
Vol 707 ◽  
pp. 74-110 ◽  
Author(s):  
Krishnendu Sinha
Keyword(s):  
Shock Wave ◽  
Numerical Simulations ◽  
High Speed ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Interaction Analysis ◽  
Direct Numerical Simulations ◽  
Homogeneous Isotropic Turbulence ◽  
Mean Flow ◽  
Linear Interaction

AbstractInteraction of turbulent fluctuations with a shock wave plays an important role in many high-speed flow applications. This paper studies the amplification of enstrophy, defined as mean-square fluctuating vorticity, in homogeneous isotropic turbulence passing through a normal shock. Linearized Navier–Stokes equations written in a frame of reference attached to the unsteady shock wave are used to derive transport equations for the vorticity components. These are combined to obtain an equation that describes the evolution of enstrophy across a time-averaged shock wave. A budget of the enstrophy equation computed using results from linear interaction analysis and data from direct numerical simulations identifies the dominant physical mechanisms in the flow. Production due to mean flow compression and baroclinic torques are found to be the major contributors to the enstrophy amplification. Closure approximations are proposed for the unclosed correlations in the production and baroclinic source terms. The resulting model equation is integrated to obtain the enstrophy jump across a shock for a range of upstream Mach numbers. The model predictions are compared with linear theory results for varying levels of vortical and entropic fluctuations in the upstream flow. The enstrophy model is then cast in the form of$k$–$\epsilon $equations and used to compute the interaction of homogeneous isotropic turbulence with normal shocks. The results are compared with available data from direct numerical simulations. The equations are further used to propose a model for the amplification of turbulent viscosity across a shock, which is then applied to a canonical shock–boundary layer interaction. It is shown that the current model is a significant improvement over existing models, both for homogeneous isotropic turbulence and in the case of complex high-speed flows with shock waves.

scholarly journals Adjoint-based parametric sensitivity analysis for swirling M-flames

2018 ◽  
Vol 859 ◽  
pp. 516-542 ◽  
Author(s):  
Calum S. Skene ◽  
Peter J. Schmid
Keyword(s):  
Sensitivity Analysis ◽  
Frequency Response ◽  
Stokes Equations ◽  
Numerical Study ◽  
Computational Cost ◽  
Base Flow ◽  
Sensitivity Analyses ◽  
Perturbation Expansions ◽  
Mean Flow ◽  
Adjoint Solutions

A linear numerical study is conducted to quantify the effect of swirl on the response behaviour of premixed lean flames to general harmonic excitation in the inlet, upstream of combustion. This study considers axisymmetric M-flames and is based on the linearised compressible Navier–Stokes equations augmented by a simple one-step irreversible chemical reaction. Optimal frequency response gains for both axisymmetric and non-axisymmetric perturbations are computed via a direct–adjoint methodology and singular value decompositions. The high-dimensional parameter space, containing perturbation and base-flow parameters, is explored by taking advantage of generic sensitivity information gained from the adjoint solutions. This information is then tailored to specific parametric sensitivities by first-order perturbation expansions of the singular triplets about the respective parameters. Valuable flow information, at a negligible computational cost, is gained by simple weighted scalar products between direct and adjoint solutions. We find that for non-swirling flows, a mode with azimuthal wavenumber $m=2$ is the most efficiently driven structure. The structural mechanism underlying the optimal gains is shown to be the Orr mechanism for $m=0$ and a blend of Orr and other mechanisms, such as lift-up, for other azimuthal wavenumbers. Further to this, velocity and pressure perturbations are shown to make up the optimal input and output showing that the thermoacoustic mechanism is crucial in large energy amplifications. For $m=0$ these velocity perturbations are mainly longitudinal, but for higher wavenumbers azimuthal velocity fluctuations become prominent, especially in the non-swirling case. Sensitivity analyses are carried out with respect to the Mach number, Reynolds number and swirl number, and the accuracy of parametric gradients of the frequency response curve is assessed. The sensitivity analysis reveals that increases in Reynolds and Mach numbers yield higher gains, through a decrease in temperature diffusion. A rise in mean-flow swirl is shown to diminish the gain, with increased damping for higher azimuthal wavenumbers. This leads to a reordering of the most effectively amplified mode, with the axisymmetric ($m=0$) mode becoming the dominant structure at moderate swirl numbers.

Analytical and Semi-empirical Models of Tornadoes and Downbursts

2021 ◽  
Author(s):  
Djordje Romanic ◽  
Horia Hangan
Keyword(s):  
Numerical Simulations ◽  
Stokes Equations ◽  
Empirical Models ◽  
Navier Stokes ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
Boundary Layer Equations ◽  
Simplifying Assumptions ◽  
The Individual ◽  
Semi Empirical

Analytical and semi-empirical models are inexpensive to run and can complement experimental and numerical simulations for risk analysis-related applications. Some models are developed by employing simplifying assumptions in the Navier-Stokes equations and searching for exact, but many times inviscid solutions occasionally complemented by boundary layer equations to take surface effects into account. Other use simple superposition of generic, canonical flows for which the individual solutions are known. These solutions are then ensembled together by empirical or semi-empirical fitting procedures. Few models address turbulent or fluctuating flow fields, and all models have a series of constants that are fitted against experiments or numerical simulations. This chapter presents the main models used to provide primarily mean flow solutions for tornadoes and downbursts. The models are organized based on the adopted solution techniques, with an emphasis on their assumptions and validity.

Stability and dynamics of the laminar wake past a slender blunt-based axisymmetric body

2011 ◽  
Vol 676 ◽  
pp. 110-144 ◽  
Author(s):  
P. BOHORQUEZ ◽  
E. SANMIGUEL-ROJAS ◽  
A. SEVILLA ◽  
J. I. JIMÉNEZ-GONZÁLEZ ◽  
C. MARTÍNEZ-BAZÁN
Keyword(s):  
Stability Analysis ◽  
Numerical Simulations ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Velocity Ratio ◽  
Reynolds Numbers ◽  
Base Flow ◽  
Direct Numerical Simulations ◽  
The Body ◽  
Stream Velocity

We investigate the stability properties and flow regimes of laminar wakes behind slender cylindrical bodies, of diameter D and length L, with a blunt trailing edge at zero angle of attack, combining experiments, direct numerical simulations and local/global linear stability analyses. It has been found that the flow field is steady and axisymmetric for Reynolds numbers below a critical value, Recs (L/D), which depends on the length-to-diameter ratio of the body, L/D. However, in the range of Reynolds numbers Recs(L/D) < Re < Reco(L/D), although the flow is still steady, it is no longer axisymmetric but exhibits planar symmetry. Finally, for Re > Reco, the flow becomes unsteady due to a second oscillatory bifurcation which preserves the reflectional symmetry. In addition, as the Reynolds number increases, we report a new flow regime, characterized by the presence of a secondary, low frequency oscillation while keeping the reflectional symmetry. The results reported indicate that a global linear stability analysis is adequate to predict the first bifurcation, thereby providing values of Recs nearly identical to those given by the corresponding numerical simulations. On the other hand, experiments and direct numerical simulations give similar values of Reco for the second, oscillatory bifurcation, which are however overestimated by the linear stability analysis due to the use of an axisymmetric base flow. It is also shown that both bifurcations can be stabilized by injecting a certain amount of fluid through the base of the body, quantified here as the bleed-to-free-stream velocity ratio, Cb = Wb/W∞.

scholarly journals Direct numerical simulations of laminar and transitional flows in diverging pipes

2019 ◽  
Vol 30 (1) ◽  
pp. 75-92 ◽  
Author(s):  
Dhanush Vittal Shenoy ◽  
Mostafa Safdari Shadloo ◽  
Jorge Peixinho ◽  
Abdellah Hadjadj
Keyword(s):  
Velocity Profile ◽  
Numerical Simulations ◽  
Stokes Equations ◽  
Spectral Element ◽  
Skin Friction Coefficient ◽  
Finite Amplitude ◽  
Direct Numerical Simulations ◽  
Recirculation Region ◽  
Cross Sectional ◽  
Content Type

Purpose Fluid flows in pipes whose cross-sectional area are increasing in the stream-wise direction are prone to separation of the recirculation region. This paper aims to investigate such fluid flow in expansion pipe systems using direct numerical simulations. The flow in circular diverging pipes with different diverging half angles, namely, 45, 26, 14, 7.2 and 4.7 degrees, are considered. The flow is fed by a fully developed laminar parabolic velocity profile at its inlet and is connected to a long straight circular pipe at its downstream to characterise recirculation zone and skin friction coefficient in the laminar regime. The flow is considered linearly stable for Reynolds numbers sufficiently below natural transition. A perturbation is added to the inlet fully developed laminar velocity profile to test the flow response to finite amplitude disturbances and to characterise sub-critical transition. Design/methodology/approach Direct numerical simulations of the Navier–Stokes equations have been solved using a spectral element method. Findings It is found that the onset of disordered motion and the dynamics of the localised turbulence patch are controlled by the Reynolds number, the perturbation amplitude and the half angle of the pipe. Originality/value The authors clarify different stages of flow behaviour under the finite amplitude perturbations and shed more light to flow physics such as existence of Kelvin–Helmholtz instabilities as well as mechanism of turbulent puff shedding in diverging pipe flows.

The analysis and modelling of dilatational terms in compressible turbulence

1991 ◽  
Vol 227 ◽  
pp. 473-493 ◽  
Author(s):  
S. Sarkar ◽  
G. Erlebacher ◽  
M. Y. Hussaini ◽  
H. O. Kreiss
Keyword(s):  
Asymptotic Analysis ◽  
Numerical Simulations ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Mixing Layer ◽  
Direct Numerical Simulations ◽  
Homogeneous Turbulence ◽  
Compressible Turbulence ◽  
Quasi Equilibrium ◽  
Compressible Mixing Layer

It is shown that the dilatational terms that need to be modelled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of the compressible velocity field, which generates these dilatational terms, is explored by asymptotic analysis of the compressible Navier-Stokes equations in the case of homogeneous turbulence. The lowest-order equations for the compressible field are solved and explicit expressions for some of the associated one-point moments are obtained. For low Mach numbers, the compressible mode has a fast timescale relative to the incompressible mode. Therefore, it is proposed that, in moderate Mach number homogeneous turbulence, the compressible component of the turbulence is in quasi-equilibrium with respect to the incompressible turbulence. A non-dimensional parameter which characterizes this equilibrium structure of the compressible mode is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

scholarly journals Parametric forcing approach to rough-wall turbulent channel flow

2012 ◽  
Vol 712 ◽  
pp. 169-202 ◽  
Author(s):  
A. Busse ◽  
N. D. Sandham
Keyword(s):  
Numerical Simulations ◽  
Channel Flow ◽  
Stokes Equations ◽  
Rough Surfaces ◽  
Turbulent Channel Flow ◽  
Rough Wall ◽  
Force Term ◽  
Shape Functions ◽  
Mean Flow ◽  
The Mean

AbstractThe effects of rough surfaces on turbulent channel flow are modelled by an extra force term in the Navier–Stokes equations. This force term contains two parameters, related to the density and the height of the roughness elements, and a shape function, which regulates the influence of the force term with respect to the distance from the channel wall. This permits a more flexible specification of a rough surface than a single parameter such as the equivalent sand grain roughness. The effects of the roughness force term on turbulent channel flow have been investigated for a large number of parameter combinations and several shape functions by direct numerical simulations. It is possible to cover the full spectrum of rough flows ranging from hydraulically smooth through transitionally rough to fully rough cases. By using different parameter combinations and shape functions, it is possible to match the effects of different types of rough surfaces. Mean flow and standard turbulence statistics have been used to compare the results to recent experimental and numerical studies and a good qualitative agreement has been found. Outer scaling is preserved for the streamwise velocity for both the mean profile as well as its mean square fluctuations in all but extremely rough cases. The structure of the turbulent flow shows a trend towards more isotropic turbulent states within the roughness layer. In extremely rough cases, spanwise structures emerge near the wall and the turbulent state resembles a mixing layer. A direct comparison with the study of Ashrafian, Andersson & Manhart (Intl J. Heat Fluid Flow, vol. 25, 2004, pp. 373–383) shows a good quantitative agreement of the mean flow and Reynolds stresses everywhere except in the immediate vicinity of the rough wall. The proposed roughness force term may be of benefit as a wall model for direct and large-eddy numerical simulations in cases where the exact details of the flow over a rough wall can be neglected.

Direct numerical simulations of roughness-induced transition in supersonic boundary layers

2012 ◽  
Vol 693 ◽  
pp. 28-56 ◽  
Author(s):  
Suman Muppidi ◽  
Krishnan Mahesh
Keyword(s):  
Numerical Simulations ◽  
Shear Layer ◽  
Boundary Layers ◽  
Plate Boundary ◽  
Direct Numerical Simulations ◽  
Transition To Turbulence ◽  
Streamwise Vortices ◽  
Mean Flow ◽  
The Mean ◽  
Near Wall

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.

Mechanics of machine milking: II. The flow-rate pattern within single pulsation cycles

1966 ◽  
Vol 33 (2) ◽  
pp. 177-191 ◽  
Author(s):  
C. C. Thiel ◽  
P. A. Clough ◽  
D. R. Westgarth ◽  
D. N. Akam
Keyword(s):  
Flow Rate ◽  
Response Rate ◽  
Sphincter Muscle ◽  
Mean Flow ◽  
Considerable Time ◽  
Flow Rates ◽  
Machine Milking ◽  
Closing Force ◽  
The Mean ◽  
Pulsation Cycle

SummaryThe milk flowing during a single pulsation cycle was collected in a circle of contiguous cups which rotated in a chamber at 1 rev/pulsation cycle just below the end of the teatcup liner. The mean flow rate during the time taken for each collecting cup to pass under the milk stream was calculated and the flow-rate curve for the milkflow period of the pulsation cycle plotted. Flow rates were measured at 130, 97, 65, 32 and 16 c/min, and also after the pulsator had been stopped with the liner open for 0·5 min (0 pulsation).It was concluded from the series of flow-rate curves at the different pulsation rates that flow rate from the teat increased in about 0·05 sec to a steady value which continued for 0·5 sec or so, and then declined over a period of about 1·5 sec to a new constant value approximately equal to that shown after milk had flowed continuously from the teat for 0·5 min.These results suggest that once the pressure difference across the streak canal during milking forces the teat sphincter open a considerable time elapses before the muscle control system responds, and that a further much longer period elapses before the full closing force of the sphincter is exerted. Thus, it would appear that at pulsation rates of about 50 c/min and above, the streak canal is closed by pressure exerted on the teat by the closing liner, the sphincter muscle playing no active part because its response rate is slow compared with the pulsation rate. At lower pulsation rates the flow rate declines during each cycle because the sphincter muscle has time to exert a closing force to a greater or lesser extent depending on the duration of the milkflow period.

Oscillatory forcing of flow through porous media. Part 1. Steady flow

2002 ◽  
Vol 465 ◽  
pp. 213-235 ◽  
Author(s):  
D. R. GRAHAM ◽  
J. J. L. HIGDON
Keyword(s):  
Porous Media ◽  
Flow Rate ◽  
Stokes Equations ◽  
Nonlinear Flow ◽  
Navier Stokes ◽  
Full Time ◽  
Inertial Effects ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
The Mean

Oscillatory forcing of a porous medium may have a dramatic effect on the mean flow rate produced by a steady applied pressure gradient. The oscillatory forcing may excite nonlinear inertial effects leading to either enhancement or retardation of the mean flow. Here, in Part 1, we consider the effects of non-zero inertial forces on steady flows in porous media, and investigate the changes in the flow character arising from changes in both the strength of the inertial terms and the geometry of the medium. The steady-state Navier–Stokes equations are solved via a Galerkin finite element method to determine the velocity fields for simple two-dimensional models of porous media. Two geometric models are considered based on constricted channels and periodic arrays of circular cylinders. For both geometries, we observe solution multiplicity yielding both symmetric and asymmetric flow patterns. For the cylinder arrays, we demonstrate that inertial effects lead to anisotropy in the effective permeability, with the direction of minimum resistance dependent on the solid volume fraction. We identify nonlinear flow phenomena which might be exploited by oscillatory forcing to yield a net increase in the mean flow rate. In Part 2, we take up the subject of unsteady flows governed by the full time-dependent Navier–Stokes equations.

Sign in / Sign up

Export Citation Format

Share Document