Abstract: Due to energy and environment context, research is looking for the use of clean and energy efficient system in cooling industry. In this regard, the ejector represents one of the promising solutions. The thermal ejector is a passive component used for thermal compression in refrigeration and cooling systems, usually activated by heat either waste or solar. The present study introduces a theoretical analysis of the cooling system which uses a gas ejector thermal compression. A theoretical model is developed and applied for the design and simulation of the ejector, as well as the whole cooling system. Besides the conservation equations of mass, energy and momentum, the gas dynamic equations, state equations, isentropic relations as well as some appropriate assumptions are applied to simulate the flow and mixing in the ejector. This model coupled with the equations of the other components (condenser, evaporator, pump, and generator) is used to analyze profiles of pressure and velocity (Mach number), as well as evaluation of the cycle cooling capacity. A FORTRAN program is developed to carry out the investigation. Properties of refrigerant R134a are calculated using real gas equations. Among many parameters, it is thought that the generator pressure is the cornerstone in the cycle, and hence considered as the key parameter in this investigation. Results show that the generator pressure has a great effect on the ejector and on the whole cooling system. At high generator pressures, strong shock waves inside the ejector are created, which lead to significant condenser pressure at the ejector exit. Additionally, at higher generator pressures, the designed system can deliver cooling capacity for high condensing pressure (hot season).
Abstract: In this paper, we show shallow water in a tin box as an analogous simulation tool for high-speed aerodynamics education and research. It is customary that we use a water tank to create shallow water flow. While a flow in a water tank is not necessarily uniform and is sometimes wavy, we can visualize a clear supercritical flow even when we move a body manually in stationary water in a simple shallow tin box. We can visualize a blunt shock wave around a moving circular cylinder together with a shock pattern around a diamond airfoil. Another interesting analogous experiment is a hydrodynamic shock tube with water and tea. We observe the contact surface clearly due to color difference of the two liquids those are invisible in the real gas dynamics experiment. We first revisit the similarities between high-speed aerodynamics and shallow water hydraulics. Several educational and research experiments are then introduced for engineering students. Shallow water experiments in a tin box simulate properly the high-speed flows.
Abstract: Controlling the flow of fluids is a challenging problem
that arises in many fields. Burgers’ equation is a fundamental
equation for several flow phenomena such as traffic, shock waves,
and turbulence. The optimal feedback control method, so-called
model predictive control, has been proposed for Burgers’ equation.
However, the model predictive control method is inapplicable to
systems whose all state variables are not exactly known. In practical
point of view, it is unusual that all the state variables of systems are
exactly known, because the state variables of systems are measured
through output sensors and limited parts of them can be only
available. In fact, it is usual that flow velocities of fluid systems
cannot be measured for all spatial domains. Hence, any practical
feedback controller for fluid systems must incorporate some type of
state estimator. To apply the model predictive control to the fluid
systems described by Burgers’ equation, it is needed to establish
a state estimation method for Burgers’ equation with limited
measurable state variables. To this purpose, we apply unscented
Kalman filter for estimating the state variables of fluid systems
described by Burgers’ equation. The objective of this study is to
establish a state estimation method based on unscented Kalman filter
for Burgers’ equation. The effectiveness of the proposed method is
verified by numerical simulations.
Abstract: The aim of this work is to analyze a flow around the axisymmetric blunt body taken into account the chemical and vibrational nonequilibrium flow. This work concerns the entry of spacecraft in the atmosphere of the planet Mars. Since the equations involved are non-linear partial derivatives, the volume method is the only way to solve this problem. The choice of the mesh and the CFL is a condition for the convergence to have the stationary solution.
Abstract: The understanding of generation and collapse of acoustic cavitation bubbles are prerequisites for application of cavitation erosion. Microbubbles generated due to rapid fluctuation of pressure induced by propagation of ultrasonic wave lead to formation of high velocity microjets and or shock waves upon collapse. Due to vast application of ultrasonic, it is important to characterize and understand cavitation collapse pressure under the radiating surface at different conditions. A comparative investigation is carried out to determine impact load and dynamic pressure distribution exerted upon bubble collapse using thin film pressure sensors. Measurements were recorded at different input conditions such as amplitude, stand-off distance, insertion depth of the horn inside the liquid and pulse on-off time of acoustic vibrations. Impact force of 2.97 N is recorded at amplitude of 108 μm and stand-off distance of 1 mm from the sensor film, whereas impulsive force as low as 0.4 N is recorded at amplitude of 12 μm and stand-off distance of 5 mm from the sensor film. The results drawn from the investigation indicated that variety of impact loads can be achieved by controlling generation and collapse of bubbles, making it suitable to use for numerous application.
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: 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: In this study, failure analysis of pipe system at a micro
hydroelectric power plant is investigated. Failure occurred at the pipe
system in the powerhouse during shut down operation of the water
flow by a valve. This locking had caused a sudden shock wave, also
called “Water-hammer effect”, resulting in noise and inside pressure
increase. After visual investigation of the effect of the shock wave on
the system, a circumference crack was observed at the pipe flange
weld region. To establish the reason for crack formation, calculations
of pressure and stress values at pipe, flange and welding seams were
carried out and concluded that safety factor was high (2.2), indicating
that no faulty design existed. By further analysis, pipe system and
hydroelectric power plant was examined. After observations it is
determined that the plant did not include a ventilation nozzle (air
trap), that prevents the system of sudden pressure increase inside the
pipes which is caused by water-hammer effect. Analyses were carried
out to identify the influence of water-hammer effect on inside
pressure increase and it was concluded that, according Jowkowsky’s
equation, shut down time is effective on inside pressure increase. The
valve closing time was uncertain but by a shut down time of even one
minute, inside pressure would increase by 7.6 bar (working pressure
was 34.6 bar). Detailed investigations were also carried out on the
assembly of the pipe-flange system by considering technical
drawings. It was concluded that the pipe-flange system was not
installed according to the instructions. Two of five weld seams were
not applied and one weld was carried out faulty. This incorrect and
inadequate weld seams resulted in; insufficient connection of the pipe
to the flange constituting a strong notch effect at weld seam regions,
increase in stress values and the decrease of strength and safety
factor.
Abstract: Traditional mechanical control systems in thrust
vectoring are efficient in rocket thrust guidance but their costs
and their weights are excessive. The fluidic injection in the nozzle
divergent constitutes an alternative procedure to achieve the goal. In
this paper, we present a 3D analytical model for fluidic injection
in a supersonic nozzle integrating an orifice. The fluidic vectoring
uses a sonic secondary injection in the divergent. As a result, the
flow and interaction between the main and secondary jet has built in
order to express the pressure fields from which the forces and thrust
vectoring are deduced. Under various separation criteria, the present
analytical model results are compared with the existing numerical
and experimental data from the literature.
Abstract: Numerical studies have been carried out using a
validated two-dimensional standard k-omega turbulence model for
the design optimization of a thrust vector control system using shock
induced self-impinging supersonic secondary double jet. Parametric
analytical studies have been carried out at different secondary
injection locations to identifying the highest unsymmetrical
distribution of the main gas flow due to shock waves, which produces
a desirable side force more lucratively for vectoring. The results from
the parametric studies of the case on hand reveal that the shock
induced self-impinging supersonic secondary double jet is more
efficient in certain locations at the divergent region of a CD nozzle
than a case with supersonic single jet with same mass flow rate. We
observed that the best axial location of the self-impinging supersonic
secondary double jet nozzle with a given jet interaction angle, built-in
to a CD nozzle having area ratio 1.797, is 0.991 times the primary
nozzle throat diameter from the throat location. We also observed
that the flexible steering is possible after invoking ON/OFF facility to
the secondary nozzles for meeting the onboard mission requirements.
Through our case studies we concluded that the supersonic self-impinging
secondary double jet at predesigned jet interaction angle
and location can provide more flexible steering options facilitating
with 8.81% higher thrust vectoring efficiency than the conventional
supersonic single secondary jet without compromising the payload
capability of any supersonic aerospace vehicle.
Abstract: Hypersonic flows around spatial vehicles during their reentry phase in planetary atmospheres are characterized by intense aerothermodynamics phenomena. The aim of this work is to analyze high temperature flows around an axisymmetric blunt body taking into account chemical and vibrational non-equilibrium for air mixture species and the no slip condition at the wall. For this purpose, the Navier-Stokes equations system is resolved by the finite volume methodology to determine the flow parameters around the axisymmetric blunt body especially at the stagnation point and in the boundary layer along the wall of the blunt body. The code allows the capture of shock wave before a blunt body placed in hypersonic free stream. The numerical technique uses the Flux Vector Splitting method of Van Leer. CFL coefficient and mesh size level are selected to ensure the numerical convergence.
Abstract: Universal modeling method well proven for industrial
compressors was applied for design of the high flow rate supersonic
stage. Results were checked by ANSYS CFX and NUMECA Fine
Turbo calculations. The impeller appeared to be very effective at
transonic flow velocities. Stator elements efficiency is acceptable at
design Mach numbers too. Their loss coefficient versus inlet flow
angle performances correlates well with Universal modeling
prediction. The impeller demonstrates ability of satisfactory operation
at design flow rate. Supersonic flow behavior in the impeller inducer
at the shroud blade to blade surface Φ des deserves additional study.
Abstract: Numerical simulation performed to investigate the behavior of the high pressure hydrogen jetting of air. High pressure hydrogen (30–40 MPa) was injected to air at atmospheric pressure through 2mm orifice. Numerical simulations were performed with Kiva3V code with 2D axisymmetric geometry. Numerical simulations showed that auto ignition of high pressure hydrogen to air are possible due to molecular diffusion. Auto ignition was predicted at hydrogen-air contact surface due to mass and energy exchange between high temperature hydrogen and air heated by shock wave.
Abstract: This work presents a numerical simulation of the interaction of an incident shock wave propagates from the left to the right with a cone placed in a tube at shock. The Mathematical model is based on a non stationary, viscous and axisymmetric flow. The Discretization of the Navier-stokes equations is carried out by the finite volume method in the integral form along with the Flux Vector Splitting method of Van Leer. Here, adequate combination of time stepping parameter, CFL coefficient and mesh size level is selected to ensure numerical convergence. The numerical simulation considers a shock tube filled with air. The incident shock wave propagates to the right with a determined Mach number and crosses the cone by leaving behind it a stationary detached shock wave in front of the nose cone. This type of interaction is observed according to the time of flow.
Abstract: Flow field around hypersonic vehicles is very
complex and difficult to simulate. The boundary layers are squeezed
between shock layer and body surface. Resolution of boundary layer,
shock wave and turbulent regions where the flow field has high
values is difficult of capture. Detached eddy simulation (DES) is a
modification of a RANS model in which the model switches to a
subgrid scale formulation in regions fine enough for LES
calculations. Regions near solid body boundaries and where the
turbulent length scale is less than the maximum grid dimension are
assigned the RANS mode of solution. As the turbulent length scale
exceeds the grid dimension, the regions are solved using the LES
mode. Therefore the grid resolution is not as demanding as pure LES,
thereby considerably cutting down the cost of the computation. In
this research study hypersonic flow is simulated at Mach 8 and
different angle of attacks to resolve the proper boundary layers and
discontinuities. The flow is also simulated in the long wake regions.
Mesh is little different than RANS simulations and it is made dense
near the boundary layers and in the wake regions to resolve it
properly. Hypersonic blunt cone cylinder body with frustrum at angle
5o and 10 o are simulated and there aerodynamics study is performed
to calculate aerodynamics characteristics of different geometries. The
results and then compared with experimental as well as with some
turbulence model (SA Model). The results achieved with DES
simulation have very good resolution as well as have excellent
agreement with experimental and available data. Unsteady
simulations are performed for DES calculations by using duel time
stepping method or implicit time stepping. The simulations are
performed at Mach number 8 and angle of attack from 0o to 10o for
all these cases. The results and resolutions for DES model found
much better than SA turbulence model.
Abstract: When the shock front (SF) hits the central electrode
axis of plasma focus device, a reflected shock wave moves radially
outwards. The current sheath (CS) results from ionization of filled
gas between two electrodes continues to compress inwards until it
hits the out-going reflected shock front. In this paper the Lagrangian
equations are solved for a parabolic shock trajectory yielding a first
and second approximation for the CS path. To determine the
accuracy of the approximation, the same problem is solved for a
straight shock.
Abstract: This paper is devoted to predict laminar and turbulent
heating rates around blunt re-entry spacecraft at hypersonic
conditions. Heating calculation of a hypersonic body is normally
performed during the critical part of its flight trajectory. The
procedure is of an inverse method, where a shock wave is assumed,
and the body shape that supports this shock, as well as the flowfield
between the shock and body, are calculated. For simplicity the
normal momentum equation is replaced with a second order pressure
relation; this simplification significantly reduces computation time.
The geometries specified in this research, are parabola and ellipsoids
which may have conical after bodies. An excellent agreement is
observed between the results obtained in this paper and those
calculated by others- research. Since this method is much faster than
Navier-Stokes solutions, it can be used in preliminary design,
parametric study of hypersonic vehicles.
Abstract: The paper presents the results of theoretical and
numerical modeling of propagation of shock waves in bubbly liquids
related to nonlinear effects (realistic equation of state, chemical
reactions, two-dimensional effects). On the basis on the Rankine-
Hugoniot equations the problem of determination of parameters of
passing and reflected shock waves in gas-liquid medium for
isothermal, adiabatic and shock compression of the gas component is
solved by using the wide-range equation of state of water in the
analitic form. The phenomenon of shock wave intensification is
investigated in the channel of variable cross section for the
propagation of a shock wave in the liquid filled with bubbles
containing chemically active gases. The results of modeling of the
wave impulse impact on the solid wall covered with bubble layer are
presented.
Abstract: In this work, we try to find the best setting
of Computational Fluid Dynamic solver available for the problems in
the field of supersonic internal flows. We used the supersonic air-toair
ejector to represent the typical problem in focus. There are
multiple oblique shock waves, shear layers, boundary layers
and normal shock interacting in the supersonic ejector making this
device typical in field of supersonic inner flows. Modeling of shocks
in general is demanding on the physical model of fluid, because
ordinary conservation equation does not conform to real conditions in
the near-shock region as found in many works. From these reasons,
we decided to take special care about solver setting in this article by
means of experimental approach of color Schlieren pictures and
pneumatic measurement. Fast pressure transducers were used to
measure unsteady static pressure in regimes with normal shock in
mixing chamber. Physical behavior of ejector in several regimes is
discussed. Best choice of eddy-viscosity setting is discussed on the
theoretical base. The final verification of the k-ω SST is done on the
base of comparison between experiment and numerical results.
Abstract: Dynamics of laser radiation – metal target interaction
in water at 1064 nm by applying Mach-Zehnder interference
technique was studied. The mechanism of generating the well
developed regime of evaporation of a metal surface and a spherical
shock wave in water is proposed. Critical intensities of the NIR for
the well developed evaporation of silver and gold targets were
determined. Dynamics of shock waves was investigated for earlier
(dozens) and later (hundreds) nanoseconds of time. Transparent
expanding plasma-vapor-compressed water object was visualized and
measured. The thickness of compressed layer of water and pressures
behind the front of a shock wave for later time delays were obtained
from the optical treatment of interferograms.