As part of the ONR sponsored Coupled Boundary Layer/Air-Sea Transfer (CBLAST) experiment, data from the NOAA WP-3D research aircraft measurements into major Hurricanes in the 2002-2004 seasons are analyzed to investigate the structure of the boundary layer. The turbulent fluxes of momentum and enthalpy are derived using the eddy correlation method. For the first time, the momentum and enthalpy fluxes were directly measured in the boundary layer of a hurricane with wind speeds up to 30 m/s. A new bulk parameterization of the momentum and enthalpy flux is developed. The vertical structure of turbulence and fluxes are presented for the entire boundary layer in the rain free region between the outer rainbands. The turbulent kinetic energy budget was estimated for the hurricane boundary layer between the outer rainbands. The universal spectra and cospectra of the wind velocity, temperature and humidity are also derived. A case study on the effects of roll vortices on the turbulent fluxes is conducted, which confirmed the existence of the boundary layer rolls and gave the first estimate of their modulation of the momentum and sensible heat flux. The CBLAST data provided an invaluable perspective on the evaluation and development of the boundary layer parameterization suited for the hurricane models. Studies on entrainment processes above of the mixed layer and turbulent transport processes induced by the inflow are recommended in the future.

Yes / Researchers have recognized that the successive growth of mid-channel bar deposits can be entertained as the raison d’être for the initiation of the braiding process, which is closely interlinked with the growth, decay, and vertical distribution of fluvial turbulent kinetic energy (TKE). Thus, focused analysis on the underlying mechanics of turbulent flow structures in the proximity of a bar deposit occurring in the middle of the channel can afford crucial scientific clues for insight into the initiating fluvial processes that give rise to braiding. In the study reported herein, a physical model of a mid-channel bar is constructed in an experimental flume to analyze the turbulence parameters in a region close to the bar. Notably, the flow velocity plays an important role in understanding the flow behavior in the scour-hole location in the upstream flow divergence zone as well as near the downstream zone of flow convergence in a mid-channel bar. Therefore, the fluctuating components of turbulent flow velocity are herein discussed and analyzed for the regions located close to the bar. In the present study, the impact of the mid-channel bar, as well as its growth in turbulent flow, on higher-order velocity fluctuation moments are investigated. For near-bed locations, the results show the dominance of ejection events in upstream zones and the dominance of sweep events at locations downstream of the mid-channel bar. In scour-hole sections, the negative value of the stream-wise flux of turbulent kinetic energy and the positive value of the vertical flux of turbulent kinetic energy indicate energy transport in downward and forward directions, respectively. The downward and forward energy transport processes lead to scouring at these locations. The maximum turbulent production rate occurs in the wake region of the bar. The high rate of turbulence production has occurred in that region, which can be ascribed to the process of shedding turbulent vortices. The results show that the impact of the presence of the bar is mainly restricted to the lower layers of flow. The turbulent dissipation rate monotonically decreases with an increase in the vertical distance from the bed. The turbulent production rate first increases and then decreases with successive increases in the vertical distance from the bed. The paper concludes with suggestions for the future potential use of the present research for the practical purpose of examining braid bar occurrences in alluvial rivers to develop an appropriate response through training measures

Yes / In a laboratory, a model of a mid-channel bar is built to study the turbulent flow structures in its vicinity. The present study on the turbulent flow structure around a mid-channel bar is based on unravelling the fluvial fluxes triggered by the bar’s 3D turbulent burst phenomenon. To this end, the three-dimensional velocity components are measured with the help of acoustic doppler velocimetry (ADV). The results indicate that the transverse component of turbulent kinetic energy cannot be neglected when analyzing turbulent burst processes, since the dominant flow is three-dimensional around the mid-channel bar. Due to the three-dimensionality of flow, the octant events approach is used for analyzing the flow in the vicinity of the mid-channel bar. The aim is to develop functional relationships between the stable movements that are modelled in the present study. To find the best Markov chain order to present experimental datasets, the zero-, first-, and second-order Markov chains are analyzed using the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). The parameter transition ratio has evolved in this research to reflect the linkage of streambed elevation changes with stable transitional movements. For a better understanding of the temporal behaviors of stable transitional movements, the residence time vs. frequency graphs are also plotted for scouring as well as for depositional regions. The study outcome herein underlines the usefulness of the octant events approach for characterizing turbulent bursts around mid-channel bar formation, which is a precursor to the initiation of braiding configuration. / This research and APC was funded by King Saud University, Riyadh, Saudi Arabia through Researchers Supporting Project number (RSP-2021/297).

