Modelling Turbulent Flow in Deformable Highly Porous Seabed and Structures

Author(s):  
Hisham Elsafti ◽  
Hocine Oumeraci

In this study, the fully-coupled and fully-dynamic Biot governing equations in the open-source geotechFoam solver are extended to account for pore fluid viscous stresses. Additionally, turbulent pore fluid flow in deformable porous media is modeled by means of the conventional eddy viscosity concept without the need to resolve all turbulence scales. A new approach is presented to account for porous media resistance to flow (solid-to-fluid coupling) by means of an effective viscosity, which accounts for tortuosity, grain shape and local turbulences induced by flow through porous media. The new model is compared to an implemented extended Darcy-Forchheimer model in the Navier-Stokes equations, which accounts for laminar, transitional, turbulent and transient flow regimes. Further, to account for skeleton deformation, the porosity and other model parameters are updated with regard to strain of geomaterials. The presented model is calibrated by means of available results of physical experiments of unidirectional and oscillatory flows.

2016 ◽  
Vol 798 ◽  
pp. 643-681 ◽  
Author(s):  
E. D. Fernández-Nieto ◽  
J. Garres-Díaz ◽  
A. Mangeney ◽  
G. Narbona-Reina

In this work we present a multilayer shallow model to approximate the Navier–Stokes equations with the ${\it\mu}(I)$-rheology through an asymptotic analysis. The main advantages of this approximation are (i) the low cost associated with the numerical treatment of the free surface of the modelled flows, (ii) the exact conservation of mass and (iii) the ability to compute two-dimensional profiles of the velocities in the directions along and normal to the slope. The derivation of the model follows Fernández-Nieto et al. (J. Comput. Phys., vol. 60, 2014, pp. 408–437) and introduces a dimensional analysis based on the shallow flow hypothesis. The proposed first-order multilayer model fully satisfies a dissipative energy equation. A comparison with steady uniform Bagnold flow – with and without the sidewall friction effect – and laboratory experiments with a non-constant normal profile of the downslope velocity demonstrates the accuracy of the numerical model. Finally, by comparing the numerical results with experimental data on granular collapses, we show that the proposed multilayer model with the ${\it\mu}(I)$-rheology qualitatively reproduces the effect of the erodible bed on granular flow dynamics and deposits, such as the increase of runout distance with increasing thickness of the erodible bed. We show that the use of a constant friction coefficient in the multilayer model leads to the opposite behaviour. This multilayer model captures the strong change in shape of the velocity profile (from S-shaped to Bagnold-like) observed during the different phases of the highly transient flow, including the presence of static and flowing zones within the granular column.


2009 ◽  
Vol 36 (7) ◽  
pp. 1200-1206 ◽  
Author(s):  
Guillermo A. Narsilio ◽  
Olivier Buzzi ◽  
Stephen Fityus ◽  
Tae Sup Yun ◽  
David W. Smith

Geophysics ◽  
1988 ◽  
Vol 53 (4) ◽  
pp. 509-518 ◽  
Author(s):  
Daniel H. Rothman

Numerical models of fluid flow through porous media can be developed from either microscopic or macroscopic properties. The large‐scale viewpoint is perhaps the most prevalent. Darcy’s law relates the chief macroscopic parameters of interest—flow rate, permeability, viscosity, and pressure gradient—and may be invoked to solve for any of these parameters when the others are known. In practical situations, however, this solution may not be possible. Attention is then typically focused on the estimation of permeability, and numerous numerical methods based on knowledge of the microscopic pore‐space geometry have been proposed. Because the intrinsic inhomogeneity of porous media makes the application of proper boundary conditions difficult, microscopic flow calculations have typically been achieved with idealized arrays of geometrically simple pores, throats, and cracks. I propose here an attractive alternative which can freely and accurately model fluid flow in grossly irregular geometries. This new method solves the Navier‐Stokes equations numerically using the cellular‐automaton fluid model introduced by Frisch, Hasslacher, and Pomeau. The cellular‐ automaton fluid is extraordinarily simple—particles of unit mass traveling with unit velocity reside on a triangular lattice and obey elementary collision rules—but is capable of modeling much of the rich complexity of real fluid flow. Cellular‐automaton fluids are applicable to the study of porous media. In particular, numerical methods can be used to apply the appropriate boundary conditions, create a pressure gradient, and measure the permeability. Scale of the cellular‐automaton lattice is an important issue; the linear dimension of a void region must be approximately twice the mean free path of a lattice gas particle. Finally, an example of flow in a 2-D porous medium demonstrates not only the numerical solution of the Navier‐Stokes equations in a highly irregular geometry, but also numerical estimation of permeability and a verification of Darcy’s law.


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.


