Large eddy simulation of shock boundary layer interaction control using micro-vortex generators

The performance of supersonic engine inlets and external aerodynamic surfaces can be critically affected by shock wave / boundary layer interactions (SBLIs), whose severe adverse pressure gradients can cause boundary layer separation. Currently such problems are avoided primarily through the use of boundary layer bleed/suction which can be a source of significant performance degradation. This study investigates a novel type of flow control device called micro-vortex generators (µVGs) which may offer similar control benefits without the bleed penalties. µVGs have the ability to alter the near-wall structure of compressible turbulent boundary layers to provide increased mixing of high speed fluid which improves the boundary layer health when subjected to flow disturbance. Due to their small size,µVGs are embedded in the boundary layer which provide reduced drag compared to the traditional vortex generators while they are cost-effective, physically robust and do not require a power source. To examine the potential of µVGs, a detailed experimental and computational study of micro-ramps in a supersonic boundary layer at Mach 3 subjected to an oblique shock was undertaken. The experiments employed a flat plate boundary layer with an impinging oblique shock with downstream total pressure measurements. The moderate Reynolds number of 3,800 based on displacement thickness allowed the computations to use Large Eddy Simulations without the subgrid stress model (LES-nSGS). The LES predictions indicated that the shock changes the structure of the turbulent eddies and the primary vortices generated from the micro-ramp. Furthermore, they generally reproduced the experimentally obtained mean velocity profiles, unlike similarly-resolved RANS computations. The experiments and the LES results indicate that the micro-ramps, whose height is h≈0.5δ, can significantly reduce boundary layer thickness and improve downstream boundary layer health as measured by the incompressible shape factor, H. Regions directly behind the ramp centerline tended to have increased boundary layer thickness indicating the significant three-dimensionality of the flow field. Compared to baseline sizes, smaller micro-ramps yielded improved total pressure recovery. Moving the smaller ramps closer to the shock interaction also reduced the displacement thickness and the separated area. This effect is attributed to decreased wave drag and the closer proximity of the vortex pairs to the wall. In the second part of the study, various types of µVGs are investigated including micro-ramps and micro-vanes. The results showed that vortices generated from µVGs can partially eliminate shock induced flow separation and can continue to entrain high momentum flux for boundary layer recovery downstream. The micro-ramps resulted in thinner downstream displacement thickness in comparison to the micro-vanes. However, the strength of the streamwise vorticity for the micro-ramps decayed faster due to dissipation especially after the shock interaction. In addition, the close spanwise distance between each vortex for the ramp geometry causes the vortex cores to move upwards from the wall due to induced upwash effects. Micro-vanes, on the other hand, yielded an increased spanwise spacing of the streamwise vortices at the point of formation. This resulted in streamwise vortices staying closer to the wall with less circulation decay, and the reduction in overall flow separation is attributed to these effects. Two hybrid concepts, named “thick-vane” and “split-ramp”, were also studied where the former is a vane with side supports and the latter has a uniform spacing along the centerline of the baseline ramp. These geometries behaved similar to the micro-vanes in terms of the streamwise vorticity and the ability to reduce flow separation, but are more physically robust than the thin vanes. Next, Mach number effect on flow past the micro-ramps (h~0.5δ) are examined in a supersonic boundary layer at M=1.4, 2.2 and 3.0, but with no shock waves present. The LES results indicate that micro-ramps have a greater impact at lower Mach number near the device but its influence decays faster than that for the higher Mach number cases. This may be due to the additional dissipation caused by the primary vortices with smaller effective diameter at the lower Mach number such that their coherency is easily lost causing the streamwise vorticity and the turbulent kinetic energy to decay quickly. The normal distance between the vortex core and the wall had similar growth indicating weak correlation with the Mach number; however, the spanwise distance between the two counter-rotating cores further increases with lower Mach number. Finally, various µVGs which include micro-ramp, split-ramp and a new hybrid concept “ramped-vane” are investigated under normal shock conditions at Mach number of 1.3. In particular, the ramped-vane was studied extensively by varying its size, interior spacing of the device and streamwise position respect to the shock. The ramped-vane provided increased vorticity compared to the micro-ramp and the split-ramp. This significantly reduced the separation length downstream of the device centerline where a larger ramped-vane with increased trailing edge gap yielded a fully attached flow at the centerline of separation region. The results from coarse-resolution LES studies show that the larger ramped-vane provided the most reductions in the turbulent kinetic energy and pressure fluctuation compared to other devices downstream of the shock. Additional benefits include negligible drag while the reductions in displacement thickness and shape factor were seen compared to other devices. Increased wall shear stress and pressure recovery were found with the larger ramped-vane in the baseline resolution LES studies which also gave decreased amplitudes of the pressure fluctuations downstream of the shock.

