Applicability of Dynamic Modeling for Turbulent Schmidt and Prandtl Numbers on Density Stratification Breakup in Several Flow Conditions

Author(s):  
Satoshi Abe ◽  
Etienne Studer ◽  
Masahiro Ishigaki ◽  
Yasuteru Sibamoto ◽  
Taisuke Yonomoto

Abstract Many experiments on density stratification breakup in several flow conditions have been performed with the large- and small-scale experimental facilities to understand the mechanism underlying hydrogen behavior in a nuclear containment vessel during a severe accident. To improve the predictability of the RANS (Reynolds-averaged Navier Stokes) approach, we implemented the dynamic modeling for turbulent Schmidt Sct and Prandtl Prt numbers. In this paper, the capability of the RANS analysis with dynamic Sct modeling is assessed with several experimental data obtained by using the MISTRA (Commissariat à l’énergie atomique et aux énergies alternatives, CEA, France), CIGMA and VIMES (Japan Atomic Energy Agency, Japan). For the quantitative assessment, the completion time of the stratification breakup, defined as when helium concentration in the upper region decreases to the same value in the lower region, is focused. The comparison study shows the good performance of the dynamic modeling for Sct and Prt. Besides, in the case with the low jet Froude number, the CFD accuracy declines significantly, because the jet upward bending is over-estimated.

1987 ◽  
Vol 174 ◽  
pp. 209-231 ◽  
Author(s):  
H. Gao ◽  
G. Metcalfe ◽  
T. Jung ◽  
R. P. Behringer

This paper first describes an apparatus for measuring the Nusselt number N versus the Rayleigh number R of convecting normal liquid 4He layers. The most important feature of the apparatus is its ability to provide layers of different heights d, and hence different aspect ratios [Gcy ]. The horizontal cross-section of each layer is circular, and [Gcy ] is defined by [Gcy ] = D/2d where D is the diameter of the layer. We report results for 2.4 [les ] [Gcy ] [les ] 16 and for Prandtl numbers Pr spanning 0.5 [lsim ] Pr [lsim ] 0.9 These results are presented in terms of the slope N1 = RcdN/dR evaluated just above the onset of convection at Rc. We find that N1 is only a slowly increasing function of [Gcy ] in the range 6 [lsim ] [Gcy ] [lsim ] 16, and that it has a value there which is quite close to 0.72. This value of N1 is in good agreement with variational calcuations by Ahlers et al. (1981) pertinent to parallel convection rolls in cylindrical geometry. Particularly for [Gcy ] [lsim ] 6, we find additional small-scale structure in N1 associated with changes in the number of convection rolls with changing [Gcy ]. An additional test of the linearzied hydrodynamics is given by measurements of Rc. We find good agreement between theory and our data for Rc.


1994 ◽  
Vol 47 (6S) ◽  
pp. S3-S13 ◽  
Author(s):  
Parviz Moin ◽  
Thomas Bewley

A brief review of current approaches to active feedback control of the fluctuations arising in turbulent flows is presented, emphasizing the mathematical techniques involved. Active feedback control schemes are categorized and compared by examining the extent to which they are based on the governing flow equations. These schemes are broken down into the following categories: adaptive schemes, schemes based on heuristic physical arguments, schemes based on a dynamical systems approach, and schemes based on optimal control theory applied directly to the Navier-Stokes equations. Recent advances in methods of implementing small scale flow control ideas are also reviewed.


Author(s):  
Nur Alom ◽  
Ujjwal K. Saha

The Savonius rotor appears to be particularly promising for the small-scale applications because of its design simplicity, good starting ability, and insensitivity to wind directions. There has been a growing interest in recent times to harness wind energy in an efficient manner by developing newer blade profiles of Savonius rotor. The overlap ratio (OR), one of the important geometric parameters, plays a crucial role in the turbine performance. In a recent study, an elliptical blade profile with a sectional cut angle (θ) of 47.5° has demonstrated its superior performance when set at an OR = 0.20. However, this value of OR is ideal for a semicircular profile, and therefore, requires further investigation to arrive at the optimum overlap ratio for the elliptical profile. In view of this, the present study attempts to make a systemic numerical study to arrive at the optimum OR of the elliptical profile having sectional cut angle, θ = 47.5°. The 2D unsteady simulation is carried out around the elliptical profile considering various overlap ratios in the range of 0.0 to 0.30. The continuity, unsteady Reynolds Averaged Navier-Stokes (URANS) equations and two equation eddy viscosity SST (Shear Stress transport) k-ω model are solved by using the commercial finite volume method (FVM) based solver ANSYS Fluent. The torque and power coefficients are calculated as a function of tip speed ratio (TSR) and at rotating conditions. The total pressure, velocity magnitude and turbulence intensity contours are obtained and analyzed to arrive at the intended objective. The numerical simulation demonstrates an improved performance of the elliptical profile at an OR = 0.15.


