New Insights from Locally Resolved Hydrodynamics in Stirred Cell Culture Reactors

The link between hydrodynamics and biological process behavior of antibody-producing mammalian cell cultures is still not fully understood. Common methods to describe dependencies refer mostly to averaged hydrodynamic parameters obtained for individual cultivation systems. In this study, cellular effects and locally resolved hydrodynamics were investigated for impellers with different spatial hydrodynamics. Therefore, the hydrodynamics, mainly flow velocity, shear rate and power input, in a single- and a three-impeller bioreactor setup were analyzed by means of CFD simulations, and cultivation experiments with antibody-producing Chinese hamster ovary (CHO) cells were performed at various agitation rates in both reactor setups. Within the three-impeller bioreactor setup, cells could be cultivated successfully at much higher agitation rates as in the single-impeller bioreactor, probably due to a more uniform flow pattern. It could be shown that this different behavior cannot be linked to parameters commonly used to describe shear effects on cells such as the mean energy dissipation rate or the Kolmogorov length scale, even if this concept is extended by locally resolved hydrodynamic parameters. Alternatively, the hydrodynamic heterogeneity was statistically quantified by means of variance coefficients of the hydrodynamic parameters fluid velocity, shear rate, and energy dissipation rate. The calculated variance coefficients of all hydrodynamic parameters were higher in the setup with three impellers than in the single impeller setup, which might explain the rather stable process behavior in multiple impeller systems due to the reduced hydrodynamic heterogeneity. Such comprehensive insights lead to a deeper understanding of the bioprocess.

Hydraulic Characteristics and Numerical Simulation of the Entrance Section of Ladder-Shaped Spillway

The ladder-shaped spillway in a certain reservoir junction is set as the engineering background in the paper. The hydraulic similarly model experiment and three-dimensional numerical simulation of hydraulic characteristics of water flow are performed. The outflow capacity, flow state analysis, velocity distribution, water surface line, pressure, and the energy dissipation rate are analyzed, and experimental results are compared with the numerical results. The conclusions demonstrate that the numerical results of the flow characteristics are very proximate to actual experimental results, the changeable law is the same, and their energy dissipation rate is basically consistent; it shows the feasibility of three-dimensional numerical simulation; the conclusions can provide the basis for the optimization about the flow state of the ladder-shaped spillway in the future.

Law of bounded dissipation and its consequences in turbulent wall flows

The dominant paradigm in turbulent wall flows is that the mean velocity near the wall, when scaled on wall variables, is independent of the friction Reynolds number $Re_\tau$ . This paradigm faces challenges when applied to fluctuations but has received serious attention only recently. Here, by extending our earlier work (Chen & Sreenivasan, J. Fluid Mech., vol. 908, 2021, p. R3) we present a promising perspective, and support it with data, that fluctuations displaying non-zero wall values, or near-wall peaks, are bounded for large values of $Re_\tau$ , owing to the natural constraint that the dissipation rate is bounded. Specifically, $\varPhi _\infty - \varPhi = C_\varPhi \,Re_\tau ^{-1/4},$ where $\varPhi$ represents the maximum value of any of the following quantities: energy dissipation rate, turbulent diffusion, fluctuations of pressure, streamwise and spanwise velocities, squares of vorticity components, and the wall values of pressure and shear stresses; the subscript $\infty$ denotes the bounded asymptotic value of $\varPhi$ , and the coefficient $C_\varPhi$ depends on $\varPhi$ but not on $Re_\tau$ . Moreover, there exists a scaling law for the maximum value in the wall-normal direction of high-order moments, of the form $\langle \varphi ^{2q}\rangle ^{{1}/{q}}_{max}= \alpha _q-\beta _q\,Re^{-1/4}_\tau$ , where $\varphi$ represents the streamwise or spanwise velocity fluctuation, and $\alpha _q$ and $\beta _q$ are independent of $Re_\tau$ . Excellent agreement with available data is observed. A stochastic process for which the random variable has the form just mentioned, referred to here as the ‘linear $q$ -norm Gaussian’, is proposed to explain the observed linear dependence of $\alpha _q$ on $q$ .

