A weakly nonlinear framework to study shock–vorticity interaction

2022 ◽  
Vol 933 ◽  
Author(s):  
Pranav Thakare ◽  
Vineeth Nair ◽  
Krishnendu Sinha
Keyword(s):  
Numerical Simulations ◽  
Acoustic Waves ◽  
Significant Contribution ◽  
Interaction Analysis ◽  
Nonlinear Effects ◽  
Normal Shock ◽  
Linear Interaction ◽  
Weakly Nonlinear ◽  
High Incidence ◽  
Vorticity Generation

Linear interaction analysis (LIA) is routinely used to study the shock–turbulence interaction in supersonic and hypersonic flows. It is based on the inviscid interaction of elementary Kovásznay modes with a shock discontinuity. LIA neglects nonlinear effects, and hence it is limited to small-amplitude disturbances. In this work, we extend the LIA framework to study the fundamental interaction of a two-dimensional vorticity wave with a normal shock. The predictions from a weakly nonlinear framework are compared with high-order accurate numerical simulations over a range of wave amplitudes ( $\epsilon$ ), incidence angles ( $\alpha$ ) and shock-upstream Mach numbers ( $M_1$ ). It is found that the nonlinear generation of vorticity at the shock has a significant contribution from the intermodal interaction between vorticity and acoustic waves. Vorticity generation is also strongly influenced by the curvature of the normal shock wave, especially for high incidence angles. Further, the weakly nonlinear analysis is able to predict the correct scaling of the nonlinear effects observed in the numerical simulations. The analysis also predicts a Mach number dependent limit for the validity of LIA in terms of the maximum possible amplitude of the upstream vorticity wave.

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 Numerical study of variable density turbulence interaction with a normal shock wave

2017 ◽  
Vol 829 ◽  
pp. 551-588 ◽  
Author(s):  
Yifeng Tian ◽  
Farhad A. Jaberi ◽  
Zhaorui Li ◽  
Daniel Livescu
Keyword(s):  
Shock Wave ◽  
Numerical Simulations ◽  
Numerical Study ◽  
Variable Density ◽  
Computational Method ◽  
Normal Shock ◽  
Turbulence Structure ◽  
Linear Interaction ◽  
Normal Shock Wave ◽  
Turbulent Statistics

Accurate numerical simulations of shock–turbulence interaction (STI) are conducted with a hybrid monotonicity-preserving–compact-finite-difference scheme for a detailed study of STI in variable density flows. Theoretical and numerical assessments of data confirm that all turbulence scales as well as the STI are well captured by the computational method. Linear interaction approximation (LIA) convergence tests conducted with the shock-capturing simulations exhibit a similar trend of converging to LIA predictions to shock-resolving direct numerical simulations (DNS). The effects of density variations on STI are studied by comparing the results corresponding to an upstream multi-fluid mixture with the single-fluid case. The results show that for the current parameter ranges, the turbulence amplification by the normal shock wave is much higher and the reduction in turbulence length scales is more significant when strong density variations exist. Turbulent mixing enhancement by the shock is also increased and stronger mixing asymmetry in the postshock region is observed when there is significant density variation. The turbulence structure is strongly modified by the shock wave, with a differential distribution of turbulent statistics in regions having different densities. The dominant mechanisms behind the variable density STI are identified by analysing the transport equations for the Reynolds stresses, vorticity, normalized mass flux and density specific volume covariance.

Linear interaction of two-dimensional free-stream disturbances with an oblique shock wave

2019 ◽  
Vol 873 ◽  
pp. 1179-1205 ◽  
Author(s):  
Zhangfeng Huang ◽  
Huilin Wang
Keyword(s):  
Shock Wave ◽  
Plane Wave ◽  
Transmission Coefficient ◽  
Acoustic Waves ◽  
Interaction Analysis ◽  
Oblique Shock Wave ◽  
Theoretical Problem ◽  
Linear Interaction ◽  
Incident Plane ◽  
Excited Wave