Yes / The aquatic vegetation patch plays a significant role on sediment net deposition in the vegetated channels. Particularly, the flow is decelerated at the leading edge of a patch that tends to induce vertical updraft, that is, a diverging flow region, in which vegetation greatly affects the pattern of sediment net deposition. This study focuses on the simulation of the sediment net deposition in the whole vegetation patch region through an innovative random displacement model, a Lagrange method, with probability-based boundary conditions, instead of the reflection or sorption boundary at the channel bottom. The probability model of deposition and resuspension is proposed according to the flow field characteristics in the different regions of the vegetation patch. The variation of the sediment deposition and resuspension with the turbulent kinetic energy is analyzed to illustrate the effect of the turbulence induced by vegetation, represented by the dimensionless turbulent kinetic energy (ψ), on the sediment deposition and resuspension. The sediment deposition predicted by the proposed model agrees well with the experimental measurements. Results show that the effect of vegetation on the sediment deposition and resuspension motions begins to prevail when the vegetation-induced ψ is larger than its threshold, ψ *. The threshold of ψ is predicted to be within 6.8–10 according to the simulation results in this study. As the turbulent kinetic energy increases, the deposition probability decreases continuously when ψ > ψ *. / ational Natural Science Foundation of China (NSFC). Grant Numbers: 52020105006, 11872285; UK Royal Society - International Exchanges Program. Grant Number: IES\R2\181122; Open Funding of State Key Laboratory of Water Resources and Hydropower Engineering Science (WRHES), Wuhan University. Grant Number: 2018HLG01

This research is intended to determine the role of riparian vegetation in stream morphology. This experiment examined the effects of riparian vegetation on boundary shear stress (BSS) by completing the following objectives: (1) evaluating the effects of streambank vegetation on near-bank velocity and turbulence; (2) determining a method for measuring BSS; and, (3) examining the effects of streambank vegetation on BSS using an existing model. A second order prototype stream, with individual reaches dominated by the three vegetation types (trees, shrubs, and grass) was modeled using a fixed-bed Froude-scale modeling technique. One model streambank of the prototype stream was constructed for each vegetation type in addition to one bank with only grain roughness. Velocity profiles were measured using an acoustic Doppler velocimeter (ADV) and a miniature propeller (MP). A flush-mounted Dantec MiniCTA system was used to measure shear stress at the streambank wall. The addition of vegetation on a sloping streambank increased the streamwise free stream velocity and decreased the near-bank streamwise velocity. The turbulence caused by the upright shrub treatment increased turbulent kinetic energy and Reynolds stresses near the streambank toe, an area susceptible to fluvial erosion. The presence of dense, semi-rigid vegetation may encourage the formation of a wider channel with a vertical streambank. The small range of CTA shear stress measurements (0.02—2.14 Pa) suggested that one estimate can describe a streambank. The law of the wall technique is not appropriate because the velocity profiles did not follow the necessary logarithmic shape. Vegetative roughness present in channels created secondary flow; turbulence characteristics more appropriately estimated BSS. The BSS model predicted velocity fields in similar distribution to that measured by the ADV and MP. BSS calculated using the ray-isovel method for both velocity measurement devices were different than the measured BSS values, likely due to distortions in the measured velocity field. In general, the predicted BSS distribution increased with water depth and decreased with increasing vegetation density. The predicted BSS at the shrub toe indicated a spike in shear stress consistent with TKE estimates. / Ph. D.

