潜入式摆动喷管两相内流场数值模拟

气相采用三维薄层近似N S方程,固相采用拉格朗日坐标系下颗粒轨道模型,计算方法采用显式ENO差分格式和分区多块网格技术,建立了固体火箭发动机内流场模拟软件。对轴对称JPL喷管进行模拟,计算结果与文献进行比较,符合较好。推广到实际三维装药型面时,发现颗粒数目增加导致计算效率降低。采用软件中的气相计算模块对潜入式摆动喷管内流场进行数值模拟,比较了喷管摆动角在发动机喉道附近引起流动变化,研究了非对称流动对喷管内表面压强的影响。

Effect of Boundary Layers on Polycrystalline Silicon Chemical Vapor Deposition in a Trichlorosilane and Hydrogen System

This paper presents the numerical investigation of the effects of momentum, thermal and species boundary layers on the characteristics of polycrystalline silicon deposition by comparing the deposition rates in three chemical vapor deposition （CVD） reactors. A two-dimensional model for the gas flow, heat transfer, and mass transfer was coupled to the gas-phase reaction and surface reaction mechanism for the deposition of polycrystalline silicon from trichlorosilane （TCS）-hydrogen system. The model was verified by comparing the simulated growth rate with the experimental and numerical data in the open literature. Computed results in the reactors indicate that the deposition characteristics are closely related to the momentum, thermal and mass boundary layer thickness. To yield higher deposition rate, there should be higher concentration of TCS gas on the substrate, and there should also be thinner boundary layer of HCl gas so that HCl gas could be pushed away from the surface of the substrate immediately.

Gustiness and coherent structure under weak wind period in atmospheric boundary layer

Statistical analysis of turbulent and gusty characteristics in the atmospheric boundary layer under weak wind period has been carried out.The data used in the analysis were from the multilevel ultrasonic anemometer-thermometers at 47 m,120 m,and 280 m levels on Beijing 325 m meteorological tower.The time series of 3D atmospheric velocity were analyzed by using conventional Fourier spectral analysis and decompose into three parts:basic mean flow(period > 10 min),gusty disturbances(1 min < period < 10 min)and turbulence fluctuations(period < 1 min).The results show that under weak mean wind condition:1)the gusty disturbances are the most strong fluctuations,contribute about 60% kinetic energy of eddy kinetic energy and 80% downward flux of momentum,although both the eddy kinetic energy and momentum transport are small in comparison with those in strong mean wind condition;2)the gusty wind disturbances are anisotropic;3)the gusty wind disturbances have obviously coherent structure,and their horizontal and vertical component are negatively correlated and make downward transport of momentum more effectively;4)the friction velocities related to turbulence and gusty wind are approximately constant with height in the surface layer.

Multi-level CO_(2) fluxes over Beijing megacity with the eddy covariance method

Based on five years of eddy covariance measurements at multiple levels(47,140,and 280 m)of Beijing's 325-m meteorological tower,the exchange process of CO_(2) fluxes between the atmosphere and urban surface were investigated.As a result of the total vehicle control policy from 2011 in Beijing,the growth rate of annual total CO_(2) flux at 140 m is 7.8% from 2008-2010 but 2.3%from 2010-2012.With the minimum vegetation cover and largest population density,the 5-yr average annual total CO_(2) flux at 140 m is largest(6.41 kg C m^(−2) yr^(−1)),compared with that at 47 m(5.78 kg C m^(−2) yr^(−1))and 280 m(3.99 kg C m^(−2) yr^(−1)).With regards to annual total CO_(2) fluxes in Beijing,vehicle numbers and population are the main controlling factors.The measured CO_(2) fluxes were highly dependent on land cover/use in the prevailing wind direction.The CO_(2) fluxes at three layers all correlated positively with road fraction,with the R2 values being 0.69,0.57,and 0.54(P<0.05),respectively.The decreasing fraction of vegetation caused an increasing of the annual total CO_(2) flux,and there was an exponential relationship between them.The annual total CO_(2) fluxes were larger with higher population density.

Boundary States and Correlation Functions of Tricritical Ising Model from Coulomb-Gas Formalism

We consider the minimal conformaJ model describing the tricritical Ising model on the disk and on the upper half plane. Using the coulomb-gas formalism we determine its consistents boundary states as well as its one-point and two-point correlation functions.

Shock wave-boundary layer interactions control by plasma aerodynamic actuation

