固体火箭发动机喷管气固两相流动的数值模拟

对颗粒相采用颗粒轨道模型,气相求解可压缩N-S方程组,计算方法采用显式Runge-Kutta时间推进法与有总变差衰减(TVD)性质的高精度MUSCL-Roe格式;自主开发了曲线坐标系下二维轴对称可压缩N-S方程组的解算器Solve2D,研究了固体火箭发动机喷管中颗粒相对流场的影响以及不同尺寸颗粒运动规律.结果表明:颗粒相对流场的影响主要表现在喷管喉部以及扩张段,和单相流场相比,沿轴线马赫数减小,且颗粒尺寸越小减少得越多;沿轴线气相温度升高,且颗粒尺寸越小温度升高越多;颗粒尺寸越小,无粒子区越小;颗粒越大与收缩段壁面碰撞越剧烈,无粒子区越大.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

An investigation of a method of underwater propulsion by direct gas injection

This report presents the results of an investigation of a method of underwater propulsion. The propelling system utilizes the energy of a small mass of expanding gas to accelerate the flow of a large mass of water through an open ended duct of proper shape and dimensions to obtain a resultant thrust. The investigation was limited to making a large number of runs on a hydroduct of arbitrary design, varying between wide limits the water flow and gas flow through the device, and measuring the net thrust caused by the introduction and expansion of the gas.

In comparison with the effective exhaust velocity of about 6,000 feet per second observed in rocket motors, this hydroduct model attained a maximum effective exhaust velocity of more than 27,000 feet per second, using nitrogen gas. Using hydrogen gas, effective exhaust velocities of 146,000 feet per second were obtained. Further investigation should prove this method of propulsion not only to be practical but very efficient.

This investigation was conducted at Project No. 1, Guggenheim Aeronautical Laboratory, California Institute of Technology, Pasadena, California.

Structure of cryogenic flames at elevated pressures

This paper presents new experimental results on cryogenic jet flames formed by a coaxial injector at a pressure of 70 bar, which approaches the pressures found in rocket engines. This element, fed with liquid oxygen and gaseous hydrogen, is placed in a square combustion chamber equipped with quartz windows. The flame is examined via spectroscopy, OH* emission, and backlighting, the aim being to provide basic information on the flame structure. It is found that some of the OH* emission is absorbed by the OH radicals present in the flame. A detailed examination of this effect is presented, in which it is shown that, for this turbulent flame, the Abel transform gives the position of the intense reaction region, whether or not absorption is signficant. The flame is attached to the oxygen injector, as at low pressure. At high pressure, flame expansion is reduced compared with low pressure and is also less dependent on the momentum flux ratio between the hydrogen and the oxygen streams. An analysis of the relevant Damköhler numbers suggests that this is because the rate of combustion is mainly controlled by large-scale turbulent mixing at high pressure, and it is dominated by jet break-up, atomization, and vaporization at low pressures. Jet break-up is particularly dependent on the momentum flux ratio. Finally, the mean volumetric heat release rates and flame surface density in the experimental facility are estimated.

激光衰减法测量固体火箭发动机尾流粒子参数

为了测量某固体火箭发动机尾流粒子参数,设计了基于激光衰减法的实验系统,并在发动机试车中采取了消除尾流辐射干扰的措施.将实验结果与参考文献的激光衰减法测量结果进行了比较,结果一致,表明所获得光强衰减信号和测量结果是可信的.在发动机燃烧室压强为9.5 MPa时,利用激光衰减方法测量喷管出口的粒子尺度d32=0.9μm,实验数据误差为20%

Transition Space

Informal causal descriptions of physical systems abound in sources such as encyclopedias, reports and user's manuals. Yet these descriptions remain largely opaque to computer processing. This paper proposes a representational framework in which such descriptions are viewed as providing partial specifications of paths in a space of possible transitions, or transition space. In this framework, the task of comprehending informal causal descriptions emerges as one of completing the specifications of paths in transition space---filling causal gaps and relating accounts of activity varied by analogy and abstraction. The use of the representation and its operations is illustrated in the context of a simple description concerning rocket propulsion.

Mapping the structure of perceptual and visual-motor abilities in healthy young adults.

The ability to quickly detect and respond to visual stimuli in the environment is critical to many human activities. While such perceptual and visual-motor skills are important in a myriad of contexts, considerable variability exists between individuals in these abilities. To better understand the sources of this variability, we assessed perceptual and visual-motor skills in a large sample of 230 healthy individuals via the Nike SPARQ Sensory Station, and compared variability in their behavioral performance to demographic, state, sleep and consumption characteristics. Dimension reduction and regression analyses indicated three underlying factors: Visual-Motor Control, Visual Sensitivity, and Eye Quickness, which accounted for roughly half of the overall population variance in performance on this battery. Inter-individual variability in Visual-Motor Control was correlated with gender and circadian patters such that performance on this factor was better for males and for those who had been awake for a longer period of time before assessment. The current findings indicate that abilities involving coordinated hand movements in response to stimuli are subject to greater individual variability, while visual sensitivity and occulomotor control are largely stable across individuals.

