Study of Gas Flow in Micronozzles Using an Unstructured DSMC Method

Author(s):  
Ehsan Roohi ◽  
Masoud Darbandi ◽  
Vahid Mirjalili

The current research uses an unstructured direct simulation Monte Carlo (DSMC) method to numerically investigate supersonic and subsonic flow behavior in micro convergent–divergent nozzle over a wide range of rarefied regimes. The current unstructured DSMC solver has been suitably modified via using uniform distribution of particles, employing proper subcell geometry, and benefiting from an advanced molecular tracking algorithm. Using this solver, we study the effects of back pressure, gas/surface interactions (diffuse/specular reflections), and Knudsen number, on the flow field in micronozzles. We show that high viscous force manifesting in boundary layers prevents supersonic flow formation in the divergent section of nozzles as soon as the Knudsen number increases above a moderate magnitude. In order to accurately simulate subsonic flow at the nozzle outlet, it is necessary to add a buffer zone to the end of nozzle. If we apply the back pressure at the outlet, boundary layer separation is observed and a region of backward flow appears inside the boundary layer while the core region of inviscid flow experiences multiple shock-expansion waves. We also show that the wall boundary layer prevents forming shocks in the divergent part. Alternatively, Mach cores appear at the nozzle center followed by bow shocks and an expansion region.

Author(s):  
Stephan Priebe ◽  
Daniel Wilkin ◽  
Andy Breeze-Stringfellow ◽  
Giridhar Jothiprasad ◽  
Lawrence C. Cheung

Abstract Shock/boundary layer interactions (SBLI) are a fundamental fluid mechanics problem relevant in a wide range of applications including transonic rotors in turbomachinery. This paper uses wall-resolved large eddy simulation (LES) to examine the interaction of normal shocks with laminar and turbulent inflow boundary layers in transonic flow. The calculations were performed using GENESIS, a high-order, unstructured LES solver. The geometry created for this study is a transonic passage with a convergent-divergent nozzle that expands the flow to the desired Mach number upstream of the shock and then introduces constant radius curvature to simulate local airfoil camber. The Mach numbers in the divergent section of the transonic passage simulate single stage commercial fan blades. The results predicted with the LES calculations show significant differences between laminar and turbulent SBLI in terms of shock structure, boundary layer separation and transition, and aerodynamic losses. For laminar flow into the shock, significant flow separation and low-frequency unsteadiness occur, while for turbulent flow into the shock, both the boundary layer loss and the low-frequency unsteadiness are reduced.


Author(s):  
Jeffrey P. Bons ◽  
Rolf Sondergaard ◽  
Richard B. Rivir

The effects of pulsed vortex generator jets on a naturally separating low pressure turbine boundary layer have been investigated experimentally. Blade Reynolds numbers in the linear turbine cascade match those for high altitude aircraft engines and industrial turbine engines with elevated turbine inlet temperatures. The vortex generator jets (30 degree pitch and 90 degree skew angle) are pulsed over a wide range of frequency at constant amplitude and selected duty cycles. The resulting wake loss coefficient vs. pulsing frequency data add to previously presented work by the authors documenting the loss dependency on amplitude and duty cycle. As in the previous studies, vortex generator jets are shown to be highly effective in controlling laminar boundary layer separation. This is found to be true at dimensionless forcing frequencies (F+) well below unity and with low (10%) duty cycles. This unexpected low frequency effectiveness is due to the relatively long relaxation time of the boundary layer as it resumes its separated state. Extensive phase-locked velocity measurements taken in the blade wake at an F+ of 0.01 with 50% duty cycle (a condition at which the flow is essentially quasi-steady) document the ejection of bound vorticity associated with a low momentum fluid packet at the beginning of each jet pulse. Once this initial fluid event has swept down the suction surface of the blade, a reduced wake signature indicates the presence of an attached boundary layer until just after the jet termination. The boundary layer subsequently relaxes back to its naturally separated state. This relaxation occurs on a timescale which is 5–6 times longer than the original attachment due to the starting vortex. Phase-locked boundary layer measurements taken at various stations along the blade chord illustrate this slow relaxation phenomenon. This behavior suggests that some economy of jet flow may be possible by optimizing the pulse duty cycle and frequency for a particular application. At higher pulsing frequencies, for which the flow is fully dynamic, the boundary layer is dominated by periodic shedding and separation bubble migration, never recovering its fully separated (uncontrolled) state.


2006 ◽  
Author(s):  
J. M. Jones ◽  
D. K. Walters

This paper presents the initial development and validation of a modified two-equation eddy-viscosity turbulence model for computational fluid dynamics (CFD) prediction of transitional and turbulent flow. The new model is based on a k-ω model framework, making it more easily implemented into existing general-purpose CFD solvers than other recently proposed model forms. The model incorporates inviscid and viscous damping functions for the eddy viscosity, as well as a production damping term, in order to reproduce the appropriate effects of laminar, transitional, and turbulent boundary layer flow. It has been implemented into a commercially available flow solver (FLUENT) and evaluated for simple attached and separated flow conditions, including 2-D flow over a flat plate and a circular cylinder. The results presented show that the new model is able to yield reasonable predictions of transitional flow behavior using a very simple modeling framework, including an appropriate response to freestream turbulence and boundary layer separation.


