The dynamics of breaking internal solitary waves on slopes

2014 ◽  
Vol 761 ◽  
pp. 360-398 ◽  
Author(s):  
Robert S. Arthur ◽  
Oliver B. Fringer
Keyword(s):  
Reynolds Number ◽  
Velocity Structure ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Mixing Efficiency ◽  
Streamwise Vorticity ◽  
Initial Wave ◽  
Reynolds Number Effects ◽  
Dimensional Simulation ◽  
Dimensional Domain

AbstractUsing direct numerical simulations (DNS), we investigate the structure and energetics of breaking internal waves on slopes. We employ a Navier–Stokes code in an idealized three-dimensional domain where an internal solitary wave of depression impinges upon a sloping bottom. Seven cases with varying initial wave amplitude and bathymetric slope, but constant wave Reynolds number $\mathit{Re}_{w}$ are considered. Volume-integrated values of dissipation and irreversible mixing are related to the density and velocity structure of the wave throughout the breaking process. The majority of dissipation (63 %) occurs along the no-slip bottom boundary. Most of the remaining dissipation (35 %) and nearly all irreversible mixing occurs in the interior after breaking, when density overturns are present at the interface. Breaking introduces three-dimensionality to the flow field that is driven by the lateral breakdown of density overturns and the lobe–cleft instability typical of gravity currents. The resulting longitudinal rolls (streamwise vorticity) increase dissipation by roughly 8 % and decrease irreversible mixing by roughly 20 % when compared with a similar two-dimensional simulation. The bulk mixing efficiency is shown to increase for larger and smaller values of the internal Iribarren number ${\it\xi}$, with a minimum for intermediate values of ${\it\xi}$ and a peak near ${\it\xi}=0.8$ for plunging breakers. This trend is explained by the degree of two-dimensionality in the flow, and agrees with previous results in the literature after accounting for Reynolds number effects. Local turbulence quantities are also calculated at ‘virtual moorings’, and a location upslope of the breakpoint but downslope of the intersection of the pycnocline and the bottom is shown to provide a signal that is most representative of the volume-integrated dissipation and mixing results.

Finite-Reynolds-Number Effects in Steady, Three-Dimensional Airway Reopening

2006 ◽  
Vol 128 (4) ◽  
pp. 573-578 ◽  
Author(s):  
Andrew L. Hazel ◽  
Matthias Heil
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Elastic Tube ◽  
Navier Stokes ◽  
Fluid Inertia ◽  
Navier Stokes Equations ◽  
Airway Reopening ◽  
Reynolds Number Effects ◽  
Steady Propagation

Motivated by the physiological problem of pulmonary airway reopening, we study the steady propagation of an air finger into a buckled elastic tube, initially filled with viscous fluid. The system is modeled using geometrically non-linear, Kirchhoff-Love shell theory, coupled to the free-surface Navier-Stokes equations. The resulting three-dimensional, fluid-structure-interaction problem is solved numerically by a fully coupled finite element method. Our study focuses on the effects of fluid inertia, which has been neglected in most previous studies. The importance of inertial forces is characterized by the ratio of the Reynolds and capillary numbers, Re∕Ca, a material parameter. Fluid inertia has a significant effect on the system’s behavior, even at relatively small values of Re∕Ca. In particular, compared to the case of zero Reynolds number, fluid inertia causes a significant increase in the pressure required to drive the air finger at a given speed.

Analysis and direct numerical simulation of the flow at a gravity-current head. Part 1. Flow topology and front speed for slip and no-slip boundaries

2000 ◽  
Vol 418 ◽  
pp. 189-212 ◽  
Author(s):  
CARLOS HÄRTEL ◽  
ECKART MEIBURG ◽  
FRIEDER NECKER
Keyword(s):  
Reynolds Number ◽  
Stagnation Point ◽  
Flow Structure ◽  
Gravity Current ◽  
Three Dimensional ◽  
Head Part ◽  
Two Dimensional ◽  
Flow Topology ◽  
Reynolds Number Effects ◽  
Dimensional Simulation

