Estimation of the turbulent kinetic energy dissipation rate from 2D-PIV measurements in a vessel stirred by an axial Mixel TTP impeller

Melting of glacier ice enhanced by bursting air bubbles

Feedbacks between ice melt, glacier flow and ocean circulation can rapidly accelerate ice loss at tidewater glaciers and alter projections of sea-level rise. At the core of these projections is a model for ice melt that neglects the fact that glacier ice contains pressurized bubbles of air due to its formation from compressed snow. Current model estimates can underpredict glacier melt at termini outside the region influenced by the subglacial discharge plume by a factor of 10-100 compared with observations. Here we use laboratory-scale experiments and theoretical arguments to show that the bursting of pressurized bubbles from glacier ice could be a source of this discrepancy. These bubbles eject air into the seawater, delivering additional buoyancy and impulses of turbulent kinetic energy to the boundary layer, accelerating ice melt. We show that real glacier ice melts 2.25 times faster than clear bubble-free ice when driven by natural convection in a laboratory setting. We extend these results to the geophysical scale to show how bubble dynamics contribute to ice melt from tidewater glaciers. Consequently, these results could increase the accuracy of modelled predictions of ice loss to better constrain sea-level rise projections globally.

Efficiency of PBT Impellers with Different Blade Cross-Sections

In this study, the mixing power and the axial and radial velocity distributions were determined for the standard PBT45-6 impeller with inclined blades and for three NACA impellers with airfoil blades. On the basis of the velocity distributions obtained by means of PIV techniques, the pumping efficiency and the size of the secondary circulation for the tested impellers were determined. The next stage was to calculate the efficiency of the impeller operation, in which the comparative criteria were the values of the mixing energy in the form of a dimensional and dimensionless criterion. The lowest mixing energy was obtained for impellers with symmetrical profiled NACA0021 blades, which was 20% lower than the similar mixing energy obtained for the standard PBT45-6 impeller.

Elastic turbulence generates anomalous flow resistance in porous media

Many energy, environmental, industrial, and microfluidic processes rely on the flow of polymer solutions through porous media. Unexpectedly, the macroscopic flow resistance often increases above a threshold flow rate in a porous medium, but not in bulk solution. The reason why has been a puzzle for over half a century. Here, by directly visualizing flow in a transparent 3D porous medium, we demonstrate that this anomalous increase is due to the onset of an elastic instability in which the flow exhibits strong spatiotemporal fluctuations reminiscent of inertial turbulence, despite the small Reynolds number. Our measurements enable us to quantitatively establish that the energy dissipated by pore-scale fluctuations generates the anomalous increase in the overall flow resistance. Because the macroscopic resistance is one of the most fundamental descriptors of fluid flow, our results both help deepen understanding of complex fluid flows and provide guidelines to inform a broad range of applications.

Effect of Flow Rate on Turbulence Dissipation Rate Distribution in a Multiphase Pump

The turbulence dissipation will cause the increment of energy loss in the multiphase pump and deteriorate the pump performance. In order to research the turbulence dissipation rate distribution characteristics in the pressurized unit of the multiphase pump, the spiral axial flow type multiphase pump is researched numerically in the present study. This research is focused on the turbulence dissipation rate distribution characteristics in the directions of inlet to outlet, hub to rim, and in the circumferential direction of the rotating impeller blades. Numerical simulation based on the RANS (Reynolds averaged Navier–Stokes equations) and the k-ω SST (Shear Stress Transport) turbulence model has been carried out. The numerical method is verified by comparing the numerical results with the experimental data. Results show that the regions of the large turbulence dissipation rate are mainly at the inlet and outlet of the rotating impeller and static impeller, while it is almost zero from the inlet to the middle of outlet in the suction surface and pressure surface of the first-stage rotating impeller blades. The turbulence dissipation rate is increased gradually from the hub to the rim of the inlet section of the first-stage rotating impeller, while it is decreased firstly and then increased on the middle and outlet sections. The turbulence dissipation rate distributes unevenly in the circumferential direction on the outlet section. The maximum value of the turbulence dissipation rate occurs at 0.9 times of the rated flow rate, while the minimum value at 1.5 times of the rated flow rate. Four turning points in the turbulence dissipation rate distribution that are the same as the number of impeller blades occur at 0.5 times the blade height at 0.9 times the rated flow rate condition. The turbulence dissipation rate distribution characteristics in the pressurized unit of the multiphase pump have been studied carefully in this paper, and the research results have an important significance for improving the performance of the multiphase pump theoretically.

Hybrid Approach for Mixing Time Characterization and Scale-Up in Geometrical Nonsimilar Stirred Vessels Equipped with Eccentric Multi-Impeller Systems—An Industrial Perspective

Multipurpose stirring and blending vessels equipped with various impeller systems are indispensable in the pharmaceutical industry because of the high flexibility necessary during multiproduct manufacturing. On the other hand, process scale-up and scale-down during process development and transfer from bench or pilot to manufacturing scale, or the design of so-called scale-down models (SDMs), is a difficult task due to the geometrical differences of used vessels. The present work comprises a hybrid approach to predict mixing times from pilot to manufacturing scale for geometrical nonsimilar vessels equipped with single top, bottom or multiple eccentrically located impellers. The developed hybrid approach is based on the experimental characterization of mixing time in the dedicated equipment and evaluation of the vessel-averaged energy dissipation rate employing computational fluid dynamics (CFD) using single-phase steady-state simulations. Obtained data are consequently used to develop a correlation of mixing time as a function of vessel filling volume and vessel-averaged energy dissipation rate, which enables the prediction of mixing times in specific vessels based on the process parameters. Predicted mixing times are in good agreement with those simulated using time-dependent CFD simulations for tested operating conditions.

