scholarly journals Numerical Simulation of Thermal-Solutal Capillary-Buoyancy Flow of Ge1–xSix Single Crystals Driven by Surface-Tension and Rotation in a Czochralski Configuration

Crystals ◽  
2019 ◽  
Vol 9 (4) ◽  
pp. 217 ◽  
Author(s):  
Jia-Jia Yu ◽  
Lu Zhang ◽  
Ting Shen ◽  
Li Zhang ◽  
You-Rong Li
Keyword(s):  
Reynolds Number ◽  
Traveling Waves ◽  
Rotation Rate ◽  
Three Dimensional ◽  
Maximum Amplitude ◽  
Reynolds Numbers ◽  
Solute Concentration ◽  
Buoyancy Flow ◽  
Constant Rate ◽  
Crucible Rotation

A series of three-dimensional numerical simulations were performed to understand the thermal-solutal capillary-buoyancy flow of Ge1-xSix melts during Czochralski crystal growth with a rotating crystal or crucible. The crystal and crucible rotation Reynolds numbers in this work are 0∼3.5 × 103 (0∼4.4 rpm) and 0∼−2.4 × 103 (0∼−1.5 rpm), respectively. Simulation results show that if the thermal capillary Reynolds number is relatively low, the flow will be steady and axisymmetric, even though the crystal or crucible rotates at a constant rate. The critical thermal capillary Reynolds number for the initiation of the three-dimensional oscillatory flow is larger than that of pure fluids. As the crystal or crucible rotation rate increases, the critical thermal capillary Reynolds number first increases and then decreases. The dominant flow pattern after the flow destabilization is azimuthal traveling waves. Furthermore, a reversed evolution from the oscillatory spoke pattern to traveling waves appears in the melt. Once the crystal or crucible rotation rate is relatively large, the traveling waves respectively evolve to rotating waves at the crystal rotation and a spindle-like pattern at the crucible rotation. In addition, the maximum amplitude of solute concentration oscillation on the free surface initially decreases, but finally rises with the crystal or crucible rotation rate increasing.

Numerical investigation of wake flow regimes behind a high-speed rotating circular cylinder in steady flow

2019 ◽  
Vol 878 ◽  
pp. 875-906
Author(s):  
Adnan Munir ◽  
Ming Zhao ◽  
Helen Wu ◽  
Lin Lu
Keyword(s):  
Reynolds Number ◽  
Circular Cylinder ◽  
Vortex Shedding ◽  
Rotation Rate ◽  
High Speed ◽  
Three Dimensional ◽  
Flow Regimes ◽  
Reynolds Numbers ◽  
Wake Flow ◽  
Rotating Circular Cylinder

Flow around a high-speed rotating circular cylinder for $Re\leqslant 500$ is investigated numerically. The Reynolds number is defined as $UD/\unicode[STIX]{x1D708}$ with $U$, $D$ and $\unicode[STIX]{x1D708}$ being the free-stream flow velocity, the diameter of the cylinder and the kinematic viscosity of the fluid, respectively. The aim of this study is to investigate the effect of a high rotation rate on the wake flow for a range of Reynolds numbers. Simulations are performed for Reynolds numbers of 100, 150, 200, 250 and 500 and a wide range of rotation rates from 1.6 to 6 with an increment of 0.2. Rotation rate is the ratio of the rotational speed of the cylinder surface to the incoming fluid velocity. A systematic study is performed to investigate the effect of rotation rate on the flow transition to different flow regimes. It is found that there is a transition from a two-dimensional vortex shedding mode to no vortex shedding mode when the rotation rate is increased beyond a critical value for Reynolds numbers between 100 and 200. Further increase in rotation rate results in a transition to three-dimensional flow which is characterized by the presence of finger-shaped (FV) vortices that elongate in the wake of the cylinder and very weak ring-shaped vortices (RV) that wrap the surface of the cylinder. The no vortex shedding mode is not observed at Reynolds numbers greater than or equal to 250 since the flow remains three-dimensional. As the rotation rate is increased further, the occurrence frequency and size of the ring-shaped vortices increases and the flow is dominated by RVs. The RVs become bigger in size and the flow becomes chaotic with increasing rotation rate. A detailed analysis of the flow structures shows that the vortices always exist in pairs and the strength of separated shear layers increases with the increase of rotation rate. A map of flow regimes on a plane of Reynolds number and rotation rate is presented.

scholarly journals Flow Instabilities of Coupled Rotation and Thermal-Solutal Capillary Convection of Binary Mixture in Czochralski Configuration

