Performance calculations of new propellant systems

The purpose of this thesis is to investigate the effect on performance and chamber temperature of adding hydrogen to a propellant system. The systems investigated are:

(1) RFNA-Aniline

(2) Nitromethane

(3) Anhydrous hydrazene-liquid oxygen

Since a systematic investigation of the performance parameters of the RFNA-Aniline system over a wide range of mixture ratios has never been made, it was decided to make these calculations, in addition to the investigations stated above.

The results of the calculations can best be summarized by a study of the figures at the end of the thesis. A few generalizations can be made. The effect of adding hydrogen in small quantities to a high temperature system is to increase the performance considerably without too much change in the chamber temperature. As more hydrogen is added, the percentage increase in performance. If hydrogen is added in large quantities, both the performance curve (effective exhaust velocity) and the chamber temperature curve flatten out.

The behavior discussed above is characteristic of hot propellant systems such as RFNA-Aniline and anhydrous hydrazene. In a low temperature system, such as nitromethane, the effect is quite different. The addition of hydrogen in small quantities causes a rapid decrease in chamber temperature, but the increase in performance is considerably less on a percentage basis. As more hydrogen is added the changes in performance and chamber temperature are almost linear.

On the stabilisation of turbulent lifted flames

This study explores the stabilisation mechanisms of turbulent lifted flames by examining the scalar dissipation rate (SDR) of both passive and reactive scalars and their cross dissipation (CDR) in the stabilisation region. DNS results of a laboratory scale hydrogen turbulent lifted flame has been used for this analysis. Various definitions of the flame leading edge (FLE) has been compared and differences are illustrated. Time and spatial averaged statistic of SDR and CDR were examined. It was found that the averaged SDR for the mixture fraction at FLE was well below the reference quenching value for stoichiometric mixture. The averaged SDR for the progress variable is in the same order of the unstrained premixed laminar flame value. It was observed that the averaged CDR changed from negative to weakly positive at FLE. The change in sign was explained by a change in the relative alignment of the gradients of mixture fraction and progress variable. It was thus evident that the CDR was a good marker for stabilisation region and an important quantity in stabilisation mechanism.

Image sequence segmentation using mixture-models with uniqueness and spatio-temporal consistency constraints

The use of mixture-model techniques for motion estimation and image sequence segmentation was discussed. The issues such as modeling of occlusion and uncovering, determining the relative depth of the objects in a scene, and estimating the number of objects in a scene were also investigated. The segmentation algorithm was found to be computationally demanding, but the computational requirements were reduced as the motion parameters and segmentation of the frame were initialized. The method provided a stable description, in whichthe addition and removal of objects from the description corresponded to the entry and exit of objects from the scene.

The internal structure of igniting turbulent sprays as revealed by complex chemistry DNS

A parametric study of spark ignition in a uniform monodisperse turbulent spray is performed with complex chemistry three-dimensional Direct Numerical Simulations in order to improve the understanding of the structure of the ignition kernel. The heat produced by the kernel increases with the amount of fuel evaporated inside the spark volume. Moreover, the heat sink by evaporation is initially higher than the heat release and can have a negative effect on ignition. With the sprays investigated, heat release occurs over a large range of mixture fractions, being high within the nominal flammability limits and finite but low below the lean flammability limit. The burning of very lean regions is attributed to the diffusion of heat and species from regions of high heat release, and from the spark, to lean regions. Two modes of spray ignition are reported. With a relatively dilute spray, nominally flammable material exists only near the droplets. Reaction zones are created locally near the droplets and have a non-premixed character. They spread from droplet to droplet through a very lean interdroplet spacing. With a dense spray, the hot spark region is rich due to substantial evaporation but the cold region remains lean. In between, a large surface of flammable material is generated by evaporation. Ignition occurs there and a large reaction zone propagates from the rich burned region to the cold lean region. This flame is wrinkled due to the stratified mixture fraction field and evaporative cooling. In the dilute spray, the reaction front curvature pdf contains high values associated with single droplet combustion, while in the dense spray, the curvature is lower and closer to the curvature associated with gaseous fuel ignition kernels. © 2011 The Combustion Institute.

