scholarly journals Turbulent Flow over Steep Steady and Unsteady Waves under Strong Wind Forcing

2018 ◽  
Vol 48 (1) ◽  
pp. 3-27 ◽  
Author(s):  
Peter P. Sullivan ◽  
Michael L. Banner ◽  
Russel P. Morison ◽  
William L. Peirson

AbstractTurbulent flow over strongly forced steep steady and unsteady waves is simulated using large-eddy simulation (LES) with time t and space x varying wave height h(x, t) imposed as a lower boundary condition. With steady waves, h(x, t) is based on measurements of incipient and active breaking waves collected in a wind-wave flume, while a numerical wave code is used to generate an unsteady evolving wave packet (group). Highly intermittent airflow separation is found in the simulations, and the results suggest separation near a wave crest occurs prior to the onset of wave breaking. The form (pressure) drag is most sensitive to the wave slope, and the form drag can contribute as much as 74% to the total stress. Wind and scalar profiles from the LES display log-linear variations above the wave surface; the LES wind profiles are in good agreement with the measurements. The momentum roughness increases as the water surface changes from wind ripples to incipient breaking to active breaking. However, the scalar roughness decreases as the wave surface becomes rougher. This highlights major differences in momentum and scalar transport over a rough wavy surface. For a rapidly evolving, strongly forced wave group, the form drag is highly correlated with the wave slope, and intermittent separation is found early in the packet evolution when the local wave slope −∂h/∂x(x, t) ≥ 0.22. The packet root-mean-square wave slope is 0.084, but the form drag fraction is 2.4 times larger than a comparably forced steady wave. Thus, a passing wave group can induce unsteadiness in the wind stress.

2000 ◽  
Vol 404 ◽  
pp. 47-85 ◽  
Author(s):  
PETER P. SULLIVAN ◽  
JAMES C. McWILLIAMS ◽  
CHIN-HOH MOENG