Direct numerical simulations are performed of gravity-current fronts in the lock-exchange configuration. The case of small density differences is considered, where the Boussinesq approximations can be adopted. The key objective of the investigation is a detailed analysis of the flow structure at the foremost part of the front, where no previous high-resolution data were available. For the simulations, high-order numerical methods are used, based on spectral and spectral-element discretizations and compact finite differences. A three-dimensional simulation is conducted of a front spreading along a no-slip boundary at a Reynolds number of about 750. The simulation exhibits all features typically observed in experimental flows near the gravity-current head, including the lobe-and-cleft structure at the leading edge. The results reveal that the flow topology at the head differs from what has been assumed previously, in that the foremost point is not a stagnation point in a translating system. Rather, the stagnation point is located below and slightly behind the foremost point in the vicinity of the wall. The relevance of this finding for the mechanism behind the lobe-and-cleft instability is discussed. In order to explore the high-Reynolds-number regime, and to assess potential Reynolds-number effects, two-dimensional simulations are conducted for Reynolds numbers up to about 30 000, for both no-slip and slip (i.e. shear-stress free) boundaries. It is shown that although quantitative Reynolds-number effects persist over the whole range examined, no qualitative changes in the flow structure at the head can be observed. A comparison of the two-dimensional results with laboratory data and the three-dimensional simulation provides evidence that a two-dimensional model is able to capture essential features of the flow at the head. The simulations also show that for the free-slip case the shape of the head agrees closely with the classical inviscid theory of Benjamin.

Direct simulation of evolution and control of three-dimensional instabilities in attachment-line boundary layers

1995 ◽  
Vol 291 ◽  
pp. 369-392 ◽  
Author(s):  
Ronald D. Joslin
Keyword(s):  
Boundary Layer ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Plate Boundary ◽  
Navier Stokes ◽  
Two Dimensional ◽  
Computationally Efficient ◽  
Simulation Results ◽  
Dimensional Simulation ◽  
Attachment Line

The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier–Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic-source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in flat-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.

The development of asymmetry and period doubling for oscillatory flow in baffled channels

1996 ◽  
Vol 328 ◽  
pp. 19-48 ◽  
Author(s):  
E. P. L. Roberts ◽  
M. R. Mackley
Keyword(s):  
Reynolds Number ◽  
Oscillatory Flow ◽  
Three Dimensional ◽  
Period Doubling ◽  
Progressive Increase ◽  
Newtonian Flow ◽  
Chaotic Flow ◽  
Spatially Periodic ◽  
Dimensional Simulation ◽  
Time Periodic

We report experimental and numerical observations on the way initially symmetric and time-periodic fluid oscillations in baffled channels develop in complexity. Experiments are carried out in a spatially periodic baffled channel with a sinusoidal oscillatory flow. At modest Reynolds number the observed vortex structure is symmetric and time periodic. At higher values the flow progressively becomes three-dimensional, asymmetric and aperiodic. A two-dimensional simulation of incompressible Newtonian flow is able to follow the flow pattern at modest oscillatory Reynolds number. At higher values we report the development of both asymmetry and a period-doubling cascade leading to a chaotic flow regime. A bifurcation diagram is constructed that can describe the progressive increase in complexity of the flow.

scholarly journals Computed Eccentricity Effects on Turbine Rim Seals at Engine Conditions With a Mainstream

1994 ◽  
Author(s):  
Zhou Guo ◽  
David L. Rhode ◽  
Fred M. Davis
Keyword(s):  
Reynolds Number ◽  
Temperature Rise ◽  
Three Dimensional ◽  
Computer Code ◽  
Interaction Effects ◽  
Navier Stokes ◽  
Radial Clearance ◽  
Turbine Wheel ◽  
Blade Root

