Abstract: In this study OpenFOAM 4.4.2 was used to investigate flow inside the coronary artery of the heart. This step is the first step of our future project, which is to include conjugate heat transfer of the heart with three main coronary arteries. Three different velocities were used as inlet boundary conditions to see the effect of velocity increase on velocity, pressure, and wall shear of the coronary artery. Also, three different fluids, namely the University of Wisconsin solution, gelatin, and blood was used to investigate the effect of different fluids on flow inside the coronary artery. A code based on Reynolds Stress Navier Stokes (RANS) equations was written and implemented with the real boundary condition that was calculated based on MRI images. In order to improve the accuracy of the current numerical scheme, hex dominant mesh is utilized. When the inlet velocity increases to 0.5 m/s, velocity, wall shear stress, and pressure increase at the narrower parts.
Abstract: Cavitation inside a diesel injector nozzle is investigated numerically in this study. The Reynolds Stress Navier Stokes set of equations (RANS) are utilized to investigate flow behavior inside the nozzle numerically. Moreover, K-ε turbulent model is found to be a better approach comparing to K-ω turbulent model. The Winklhofer rectangular shape nozzle is also simulated in order to verify the current numerical scheme, and with the mass flow rate approach, the current solution is verified. Afterward, a six-hole real size nozzle was simulated and it was found that among the different fuels used in this study with the same condition, diesel fuel provides the largest length of cavitation. Also, it was found that at the same boundary condition, rapeseed methyl ester (RME) fuel leads to the highest value of discharge coefficient and mass flow rate.
Abstract: The influence of the geometric parameters of trapezoidal labyrinth channel on the emitter discharge is investigated in this work. The impact of the dentate angle, the dentate spacing, and the dentate height are studied among the geometric parameters of the labyrinth channel. Numerical simulations of the water flow movement are performed according to central cubic composite design using Commercial codes GAMBIT and FLUENT. Inlet pressure of the dripper is set up to be 1 bar. The objective of this paper is to derive a mathematical model of the emitter discharge depending on the dentate angle, the dentate spacing, the dentate height of the labyrinth channel. As a result, the obtained mathematical model is a second-order polynomial reporting 2-way interactions among the geometric parameters. The dentate spacing has the most important and positive influence on the emitter discharge, followed by the simultaneous impact of the dentate spacing and the dentate height. The dentate angle in the observed interval has no significant effect on the emitter discharge. The obtained model can be used as a basis for a future emitter design.
Abstract: Turbulence modelling is still evolving, and efforts are on to improve and develop numerical methods to simulate the real turbulence structures by using the empirical and experimental information. The monotonically integrated large eddy simulation (MILES) is an attractive approach for modelling turbulence in high Re flows, which is based on the solving of the unfiltered flow equations with no explicit sub-grid scale (SGS) model. In the current work, this approach has been used, and the action of the SGS model has been included implicitly by intrinsic nonlinear high-frequency filters built into the convection discretization schemes. The MILES solver is developed using the opensource CFD OpenFOAM libraries. The role of flux limiters schemes namely, Gamma, superBee, van-Albada and van-Leer, is studied in predicting turbulent statistical quantities for a fully developed channel flow with a friction Reynolds number, ReT = 180, and compared the numerical predictions with the well-established Direct Numerical Simulation (DNS) results for studying the wall generated turbulence. It is inferred from the numerical predictions that Gamma, van-Leer and van-Albada limiters produced more diffusion and overpredicted the velocity profiles, while superBee scheme reproduced velocity profiles and turbulence statistical quantities in good agreement with the reference DNS data in the streamwise direction although it deviated slightly in the spanwise and normal to the wall directions. The simulation results are further discussed in terms of the turbulence intensities and Reynolds stresses averaged in time and space to draw conclusion on the flux limiter schemes performance in OpenFOAM context.
Abstract: Eddy viscosity models in turbulence modeling can be mainly classified as linear and nonlinear models. Linear formulations are simple and require less computational resources but have the disadvantage that they cannot predict actual flow pattern in complex geophysical flows where streamline curvature and swirling motion are predominant. A constitutive equation of Reynolds stress anisotropy is adopted for the formulation of eddy viscosity including all the possible higher order terms quadratic in the mean velocity gradients, and a simplified model is developed for actual oceanic flows where only the vertical velocity gradients are important. The new model is incorporated into the one dimensional General Ocean Turbulence Model (GOTM). Two realistic oceanic test cases (OWS Papa and FLEX' 76) have been investigated. The new model predictions match well with the observational data and are better in comparison to the predictions of the two equation k-epsilon model. The proposed model can be easily incorporated in the three dimensional Princeton Ocean Model (POM) to simulate a wide range of oceanic processes. Practically, this model can be implemented in the coastal regions where trasverse shear induces higher vorticity, and for prediction of flow in estuaries and lakes, where depth is comparatively less. The model predictions of marine turbulence and other related data (e.g. Sea surface temperature, Surface heat flux and vertical temperature profile) can be utilized in short term ocean and climate forecasting and warning systems.
Abstract: The paper presents the results of a detailed assessment of several modern Reynolds Averaged Navier-Stokes (RANS) turbulence models for prediction of C3X vane film cooling at various injection regimes. Three models are considered, namely the Shear Stress Transport (SST) model, the modification of the SST model accounting for the streamlines curvature (SST-CC), and the Explicit Algebraic Reynolds Stress Model (EARSM). It is shown that all the considered models face with a problem in prediction of the adiabatic effectiveness in the vicinity of the cooling holes; however, accounting for the Reynolds stress anisotropy within the EARSM model noticeably increases the solution accuracy. On the other hand, further downstream all the models provide a reasonable agreement with the experimental data for the adiabatic effectiveness and among the considered models the most accurate results are obtained with the use EARMS.
Abstract: The turbulent structures in the wake (x/d =10 to 60) of a screen cylinder have been educed to understand the roles of the various structures as evolving downstream by comparing with those obtained in a solid circular cylinder wake at Reynolds number, Re of 7000. Using a wavelet multiresolution technique, the flow structures are decomposed into a number of wavelet components based on their central frequencies. It is observed that in the solid cylinder wake, large-scale structures (of frequencyf0 and 1.2 f0) make the largest contribution to the Reynolds stresses although they start to lose their roles significantly at x/d> 20. In the screen cylinder wake, the intermediate-scale structures (2f0 and 4f0) contribute the most to the Reynolds stresses atx/d =10 before being taken over by the large-scale structures (f0) further downstream.
Abstract: Among the various cooling processes in industrial
applications such as: electronic devices, heat exchangers, gas
turbines, etc. Gas turbine blades cooling is the most challenging one.
One of the most common practices is using ribbed wall because of
the boundary layer excitation and therefore making the ultimate
cooling. Vortex formation between rib and channel wall will result in
a complicated behavior of flow regime. At the other hand, selecting
the most efficient method for capturing the best results comparing to
experimental works would be a fascinating issue. In this paper 4
common methods in turbulence modeling: standard k-e, rationalized
k-e with enhanced wall boundary layer treatment, k-w and RSM
(Reynolds stress model) are employed to a square ribbed channel to
investigate the separation and thermal behavior of the flow in the
channel. Finally all results from different methods which are used in
this paper will be compared with experimental data available in
literature to ensure the numerical method accuracy.
Abstract: Numerical analysis of flow characteristics and
separation efficiency in a high-efficiency cyclone has been performed.
Several models based on the experimental observation for a design
purpose were proposed. However, the model is only estimated the
cyclone's performance under the limited environments; it is difficult to
obtain a general model for all types of cyclones. The purpose of this
study is to find out the flow characteristics and separation efficiency
numerically. The Reynolds stress model (RSM) was employed instead
of a standard k-ε or a k-ω model which was suitable for isotropic
turbulence and it could predict the pressure drop and the Rankine
vortex very well. For small particles, there were three significant
components (entrance of vortex finder, cone, and dust collector) for
the particle separation. In the present work, the particle re-entraining
phenomenon from the dust collector to the cyclone body was observed
after considerable time. This re-entrainment degraded the separation
efficiency and was one of the significant factors for the separation
efficiency of the cyclone.
Abstract: Linear stochastic estimation and quadratic stochastic
estimation techniques were applied to estimate the entire velocity
flow-field of an open cavity with a length to depth ratio of 2. The
estimations were done through the use of instantaneous velocity
magnitude as estimators. These measurements were obtained by
Particle Image Velocimetry. The predicted flow was compared
against the original flow-field in terms of the Reynolds stresses and
turbulent kinetic energy. Quadratic stochastic estimation proved to be
more superior than linear stochastic estimation in resolving the shear
layer flow. When the velocity fluctuations were scaled up in the
quadratic estimate, both the time-averaged quantities and the
instantaneous cavity flow can be predicted to a rather accurate extent.
Abstract: The objective of this research was to find the diffusion properties of vehicles on the road by using the V-Sphere Code. The diffusion coefficient and the size of the height of the wake were estimated with the LES option and the third order MUSCL scheme. We evaluated the code with the changes in the moments of Reynolds Stress along the mean streamline. The results show that at the leading part of a bluff body the LES has some advantages over the RNS since the changes in the strain rates are larger for the leading part. We estimated that the diffusion coefficient with the computed Reynolds stress (non-dimensional) was about 0.96 times the mean velocity.
Abstract: Turbulence studies were made in the wake of a rotating
circular cylinder in a uniform free stream. The interest was to
examine the turbulence properties at the suppression of periodicity in
vortex formation process. An experimental study of the turbulent near
wake of a rotating circular cylinder was made at a Reynolds number
of 9000 for velocity ratios, λ between 0 and 2.7. Hot-wire
anemometry and particle image velocimetry results indicate that the
rotation of the cylinder causes significant changes in the vortical
activities. The turbulence quantities are getting smaller as λ increases
due to suppression of coherent vortex structures.
Abstract: A numerical investigation of the effects of nanosecond
barrier discharge on the stability of a two-dimensional free shear layer
is performed. The computations are carried out using a compressible
Navier-Stokes algorithm coupled with a thermodynamic model of the
discharge. The results show that significant increases in the shear
layer-s momentum thickness and Reynolds stresses occur due to
actuation. Dependence on both frequency and amplitude of actuation
are considered, and a comparison is made of the computed growth
rates with those predicted by linear stability theory. Amplitude and
frequency ranges for the efficient promotion of shear-layer instabilities
are identified.
Abstract: This paper analytically investigates the 3D flow
pattern at the confluences of two rectangular channels having 900
angles using Navier-Stokes equations based on Reynolds Stress
Turbulence Model (RSM). The equations are solved by the Finite-
Volume Method (FVM) and the flow is analyzed in terms of steadystate
(single-phased) conditions. The Shumate experimental findings
were used to test the validity of data. Comparison of the simulation
model with the experimental ones indicated a close proximity
between the flow patterns of the two sets. Effects of the discharge
ratio on separation zone dimensions created in the main-channel
downstream of the confluence indicated an inverse relation, where a
decrease in discharge ratio, will entail an increase in the length and
width of the separation zone. The study also found the model as a
powerful analytical tool in the feasibility study of hydraulic
engineering projects.
Abstract: This paper looks into detailed investigation of
thermal-hydraulic characteristics of the flow field in a fuel rod
model, especially near the spacer. The area investigate represents a
source of information on the velocity flow field, vortex, and on the
amount of heat transfer into the coolant all of which are critical for
the design and improvement of the fuel rod in nuclear power plants.
The flow field investigation uses three-dimensional Computational
Fluid Dynamics (CFD) with the Reynolds stresses turbulence model
(RSM). The fuel rod model incorporates a vertical annular channel
where three different shapes of spacers are used; each spacer shape is
addressed individually. These spacers are mutually compared in
consideration of heat transfer capabilities between the coolant and
the fuel rod model. The results are complemented with the calculated
heat transfer coefficient in the location of the spacer and along the
stainless-steel pipe.
Abstract: The present paper represents a methodology for
investigating flow characteristics near orifice plate by using a
commercial computational fluid dynamics code. The flow
characteristics near orifice plate which is located in the auxiliary
feedwater system were modeled via three different levels of grid and
four different types of Reynolds Averaged Navier-Stokes (RANS)
equations with proper near-wall treatment. The results from CFD code
were compared with experimental data in terms of differential pressure
through the orifice plate. In this preliminary study, the Realizable k-ε
and the Reynolds stress models with enhanced wall treatment were
suitable to analyze flow characteristics near orifice plate, and the
results had a good agreement with experimental data.
Abstract: In this research the separation efficiency of deoiling hydrocyclone is evaluated using three-dimensional simulation of multiphase flow based on Eulerian-Eulerian finite volume method. The mixture approach of Reynolds Stress Model is also employed to capture the features of turbulent multiphase swirling flow. The obtained separation efficiency of Colman's design is compared with available experimental data and showed that the separation curve of deoiling hydrocyclones can be predicted using numerical simulation.
Abstract: Circle grid space filling plate is a flow conditioner with a fractal pattern and used to eliminate turbulence originating from pipe fittings in experimental fluid flow applications. In this paper, steady state, incompressible, swirling turbulent flow through circle grid space filling plate has been studied. The solution and the analysis were carried out using finite volume CFD solver FLUENT 6.2. Three turbulence models were used in the numerical investigation and their results were compared with the pressure drop correlation of BS EN ISO 5167-2:2003. The turbulence models investigated here are the standard k-ε, realizable k-ε, and the Reynolds Stress Model (RSM). The results showed that the RSM model gave the best agreement with the ISO pressure drop correlation. The effects of circle grids space filling plate thickness and Reynolds number on the flow characteristics have been investigated as well.
Abstract: For the first time since 1940 and presentation of
theodorson-s theory, distribution of thrust, torque and efficiency
along the blade of a counter rotating propeller axial fan was studied
with a novel method in this research. A constant chord, constant
pitch symmetric fan was investigated with Reynolds Stress
Turbulence method in this project and H.E.S. method was utilized to
obtain distribution profiles from C.F.D. tests outcome. C.F.D. test
results were validated by estimation from Playlic-s analytical method.
Final results proved ability of H.E.S. method to obtain distribution
profiles from C.F.D test results and demonstrated interesting facts
about effects of solidity and differences between distributions in front
and rear section.