The problem of interaction between disturbances and shock waves was solved by a theoretical approach called linear interaction analysis in the mid-twentieth century. More recently, great progress has been made in analysing shock–turbulence interactions by direct numerical simulation. However, an unsolved theoretical problem remains: What happens when no acoustic waves are stimulated behind the shock wave? The concept of a damped wave is introduced, which is a type of excited plane wave. Based on this, the dispersion and amplitude relationships between any incident plane wave and resulting stimulated waves are constructed analytically, systematically and comprehensively. The physical essence of damped waves and the existence of critical angles are clarified. It is demonstrated that a damped wave is a complex number space solution to the acoustic dispersion relationship under certain conditions. It acts as a bridge connecting fast and slow acoustic waves at the position where the $x$ component of the group velocity is zero. There are two critical angles that can excite fast and slow acoustic waves, which determine the conditions that stimulate a damped wave. Our results show good agreement with theoretical and simulation results. The contribution of each excited wave to the transmission coefficient is evaluated, the distribution of the transmission coefficient is analysed and application to an engineering wedge model is performed.

Simulation of multiwavelength conditions in laser picosecond ultrasonics

SIMULATION ◽  
2021 ◽  
pp. 003754972199645
Author(s):  
Philippe Babilotte
Keyword(s):  
Numerical Simulations ◽  
Acoustic Signal ◽  
Acoustic Waves ◽  
Bulk Material ◽  
Deformation Field ◽  
Optical Pump ◽  
Pump And Probe ◽  
Oriented Object ◽  
Numerical Approaches ◽  
Picosecond Ultrasonics

Complete numerical simulations are given under SciLab® and MATLAB® coding environments, concerning propagative acoustic wavefronts, for laser picosecond ultrasonics under multiwavelength conditions. Simulations of the deformation field and its propagation into bulk material are given under different wavelength configurations for optical pump and probe beams, which are used to generate and to detect the acoustic signal. Complete insights concerning the dynamics of the acoustic waves are given, considering the absence of carrier diffusions into the material. Several numerical approaches are proposed concerning both the functions introduced to simulate the wavefront ( Heaviside or error) and the coding approach (linear/vectorized/ Oriented Object Programming), under the pure thermo-elastic approach.

Predicting Thermoacoustic Instability in an Industrial Gas Turbine Combustor: Combining a Low Order Network Model With Flame LES

2017 ◽  
Author(s):  
Y. Xia ◽  
A. S. Morgans ◽  
W. P. Jones ◽  
J. Rogerson ◽  
G. Bulat ◽  
...  
Keyword(s):  
Gas Turbine ◽  
Acoustic Waves ◽  
Gas Turbine Combustor ◽  
Network Modelling ◽  
Weakly Nonlinear ◽  
Thermoacoustic Instability ◽  
Reduced Chemistry ◽  
Flame Describing Function ◽  
Industrial Gas Turbine ◽  
Low Order

The thermoacoustic modes of a full scale industrial gas turbine combustor have been predicted numerically. The predictive approach combines low order network modelling of the acoustic waves in a simplified geometry, with a weakly nonlinear flame describing function, obtained from incompressible large eddy simulations of the flame region under upstream forced velocity perturbations, incorporating reduced chemistry mechanisms. Two incompressible solvers, each employing different numbers of reduced chemistry mechanism steps, are used to simulate the turbulent reacting flowfield to predict the flame describing functions. The predictions differ slightly between reduced chemistry approximations, indicating the need for more involved chemistry. These are then incorporated into a low order thermoacoustic solver to predict thermoacoustic modes. For the combustor operating at two different pressures, most thermoacoustic modes are predicted to be stable, in agreement with the experiments. The predicted modal frequencies are in good agreement with the measurements, although some mismatches in the predicted modal growth rates and hence modal stabilities are observed. Overall, these findings lend confidence in this coupled approach for real industrial gas turbine combustors.

scholarly journals Anomalous width variation of rarefactive ion acoustic solitary waves in the context of auroral plasmas

2004 ◽  
Vol 11 (2) ◽  
pp. 219-228 ◽  
Author(s):  
S. S. Ghosh ◽  
G. S. Lakhina
Keyword(s):  
Solitary Waves ◽  
Nonlinear Theory ◽  
Acoustic Waves ◽  
Large Amplitude ◽  
Magnetized Plasma ◽  
Amplitude Variation ◽  
Fully Nonlinear ◽  
Weakly Nonlinear ◽  
Ion Acoustic Solitary Waves ◽  
Ion Acoustic

