Perturbation theory and numerical modelling of weakly and moderately nonlinear dynamics of the incompressible Richtmyer–Meshkov instability

2014 ◽  
Vol 751 ◽  
pp. 432-479 ◽  
Author(s):  
A. L. Velikovich ◽  
M. Herrmann ◽  
S. I. Abarzhi

AbstractA study of incompressible two-dimensional (2D) Richtmyer–Meshkov instability (RMI) by means of high-order perturbation theory and numerical simulations is reported. Nonlinear corrections to Richtmyer’s impulsive formula for the RMI bubble and spike growth rates have been calculated for arbitrary Atwood number and an explicit formula has been obtained for it in the Boussinesq limit. Conditions for early-time acceleration and deceleration of the bubble and the spike have been elucidated. Theoretical time histories of the interface curvature at the bubble and spike tip and the profiles of vertical and horizontal velocities have been calculated and favourably compared to simulation results. In our simulations we have solved 2D unsteady Navier–Stokes equations for immiscible incompressible fluids using the finite volume fractional step flow solver NGA developed by Desjardins et al. (J. Comput. Phys., vol. 227, 2008, pp. 7125–7159) coupled to the level set based interface solver LIT (Herrmann, J. Comput. Phys., vol. 227, 2008, pp. 2674–2706). We study the impact of small amounts of viscosity on the flow dynamics and compare simulation results to theory to discuss the influence of the theory’s ideal inviscid flow assumption.

2020 ◽  
Vol 67 ◽  
pp. 100-119 ◽  
Author(s):  
Laurent Boudin ◽  
Céline Grandmont ◽  
Bérénice Grec ◽  
Sébastien Martin ◽  
Amina Mecherbet ◽  
...  

In this paper, we propose a coupled fluid-kinetic model taking into account the radius growth of aerosol particles due to humidity in the respiratory system. We aim to numerically investigate the impact of hygroscopic effects on the particle behaviour. The air flow is described by the incompressible Navier-Stokes equations, and the aerosol by a Vlasov-type equation involving the air humidity and temperature, both quantities satisfying a convection-diffusion equation with a source term. Conservations properties are checked and an explicit time-marching scheme is proposed. Twodimensional numerical simulations in a branched structure show the influence of the particle size variations on the aerosol dynamics.


Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 168 ◽  
Author(s):  
Agostino Lauria ◽  
Giancarlo Alfonsi ◽  
Ali Tafarojnoruz

Ski jump spillways are frequently implemented to dissipate energy from high-speed flows. The general feature of this structure is to transform the spillway flow into a free jet up to a location where the impact of the jet creates a plunge pool, representing an area for potential erosion phenomena. In the present investigation, several tests with different ski jump bucket angles are executed numerically by means of the OpenFOAM® digital library, taking advantage of the Reynolds-averaged Navier–Stokes equations (RANS) approach. The results are compared to those obtained experimentally by other authors as related to the jet length and shape, obtaining physical insights into the jet characteristics. Particular attention is given to the maximum pressure head at the tailwater. Simple equations are proposed to predict the maximum dynamic pressure head acting on the tailwater, as dependent upon the Froude number, and the maximum pressure head on the bucket. Results of this study provide useful suggestions for the design of ski jump spillways in dam construction.


2020 ◽  
Vol 8 (2) ◽  
pp. 87 ◽  
Author(s):  
Paran Pourteimouri ◽  
Kourosh Hejazi

An integrated two-dimensional vertical (2DV) model was developed to investigate wave interactions with permeable submerged breakwaters. The integrated model is capable of predicting the flow field in both surface water and porous media on the basis of the extended volume-averaged Reynolds-averaged Navier–Stokes equations (VARANS). The impact of porous medium was considered by the inclusion of the additional terms of drag and inertia forces into conventional Navier–Stokes equations. Finite volume method (FVM) in an arbitrary Lagrangian–Eulerian (ALE) formulation was adopted for discretization of the governing equations. Projection method was utilized to solve the unsteady incompressible extended Navier–Stokes equations. The time-dependent volume and surface porosities were calculated at each time step using the fraction of a grid open to water and the total porosity of porous medium. The numerical model was first verified against analytical solutions of small amplitude progressive Stokes wave and solitary wave propagation in the absence of a bottom-mounted barrier. Comparisons showed pleasing agreements between the numerical predictions and analytical solutions. The model was then further validated by comparing the numerical model results with the experimental measurements of wave propagation over a permeable submerged breakwater reported in the literature. Good agreements were obtained for the free surface elevations at various spatial and temporal scales, velocity fields around and inside the obstacle, as well as the velocity profiles.


1983 ◽  
Vol 50 (2) ◽  
pp. 265-269
Author(s):  
D. Nixon

The perturbation theory for transonic flow is further developed for solutions of the Navier-Stokes equations in two dimensions or for experimental results. The strained coordinate technique is used to treat changes in location of any shock waves or large gradients.


2020 ◽  
Vol 8 (11) ◽  
pp. 903
Author(s):  
Sixtine Neuvéglise ◽  
Gaële Perret ◽  
Hassan Smaoui ◽  
François Marin ◽  
Philippe Sergent