Crystals ◽  
2019 ◽  
Vol 9 (2) ◽  
pp. 72 ◽  
Author(s):  
Chunmei Wu ◽  
Bo Yuan ◽  
Yourong Li
Keyword(s):  
Reynolds Number ◽  
Binary Mixture ◽  
Rotation Rate ◽  
Oscillatory Flow ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Basic Flow ◽  
Flow Instabilities ◽  
Flow Transitions ◽  
Rotational Waves

In order to understand the flow instabilities of coupled rotation and thermal-solutal capillary convection of binary mixture in a Czochralski configuration subjected to simultaneous radial thermal and solutal gradients, a series of three-dimensional direct numerical simulation have been conducted. The capillary ratio of the silicon-germanium mixture is −0.2. The rotation Reynolds numbers of crystal and crucible, Res and Rec range from 0 to 3506 and 0 to 1403, respectively. Results show that the basic flow is axisymmetric and steady. It has rich flow structures in the meridian plane, depending on the competitions among the driving forces. With the increase of thermocapillary and rotation Reynolds numbers, the basic flow will transit to three dimensional oscillatory flow. For different combination of rotation rate and thermocapillary Reynolds number, the oscillatory flow can be displayed as spoke patterns which is steady in time but oscillate in space, spoke patterns propagate in azimuthal direction, rotational waves or coexistence of spokes and rotational waves. The crucible rotation has an inhibitory effect on the flow instability, inducing the monotonically increase of critical value for flow transitions, however, for crystal rotation, the critical thermocapillary Reynolds number increases at first and then decreases. When the rotation rate is large, two flow transitions are captured.

Bifurcation Characteristics of Flows in Rectangular Sudden Expansion Channels

2005 ◽  
Author(s):  
Francine Battaglia ◽  
George Papadopoulos
Keyword(s):  
Reynolds Number ◽  
Laminar Flow ◽  
Critical Reynolds Number ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Expansion Ratio ◽  
Sudden Expansion ◽  
Two Dimensional ◽  
Effective Expansion ◽  
Three Dimensional Simulations

The effect of three-dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio of in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations fail to capture completely the total expansion effect on the flow, which couples both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When the two-dimensional geometry was consistent with the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations was compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, helping to sustain laminar flow symmetry to higher Reynolds numbers in comparison to nominally two-dimensional double-expansion geometries. Lastly, and most importantly, when the logarithm of the critical Reynolds number from all these studies was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.

scholarly journals Nonlinear amplification in hydrodynamic turbulence

2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung
Keyword(s):  
Reynolds Number ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Developed Turbulence ◽  
Advection Term ◽  
Nonlinear Amplification

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.

Experimental characterization of three-dimensional corner flows at low Reynolds numbers

2012 ◽  
Vol 707 ◽  
pp. 37-52 ◽  
Author(s):  
J. Sznitman ◽  
L. Guglielmini ◽  
D. Clifton ◽  
D. Scobee ◽  
H. A. Stone ◽  
...  
Keyword(s):  
Reynolds Number ◽  
Numerical Solutions ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Secondary Flows ◽  
Channel Cross Section ◽  
Experimental Characterization ◽  
Low Reynolds Numbers ◽  
Low Reynolds

AbstractWe investigate experimentally the characteristics of the flow field that develops at low Reynolds numbers ($\mathit{Re}\ll 1$) around a sharp $9{0}^{\ensuremath{\circ} } $ corner bounded by channel walls. Two-dimensional planar velocity fields are obtained using particle image velocimetry (PIV) conducted in a towing tank filled with a silicone oil of high viscosity. We find that, in the vicinity of the corner, the steady-state flow patterns bear the signature of a three-dimensional secondary flow, characterized by counter-rotating pairs of streamwise vortical structures and identified by the presence of non-vanishing transverse velocities (${u}_{z} $). These results are compared to numerical solutions of the incompressible flow as well as to predictions obtained, for a similar geometry, from an asymptotic expansion solution (Guglielmini et al., J. Fluid Mech., vol. 668, 2011, pp. 33–57). Furthermore, we discuss the influence of both Reynolds number and aspect ratio of the channel cross-section on the resulting secondary flows. This work represents, to the best of our knowledge, the first experimental characterization of the three-dimensional flow features arising in a pressure-driven flow near a corner at low Reynolds number.

Secondary instability and subcritical transition of the leading-edge boundary layer

2016 ◽  
Vol 792 ◽  
pp. 682-711 ◽  
Author(s):  
Michael O. John ◽  
Dominik Obrist ◽  
Leonhard Kleiser
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
High Speed ◽  
Three Dimensional ◽  
Leading Edge ◽  
Reynolds Numbers ◽  
Secondary Instability ◽  
Bypass Transition ◽  
Wall Suction ◽  
Transition Mechanism

