Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.

Relativistic Continuum Random Phase Approximation and Applications I. Formalism

A fully consistent relativistic continuum random phase approximation (RCRPA) is constructed in terms of the Green's function technique. In this method the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single particle Green's function, which includes also the negative states in the Dirac sea in the nose aapproximation. The theoretical formalism of RCRPA and numerical details are presented. The single particle Green's function is calculated numerically by a proper product of regular and irregular solutions of the Dirac equation. The numerical details and the formalism of RCRPA in the momentum representation are presented.

Investigation of Chemical Bond Characteristics, Thermal Expansion Coefficients and Bulk Moduli of alpha-R2MoO6 and R2Mo2O7 (R = rare earths) by Using a Dielectric Chemical Bond Method

Theoretical researches are performed on the alpha-R2MoO6 (R = Y, Gd, Tb Dy, Ho, Er, Tm and Yb) and pyrochlore-type R2Mo2O7 (R = Y, Nd, Sm, Gd, Tb and Dy) rare earth molybdates by using chemical bond theory of dielectric description. The chemical bonding characteristics and their relationship with thermal expansion property and compressibility are explored. The calculated values of linear thermal expansion coefficient (LTEC) and bulk modulus agree well with the available experimental values. The calculations reveal that the LTECs and the bulk moduli do have linear relationship with the ionic radii of the lanthanides: the LTEC decreases from 6.80 to 6.62 10(-6)/K and the bulk modulus increases from 141 to 154 GPa when R goes in the order Gd, Tb Dy, Ho, Er, Tm, and Yb in the alpha-R2MoO6 series; while in the R2Mo2O7 series, the LTEC ranges from 6.80 to 6.61 10(-6)/K and the bulk modulus ranges from 147 to 163 GPa when R varies in the order Nd, Sm, Gd, Tb and Dy.

First Principles Investigation on the Ultra-Incompressible and Hard TaN

The structural, electronic, and mechanical properties of TaN were investigated by use of the density functional theory (DFT). Eight structures were considered, i.e.. hexagonal WC TaN, NiAs, wurtzite, and CoSn structures. cubic NaCl. zinc-blende and CsCl structures. The results indicate that TaN in TaN-type structure is the most stable at ambient conditions among the considered structures. Above 5 GPa, TaN in WC-type structure becomes energetically the most stable phase. They are also stable both thermodynamically and mechanically. TaN in WC-type has the largest shear Modulus 243 GPa and large bulk modulus 337 GPa among the considered structures. The Volume compressibility is slightly larger than diamond, but smaller than c-BN at pressures from 0 to 100 GPa. The compressibility along the c axis is smaller than the linear compressibility of both diamond and c-BN.

Calculation of the bulk modulus of simple and complex crystals with the chemical bond method

The relation between the lattice energies and the bulk moduli on binary inorganic crystals was studied, and the concept of lattice energy density is introduced. We find that the lattice energy densities are in good linear relation with the bulk moduli in the same type of crystals, the slopes of fitting lines for various types of crystals are related to the valence and coordination number of cations of crystals, and the empirical expression of calculated slope is obtained. From crystal structure, the calculated results are in very good agreement with the experimental values. At the same time, by means of the dielectric theory of the chemical bond and the calculating method of the lattice energy of complex crystals, the estimative method of the bulk modulus of complex crystals was established reasonably, and the calculated results are in very good agreement with the experimental values.

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.

Abnormal pressure dependence of the phase boundaries in TL/PEO/P(EO-b-DMS) ternary mixtures

