Abstract: Detailed numerical calculations are illustrated in our investigation for unsteady natural convection heat and mass transfer of non-Newtonian Casson fluid along a vertical wavy surface. The surface of the plate is kept at a constant temperature and uniform concentration. To transform the complex wavy surface to a flat plate, a simple coordinate transformation is employed. The resulting partial differential equations are solved using the fully implicit finite difference method with SUR procedure. Flow and heat transfer characteristics are investigated for a wide range of values of the Casson parameter, the dimensionless time parameter, the buoyancy ratio and the amplitude-wavelength parameter. It is found that, the variations of the Casson parameter have significant effects on the fluid motion, heat and mass transfer. Also, the maximum and minimum values of the local Nusselt and Sherwood numbers increase by increase either the Casson parameter or the buoyancy ratio.
Abstract: Supersonic wind tunnel nozzles are generally capable of producing a constant Mach number flow in the test section of the wind tunnel. As a result, most of the supersonic vehicles are widely designed using steady state flow characteristics which may have errors while facing unsteady situations. This study aims to explore the possibility of varying the Mach number of the flow during wind tunnel operation. The nozzle walls are restricted to be inflexible for cooling near the throat due to high stagnation temperature requirement of the flow to simulate the conditions as experienced by the vehicle. Two simple independent mechanisms, rotation and translation of nozzle walls have been analyzed and the nozzle ranges have been optimized to vary the Mach number from Mach 2 to Mach 5 using minimum number of nozzles in the wind tunnel.
Abstract: The prediction of the flow around marine propellers and vessel hulls propeller interaction is one of the challenges of Computational fluid dynamics (CFD). The CFD has emerged as a potential tool in recent years and has promising applications. The objective of the current study is to predict the hydrodynamic performances of HSP marine propeller in open water and behind a vessel. The unsteady 3-D flow was modeled numerically along with respectively the K-ω standard and K-ω SST turbulence models for steady and unsteady cases. The hydrodynamic performances such us a torque and thrust coefficients and efficiency show good agreement with the experiment results.
Abstract: Complicated unsteady flow in axial turbines produces high-frequency unsteady aerodynamic exciting force, which threatens the safe operation of turbines. This paper illustrates how negative-bowed stator reduces the rotor unsteady aerodynamic exciting force by unsteady flow field. With the support of three-dimensional viscous compressible Navier-Stokes equation, the single axial turbines with 0, -10 and -20 degree bowed stator are comparably investigated, aiming to identify the flow field structure difference caused by various negative-bowed degrees. The results show that negative-bowed stator strengthens the turbulence kinetic energy, which is further strengthened with the increase of negative-bowed degree. Meanwhile, the flow phenomenon including stator wakes and passage vortex is shown. In addition, the interaction of upstream negative-bowed wakes contributes to the reduction of unsteady blade load fluctuation. Furthermore, the aerodynamic exciting force decreases with the increasing negative bowed degree, while the efficiency is correspondingly reduced. This paper provides the reference for the alleviation of the harmful impact caused by unsteady interaction with the method of wake control.
Abstract: Additive manufacturing (AM) is increasingly crucial in biomedical and aerospace industries. As a recently developed AM technique, selective laser melting (SLM) has become a commercial method for various manufacturing processes. However, the molten pool configuration during SLM of metal powders is a decisive issue for the product quality. It is very important to investigate the heat transfer characteristics during the laser heating process. In this work, the finite element method (FEM) software ANSYS® (work bench module 16.0) was used to predict the unsteady temperature distribution for resolving molten pool dimensions with consideration of temperature-dependent thermal physical properties of TiAl6V4 at different laser powers and scanning speeds. The simulated results of the temperature distributions illustrated that the ratio of laser power to scanning speed can greatly influence the size of molten pool of titanium alloy powder for SLM development.
Abstract: A local nonlinear stability analysis using a eight-mode
expansion is performed in arriving at the coupled amplitude equations
for Rayleigh-Bénard-Brinkman convection (RBBC) in the presence
of LTNE effects. Streamlines and isotherms are obtained in the
two-dimensional unsteady finite-amplitude convection regime. The
parameters’ influence on heat transport is found to be more
pronounced at small time than at long times. Results of the
Rayleigh-Bénard convection is obtained as a particular case of
the present study. Additional modes are shown not to significantly
influence the heat transport thus leading us to infer that five minimal
modes are sufficient to make a study of RBBC. The present problem
that uses rolls as a pattern of manifestation of instability is a needed
first step in the direction of making a very general non-local study of
two-dimensional unsteady convection. The results may be useful in
determining the preferred range of parameters’ values while making
rheometric measurements in fluids to ascertain fluid properties such
as viscosity. The results of LTE are obtained as a limiting case of
the results of LTNE obtained in the paper.
Abstract: Rayleigh-B´enard convection of a nanoliquid in shallow, square and tall enclosures is studied using the Khanafer-Vafai-Lightstone single-phase model. The thermophysical properties of water, copper, copper-oxide, alumina, silver and titania at 3000 K under stagnant conditions that are collected from literature are used in calculating thermophysical properties of water-based nanoliquids. Phenomenological laws and mixture theory are used for calculating thermophysical properties. Free-free, rigid-rigid and rigid-free boundary conditions are considered in the study. Intractable Lorenz model for each boundary combination is derived and then reduced to the tractable Ginzburg-Landau model. The amplitude thus obtained is used to quantify the heat transport in terms of Nusselt number. Addition of nanoparticles is shown not to alter the influence of the nature of boundaries on the onset of convection as well as on heat transport. Amongst the three enclosures considered, it is found that tall and shallow enclosures transport maximum and minimum energy respectively. Enhancement of heat transport due to nanoparticles in the three enclosures is found to be in the range 3% - 11%. Comparison of results in the case of rigid-rigid boundaries is made with those of an earlier work and good agreement is found. The study has limitations in the sense that thermophysical properties are calculated by using various quantities modelled for static condition.
Abstract: The purpose of this work is to simulate the flow at the exit of Vulcan 1 engine of European launcher Ariane 5. The geometry of the propellant nozzle is already determined using the characteristics method. The pressure in the outlet section of the nozzle is less than atmospheric pressure on the ground, causing the existence of oblique and normal shock waves at the exit. During the rise of the launcher, the atmospheric pressure decreases and the shock wave disappears. The code allows the capture of shock wave at exit of nozzle. The numerical technique uses the Flux Vector Splitting method of Van Leer to ensure convergence and avoid the calculation instabilities. The Courant, Friedrichs and Lewy coefficient (CFL) and mesh size level are selected to ensure the numerical convergence. The nonlinear partial derivative equations system which governs this flow is solved by an explicit unsteady numerical scheme by the finite volume method. The accuracy of the solution depends on the size of the mesh and also the step of time used in the discretized equations. We have chosen in this study the mesh that gives us a stationary solution with good accuracy.
Abstract: Numerical investigations are performed to analyze the flow behavior over NACA0015 and to evaluate the efficiency of synthetic jet as active control device. The second objective of this work is to investigate the influence of momentum coefficient of synthetic jet on the flow behaviour. The unsteady Reynolds-averaged Navier-Stokes equations of the turbulent flow are solved using, k-ω SST provided by ANSYS CFX-CFD code. The model presented in this paper is a comprehensive representation of the information found in the literature. Comparison of obtained numerical flow parameters with the experimental ones shows that the adopted computational procedure reflects nearly the real flow nature. Also, numerical results state that use of synthetic jets devices has positive effects on the flow separation, and thus, aerodynamic performance improvement of NACA0015 airfoil. It can also be observed that the use of synthetic jet increases the lift coefficient about 13.3% and reduces the drag coefficient about 52.7%.
Abstract: The unsteady supersonic jet formed by a shock tube with a small high-pressure chamber was used as a simple alternative model for pulsed laser ablation. Understanding the vortex ring formed by the shock wave is crucial in clarifying the behavior of unsteady supersonic jet discharged from an elliptical cell. Therefore, this study investigated the behavior of vortex rings and a jet. The experiment and numerical calculation were conducted using the schlieren method and by solving the axisymmetric two-dimensional compressible Navier–Stokes equations, respectively. In both, the calculation and the experiment, laser ablation is conducted for a certain duration, followed by discharge through the exit. Moreover, a parametric study was performed to demonstrate the effect of pressure ratio on the interaction among vortex rings and the supersonic jet. The interaction between the supersonic jet and the vortex rings increased the velocity of the supersonic jet up to the magnitude of the velocity at the center of the vortex rings. The interaction between the vortex rings increased the velocity at the center of the vortex ring.
Abstract: Numerical simulations of vortex-induced vibration of a three-dimensional flexible tube under uniform turbulent flow are calculated when Reynolds number is 1.35×104. In order to achieve the vortex-induced vibration, the three-dimensional unsteady, viscous, incompressible Navier-Stokes equation and LES turbulence model are solved with the finite volume approach, the tube is discretized according to the finite element theory, and its dynamic equilibrium equations are solved by the Newmark method. The fluid-tube interaction is realized by utilizing the diffusion-based smooth dynamic mesh method. Considering the vortex-induced vibration system, the variety trends of lift coefficient, drag coefficient, displacement, vertex shedding frequency, phase difference angle of tube are analyzed under different frequency ratios. The nonlinear phenomena of locked-in, phase-switch are captured successfully. Meanwhile, the limit cycle and bifurcation of lift coefficient and displacement are analyzed by using trajectory, phase portrait, and Poincaré sections. The results reveal that: when drag coefficient reaches its minimum value, the transverse amplitude reaches its maximum, and the “lock-in” begins simultaneously. In the range of lock-in, amplitude decreases gradually with increasing of frequency ratio. When lift coefficient reaches its minimum value, the phase difference undergoes a suddenly change from the “out-of-phase” to the “in-phase” mode.
Abstract: In aerovehicles context, the flow around an Ahmed
body profile is simulated using the velocity-vorticity formulation of
the Navier-Stokes equations, associated to a penalization method for
solids and Large Eddy Simulation for turbulence. The study focuses
both on the ground influence on the flow and on the dissymetry of
the wake, observed for a ground clearance greater than 10% of the
body height H. Unsteady and mean flows are presented and analyzed.
POD study completes the analysis and gives information on the most
energetic structures of the flow.
Abstract: Present study investigates the effect of unsteady wakes on heat transfer in blade tip. Heat/mass transfer was measured in blade tip region depending on a variety of strouhal number by naphthalene sublimation technique. Naphthalene sublimation technique measures heat transfer using a heat/mass transfer analogy. Experiments are performed in linear cascade which is composed of five turbine blades and rotating rods. Strouhal number of inlet flow are changed ranging from 0 to 0.22. Reynolds number is 100,000 based on 11.4 m/s of outlet flow and axial chord length. Three different squealer tip geometries such as base squealer tip, vertical rib squealer tip, and camber line squealer tip are used to study how unsteady wakes affect heat transfer on a blade tip. Depending on squealer tip geometry, different flow patterns occur on a blade tip. Also, unsteady wakes cause reduced tip leakage flow and turbulent flow. As a result, as strouhal number increases, heat/mass transfer coefficients decrease due to the reduced leakage flow. As strouhal number increases, heat/ mass transfer coefficients on a blade tip increase in vertical rib squealer tip.
Abstract: Relative permeabilities are practical factors that are used to correct the single phase Darcy’s law for application to multiphase flow. For effective characterisation of large-scale multiphase flow in hydrocarbon recovery, relative permeability and capillary pressures are used. These parameters are acquired via special core flooding experiments. Special core analysis (SCAL) module of reservoir simulation is applied by engineers for the evaluation of these parameters. But, core flooding experiments in shale core sample are expensive and time consuming before various flow assumptions are achieved for instance Darcy’s law. This makes it imperative for the application of coreflooding simulations in which various analysis of relative permeabilities and capillary pressures of multiphase flow can be carried out efficiently and effectively at a relative pace. This paper presents a Sendra software simulation of core flooding to achieve to relative permeabilities and capillary pressures using different correlations. The approach used in this study was three steps. The first step, the basic petrophysical parameters of Marcellus shale sample such as porosity was determined using laboratory techniques. Secondly, core flooding was simulated for particular scenario of injection using different correlations. And thirdly the best fit correlations for the estimation of relative permeability and capillary pressure was obtained. This research approach saves cost and time and very reliable in the computation of relative permeability and capillary pressures at steady or unsteady state, drainage or imbibition processes in oil and gas industry when compared to other methods.
Abstract: The damage caused by surge waves generated in water
bodies by flow-like landslides can be very high in terms of human
lives and economic losses. The complicated phenomena occurred in
this highly unsteady process are difficult to model because three
interacting phases: air, water and sediment are involved. The problem
therefore is challenging since the effects of non-Newtonian fluid
describing the rheology of the flow-like landslides, multi-phase flow
and free surface have to be included in the simulation. In this work, the
commercial computational fluid dynamics (CFD) package FLUENT is
used to model the surge waves due to flow-like landslides. The
comparison between the numerical results and experimental data
reported in the literature confirms the accuracy of the method.
Abstract: The power generated by eight point absorber type wave energy converters each having a different buoy are calculated in order to investigate the performances of buoys in this study. The calculations are carried out by modeling three different sea states observed in two different locations in the Black Sea. The floats analyzed in this study have two basic geometries and four different draft/radius (d/r) ratios. The buoys possess the shapes of a semi-ellipsoid and a semi-elliptic paraboloid. Additionally, the draft/radius ratios range from 0.25 to 1 by an increment of 0.25. The radiation forces acting on the buoys due to the oscillatory motions of these bodies are evaluated by employing a 3D panel method along with a distribution of 3D pulsating sources in frequency domain. On the other hand, the wave forces acting on the buoys which are taken as the sum of Froude-Krylov forces and diffraction forces are calculated by using linear wave theory. Furthermore, the wave energy converters are assumed to be taut-moored to the seabed so that the secondary body which houses a power take-off system oscillates with much smaller amplitudes compared to the buoy. As a result, it is assumed that there is not any significant contribution to the power generation from the motions of the housing body and the only contribution to power generation comes from the buoy. The power take-off systems of the wave energy converters are high pressure oil hydraulic systems which are identical in terms of their characteristic parameters. The results show that the power generated by wave energy converters which have semi-ellipsoid floats is higher than that of those which have semi elliptic paraboloid floats in both locations and in all sea states. It is also determined that the power generated by the wave energy converters follow an unsteady pattern such that they do not decrease or increase with changing draft/radius ratios of the floats. Although the highest power level is obtained with a semi-ellipsoid float which has a draft/radius ratio equal to 1, other floats of which the draft/radius ratio is 0.25 delivered higher power that the floats with a draft/radius ratio equal to 1 in some cases.
Abstract: The exploitation of flow pulsation in micro- and
mini-channels is a potentially useful technique for enhancing cooling
of high-end photonics and electronics systems. It is thought that
pulsation alters the thickness of the hydrodynamic and thermal
boundary layers, and hence affects the overall thermal resistance
of the heat sink. Although the fluid mechanics and heat transfer
are inextricably linked, it can be useful to decouple the parameters
to better understand the mechanisms underlying any heat transfer
enhancement. Using two-dimensional, two-component particle image
velocimetry, the current work intends to characterize the heat transfer
mechanisms in pulsating flow with a mean Reynolds number of
48 by experimentally quantifying the hydrodynamics of a generic
liquid-cooled channel geometry. Flows circulated through the test
section by a gear pump are modulated using a controller to achieve
sinusoidal flow pulsations with Womersley numbers of 7.45 and
2.36 and an amplitude ratio of 0.75. It is found that the transient
characteristics of the measured velocity profiles are dependent on the
speed of oscillation, in accordance with the analytical solution for
flow in a rectangular channel. A large velocity overshoot is observed
close to the wall at high frequencies, resulting from the interaction
of near-wall viscous stresses and inertial effects of the main fluid
body. The steep velocity gradients at the wall are indicative of
augmented heat transfer, although the local flow reversal may reduce
the upstream temperature difference in heat transfer applications.
While unsteady effects remain evident at the lower frequency, the
annular effect subsides and retreats from the wall. The shear rate at
the wall is increased during the accelerating half-cycle and decreased
during deceleration compared to steady flow, suggesting that the flow
may experience both enhanced and diminished heat transfer during
a single period. Hence, the thickness of the hydrodynamic boundary
layer is reduced for positively moving flow during one half of the
pulsation cycle at the investigated frequencies. It is expected that the
size of the thermal boundary layer is similarly reduced during the
cycle, leading to intervals of heat transfer enhancement.
Abstract: This paper presents the performance characteristics of
Darrieus-type vertical axis wind turbine (VAWT) with NACA airfoil
blades. The performance of Darrieus-type VAWT can be
characterized by torque and power. There are various parameters
affecting the performance such as chord length, helical angle, pitch
angle and rotor diameter. To estimate the optimum shape of Darrieustype
wind turbine in accordance with various design parameters, we
examined aerodynamic characteristics and separated flow occurring
in the vicinity of blade, interaction between flow and blade, and
torque and power characteristics derived from it. For flow analysis,
flow variations were investigated based on the unsteady RANS
(Reynolds-averaged Navier-Stokes) equation. Sliding mesh algorithm
was employed in order to consider rotational effect of blade. To
obtain more realistic results we conducted experiment and numerical
analysis at the same time for three-dimensional shape. In addition,
several parameters (chord length, rotor diameter, pitch angle, and
helical angle) were considered to find out optimum shape design and
characteristics of interaction with ambient flow. Since the NACA
airfoil used in this study showed significant changes in magnitude of
lift and drag depending on an angle of attack, the rotor with low drag,
long cord length and short diameter shows high power coefficient in
low tip speed ratio (TSR) range. On the contrary, in high TSR range,
drag becomes high. Hence, the short-chord and long-diameter rotor
produces high power coefficient. When a pitch angle at which airfoil
directs toward inside equals to -2° and helical angle equals to 0°,
Darrieus-type VAWT generates maximum power.
Abstract: This paper discusses the applicability of the numerical model for a damage prediction method of the accidental hydrogen explosion occurring in a hydrogen facility. The numerical model was based on an unstructured finite volume method (FVM) code “NuFD/FrontFlowRed”. For simulating unsteady turbulent combustion of leaked hydrogen gas, a combination of Large Eddy Simulation (LES) and a combustion model were used. The combustion model was based on a two scalar flamelet approach, where a G-equation model and a conserved scalar model expressed a propagation of premixed flame surface and a diffusion combustion process, respectively. For validation of this numerical model, we have simulated the previous two types of hydrogen explosion tests. One is open-space explosion test, and the source was a prismatic 5.27 m3 volume with 30% of hydrogen-air mixture. A reinforced concrete wall was set 4 m away from the front surface of the source. The source was ignited at the bottom center by a spark. The other is vented enclosure explosion test, and the chamber was 4.6 m × 4.6 m × 3.0 m with a vent opening on one side. Vent area of 5.4 m2 was used. Test was performed with ignition at the center of the wall opposite the vent. Hydrogen-air mixtures with hydrogen concentrations close to 18% vol. were used in the tests. The results from the numerical simulations are compared with the previous experimental data for the accuracy of the numerical model, and we have verified that the simulated overpressures and flame time-of-arrival data were in good agreement with the results of the previous two explosion tests.
Abstract: This paper outlines the development of an
experimental technique in quantifying supersonic jet flows, in an
attempt to avoid seeding particle problems frequently associated with
particle-image velocimetry (PIV) techniques at high Mach numbers.
Based on optical flow algorithms, the idea behind the technique
involves using high speed cameras to capture Schlieren images of the
supersonic jet shear layers, before they are subjected to an adapted
optical flow algorithm based on the Horn-Schnuck method to
determine the associated flow fields. The proposed method is capable
of offering full-field unsteady flow information with potentially
higher accuracy and resolution than existing point-measurements or
PIV techniques. Preliminary study via numerical simulations of a
circular de Laval jet nozzle successfully reveals flow and shock
structures typically associated with supersonic jet flows, which serve
as useful data for subsequent validation of the optical flow based
experimental results. For experimental technique, a Z-type Schlieren
setup is proposed with supersonic jet operated in cold mode,
stagnation pressure of 4 bar and exit Mach of 1.5. High-speed singleframe
or double-frame cameras are used to capture successive
Schlieren images. As implementation of optical flow technique to
supersonic flows remains rare, the current focus revolves around
methodology validation through synthetic images. The results of
validation test offers valuable insight into how the optical flow
algorithm can be further improved to improve robustness and
accuracy. Despite these challenges however, this supersonic flow
measurement technique may potentially offer a simpler way to
identify and quantify the fine spatial structures within the shock shear
layer.