blips, coils, and pedals: instabilities of buoyant jets in fluids with stable density stratification...
TRANSCRIPT
![Page 1: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/1.jpg)
Blips, Coils, and Pedals: Instabilities of buoyant jets in
fluids with stable density stratification
Sam Schofield
Juan Restrepo
Raymond Goldstein
Program in Applied MathematicsThe University of Arizona
![Page 2: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/2.jpg)
Abstract
In biological flow experiments, algae gather in a negatively buoyant plume that breaks apart into round blips as it descends. The question of whether the motion of the individual cells is essential to the formation of the instability was posed. Recent three-dimensional and Hele-Shaw experiments have demonstrated that similar behaviors are found in dense jets descending into stably stratified fluids. These experiments have shown a rich array of dynamics including buckling and overturning, blip formation, shear type instabilities in the entrained conduit and multiple jet interactions. However, the flows are quite complex requiring viscous, buoyancy, convection and stratification effects to be considered together. Direct numerical simulation of the Boussinesq equations with high order spectral element methods allows the detailed aspects of the flows to be captured and understood without resorting to excessive approximations. We show that these behaviors are not tied to cell motions in biological flows, surface tensions effects, or Hele-Shaw cells and can be understood in terms of the parameter regimes being considered. Furthermore, we find a vortex street behavior that may be the pedal breakdown observed in other liquid into liquid jet experiments.
![Page 3: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/3.jpg)
Overview
We consider the behavior of a dense jet injected into a fluid with a stable density stratification
Example application areas: • Atmospheric circulation• Ocean models• Effluent discharges• Smokestacks• Bacterial Swarms
![Page 4: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/4.jpg)
Historical Developments
• A.J. Reynolds, JFM Feb 15, 1962 “Observations of a liquid into liquid jet”– Observed (a) shearing puffs, (b) symmetric condensations, (c) sinuous
undulations, (d) pedal breakdown and (e) confused breakdown depending on the injection rate.
– Breakdown was observed a Reynolds number around 300 d
QVdRe
4
![Page 5: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/5.jpg)
Hele-Shaw Experiments
by A.L. Pesci, et al. PRL 2003
• A dense jet is injected into a tank with a stable density stratification.
• The jet is injected at the top from a needle with .24mm diameter at velocities of 0.02cm/s to 1 cm/s. (Re = 0.03 to 3)
• Visualized using the difference in index of refraction from concentration variations
Salinity
y
30 cm
30 cm
0.24mm diameter nozzle
Salinity gradient 0.04M/l-cm
3 mm
![Page 6: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/6.jpg)
Hele-Shaw Experiments - Buckling Jet.mov.lnk
![Page 7: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/7.jpg)
Hele-Shaw Experiments – Conduit
![Page 8: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/8.jpg)
Mathematical Formulation
•The small but nonzero Reynolds number does not permit the use of boundary layer approximations. T•There is not a similarity structure to the governing equations making analysis difficult.•Viscous, convection, and buoyancy effects must all be considered together.
LU
Re 0
2
30
LSSg
Gr j
DSc
Reynolds number
Grasshof number
Schmidt number
![Page 9: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/9.jpg)
Mathematical Formulation and Direct Numerical Simulation (DNS)
The Boussinesq approximation to the incompressible Navier-Stokes equations was used:
where u is velocity, S is salinity, κ is diffusivity, ν is kinematic viscosity, ρ is the density.
DNS was done with a high-order spectral element method. 4740 elements were used with ~330,000 grid points.
The relevant parameters are the injection rate, the density of injected jet, the density stratification, and the viscosity.
SDSt
S
SSgpt
u
u
uuuu
0
1
![Page 10: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/10.jpg)
Direct Numerical Simulation
![Page 11: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/11.jpg)
Spectral Element Method
Example (Maday and Patera, 1989)
1D Helmholtz problem
11
0
011
0
2
,x
)x(p
)(u)(u
)x(fu)u)x(p(
Ref:
Maday, Y. and Patera, A.T. “Spectral element methods for the Navier-Stokes equations,” in State of the
Art Surveys in Computational Mechanics, A.K. Noor, ed., ASME, New York, pp. 71-143 (1989)
Deville, M.O, Fischer, P.F, and Mund E.H. High order methods for incompressible fluid flow. Cambridge
University Press, 2002
![Page 12: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/12.jpg)
Look for a solution in the Sobolev Space
Recast the equation as a variational problem
The domain can be mapped back to a reference domain (i.e. [-1,1]). This is usually done at the discrete level using an affine transformation or using the Gordon-Hall algorithm (for more complex deformation in higher dimensions)
0)(),('),( 2210 vLvLvH
dx)x(v)x(fdx)x(v)x(udx)x(v)x(u)x(p),(Hv 210
dxxvxudxxvxuxpvua
vfvuaHv
)()()()()(),(
),(),(),(2
10
Weak formulationWeak formulationDiscretizationQuadrature RulesBasis FunctionsTensor Products
![Page 13: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/13.jpg)
Discretization
Split Ω into subintervals (spectral elements)
Now look for a solution in this subspace Xh
k1
on degree of sPolynomial
N)(P
)K,N(h
kN
K
kk
dx)x(v)x(fdx)x(v)x(udx)x(v)x(u)x(p,Xv hhhhhh2
Weak formulationDiscretizationQuadrature RulesBasis FunctionsTensor Products
![Page 14: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/14.jpg)
Gauss-Lobatto Quadrature
)(v)(udr)r(v)r(u]),,([Pv,u i
N
iiiN
0
1
1
11
Approximate the inner product (u,v) with a finite sum.
Select the quadrature points as the Gauss-Lobatto or Gauss-Lobatto-Legendre points to integrate exactly polynomials of degree 2N-1.
011110
012
)(L;N,,i,,
)(dr)r(),(P
iNN
N
iiiN
Weak formulationDiscretizationQuadrature RulesBasis FunctionsTensor Products
![Page 15: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/15.jpg)
Basis for Solution Space
Choose a basis for Xh, such as Lagrangian Interpolants where
Express functions in Xh in terms of basis functions
Enforce continuity
Boundary conditions
ijiih )(
1,,1,10 NkWW kk
N
010 K
NWW
N
ii
ki
k rhwrWh
0
)()(
Note the choice of Gauss-Lobatto Lagrangian interpolants result in coupling only at the endpoints and a diagonal mass matrix.
)(h)(hdx)x(h)x(h in
N
iiminm
0
1
1
Weak formulationDiscretizationQuadrature RulesBasis FunctionsTensor Products
![Page 16: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/16.jpg)
Tensor Product MatricesWeak formulationDiscretizationQuadrature RulesBasis FunctionsTensor Products
For higher dimensions, use tensor product of GLL points
Deville, M.O, Fischer, P.F, and Mund E.H. High order methods for incompressible fluid flow. Cambridge
University Press, 2002
![Page 17: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/17.jpg)
Tensor Product MatricesWeak formulationDiscretizationQuadrature RulesBasis FunctionsTensor Products
Operators also take a tensor product form,
x
x
x
x
x
D
D
D
D
DID
ˆ
ˆ
ˆ
ˆ
ˆx
Evaluation of the action of operator is efficient 1dnO
![Page 18: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/18.jpg)
Jet
![Page 19: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/19.jpg)
Plume
![Page 20: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/20.jpg)
Flow zones for buoyant jets in stratified fluids(Tenner and Gebhart ’71)
1. Jet core
2. Jet-cell interface
3. Toroidal cell
4. Shroud
5. Far Field
6. Cloud
7. Shroud Base
A. R. Tenner and Gebhart, B. “Laminar and axysmmetric jets in a stably stratified environment” Int. J. Heat and Mass Transfer, 14, 1971
![Page 21: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/21.jpg)
![Page 22: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/22.jpg)
![Page 23: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/23.jpg)
![Page 24: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/24.jpg)
Stability
xS
xSSL
LULxSSg js
sj s
)0(Re
)0(Gr 0
2
3
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
-8 -7 -6 -5 -4 -3 -2 -1 0
-log(Gr/Re^2)
Li/Ls
Blip
Blip/Sine
Sine
![Page 25: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/25.jpg)
![Page 26: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/26.jpg)
![Page 27: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/27.jpg)
Conduit Instability
![Page 28: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/28.jpg)
Blips1
0
0
001.01
/04.0
100
cmx
scmU
SC
14.4s 14.6s 14.8s 15.0s 15.2s
Salinity.mov
![Page 29: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/29.jpg)
Transitory Blips/Sines
![Page 30: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/30.jpg)
Buckling Instability1
0
0
001.01
/50.0
100
cmx
scmU
SC
5.6s 5.8s 6.0s 6.2s 6.4s
Salinity.mov
![Page 31: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/31.jpg)
Vortex Street (Sc=7)
![Page 32: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/32.jpg)
Planar Coiling (Sc=7)
3.80s 3.85s 3.90s 3.95s 4.0s3.75s3.70s
![Page 33: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/33.jpg)
Planar Coiling (Sc=100)
4.3s 4.4s 4.5s 4.6s4.2s 4.7s
![Page 34: Blips, Coils, and Pedals: Instabilities of buoyant jets in fluids with stable density stratification Sam Schofield Juan Restrepo Raymond Goldstein Program](https://reader035.vdocument.in/reader035/viewer/2022062308/56649f565503460f94c7a6e9/html5/thumbnails/34.jpg)
Acknowledgements
• Krell Institute, DOE Computational Science Graduate Fellowship• Dr. Paul Fischer, Argonne National Laboratory
Ongoing Research
• Numerical bifurcation analysis• The stability impact of the evolution of velocity profile along the jet axis. • Analytical/reduced models• Interaction of the jet buckling instability with conduit.• Multiple jets interactions.