scholarly journals Sustained gravity currents in a channel

2016 ◽  
Vol 798 ◽  
pp. 853-888 ◽  
Author(s):  
Andrew J. Hogg ◽  
Mohamad M. Nasr-Azadani ◽  
Marius Ungarish ◽  
Eckart Meiburg

Gravitationally driven motion arising from a sustained constant source of dense fluid in a horizontal channel is investigated theoretically using shallow-layer models and direct numerical simulations of the Navier–Stokes equations, coupled to an advection–diffusion model of the density field. The influxed dense fluid forms a flowing layer underneath the less dense fluid, which initially filled the channel, and in this study its speed of propagation is calculated; the outflux is at the end of the channel. The motion, under the assumption of hydrostatic balance, is modelled using a two-layer shallow-water model to account for the flow of both the dense and the overlying less dense fluids. When the relative density difference between the fluids is small (the Boussinesq regime), the governing shallow-layer equations are solved using analytical techniques. It is demonstrated that a variety of flow-field patterns are feasible, including those with constant height along the length of the current and those where the height varies continuously and discontinuously. The type of solution realised in any scenario is determined by the magnitude of the dimensionless flux issuing from the source and the source Froude number. Two important phenomena may occur: the flow may be choked, whereby the excess velocity due to the density difference is bounded and the height of the current may not exceed a determined maximum value, and it is also possible for the dense fluid to completely displace all of the less dense fluid originally in the channel in an expanding region close to the source. The onset and subsequent evolution of these types of motions are also calculated using analytical techniques. The same range of phenomena occurs for non-Boussinesq flows; in this scenario, the solutions of the model are calculated numerically. The results of direct numerical simulations of the Navier–Stokes equations are also reported for unsteady two-dimensional flows in which there is an inflow of dense fluid at one end of the channel and an outflow at the other end. These simulations reveal the detailed mechanics of the motion and the bulk properties are compared with the predictions of the shallow-layer model to demonstrate good agreement between the two modelling strategies.

Author(s):  
Giacomo Busco ◽  
Yassin A. Hassan

The highly turbulent flow inside a pressurized water reactor makes unpractical the use of scale resolving simulations, due to the large number of space and time turbulent structures. The high computational cost associated with typical large eddies simulations or direct numerical simulations techniques is unsuitable due to the large spatiotemporal resolution required. Partially averaged Navier-Stokes turbulence model is presented as bridging model between Reynolds averaged Navier-Stokes equations and direct numerical simulations. As filtered representation of the Navier-Stokes equations, the model is able to continuously shift its energy-based filter, inside the turbulence spectrum, being able to resolve the turbulent scales of interest. The choice of energy based cut-off filters gives the chance to directly impose the degree of needed resolution, where the most important large scales unsteadiness are resolved at minimal computational expenses. The partially averaged Navier-Stokes modelling approach has been tested for a Reynolds number of 14,000, inside a 5 × 5 fuel bundle, with a single spacer grid and split-type mixing vanes. Four different filters have been tested, whose resolution ranged from Reynolds averaged Navier-Stokes and large eddy simulation. A comparison with large eddy simulation will be presented. The results show that the partially averaged Navier-Stokes modeling produces results comparable to those of large eddy simulation when the appropriate cut-off energy filter is chosen. The turbulence models results will be compared with the available particle image velocimetry experimental data.


Volume 3 ◽  
2004 ◽  
Author(s):  
Gre´tar Tryggvason ◽  
Asghar Esmaeeli

Numerical simulations of boiling flows are discussed. The change of phase from liquid to vapor and vice-versa usually takes place in a highly unsteady manner where the phase boundary is very convoluted. Direct numerical simulations therefore require the accurate solution of the Navier-Stokes equations and the energy equation in each phase and the correct incorporation of the unsteady phase boundary. Such simulations, where the motion of an unsteady phase boundary is followed for a sufficiently long time to allow computation of average heat transfer are very recent. Here, we will describe one method that has been used successfully to simulate boiling flows and show a few examples of studies using the method.


2013 ◽  
Vol 716 ◽  
Author(s):  
M. Pradas ◽  
S. Kalliadasis ◽  
P.-K. Nguyen ◽  
V. Bontozoglou

AbstractWe examine pulse interaction and bound-state formation in interfacial turbulence using the problem of a falling liquid film as a model system. We perform direct numerical simulations (DNSs) of the full Navier–Stokes equations and associated wall and free-surface boundary conditions and we examine both analytically and numerically a low-dimensional (LD) model for the film. For a two-pulse system, DNSs reveal the existence of very rich and complex pulse interactions, characterized by either overdamped, underdamped or self-sustained oscillating behaviours, all of them found to be in excellent agreement with LD results. Having demonstrated the reliability of the LD model for two-pulse systems/smaller domains, we use it to investigate larger domains with many interacting pulses, where DNSs are computationally expensive. We demonstrate that such systems are likely to be dominated by a self-sustained oscillatory dynamics.


Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 756
Author(s):  
Federico Lluesma-Rodríguez ◽  
Francisco Álcantara-Ávila ◽  
María Jezabel Pérez-Quiles ◽  
Sergio Hoyas

One numerical method was designed to solve the time-dependent, three-dimensional, incompressible Navier–Stokes equations in turbulent thermal channel flows. Its originality lies in the use of several well-known methods to discretize the problem and its parallel nature. Vorticy-Laplacian of velocity formulation has been used, so pressure has been removed from the system. Heat is modeled as a passive scalar. Any other quantity modeled as passive scalar can be very easily studied, including several of them at the same time. These methods have been successfully used for extensive direct numerical simulations of passive thermal flow for several boundary conditions.


1974 ◽  
Vol 96 (4) ◽  
pp. 394-400 ◽  
Author(s):  
V. A. Marple ◽  
B. Y. H. Liu ◽  
K. T. Whitby

The flow field in an inertial impactor was studied experimentally with a water model by means of a flow visualization technique. The influence of such parameters as Reynolds number and jet-to-plate distance on the flow field was determined. The Navier-Stokes equations describing the laminar flow field in the impactor were solved numerically by means of a finite difference relaxation method. The theoretical results were found to be in good agreement with the empirical observations made with the water model.


Author(s):  
Djordje Romanic ◽  
Horia Hangan

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.


Author(s):  
Jaromi´r Hora´cˇek ◽  
Miloslav Feistauer ◽  
Petr Sva´cˇek

The contribution deals with the numerical simulation of the flutter of an airfoil with three degrees of freedom (3-DOF) for rotation around an elastic axis, oscillation in the vertical direction and rotation of a flap. The finite element (FE) solution of two-dimensional (2-D) incompressible Navier-Stokes equations is coupled with a system of nonlinear ordinary differential equations describing the airfoil vibrations with large amplitudes taking into account the nonlinear mass matrix. The time-dependent computational domain and a moving grid are treated by the Arbitrary Lagrangian-Eulerian (ALE) method and a suitable stabilization of the FE discretization is applied. The developed method was successfully tested by the classical flutter computation of the critical flutter velocity using NASTRAN program considering the linear model of vibrations and the double-lattice aerodynamic theory. The method was applied to the numerical simulations of the post flutter regime in time domain showing Limit Cycle Oscillations (LCO) due to nonlinearities of the flow model and vibrations with large amplitudes. Numerical experiments were performed for the airfoil NACA 0012 respecting the effect of the air space between the flap and the main airfoil.


Sign in / Sign up

Export Citation Format

Share Document