Abstract. The presence of dynamic, large amplitude solitary waves in the auroral regions of space is well known. Since their velocities are of the order of the ion acoustic speed, they may well be considered as being generated from the nonlinear evolution of ion acoustic waves. However, they do not show the expected width-amplitude correlation for K-dV solitons. Recent POLAR observations have actually revealed that the low altitude rarefactive ion acoustic solitary waves are associated with an increase in the width with increasing amplitude. This indicates that a weakly nonlinear theory is not appropriate to describe the solitary structures in the auroral regions. In the present work, a fully nonlinear analysis based on Sagdeev pseudopotential technique has been adopted for both parallel and oblique propagation of rarefactive solitary waves in a two electron temperature multi-ion plasma. The large amplitude solutions have consistently shown an increase in the width with increasing amplitude. The width-amplitude variation profile of obliquely propagating rarefactive solitary waves in a magnetized plasma have been compared with the recent POLAR observations. The width-amplitude variation pattern is found to fit well with the analytical results. It indicates that a fully nonlinear theory of ion acoustic solitary waves may well explain the observed anomalous width variations of large amplitude structures in the auroral region.

scholarly journals New Explicit Solutions to the Fractional-Order Burgers’ Equation

2021 ◽  
Vol 2021 ◽  
pp. 1-11
Author(s):  
M. Hafiz Uddin ◽  
Mohammad Asif Arefin ◽  
M. Ali Akbar ◽  
Mustafa Inc
Keyword(s):  
Differential Equations ◽  
Rational Function ◽  
Acoustic Waves ◽  
Burgers Equation ◽  
Three Dimensional ◽  
Wall Friction ◽  
Weakly Nonlinear ◽  
Bubbly Liquids ◽  
Fractional Differential ◽  
Accuracy Of Results

The closed-form wave solutions to the time-fractional Burgers’ equation have been investigated by the use of the two variables G ′ / G , 1 / G -expansion, the extended tanh function, and the exp-function methods translating the nonlinear fractional differential equations (NLFDEs) into ordinary differential equations. In this article, we ascertain the solutions in terms of tanh , sech , sinh , rational function, hyperbolic rational function, exponential function, and their integration with parameters. Advanced and standard solutions can be found by setting definite values of the parameters in the general solutions. Mathematical analysis of the solutions confirms the existence of different soliton forms, namely, kink, single soliton, periodic soliton, singular kink soliton, and some other types of solitons which are shown in three-dimensional plots. The attained solutions may be functional to examine unidirectional propagation of weakly nonlinear acoustic waves, the memory effect of the wall friction through the boundary layer, bubbly liquids, etc. The methods suggested are direct, compatible, and speedy to simulate using algebraic computation schemes, such as Maple, and can be used to verify the accuracy of results.

scholarly journals “Extraordinary” modulation instability in optics and hydrodynamics

2021 ◽  
Vol 118 (14) ◽  
pp. e2019348118
Author(s):  
Guillaume Vanderhaegen ◽  
Corentin Naveau ◽  
Pascal Szriftgiser ◽  
Alexandre Kudlinski ◽  
Matteo Conforti ◽  
...  
Keyword(s):  
Wave Propagation ◽  
Stability Analysis ◽  
Linear Stability ◽  
Nonlinear Theory ◽  
Linear Stability Analysis ◽  
Water Waves ◽  
Modulation Instability ◽  
Nonlinear Effects ◽  
Weakly Nonlinear ◽  
Weakly Nonlinear Theory

The classical theory of modulation instability (MI) attributed to Bespalov–Talanov in optics and Benjamin–Feir for water waves is just a linear approximation of nonlinear effects and has limitations that have been corrected using the exact weakly nonlinear theory of wave propagation. We report results of experiments in both optics and hydrodynamics, which are in excellent agreement with nonlinear theory. These observations clearly demonstrate that MI has a wider band of unstable frequencies than predicted by the linear stability analysis. The range of areas where the nonlinear theory of MI can be applied is actually much larger than considered here.

scholarly journals Precessional instability of a fluid cylinder

2011 ◽  
Vol 666 ◽  
pp. 104-145 ◽  
Author(s):  
ROMAIN LAGRANGE ◽  
PATRICE MEUNIER ◽  
FRANÇOIS NADAL ◽  
CHRISTOPHE ELOY
Keyword(s):  
Reynolds Number ◽  
Nonlinear Effects ◽  
Intermittent Flow ◽  
Mean Flow ◽  
Weakly Nonlinear ◽  
Weakly Nonlinear Theory ◽  
Aspect Ratios ◽  
Image Velocimetry ◽  
Precessional Motion ◽  
Weakly Nonlinear Model

