Flow past a sphere moving vertically in a stratified diffusive fluid

2000 ◽  
Vol 417 ◽  
pp. 211-236 ◽  
Author(s):  
C. R. TORRES ◽  
H. HANAZAKI ◽  
J. OCHOA ◽  
J. CASTILLO ◽  
M. VAN WOERT

Numerical studies are described of the flows generated by a sphere moving vertically in a uniformly stratified fluid. It is found that the axisymmetric standing vortex usually found in homogeneous fluids at moderate Reynolds numbers (25 [les ] Re [les ] 200) is completely collapsed by stable stratification, generating a strong vertical jet. This is consistent with our experimental visualizations. For Re = 200 the complete collapse of the vortex occurs at Froude number F ≃ 19, and the critical Froude number decreases slowly as Re increases. The Froude number and the Reynolds number are here defined by F = W/Na and Re = 2Wa/v, with W being the descent velocity of the sphere, N the Brunt–Väisälä frequency, a the radius of the sphere and v the kinematic viscosity coefficient. The inviscid processes, including the generation of the vertical jet, have been investigated by Eames & Hunt (1997) in the context of weak stratification without buoyancy effects. They showed the existence of a singularity of vorticity and density gradient on the rear axis of the flow and also the impossibility of realizing a steady state. When there is no density diffusion, all the isopycnal surfaces which existed initially in front of the sphere accumulate very near the front surface because of density conservation and the fluid in those thin layers generates a rear jet when returning to its original position. In the present study, however, the fluid has diffusivity and the buoyancy effects also exist. The density diffusion prevents the extreme piling up of the isopycnal surfaces and allows the existence of a steady solution, preventing the generation of a singularity or a jet. On the other hand, the buoyancy effect works to increase the vertical velocity to the rear of the sphere by converting the potential energy to vertical kinetic energy, leading to the formation of a strong jet. We found that the collapse of the vortex and the generation of the jet occurs at much weaker stratifications than those necessary for the generation of strong lee waves, showing that jet formation is independent of the internal waves. At low Froude numbers (F [les ] 2) the lee wave patterns showed good agreement with the linear wave theory and the previous experiments by Mowbray & Rarity (1967). At very low Froude numbers (F [les ] 1) the drag on a sphere increases rapidly, partly due to the lee wave drag but mainly due to the large velocity of the jet. The jet causes a reduction of the pressure on the rear surface of the sphere, which leads to the increase of pressure drag. High velocity is induced also just outside the boundary layer of the sphere so that the frictional drag increases even more significantly than the pressure drag.

2010 ◽  
Vol 67 (7) ◽  
pp. 2260-2278 ◽  
Author(s):  
Stephen D. Eckermann ◽  
John Lindeman ◽  
Dave Broutman ◽  
Jun Ma ◽  
Zafer Boybeyi

Abstract Fully nonlinear mesoscale model simulations are used to investigate the momentum fluxes of gravity waves that emerge at a “far-field” height of 6 km from steady unsheared flow over both an axisymmetric and elliptical obstacle for nondimensional mountain heights ĥm = Fr−1 in the range 0.1–5, where Fr is the surface Froude number. Fourier- and Hilbert-transform diagnostics of model output yield local estimates of phase-averaged momentum flux, while area integrals of momentum flux quantify the amount of surface pressure drag that translates into far-field gravity waves, referred to here as the “wave drag” component. Estimates of surface and wave drag are compared to parameterization predictions and theory. Surface dynamics transition from linear to high-drag (wave breaking) states at critical inverse Froude numbers Frc−1 predicted to within 10% by transform relations. Wave drag peaks at Frc−1 < ĥm ≲ 2, where for the elliptical obstacle both surface and wave drag vacillate owing to cyclical buildup and breakdown of waves. For the axisymmetric obstacle, this occurs only at ĥm = 1.2. At ĥm ≳ 2–3 vacillation abates and normalized pressure drag assumes a common normalized form for both obstacles that varies approximately as ĥm−1.3. Wave drag in this range asymptotes to a constant absolute value that, despite its theoretical shortcomings, is predicted to within 10%–40% by an analytical relation based on linear clipped-obstacle drag for a Sheppard-based prediction of dividing streamline height. Constant wave drag at ĥm ∼ 2–5 arises despite large variations with ĥm in the three-dimensional morphology of the local wave momentum fluxes. Specific implications of these results for the parameterization of subgrid-scale orographic drag in global climate and weather models are discussed.


Author(s):  
Callum J. Shakespeare ◽  
Brian K. Arbic ◽  
Andrew McC. Hogg

