Abstract: In this study, the dispersion of heavy particles line in
an isotropic and incompressible three-dimensional turbulent flow has
been studied using the Kinematic Simulation techniques to find out
the evolution of the line fractal dimension. The fractal dimension of
the line is found in the case of different particle gravity (in practice,
different values of particle drift velocity) in the presence of small
particle inertia with a comparison with that obtained in the diffusion
case of material line at the same Reynolds number. It can be
concluded for the dispersion of heavy particles line in turbulent flow
that the particle gravity affect the fractal dimension of the line for
different particle gravity velocities in the range 0.2 < W < 2. With
the increase of the particle drift velocity, the fractal dimension of the
line decreases which may be explained as the particles pass many
scales in their journey in the direction of the gravity and the particles
trajectories do not affect by these scales at high particle drift
velocities.
Abstract: Coronary artery bypass grafts (CABG) are widely
studied with respect to hemodynamic conditions which play
important role in presence of a restenosis. However, papers which
concern with constitutive modeling of CABG are lacking in the
literature. The purpose of this study is to find a constitutive model for
CABG tissue. A sample of the CABG obtained within an autopsy
underwent an inflation–extension test. Displacements were
recoredered by CCD cameras and subsequently evaluated by digital
image correlation. Pressure – radius and axial force – elongation
data were used to fit material model. The tissue was modeled as onelayered
composite reinforced by two families of helical fibers. The
material is assumed to be locally orthotropic, nonlinear,
incompressible and hyperelastic. Material parameters are estimated
for two strain energy functions (SEF). The first is classical
exponential. The second SEF is logarithmic which allows
interpretation by means of limiting (finite) strain extensibility.
Presented material parameters are estimated by optimization based
on radial and axial equilibrium equation in a thick-walled tube. Both
material models fit experimental data successfully. The exponential
model fits significantly better relationship between axial force and
axial strain than logarithmic one.
Abstract: A new numerical method for solving the twodimensional,
steady, incompressible, viscous flow equations on a
Curvilinear staggered grid is presented in this paper. The proposed
methodology is finite difference based, but essentially takes
advantage of the best features of two well-established numerical
formulations, the finite difference and finite volume methods. Some
weaknesses of the finite difference approach are removed by
exploiting the strengths of the finite volume method. In particular,
the issue of velocity-pressure coupling is dealt with in the proposed
finite difference formulation by developing a pressure correction
equation in a manner similar to the SIMPLE approach commonly
used in finite volume formulations. However, since this is purely a
finite difference formulation, numerical approximation of fluxes is
not required. Results obtained from the present method are based on
the first-order upwind scheme for the convective terms, but the
methodology can easily be modified to accommodate higher order
differencing schemes.
Abstract: The numerical simulation of the slip effect via
vicoelastic fluid for 4:1 contraction problem is investigated with
regard to kinematic behaviors of streamlines and stress tensor by
models of the Navier-Stokes and Oldroyd-B equations. Twodimensional
spatial reference system of incompressible creeping flow
with and without slip velocity is determined and the finite element
method of a semi-implicit Taylor-Galerkin pressure-correction is
applied to compute the problem of this Cartesian coordinate system
including the schemes of velocity gradient recovery method and the
streamline-Upwind / Petrov-Galerkin procedure. The slip effect at
channel wall is added to calculate after each time step in order to
intend the alteration of flow path. The result of stress values and the
vortices are reduced by the optimum slip coefficient of 0.1 with near
the outcome of analytical solution.
Abstract: A 3D simulation study for an incompressible
slip flow around a spherical aerosol particle was performed.
The full Navier-Stokes equations were solved and the velocity
jump at the gas-particle interface was treated numerically by
imposition of the slip boundary condition. Analytical solution
to the Stokesian slip flow past a spherical particle was used as
a benchmark for code verification, and excellent agreement
was achieved. The Simulation results showed that in addition
to the Knudsen number, the Reynolds number affects the slip
correction factor. Thus, the Cunningham-based slip corrections
must be augmented by the inclusion of the effect of
Reynolds number for application to Lagrangian tracking of
fine particles. A new expression for the slip correction factor
as a function of both Knudsen number and Reynolds number
was developed.
Abstract: In this study a two dimensional axisymmetric, steady state and incompressible laminar flow in a rotating single disk is numerically investigated. The finite volume method is used for solving the momentum equations. The numerical model and results
are validated by comparing it to previously reported experimental data for velocities, angles and moment coefficients. It is
demonstrated that increasing the axial distance increases the value of axial velocity and vice versa for tangential and total velocities. However, the maximum value of nondimensional radial velocity
occurs near the disk wall. It is also found that with increase rotational Reynolds number, moment coefficient decreases.
Abstract: A conjugate heat transfer for steady two-dimensional
mixed convection with magnetic hydrodynamic (MHD) flow of an
incompressible quiescent fluid over an unsteady thermal forming
stretching sheet has been studied. A parameter, M, which is used to
represent the dominance of the magnetic effect has been presented in
governing equations. The similar transformation and an implicit
finite-difference method have been used to analyze the present
problem. The numerical solutions of the flow velocity distributions,
temperature profiles, the wall unknown values of f''(0) and '(θ (0) for
calculating the heat transfer of the similar boundary-layer flow are
carried out as functions of the unsteadiness parameter (S), the Prandtl
number (Pr), the space-dependent parameter (A) and
temperature-dependent parameter (B) for heat source/sink and the
magnetic parameter (M). The effects of these parameters have also
discussed. At the results, it will produce greater heat transfer effect
with a larger Pr and M, S, A, B will reduce heat transfer effects. At
last, conjugate heat transfer for the free convection with a larger G has
a good heat transfer effect better than a smaller G=0.
Abstract: Pressure driven microscale gas flow-separation has
been investigated by solving the compressible Navier-Stokes (NS)
system of equations. A two dimensional explicit finite volume (FV)
compressible flow solver has been developed using modified
advection upwind splitting methods (AUSM+) with no-slip/first
order Maxwell-s velocity slip conditions to predict the flowseparation
behavior in microdimensions. The effects of scale-factor
of the flow geometry and gas species on the microscale gas flowseparation
have been studied in this work. The intensity of flowseparation
gets reduced with the decrease in scale of the flow
geometry. In reduced dimension, flow-separation may not at all be
present under similar flow conditions compared to the larger flow
geometry. The flow-separation patterns greatly depend on the
properties of the medium under similar flow conditions.
Abstract: The steady incompressible flow has been solved in cylindrical coordinates in both vapour region and wick structure. The governing equations in vapour region are continuity, Navier-Stokes and energy equations. These equations have been solved using SIMPLE algorithm. For study of parameters variation on heat pipe operation, a benchmark has been chosen and the effect of changing one parameter has been analyzed when the others have been fixed.
Abstract: In this paper, one-dimensional analysis of flow in a
single-stage gas gun is conducted. The compressible inviscid flow
equations are numerically solved by the second-order Roe TVD
method, by using moving boundaries. For investigation of real gas
effect the Noble-Able equation is applied. The numerical results are
compared with the experimental data to validate the numerical
scheme. The results show that with using the Noble-Able equation,
the muzzle velocity decreases.
Abstract: This article presents the boundary conditions for the problem of turbulent supersonic gas flow in a plane channel with a perpendicular injection jets. The non-reflection boundary conditions for direct modeling of compressible viscous gases are studied. A formulation using the NSCBC (Navier- Stocks characteristic boundary conditions) through boundaries is derived for the subsonic inflow and subsonic non-reflection outflow situations. Verification of the constructed algorithm of boundary conditions is carried out by solving a test problem of perpendicular sound of jets injection into a supersonic gas flow in a plane channel.
Abstract: The paper presents a one-dimensional transient
mathematical model of compressible thermal multi-component gas
mixture flows in pipes. The set of the mass, momentum and enthalpy
conservation equations for gas phase is solved. Thermo-physical
properties of multi-component gas mixture are calculated by solving
the Equation of State (EOS) model. The Soave-Redlich-Kwong
(SRK-EOS) model is chosen. Gas mixture viscosity is calculated on
the basis of the Lee-Gonzales-Eakin (LGE) correlation. Numerical
analysis on rapid decompression in conventional dry gases is
performed by using the proposed mathematical model. The model is
validated on measured values of the decompression wave speed in
dry natural gas mixtures. All predictions show excellent agreement
with the experimental data at high and low pressure. The presented
model predicts the decompression in dry natural gas mixtures much
better than GASDECOM and OLGA codes, which are the most
frequently-used codes in oil and gas pipeline transport service.
Abstract: In the present paper, a numerical investigation has
been carried out to classify and clarify the effects of paramount
parameters on turbulent impinging slot jets. The effects of nozzle-s
exit turbulent intensity, distance between nozzle and impinging plate
are studied at Reynolds number 5000 and 20000. In addition, the
effect of Mach number that is varied between 0.3-0.8 at a constant
Reynolds number 133000 is investigated to elucidate the effect of
compressibility in impinging jet upon a flat plate. The wall that is
located at the same level with nozzle-s exit confines the flow. A
compressible finite volume solver is implemented for simulation the
flow behavior. One equation Spalart-Allmaras turbulent model is
used to simulate turbulent flow at this study. Assessment of the
Spalart-Allmaras turbulent model at high nozzle to plate distance,
and giving enough insights to characterize the effect of Mach number
at high Reynolds number for the complex impinging jet flow are the
remarkable results of this study.
Abstract: In this paper, the local grid refinement is focused by
using a nested grid technique. The Cartesian grid numerical method is
developed for simulating unsteady, viscous, incompressible flows
with complex immersed boundaries. A finite volume method is used in
conjunction with a two-step fractional-step procedure. The key aspects
that need to be considered in developing such a nested grid solver are
imposition of interface conditions on the inter-block and accurate
discretization of the governing equation in cells that are with the
inter-block as a control surface. A new interpolation procedure is
presented which allows systematic development of a spatial
discretization scheme that preserves the spatial accuracy of the
underlying solver. The present nested grid method has been tested by
two numerical examples to examine its performance in the two
dimensional problems. The numerical examples include flow past a
circular cylinder symmetrically installed in a Channel and flow past
two circular cylinders with different diameters. From the numerical
experiments, the ability of the solver to simulate flows with
complicated immersed boundaries is demonstrated and the nested grid
approach can efficiently speed up the numerical solutions.
Abstract: Plasma Wind Tunnels (PWT) are extensively used for screening and qualification of re-entry Thermel Protection System (TPS) materials. Proper design of a supersonic diffuser for plasma wind tunnel is of importance for achieving good pressurerecovery (thereby reducing vacuum pumping requirement & run time costs) and isolating downstream stream fluctuations from propagating costs) and isolating downstream stream fluctuationnts the details of a rapid design methodology successfully employed for designing supersonic diffuser for high power (several megawatts)plasma wind tunnels and numerical performance analysis of a diffuser configuration designed for one megawatt power rated plasma wind tunnel(enthalpy ~ 30 MJ/kg) using FLUENT 6.3® solver for different diffuser operating sub-atmospheric back-pressures.
Abstract: Gas hydrates can agglomerate and block multiphase oil and gas pipelines when water is present at hydrate forming conditions. Using "Cold Flow Technology", the aim is to condition gas hydrates so that they can be transported as a slurry mixture without a risk of agglomeration. During the pipeline shut down however, hydrate particles may settle in bends and build hydrate plugs. An experimental setup has been designed and constructed to study the flow of such plugs at start up operations. Experiments have been performed using model fluid and model hydrate particles. The propagations of initial plugs in a bend were recorded with impedance probes along the pipe. The experimental results show a dispersion of the plug front. A peak in pressure drop was also recorded when the plugs were passing the bend. The evolutions of the plugs have been simulated by numerical integration of the incompressible mass balance equations, with an imposed mixture velocity. The slip between particles and carrier fluid has been calculated using a drag relation together with a particle-fluid force balance.
Abstract: An analysis is made of the flow of an incompressible viscoelastic fluid (of small memory) over a porous plate subject to suction or blowing. It is found that velocity at a point increases with increase in the elasticity in the fluid. It is also shown that wall shear stress depends only on suction and is also independent of the material of fluids. No steady solution for velocity distribution exists when there is blowing at the plate. Temperature distribution in the boundary layer is determined and it is found that temperature at a point decreases with increase in the elasticity in the fluid.
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 one-step conservative level set method, combined with a global mass correction method, is developed in this study to simulate the incompressible two-phase flows. The present framework do not need to solve the conservative level set scheme at two separated steps, and the global mass can be exactly conserved. The present method is then more efficient than two-step conservative level set scheme. The dispersion-relation-preserving schemes are utilized for the advection terms. The pressure Poisson equation solver is applied to GPU computation using the pCDR library developed by National Center for High-Performance Computing, Taiwan. The SMP parallelization is used to accelerate the rest of calculations. Three benchmark problems were done for the performance evaluation. Good agreements with the referenced solutions are demonstrated for all the investigated problems.
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.