scholarly journals Deep-water sediment wave formation: linear stability analysis of coupled flow/bed interaction

2011 ◽  
Vol 680 ◽  
pp. 435-458 ◽  
Author(s):  
L. LESSHAFFT ◽  
B. HALL ◽  
E. MEIBURG ◽  
B. KNELLER
Keyword(s):  
Boundary Layer ◽  
Stability Analysis ◽  
Internal Wave ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Stokes Equations ◽  
Wave Mode ◽  
Bottom Current ◽  
Field Observations ◽  
Sediment Bed

A linear stability analysis is carried out for the interaction of an erodible sediment bed with a sediment-laden, stratified flow above the bed, such as a turbidity or bottom current. The fluid motion is described by the full, two-dimensional Navier–Stokes equations in the Boussinesq approximation, while erosion is modelled as a diffusive flux of particles from the bed into the fluid. The stability analysis shows the existence of both Tollmien–Schlichting and internal wave modes in the stratified boundary layer. For the internal wave mode, the stratified boundary layer acts as a wave duct, whose height can be determined analytically from the Brunt–Väisälä frequency criterion. Consistent with this criterion, distinct unstable perturbation wavenumber regimes exist for the internal wave mode, which are associated with different numbers of pressure extrema in the wall-normal direction. For representative turbidity current parameters, the analysis predicts unstable wavelengths that are consistent with field observations. As a key condition for instability to occur, the base flow velocity boundary layer needs to be thinner than the corresponding concentration boundary layer. For most of the unstable wavenumber ranges, the phase relations between the sediment bed deformation and the associated wall shear stress and concentration perturbations are such that the sediment waves migrate in the upstream direction, which again is consistent with field observations.

Variable density and viscosity, miscible displacements in horizontal Hele-Shaw cells. Part 1. Linear stability analysis

2013 ◽  
Vol 721 ◽  
pp. 268-294 ◽  
Author(s):  
L. Talon ◽  
N. Goyal ◽  
E. Meiburg
Keyword(s):  
Growth Rate ◽  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Stokes Equations ◽  
Variable Density ◽  
Local Concentration ◽  
Miscible Displacements ◽  
Instability Modes ◽  
Finger Tip

AbstractA computational investigation of variable density and viscosity, miscible displacements in horizontal Hele-Shaw cells is presented. As a first step, two-dimensional base states are obtained by means of simulations of the Stokes equations, which are nonlinear due to the dependence of the viscosity on the local concentration. Here, the vertical position of the displacement front is seen to reach a quasisteady equilibrium value, reflecting a balance between viscous and gravitational forces. These base states allow for two instability modes: first, there is the familiar tip instability driven by the unfavourable viscosity contrast of the displacement, which is modulated by the presence of density variations in the gravitational field; second, a gravitational instability occurs at the unstably stratified horizontal interface along the side of the finger. Both of these instability modes are investigated by means of a linear stability analysis. The gravitational mode along the side of the finger is characterized by a wavelength of about one half to one full gap width. It becomes more unstable as the gravity parameter increases, even though the interface is shifted closer to the wall. The growth rate is largest far behind the finger tip, where the interface is both thicker, and located closer to the wall, than near the finger tip. The competing influences of interface thickness and wall proximity are clarified by means of a parametric stability analysis. The tip instability mode represents a gravity-modulated version of the neutrally buoyant mode. The analysis shows that in the presence of density stratification its growth rate increases, while the dominant wavelength decreases. This overall destabilizing effect of gravity is due to the additional terms appearing in the stability equations, which outweigh the stabilizing effects of gravity onto the base state.

Linear Stability Analysis and Application of a New Solution Method of the Elastodynamic Equations Suitable for a Unified Fluid-Structure-Interaction Approach

2008 ◽  
Vol 130 (3) ◽  
Author(s):  
C. G. Giannopapa ◽  
G. Papadakis
Keyword(s):  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Solution Method ◽  
Stokes Equations ◽  
Fluid Structure Interaction ◽  
Analytic Solutions ◽  
Computational Domain ◽  
Fluid Structure ◽  
Structure Interaction

In the conventional approach for fluid-structure-interaction problems, the fluid and solid components are treated separately and information is exchanged across their interface. According to the conventional terminology, the current numerical methods can be grouped in two major categories: partitioned methods and monolithic methods. Both methods use separate sets of equations for fluid and solid that have different unknown variables. A unified solution method has been presented in the previous work of Giannopapa and Papadakis (2004, “A New Formulation for Solids Suitable for a Unified Solution Method for Fluid-Structure Interaction Problems,” ASME PVP 2004, San Diego, CA, July, PVP Vol. 491–1, pp. 111–117), which is different from these methods. The new approach treats both fluid and solid as a single continuum; thus, the whole computational domain is treated as one entity discretized on a single grid. Its behavior is described by a single set of equations, which are solved fully implicitly. In this paper, the elastodynamic equations are reformulated so that they contain the same unknowns as the Navier–Stokes equations, namely, velocities and pressure. Two time marching and one spatial discretization scheme, widely used for fluid equations, are applied for the solution of the reformulated equations for solids. Using linear stability analysis, the accuracy and dissipation characteristics of the resulting difference equations are examined. The aforementioned schemes are applied to a transient structural problem (beam bending) and the results compare favorably with available analytic solutions and are consistent with the conclusions of the stability analysis. A parametric investigation using different meshes, time steps, and beam dimensions is also presented. For all cases examined, the numerical solution was stable and robust and therefore is suitable for the next stage of application to full fluid-structure-interaction problems.

A Linear Stability Analysis of Incompressible Coaxial Jets Using an Accurate Boundary-Layer Approximation as Base Flow

2015 ◽  
Author(s):  
Helio Ricardo de Aguiar Quintanilha Júnior ◽  
Leonardo Santos de Brito Alves ◽  
Oberdan Miguel Rodrigues de Souza ◽  
Marcio Teixeira de Mendonça
Keyword(s):  
Boundary Layer ◽  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Base Flow ◽  
Boundary Layer Approximation ◽  
Coaxial Jets

Non-parallel linear stability analysis of the vertical boundary layer in a differentially heated cavity

1997 ◽  
Vol 352 ◽  
pp. 265-281 ◽  
Author(s):  
A. M. H. BROOKER ◽  
J. C. PATTERSON ◽  
S. W. ARMFIELD
Keyword(s):  
Boundary Layer ◽  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Navier Stokes ◽  
Parabolized Stability Equations ◽  
The Stability ◽  
Flow Quantity ◽  
Energy Equations ◽  
Made In

A non-parallel linear stability analysis which utilizes the assumptions made in the parabolized stability equations is applied to the buoyancy-driven flow in a differentially heated cavity. Numerical integration of the complete Navier–Stokes and energy equations is used to validate the non-parallel theory by introducing an oscillatory heat input at the upstream end of the boundary layer. In this way the stability properties are obtained by analysing the evolution of the resulting disturbances. The solutions show that the spatial growth rate and wavenumber are highly dependent on the transverse location and the disturbance flow quantity under consideration. The local solution to the parabolized stability equations accurately predicts the wave properties observed in the direct simulation whereas conventional parallel stability analysis overpredicts the spatial amplification and the wavenumber.

Sign in / Sign up

Export Citation Format

Share Document