Complex chemistry DNS of n-heptane spray autoignition at high pressure and intermediate temperature conditions

Direct Numerical Simulations (DNS) of turbulent n-heptane sprays autoigniting at high pressure (P=24bar) and intermediate air temperature (Tair=1000K) have been performed to investigate the physical mechanisms present under conditions where low-temperature chemistry is expected to be important. The initial turbulence in the carrier gas, the global equivalence ratio in the spray region, and the initial droplet size distribution of the spray were varied. Results show that spray ignition exhibits a spotty nature, with several kernels developing independently in those regions where the mixture fraction is close to its most reactive value ξMR (as determined from homogeneous reactor calculations) and the scalar dissipation rate is low. Turbulence reduces the ignition delay time as it promotes mixing between air and the fuel vapor, eventually resulting in lower values of scalar dissipation. High values of the global equivalence ratio are responsible for a larger number of ignition kernels, due to the higher probability of finding regions where ξ=ξMR. Spray polydispersity results in the occurrence of ignition over a wider range of mixture fraction values. This is a consequence of the inhomogeneities in the mixing field that characterize these sprays, where poorly mixed rich spots are seen to alternate with leaner ones which are well-mixed. The DNS simulations presented in this work have also been used to assess the applicability of the Conditional Moment Closure (CMC) method to the simulation of spray combustion. CMC is found to be a valid method for capturing spray autoignition, although care should be taken in the modelling of the unclosed terms appearing in the CMC equations. © 2013 The Combustion Institute.

Low-order modelling of the response of ducted flames in annular geometries

The adoption of lean premixed prevaporised combustion systems can reduce NOx emissions from gas turbines, but unfortunately also increases their susceptibility to thermoacoustic instabilities. Initially, acoustic waves can produce heat release fluctuations by a variety of mechanisms, often by perturbing the equivalence ratio. If correctly phased, heat release fluctuations can subsequently generate more acoustic waves, which at high amplitude can result in significant structural damage to the combustor. The prediction of this phenomenon is of great industrial interest. In previous work, we have coupled a physics based, kinematic model of the flame with a network model to provide the planar acoustic response necessary to close the feedback loop and predict the onset and amplitude of thermoacoustic instabilities in a lab-scale, axisymmetric single burner combustor. The advantage of a time domain approach is that the modal interaction, the influence of harmonics, and flame saturation can be investigated. This paper extends this approach to more realistic, annular geometries, where both planar and circumferential modes must be considered. In lean premixed prevaporised combustors, fluctuations in equivalence ratio have been shown to be a dominant cause of unsteady combustion. These can occur, for example, due to velocity perturbations in the premix ducts, which can lead to equivalence ratio fluctuations at the fuel injectors, which are subsequently convected downstream to the flame surfaces. Here, they can perturb the heat release by locally altering the flame speed, enthalpy of combustion, and, indirectly, the flame surface area. In many gas turbine designs, particularly aeroengines, the geometries are composed of a ring of premix ducts linking a plenum and an annular combustor. The most unstable modes are often circumferential modes. The network model is used to characterise the flow response of the geometry to heat fluctuations at an appropriate location, such as the fuel injectors. The heat release at each flame holder is determined in the time domain using the kinematic flame model derived, as a function of the flow perturbations in the premix duct. This approach is demonstrated for an annular ring of burners on a in a simple geometry. The approach is then extended to an industrial type gas turbine combustor, and used to predict the limit cycle amplitudes. Copyright © 2012 by ASME.

Modelling of turbulent lifted jet flames using flamelets: A priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z-c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z-c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences. © 2014 Taylor & Francis.

Low-temperature heat capacities and derived thermodynamic functions of cyclohexane

