Computational Analysis of Flow through a Multiple Nozzle Driven Laser Cavity and Diffuser

Computations have been carried out for simulating supersonic flow through a set of converging-diverging nozzles with their expanding jets forming a laser cavity and flow patterns through diffusers, past the cavity. A thorough numerical investigation with 3-D RANS code is carried out to capture the flow distribution which comprises of shock patterns and multiple supersonic jet interactions. The analysis of pressure recovery characteristics during the flow through the diffusers is an important parameter of the simulation and is critical for the performance of the laser device. The results of the computation have shown a close agreement with the experimentally measured parameters as well as other established results indicating that the flow analysis done is found to be satisfactory.

Gibbs energy of formation of nickel orthosilicate (Ni2SiO4)

Two solid state galvanic cells:Pt, Ni + Ni2Si04 + Si02/(Y203)Zr02/Ni + + NiO, Pt (1) and Pt, Ni + NizSiOj + Si02/CaF2/Ni + + NiO, Pt (11) have been employed for the determination of the Gibbs' energy of formation of nickel orthosilicate(Ni2Si04) from nickel oxide and quartz. The emf of cell (I) was reversible and reproducible in the temperature range 925 to 1375K whereas emf of cell (11) drifted with time and changed polarity. From the results of cell (I), the Gibbs' energy of formation of nickel silicate is obtained as,2Ni0 (r.s.) + Si02 (quartz) + Ni2Si04 (olivine)Gibbs' energy of formation of the spinel form of Ni2Si04 is obtained by combining the data for olivine obtained in this study with high pressure data on olivine to spinel transition reported in the literature. The complex time dependence of the emf of cell (11) can be rationalised on the basis of formation of calcium silicates from calcium oxide, generally present as an impurity in the calcium fluoride electrolyte, and silica. The emf of cell (11) is shown to be the function of the activity of calcium oxide at the electrolyte/ electrode interface. The results provide strong evidence against the recent suggestion of mixed anionic conduction in calcium fluoride.

Solubility and activity of oxygen in liquid nickel in equilibrium with -Al2O3 and NiO.(1+x)Al2O3