The main function of the heart is to pump blood throughout the cardiovascular system by generating pressure differences created through volume changes. Although the main purpose of the heart and vessels is to lead the flowing blood throughout the body, clinical assessments of cardiac function are usually based on morphology, approximating the flow features by viewing the motion of the myocardium and vessels. Measurement of three-directional, three-dimensional and time-resolved velocity (4D Flow) data is feasible using magnetic resonance (MR). The focus of this thesis is the development and application of methods that facilitate the analysis of larger groups of data in order to increase our understanding of intracardiac flow patterns and take the 4D flow technique closer to the clinical setting. In the first studies underlying this thesis, a pathline based method for analysis of intra ventricular blood flow patterns has been implemented and applied. A pathline is integrated from the velocity data and shows the path an imaginary massless particle would take through the data volume. This method separates the end-diastolic volume (EDV) into four functional components, based on the position for each individual pathline at end-diastole (ED) and end-systole (ES). This approach enables tracking of the full EDV over one cardiac cycle and facilitates calculation of parameters such as e.g. volumes and kinetic energy (KE). Besides blood flow, pressure plays an important role in the cardiac dynamics. In order to study this parameter in the left ventricle, the relative pressure field was computed using the pressure Poisson equation. A comprehensive presentation of the pressure data was obtained dividing the LV blood pool into 17 pie-shaped segments based on a modification of the standard seventeen segment model. Further insight into intracardiac blood flow dynamics was obtained by studying the turbulent kinetic energy (TKE) in the LV. The methods were applied to data from a group of healthy subjects and patients with dilated cardiomyopathy (DCM). DCM is a pathological state where the cardiac function is impaired and the left ventricle or both ventricles are dilated. The validation study of the flow analysis method showed that a reliable user friendly tool for intra ventricular blood flow analysis was obtained. The application of this tool also showed that roughly one third of the blood that enters the LV, directly leaves the LV again in the same heart beat. The distribution of the four LV EDV components was altered in the DCM group as compared to the healthy group; the component that enters and leaves the LV during one cardiac cycle (Direct Flow) was significantly larger in the healthy subjects. Furthermore, when the kinetic energy was normalized by the volume for each component, at time of ED, the Direct Flow had the highest values in the healthy subjects. In the DCM group, however, the Retained Inflow and Delayed Ejection Flow had higher values. The relative pressure field showed to be highly heterogeneous, in the healthy heart. During diastole the predominate pressure differences in the LV occur along the long axis from base to apex. The distribution and variability of 3D pressure fields differ between early and late diastolic filling phases, but common to both phases is a relatively lower pressure in the outflow segment. In the normal LV, TKE values are low. The highest TKE values can be seen during early diastole and are regionally distributed near the basal LV regions. In contrast, in a heterogeneous group of DCM patients, total diastolic and late diastolic TKE values are higher than in normals, and increase with the LV volume. In conclusion, in this thesis, methods for analysis of multidirectional intra cardiac velocity data have been obtained. These methods allow assessment of data quality, intra cardiac blood flow patterns, relative pressure fields, and TKE. Using these methods, new insights have been obtained in intra cardiac blood flow dynamics in health and disease. The work underlying this thesis facilitates assessment of data from a larger population of healthy subjects and patients, thus bringing the 4D Flow MRI technique closer to the clinical setting.

Flow mechanisms of compound meandering channels are recognised to be far more complicated than compound straight channels. The compound meandering channels are mainly characterised by the continuous variation of mean and turbulent flow parameters along a meander wavelength; the existence of horizontal shear layer at the bankfull level and the presence of strong helical secondary flow circulations in the streamwise direction. The secondary flow circulations are very important as they govern the advection of flow momentum, distort isovels, and influence bed shear stress, thus producing a complicated and fully three-dimensional turbulent flow structures. A great deal of experiments has been conducted in the past, which explains flow mechanisms, mixing patterns and the behaviour of secondary flow circulations. However, a complete understanding of secondary flow structures still remains far from conclusive mainly because the secondary flow structures are influenced by the host of geometrical and flow parameters, which are yet to be investigated in detail. The three-dimensional Reynolds-averaged Navier-Stokes and continuity equations were solved using a standard Computational Fluid Dynamics solver to predict mean velocity, secondary flow and turbulent kinetic energy. Five different flow cases of various model scales and relative depths were considered. Detailed analyses of the measured and predicted flow variables were carried out to understand mean flow mechanisms and turbulent secondary flow structures in compound meandering channels. The streamwise vorticity equation was used to quantify the complex and three-dimensional behaviour of secondary flow circulations in terms of their generation, development and decay along the half-meander wavelength. The turbulent kinetic energy equation was used to understand energy expense mechanisms of secondary flow circulations. The strengths of secondary flow circulations were calculated and compared for different flow cases considered. The main findings from this research are as follows. The shearing of the main channel flow as the floodplain flow plunges into and over the main channel influences the mean and turbulent flow structures particularly in the crossover region. The horizontal shear layer at the inner bankfull level generates secondary flow circulations. As the depth of flow increases, the point of generation of secondary flow circulations moves downstream. The secondary shear stress significantly contributes towards the generation of streamwise vorticity and the production of turbulent kinetic energy. The rate of turbulence kinetic energy production was found to be higher than the rate of its dissipation in the crossover region, which demonstrates that the turbulence extracts more energy from the mean flu\\' than what is actually dissipated. This also implies that, in the crossover region, the turbulence is always advected downstream by the mean and secondary flows, The strength of geometry induced secondary flow circulation increases with the increase in the relative depth.