The leading-edge boundary layer (LEBL) in the front part of swept airplane wings is prone to three-dimensional subcritical instability, which may lead to bypass transition. The resulting increase of airplane drag and fuel consumption implies a negative environmental impact. In the present paper, we present a temporal biglobal secondary stability analysis (SSA) and direct numerical simulations (DNS) of this flow to investigate a subcritical transition mechanism. The LEBL is modelled by the swept Hiemenz boundary layer (SHBL), with and without wall suction. We introduce a pair of steady, counter-rotating, streamwise vortices next to the attachment line as a generic primary disturbance. This generates a high-speed streak, which evolves slowly in the streamwise direction. The SSA predicts that this flow is unstable to secondary, time-dependent perturbations. We report the upper branch of the secondary neutral curve and describe numerous eigenmodes located inside the shear layers surrounding the primary high-speed streak and the vortices. We find secondary flow instability at Reynolds numbers as low as$Re\approx 175$, i.e. far below the linear critical Reynolds number$Re_{crit}\approx 583$of the SHBL. This secondary modal instability is confirmed by our three-dimensional DNS. Furthermore, these simulations show that the modes may grow until nonlinear processes lead to breakdown to turbulent flow for Reynolds numbers above$Re_{tr}\approx 250$. The three-dimensional mode shapes, growth rates, and the frequency dependence of the secondary eigenmodes found by SSA and the DNS results are in close agreement with each other. The transition Reynolds number$Re_{tr}\approx 250$at zero suction and its increase with wall suction closely coincide with experimental and numerical results from the literature. We conclude that the secondary instability and the transition scenario presented in this paper may serve as a possible explanation for the well-known subcritical transition observed in the leading-edge boundary layer.

Numerical Investigation on Turbulence Statistics and Heat Transfer of a Circular Jet Impinging on a Roughened Flat Plate

2021 ◽  
Vol 13 (4) ◽  
Author(s):  
Abdulrahman Alenezi ◽  
Abdulrahman Almutairi ◽  
Hamad Alhajeri ◽  
Abdulaziz Gamil ◽  
Faisal Alshammari
Keyword(s):  
Heat Transfer ◽  
Reynolds Number ◽  
Numerical Study ◽  
Roughness Element ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Target Plate ◽  
Cross Sectional ◽  
Flow Physics ◽  
Baseline Case

Abstract A detailed heat transfer numerical study of a three-dimensional impinging jet on a roughened isothermal surface is presented and is investigated from flow physics vantage point under the influence of different parameters. The effects of the Reynolds number, roughness location, and roughness dimension on the flow physics and heat transfer parameters are studied. Additionally, the relations between average heat transfer coefficient (AHTC) and flow physics including pressure, wall shear and flow vortices with thermodynamic nonequilibrium are offered. This paper studies the effect of varying both location and dimension of the roughness element which took the shape of square cross-sectional continuous ribs to deliver a favorable trade-off between total pressure loss and heat transfer rate. The roughness element was tested for three different radial locations (R/D) = 1, 1.5, and 2 and at each location its height (i.e., width) (e) was changed from 0.25 to 1 mm in incremental steps of 0.25. The study used a jet angle (α) of 90 deg, jet-to-target distance (H/D = 6), and Re ranges from 10,000 to 50,000, where H is the vertical distance between the target plate and jet exit. The results show that the AHTC can be significantly affected by changing the geometry and dimensions of the roughness element. This variation can be either an augmentation of, or decrease in, the (HTC) when compared with the baseline case. An enhancement of 12.9% in the AHTC was achieved by using optimal location and dimensions of the roughness element at specific Reynolds number. However, a diminution between 10% and 30% in (AHTC) was attained by the use of rib height e = 1 mm at Re = 50k. The variation of both rib location and height showed better contribution in increasing heat transfer for low-range Reynolds numbers.

Three-dimensional instabilities in the boundary-layer flow over a long rectangular plate

2011 ◽  
Vol 681 ◽  
pp. 411-433 ◽  
Author(s):  
HEMANT K. CHAURASIA ◽  
MARK C. THOMPSON
Keyword(s):  
Reynolds Number ◽  
Shear Layer ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Vortical Flow ◽  
Forced Flow ◽  
Recirculation Region ◽  
Two Dimensional ◽  
Optimal Perturbation ◽  
Dimensional Instability