A previously verified axisymmetric Navier-Stokes computer code was extended for three-dimensional computation of eccentric rim seals of almost any configuration. All compressibility and thermal/momentum interaction effects are completely, included, and the temperature, pressure and Reynolds number of the mainstream, coolant stream and turbine wheel are fixed at actual engine conditions. Regardless of the seal eccentricity, both ingress and egress are found between θ = −30° and 100°. which encompasses the location of maximum radial clearance at θ = 0°. All other θ locations within the rim seal show only egress, as does the concentric basecase for all circumferential locations. Further, the maximum ingress occurs near θ = 30° for all eccentricities. This is found to produce a blade root/retainer temperature rise from the concentric case of 390 percent at 50 percent eccentricity and a 77 percent rise at 7.5 percent eccentricity. In addition, the nature of an increased eccentricity causing a decreased seal effectiveness is examined, along with the corresponding increase of cavity-averaged temperature.

scholarly journals Nonlinear amplification in hydrodynamic turbulence

2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung
Keyword(s):  
Reynolds Number ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Developed Turbulence ◽  
Advection Term ◽  
Nonlinear Amplification

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.

scholarly journals U-shaped fairings suppress vortex-induced vibrations for cylinders in cross-flow

2015 ◽  
Vol 782 ◽  
pp. 300-332 ◽  
Author(s):  
Fangfang Xie ◽  
Yue Yu ◽  
Yiannis Constantinides ◽  
Michael S. Triantafyllou ◽  
George Em Karniadakis
Keyword(s):  
Reynolds Number ◽  
Drag Force ◽  
Three Dimensional ◽  
Cross Flow ◽  
Wake Flow ◽  
Streamwise Vorticity ◽  
Combined System ◽  
Force Coefficient ◽  
Vortex Induced Vibrations ◽  
Adjacent Segments

We employ three-dimensional direct and large-eddy numerical simulations of the vibrations and flow past cylinders fitted with free-to-rotate U-shaped fairings placed in a cross-flow at Reynolds number $100\leqslant \mathit{Re}\leqslant 10\,000$. Such fairings are nearly neutrally buoyant devices fitted along the axis of long circular risers to suppress vortex-induced vibrations (VIVs). We consider three different geometric configurations: a homogeneous fairing, and two configurations (denoted A and AB) involving a gap between adjacent segments. For the latter two cases, we investigate the effect of the gap on the hydrodynamic force coefficients and the translational and rotational motions of the system. For all configurations, as the Reynolds number increases beyond 500, both the lift and drag coefficients decrease. Compared to a plain cylinder, a homogeneous fairing system (no gaps) can help reduce the drag force coefficient by 15 % for reduced velocity $U^{\ast }=4.65$, while a type A gap system can reduce the drag force coefficient by almost 50 % for reduced velocity $U^{\ast }=3.5,4.65,6$, and, correspondingly, the vibration response of the combined system, as well as the fairing rotation amplitude, are substantially reduced. For a homogeneous fairing, the cross-flow amplitude is reduced by about 80 %, whereas for fairings with a gap longer than half a cylinder diameter, VIVs are completely eliminated, resulting in additional reduction in the drag coefficient. We have related such VIV suppression or elimination to the features of the wake flow structure. We find that a gap causes the generation of strong streamwise vorticity in the gap region that interferes destructively with the vorticity generated by the fairings, hence disorganizing the formation of coherent spanwise cortical patterns. We provide visualization of the incoherent wake flow that leads to total elimination of the vibration and rotation of the fairing–cylinder system. Finally, we investigate the effect of the friction coefficient between cylinder and fairing. The effect overall is small, even when the friction coefficients of adjacent segments are different. In some cases the equilibrium positions of the fairings are rotated by a small angle on either side of the centreline, in a symmetry-breaking bifurcation, which depends strongly on Reynolds number.

scholarly journals A Review of Passive Micromixers with a Comparative Analysis

Micromachines ◽  
2020 ◽  
Vol 11 (5) ◽  
pp. 455 ◽  
Author(s):  
Wasim Raza ◽  
Shakhawat Hossain ◽  
Kwang-Yong Kim
Keyword(s):  
Reynolds Number ◽  
Comparative Analysis ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Number Range ◽  
Pressure Drops ◽  
Convection Diffusion ◽  
Wide Range ◽  
Quantitative Analyses ◽  
Mixing Costs