2010 ◽  
Vol 659 ◽  
pp. 289-317 ◽  
Author(s):  
PAYAM AGHSAEE ◽  
LEON BOEGMAN ◽  
KEVIN G. LAMB

The breaking of fully nonlinear internal solitary waves of depression shoaling upon a uniformly sloping boundary in a smoothed two-layer density field was investigated using high-resolution two-dimensional simulations. Our simulations were limited to narrow-crested waves, which are more common than broad-crested waves in geophysical flows. The simulations were performed for a wide range of boundary slopes S ∈ [0.01, 0.3] and wave slopes extending the parameter range to weaker slopes than considered in previous laboratory and numerical studies. Over steep slopes (S ≥ 0.1), three distinct breaking processes were observed: surging, plunging and collapsing breakers which are associated with reflection, convective instability and boundary-layer separation, respectively. Over mild slopes (S ≤ 0.05), nonlinearity varies gradually and the wave fissions into a train of waves of elevation as it passes through the turning point where solitary waves reverse polarity. The dynamics of each breaker type were investigated and the predominance of a particular mechanism was associated with a relative developmental time scale. The breaking location was modelled as a function of wave amplitude (a), characteristic wave length and the isopycnal length along the slope. The breaker type was characterized in wave slope (Sw = a/Lw, where Lw is a measure of half of the wavelength) versus S space, and the reflection coefficient (R), modelled as a function of the internal Iribarren number, was in agreement with other studies. The effects of grid resolution and wave Reynolds number (Rew) on R, boundary-layer separation and the evolution of global instability were studied. High Reynolds numbers (Rew ~ 104) were found to trigger a global instability, which modifies the breaking process relative to the lower Rew case, but not necessarily the breaking location, and results in a ~ 10 % increase in R, relative to the Rew ~ 103 case.


2010 ◽  
Vol 14 (suppl.) ◽  
pp. 183-196
Author(s):  
Slobodan Savic ◽  
Branko Obrovic ◽  
Milan Despotovic ◽  
Dusan Gordic

This paper studies the influence of the magnetic field on the planar laminar steady flow of the ionized gas in the boundary layer. The present outer magnetic field is homogenous and perpendicular to the body within the fluid. The gas of the same physical characteristics as the gas in the main flow is injected (ejected) through the contour of the body. The governing boundary layer equations for different forms of the electroconductivity variation law are transformed, brought to a generalized form and solved numerically in a four-parametric approximation. It has been determined that the magnetic field, through the magnetic parameter, has a great influence on certain quantities and characteristics of the boundary layer. It has also been shown that this parameter has an especially significant influence on the non-dimensional friction function, and hence the boundary layer separation point.


2001 ◽  
Vol 124 (1) ◽  
pp. 77-85 ◽  
Author(s):  
Jeffrey P. Bons ◽  
Rolf Sondergaard ◽  
Richard B. Rivir

The effects of pulsed vortex generator jets on a naturally separating low-pressure turbine boundary layer have been investigated experimentally. Blade Reynolds numbers in the linear turbine cascade match those for high-altitude aircraft engines and industrial turbine engines with elevated turbine inlet temperatures. The vortex generator jets (30 deg pitch and 90 deg skew angle) are pulsed over a wide range of frequency at constant amplitude and selected duty cycles. The resulting wake loss coefficient versus pulsing frequency data add to previously presented work by the authors documenting the loss dependency on amplitude and duty cycle. As in the previous studies, vortex generator jets are shown to be highly effective in controlling laminar boundary layer separation. This is found to be true at dimensionless forcing frequencies F+ well below unity and with low (10 percent) duty cycles. This unexpected low-frequency effectiveness is due to the relatively long relaxation time of the boundary layer as it resumes its separated state. Extensive phase-locked velocity measurements taken in the blade wake at an F+ of 0.01 with 50 percent duty cycle (a condition at which the flow is essentially quasi-steady) document the ejection of bound vorticity associated with a low-momentum fluid packet at the beginning of each jet pulse. Once this initial fluid event has swept down the suction surface of the blade, a reduced wake signature indicates the presence of an attached boundary layer until just after the jet termination. The boundary layer subsequently relaxes back to its naturally separated state. This relaxation occurs on a timescale which is five to six times longer than the original attachment due to the starting vortex. Phase-locked boundary layer measurements taken at various stations along the blade chord illustrate this slow relaxation phenomenon. This behavior suggests that some economy of jet flow may be possible by optimizing the pulse duty cycle and frequency for a particular application. At higher pulsing frequencies, for which the flow is fully dynamic, the boundary layer is dominated by periodic shedding and separation bubble migration, never recovering its fully separated (uncontrolled) state.