AbstractInternal waves generated at the seafloor propagate through the interior of the ocean, driving mixing where they break and dissipate. However, existing theories only describe these waves in two limiting cases. In one limit, the presence of an upper boundary permits bottom-generated waves to reflect from the ocean surface back to the seafloor, and all the energy flux is at discrete wavenumbers corresponding to resonant modes. In the other limit, waves are strongly dissipated such that they do not interact with the upper boundary and the energy flux is continuous over wavenumber. Here, a novel linear theory is developed for internal tides and lee waves that spans the parameter space in between these two limits. The linear theory is compared with a set of numerical simulations of internal tide and lee wave generation at realistic abyssal hill topography. The linear theory is able to replicate the spatially-averaged kinetic energy and dissipation of even highly non-linear wave fields in the numerical simulations via an appropriate choice of the linear dissipation operator, which represents turbulent wave breaking processes.


2012 ◽  
Vol 708 ◽  
pp. 250-278 ◽  
Author(s):  
Nicolas Grisouard ◽  
Oliver Bühler

AbstractWe present a theoretical and numerical study of the effective mean force exerted on an oceanic mean flow due to the presence of small-amplitude internal waves that are forced by the oscillatory flow of a barotropic tide over undulating topography and are also subject to dissipation. This extends the classic lee-wave drag problem of atmospheric wave–mean interaction theory to a more complicated oceanographic setting, because now the steady lee waves are replaced by oscillatory internal tides and, most importantly, because now the three-dimensional oceanic mean flow is defined by time averaging over the fast tidal cycles rather than by the zonal averaging familiar from atmospheric theory. Although the details of our computation are quite different, we recover the main action-at-a-distance result from the atmospheric setting, namely that the effective mean force that is felt by the mean flow is located in regions of wave dissipation, and not necessarily near the topographic wave source. Specifically, we derive an explicit expression for the effective mean force at leading order using a perturbation series in small wave amplitude within the framework of generalized Lagrangian-mean theory, discuss in detail the range of situations in which a strong, secularly growing mean-flow response can be expected, and then compute the effective mean force numerically in a number of idealized examples with simple topographies.


2014 ◽  
Vol 44 (9) ◽  
pp. 2409-2432 ◽  
Author(s):  
H. W. Wijesekera ◽  
E. Jarosz ◽  
W. J. Teague ◽  
D. W. Wang ◽  
D. B. Fribance ◽  
...  

Abstract Pressure differences across topography generate a form drag that opposes the flow in the water column, and viscous and pressure forces acting on roughness elements of the topographic surface generate a frictional drag on the bottom. Form drag and bottom roughness lengths were estimated over the East Flower Garden Bank (EFGB) in the Gulf of Mexico by combining an array of bottom pressure measurements and profiles of velocity and turbulent kinetic dissipation rates. The EFGB is a coral bank about 6 km wide and 10 km long located at the shelf edge that rises from 100-m water depth to about 18 m below the sea surface. The average frictional drag coefficient over the entire bank was estimated as 0.006 using roughness lengths that ranged from 0.001 cm for relatively smooth portions of the bank to 1–10 cm for very rough portions over the corals. The measured form drag over the bank showed multiple time-scale variability. Diurnal tides and low-frequency motions with periods ranging from 4 to 17 days generated form drags of about 2000 N m−1 with average drag coefficients ranging between 0.03 and 0.22, which are a factor of 5–35 times larger than the average frictional drag coefficient. Both linear wave and quadratic drag laws have similarities with the observed form drag. The form drag is an important flow retardation mechanism even in the presence of the large frictional drag associated with coral reefs and requires parameterization.


2016 ◽  
Vol 97 ◽  
pp. 109-128 ◽  
Author(s):  
David S. Trossman ◽  
Brian K. Arbic ◽  
James G. Richman ◽  
Stephen T. Garner ◽  
Steven R. Jayne ◽  
...  
Keyword(s):  

2013 ◽  
Vol 72 ◽  
pp. 119-142 ◽  
Author(s):  
David S. Trossman ◽  
Brian K. Arbic ◽  
Stephen T. Garner ◽  
John A. Goff ◽  
Steven R. Jayne ◽  
...  

2008 ◽  
Vol 24 (1) ◽  
pp. 35-42 ◽  
Author(s):  
Per-Ludvik Kjendlie ◽  
Robert Keig Stallman

The aims of this study were to compare drag in swimming children and adults, quantify technique using the technique drag index (TDI), and use the Froude number (Fr) to study whether children or adults reach hull speed at maximal velocity (vmax). Active and passive drag was measured by the perturbation method and a velocity decay method, respectively, including 9 children aged 11.7 ± 0.8 and 13 adults aged 21.4 ± 3.7. The children had significantly lower active (kAD) and passive drag factor (kPD) compared with the adults. TDI (kAD/kPD) could not detect any differences in swimming technique between the two groups, owing to the adults swimming maximally at a higher Fr, increasing the wave drag component, and masking the effect of better technique. The children were found not to reach hull speed atvmax, and their Fr were 0.37 ± 0.01 vs. the adults 0.42 ± 0.01, indicating adults’ larger wave-making component of resistance atvmaxcompared with children. Fr is proposed as an evaluation tool for competitive swimmers.


