Study on the Fuel Air Mixing Induced by a Shock Wave Propagating into a H2-Air Interface

2nd-order upwind TVD scheme was used to solve the laminar, fully Navier-Stokes equations. The numerical simulations were done on the propagation of a shock wave with Ma _S = 2 and 4 into a hydrogen and air mixture in a duct and a duct with a rearward step. The results indicate that a swirling vortex may be generated in the lopsided interface behind the moving shock. Meanwhile, the complex shock system is also formed in this shear flow region. A large swirling vortex is produced and the fuel mixing can be enhanced by a shock wave at low Mach number. But in a duct with a rearward step, the shock almost disappears in hydrogen for Ma _S = 2. The shock in hydrogen will become strong if Ma _S is large. Similar to the condition of a shock moving in a duct full of hydrogen and air, a large vortex can be formed in the shear flow region. The large swirling vortex even gets through the reflected shock and impacts on the lower wall. Then, the distribution of hydrogen behind the rearward step is divided into two regions. The transition from regular reflection to Mach reflection was observed as well in case Ma _S = 4.     

STUDY ON THE FUEL AIR MIXING INDUCED BY A SHOCK WAVE PROPAGATING INTO A H_2-AIR INTERFACE

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｈｅａｒｌａｙｅｒｆｌｏｗｉｓａｎｉｍｐｏｒｔａｎｔｐｒｏｂｌｅｍｉｎｓｔｕｄｙｉｎｇｃｏｍｐｒｅｓｓｉｂｌｅｆｌｏｗａｎｄｓｕｐｅｒｓｏｎｉｃｃｏｍｂｕｓｔｉｏｎ ,ｓｕｃｈａｓｆｕｅｌｄｉｆｆｕｓｉｏｎａｎｄｍｉｘｉｎｇｉｎａＳｃｒａｍｊｅｔ[1].Ｒｏｓｈｋｏ[2 ]ｓｔｕｄｉｅｄｅｘｐｅｒｉｍｅｎｔａｌｌｙｔｈｅｐｈｅｎｏｍｅｎｏｎｏｆｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｕｂｓｏｎｉｃｓｈｅａｒｆｌｏｗａｎｄｆｏｕｎｄｌａｒｇｅ ,ｃｏｈ…     

Analysis of Flow and Heat Transfer in a Vertical Rectangular Duct Using a Non-Darcy Model

Numerical investigation of steady natural convection flow through a fluid-saturated porous medium in a vertical rectangular duct is investigated. The Darcy-Forchheimer-Brinkman model is used to represent the fluid transport within the porous medium. One of the vertical walls of the duct is cooled to a constant temperature, while the other wall is heated to constant but different temperature. The other two sides of the duct are insulated. The finite difference method of second-order accuracy is used to solve the non-dimensional governing equations. The results are presented graphically to show the effects of the Darcy number, inertial parameter, Grashof number, Brinkman number, aspect ratio, and viscosity ratio. It is found that an increase in the Darcy number and inertial parameter reduces the flow intensity whereas an increase in the Grashof number, Brinkman number, aspect ratio, and viscosity ratio increases the flow intensity.     

2D numerical simulation of passive autocatalytic recombiner for hydrogen mitigation

Resolving hydrogen related safety issues, pertaining to nuclear reactor safety has been an important area of research world over for the past decade. The studies on hydrogen transport behavior and development of hydrogen mitigation systems are still being pursued actively in various research labs, including Bhabha Atomic Research Centre (BARC), in India. The passive autocatalytic recombiner (PAR) is one of such hydrogen mitigating device consisting of catalyst surfaces arranged in an open-ended enclosure. In the plate type recombiner design sheets made of stainless steel and coated with platinum catalyst material are arranged in parallel inside a flow channel. The catalyst elements are exposed to a constant flow of a mixture of air, hydrogen and steam, a catalytic reaction occurs spontaneously at the catalyst surfaces and the heat of reaction produces natural convection flow through the enclosure. Numerical simulation and experiments are required for an in-depth knowledge of such plate type PAR. Specific finite volume based in-house 2D computational fluid dynamics (CFD) code has been developed to model and analyse the working of these recombiners and has been used to simulate one literature quoted experiment. The validation results were in good agreement against literature quoted German REKO experiments. Parametric study has been performed for particular recombiner geometry for various inlet conditions. Salient features of the simplified CFD model developed at BARC and results of the present model calculations are presented in this paper.     