Author(s):  
Shuzhen Hu ◽  
Yanfeng Zhang ◽  
Xue Feng Zhang ◽  
Edward Vlasic

The inter-turbine transition duct (ITD) of a gas turbine engine has significant potential for engine weight reduction and/or aerodynamic performance improvement. This potential arises because very little is understood of the flow behavior in the duct in relation to the hub and casing shapes and the flow entering the duct (e.g., swirl angle, turbulence intensity, periodic unsteadiness and blade tip vortices from upstream HP turbine blade rows). In this study, the flow development in an ITD with different inlet swirl distributions was investigated experimentally and numerically. The current paper, which is the first part of a two-part paper, presents the investigations of the influences of the casing swirl variations on the flow physics in the ITD. The results show a fair agreement between the predicted and experimental data. The radial pressure gradient at the first bend of ITD drives the low momentum hub boundary layer and wake flow radially, which results in a pair of hub counter-rotating vortices. Furthermore, the radially moving low momentum wake flow feeds into the casing region and causes 3D casing boundary layer. At the second bend, the reversed radial pressure gradient together with the 3D casing boundary layer generates a pair of casing counter-rotating vortices. Due to the local adverse pressure gradient, 3D boundary layer separation occurs on both the casing and hub at the second bend and the exit of the ITD, respectively. The casing 3D separation enhances the 3D features of the casing boundary layer as well as the existing casing counter-rotating vortices. With increasing casing swirl angle, the casing 3D boundary layer separation is delayed and the casing counter-rotating vortices are weakened. On the other hand, although the hub swirls are kept constant, the hub counter-rotating vortices get stronger with the increasing inlet swirl gradient. The total pressure coefficients within the ITD are significantly redistributed by the casing and hub counter-rotating vortices.


2012 ◽  
Vol 694 ◽  
pp. 191-224 ◽  
Author(s):  
Satoshi Taguchi ◽  
Kazuo Aoki

AbstractA rarefied gas flow thermally induced around a heated (or cooled) flat plate, contained in a vessel, is considered in two different situations: (i) both sides of the plate are simultaneously and uniformly heated (or cooled); and (ii) only one side of the plate is uniformly heated. The former is known as the thermal edge flow and the latter, typically observed in the Crookes radiometer, may be called the radiometric flow. The steady behaviour of the gas induced in the container is investigated on the basis of the Bhatnagar–Gross–Krook (BGK) model of the Boltzmann equation and the diffuse reflection boundary condition by means of an accurate finite-difference method. The flow features are clarified for a wide range of the Knudsen number, with a particular emphasis placed on the structural similarity between the two flows. The limiting behaviour of the flow as the Knudsen number tends to zero (and thus the system approaches the continuum limit) is investigated for both flows. The detailed structure of the normal stress on the plate as well as the cause of the radiometric force (the force acting on the plate from the hotter to the colder side) is also clarified for the present infinitely thin plate.


2004 ◽  
Vol 127 (4) ◽  
pp. 827-834 ◽  
Author(s):  
Masayoshi Shimo ◽  
James V. Canino ◽  
Stephen D. Heister

Oil flow behavior within a lubrication system utilized in a turbofan engine has been studied using a two-dimensional model for a seal runner and a sump wall. A two-dimensional, axisymmetric boundary layer method is utilized to derive a model for the film on the seal runner. An integral method analysis of boundary layer equations are utilized to derive a model for the film behavior on the outer wall of the sump neglecting variations in the axial direction. Parametric studies for oil flow at the seal runner and one the sump wall have been generated for a wide range of oil film properties and ambient conditions.


2019 ◽  
Vol 2019 ◽  
pp. 1-14
Author(s):  
Mohamed Hssikou ◽  
Jamal Baliti ◽  
Mohammed Alaoui

A Maxwell gas confined within a micro cavity with nonisothermal walls is investigated in the slip and early transition regimes using the classical and extended continuum theories. The vertical sides of the cavity are kept at the uniform and environmental temperature T0, while the upper and bottom ones are linearly heated in opposite directions from the cold value T0 to the hot one TH. The gas flow is, therefore, induced only by the temperature gradient created along the longitudinal walls. The problem is treated from a macroscopic point of view by solving numerically the so-called regularized 13-moment equations (R13) recently developed as an extension of Grad 13-moment theory to the third order of the Knudsen number powers in the Chapman-Enskog expansion. The gas macroscopic properties obtained by this method are compared with the classical continuum theory results (NSF) using the first and second order of velocity slip and temperature jump boundary conditions. The gas flow behavior is studied as a function of the Knudsen number (Kn), nonlinear effects, for different heating rates T0/TH. The micro cavity aspect ratio effect is also evaluated on the flow fields in this study.


Sign in / Sign up

Export Citation Format

Share Document