The cloud-point temperatures (T-cl's) of both binary poly(ethylene oxide) (PEO)-poly(ethylene oxide-b-dimethylsiloxane) [P(EO-b-DMS)] and ternary[toluene/PEO/P(EO-b-DMS)] systems were determined by light scattering measurements at atmospheric pressure. The phase separation behavior upon cooling in the ternary system has been investigated at atmospheric pressure and under high pressure and compared to the phase behavior in the binary system. The phase transition temperatures have been obtained for all of the samples. As a result, the pressure induces compatibility in the binary mixtures, but for the ternary system, pressure not only can induce mixing but also can induce phase separation.

Study of PVT behavior of biaxially oriented polypropylene resins in the melting state

The PVT data of five kinds of biaxially oriented polypropylene (BOPP) Resins was measured by the PVT-100 apparatus. Thermal expansion coefficients (alpha) and isothermal compressibility (beta) were evaluated from Tait equation in the melting state and then compared with those fitted with the value of experiment. The results showed that it was reasonable to calculate alpha and beta with Tait equation in the melting state. At the same time, it was found that thermal expansion coefficients, isothermal compressibility and the melting temperature (T-m) of one of BOPP melts (S28C) were lower than those of the others in the same test conditions, indicating that the volume deformation of S28C resin is' less so that it could be realized to avoid arising surface defects of the film (biaxially oriented polypropylene film) due to. contracting, thereby decrease damage to the film in the subsequence process. Accordingly superior processing properties of S28C resin are confirmed from PVT. speciality.

Kinetic isotope effects of proton transfer reaction for amines by ultrasonic relaxation methods: Unexpected evidence from the results in ethylamine solutions

Ultrasonic absorption coefficients for ethylamine in heavy water (D2O) and in light water (H2O) have been measured in the frequency range from 0.8 to 220 MHz at 25 degrees C. A single relaxational process has been observed in these two kinds of solutions. From the concentration dependence of the ultrasonic relaxation parameters, and following the reaction mechanism proposed by Eigen et al. for ethylamine in H2O, the causes of the relaxations have been attributed to a perturbation of an equilibrium associated with a deuteron or proton transfer reaction. The rate and equilibrium constants have been estimated from deuterioxide or hydroxide ion concentration dependence of the relaxation frequency, and the kinetic isotope effects have been determined. In addition, the standard volume changes of the reactions have been calculated from the concentration dependence of the maximum absorption per wavelength, and the adiabatic compressibility has also been determined from the density and sound velocity for ethylamine in D2O and in H2O, respectively. These results are compared with those for propylamine and butylamine and are discussed in relation to the different kinetic properties between D2O and H2O, the reaction radii derived by Debye theory, and the structural properties of the reaction intermediate.

Ultrasonic relaxation kinetics on fast deuteron transfer reaction

Propylamine has been selected to investigate the isotope effect of a fast deuteron transfer reaction by ultrasonic relaxation method. Ultrasonic absorption coefficients of propylamine in heavy water (D2O) at 25 degrees C in the concentration range from 0.0107 to 0.6300 mol dm(-3) have been measured by pulse and resonance methods over the frequency range from 0.8 to 220 MHz. A Debye-type single relaxation absorption has been observed in the solution. From the dependence of the ultrasonic relaxation parameters on the concentration and solution pH, the source of the observed relaxation has been attributed to a perturbation of the chemical equilibrium associated with the deuteron transfer reaction. The rate and equilibrium constants have been determined by the measurement of the deuteroxyl ion concentration dependence of the relaxation frequency. Also the standard volume change of the reaction has been determined from the concentration dependence of the maximum absorption per wavelength and the adiabatic compressibility has been calculated from the density and the sound velocity in the solution. These results have then been compared with those obtained for propylamine in light water (H2O). The forward rate constant is greater and the reverse rate constant is smaller in DO than in H2O. The standard volume change for deuteron transfer is greater than that for proton transfer reaction, and the adiabatic compressibility shows a similar trend. These data support an argument that there exists a stronger hydrogen bond in D2O than in H2O. The difference of the stability in the intermediate states, R-ND3+... OD- and R-NH3+... OH-, has also been considered from the results of the isotope effects.

STATISTICAL THERMODYNAMICS OF POLYDISPERSE POLYMERS

A statistical thermodynamics theory of polydisperse polymers based on a lattice model of fluids is formulated. Pure polydisperse polymer can be completely characterized by three scale factors and the molecular weight distribution of the system. The equation of state does not satisfy a simple corresponding-states principle, except for a polymer fluid of sufficiently high molecular weight. The relationships between thermal expansion coefficient alpha and isothermal compressibility beta with reduced variables are also predicted.

深水滑翔机器人耐压壳体结构优化设计

针对深水滑翔机器人耐压壳体在深水中压缩变形和海水密度随水深增大而变化等因素造成驱动浮力变化问题,利用压缩变形和密度变化对驱动浮力影响的互补特性,归纳了载体耐压壳体结构优化设计方法。保证耐压壳体抗压条件和质量最小前提下采用结构优化设计方法,降低耐压壳体压缩变形和海水密度变化对驱动浮力的影响,提高系统能源的利用率。

天津滨海新区软土不同含砂量的强度特征研究

Based on the features of soft soil in Tianjing Coastal New Developing Area, this kind of soil with different content of sand was researched systematically, according to the indoor experiment, about its characteristics of strength and deformation. The main results are summarized in the following: Firstly, on the basis of geological engineering investigation, the systemic experiments about the physical characteristics were conducted. The test soil samples were taken from the gray and gray-yellow silty soft soil which was formatted by near-shore marine sediment and marine-continental interactive sediment. The original condition of the sample soil was in saturation and the basic indexes are: liquid limit36.1%, plastic limit 18.8%, plasticity index. Then, the condensation characteristics of the soft soil were analyzed through high-pressure consolidation tests. The results show that，in various loading serials, the coefficient of compressibility under P=100kPa and 200kPa are all larger than 0.5MPa-1. So the sample soil is a kind of high-compressibility soil. Secondly, triaxial strength of undisturbed soil and remoulded soil was researched by using triaxial test. The types of stress-strain curve of both undisturbed and remoulded soil are the stress stabilization and softening type, which show the specific plastic character. Furthermore, the cohesion and friction angle of undisturbed soil changes, when the ambient pressure is different, instead of a stable value for all time; the cohesion and friction angle of remoulded soil changes with the compactness and sand-carrying capacity which is wholly higher than undisturbed soil. At last, the stress-strain results of both undisturbed and remoulded soil were normalized by using the ambient pressure as the normalization factor. The results show that, there are all some normalization characters in both undisturbed and remoulded soil, however, the feature of normalization of undisturbed soil is worse than the remoulded ones. The main reason is that the undisturbed samples are worse in equality and the unavoidable disturb through the process of sampling and experiments will also make them can not put up good normalization. Therefore, it is feasible to normalize the soil in Tianjing Coastal New Developing Area with the ambient pressure as normalization factor.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.