Flow characteristics of the Rushton and pitched blade turbines in turbulent and laminar mixing

The velocity field and concentration distribution of three radial and axial impellers were studied in miscible liquid blending in a stirred tank reactor. To obtain a better insight into the applicability of each impeller system, the influence of velocity profiles on the generated concentration profiles was studied in detail for each impeller. Particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF) methods were used as the measurement techniques. In the turbulent regime, the Rushton turbine provided a homogeneous concentration field throughout the reactor by generating high radial and axial velocities in the bottom and upper zones of the tank. Operation comparisons of the up and down-pumping pitched blade turbines illustrated that the up-pumping turbine had a better performance in overall because of producing high turbulence in the upper half of the reactor. In this region, the radial concentration profile approached the final homogenized value in a shorter time. In the laminar regime, all three impellers acted in a similar way and produced a radial flow throughout the reactor.

Numerical Simulation of Biomass Growth in OKTOP®9000 Reactor at Industrial Scale

Computational fluid dynamics is a powerful method for scale-up of reactors although it is still challenging to fully embrace hydrodynamics and biological complexities. In this article, an aerobic fermentation of Pichia pastoris cells is modeled in a batch OKTOP®9000 reactor. The 800 m³ industrial scale reactor is equipped with a radial impeller, designed by Outotec Oy for gas dispersion in the draft tube reactor. Measured N_p of the impeller is used in hydrodynamics validation. The resolved energy dissipation rate is compensated, and its influence on mass transfer is analyzed and discussed. Gas–liquid drag force is modified to simulate effects of liquid turbulence and bubble swarms. Resolved steady state multiphase hydrodynamics is used to simulate the fermentation process. Temporal evolution of species concentrations is compared to experimental data measured in a small copy of the reactor at lab scale (14 L). The effect of oxygenation on the P. pastoris cells cultivation is considered.

Computational fluid dynamics benchmark dataset of airflow in tracheas

Computational Fluid Dynamics (CFD) is fast becoming a useful tool to aid clinicians in pre-surgical planning through the ability to provide information that could otherwise be extremely difficult if not impossible to obtain. However, in order to provide clinically relevant metrics, the accuracy of the computational method must be sufficiently high. There are many alternative methods employed in the process of performing CFD simulations within the airways, including different segmentation and meshing strategies, as well as alternative approaches to solving the Navier–Stokes equations. However, as in vivo validation of the simulated flow patterns within the airways is not possible, little exists in the way of validation of the various simulation techniques. The data presented here consists of very highly resolved flow data. The degree of resolution is compared to the highest necessary resolutions of the Kolmogorov length and time scales. Therefore this data is ideally suited to act as a benchmark case to which cheaper computational methods may be compared. A dataset and solution setup for one such more efficient method, large eddy simulation (LES), is also presented.

Characterization of Turbulent Properties in the EPA Baffled Flask for Dispersion Effectiveness Testing

The baffled flask test (BFT) has been proposed by United States Environmental Protection Agency to be adopted as the official standard protocol for testing dispersant effectiveness. The mixing energy in the baffled flask is investigated in this paper. Particle image velocimetry (PIV) was used to measure the water velocity in the flask placed at an orbital shaker that was rotated at seven rotation speeds: 100, 125, 150, 160, 170, 200, and 250 rpm. Two dimensional velocity fields in large and small vertical cross sections of the flask for each rotation speed were obtained. The one-dimensional (1D) energy spectra indicates the existence of inertial subrange. The estimated average energy dissipation rates were in the range 7.65×10^-3 to 4 W/kg for rotation speeds of Ω=100-250 rpm, of which it is larger than the one estimated by prior studies using single-point velocity measurement techniques for Ω=100 and 200 rpm. Factors such as instruments used, velocity components measured, and different analysis methods could contribute to the discrepancies in the results. The Kolmogorov scale estimated in this study for all seven rotation speeds approached the size of oil droplets observed at sea, which is 50-400 μm. The average energy dissipation rate, ε and Kolmogorov microscale, η, in the flasks were correlated to the rotation speed, and it was found that ε ¯ = 9.0 × 10 - 5 Exp (0.043Ω) with R ² = 0.97 and η ¯ = 1 , 463 Exp (-0.015Ω) with R ² = 0.98.

Dynamics of airflow in a short inhalation

During a rapid inhalation, such as a sniff, the flow in the airways accelerates and decays quickly. The consequences for flow development and convective transport of an inhaled gas were investigated in a subject geometry extending from the nose to the bronchi. The progress of flow transition and the advance of an inhaled non-absorbed gas were determined using highly resolved simulations of a sniff 0.5 s long, 1 l s⁻¹ peak flow, 364 ml inhaled volume. In the nose, the distribution of airflow evolved through three phases: (i) an initial transient of about 50 ms, roughly the filling time for a nasal volume, (ii) quasi-equilibrium over the majority of the inhalation, and (iii) a terminating phase. Flow transition commenced in the supraglottic region within 20 ms, resulting in large-amplitude fluctuations persisting throughout the inhalation; in the nose, fluctuations that arose nearer peak flow were of much reduced intensity and diminished in the flow decay phase. Measures of gas concentration showed non-uniform build-up and wash-out of the inhaled gas in the nose. At the carina, the form of the temporal concentration profile reflected both shear dispersion and airway filling defects owing to recirculation regions.