El hombre metrosexual como Narciso contemporáneo: el papel del discurso publicitario actual en la construcción de la masculinidad

Existen diferentes definiciones del concepto de masculinidad. Por un lado están las normativas, que la entienden como lo que los hombres deberían ser y, por otro, los enfoques semióticos, que la definen como no-feminidad. Partiendo de este último enfoque y teniendo en cuenta que la feminidad ha evolucionado en los últimos años junto con la sociedad desdibujando la frontera entre lo que forma parte del rol femenino y lo que no, se puede afirmar que han surgido a partir de los renovados modelos de feminidad, diferentes y nuevas masculinidades entre las que se encuentra la del hombre metrosexual. Es en este tipo de masculinidad en concreto en la que se centra la presente investigación y en cómo se construye la narración de dicha masculinidad en el discurso publicitario actual recurriendo a la figura mítica de Narciso. Así, a partir de una revisión bibliográfica y de un análisis de contenido, se pretende analizar la influencia de la publicidad en la aparición del hombre metrosexual como nueva masculinidad en respuesta a la transformación de la feminidad, a la par que confirmar la presencia de la figura de Narciso como recurso en la construcción de la masculinidad en la publicidad actual.

Fe XIII emission lines in active region spectra obtained with the Solar Extreme-Ultraviolet Research Telescope and Spectrograph

Recent fully relativistic calculations of radiative rates and electron impact excitation cross-sections for FeXIII are used to generate emission-line ratios involving 3s23p2-3s3p3 and 3s23p2-3s23p3d transitions in the 170-225 and 235-450 Å wavelength ranges covered by the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS). A comparison of these line ratios with SERTS active region observations from rocket flights in 1989 and 1995 reveals generally very good agreement between theory and experiment. Several new FeXIII emission features are identified, at wavelengths of 203.79, 259.94, 288.56 and 290.81 Å. However, major discrepancies between theory and observation remain for several FeXIII transitions, as previously found by Landi and others, which cannot be explained by blending. Errors in the adopted atomic data appear to be the most likely explanation, in particular for transitions which have 3s23p3d1D2 as their upper level. The most useful FeXIII electron-density diagnostics in the SERTS spectral regions are assessed, in terms of the line pairs involved being (i) apparently free of atomic physics problems and blends, (ii) close in wavelength to reduce the effects of possible errors in the instrumental intensity calibration, and (iii) very sensitive to changes in Ne over the range 108-1011cm-3. It is concluded that the ratios which best satisfy these conditions are 200.03/202.04 and 203.17/202.04 for the 170-225 Å wavelength region, and 348.18/320.80, 348.18/368.16, 359.64/348.18 and 359.83/368.16 for 235-450 Å.

SiIX emission lines in spectra obtained with the Solar EUV Research Telescope and Spectrograph (SERTS)

Theoretical electron-density-sensitive emission line ratios involving 2s(2)2p(2)-2s2p(3) transitions in Si IX between 223 and 350 Angstrom are presented. A comparison of these with an extensive dataset of solar-active-region, quiet-Sun, subflare and off-limb observations, obtained during rocket flights by the Solar EUV Research Telescope and Spectrograph (SERTS), reveals generally very good agreement between theory and experiment. This provides support for the accuracy of the line- ratio diagnostics, and hence the atomic data on which they are based. In particular, the density-sensitive intensity ratio I (258.10 Angstrom)/ I (349.87 Angstrom) offers an especially promising diagnostic for studies of coronal plasmas, as it involves two reasonably strong emission lines and varies by more than an order of magnitude over the useful density range of 10(9)-10(11) cm(-3). The 2s(2)2p(2) S-1(0) - 2s2p(3) P-1(1) transition at 259.77 Angstrom is very marginally identified for the first time in the SERTS database, although it has previously been detected in solar flare observations.

Emission lines of Na-like ions in spectra obtained with the Solar EUV Research Telescope and Spectrograph (SERTS)

Theoretical emission-line ratios involving transitions in the 236-412 Angstrom wavelength range are presented for the Na-like ions Ar viii, Cr xiv, Mn xv, Fe xvi, Co xvii, Ni xviii and Zn xx. A comparison of these with an extensive data set of the solar active region, quiet-Sun, subflare and off-limb observations, obtained during rocket flights by the Solar EUV Research Telescope and Spectrograph (SERTS), reveals generally very good agreement between theory and experiment. This indicates that most of the Na-like ion lines are reliably detected in the SERTS observations, and hence may be employed with confidence in solar spectral analyses. However, the features in the SERTS spectra at 236.34 and 300.25 Angstrom, originally identified as the Ni xviii 3p (2) P-3/2 -3d (2) D- 3/2 and Cr xiv 3p (2) P-3/2 -3d (2) D-5/2 transitions, respectively, are found to be due to emission lines of Ar xiii (236.34 Angstrom) and possibly S v or Ni vi (300.25 Angstrom). The Co xvii 3s (2) S-3p (2) P-3/2 line at 312.55 Angstrom is always badly blended with an Fe xv feature at the same wavelength, but Mn xv 3s (2) S-3p (2) P-1/2 at 384.75 Angstrom may not always be as affected by second-order emission from Fe xii 192.37 Angstrom as previously thought. On the other hand, we find that the Zn xx 3s (2) S-3p (2) P-3/2 transition can sometimes make a significant contribution to the Zn xx/Fe xiii 256.43- Angstrom blend, and hence care must be taken when using this feature as an Fe xiii electron density diagnostic. A line in the SERTS-89 active region spectrum at 265.00 Angstrom has been re-assessed, and we confirm its identification as the Fe xvi 3p (2) P-3/2 -3d (2) D-3/2 transition.