Detached eddy simulation of unsteady turbulent flows in the draft tube of a bulb turbine

Les aspirateurs de turbines hydrauliques jouent un rôle crucial dans l’extraction de l’énergie disponible. Dans ce projet, les écoulements dans l’aspirateur d’une turbine de basse chute ont été simulés à l’aide de différents modèles de turbulence dont le modèle DDES, un hybride LES/RANS, qui permet de résoudre une partie du spectre turbulent. Déterminer des conditions aux limites pour ce modèle à l’entrée de l’aspirateur est un défi. Des profils d’entrée 1D axisymétriques et 2D instationnaires tenant compte des sillages et vortex induits par les aubes de la roue ont notamment été testés. Une fluctuation artificielle a également été imposée, afin d’imiter la turbulence qui existe juste après la roue. Les simulations ont été effectuées pour deux configurations d’aspirateur du projet BulbT. Pour la deuxième, plusieurs comparaisons avec des données expérimentales ont été faites pour deux conditions d’opération, à charge partielle et dans la zone de baisse rapide du rendement après le point de meilleur rendement. Cela a permis d’évaluer l’efficacité et les lacunes de la modélisation turbulente et des conditions limites à travers leurs effets sur les quantités globales et locales. Les résultats ont montrés que les structures tourbillonnaires et sillages sortant de la roue sont adéquatement résolus par les simulations DDES de l’aspirateur, en appliquant les profils instationnaires bidimensionnels et un schéma de faible dissipation pour le terme convectif. En outre, les effets de la turbulence artificielle à l’entrée de l’aspirateur ont été explorés à l’aide de l’estimation de l’intermittence du décollement, de corrélations en deux points, du spectre d’énergie et du concept de structures cohérentes lagrangiennes. Ces analyses ont montré que les détails de la dynamique de l’écoulement et de la séparation sont modifiés, ainsi que les patrons des lignes de transport à divers endroits de l’aspirateur. Cependant, les quantités globales comme le coefficient de récupération de l’aspirateur ne sont pas influencées par ces spécificités locales.

DopSCAT: A mission concept for simultaneous measurements of marine winds and surface currents

A radar scatterometer operates by transmitting a pulse of microwave energy toward the ocean's surface and measuring the normalized (per-unit-surface) radar backscatter coefficient (σ°). The primary application of scatterometry is the measurement of near-surface ocean winds. By combining σ° measurements from different azimuth angles, the 10 m vector wind can be determined through a Geophysical Model Function (GMF), which relates wind and backscatter. This paper proposes a mission concept for the measurement of both oceanic winds and surface currents, which makes full use of earlier C-band radar remote sensing experience. For the determination of ocean currents, in particular, the novel idea of using two chirps of opposite slope is introduced. The fundamental processing steps required to retrieve surface currents are given together with their associated accuracies. A detailed description of the mission proposal and comparisons between real and retrieved surface currents are presented. The proposed ocean Doppler scatterometer can be used to generate global surface ocean current maps with accuracies better than 0.2 m/s at a spatial resolution better than 25 km (i.e., 12.5 km spatial sampling) on a daily basis. These maps will allow gaining some insights on the upper ocean mesoscale dynamics. The work lies at a frontier, given that the present inability to measure ocean currents from space in a consistent and synoptic manner represents one of the greatest weaknesses in ocean remote sensing.

Origin and fate of surface drift in the oceanic convergence zones of the eastern Pacific