Combined forced and free flow in a vertical rectangular duct with prescribed wall heat flux

In this paper, combined forced and free convection is studied in a vertical rectangular duct with a prescribed uniform wall heat flux (H2 boundary condition). A different heat flux value for each plane wall is considered; the condition of a uniform wall heat flux throughout the duct results as a special case. The local momentum and energy balance equations are written in a dimensionless form and solved numerically, by means of a Galerkin finite element method. The numerical solution gives the dimensionless velocity and temperature distributions, together with the values of the Fanning friction factor, of the Nusselt number, of the momentum flux correction factor and of the kinetic energy correction factor. These dimensionless parameters are reported as functions of the aspect ratio and of the ratio between the Grashof number, Gr, and the Reynolds number, Re. The threshold values of Gr/Re for the onset of flow reversal are evaluated.     

Detonation combustion of hydrogen in an axisymmetric channel with a central body

The problem of initiation and stabilization of detonation combustion of a hydrogen–air mixture injected into an axisymmetric channel with a finite-length central body in a flow with a Mach number M₀ = 5–9 is solved. It is numerically demonstrated that the presence of the central body both in a convergent–divergent nozzle and in an expanding channel leads to stabilization of detonation combustion of a stoichiometric hydrogen–air mixture at free-stream Mach numbers M₀ > 7. Various channel configurations that ensure different values of thrust generated by detonation combustion of a stoichiometric hydrogen–air mixture are compared.     

Simultaneous heat and mass transfer inside a vertical tube in evaporating a heated falling alcohols liquid film into a stream of dry air

A numerical study of the evaporation in mixed convection of a pure alcohol liquid film: ethanol and methanol was investigated. It is a turbulent liquid film falling on the internal face of a vertical tube. A laminar flow of dry air enters the vertical tube at constant temperature in the downward direction. The wall of the tube is subjected to a constant and uniform heat flux. The model solves the coupled parabolic governing equations in both phases including turbulent liquid film together with the boundary and interfacial conditions. The systems of equations obtained by using an implicit finite difference method are solved by TDMA method. A Van Driest model is adopted to simulate the turbulent liquid film flow. The influence of the inlet liquid flow, Reynolds number in the gas flow and the wall heat flux on the intensity of heat and mass transfers are examined. A comparison between the results obtained for studied alcohols and water in the same conditions is made.     

Aeroacoustical coupling in a ducted shallow cavity and fluid/structure effects on a steam line

A pure tone phenomenon has been observed at 460 Hz in a piping steam line. The acoustical energy has been identified to be generated in an open gate valve and to be of cavity noise type. This energy is then transmitted to the main pipe by fluid/structure coupling. The objectives here are to display the mechanism of the flow acoustic coupling in the cavity and in the duct through an aeroacoustical analysis and to understand the way of energy transfer from the fluid to the main pipe through a vibroacoustical analysis. Concerning the first objective, an experimental study by means of 2/7 scale models in air is analysed by means of numerical flow simulation. The flow acoustic phenomena are modelled by computing the Euler equations. Two different computations are carried out: in the first one, a pure Euler modelling is used, in the second one, a boundary layer obtained from experimental data is introduced in the computation in order to have a realistic flow profile upstream the cavity. The boundary layer flow profile appears to be essential to recover the experimentally observed coupling between the shear-layer instability and the acoustical transverse mode of the pipe. The numerical results confirm that the second aerodynamic mode is responsible for the oscillation. While the predicted frequency agrees about 1% with the scale model experiments, the predicted amplitude is approximately 15 dB too low. For the second objective, fluid/structure coupling in the main pipe is studied using two fully coupled methods. The first method consists in a modal analysis of the line using a fluid–structure finite element model. The second one is based on the analysis of dispersion diagrams derived from the local equations of cylindrical shells filled with fluid. The way of energy transfer in transverse acoustical waves coupled with flexion-ovalization deformations of the pipe is highlighted using both methods. The dispersion diagrams allow a fast and accurate analysis. The modal analysis using a finite-element model may complete the first one with quantitative data. The link between the fluid/acoustic and the fluid/structure analysis is then the excitation of the transverse acoustical mode of the duct.     