Author(s):  
Jun Ishikawa ◽  
Tomoyuki Sugiyama ◽  
Yu Maruyama

The Japan Atomic Energy Agency (JAEA) is pursuing the development and application of the methodologies on fission product (FP) chemistry for source term analysis by using the integrated severe accident analysis code THALES2. In the present study, models for the eutectic interaction of boron carbide (B4C) with steel and the B4C oxidation were incorporated into THALES2 code and applied to the source term analyses for a boiling water reactor (BWR) with Mark-I containment vessel (CV). Two severe accident sequences with drywell (D/W) failure by overpressure initiated by loss of core coolant injection (TQUV sequence) and long-term station blackout (TB sequence) were selected as representative sequences. The analyses indicated that a much larger amount of species from the B4C oxidation was produced in TB sequence than TQUV sequence. More than a half of carbon dioxide (CO2) produced by the B4C oxidation was predicted to dissolve into the water pool of the suppression chamber (S/C), which could largely influence pH of the water pool and consequent formation and release of volatile iodine species.


GeoArabia ◽  
2003 ◽  
Vol 8 (1) ◽  
pp. 47-86 ◽  
Author(s):  
Jürgen Grötsch ◽  
Omar Suwaina ◽  
Ghiath Ajlani ◽  
Ahmed Taher ◽  
Reyad El-Khassawneh ◽  
...  

ABSTRACT A 3-D geological model of the Kimmeridgian-Tithonian Manifa, Hith, Arab, and Upper Diyab formations in the area of the onshore Central Abu Dhabi Ridge was based on a newly established sequence stratigraphic, sedimentologic, and diagenetic model. It was part of an inter-disciplinary study of the large sour-gas reserves in Abu Dhabi that are mainly hosted by the Arab Formation. The model was used for dynamic evaluations and recommendations for further appraisal and development planning in the studied field. Fourth-order aggradational and progradational cycles are composed of small-scale fifth-order shallowing-upward cycles, mostly capped by anhydrite within the Arab-ABC. The study area is characterized by a shoreline progradation of the Arab Formation toward the east-northeast marked by high-energy oolitic/bioclastic grainstones of the Upper Arab-D and the Asab Oolite. The Arab-ABC, Hith, and Manifa pinch out toward the northeast. The strongly bioturbated Lower Arab-D is an intrashelf basinal carbonate ramp deposit, largely time-equivalent to the Arab-ABC. The deposition of the Manifa Formation over the Arab Formation was a major back-stepping event of the shallow-water platform before the onset of renewed progradation in the Early Cretaceous. Well productivity in the Arab-ABC is controlled mainly by thin, permeable dolomitic streaks in the fifth-order cycles at the base of the fourth-order cycles. This has major implications for reservoir management, well completion and stimulation, and development planning. Good reservoir properties have been preserved in the early diagenetic dolomitic streaks. In contrast, the reservoir properties of the Upper Arab-D oolitic/bioclastic grainstones deteriorate with depth due to burial diagenesis. A rock-type scheme was established because complex diagenetic overprinting prevented the depositional facies from being directly related to petrophysical properties. Special core analysis and the attribution of saturation functions to static and dynamic models were made on a cell-by-cell basis using the scheme and honoring the 3-D depositional facies and property model. The results demonstrated the importance of integrating sedimentological analysis and diagenesis with rock typing and static and dynamic modeling so as to enhance the predictive capabilities of subsurface models.


1984 ◽  
Vol 106 (3) ◽  
pp. 552-557 ◽  
Author(s):  
K. V. Rao ◽  
B. F. Armaly ◽  
T. S. Chen