In this paper, the instability of a fluid inside a precessing cylinder is addressed theoretically and experimentally. The precessional motion forces Kelvin modes in the cylinder, which can become resonant for given precessional frequencies and cylinder aspect ratios. When the Reynolds number is large enough, these forced resonant Kelvin modes eventually become unstable. A linear stability analysis based on a triadic resonance between a forced Kelvin mode and two additional free Kelvin modes is carried out. This analysis allows us to predict the spatial structure of the instability and its threshold. These predictions are compared to the vorticity field measured by particle image velocimetry with an excellent agreement. When the Reynolds number is further increased, nonlinear effects appear. A weakly nonlinear theory is developed semi-empirically by introducing a geostrophic mode, which is triggered by the nonlinear interaction of a free Kelvin mode with itself in the presence of viscosity. Amplitude equations are obtained coupling the forced Kelvin mode, the two free Kelvin modes and the geostrophic mode. They show that the instability saturates to a fixed point just above threshold. Increasing the Reynolds number leads to a transition from a steady saturated regime to an intermittent flow in good agreement with experiments. Surprisingly, this weakly nonlinear model still gives a correct estimate of the mean flow inside the cylinder even far from the threshold when the flow is turbulent.

Finite-amplitude internal wavepacket dispersion and breaking

2001 ◽  
Vol 429 ◽  
pp. 343-380 ◽  
Author(s):  
BRUCE R. SUTHERLAND
Keyword(s):  
Numerical Simulations ◽  
Internal Waves ◽  
Large Amplitude ◽  
Finite Amplitude ◽  
Critical Amplitude ◽  
Fully Nonlinear ◽  
Weakly Nonlinear ◽  
The Stability ◽  
The Waves ◽  
Horizontal Wavelength

The evolution and stability of two-dimensional, large-amplitude, non-hydrostatic internal wavepackets are examined analytically and by numerical simulations. The weakly nonlinear dispersion relation for horizontally periodic, vertically compact internal waves is derived and the results are applied to assess the stability of weakly nonlinear wavepackets to vertical modulations. In terms of Θ, the angle that lines of constant phase make with the vertical, the wavepackets are predicted to be unstable if [mid ]Θ[mid ] < Θc, where Θc = cos−1 (2/3)1/2 ≃ 35.3° is the angle corresponding to internal waves with the fastest vertical group velocity. Fully nonlinear numerical simulations of finite-amplitude wavepackets confirm this prediction: the amplitude of wavepackets with [mid ]Θ[mid ] > Θc decreases over time; the amplitude of wavepackets with [mid ]Θ[mid ] < Θc increases initially, but then decreases as the wavepacket subdivides into a wave train, following the well-known Fermi–Pasta–Ulam recurrence phenomenon.If the initial wavepacket is of sufficiently large amplitude, it becomes unstable in the sense that eventually it convectively overturns. Two new analytic conditions for the stability of quasi-plane large-amplitude internal waves are proposed. These are qualitatively and quantitatively different from the parametric instability of plane periodic internal waves. The ‘breaking condition’ requires not only that the wave is statically unstable but that the convective instability growth rate is greater than the frequency of the waves. The critical amplitude for breaking to occur is found to be ACV = cot Θ (1 + cos2 Θ)/2π, where ACV is the ratio of the maximum vertical displacement of the wave to its horizontal wavelength. A second instability condition proposes that a statically stable wavepacket may evolve so that it becomes convectively unstable due to resonant interactions between the waves and the wave-induced mean flow. This hypothesis is based on the assumption that the resonant long wave–short wave interaction, which Grimshaw (1977) has shown amplifies the waves linearly in time, continues to amplify the waves in the fully nonlinear regime. Using linear theory estimates, the critical amplitude for instability is ASA = sin 2Θ/(8π2)1/2. The results of numerical simulations of horizontally periodic, vertically compact wavepackets show excellent agreement with this latter stability condition. However, for wavepackets with horizontal extent comparable with the horizontal wavelength, the wavepacket is found to be stable at larger amplitudes than predicted if Θ [lsim ] 45°. It is proposed that these results may explain why internal waves generated by turbulence in laboratory experiments are often observed to be excited within a narrow frequency band corresponding to Θ less than approximately 45°.

Sign in / Sign up

Export Citation Format

Share Document