2002 ◽  
Vol 465 ◽  
pp. 213-235 ◽  
Author(s):  
D. R. GRAHAM ◽  
J. J. L. HIGDON

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.


2016 ◽  
Vol 26 (05) ◽  
pp. 1650086
Author(s):  
Tingting Tang ◽  
Zhiyong Li ◽  
J. M. McDonough ◽  
P. D. Hislop

In this paper, a discrete dynamical system (DDS) is derived from the generalized Navier–Stokes equations for incompressible flow in porous media via a Galerkin procedure. The main difference from the previously studied poor man’s Navier–Stokes equations is the addition of forcing terms accounting for linear and nonlinear drag forces of the medium — Darcy and Forchheimer terms. A detailed numerical investigation focusing on the bifurcation parameters due to these additional terms is provided in the form of regime maps, time series, power spectra, phase portraits and basins of attraction, which indicate system behaviors in agreement with expected physical fluid flow through porous media. As concluded from the previous studies, this DDS can be employed in subgrid-scale models of synthetic-velocity form for large-eddy simulation of turbulent flow through porous media.


2003 ◽  
Author(s):  
Tien-Chien Jen ◽  
Tuan-Zhou Yan ◽  
S. H. Chan

A three-dimensional computational model is developed to analyze fluid flow in a semi-porous channel. In order to understand the developing fluid flow and heat transfer process inside the semi-porous channels, the conventional Navier-Stokes equations for gas channel, and volume-averaged Navier-Stokes equations for porous media layer are adopted individually in this study. Conservation of mass, momentum and energy equations are solved numerically in a coupled gas and porous media domain in a channel using the vorticity-velocity method with power law scheme. Detailed development of axial velocity, secondary flow and temperature fields at various axial positions in the entrance region are presented. The friction factor and Nusselt number are presented as a function of axial position, and the effects of the size of porous media inside semi-porous channel are also analyzed in the present study.


Author(s):  
Kah Joon Yong ◽  
Max Meindl ◽  
Wolfgang Polifke ◽  
Camilo F. Silva

Abstract This study investigates the effect of partial acoustic reflection at inlet or outlet of a combustor on thermoacoustic stability. Parametric maps of the thermoacoustic spectrum are utilized for this purpose, which represent frequencies and growth rates of eigenmodes for a wide range of model parameters. It is found that a decrease of the acoustic reflection at the boundaries does not always imply an increase in the stability margin of the thermoacoustic system. As a matter of fact, a reduction in the acoustic reflection may sometimes destabilize a thermoacoustic mode. Additionally, we show that perturbed passive thermoacoustic modes may become ITA modes in the fully anechoic case. We briefly discuss the mode definitions ‘acoustic’ and ‘intrinsic’ commonly found in the literature. The computational analysis is based on a state-space formulation of the Linearized Navier-Stokes Equations (LNSE) with discontinuous Galerkin discretization. This approach allows to describe the thermoacoustic system as a linear combination of internal acoustics, flame dynamics and acoustic boundaries. Such a segregation grants a clear analysis of the respective effects of the individual subsystems on the general stability of the system, expressed in terms of adjoint-based eigenvalue sensitivity. The state-space formulation of the LNSE proposed in this paper offers a powerful and flexible framework to carry out thermoacoustic studies of combustors with arbitrary geometry and acoustic boundary conditions.


2010 ◽  
Vol 644 ◽  
pp. 97-106 ◽  
Author(s):  
B. U. FELDERHOF

The flow of a viscous compressible fluid in a circular tube generated by a sudden impulse at a point on the axis is studied on the basis of the linearized Navier–Stokes equations. A no-slip boundary condition is assumed to hold on the wall of the tube. Owing to the finite velocity of sound the flow behaviour differs qualitatively from that of an incompressible fluid. The flow velocity and the pressure disturbance at any fixed point different from the source point vanish at short time and decay at long times with a t−3/2 power law.


Sign in / Sign up

Export Citation Format

Share Document