scholarly journals Numerical solution of the Navier-Stokes equations for incompressible flow in porous media with free surface boundary

2005 ◽  
Vol 8 ◽  
pp. 225-231 ◽  
Author(s):  
Nenad JACIMOVIC ◽  
Takashi HOSODA ◽  
Kiyoshi KISHIDA ◽  
Marko IVETIC
1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.


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.


2010 ◽  
Vol 658 ◽  
pp. 33-62 ◽  
Author(s):  
XIN GUO ◽  
LIAN SHEN

Direct numerical simulation is performed for the interaction between a deformable free surface and the homogeneous isotropic turbulent flow underneath. The Navier–Stokes equations subject to fully nonlinear free-surface boundary conditions are simulated by using a pseudospectral method in the horizontal directions and a finite-difference method in the vertical direction. Statistically, steady turbulence is generated by using a linear forcing method in the bulk flow below. Through investigation of cases of different Froude and Weber numbers, the present study focuses on the effect of surface deformation of finite amplitude. It is found that the motion of the free surface is characterized by propagating waves and turbulence-generated surface roughness. Statistics of the turbulence field near the free surface are analysed in detail in terms of fluctuations of velocity, fluctuations of velocity gradients and strain rates and the energy budget for horizontal and vertical turbulent motions. Our results illustrate the effects of surface blockage and vanishing shear stress on the anisotropy of the flow field. Using conditional averaging analysis, it is shown that splats and antisplats play an essential role in energy inter-component exchange and vertical transport.


AIAA Journal ◽  
2000 ◽  
Vol 38 ◽  
pp. 1603-1614
Author(s):  
Martin Scholtysik ◽  
Bernhard Mueller ◽  
Torstein K. Fannelop

Author(s):  
B. Elie ◽  
G. Reliquet ◽  
P.-E. Guillerm ◽  
O. Thilleul ◽  
P. Ferrant ◽  
...  

This paper compares numerical and experimental results in the study of the resonance phenomenon which appears between two side-by-side fixed barges for different sea-states. Simulations were performed using SWENSE (Spectral Wave Explicit Navier-Stokes Equations) approach and results are compared with experimental data on two fixed barges with different headings and bilges. Numerical results, obtained using the SWENSE approach, are able to predict both the frequency and the magnitude of the RAO functions.


Sign in / Sign up

Export Citation Format

Share Document