scholarly journals Direct numerical simulation of turbulent mixing

Author(s):  
V. P. Statsenko ◽  
Yu. V. Yanilkin ◽  
V. A. Zhmaylo

The results of three-dimensional numerical simulations of turbulent flows obtained by various authors are reviewed. The paper considers the turbulent mixing (TM) process caused by the development of the main types of instabilities: those due to gravitation (with either a fixed or an alternating-sign acceleration), shift and shock waves. The problem of a buoyant jet is described as an example of the mixed-type problem. Comparison is made with experimental data on the TM zone width, profiles of density, velocity and turbulent energy and degree of homogeneity.

Author(s):  
Allen E. Badeau ◽  
Ismail B. Celik

The objective of this study is to investigate vertical buoyant jets in an enclosure using Large Eddy Simulation (LES) methods with no sub-grid scale model. This type of methodology is called Implicit Turbulent Modeling (ITM). Two different boundary conditions are applied at the inlet, being a uniform and periodic forcing velocity distribution. To accomplish this goal, a numerical solver was written, named DREAM®, which is capable of solving three dimensional, transient flows using an accurate monotonic upwinding scheme. The three-dimensional Navier-Stokes equations are solved in Cartesian coordinates, with the control volume approach being implemented on a staggered grid. The numerical scheme uses a fractional time step method, Crank-Nicolson, with the overall spatial and temporal accuracy being second order. In ITM simulations, there is no explicit subgrid-scale model (SGS) used for the modeling of the small scale vortical structures. ITM simulations assume that through strict conservation of the fluxing quantities in and out of the cell, the grid resolution is fully capable of capturing the important scales of the flow. The control volume averaging techniques used in the ITM methods acts as an implicit subgrid-scale model, and the resolvable scales of the flow are only dependent on the grid resolution within the domain. The available experimental data, as well as simulations that used SGS models, compare favorably to the ITM simulations from DREAM® in most cases as long as an “adequate” grid resolution is maintained. Results show that the density stratification tends to accelerate the jet and increase the amount of turbulence present within the flow. Perturbation of the inlet boundary condition ensures a sooner onset of turbulence, which is faster than the non-perturbed inlet boundary condition. A similarity solution is achieved at approximately 8 and 13 inlet diameters downstream of the inlet for the perturbed and uniform inlet boundary condition. Comparison between the vertical buoyant jet simulations to the available experimental data shows good agreement for the jet width and buoyant path centerline locations based on the internal densimetric Froude number. The application of these methods to immiscible fluids shows a new dimension to ITM and allows for a high resolution of the resulting flow field without the need for an explicit SGS model.


Kerntechnik ◽  
2021 ◽  
Vol 86 (3) ◽  
pp. 210-216
Author(s):  
M. P. Sharma ◽  
A. Moharana

Abstract Subchannel analysis codes are widely used for the thermal-hydraulic design of nuclear reactor rod bundle. The effectiveness of subchannel analysis codes depends on turbulent mixing between these subchannels. Turbulent mixing has no direct contribution to the axial mass flow rate through subchannel but it will cause exchange of momentum and energy between the neighboring subchannels. Thus, it is important to evaluate the turbulent mixing coefficient for reactor rod bundle as it is a significant factor in the lateral energy and momentum equation for subchannel analysis codes like COBRA IIIC, COBRA-IV and MATRA LMR-FB. With the rapid developments in computational fluid dynamics and computer performance, three-dimensional analyses of turbulent flows occurring in the nuclear rod bundle have become more prominent. Several numerical analyses have already been attempted to investigate the flow behavior in rod bundles of different reactors. Much of these are dedicated to find out the structure of turbulence in rod bundle but a few analyses has been done to evaluate the magnitude of the turbulent mixing coefficient. In view of this, CFD analyses were carried out to determine the turbulent mixing coefficient in the simulated sub-channels of the reactor rod bundle. Previous studies on the structure of turbulence reveals that it is highly anisotropic. Hence, the Reynolds Stress Model (RSM), finer mesh and near wall distance ( y + ≤ 2) is required to capture turbulent mixing phenomena. The validation of results is done by comparing with subchannel mixing experiments.


2016 ◽  
Vol 804 ◽  
pp. 130-161 ◽  
Author(s):  
M. MacDonald ◽  
L. Chan ◽  
D. Chung ◽  
N. Hutchins ◽  
A. Ooi

