Physical oceanography and processed 2 minutes-averaged continuous VM-ADCP profiles during POLARSTERN cruise ANT-XXVIII/3

The physical and biological carbon pumps in the different hydrographic and biogeochemical regimes of the Atlantic Sector of the Southern Ocean are controlled by a series of coupled physical, chemical and biological processes and a project named Eddy-Pump was designed to study them. The Eddy Pump field campaign was carried out during RV Polarstern Cruise ANT-XXVIII/3 between January and March 2012. Particular emphasis was laid on the differences which occur along the axis of the Antarctic Circumpolar Current (ACC) with its associated mesoscale eddy field. The study sites were selected in order to represent (1) the central ACC with its regular separation in different frontal jets, investigated by a meridional transect along 10°E; (2) a large-scale bloom west of the Mid-Atlantic Ridge which lasted several months with conspicuous chlorophyll-poor waters to its immediate east studied by a three-dimensional mesoscale survey centred at 12°40'W; and (3) the Georgia Basin north of the island of South Georgia, which regularly features an extended and dense phytoplankton bloom, was investigated by a mesoscale survey centred at 38°12'W. While Eddy-Pump represents an interdisciplinary project by design, we here focus on describing the variable physical environment within which the different biogeochemical regimes developed. For describing the physical environment we use measurements of temperature, salinity and density, of mixed-layer turbulence parameters, of dynamic heights and horizontal current vectors, and of flow trajectories obtained from surface drifters and submerged floats. This serves as background information for the analyses of biological and chemical processes and of biogeochemical fluxes addressed by other papers in this issue. The section along 10°E between 44°S and 53°S showed a classical ACC structure with well-known hydrographic fronts, the Subantarctic Front (SAF) at 46.5°S, the Antarctic Polar Front (APF) split in two, at 49.25°S and 50.5°S, and the Southern Polar Front (SPF) at 52.5°S. Each front was associated with strong eastward flows. The West Mid-Atlantic Ridge Survey showed a weak and poorly resolved meander structure between the APF and the SPF. During the first eight days of the survey the oceanographic conditions at the Central Station at 12°40'W remained reasonably constant. However after that, conditions became more variable in the thermocline with conspicuous temperature inversions and interleavings and also a decrease in temperature in the surface layer. At the very end of the period of observation the conditions in the thermocline returned to being similar to those observed during the early part of the period with however the mixed layer temperature raised. The period of enhanced thermohaline variability was accompanied by increased currents. The Georgia Basin Survey showed a very strong zonal jet at its northern edge which connects to a large cyclonic meander that itself joins an anticyclonic eddy in the southeastern quadrant. The water mass contrasts in this survey were stronger than in the West Mid-Atlantic Ridge Survey, but similar to those met along 10°E with the exception that the warm and saline surface water typical of the northern side of the SAF was not covered by the Georgia Basin Survey. Mixed layers found during Eddy-Pump were typically deep, but varied between the three survey areas; the mean depths and standard variations of the mixed layer along the 10°E were 77.2±24.7 m, at the West Mid-Atlantic Ridge 66.7±17.7 m, and in the Georgia Basin 36.8±10.7 m.

Upper-Pleistocene terrace deposits in Mediterranean climate: geomorphological and source-rock control on mineral and geochemical signatures (Betic Cordillera, SE Spain)

Mineral and chemical composition of alluvial Upper-Pleistocene deposits from the Alto Guadalquivir Basin (SE Spain) were studied as a tool to identify sedimentary and geomorphological processes controlling its formation. Sediments located upstream, in the north-eastern sector of the basin, are rich in dolomite, illite, MgO and KB2BO. Downstream, sediments at the sequence base are enriched in calcite, smectite and CaO, whereas the upper sediments have similar features to those from upstream. Elevated rare-earth elements (REE) values can be related to low carbonate content in the sediments and the increase of silicate material produced and concentrated during soil formation processes in the neighbouring source areas. Two mineralogical and geochemical signatures related to different sediment source areas were identified. Basal levels were deposited during a predominantly erosive initial stage, and are mainly composed of calcite and smectite materials enriched in REE coming from Neogene marls and limestones. Then the deposition of the upper levels of the alluvial sequences, made of dolomite and illitic materials depleted in REE coming from the surrounding Sierra de Cazorla area took place during a less erosive later stage of the fluvial system. Such modification was responsible of the change in the mineralogical and geochemical composition of the alluvial sediments.

Modeling of Turbulent Natural Gas and Biogas Flames of the Delft Jet-in-Hot-Coflow Burner: Effects of Coflow Temperature, Fuel Temperature and Fuel Composition on the Flame Lift-Off Height