The limiting solubility of oxygen in liquid nickel in equilibrium withα-alumina and nickel aluminate has been measured by inert gas fusion analysis of suction samples in the temperature range 1730 to 1975 K. The corresponding oxygen potential has been monitored by a solid electrolyte cell consisting of calcia stabilized zirconia as the electrolyte and Mo + MoO2 as the reference electrode. The results can be summarized by the following equations: log(at. pct O) = \frac - 10,005T + 4.944 ( ±0.015)log(atpctO)=T−10005+4944(0015) % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn DmO2 /4.606RT = log P O2 1/2 = \frac - 13,550T + 4.411 ( ±0.009)O24606RT=logPO212=T−13550+4411(0009) From simultaneous measurements of the potential and concentration of oxygen in melts, not in thermodynamic equilibrium with alumina and aluminate phases, information on the composition dependence of the activity coefficient and the standard free energy of solution of oxygen is obtained. For the reaction, $\frac{1}{2} O_2 \to \underset{\raise0.3em\hbox{$Missing close brace ΔG o = -72,930 - 7.11T (±840) J gr.at.–1 = + 0.216 at. pct OlogfO=T−500+0216atpctO where the standard state for dissolved oxygen is that which makes the value of activity equal to the concentration (in at. pct) in the limit as concentration approaches zero. The oxygen solubility in liquid nickel in equilibrium with solid NiO, evaluated from thermodynamic data, is compared with information reported in the literature. Implications of the results to the deoxidation equilibria of aluminum in nickel are discussed.

Kinetics of Sulphation of Nickel and Cobalt Ferrites with Na$O, Melt Containing F'e3+ Ions

It has been experimentally established that nickel and cobalt can be extracted from their ferrites with sodium sulphate melt containing femc ions. The kinetics of extraction from synthetic ferrites using a melt of sodium and ferric sulphates of eutectic composition has been studied as a function of the particle size of the ferrite and temperature in the range 900 to 1073 K. The divalent ions in the ferrite exchange with the ferric ion in the melt, leaving a residue of hematite.The rate of reaction conforms to the Crank-Ginstling-Brounshtein diffusion model. The reaction rate is governed by the counter-diffusion of ~ e an~d ~+i ' +(or co2+) ions in the hematite lattice. Analytical expressions for the rate constants have been derived from the experimental data as a function of particle size and temperature. The activation energy for the extraction of nickel from nickel ferrite is 154(+10) kJ mol-' and the corresponding value for cobalt is 142(+10)kJ mol;'. In sulphation roasting of minerals containing nickel, the yield of nickel is generally limited to 75% due to the formation of insoluble ferrites. The use of melts based on sodium sulphate provides a possible route for enhancing the recovery of nickel to approximately 98%.

A new matched thermochemical diagram for vacuum refining of nickel and cupronickel alloys

Sulfur and oxygen dissolved in nickel and cupronickel melts can be remwed as gaseous oxides of sulfur by a vacuum treatment. Presented in this paper is a new matched thermcxhemical disgran~ that permit.. direct evaluation of the equilibrium partial pressure of SO, as a function of temperature wer an alloy of specified compition. The matched thermochemical diagram consists of a central plot which shows the integral Gibbs' energy of mixing for the binary system SO, at different temperatures. The central plot is flanked on either side by terminal plots of the chemical potentials of oxygen and sulfur, as functions of temperature, for different alloy compositions. By projecting the chemical wtentials of oxygen and sulfur from the terminal lots on to the central diagram, ihe equilibrium partial pressure of S0,can be directly ;cad on the nomograms on the central plot at different temperatures. The matched therrnochemical diagrams are useful in assuring the efficiency of vacuum refining.

Production of niobium powder by electronically mediated reaction (EMR) using calcium as a reductant

Explored in this study is an electronically mediated reaction (EMR) route for the production of niobium powder using calcium as a reductant for niobium oxide (Nb2O5). Feed material, Nb2O5, and reductant calcium alloy containing aluminum and nickel were charged into electronically isolated locations in a molten salt (e.g. CaCl2) at 1173 K. The current flow through an external path between the feed and reductant locations was monitored. A current approximately 0.4 A was measured during the reaction in the external circuit connecting cathode and anode location. Niobium powder with low aluminum and nickel content was obtained although liquid Ca–Al–Ni alloy was used as the reductant. This clearly demonstrates that niobium metal powder can be produced by an electronically mediated reaction (EMR), without direct physical contact between feed (Nb2O5) and reductant (calcium). Mechanism of calciothermic reduction of Nb2O5 in the molten salt is discussed using an isothermal chemical potential diagram.

Analysis of Lined Ducts With Mean Flow, With Application to Dissipative Mufflers

Mufflers with at least one acoustically absorptive duct are generally called dissipative mufflers. Generally, for want of systems approach, these mufflers are characterized by transmission loss of the lined duct with overriding corrections for the terminations, mean flow, etc. In this article, it is proposed that dissipative duct should be integrated with other muffler elements, source impedance and radiation impedance, by means of transfer matrix approach. Towards this end, the transfer matrix for rectangular duct with mean flow has been derived here, for the least attenuated mode. Mean flow introduces a coupling between transverse wave numbers and axial wave number, the evaluation of which therefore calls for simultaneous solution of two or three transcendental equations. This is done by means of a Newton-Raphson iteration scheme, which is illustrated here for square ducts lined with porous ceramic tiles.

Transfer matrix modeling of hyperbolic and parabolic ducts with incompressible mean flow

A general differential equation for the propagation of sound in a variable area duct or nozzle carrying incompressible mean flow (of low Mach number) is derived and solved for hyperbolic and parabolic shapes. Expressions for the state variables of acoustic pressure and acoustic mass velocity of the shapes are derived. Self‐consistent expressions for the four‐pole parameters are developed. The conical, exponential, catenoidal, sine, and cosine ducts are shown to be special cases of hyperbolic ducts. Finally, it is shown that if the mean flow in computing the transmission loss of the mufflers involving hyperbolic and parabolic shapes was not neglected, little practical benefit would be derived.

The generalized Onsager model for the secondary flow in a high-speed rotating cylinder

The generalizations of the Onsager model for the radial boundary layer and the Carrier-Maslen model for the end-cap axial boundary layer in a high-speed rotating cylinder are formulated for studying the secondary gas flow due to wall heating and due to insertion of mass, momentum and energy into the cylinder. The generalizations have wider applicability than the original Onsager and Carrier-Maslen models, because they are not restricted to the limit A >> 1, though they are restricted to the limit R e >> 1 and a high-aspect-ratio cylinder whose length/diameter ratio is large. Here, the stratification parameter A = root m Omega(2)R(2)/2k(B)T). This parameter A is the ratio of the peripheral speed, Omega R, to the most probable molecular speed, root 2k(B)T/m, the Reynolds number Re = rho w Omega R(2)/mu, where m is the molecular mass, Omega and R are the rotational speed and radius of the cylinder, k(B) is the Boltzmann constant, T is the gas temperature, rho(w) is the gas density at wall, and mu is the gas viscosity. In the case of wall forcing, analytical solutions are obtained for the sixth-order generalized Onsager equations for the master potential, and for the fourth-order generalized Carrier-Maslen equation for the velocity potential. For the case of mass/momentum/energy insertion into the flow, the separation-of-variables procedure is used, and the appropriate homogeneous boundary conditions are specified so that the linear operators in the axial and radial directions are self-adjoint. The discrete eigenvalues and eigenfunctions of the linear operators (sixth-order and second-order in the radial and axial directions for the Onsager equation, and fourth-order and second-order in the axial and radial directions for the Carrier-Maslen equation) are determined. These solutions are compared with direct simulation Monte Carlo (DSMC) simulations. The comparison reveals that the boundary conditions in the simulations and analysis have to be matched with care. The commonly used `diffuse reflection' boundary conditions at solid walls in DSMC simulations result in a non-zero slip velocity as well as a `temperature slip' (gas temperature at the wall is different from wall temperature). These have to be incorporated in the analysis in order to make quantitative predictions. In the case of mass/momentum/energy sources within the flow, it is necessary to ensure that the homogeneous boundary conditions are accurately satisfied in the simulations. When these precautions are taken, there is excellent agreement between analysis and simulations, to within 10 %, even when the stratification parameter is as low as 0.707, the Reynolds number is as low as 100 and the aspect ratio (length/diameter) of the cylinder is as low as 2, and the secondary flow velocity is as high as 0.2 times the maximum base flow velocity. The predictions of the generalized models are also significantly better than those of the original Onsager and Carrier-Maslen models, which are restricted to thin boundary layers in the limit of high stratification parameter.

Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 1. DNS and fluctuating force model

The fluctuating force model is developed and applied to the turbulent flow of a gas-particle suspension in a channel in the limit of high Stokes number, where the particle relaxation time is large compared to the fluid correlation time, and low particle Reynolds number where the Stokes drag law can be used to describe the interaction between the particles and fluid. In contrast to the Couette flow, the fluid velocity variances in the different directions in the channel are highly non-homogeneous, and they exhibit significant variation across the channel. First, we analyse the fluctuating particle velocity and acceleration distributions at different locations across the channel. The distributions are found to be non-Gaussian near the centre of the channel, and they exhibit significant skewness and flatness. However, acceleration distributions are closer to Gaussian at locations away from the channel centre, especially in regions where the variances of the fluid velocity fluctuations are at a maximum. The time correlations for the fluid velocity fluctuations and particle acceleration fluctuations are evaluated, and it is found that the time correlation of the particle acceleration fluctuations is close to the time correlations of the fluid velocity in a `moving Eulerian' reference, moving with the mean fluid velocity. The variances of the fluctuating force distributions in the Langevin simulations are determined from the time correlations of the fluid velocity fluctuations and the results are compared with direct numerical simulations. Quantitative agreement between the two simulations are obtained provided the particle viscous relaxation time is at least five times larger than the fluid integral time.

Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 2. Comparison of fluctuating force simulations and experiments

The particle and fluid velocity fluctuations in a turbulent gas-particle suspension are studied experimentally using two-dimensional particle image velocimetry with the objective of comparing the experiments with the predictions of fluctuating force simulations. Since the fluctuating force simulations employ force distributions which do not incorporate the modification of fluid turbulence due to the particles, it is of importance to quantify the turbulence modification in the experiments. For experiments carried out at a low volume fraction of 9.15 x 10(-5) (mass loading is 0.19), where the viscous relaxation time is small compared with the time between collisions, it is found that the gas-phase turbulence is not significantly modified by the presence of particles. Owing to this, quantitative agreement is obtained between the results of experiments and fluctuating force simulations for the mean velocity and the root mean square of the fluctuating velocity, provided that the polydispersity in the particle size is incorporated in the simulations. This is because the polydispersity results in a variation in the terminal velocity of the particles which could induce collisions and generate fluctuations; this mechanism is absent if all of the particles are of equal size. It is found that there is some variation in the particle mean velocity very close to the wall depending on the wall-collision model used in the simulations, and agreement with experiments is obtained only when the tangential wall-particle coefficient of restitution is 0.7. The mean particle velocity is in quantitative agreement for locations more than 10 wall units from the wall of the channel. However, there are systematic differences between the simulations and theory for the particle concentrations, possibly due to inadequate control over the particle feeding at the entrance. The particle velocity distributions are compared both at the centre of the channel and near the wall, and the shape of the distribution function near the wall obtained in experiments is accurately predicted by the simulations. At the centre, there is some discrepancy between simulations and experiment for the distribution of the fluctuating velocity in the flow direction, where the simulations predict a bi-modal distribution whereas only a single maximum is observed in the experiments, although both distributions are skewed towards negative fluctuating velocities. At a much higher particle mass loading of 1.7, where the time between collisions is smaller than the viscous relaxation time, there is a significant increase in the turbulent velocity fluctuations by similar to 1-2 orders of magnitude. Therefore, it becomes necessary to incorporate the modified fluid-phase intensity in the fluctuating force simulation; with this modification, the mean and mean-square fluctuating velocities are within 20-30% of the experimental values.

Flow-acoustic Characterisation of a Cavity-based Combustor Configuration

This study concerns the flow-acoustic characterisation of a cavity-based combustor configuration. A well-validated numerical tool has been used to simulate the unsteady, two-dimensional reacting flow. Initially, a conventional flow over a cavity with dimensions and conditions corresponding to a compact cavity combustor was studied. Cavity mass injections in the form of fuel and air injections required for trapped vortex formation were then employed and the resonance features of this configuration were studied. The results indicate that the cavity depth mode resonance mechanism is dominant at the conditions studied in this work and that the oscillation frequencies do not change with cavity air injection. This observation is important since it implies that the only important variable which can alter resonant frequencies is the cavity depth. With combustion, the pressure oscillation amplitude was observed to increases significantly due to periodic entrainment of the cavity air jet and fluctuation of fuel-air mixture composition to produce highly fluctuating heat-release rates. The underlying mechanisms of the unsteady flow in the cavity combustor identified in this study indicate the strong dependence of the acoustics on the cavity injection strategies.