This study investigates the structure and intensity of the surface pathways connecting to and from the central areas of the large-scale convergence regions of the eastern Pacific Ocean. Surface waters are traced with numerical Lagrangian particles transported in the velocity field of three different ocean models with horizontal resolutions that range from ¼° to 1/32°. The connections resulting from the large-scale convergent Ekman dynamics agree qualitatively but are strongly modulated by eddy variability that introduces meridional asymmetry in the amplitude of transport. Lagrangian forward-in-time integrations are used to analyze the fate of particles originating from the central regions of the convergence zones and highlight specific outflows not yet reported for the southeastern Pacific when using the currents at the highest resolutions (1/12° and 1/32°). The meridional scales of these outflows are comparable to the characteristic width of the fine-scale striation of mean currents.

Numerical study of circulation in the Caspian Sea and its effects on the southern coastal currents

In this study, numerical simulation of the Caspian Sea circulation was performed using COHERENS three-dimensional numerical model and field data. The COHERENS three-dimensional model and FVCOM were performed under the effect of the wind driven force, and then the simulation results obtained were compared. Simulation modeling was performed at the Caspian Sea. Its horizontal grid size is approximately equal to 5 Km and 30 sigma levels were considered. The numerical simulation results indicate that the winds' driven-forces and temperature gradient are the most important driving force factors of the Caspian circulation pattern. One of the effects of wind-driven currents was the upwelling phenomenon that was formed in the eastern shores of the Caspian Sea in the summer. The simulation results also indicate that this phenomenon occurred at a depth less than 40 meters, and the vertical velocity in July and August was 10 meters and 7 meters respectively. During the upwelling phenomenon period the temperatures on the east coast compared to the west coast were about 5°C lower. In autumn and winter, the warm waters moved from the south east coast to the north and the cold waters moved from the west coast of the central Caspian toward the south. In the subsurface and deep layers, these movements were much more structured and caused strengthening of the anti-clockwise circulation in the area, especially in the central area of Caspian. The obtained results of the two models COHERENS and FVCOM performed under wind driven-force show a high coordination of the two models, and so the wind current circulation pattern for both models is almost identical.

How essential are Argo observations to constrain a global ocean data assimilation system?