A detailed numerical study of the separating and reattaching flow over a square leading-edge plate is presented, examining the instability modes governing transition from two- to three-dimensional flow. Under the influence of background noise, experiments show that the transition scenario typically is incompletely described by either global stability analysis or the transient growth of dominant optimal perturbation modes. Instead two-dimensional transition effectively can be triggered by the convective Kelvin–Helmholtz (KH) shear-layer instability; although it may be possible that this could be described alternatively in terms of higher-order optimal perturbation modes. At least in some experiments, observed transition occurs by either: (i) KH vortices shedding downstream directly and then almost immediately undergoing three-dimensional transition or (ii) at higher Reynolds numbers, larger vortical structures are shed that are also three-dimensionally unstable. These two paths lead to distinctly different three-dimensional arrangements of vortical flow structures. This paper focuses on the mechanisms underlying these three-dimensional transitions. Floquet analysis of weakly periodically forced flow, mimicking the observed two-dimensional quasi-periodic base flow, indicates that the two-dimensional vortex rollers shed from the recirculation region become globally three-dimensionally unstable at a Reynolds number of approximately 380. This transition Reynolds number and the predicted wavelength and flow symmetries match well with those of the experiments. The instability appears to be elliptical in nature with the perturbation field mainly restricted to the cores of the shed rollers and showing the spatial vorticity distribution expected for that instability type. Indeed an estimate of the theoretical predicted wavelength is also a good match to the prediction from Floquet analysis and theoretical estimates indicate the growth rate is positive. Fully three-dimensional simulations are also undertaken to explore the nonlinear development of the three-dimensional instability. These show the development of the characteristic upright hairpins observed in the experimental dye visualisations. The three-dimensional instability that manifests at lower Reynolds numbers is shown to be consistent with an elliptic instability of the KH shear-layer vortices in both symmetry and spanwise wavelength.

Simulations of Droplet Collisions in Shear Flow

2012 ◽  
Author(s):  
Orest Shardt ◽  
J. J. Derksen ◽  
Sushanta K. Mitra
Keyword(s):  
Reynolds Number ◽  
Shear Flow ◽  
Capillary Number ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Simple Shear Flow ◽  
Number Range ◽  
Droplet Collisions ◽  
Critical Capillary Number ◽  
Boltzmann Method

When droplets collide in a shear flow, they may coalesce or remain separate after the collision. At low Reynolds numbers, droplets coalesce when the capillary number does not exceed a critical value. We present three-dimensional simulations of droplet coalescence in a simple shear flow. We use a free-energy lattice Boltzmann method (LBM) and study the collision outcome as a function of the Reynolds and capillary numbers. We study the Reynolds number range from 0.2 to 1.4 and capillary numbers between 0.1 and 0.5. We determine the critical capillary number for the simulations (0.19) and find that it is does not depend on the Reynolds number. The simulations are compared with experiments on collisions between confined droplets in shear flow. The critical capillary number in the simulations is about a factor of 25 higher than the experimental value.

Three-dimensional dynamics and transition to turbulence in the wake of bluff objects

1992 ◽  
Vol 238 ◽  
pp. 1-30 ◽  
Author(s):  
George Em Karniadakis ◽  
George S. Triantafyllou
Keyword(s):  
Reynolds Number ◽  
Three Dimensional ◽  
Period Doubling ◽  
Reynolds Numbers ◽  
Transition To Turbulence ◽  
Two Dimensional ◽  
Near Wake ◽  
Vortex Filaments ◽  
Average Flow ◽  
Time Average

The wakes of bluff objects and in particular of circular cylinders are known to undergo a ‘fast’ transition, from a laminar two-dimensional state at Reynolds number 200 to a turbulent state at Reynolds number 400. The process has been documented in several experimental investigations, but the underlying physical mechanisms have remained largely unknown so far. In this paper, the transition process is investigated numerically, through direct simulation of the Navier—Stokes equations at representative Reynolds numbers, up to 500. A high-order time-accurate, mixed spectral/spectral element technique is used. It is shown that the wake first becomes three-dimensional, as a result of a secondary instability of the two-dimensional vortex street. This secondary instability appears at a Reynolds number close to 200. For slightly supercritical Reynolds numbers, a harmonic state develops, in which the flow oscillates at its fundamental frequency (Strouhal number) around a spanwise modulated time-average flow. In the near wake the modulation wavelength of the time-average flow is half of the spanwise wavelength of the perturbation flow, consistently with linear instability theory. The vortex filaments have a spanwise wavy shape in the near wake, and form rib-like structures further downstream. At higher Reynolds numbers the three-dimensional flow oscillation undergoes a period-doubling bifurcation, in which the flow alternates between two different states. Phase-space analysis of the flow shows that the basic limit cycle has branched into two connected limit cycles. In physical space the period doubling appears as the shedding of two distinct types of vortex filaments.Further increases of the Reynolds number result in a cascade of period-doubling bifurcations, which create a chaotic state in the flow at a Reynolds number of about 500. The flow is characterized by broadband power spectra, and the appearance of intermittent phenomena. It is concluded that the wake undergoes transition to turbulence following the period-doubling route.

Sign in / Sign up

Export Citation Format

Share Document