Energy dissipation performance of the trapezoidal stepped spillway

Stepped spillways are a more effective type of spillway in energy dissipation than conventional chute channels. Therefore, the dimensions of the energy breaker at the downstream of the stepped spillways are lower. It is an alternative especially for the downstream pool that cannot be built in sufficient length due to the terrain conditions. In this study, the energy dissipation performance of the trapezoidal stepped spillways was investigated numerically by using Flow3D software. Four different models and three different discharges were utilized for this aim. According to the results, the trapezoidal stepped spillway is more effective up to 30% than classical stepped spillways in energy dissipation. The depth of the trapezoidal step and the bottom base length of the trapezoid significantly affected the energy dissipation rate for the trapezoidal stepped spillway.

Determination of hydraulic resistance of channels using spectral geometry methods

In this paper, we propose a new method for calculation of hydraulic resistance of channels with constant cross-section. The method is based on the obtained estimates for the average energy dissipation rate in a turbulent flow. The first part of the paper is devoted to theoretical justification of the method. The second part is devoted to calculation of hydraulic resistance of various channels using the abovementioned method and comparison of these values with the known results. The proposed method allows for calculation of hydraulic resistance of various channels with sufficiently high accuracy and is based only on the information about the channel geometry.

Numerical evidence of anomalous energy dissipation in incompressible Euler flows: towards grid-converged results for the inviscid Taylor–Green problem

The well-known energy dissipation anomaly in the inviscid limit, related to velocity singularities according to Onsager, still needs to be demonstrated by numerical experiments. The present work contributes to this topic through high-resolution numerical simulations of the inviscid three-dimensional Taylor–Green vortex problem using a novel high-order discontinuous Galerkin discretisation approach for the incompressible Euler equations. The main methodological ingredient is the use of a discretisation scheme with inbuilt dissipation mechanisms, as opposed to discretely energy-conserving schemes, which – by construction – rule out the occurrence of anomalous dissipation. We investigate effective spatial resolution up to $8192^3$ (defined based on the $2{\rm \pi}$ -periodic box) and make the interesting phenomenological observation that the kinetic energy evolution does not tend towards exact energy conservation for increasing spatial resolution of the numerical scheme, but that the sequence of discrete solutions seemingly converges to a solution with non-zero kinetic energy dissipation rate. Taking the fine-resolution simulation as a reference, we measure grid-convergence with a relative $L^2$ -error of $0.27\,\%$ for the temporal evolution of the kinetic energy and $3.52\,\%$ for the kinetic energy dissipation rate against the dissipative fine-resolution simulation. The present work raises the question of whether such results can be seen as a numerical confirmation of the famous energy dissipation anomaly. Due to the relation between anomalous energy dissipation and the occurrence of singularities for the incompressible Euler equations according to Onsager's conjecture, we elaborate on an indirect approach for the identification of finite-time singularities that relies on energy arguments.

Percampuran vertikal di Perairan Laut Maluku dan Talaud pada bulan Februari 2021

<strong>Vertical mixing in the northern Maluku Sea and Talaud Waters in February 2021. </strong>The spatial variability of water mass mixing in the northern Maluku Sea and Talaud waters are presented based on the results of Eastern Indonesia Expedition (EIT) 2021 using RV Baruna Jaya VIII-LIPI. The turbulent kinetic energy dissipation rate was obtained using the Kunze-Williams-Briscoe (KWB) Method calculated from CTD (Conductivity, Temperature, Depth) and LADCP (Lowered Acoustic Doppler Current Profiler) datasets. We found the dissipation rate in the core layer of North Pacific Subtropical Water (NPSW) and North Pacific Intermediate Water (NPIW) are in the order of 10^-6 W/kg and 10^-8 W/kg, respectively. The KWB Method used in this study is also proven comparable with the Thorpe Method.

Statistical Analysis of Dynamic Subgrid Modeling Approaches in Large Eddy Simulation