Le décrochage est un phénomène aérodynamique instationnaire susceptible d'apparaître sur de nombreux profils aérodynamiques. Il résulte d'un décollement important de l'écoulement vis-à-vis de la paroi de l'aile et dégrade considérablement les performances de vol. Sur certains profils de pales d'hélicoptères, d'éoliennes ou de rotors, ce phénomène se produit dans des conditions d'utilisation normales et justifie la recherche de méthodes de modélisation accessibles industriellement. Le décrochage est initié au bord d'attaque par l'apparition d'une petite région de recirculation de fluide appelée bulbe de décollement laminaire où l'écoulement transitionne de l'état laminaire vers l'état turbulent. Ce phénomène encore mal connu met en jeu transition et écoulements hors équilibre auxquels les outils de modélisation RANS habituellement employés ne sont pas adaptés. Dans cette étude, un bulbe transitionnel typique d'un écoulement de bord d'attaque de pale d'hélicoptère (profil OA209 à un nombre de Reynolds Rec∞=1.8x106 et 15° d'incidence) est isolé sur une plaque plane. Une simulation DNS de cet écoulement est réalisée à l'aide du logiciel FUNk de l'ONERA afin de servir de base de données pour l'amélioration des modèles RANS. L'évolution des bilans de l'équation de transport de l'énergie cinétique turbulente ainsi que les principales hypothèses RANS (isotropie de la turbulence, Boussinesq, équilibre production/dissipation) sont analysées. Une étude des principaux modèles RANS développés dans le logiciel elsA de l'ONERA est ensuite réalisée en pondérant les grandeurs turbulentes par une fonction de transition reproduisant l'intermittence de la turbulence. Le modèle k-ω de Wilcox couplé à une fonction de transition optimisée a donné les résultats les plus proches de la DNS et a donc été l'objet d'une analyse plus approfondie, notamment une évaluation des principales équations bilans et une application de ce modèle et de sa méthode de transition à un cas de transition naturelle de plaque plane. / The stall is an unsteady aerodynamic phenomenon that may occur on many aerodynamic profiles. It consists in a large separation of the flow from the wall of the wing and significantly deteriorates the flight performances. On some blade profiles such as helicopters, turbines or rotors, this phenomenon occurs under normal conditions of use and justifies the research of industrially accessible modeling methods. The stall is initiated at the leading edge by the appearance of a small region of fluid recirculation called a “laminar separation bubble” where the flow transitions from the laminar to the turbulent state. This still poorly understood phenomenon involves transition and non-equilibrium flows for which commonly used RANS modeling tools are not suitable. In this study, a transitional bubble typical of an helicopter leading edge flow (OA209 profile at a Reynolds number Rec∞=1.8x106 and 15° of incidence) is reproduced on a flat plate. A DNS simulation of this flow is performed using the ONERA FUNk software to serve as a database for RANS models improvements. The evolution of turbulent kinetic energy budgets as well as the main RANS assumptions (isotropy of turbulence, Boussinesq hypothesis, production/dissipation balance) are analyzed. The main RANS models developed in the ONERA elsA software are then studied by weighting the turbulent quantities with a transition function reproducing the intermittency of the turbulence. The k-ω Wilcox model coupled with an optimized transition function gave the best results and was therefore kept for a more in-depth analysis, including an assessment of the main budgets and an application of this model and its transition method to a natural transition test case on a flat plate.

The reliable prediction of turbulent non-equilibrium flows is of high academic and industrial interest in several engineering fields. Most turbulent flows are often predicted using single-time-scale Reynolds-Averaged-Navier-Stokes (RANS) turbulence models which assume the flows can be modelled through a single time or length scale which is an admittedly incorrect assumption. Therefore they are not expected to capture the lag in the response of the turbulence in non-equilibrium flows. In attempts to improve prediction of these flows, by taking into consideration some features of the turbulent kinetic energy spectrum, the multiple-time-scale models arose. A number of two-scale models have been proposed, but so far their use has been rather limited.This work thus focusses on the development of two-time-scale approaches. Two two-time-scale linear-eddy-viscosity models, referred to as NT1 and NT2 models, have been developed and the initial stages of the development of two-time-scale non-linear-eddy-viscosity models are also reported. The models' coefficients have been determined through asymptotic analysis of decaying grid turbulence, homogeneous shear flows and the flow in a boundary layer in local equilibrium. Three other important features of these models are that there is consistent partition of the large and the small scales for all above limiting cases, model sensitivity to the partition and production rate ratios and sensitivity of the eddy viscosity sensitive to the mean strain rates.The models developed have been tested through computations of a wide range of flows such as homogeneous shear and normally strained flows, fully developed channel flows, zero-pressure-gradient, adverse-pressure-gradient, favourable-pressure-gradient and oscillatory boundary layer flows, fully developed oscillatory and ramp up pipe flows and steady and pulsated backward-facing-step flows.The proposed NT1 and NT2 two-scale models have been shown to perform well in all test cases, being, among the benchmarked models tested, the models which best performed in the wide range of dimensionless shear values of homogeneous shear flows, the only linear-eddy-viscosity models which predicted well the turbulent kinetic energy in the normally strained cases and the only models which showed satisfactory sensitivity in predicting correctly the reattachment point in the unsteady backward facing step cases with different forcing frequencies. Although the development of the two-time-scale non-linear-eddy-viscosity models is still in progress, the interim versions proposed here have resulted in predictions of the Reynolds normal stresses similar to those of much more complex models in all test cases studied and in predictions of the turbulent kinetic energy in normally strained flows which are better than those of the other models tested in this study.

No description available.