Observing system experiments (OSEs) are carried out over a 1-year period to quantify the impact of Argo observations on the Mercator Ocean 0.25° global ocean analysis and forecasting system. The reference simulation assimilates sea surface temperature (SST), SSALTO/DUACS (Segment Sol multi-missions dALTimetrie, d'orbitographie et de localisation précise/Data unification and Altimeter combination system) altimeter data and Argo and other in situ observations from the Coriolis data center. Two other simulations are carried out where all Argo and half of the Argo data are withheld. Assimilating Argo observations has a significant impact on analyzed and forecast temperature and salinity fields at different depths. Without Argo data assimilation, large errors occur in analyzed fields as estimated from the differences when compared with in situ observations. For example, in the 0–300 m layer RMS (root mean square) differences between analyzed fields and observations reach 0.25 psu and 1.25 °C in the western boundary currents and 0.1 psu and 0.75 °C in the open ocean. The impact of the Argo data in reducing observation–model forecast differences is also significant from the surface down to a depth of 2000 m. Differences between in situ observations and forecast fields are thus reduced by 20 % in the upper layers and by up to 40 % at a depth of 2000 m when Argo data are assimilated. At depth, the most impacted regions in the global ocean are the Mediterranean outflow, the Gulf Stream region and the Labrador Sea. A significant degradation can be observed when only half of the data are assimilated. Therefore, Argo observations matter to constrain the model solution, even for an eddy-permitting model configuration. The impact of the Argo floats' data assimilation on other model variables is briefly assessed: the improvement of the fit to Argo profiles do not lead globally to unphysical corrections on the sea surface temperature and sea surface height. The main conclusion is that the performance of the Mercator Ocean 0.25° global data assimilation system is heavily dependent on the availability of Argo data.

Release currents of IP₃ receptor channel clusters and concentration profiles

We simulate currents and concentration profiles generated by Ca2+ release from the endoplasmic reticulum (ER) to the cytosol through IP3 receptor channel clusters. Clusters are described as conducting pores in the lumenal membrane with a diameter from 6 nm to 36 nm. The endoplasmic reticulum is modeled as a disc with a radius of 1–12 mm and an inner height of 28 nm. We adapt the dependence of the currents on the trans Ca2+ concentration (intralumenal) measured in lipid bilayer experiments to the cellular geometry. Simulated currents are compared with signal mass measurements in Xenopus oocytes. We find that release currents depend linearly on the concentration of free Ca2+ in the lumen. The release current is approximately proportional to the square root of the number of open channels in a cluster. Cytosolic concentrations at the location of the cluster range from 25 μM to 170 μM. Concentration increase due to puffs in a distance of a few micrometers from the puff site is found to be in the nanomolar range. Release currents decay biexponentially with timescales of < 1 s and a few seconds. Concentration profiles decay with timescales of 0.125–0.250 s upon termination of release.

Hydrothermal-vent alvinellid polychaete dispersal in the eastern Pacific .1. Influence of vent site distribution, bottom currents, and biological patterns

Deep-sea hydrothermal-vent habitats are typically linear, discontinuous, and short-lived. Some of the vent fauna such as the endemic polychaete family Alvinellidae are thought to lack a planktotrophic larval stage and therefore not to broadcast-release their offspring. The genetic evidence points to exchanges on a scale that seems to contradict this type of reproductive pattern. However, the rift valley may topographically rectify the bottom currents, thereby facilitating the dispersal of propagules between active vent sites separated in some cases by 10s of kilometers or more along the ridge axis. A propagule flux model based on a matrix of intersite distances, long-term current-meter data, and information on the biology and ecology of Alvinellidae was developed to test this hypothesis. Calculations of the number of migrants exchanged between two populations per generation (N-m) allowed comparisons with estimates obtained from genetic studies. N, displays a logarithmic decrease with increasing dispersal duration and reaches the critical value of 1 after 8 d when the propagule Aux model was run in standard conditions. At most, propagule traveling time cannot reasonably exceed 15-30 d, according to the model, whereas reported distances between sites would require longer lasting dispersal abilities. Two nonexclusive explanations are proposed. First, some aspects of the biology of Alvinellidae have been overlooked and long-distance dispersal does occur. Second, such dispersal never occurs in Alvinellidae, but the spatial-temporal dynamics of vent sites over geological timescales allows short-range dispersal processes to maintain gene flow.

Tide-topography interaction along the eastern Brazilian shelf

In the Tropics, continental shelves governed by western boundary currents are considered to be among the least productive ocean margins in the world, unless eddy-induced shelf-edge upwelling becomes significant. The eastern Brazilian shelf in the Southwest Atlantic is one of these, and since the slight nutrient input from continental sources is extremely oligotrophic. It is characterized by complex bathymetry with the presence of shallow banks and seamounts. In this work, a full three-dimensional nonlinear primitive equation ocean model is used to demonstrate that the interaction of tidal currents and the bottom topography of the east Brazil continental shelf is capable of producing local upwelling of South Atlantic Central Water, bringing nutrients up from deep waters to the surface layer. Such upper layer enrichment is found to be of significance in increasing local primary productivity. (c) 2005 Elsevier Ltd. All rights reserved.

Pelvic-Floor Rehabilitation, Part 2 : Pelvic-Floor Reeducation With Interferential Currents and Exercise in the Treatment of Genuine Stress Incontinence in Postpartum Women : A Cohort Study

Background and Purpose. This descriptive cohort study investigated a physical therapy program of pelvic-floor neuromuscular electrostimulation (NMES) combined with exercises, with the aim of developing a simple, inexpensive, and conservative treatment for postpartum genuine stress incontinence (GSI). Subjects. Eight female subjects with urodynamically established GSI persisting more than 3 months after delivery participated in the study. The subjects ranged in age from 24 to 37 years (X̅=32, SD=4.2). Methods. This was a descriptive multiple-subject cohort study. Each subject received a total of nine treatment sessions during 3 consecutive weeks, consisting of two 15-minute sessions of NMES followed by a 15-minute pelvic-floor muscle exercise program. Patients also practiced daily pelvic-floor exercises during the 3-week treatment period. The treatment intervention was measured using three separate variables. Maximum muscle contractions (pretraining, during training, and posttraining) were measured indirectly as pressure, using perineometry. Urine loss pretraining and posttraining was measured by means of a Pad test. Self-reported frequency of incontinence was recorded daily throughout the period of the study, using a diary. Data were analyzed using a one-way repeated measures analysis of variance (ANOVA), a Wilcoxon signed-ranks test, and a Friedman two-way ANOVA by ranks. Results. The results indicated that maximum pressure generated by pelvic-floor contractions was greater and both the quantity of urine loss and the frequency of incontinence were lower following the implementation of the physical therapy program. Five subjects became continent, and three others improved. A follow-up survey 1 year later confirmed the consistency of these results. Conclusion and Discussion. The results suggest that the proposed physical therapy program may influence postpartum GSI. Further studies are needed to validate this simple, inexpensive, and conservative physical therapy protocol.

Pelvic-Floor Rehabilitation, Part 1 : Comparison of Two Surface Electrode Placements During Stimulation of the Pelvic-Floor Musculature in Women Who Are Continent Using Bipolar Interferential Currents

Background and Purpose. Electrical stimulation of the pelvic floor is used as an adjunct in the conservative treatment of urinary incontinence. No consensus exists, however, regarding electrode placements for optimal stimulation of the pelvic-floor musculature. The purpose of this study was to compare two different bipolar electrode placements, one suggested by Laycock and Green (L2) the other by Dumoulin (D2), during electrical stimulation with interferential currents of the pelvic-floor musculature in continent women, using a two-group crossover design. Subjects. Ten continent female volunteers, ranging in age from 20 to 39 years (X̅=27.3, SD=5.6), were randomly assigned to one of two study groups. Methods. Each study group received neuromuscular electrical stimulation (NMES) of the pelvic-floor musculature using both electrode placements, the order of application being reversed for each group. Force of contraction was measured as pressure (in centimeters of water [cm H2O]) exerted on a vaginal pressure probe attached to a manometer. Data were analyzed using a two-way, mixed-model analysis of variance. Results. No difference in pressure was observed between the two electrode placements. Differences in current amplitude were observed, with the D2 electrode placement requiring less current amplitude to produce a maximum recorded pressure on the manometer. Subjective assessment by the subjects revealed a preference for the D2 electrode placement (7 of 10 subjects). Conclusion and Discussion. The lower current amplitudes required with the D2 placement to obtain recordings comparable to those obtained with the L2 technique suggest a more comfortable stimulation of the pelvic-floor muscles. The lower current amplitudes required also suggest that greater increases in pressure might be obtained with the D2 placement by increasing the current amplitude while remaining within the comfort threshold. These results will help to define treatment guidelines for a planned clinical study investigating the effects of NMES and exercise in the treatment of urinary stress incontinence in women postpartum.

Subfilter-scale Stress Modelling for Large Eddy Simulations

A subfilter-scale (SFS) stress model is developed for large-eddy simulations (LES) and is tested on various benchmark problems in both wall-resolved and wall-modelled LES. The basic ingredients of the proposed model are the model length-scale, and the model parameter. The model length-scale is defined as a fraction of the integral scale of the flow, decoupled from the grid. The portion of the resolved scales (LES resolution) appears as a user-defined model parameter, an advantage that the user decides the LES resolution. The model parameter is determined based on a measure of LES resolution, the SFS activity. The user decides a value for the SFS activity (based on the affordable computational budget and expected accuracy), and the model parameter is calculated dynamically. Depending on how the SFS activity is enforced, two SFS models are proposed. In one approach the user assigns the global (volume averaged) contribution of SFS to the transport (global model), while in the second model (local model), SFS activity is decided locally (locally averaged). The models are tested on isotropic turbulence, channel flow, backward-facing step and separating boundary layer. In wall-resolved LES, both global and local models perform quite accurately. Due to their near-wall behaviour, they result in accurate prediction of the flow on coarse grids. The backward-facing step also highlights the advantage of decoupling the model length-scale from the mesh. Despite the sharply refined grid near the step, the proposed SFS models yield a smooth, while physically consistent filter-width distribution, which minimizes errors when grid discontinuity is present. Finally the model application is extended to wall-modelled LES and is tested on channel flow and separating boundary layer. Given the coarse resolution used in wall-modelled LES, near the wall most of the eddies become SFS and SFS activity is required to be locally increased. The results are in very good agreement with the data for the channel. Errors in the prediction of separation and reattachment are observed in the separated flow, that are somewhat improved with some modifications to the wall-layer model.