We investigate rough-wall turbulent flows through direct numerical simulations of flow over three-dimensional transitionally rough sinusoidal surfaces. The roughness Reynolds number is fixed at $k^{+}=10$, where $k$ is the sinusoidal semi-amplitude, and the sinusoidal wavelength is varied, resulting in the roughness solidity $\unicode[STIX]{x1D6EC}$ (frontal area divided by plan area) ranging from 0.05 to 0.54. The high cost of resolving both the flow around the dense roughness elements and the bulk flow is circumvented by the use of the minimal-span channel technique, recently demonstrated by Chung et al. (J. Fluid Mech., vol. 773, 2015, pp. 418–431) to accurately determine the Hama roughness function, $\unicode[STIX]{x0394}U^{+}$. Good agreement of the second-order statistics in the near-wall roughness-affected region between minimal- and full-span rough-wall channels is observed. In the sparse regime of roughness ($\unicode[STIX]{x1D6EC}\lesssim 0.15$) the roughness function increases with increasing solidity, while in the dense regime the roughness function decreases with increasing solidity. It was found that the dense regime begins when $\unicode[STIX]{x1D6EC}\gtrsim 0.15{-}0.18$, in agreement with the literature. A model is proposed for the limit of $\unicode[STIX]{x1D6EC}\rightarrow \infty$, which is a smooth wall located at the crest of the roughness elements. This model assists with interpreting the asymptotic behaviour of the roughness, and the rough-wall data presented in this paper show that the near-wall flow is tending towards this modelled limit. The peak streamwise turbulence intensity, which is associated with the turbulent near-wall cycle, is seen to move further away from the wall with increasing solidity. In the sparse regime, increasing $\unicode[STIX]{x1D6EC}$ reduces the streamwise turbulent energy associated with the near-wall cycle, while increasing $\unicode[STIX]{x1D6EC}$ in the dense regime increases turbulent energy. An analysis of the difference of the integrated mean momentum balance between smooth- and rough-wall flows reveals that the roughness function decreases in the dense regime due to a reduction in the Reynolds shear stress. This is predominantly due to the near-wall cycle being pushed away from the roughness elements, which leads to a reduction in turbulent energy in the region previously occupied by these events.


Author(s):  
Rudolf A. Izmaylov ◽  
Henry D. Lopulalan ◽  
Gertruida S. Norimarna

Based on experimental investigation of unsteady flow phenomena in centrifugal compressor, conducted in LPI (SPbSPU), including the measurements of high frequency pressure and hot-wire probes, numerical calculation was conducted with the aid of commercial ANSYS CFX12. The object was a model of industrial compressor, tested in full range of mass flows (from maximum to the onset of surge). Unsteady turbulent flows were calculated with URANS, turbulence model - SST-Menter. Region between the impeller and diffuser was calculated with Transient Method option. Time discretization Δ = 12 μsec. Sequential 20 revolutions of rotors were calculated, also LES model were used. Unsteady calculation results demonstrated reasonable agreement with experimental data both in meridional, and radial sections concerning the character of “jet – wakes” propagation, distribution of velocities in diffuser, and clearly showed three-dimensional flow conditions and reverse zones on the diffuser walls. “Frozen” velocity distributions during rotating stall demonstrated possibility of unsteady flow prediction with the aid of ANSYS CFX12. CFD results comparison with experimental data showed reasonable agreement for distributions of total pressures in diffuser, as well total and static characteristics of impeller and diffuser.


Water ◽  
2021 ◽  
Vol 13 (8) ◽  
pp. 1079
Author(s):  
Lena Mahl ◽  
Patrick Heneka ◽  
Martin Henning ◽  
Roman B. Weichert

The efficiency of a fishway is determined by the ability of immigrating fish to follow its attraction flow (i.e., its jet) to locate and enter the fishway entrance. The hydraulic characteristics of fishway entrance jets can be simplified using findings from widely investigated surface jets produced by shaped nozzles. However, the effect of the different boundary conditions of fishway entrance jets (characterized by vertical entrance slots) compared to nozzle jets must be considered. We investigate the downstream propagation of attraction jets from the vertical slot of a fishway entrance into a quiescent tailrace, considering the following boundary conditions not considered for nozzle jets: (1) slot geometry, (2) turbulence characteristics of the approach flow to the slot, and (3) presence of a lateral wall downstream of the slot. We quantify the effect of these boundary conditions using three-dimensional hydrodynamic-numeric flow simulations with DES and RANS turbulence models and a volume-of-fluid method (VoF) to simulate the free water surface. In addition, we compare jet propagation with existing analytical methods for describing jet propagations from nozzles. We show that a turbulent and inhomogeneous approach flow towards a vertical slot reduces the propagation length of the slot jet in the tailrace due to increased lateral spreading compared to that of a jet produced by a shaped nozzle. An additional lateral wall in the tailrace reduces lateral spreading and significantly increases the propagation length. For highly turbulent flows at fishway entrances, the RANS model tends to overestimate the jet propagation compared to the transient DES model.


Author(s):  
Rahid Zaman ◽  
Yujiang Xiang ◽  
Jazmin Cruz ◽  
James Yang