1968 ◽  
Vol 32 (3) ◽  
pp. 549-567 ◽  
Author(s):  
John W. Miles

The lee-wave amplitudes and wave drag for a thin barrier in a two-dimensional stratified flow in which the upstream dynamic pressure and density gradient are constant (Long's model) are determined as functions of barrier height and Froude number for a channel of finite height and for a half-space. Variational approximations to these quantities are obtained and validated by comparison with the earlier results of Drazin & Moore (1967) for the channel and with the results of an exact solution for the half-space, as obtained through separation of variables. An approximate solution of the integral equation for the channel also is obtained through a reduction to a singular integral equation of potential theory. The wave drag tends to increase with decreasing wind speed, but it seems likely that the flow is unstable in the region of high drag. The maximum attainable drag coefficient consistent with stable lee-wave formation appears to be roughly two and almost certainly less than three.


1969 ◽  
Vol 35 (3) ◽  
pp. 497-525 ◽  
Author(s):  
John W. Miles ◽  
Herbert E. Huppert

A two-dimensional stratified flow over an obstacle in a half space is considered on the assumptions that the upstream dynamic pressure and density gradient are constant (Long's model). A general solution of the resulting boundary-value problem is established in terms of an assumed distribution of dipole sources. Asymptotic solutions for prescribed bodies are established for limiting values of the slenderness ratio ε (height/breadth) of the obstacle and the reduced frequency k (inverse Froude number based on the obstacle breadth) as follows: (i) ε → 0 withkfixed; (ii)k→ 0 with ε fixed; (iii)k→ ∞ withkεfixed. The approximation (i) is deveoped to both first (linearized theory) and second order in ε in terms of Fourier integrals. The approximation (ii), which constitutes a modification of Rayleigh-scattering theory, is obtained by the method of matched asymptotic expansions and depends essentially on thedipole form(which is proportional to the sum of the displaced and virtual masses) of the obstacle with respect to a uniform flow. A simple approximation to this dipole form is proposed and validated by a series of examples in an appendix. The approximation (iii) is obtained through the reduction of the original integral equation to a singular integral equation of Hilbert's type that is solved by the techniques of function theory. A composite approximation to the lee-wave field that is valid in each of the limits (i)-(iii) also is obtained. The approximation (iii) yields an estimate of the maximum value ofkεfor which completely stable lee-wave formation behind a slender obstacle is possible. The differential and total scattering cross-sections and the wave drag on the obstacle are related to the power spectrum of the dipole density. It is shown that the drag is invariant under a reversal of the flow in the limits (i) and (ii), but only for a symmetric obstacle in the limit (iii). The results are applied to a semi-ellipse, an asymmetric generalization thereof, the Witch of Agnesi (Queney's mountain), and a rectangle. The approximate results for the semi-ellipse are compared with the more accurate results obtain by Huppert & Miles (1969). It appears from this comparison that the approximate solutions should be adequate for any slender obstacle within the range ofkεfor which completely stable lee-wave formation is possible. The extension to obstacles in a channel of finite height is considered in an appendix.


2013 ◽  
Vol 70 (9) ◽  
pp. 2930-2947 ◽  
Author(s):  
Miguel A. C. Teixeira ◽  
José Luis Argaín ◽  
Pedro M. A. Miranda

Abstract The drag produced by 2D orographic gravity waves trapped at a temperature inversion and waves propagating in the stably stratified layer existing above are explicitly calculated using linear theory, for a two-layer atmosphere with neutral static stability near the surface, mimicking a well-mixed boundary layer. For realistic values of the flow parameters, trapped-lee-wave drag, which is given by a closed analytical expression, is comparable to propagating-wave drag, especially in moderately to strongly nonhydrostatic conditions. In resonant flow, both drag components substantially exceed the single-layer hydrostatic drag estimate used in most parameterization schemes. Both drag components are optimally amplified for a relatively low-level inversion and Froude numbers Fr ≈ 1. While propagating-wave drag is maximized for approximately hydrostatic flow, trapped-lee-wave drag is maximized for l2a = O(1) (where l2 is the Scorer parameter in the stable layer and a is the mountain width). This roughly happens when the horizontal scale of trapped lee waves matches that of the mountain slope. The drag behavior as a function of Fr for l2H = 0.5 (where H is the inversion height) and different values of l2a shows good agreement with numerical simulations. Regions of parameter space with high trapped-lee-wave drag correlate reasonably well with those where lee-wave rotors were found to occur in previous nonlinear numerical simulations including frictional effects. This suggests that trapped-lee-wave drag, besides giving a relevant contribution to low-level drag exerted on the atmosphere, may also be useful to diagnose lee-rotor formation.


Sign in / Sign up

Export Citation Format

Share Document