MHD non-Darcian flow through a non-isothermal vertical surface embedded in a porous medium with radiation

The effect of MHD on steady two-dimensional laminar mixed flow about a vertical porous surface is numerically analyzed. Also the effects of radiation and heat generation and absorption are considered. A power law variation of temperature along the vertical wall is assumed. The nonlinear boundary-layer equations were transformed and the resulting differential equations were solved by an implicit finite difference scheme (Keller box method). Numerical results for the velocity distribution and the temperature distribution are presented for various values of Prandtl number Pr, magnetic parameter, porous medium parameter and internal heat generation or absorption coefficient. Further validation with previous works is carried out.     

Two-dimensional problem of the impact of a vertical wall on a layer of a partially aerated liquid

The two-dimensional unsteady problem of the impact of a vertical wall on a layer of a liquid which is mixed with air near the wall and does not contain air bubbles away from the wall is solved in a linear approximation. The gas-liquid mixture is modeled by a homogeneous, ideal, and weakly compressible medium with a reduced sound velocity dependent on the air concentration in the gas-liquid mixture. Outside the gas-liquid layer, the liquid is considered ideal and incompressible. During the initial stage of the impact, the liquid flow and the hydrodynamic pressure are determined using the linear theory of the potential motion of an inhomogeneous liquid. The dependence of the amplitude of the impact pressure along the wall on the air concentration in the gas-liquid layer and on the thickness of this layer is investigated. For a small relative thickness of the layer, the thin-layer approximation is used. It is shown that the solution of the original problem tends to the approximate solution as the thickness of the layer decreases. It is shown that the presence of the gas-liquid layer leads to wall pressure oscillations. Estimates are obtained for the pressure amplitude and the oscillation period. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 34–46, September–October, 2006.     

Finite element analysis of air flow around permeable sand fences

A numerical study of the turbulent air flow in a trench trap and the turbulent flow around a permeable sand fence is reported in this paper. The two-dimensional modified k–ε turbulence model proposed by Kato and Launder is used to predict the turbulent characteristics of the air flow. The discretization method for the governing equations is the three-step Taylor/Galerkin finite element method proposed by the authors. For the flow in a trench trap the numerical results are compared with experimental data obtained under realistic conditions using a large wind tunnel. For the air flow around a permeable sand fence a pressure loss model is used to represent the effect of the porosity of the fence on the flow field. © 1997 John Wiley & Sons, Ltd.     

Numerical study of condensing a small concentration of vapour inside a vertical tube

The purpose of this study is to analyse the combined heat and mass transfer of liquid film condensation from a small steam–air mixtures flowing downward along a vertical tube. Both liquid and gas stream are approached by two coupled laminar boundary layer. An implicit finite difference method is employed to solve the coupled governing equations for liquid film and gas flow together with the interfacial matching conditions. The effects of a wide range of changes of three independent variables (inlet pressure, inlet Reynolds number and wall temperature) on the concentration at exit tube, local Nusselt and Sherwood numbers, film thickness, accumulated condensate rate and temperature are carefully examined. The numerical results indicate that in the case of condensing a small concentration of vapours from a mixture, the resistance to heat and mass transfer by non-condensable gas becomes very intense. The comparisons of average Nusselt number and local condensate heat transfer coefficient with the literature results are in good agreement.     

Fully developed magneto convection flow in a vertical rectangular duct

An analysis is performed to study the MHD free convection flow in a vertical rectangular duct for laminar and fully developed regime taking into consideration the effects of Ohmic heating and viscous dissipation. Numerical solutions are found using finite difference method of second-order accuracy. The effects of various physical parameters such as Hartmann number, aspect ratio, buoyancy parameter and circuit parameter are presented graphically. It is found that as Hartmann number, buoyancy parameter and aspect ratio increase, the upward and downward flow rates are increased for open circuit but decrease for short circuit.     