A wide range of existing passive micromixers are reviewed, and quantitative analyses of ten typical passive micromixers were performed to compare their mixing indices, pressure drops, and mixing costs under the same axial length and flow conditions across a wide Reynolds number range of 0.01–120. The tested micromixers were selected from five types of micromixer designs. The analyses of flow and mixing were performed using continuity, Navier-Stokes and convection-diffusion equations. The results of the comparative analysis were presented for three different Reynolds number ranges: low-Re (Re ≤ 1), intermediate-Re (1 < Re ≤ 40), and high-Re (Re > 40) ranges, where the mixing mechanisms are different. The results show a two-dimensional micromixer of Tesla structure is recommended in the intermediate- and high-Re ranges, while two three-dimensional micromixers with two layers are recommended in the low-Re range due to their excellent mixing performance.

Dielectrophoretic Mixing With Novel Electrode Geometry

2009 ◽  
Author(s):  
G. Naga Siva Kumar ◽  
Sushanta K. Mitra ◽  
Subir Bhattacharjee
Keyword(s):  
Electric Field ◽  
Cell Activation ◽  
Colloidal Particles ◽  
Colloidal Suspension ◽  
Three Dimensional ◽  
Colloidal Suspensions ◽  
Navier Stokes ◽  
Mixing Efficiency ◽  
Computational Domain ◽  
Element Analysis

Electrokinetic mixing of analytes at micro-scale is important in several biochemical applications like cell activation, DNA hybridization, protein folding, immunoassays and enzyme reactions. This paper deals with the modeling and numerical simulation of micromixing of two different types of colloidal suspensions based on principle of dielectrophoresis (DEP). A mathematical model is developed based on Laplace, Navier-Stokes, and convection-diffusion-migration equations to calculate electric field, velocity, and concentration distributions, respectively. Mixing of two colloidal suspensions is simulated in a three-dimensional computational domain using finite element analysis considering dielectrophoretic, gravitational and convective (advective)–diffusive forces. Phase shifted AC signal is applied to the alternating electrodes for achieving the mixing of two different colloidal suspensions. The results indicate that the electric field and DEP forces are maximum at the edges of the electrodes and become minimum elsewhere. As compared to curved edges, straight edges of electrodes have lower electric field and DEP forces. The results also indicate that DEP force decays exponentially along the height of the channel. The effect of DEP forces on the concentration profile is studied. It is observed that, the concentration of colloidal particles at the electrodes edges is very less compared to elsewhere. Mixing of two colloidal suspensions due to diffusion is observed at the interface of the two suspensions. The improvement in mixing after applying the repulsive DEP forces on the colloidal suspension is observed. Most of the mixing takes place across the slant edges of the triangular electrodes. The effect of electrode pairs and the mixing length on degree of mixing efficiency are also observed.

scholarly journals Three-Dimensional Simulation of Fluid–Structure Interaction Problems Using Monolithic Semi-Implicit Algorithm

Fluids ◽  
2019 ◽  
Vol 4 (2) ◽  
pp. 94 ◽  
Author(s):  
Cornel Marius Murea
Keyword(s):  
Stokes Equations ◽  
Fluid Structure Interaction ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Arbitrary Lagrangian Eulerian ◽  
Time Step ◽  
Fluid Structure ◽  
Structure Interaction ◽  
Dimensional Simulation ◽  
Updated Lagrangian

A monolithic semi-implicit method is presented for three-dimensional simulation of fluid–structure interaction problems. The updated Lagrangian framework is used for the structure modeled by linear elasticity equation and, for the fluid governed by the Navier–Stokes equations, we employ the Arbitrary Lagrangian Eulerian method. We use a global mesh for the fluid–structure domain where the fluid–structure interface is an interior boundary. The continuity of velocity at the interface is automatically satisfied by using globally continuous finite element for the velocity in the fluid–structure mesh. The method is fast because we solve only a linear system at each time step. Three-dimensional numerical tests are presented.

Sign in / Sign up

Export Citation Format

Share Document