Turbulent flow over idealized water waves with varying wave slope ak and wave age c/u∗ is investigated using direct numerical simulations at a bulk Reynolds number Re = 8000. In the present idealization, the shape of the water wave and the associated orbital velocities are prescribed and do not evolve dynamically under the action of the wind. The results show that the imposed waves significantly influence the mean flow, vertical momentum fluxes, velocity variances, pressure, and form stress (drag). Compared to a stationary wave, slow (fast) moving waves increase (decrease) the form stress. At small c/u∗, waves act similarly to increasing surface roughness zo resulting in mean vertical velocity profiles with shorter buffer and longer logarithmic regions. With increasing wave age, zo decreases so that the wavy lower surface is nearly as smooth as a flat lower boundary. Vertical profiles of turbulence statistics show that the wave effects depend on wave age and wave slope but are confined to a region kz < 1 (where k is the wavenumber of the surface undulation and z is the vertical coordinate). The turbulent momentum flux can be altered by as much as 40% by the waves. A region of closed streamlines (or cat's-eye pattern) centred about the critical layer height was found to be dynamically important at low to moderate values of c/u∗. The wave-correlated velocity and flux fields are strongly dependent on the variation of the critical layer height and to a lesser extent the surface orbital velocities. Above the critical layer zcr the positions of the maximum and minimum wave-correlated vertical velocity ww occur upwind and downwind of the peak in zcr, like a stationary surface. The wave-correlated flux uwww is positive (negative) above (below) the critical layer height.


2014 ◽  
Vol 44 (2) ◽  
pp. 742-763 ◽  
Author(s):  
Yevgenii Rastigejev ◽  
Sergey A. Suslov

Abstract In-depth understanding and accurate modeling of the interaction between ocean spray and a turbulent flow under high wind conditions is essential for improving the intensity forecasts of hurricanes and severe storms. Here, the authors consider the E–ε closure for a turbulent flow model that accounts for the effects of the variation of turbulent energy and turbulent mixing length caused by spray stratification. The obtained analytical and numerical solutions show significant differences between the current E–ε model and the lower-order turbulent kinetic energy (TKE) model considered previously. It is shown that the reduction of turbulent energy and mixing length above the wave crest level, where the spray droplets are generated, that is not accounted for by the TKE model results in a significant suppression of turbulent mixing in this near-wave layer. In turn, suppression of turbulence causes an acceleration of flow and a reduction of the drag coefficient that is qualitatively consistent with field observations if spray is fine (even if its concentration is low) or if droplets are large but their concentration is sufficiently high. In the latter case, spray inertia may become important. This effect is subsequently examined. It is shown that spray inertia leads to the reduction of wind velocity in the close proximity of the wave surface relative to the reference logarithmic profile. However, at higher altitudes the suppression of flow turbulence by the spray still results in the wind acceleration and the reduction of the local drag coefficient.


2004 ◽  
Vol 5 (6) ◽  
pp. 1131-1146 ◽  
Author(s):  
H. Richter ◽  
A. W. Western ◽  
F. H. S. Chiew

Abstract Numerical Weather Prediction (NWP) and climate models are sensitive to evapotranspiration at the land surface. This sensitivity requires the prediction of realistic surface moisture and heat fluxes by land surface models that provide the lower boundary condition for the atmospheric models. This paper compares simulations of a stand-alone version of the European Centre for Medium-Range Weather Forecasts (ECMWF) land surface scheme, or the Viterbo and Beljaars scheme (VB95), with various soil and vegetation parameter sets against soil moisture observations across the Murrumbidgee River catchment in southeast Australia. The study is, in part, motivated by the adoption of VB95 as the operational land surface scheme by the Australian Bureau of Meteorology in 1999. VB95 can model the temporal fluctuations in soil moisture, and therefore the moisture fluxes, fairly realistically. The monthly model latent heat flux is also fairly insensitive to soil or vegetation parameters. The VB95 soil moisture is sensitive to the soil and, to a lesser degree, the vegetation parameters. The model exhibits a significant (generally wet) bias in the absolute soil moisture that varies spatially. The use of the best Australia-wide available soils and vegetation information did not improve VB95 simulations consistently, compared with the original model parameters. Comparisons of model and observed soil moistures revealed that more realistic soil parameters are needed to reduce the model soil moisture bias. Given currently available continent-wide soils parameters, any initialization of soil moisture with observed values would likely result in significant flux errors. The soil moisture bias could be largely eliminated by using soil parameters that were derived directly from the actual soil moisture observations. Such parameters, however, are only available at very few point locations.


1976 ◽  
Vol 1 (15) ◽  
pp. 25 ◽  
Author(s):  
Edward B. Thornton ◽  
James J. Galvin ◽  
Frank L. Bub ◽  
David P. Richardson

The sight and sound of breaking waves and surf is so familiar and enjoyable that we tend to forget how little we really understand about them. Why is it, that compared to other branches of wave studies our knowledge of breaking waves is so empirical and inexact? The reason must lie partly in the difficulty of finding a precise mathematical description of a fluid flow that is in general nonlinear and time-dependent. The fluid accelerations can no longer be assumed t o be small compared t o gravity, as in Stokes's theory for periodic waves and the theory of cnoidal waves in shallow water, nor is the particle velocity any longer small compared to the phase velocity. The aim of this paper is to bring together s ome recent contributions to the calculation both of steep symmetric waves and of time-dependent surface waves. These have a bearing on the behaviour of whitecaps in deep water and of surf in the breaker zone . Since spilling breakers in gently shoaling water closely resemble solitary waves, we begin with the description of solitary waves of limiting amplitude, then discuss steep waves of arbitrary height. The observed intermittency of whitecaps is discussed in terms of the energy maximum, as a function of wave steepness, In Sections 6 and 7 a simpler description of steady symmetric waves is proposed, using an asymptotic expression for the flow near the wave crest. Finally we describe a new numerical technique (MEL, or mixed Eulerian-Lagrangian) with which it has been found possible to follow the development of periodic waves past the point when overturning takes place. Measurement of waves, and vertical and horizontal water particle velocities were made of spilling, plunging and surging breakers at sandy beaches in the vicinity of Monterey, California. The measured breaking waves, derived characteristically from swell-type waves, can be described as highly nonlinear. Spectra and cross spectra were calculated between waves and velocities. Secondary waves were noted visually and by the strong harmonics in the spectra. The strength of the harmonics is related to the beach steepness, wave height and period. The phase difference between waves and horizontal velocities indicates the unstable crest of the wave leads the velocities on the average by 5-20 degrees. Phase measurements between wave gauges in a line perpendicular to the shore show breaking waves to be frequency nondispersive indicating phase-coupling of the various wave components. The coherence squared values between the sea surface elevation and the horizontal water particle velocity were high in all runs, ranging above 0.8 at the peak of the spectra. The high coherence suggests that most of the motion in the body of breaking waves is wave-induced and not turbulent.


Author(s):  
François Lott ◽  
Bruno Deremble ◽  
Clément Soufflet

AbstractThe non-hydrostatic version of the mountain flow theory presented in Part I is detailed. In the near neutral case, the surface pressure decreases when the flow crosses the mountain to balance an increase in surface friction along the ground. This produces a form drag which can be predicted qualitatively. When stratification increases, internal waves start to control the dynamics and the drag is due to upward propagating mountain waves as in part I. The reflected waves nevertheless add complexity to the transition. First, when stability increases, upward propagating waves and reflected waves interact destructively and low drag states occur. When stability increases further, the interaction becomes constructive and high drag state are reached. In very stable cases the reflected waves do not affect the drag much. Although the drag gives a reasonable estimate of the Reynolds stress, its sign and vertical profile are profoundly affected by stability. In the near neutral case the Reynolds stress in the flow is positive, with maximum around the top of the inner layer, decelerating the large-scale flow in the inner layer and accelerating it above. In the more stable cases, on the contrary, the large-scale flow above the inner layer is decelerated as expected for dissipated mountain waves. The structure of the flow around the mountain is also strongly affected by stability: it is characterized by non separated sheltering in the near neutral cases, by upstream blocking in the very stable case, and at intermediate stability by the presence of a strong but isolated wave crest immediately downstream of the ridge.


Author(s):  
Simon SW Li ◽  
Daniel HK Chow

This study modified an electromyography-assisted optimization approach for predicting lumbar spine loading while walking with backpack loads. The modified-electromyography-assisted optimization approach eliminated the electromyography measurement at maximal voluntary contraction and adopted a linear electromyography–force relationship. Moreover, an optimal lower boundary condition for muscle gain was introduced to constrain the trunk muscle co-activation. Anthropometric information of 10 healthy young men as well as their kinematic, kinetic, and electromyography data obtained while walking with backpack loads were used as inputs in this study. A computational algorithm was used to find and analyse the sensitivity of the optimal lower boundary condition for achieving minimum deviation of the modified-electromyography-assisted optimization approach from the electromyography-assisted optimization approach for predicting lumbosacral joint compression force. Results validated that the modified-electromyography-assisted optimization approach (at optimal lower boundary condition of 0.92) predicted on average, a non-significant deviation in peak lumbosacral joint compression force of −18 N, a standard error of 9 N, and a root mean square difference in force profile of 73.8 N. The modified-electromyography-assisted optimization approach simplified the experimental process by eliminating the electromyography measurement at maximal voluntary contraction and provided comparable estimations for lumbosacral joint compression force that is also applicable to patients or individuals having difficulty in performing the maximal voluntary contraction activity.


2018 ◽  
Vol 2018 ◽  
pp. 1-17
Author(s):  
Jun Liu ◽  
Hongliang Zhao ◽  
Simon X. Yang ◽  
Qingyou Liu ◽  
Guorong Wang

The landing string is an important component of deepwater riserless drilling systems. Determination of the dynamic characteristics of the landing string plays an essential role in its design for ensuring its safe operation. In this paper, a dynamic model is developed to investigate the dynamic response characteristics of a landing string, where a landing string in a marine environment is modeled as a flexible slender tube undergoing coupled transverse and axial motions. The heaving motion of the drilling platform is taken as the upper boundary condition and the motion of the drilling bit caused by the interaction between the rock and the bit as the lower boundary condition. A semiempirical Morison equation is used to simulate the effect of the load imposed by the marine environment. The dynamic model, which is nonlinearly coupled and multibody, is discretized by a finite element method and solved by the Newmark technique. Using the proposed model, the dynamic responses of the displacement, axial force, and moment in the landing string are investigated in detail to find out the influences of driving depth of surface catheter, platform motion, bit movement, and marine environment on the dynamical characteristics of the landing string.


Author(s):  
Yanfei Deng ◽  
Jianmin Yang ◽  
Longfei Xiao

In the last few decades, the hydrodynamic performance of offshore structures has been widely studied to ensure their safety as well as to achieve an economical design. However, an increasing number of reported accidents due to rough ocean waves call for in-depth investigations on the loads and motions of offshore structures, particularly the effect of freak waves. The aim of this paper is to determine the sea conditions that may cause the maximum motion responses of offshore structures, which have a significant effect on the loads of mooring systems because of their tight relationship. As a preliminary step, the response amplitude operators (RAOs) of a semisubmersible platform of 500 meters operating depth are obtained with the frequency-domain analysis method. Subsequently, a series of predetermined extreme wave sequences with different wave group characteristics, such as the maximum crest amplitude and the time lag between successive high waves, are adopted to calculate the hydrodynamic performance of the semisubmersible with mooring systems in time-domain. The paper shows that the maximum motion responses not only depend on the largest wave crest amplitude but also the time lags between successive giant waves. This paper will provide an important reference for future designs which could consider the most dangerous wave environment.


2001 ◽  
Vol 25 (4) ◽  
pp. 483-511 ◽  
Author(s):  
Gareth Roberts

This paper presents a review of the application of Bi-directional Reflectance Distribution Function (BRDF) models in the inference of land surface parameters at regional and global scales using remotely sensed data. Information on land surface parameters, such as Leaf Area Index (LAI), fraction of Absorbed Photosynthetically Active Radiation (fAPAR), aerodynamic surface roughness and albedo, are valuable for understanding the transfer of energy and mass between terrestrial ecosystems and the atmosphere (e.g., carbon, nitrogen and methane cycling) and for ingestion into the lower boundary condition of global circulation models (GCM)s. Conventional techniques for acquiring information on land surface parameters do not account for or utilize the directional nature of surface reflectance. This paper reviews empirical, semi-empirical and, to a lesser extent, physical BRDF models that describe the surface BRDF. In each case examples are given of their application in inferring land surface parameters. The review concludes by discussing the future prospects of BRDF modelling using spaceborne sensors.


2016 ◽  
Vol 29 (15) ◽  
pp. 5643-5659 ◽  
Author(s):  
Maria A. A. Rugenstein ◽  
Jonathan M. Gregory ◽  
Nathalie Schaller ◽  
Jan Sedláček ◽  
Reto Knutti

Abstract In radiative forcing and climate feedback frameworks, the initial stratospheric and tropospheric adjustments to a forcing agent can be treated as part of the forcing and not as a feedback, as long as the average global surface temperature response is negligible. Here, a very large initial condition ensemble of the Community Earth System Model is used to analyze how the ocean shapes the fast response to radiative forcing. It is shown that not only the stratosphere and troposphere but also the ocean adjusts. This oceanic adjustment includes meridional ocean heat transport convergence anomalies, which are locally as large as the surface heat flux anomalies, and an increase of the Atlantic meridional overturning circulation. These oceanic adjustments set the lower boundary condition for the atmospheric response of the first few years, in particular, the shortwave cloud radiative effect. This cloud adjustment causes a nonlinear relationship between global energy imbalance and temperature. It proceeds with a characteristic time scale of a few years in response to the forcing rather than scaling nonlinearly with global mean temperature anomaly. It is proposed that even very short time scales are treated as a fully coupled problem and encourage other modeling groups to investigate whether our description also suits their models’ behavior. A definition of the forcing term (“virtual forcing”) including oceanic adjustment processes is introduced and serves as an interpretive idea for longer time scales.


Sign in / Sign up

Export Citation Format

Share Document