This study demonstrates the potential for shock wave-boundary layer interaction control in air by plasma aerodynamic actuation.Experimental investigations on shock wave-boundary layer interactions control by plasma aerodynamic actuation are conducted in a Mach 3 in-draft air tunnel.Schlieren imaging shows that the discharges cause the oblique shock to move forward.Schlieren imaging and static pressure probes also show that separation phenomenon shifts backward and the size of separation is enlarged when plasma aerodynamic actuation is applied.The intensity of shock wave is weakened through wall pressure probe.Furthermore,numerical investigations on shock wave-boundary layer interactions control are conducted with plasma aerodynamic actuation.The discharge is modeled as a steady volumetric heat source which is integrated into the energy equation.The input energy level is about 7 kW through discharge process.Results show that the separation phenomenon shifts backward and the intensity of shock is reduced with plasma actuation.These numerical results are consistent with the experimental results.

Estimation of key surface parameters in semi-arid region and their impacts on improvement of surface fluxes simulation

Uncertainties in some key parameters in land surface models severely restrict the improvement of model capacity for successful simulation of surface-atmosphere interaction. These key parameters are related to soil moisture and heat transfer and phy- sical processes in the vegetation canopy as well as other important aerodynamic processes. In the present study, measurements of surface-atmosphere interaction at two observation stations that are located in the typical semi-arid region of China, Tongyu Station in Jilin Province and Yuzhong Station in Gansu Province, are combined with the planetary boundary layer theory to estimate the value of two key aerodynamic parameters, i.e., surface roughness length zorn and excess resistance κB-1. Multiple parameterization schemes have been used in the study to obtain values for surface roughness length and excess resistance κB-1 at the two stations. Results indicate that Zorn has distinct seasonal and inter-annual variability. For the type of surface with low-height vegetation, there is a large difference between the default value of Zorn in the land surface model and that obtained from this study, κB-1 demonstrates a significant diurnal variation and seasonal variability. Using the modified scheme for the estimation of Zom and κB-1 in the land surface model, it is found that simulations of sensible heat flux over the semi-arid region have been greatly improved. These results suggest that it is necessary to further evaluate the default values of various parameters used in land surface models based on field measurements. The approach to combine field measurements with atmospheric boundary layer theory to retrieve realistic values for key parameters in land surface models presents a great potential in the improvement of modeling studies of surface-atmosphere interaction.

Numerical Investigation of Influence of Rotor/Stator Interaction on Blade Boundary Layer Flow in a Low Speed Compressor

Numerical method was applied to the unsteady flow simulation at the mid span of a two-stage low speed compressor,and the blade boundary layer flow under rotor/stator interaction was investigated.By the model of wake/boundary layer interaction provided in this paper,the simulated blade frictional force and the boundary layer turbulent kinetic energy,the influence of wake/potential flow interaction on the blade boundary layer flow was analyzed in detail.The results show that under the condition of rotor/stator interaction,the wake is able to induce the stator laminar boundary layer flow to develop into turbulent flow within a certain range of wake interaction.In the stator suction boundary layer,an undisturbed region occurs behind the rotor wake,which extends the laminar flow range,and the wake with high turbulent intensity has the capability to control the boundary layer separation under adverse pressure gradient.

Aerothermal characteristics of bleed slot in hypersonic flows

Two types of flow configurations with bleed their aerodynamic thermal loads and related in two-dimensional hypersonic flows flow structures at choked conditions. are numerically examined to investigate One is a turbulent boundary layer flow without shock impingement where the effects of the slot angle are discussed, and the other is shock wave boundary layer in- teractions where the effects of slot angle and slot location relative to shock impingement point are surveyed. A key separation is induced by bleed barrier shock on the upstream slot wall, resulting in a localized maximum heat flux at the reattachment point. For slanted slots, the dominating flow patterns are not much affected by the change in slot angle, but vary dramatically with slot location relative to the shock impingement point. Different flow structures are found in the case of normal slot, such as a flow pattern similar to typical Laval nozzle flow, the largest separation bubble which is almost independent of the shock position. Its larger detached distance results in 20% lower stagnation heat flux on the downstream slot corner, but with much wider area suffering from severe thermal loads. In spite of the complexity of the flow patterns, it is clearly revealed that the heat flux generally rises with the slot location moving downstream, and an increase in slot angle from 20° to 40° reduces 50% the heat flux peak at the reattachment point in the slot passage. The results further indicate that the bleed does not raise the heat flux around the slot for all cases except for the area around the downstream slot corner. Among all bleed configurations, the slot angle of 40° located slightly upstream of the incident shock is regarded as the best.