In this study, the three-dimensional (3D) asymmetric maximum weight lifting is predicted using an inverse-dynamics-based optimization method considering dynamic joint torque limits. The dynamic joint torque limits are functions of joint angles and angular velocities, and imposed on the hip, knee, ankle, wrist, elbow, shoulder, and lumbar spine joints. The 3D model has 40 degrees of freedom (DOFs) including 34 physical revolute joints and 6 global joints. A multi-objective optimization (MOO) problem is solved by simultaneously maximizing box weight and minimizing the sum of joint torque squares. A total of 12 male subjects were recruited to conduct maximum weight box lifting using squat-lifting strategy. Finally, the predicted lifting motion, ground reaction forces, and maximum lifting weight are validated with the experimental data. The prediction results agree well with the experimental data and the model’s predictive capability is demonstrated. This is the first study that uses MOO to predict maximum lifting weight and 3D asymmetric lifting motion while considering dynamic joint torque limits. The proposed method has the potential to prevent individuals’ risk of injury for lifting.


2020 ◽  
Vol 21 (20) ◽  
pp. 7702 ◽  
Author(s):  
Sofya I. Scherbinina ◽  
Philip V. Toukach

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.


Author(s):  
Stephan Uhkoetter ◽  
Stefan aus der Wiesche ◽  
Michael Kursch ◽  
Christian Beck

The traditional method for hydrodynamic journal bearing analysis usually applies the lubrication theory based on the Reynolds equation and suitable empirical modifications to cover turbulence, heat transfer, and cavitation. In cases of complex bearing geometries for steam and heavy-duty gas turbines this approach has its obvious restrictions in regard to detail flow recirculation, mixing, mass balance, and filling level phenomena. These limitations could be circumvented by applying a computational fluid dynamics (CFD) approach resting closer to the fundamental physical laws. The present contribution reports about the state of the art of such a fully three-dimensional multiphase-flow CFD approach including cavitation and air entrainment for high-speed turbo-machinery journal bearings. It has been developed and validated using experimental data. Due to the high ambient shear rates in bearings, the multiphase-flow model for journal bearings requires substantial modifications in comparison to common two-phase flow simulations. Based on experimental data, it is found, that particular cavitation phenomena are essential for the understanding of steam and heavy-duty type gas turbine journal bearings.


1984 ◽  
Vol 62 (3) ◽  
pp. 596-600 ◽  
Author(s):  
R. G. Barradas ◽  
D. S. Nadezhdin

The cathodic reduction of the lead monoxide layer formed on lead in 30% aqueous H2SO4 under anodic oxidation at 0.6 V (vs. Hg/HgSO4 reference electrode) was investigated by linear sweep voltammetry, potential step and admittance measurements. The experimental data were analyzed respectively in terms of thin-layer electrochemistry, electrocrystallisation, and changes of resistance of the PbO layer under reduction. The results seem to be best interpreted from the theory of three-dimensional electrocrystallisation as PbO is reduced to Pb. At sub-zero temperatures the PbO peak observed on our voltammograms and potentiostatic current time transients reveals the splitting of the curves into two peaks, which may be a result of reduction of the same material but of different phases, namely, orthorhombic and tetragonal PbO.


1962 ◽  
Vol 99 (6) ◽  
pp. 558-569 ◽  
Author(s):  
Peter J. Wyllie

AbstractBowen's petrogenetic grid is a PT projection containing univariant curves for decarbonation, dehydration, and solid-solid reactions, with vapour pressure (Pf) equal to total pressure (Ps). Analysis of experimental data in the system MgO–CO2–H2O leads to an expansion of this grid. Three of the important variables in metamorphism when Pf = Ps are P, T, and variation of the pore fluid composition between H2O and CO2. These can be illustrated in a three-dimensional petrogenetic model; one face is a PT plane for reactions occurring with pure H2O, and the opposite face is a similar plane for reactions with pure CO2; these are separated by an axis for pore fluid composition varying between H2O and CO2. Superposition of the PT faces of the model provides the petrogenetic grid. The reactions within the model are represented by divariant surfaces, which may meet along univariant lines. For dissociation reactions, the surfaces curve towards lower temperatures as the proportion of non-reacting volatile increases, and solid-solid reaction surfaces are parallel to the vapour composition axis and perpendicular to the PT axes. The relative temperatures of reactions and the lines of intersections of the surfaces can be illustrated in isobaric sections. Isobaric sections are used to illustrate reactions proceeding at constant pressure with (1) pore fluid composition remaining constant during the reaction, with temperature increasing (2) pore fluid composition changing during the reaction, with temperature increasing, and (3) pore fluid changing composition at constant temperature. The petrogenetic model provides a convenient framework for a wide range of experimental data.


Sign in / Sign up

Export Citation Format

Share Document