Abstract: In this paper, the differential quadrature method is applied to simulate natural convection in an inclined cubic cavity using velocity-vorticity formulation. The numerical capability of the present algorithm is demonstrated by application to natural convection in an inclined cubic cavity. The velocity Poisson equations, the vorticity transport equations and the energy equation are all solved as a coupled system of equations for the seven field variables consisting of three velocities, three vorticities and temperature. The coupled equations are simultaneously solved by imposing the vorticity definition at boundary without requiring the explicit specification of the vorticity boundary conditions. Test results obtained for an inclined cubic cavity with different angle of inclinations for Rayleigh number equal to 103, 104, 105 and 106 indicate that the present coupled solution algorithm could predict the benchmark results for temperature and flow fields. Thus, it is convinced that the present formulation is capable of solving coupled Navier-Stokes equations effectively and accurately.
Abstract: A parallel computational fluid dynamics code has been
developed for the study of aerodynamic heating problem in hypersonic
flows. The code employs the 3D Navier-Stokes equations as the basic
governing equations to simulate the laminar hypersonic flow. The cell
centered finite volume method based on structured grid is applied for
spatial discretization. The AUSMPW+ scheme is used for the inviscid
fluxes, and the MUSCL approach is used for higher order spatial
accuracy. The implicit LU-SGS scheme is applied for time integration
to accelerate the convergence of computations in steady flows. A
parallel programming method based on MPI is employed to shorten
the computing time. The validity of the code is demonstrated by
comparing the numerical calculation result with the experimental data
of a hypersonic flow field around a blunt body.
Abstract: A general purpose viscous flow solver Ansys CFX
was used to solve the unsteady three-dimensional (3D) Reynolds
Averaged Navier-Stokes Equation (RANSE) for simulating a 3D
numerical viscous wave tank. A flap-type wave generator was
incorporated in the computational domain to generate the desired
incident waves. Authors have made effort to study the physical
behaviors of Flap type wave maker with governing parameters.
Dependency of the water fill depth, Time period of oscillations and
amplitude of oscillations of flap were studied. Effort has been made
to establish relations between parameters. A validation study was
also carried out against CFD methodology with wave maker theory.
It has been observed that CFD results are in good agreement with
theoretical results. Beaches of different slopes were introduced to
damp the wave, so that it should not cause any reflection from
boundary. As a conclusion this methodology can simulate the
experimental wave-maker for regular wave generation for different
wave length and amplitudes.
Abstract: In this study, the numerical solution of unsteady flow
between two concentric rotating spheres with suction and blowing at
their boundaries is presented. The spheres are rotating about a
common axis of rotation while their angular velocities are constant.
The Navier-Stokes equations are solved by employing the finite
difference method and implicit scheme. The resulting flow patterns
are presented for various values of the flow parameters including
rotational Reynolds number Re , and a blowing/suction Reynolds
number Rew . Viscous torques at the inner and the outer spheres are
calculated, too. It is seen that increasing the amount of suction and
blowing decrease the size of eddies generated in the annulus.
Abstract: This research paper presents the CFD analysis of
oscillating airfoil during pitch cycle. Unsteady subsonic flow is
simulated for pitching airfoil at Mach number 0.283 and Reynolds
number 3.45 millions. Turbulent effects are also considered for this
study by using K-ω SST turbulent model. Two-dimensional unsteady
compressible Navier-Stokes code including two-equation turbulence
model and PISO pressure velocity coupling is used. Pressure based
implicit solver with first order implicit unsteady formulation is used.
The simulated pitch cycle results are compared with the available
experimental data. The results have a good agreement with the
experimental data. Aerodynamic characteristics during pitch cycles
have been studied and validated.
Abstract: A numerical study on the influence of electroosmotic flow on analyte preconcentration by isotachophoresis ( ITP) is made. We consider that the double layer induced electroosmotic flow ( EOF) counterbalance the electrophoretic velocity and a stationary ITP stacked zones results. We solve the Navier-Stokes equations coupled with the Nernst-Planck equations to determine the local convective velocity and the preconcentration dynamics of ions. Our numerical algorithm is based on a finite volume method along with a secondorder upwind scheme. The present numerical algorithm can capture the the sharp boundaries of step-changes ( plateau mode) or zones of steep gradients ( peak mode) accurately. The convection of ions due to EOF reduces the resolution of the ITP transition zones and produces a dispersion in analyte zones. The role of the electrokinetic parameters which induces dispersion is analyzed. A one-dimensional model for the area-averaged concentrations based on the Taylor-Aristype effective diffusivity is found to be in good agreement with the computed solutions.
Abstract: The purpose of this paper is to elucidate the flow unsteady behavior for moving plug in convergent-divergent variable thrust nozzle. Compressible axisymmetric Navier-Stokes equations are used to study this physical phenomenon. Different velocities are set for plug to investigate the effect of plug movement on flow unsteadiness. Variation of mass flow rate and thrust are compared under two conditions: First, the plug is placed at different positions and flow is simulated to reach the steady state (quasi steady simulation) and second, the plug is moved with assigned velocity and flow simulation is coupled with plug movement (unsteady simulation). If plug speed is high enough and its movement time scale is at the same order of the flow time scale, variation of the mass flow rate and thrust level versus plug position demonstrate a vital discrepancy under the quasi steady and unsteady conditions. This phenomenon should be considered especially from response time viewpoints in thrusters design.
Abstract: In the present study, the surface temperature history of the adaptor part in a two-stage supersonic launch vehicle is accurately predicted. The full Navier-Stokes equations are used to estimate the aerodynamic heat flux and the one-dimensional heat conduction in solid phase is used to compute the temperature history. The instantaneous surface temperature is used to improve the applied heat flux, to improve the accuracy of the results.
Abstract: Different variants for buoyancy-affected terms in k-ε turbulence model have been utilized to predict the flow parameters more accurately, and investigate applicability of alternative k-ε turbulence buoyant closures in numerical simulation of a horizontal gravity current. The additional non-isotropic turbulent stress due to buoyancy has been considered in production term, based on Algebraic Stress Model (ASM). In order to account for turbulent scalar fluxes, general gradient diffusion hypothesis has been used along with Boussinesq gradient diffusion hypothesis with a variable turbulent Schmidt number and additional empirical constant c3ε.To simulate buoyant flow domain a 2D vertical numerical model (WISE, Width Integrated Stratified Environments), based on Reynolds- Averaged Navier-Stokes (RANS) equations, has been deployed and the model has been further developed for different k-ε turbulence closures. Results are compared against measured laboratory values of a saline gravity current to explore the efficient turbulence model.
Abstract: The exact solutions of the equations describing the steady plane motion of an incompressible fluid of variable viscosity for an arbitrary state equation are determined in the (ξ,ψ) − or (η,ψ )- coordinates where ψ(x,y) is the stream function, ξ and η are the parts of the analytic function, ϖ =ξ( x,y )+iη( x,y ). Most of the solutions involve arbitrary function/ functions indicating
that the flow equations possess an infinite set of solutions.
Abstract: The purpose of this study was to explore the complex
flow structure a novel active-type micromixer that based on concept of
Wankle-type rotor. The characteristics of this micromixer are two
folds; a rapid mixing of reagents in a limited space due to the
generation of multiple vortices and a graduate increment in dynamic
pressure as the mixed reagents is delivered to the output ports.
Present micro-mixer is consisted of a rotor with shape of triangle
column, a blending chamber and several inlet and outlet ports. The
geometry of blending chamber is designed to make the rotor can be
freely internal rotated with a constant eccentricity ratio. When the
shape of the blending chamber and the rotor are fixed, the effects of
rotating speed of rotor and the relative locations of ports on the mixing
efficiency are numerical studied. The governing equations are
unsteady, two-dimensional incompressible Navier-Stokes equation
and the working fluid is the water. The species concentration equation
is also solved to reveal the mass transfer process of reagents in various
regions then to evaluate the mixing efficiency.
The dynamic mesh technique was implemented to model the
dynamic volume shrinkage and expansion of three individual
sub-regions of blending chamber when the rotor conducted a complete
rotating cycle. Six types of ports configuration on the mixing
efficiency are considered in a range of Reynolds number from 10 to
300. The rapid mixing process was accomplished with the multiple
vortex structures within a tiny space due to the equilibrium of shear
force, viscous force and inertial force. Results showed that the highest
mixing efficiency could be attained in the following conditions: two
inlet and two outlet ports configuration, that is an included angle of 60
degrees between two inlets and an included angle of 120 degrees
between inlet and outlet ports when Re=10.