The structure of a turbulent non-premixed flame of a biogas fuel in a hot and diluted coflow mimicking moderate and intense low dilution (MILD) combustion is studied numerically. Biogas fuel is obtained by dilution of Dutch natural gas (DNG) with CO2. The results of biogas combustion are compared with those of DNG combustion in the Delft Jet-in-Hot-Coflow (DJHC) burner. New experimental measurements of lift-off height and of velocity and temperature statistics have been made to provide a database for evaluating the capability of numerical methods in predicting the flame structure. Compared to the lift-off height of the DNG flame, addition of 30 % carbon dioxide to the fuel increases the lift-off height by less than 15 %. Numerical simulations are conducted by solving the RANS equations using Reynolds stress model (RSM) as turbulence model in combination with EDC (Eddy Dissipation Concept) and transported probability density function (PDF) as turbulence-chemistry interaction models. The DRM19 reduced mechanism is used as chemical kinetics with the EDC model. A tabulated chemistry model based on the Flamelet Generated Manifold (FGM) is adopted in the PDF method. The table describes a non-adiabatic three stream mixing problem between fuel, coflow and ambient air based on igniting counterflow diffusion flamelets. The results show that the EDC/DRM19 and PDF/FGM models predict the experimentally observed decreasing trend of lift-off height with increase of the coflow temperature. Although more detailed chemistry is used with EDC, the temperature fluctuations at the coflow inlet (approximately 100K) cannot be included resulting in a significant overprediction of the flame temperature. Only the PDF modeling results with temperature fluctuations predict the correct mean temperature profiles of the biogas case and compare well with the experimental temperature distributions.

Comparison of DES predictions of turbulent flow and heat transfer in a 45° ribbed duct

Numerical predictions of the turbulent flow and heat transfer of a stationary duct with square ribs 45° angled to the main flow direction are presented. The rib height to channel hydraulic diameter is 0.1, the rib pitch to rib height is 10. The calculations have been carried out for a bulk Reynolds number of 50,000. The flows generated by ribs are dominated by separating and reattaching shear layers with vortex shedding and secondary flows in the cross-section. The hybrid RANS-LES approach is adopted to simulate such flows at a reasonable computation cost. The capability of the various versions of DES method, depending the RANS model, such as DES-SA, DES-RKE, DES-SST, have been compared and validated against the experiment. The significant effect of RANS model on the accuracy of the DES prediction has been shown. The DES-SST method, which was able to reproduce the correct physics of flow and heat transfer in a ribbed duct showed better performance than others.

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.

Desenvolvimento de um método para correção de curvas de desempenho em bombas centrífugas submersas operando com fluidos viscosos

Electrical Submersible Pump (ESP) is used as an artificial lift technique. However, pumping viscous oil is generally associated with low Reynolds number flows. This condition leads to a performance degradation respect to the performance expected from the regular operation with water that most of the centrifugal pumps are originally designed for. These issues are considered in this investigation through a numerical study of the flow in two different multistage, semi-axial type ESPs. This investigation is carried out numerically using a Computational Fluid Dynamics (CFD) package, where the transient RANS equations are solved numerically. The turbulence is modeled using the SST model. Head curves for several operating conditions are compared with manufacturer’s curves and experimental data for a three-stage ESP, showing good agreement for a wide range of fluid viscosities and rotational speeds. Dimensionless numbers (n, n, n e Re) are used to investigate performance degradation of the ESPs. In addition, flow phenomena through the impellers of the ESPs are investigated using flow field from numerical results. Results show that performance degradation is directly related to rotational Reynolds number, Re. In addition, it was verified that performance degradation occurs for constant normalized specific speedn, which shows that performance degradation occurs similarly for different centrifugal pumps. Moreover, experimental data and numerical results agreed with a correlation from literature between head and flow correction factors proposed by Stepanoff (1967). A definition of modified Reynolds number was proposed and relates the head correction factor to viscosity. A correlation between head correction factor and the modified Reynolds number was proposed, which agreed well with numerical and experimental data. Then, a method to predict performance degradation based on the previous correlations was proposed. This method was compared with others from literature. In general, results and conclusions from this work can also be useful to bring more information about the flow of highly viscous fluids in pumps, especially in semi-axial, multistage ESPs.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.

Pressure distribution on a body a revolution in supersonic flight, by Robert R. Reynolds.

"This report appears as an appendix to the National Applied Mathematics Laboratories Quarterly report of projects and publications, January through March 1952."

Heat transfer rates for cross flow of water through a tube bank at high Reynolds numbers.

Mode of access: Internet.

Aerosol transport in a bifurcating tube at cyclic flow

Branching tube flow is a common feature of many fields of science and technology, and occurs both in animate and inanimate systems [1]. The transport of aerosol particles is of particular importance in industrial flow networks but also for the respiratory tree [2]. In this analysis a 3-D numerical study is performed to investigate transport and deposition of aerosol particles in branching tubes. Bifurcation tubes designed according to Hess-Murray law [3] but with different branching angles are analyzed. This study covers cyclic flow conditions at frequencies of 0.25 Hz, 0.50 Hz and 0.75 Hz, Stokes numbers ranging between 0.03 and 0.25, and Reynolds numbers up to 3000.

Two-dimensional Stokes flow driven by elliptical paddles

A fast and accurate numerical technique is developed for solving the biharmonic equation in a multiply connected domain, in two dimensions. We apply the technique to the computation of slow viscous flow (Stokes flow) driven by multiple stirring rods. Previously, the technique has been restricted to stirring rods of circular cross section; we show here how the prior method fails for noncircular rods and how it may be adapted to accommodate general rod cross sections, provided only that for each there exists a conformal mapping to a circle. Corresponding simulations of the flow are described, and their stirring properties and energy requirements are discussed briefly. In particular the method allows an accurate calculation of the flow when flat paddles are used to stir a fluid chaotically.