Mixed convection in a vertical channel with discrete heat sources at the wall

The mixed convection in a vertical plane-parallel channel with two heat sources of finite dimensions located at the wall is analyzed on the basis of a two-dimensional numerical simulation. The effect of the distance between the heat sources on the flow pattern and the temperature field is studied. Calculations are performed on the Grashof and Reynolds number ranges from 0–10⁵ and 0–10, respectively, at a Prandtl number of 0.7. The mathematical model is based on the time-dependent Navier-Stokes equations in the Boussinesq approximation. The problem is solved by the finite element method.     

Time-dependent liquid metal flows with free convection and a deformable free surface

The finite element method is employed to investigate time-dependent liquid metal flows with free convection, free surfaces and Marangoni effects. The liquid circulates in a two-dimensional shallow trough with differentially heated vertical walls. The spatial formulation incorporates mixed Lagrangian approximations to the velocity, pressure, temperature and free surface position. The time integration is performed with the backward Euler and trapezoid rule methods with step size control. The Galerkin method is used to reduce the problem to a set of non-linear equations which are solved with the Newton–Raphson method. Calculations are performed for conditions relevant to the electron beam vaporization of refractory metals. The Prandtl number is 0·015 and Grashof number are in the transition range between laminar and turbulent flow. The results reveal the effects of flow intensity, surface tension gradients, mesh refinement and time integration strategy.     

Dispersed multiphase flow with air-driven runback of a liquid layer at a moving boundary

Multiphase flow with impinging droplets on an icing surface with a flowing supercooled surface layer is investigated. The air-assisted flowing layer is modelled with the cross-phase shear stresses imparted at the moving liquid/air interface. Runback and runoff of the surface layer are predicted by mass flow across the boundaries between adjacent elements in the numerical formulation. This liquid runoff is determined by coupled heat and momentum balances for the unfrozen water layer. The numerical analysis is developed with a control-volume-based finite element method (CVFEM). An Eulerian formulation with volume averaging is developed to accommodate the near-wall elements containing both dispersed and continuous phases. The predicted results are successfully validated through comparisons with analytical solutions and measured data.     

Upwind adaptive finite element investigations of the two-dimensional reactive interaction of supersonic gaseous jets

We describe an upwind finite element method aimed at numerically simulating the two-dimensional transonic flow of a reactive gaseous mixture. The method uses in particular a triangular finite element mesh, with an adaptive procedure based on mesh refinement by triangle division, and an upwind non-oscillatory scheme based on an approximate Riemann solver for the evaluation of the convective terms for all species. Results concerning the reactive interaction of two supersonic gaseous jets are presented.     

Synthesis of diamond structures from the jet of the H<Subscript>2</Subscript> + CH<Subscript>4</Subscript> mixture in a cocurrent axisymmetric hydrogen flow

The flow of a hydrogen–methane mixture through heated coaxial cylindrical tungsten channels with a built-in tungsten wire is studied by the Direct Simulation Monte Carlo method. The purpose of the study is further development of the gas-phase method of deposition of diamond structures. The axial distributions of the concentrations of the components of the hydrogen–carbon mixture are calculated by means of solving a system of chemical kinetics equations. A series of experiments on deposition of diamond structures from various flows of the hydrogen–methane mixture is performed. The calculated results are compared with the experimental data. Based on these comparisons, it is concluded that numerical optimization of operation modes of gas-dynamic reactors can be used for deposition of diamond structures.     

Numerical study on the two-phase flow distribution in a T-junction

The main objective of this paper is to investigate the ability of a two-dimensional two-fluid computer code to predict the phase separation in a T-junction. A new semi-implicit numerical scheme is developed for solving the two-fluid model equations. Special attention is directed to the modelling of the constitutive for the interfacial friction term. Detailed distribution of void fraction, pressure and velocities are obtained for an air–water mixture in a vertical tee. Good agreement was obtained between the computer code results and the experimental data for the phase separation in the T-junction.