The low-temperature heat capacities of cyclohexane were measured in the temperature range from 78 to 350 K by means of an automatic adiabatic calorimeter equipped with a new sample container adapted to measure heat capacities of liquids. The sample container was described in detail. The performance of this calorimetric apparatus was evaluated by heat capacity measurements on water. The deviations of experimental heat capacities from the corresponding smoothed values lie within +/-0.3%, while the inaccuracy is within +/-0.4%, compared with the reference data in the whole experimental temperature range. Two kinds of phase transitions were found at 186.065 and 279.684 K corresponding solid-solid and solid-liquid phase transitions, respectively. The entropy and enthalpy of the phase transition, as well as the thermodynamic functions {H-(T)- H-298.15 K} and {S-(T)-S-298.15 K}, were derived from the heat capacity data. The mass fraction purity of cyclohexane sample used in the present calorimetric study was determined to be 99.9965% by fraction melting approach.

Kinetic and electrochemical properties of icosahedral quasicrystalline Ti45Zr35Ni17Cu3 powder

Kinetic and electrochemical properties of icosahedral quasicrystalline Ti45Zr35Ni17Cu3 alloy powder as negative electrode material of Ni-MH battery have been investigated at different temperatures. The calculated results show that the apparent activation enthalpy of the charge-transfer reaction is 43.89 kJ mol(-1), and the activation energy of hydrogen diffusion is 21.03 kJ mol(-1). The exchange current density and the diffusion coefficient of hydrogen in the bulky electrode increase with increasing temperature, indicating that increasing temperature is beneficial to charge-transfer reaction and hydrogen diffusion. As a result, the maximum discharge capacity, activation property and high-rate dischargeability are greatly improved with increasing temperature. However, the charge retention and the cycling stability degrade with the increase of the temperature.

Crystallization behavior of carbon nanotubes-filled polyamide 1010

The nanocomposites of polyamide1010 (PA1010) filled with carbon nanotubes (CNTs) were prepared by melt mixing techniques. The isothermal melt-crystallization kinetics and nonisothermal crystallization behavior of CNTs/PA1010 nanocomposites were investigated by differential scanning calorimetry. The peak temperature, melting point, half-time of crystallization, enthalpy of crystallization, etc. were measured. Two stages of crystallization are observed, including primary crystallization and secondary crystallization. The isothermal crystallization was also described according to Avrami's approach. It has been shown that the addition of CNTs causes a remarkable increase in the overall crystallization rate of PA1010 and affects the mechanism of nucleation and growth of PA1010 crystals. The analysis of kinetic data according to nucleation theories shows that the increment in crystallization rate of CNTs/PA1010 composites results from the decrease in lateral surface free energy.

Pressure effects on the thermodynamics of trans-decahydronaphthalene/polystyrene polymer solutions: Application of the Sanchez-Lacombe lattice fluid theory