Laminar mixed forced and free convection from a line thermal source imbedded at the leading edge of an adiabatic vertical surface is analytically investigated for the cases of buoyancy assisting and buoyancy opposing flow conditions. Temperature and velocity distributions in the boundary layer adjacent to the adiabatic surface are presented for the entire range of the buoyancy parameter ξ (x) = Grx/Rex5/2 from the pure forced (ξ(x) = 0) to the pure free (ξ(x) = ∞) convection regime for fluids having Prandtl numbers of 0.7 and 7.0. For buoyancy-assisting flow, the velocity overshoot, the temperature, and the wall shear stress increase as the plume’s strength increases. On the other hand, the velocity overshoot, the wall shear stress, and the temperature decrease as the free-stream velocity increases. For buoyancy opposing flow, the velocity and wall shear stress decrease but the temperature increases as the plume’s strength increases.


2018 ◽  
Vol 856 ◽  
Author(s):  
M. Borgnino ◽  
G. Boffetta ◽  
F. De Lillo ◽  
M. Cencini

We study the dynamics and the statistics of dilute suspensions of gyrotactic swimmers, a model for many aquatic motile microorganisms. By means of extensive numerical simulations of the Navier–Stokes equations at different Reynolds numbers, we investigate preferential sampling and small-scale clustering as a function of the swimming (stability and speed) and shape parameters, considering in particular the limits of spherical and rod-like particles. While spherical swimmers preferentially sample local downwelling flow, for elongated swimmers we observe a transition from downwelling to upwelling regions at sufficiently high swimming speed. The spatial distribution of both spherical and elongated swimmers is found to be fractal at small scales in a wide range of swimming parameters. The direct comparison between the different shapes shows that spherical swimmers are more clusterized at small stability and speed numbers, while for large values of the parameters elongated cells concentrate more. The relevance of our results for phytoplankton swimming in the ocean is briefly discussed.


2018 ◽  
Vol 140 (6) ◽  
Author(s):  
Hamed Avari ◽  
Kem A. Rogers ◽  
Eric Savory

The parallel plate flow chamber (PPFC) has gained popularity due to its applications in fields such as biological tissue engineering. However, most of the studies using PPFC refer to theoretical relations for estimating the wall shear stress (WSS) and, hence, the accuracy of such quantifications remains elusive for anything other than steady laminar flow. In the current study, a laser Doppler velocimetry (LDV) method was used to quantify the flow in a PPFC (H = 1.8 mm × W = 17.5 mm, Dh = 3.26 mm, aspect ratio = 9.72) under steady Re = 990, laminar pulsatile (carotid Re0-mean = 282 as well as a non-zero-mean sinusoidal Re0-mean = 45 pulse) and low-Re turbulent Re = 2750 flow conditions. A mini-LDV probe was applied, and the absolute location of the LDV measuring volume with the respect to the wall was determined using a signal monitoring technique with uncertainties being around ±27 μm. The uniformity of the flow across the span of the channel, as well as the WSS assessment for all the flow conditions, was measured with the uncertainties all being less than 16%. At least two points within the viscous sublayer of the low-Re turbulent flow were measured (with the y+ for the first point < 3) and the WSS was determined using two methods with the differences between the two methods being within 5%. This paper for the first time presents the experimental determination of WSS using LDV in a small-scale PPFC under various flow conditions, the challenges associated with each condition, and a comparison between the cases. The present data will be useful for those conducting biological or numerical modeling studies using such devices.


2021 ◽  
Author(s):  
Michael Kurgansky

&lt;p&gt;The linear and non-linear inertial stability of the Kolmogorov flow in a rotating viscous fluid of uniform density is investigated. A necessary condition for instability is the violation of the criterion of non-viscous inertial stability, and the sufficient condition of instability is formulated in terms of the Reynolds criterion. The existence of stable secondary stationary regimes in the problem is shown, developing in a context of loss of stability of the main flow and having the shape of rolls (cloud streets in the atmosphere) oriented along it. Stable density stratification is taken into account when the direction of gravity coincides with the direction of rotation of the fluid. In this case, the necessary condition for the inertial instability of the main flow remains the same, but the critical Reynolds number for the instability depends on two additional dimensionless parameters that appear in the problem: the stratification parameter and the Prandtl number. The case of Prandtl numbers less than or equal to unity has been studied in greater detail, when there is a secondary stationary regime, which can be unstable - in contrast to the case of a fluid that is uniform in density - and density stratification is a destabilizing factor.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document