In large eddy simulation (LES) of turbulent flows, dynamic subgrid models would account for an average cascade of kinetic energy from the largest to the smallest scales of the flow. Yet, it is unclear which of the most critical dynamical processes can ensure the criterion mentioned above. Furthermore, evidence of vortex stretching being the primary mechanism of the cascade is not out of the question. In this article, we study essential statistical characteristics of vortex stretching. Our numerical results demonstrate that vortex stretching rate provides the energy dissipation rate necessary for modeling subgrid-scale turbulence. We have compared the interaction of subgrid stresses with the filtered quantities among four models using invariants of the velocity gradient tensor. The individual and the joint probability of vortex stretching and strain amplification show that vortex stretching rate is highly correlated with the energy cascade rate. Sheet-like flow structures are correlated with viscous dissipation, and vortex tubes are more stretched than compressed. The overall results indicate that the stretching mechanism extracts energy from the large-scale straining motion and passes it onto small-scale stretched vortices.

Stochastic models for the droplet motion and evaporation in under-resolved turbulent flows at a large Reynolds number

In this work we introduce a Lagrangian stochastic model for particle motion and evaporation to be used in large-eddy simulations (LES) of turbulent liquid sprays. Effects of small-scale intermittency, usually under-resolved in LES, are explicitly included via modelling of the energy dissipation rate seen by a droplet along its trajectory. Namely, the dissipation rate is linked to the norm of the droplet sub-filtered acceleration which is included in the droplet motion equation. This norm, along with the direction of the droplet sub-filtered acceleration, is simulated as a stochastic process. With increasing Reynolds number, the distribution of the sub-filtered acceleration develops longer tails, with a slower decay in auto-correlation functions of the norm and direction of this acceleration. The stochastic models are specified for particles larger and smaller the Kolmogorov length scale. The assumption of the droplet evaporation model is similar, i.e. the evaporation rate is strongly enhanced when a droplet is subjected to very localized zones of intense velocity gradients. Thereby, the overall evaporation process is assumed to be a succession of two steady-state sub-processes with equal intensities, i.e. evaporation and vapour mixing. Then the stochastic properties of the overall evaporation rate are also controlled by fluctuations of the energy dissipation rate along the droplet path, and with increasing Reynolds number, the intensity of fluctuations of this rate is also increasing. The assessment of the presented stochastic models in LES of high-speed non-evaporating and evaporating sprays show the accurate prediction of experimental data on relatively coarser grids along with a remarkably weaker sensitivity to the grid spacing. The joint statistics and Voronoi tessellations exhibit strong intermittency of evaporation rate. The intensity of turbulence along the droplet pathway substantially promotes the vaporization rate.

On the applicability of Taylor's hypothesis, including small sampling velocities

Taylor's hypothesis, or the frozen turbulence approximation, can be used to estimate also the specific energy dissipation rate $\epsilon$ by comparing experimental results with the Kolmogorov–Obukhov expression. The hypothesis assumes that a frequency detected by an instrument moving with a constant large velocity $V$ can be related to a wavenumber by $\omega = k V$ . It is, however, not obvious how large the translational velocity has to be in order to make the hypothesis valid, or at least applicable with some acceptable uncertainty. Using the space–time-varying structure function for homogeneous and isotropic conditions, this question is addressed in the present study with emphasis on small velocities $V$ . The structure function is obtained using results from numerical solutions of the Navier–Stokes equation. Particular attention is given to the $V$ variation of the estimated specific energy dissipation, $\epsilon _{est}$ , compared with the actual value, $\epsilon$ , used in the numerical calculations. In contrast to previous studies, the results emphasize velocities $V$ less than or comparable to the one-component root-mean-square velocity, $u_{rms}$ . We find that $\epsilon$ can be determined to an acceptable accuracy for $V \geq 0.3\,u_{rms}$ . A simple analytical model is suggested to explain the main features of the observations, both Eulerian and Lagrangian. The model assumes that the observed time variations are solely due to eddies moving past the observer, thus ignoring eddy deformation and intermittency effects. In spite of these simplifications, the analysis accounts for most of the numerical results when also eddy-size-dependent velocities are accounted for.