A comparison of theoretical Si VIII emission line ratios with observations from SERTS

Recent R-matrix calculations of electron impact excitation rates in N-like Si VIII are used to derive theoretical emission line intensity ratios involving 2s(2)2p(3)-2s2p(4) transitions in the 216 -320 Angstrom wavelength range. A comparison of these with an extensive dataset of solar active region, quiet- Sun, sub-flare and off-limb observations, obtained during rocket flights of the Solar EUV Research Telescope and Spectrograph (SERTS), indicates that the ratio R-1 = I(216.94 Angstrom)/I(319.84 Angstrom) may provide a usable electron density diagnostic for coronal plasmas. The ratio involves two lines of comparable intensity, and varies by a factor of about 5 over the useful density range of 10(8)-10(11) cm(-3). However R-2 = I(276.85 Angstrom)/I(319.84 Angstrom) and R-3 = I(277.05 Angstrom)/I(319.84 Angstrom) show very poor agreement between theory and observation, due to the severe blending of the 276.85 and 277.05 Angstrom lines with Si VII and Mg VII transitions, respectively, making the ratios unsuitable as density diagnostics. The 314.35 Angstrom feature of Si VIII also appears to be blended, with the other species contributing around 20% to the total line flux.

Extreme-ultraviolet emission lines of S X in solar flare and active region spectra

R-matrix calculations of electron impact excitation rates in N- like S x are used to derive theoretical emission-line intensity ratios involving 2s(2)2p(3)-2s2p(4) transitions in the 189-265 Angstrom wavelength range. A comparison of these with observational data for solar flares and active regions, obtained with the Naval Research Laboratory's S082A spectrograph on board Skylab and the Solar EUV Rocket Telescope and Spectrograph, reveals that many of the S x lines in the spectra are badly blended with emission features from other species. However, the intensity ratios I(228.70 Angstrom)/I(264.24 Angstrom) and I(228.70 Angstrom)/I(259.49 Angstrom) are found to provide useful electron density diagnostics for flares, although the latter cannot be employed for active regions, because of blending of the 259.49 Angstrom line with an unidentified transition in these solar features.

A comparison of theoretical MgVI emission line strengths with active-region observations from serts

R-matrix calculations of electron impact excitation rates in N-like Mg vi are used to derive theoretical electron-density-sensitive emission line ratios involving 2s(2)2p(3)-2s2p(4) transitions in the 269-403 Angstrom wavelength range. A comparison of these with observations of a solar active region, obtained during the 1989 flight of the Solar EUV Rocket Telescope and Spectrograph (SERTS), reveals good agreement between theory and observation for the 2s(2)2p(3) S-4-2s2p(4) P-4 transitions at 399.28, 400.67, and 403.30 Angstrom, and the 2s(2)2p(3) P-2-2s2p(4) D-2 lines at 387.77 and 387.97 Angstrom. However, intensities for the other lines attributed to Mg vi in this spectrum by various authors do not match the present theoretical predictions. We argue that these discrepancies are not due to errors in the adopted atomic data, as previously suggested, but rather to observational uncertainties or mis-identifications. Some of the features previously identified as Mg vi lines in the SERTS spectrum, such as 291.36 and 293.15 Angstrom, are judged to be noise, while others (including 349.16 Angstrom) appear to be blended.

Electron densities in the coronae of the Sun and Procyon from extreme-ultraviolet emission line ratios in Fe XI

New R-matrix calculations of electron impact excitation rates for Fe XI are used to determine theoretical emission line ratios applicable to solar and stellar coronal observations. These are subsequently compared to solar spectra of the quiet Sun and an active region made by the Solar EUV Rocket Telescope and Spectrograph (SERTS-95), as well as Skylab observations of two flares. Line blending is identified, and electron densities of 10(9.3), 10(9.7), greater than or equal to 10(10.8), and greater than or equal to 10(11.3) cm(-3) are found for the quiet Sun, active region, and the two flares, respectively. Observations of the F5 IV-V star Procyon, made with the Extreme Ultraviolet Explorer (EUVE) satellite, are compared and contrasted with the solar observations. It is confirmed that Procyon's average coronal conditions are very similar to those seen in the quiet Sun, with N-e = 10(9.4) cm(-3). In addition, although the quiet Sun is the closest solar analog to Procyon, we conclude that Procyon's coronal temperatures are slightly hotter than solar. A filling factor of 25(-12)(+38)% was derived for the corona of Procyon.