This paper studies the behaviour of a quayside floater oscillating in front of a vertical dike. In order to study the floater motion and the impact of the dike on the floater, a linear analytical model based on 2D potential flow theory in intermediate water depth conditions and a numerical model resolving 2D Navier–Stokes equations are developed. Physical tests performed for different floater dimensions in a wave tank are used as references for the analytical and numerical models. The comparison of the results obtained analytically, numerically and experimentally leads to the validity domain of the potential model. A correction of this model is proposed, based on the optimization of the radiated coefficients, and a quadratic drag term is added according to Morison equation. The impact of the different parameters of the system on floater behaviour is considered. Results show that the draft has the most important impact on floater motion.


1985 ◽  
Vol 154 ◽  
pp. 357-375 ◽  
Author(s):  
J. A. C. Humphrey ◽  
H. Iacovides ◽  
B. E. Launder

The paper reports numerical solutions to a semi-elliptic truncation of the Navier–Stokes equations for the case of developing laminar flow in circular-sectioned bends over a range of Dean numbers. The ratios of bend radius to pipe radius are 7:1 and 20:1, corresponding with the configurations examined experimentally by Talbot and his co-workers in recent years. The semi-elliptic treatment facilitates a much finer grid than has been possible in earlier studies. Numerical accuracy has been further improved by assuming radial equilibrium over a thin sublayer immediately adjacent to the wall and by re-formulating the boundary conditions at the pipe centre.Streamwise velocity profiles at Dean numbers of 183 and 565 are in excellent agreement with laser-Doppler measurements by Agrawal, Talbot & Gong (1978). Good, albeit less complete, accord is found with the secondary velocities, though the differences that exist may be mainly due to the difficulty of making these measurements. The paper provides new information on the behaviour of the streamwise shear stress around the inner line of symmetry. Upstream of the point of minimum shear stress, our numerical predictions display a progressive shift towards the result of Stewartson, Cebici & Chang (1980) as the Dean number is successively raised. Downstream of the minimum, however, in contrast with the monotonic approach to an asymptotic level reported by Stewartson, the numerical solutions display a damped oscillatory behaviour reminiscent of those from Hawthorne's (1951) inviscid-flow calculations. The amplitude of the oscillation grows as the Dean number is raised.


1997 ◽  
Vol 339 ◽  
pp. 199-211 ◽  
Author(s):  
R. S. HEEG ◽  
N. RILEY

In this paper we present the results from numerical calculations, based upon the Navier–Stokes equations at relatively high Reynolds number, of the formation of a vortex ring when fluid is ejected from a circular tube. Our results are compared with the experiments of Didden (1979), and the inviscid flow calculations of Nitsche & Krasny (1994). Reasonable agreement is achieved except for the rate of shedding of circulation during the initial stages of ring formation. The theoretically predicted rate of shedding is substantially higher than that predicted by Didden. By contrast the inviscid theory predicts an anomalously high rate of initial shedding. We offer explanations for both of these apparent discrepancies.


1995 ◽  
Vol 303 ◽  
pp. 215-232 ◽  
Author(s):  
H. M. Badr ◽  
S. C. R. Dennis ◽  
S. Kocabiyik ◽  
P. Nguyen

The transient flow field caused by an infinitely long circular cylinder placed in an unbounded viscous fluid oscillating in a direction normal to the cylinder axis, which is at rest, is considered. The flow is assumed to be started suddenly from rest and to remain symmetrical about the direction of motion. The method of solution is based on an accurate procedure for integrating the unsteady Navier–Stokes equations numerically. The numerical method has been carried out for large values of time for both moderate and high Reynolds numbers. The effects of the Reynolds number and of the Strouhal number on the laminar symmetric wake evolution are studied and compared with previous numerical and experimental results. The time variation of the drag coefficients is also presented and compared with an inviscid flow solution for the same problem. The comparison between viscous and inviscid flow results shows a better agreement for higher values of Reynolds and a Strouhal numbers. The mean flow for large times is calculated and is found to be in good agreement with previous predictions based on boundary-layer theory.


Author(s):  
Hadi Karrabi ◽  
Mohsen Rezasoltani

An investigation to understand the impact of twisted, leaned and bowed blades on the performance of axial turbine was undertaken. A CFD code, which solves the Reynolds-averaged Navier–Stokes equations, was used to compute the complex flow field of axial turbine. The code was validated against existing Hannover turbine experimental data. Numerical data showed good agreement with measured data. Finally, the effect of geometry changes, focusing on blade lean, twist and bow, on the Avon turbine blade performance, was analyzed. Results show that twisted blade affects performance significantly. Leaned and bowed blade has minor effect on performance.


2013 ◽  
Vol 378 ◽  
pp. 418-423
Author(s):  
Gang Liu ◽  
Jia Wu ◽  
Wei Li

The three-dimensional construct of concentration field in an oscillatory flow reactor (OFR) containing periodically spaced conic ring baffles was investigated by numerical simulation employing Reynolds-averaged Navier-Stokes equations. The computation covered a range of Oscillatory Reynolds number (Reo) from 623.32 to 3116.58 at Strouhal number (St) 0.995 and 1.99. The contour of concentration field showed that the concentration in the most part of the channel is relative uniform and a small retention area is found below the conic ring baffles, which means a region of relative poor mixing. In addition, the turbulent diffusion coefficient calculated from the simulation results implied the greater oscillatory amplitude and oscillatory frequency superimposed to the fluid, the stronger is the turbulence intensity.


Sign in / Sign up

Export Citation Format

Share Document