The cloud-point temperatures (T-cl's) of trans-decahydronaphthalene(TD)/polystyrene (PS, (M) over bar (w) = 270 000) solutions were determined by light scattering measurements over a range of temperatures (1-16degreesC), pressures (100-900 bar), and compositions (4.2-21.6 vol.-% polymer). The system phase separates upon cooling and T-cl was found to increase with rising pressure for constant composition. In the absence of special effects, this finding indicates positive excess volume for the mixing. Special attention was paid to the demixing temperatures as a function of pressure for different polymer solutions and the plots in the T-phi plane (where phi signifies volume fractions). The cloud-point curves of polymer solutions under different pressures were observed for different compositions, which demonstrated that pressure has a greater effect on the TD/PS solutions when far from the critical point as opposed to near the critical point. The Sanchez-Lacombe lattice fluid theory (SLLFT) was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalpy of mixing, and the volume changes of mixing. The calculated results show that modified PS scaling parameters can describe the thermodynamics of the TD/PS system well. Moreover the SLLFT describes the experimental results well.

Study on a carbosilane liquid crystalline (LC) dendrimer of the first generation - Containing twelve 4-butoxyazobenzene mesogenic groups in periphery

The divergent synthesis of a new carbosilane liquid-crystalline (LC) dendrimer of the first generation (D1) is described. Twelve 4-butoxyazobenzene groups are used as mesogenic fragments and attached in the periphery of the molecule. Structure and properties of D1 were characterized by element analysis, H-1 NMR, MALDI-TOF-MS, IR, UV-Vis, polarizing optical micrograph, DSC and WAXD. It is argued that mesophase of nematic type is realized. It is shown that the mesophase type of the dendrimer essentially depends on the chemical nature of the mesogenic groups. Phase behavior of D1 is K82N1331132N67K. The melting point of D1 is 30similar to43 degreesC lower than that of M5, its clearing temperature is 9 similar to 11 degreesC higher than that of M5 and its mesophase region is enlarged by 39 similar to 54 degreesC compared to that of M5. Eight extinguished brushes emanating from a stationary point are observed, corresponding to the high-strength disclination of S = + 2 of dendrimer. The clearing enthalpy of D1 is smaller than the value that is commonly found for phase transition n-i in LC and LC polymers. This may be due to the presence of branched dendrimer cores which cannot be easily deformed to fit into the anisotropic LC phase structure.

Theoretical estimation of thermodynamic properties of the system PS/PPO on the basis of modified combining rule of Sanchez-Lacombe lattice fluid model

The three scaling parameters described in Sanchez-Lacombe lattice fluid theory (SLLFT), T*, P* and rho* of pure polystyrene (PS), pure poly(2,6-dimethyl-1,4-phenylene oxide) (PPO) and their mixtures are obtained by fitting corresponding experimental pressure volume-temperature data with equation-of-state of SLLFT. A modified combining rule in SLLFT used to match the volume per mer, v* of the PS/PPO mixtures was advanced and the enthalpy of mixing and Flory-Huggins (FH) interaction parameter were calculated using the new rule. It is found that the difference between the new rule and the old one presented by Sanchez and Lacombe is quite small in the calculation of the enthalpy of mixing and FH interaction parameter and the effect of volume-combining rule on the calculation of thermodynamic properties is much smaller than that of energy-combining rule. But the relative value of interaction parameter changes much due to the new volume-based combining rule. This effect can affect the position of phase diagram very much, which is reported elsewhere [Macromolecules 34 (2001) 6291]

Effect of PMR-POI on the melt behavior and isothermal crystallization of poly(phenylene sulfide)

The crystallization behavior of neat PPS and PPS in blends with PMR-POI prepared by melt mixing were investigated by differential scanning calorimetry (DSC). It was found that POI was an effective nucleation agent of the crystallization for PPS. The enthalpy of crystallization of PPS in the blends increased compared with that of neat PPS. During isothermal crystallization from melt, the dependence of relative degree of crystallinity on time was described by the Avrami equation. It has been shown that the addition of POI causes an increase in the overall crystallization rate of PPS; it also changed the mechanism of nucleation of the PHB crystals from homogeneous nucleation to heterogeneous nucleation. The equilibrium melting temperature of PPS and PPS/POI blends were determined. The analysis of kinetic data according to nucleation theories shows that the increase in crystallization rate of PPS in the composite is due to the decrease in surface energy of the extremity surface.

Thermodynamics of PMMA/SAN blends: Application of the Sanchez-Lacombe lattice fluid theory

Polymer blends of poly(methyl methacrylate) (PMMA) and poly(styrene-co-acrylonitrile) (SAN) with an acrylonitrile content of about 30 wt % were prepared by means of solution-casting and characterized by virtue of pressure-volume-temperature (PVT) dilatometry. The Sanchez-Lacombe (SL) lattice fluid theory was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalpy of the mixing, the volume change of the mixing, and the combinatorial and vacancy entropies of the mixing for the PMMA/SAN system. A new volume-combining rule was used to evaluate the close-packed volume per mer, upsilon*, of the PMMA/SAN blends. The calculated results showed that the new and the original volume-combining rules had a slight influence on the FH interaction parameter, the enthalpy of the mixing, and the combinatorial entropy of the mixing. Moreover, the spinodals and the binodals calculated with the SL theory by means of the new volume-combining rule could coincide with the measured data for the PMMA/SAN system with a lower critical solution temperature, whereas those obtained by means of the original one could not.