In this work, water adsorbed in single-walled carbon nanotube (SWCNT) is simulated by the Grand Canonical Monte Carlo (GCMC) method. The effects of pore length, pore size and surface strength on water adsorption behaviors at 298 K are systematically studied, via the characterization of isotherms, local density distribution and isosteric heat. There are hysteresis loops observed for water adsorbed and desorbed in SWCNTs with relatively large pore radii (PR=0.8, 1.0 nm), and the hysteresis loop disappears as the pore radius decreases to 0.55 nm. In addition, the water molecule packing manners are in the arrangement form of single chain, double helix chain and water clusters due to the pore size effects. When the pore length is in the range of 4~8 nm, the initial adsorption pressure becomes smaller and smaller as the pore length increases, but this effect rule gradually disappears with increasing pore length to 10 nm. Finally, the initial pressure of water adsorption decreases with surface strength; and when surface strengths are 20, 28, 32 K, the capillary evaporation phase transitions completed instantaneously. While the surface strength increases to 40 K, the capillary evaporation phase transition is shown gradually desorbed steps. As surface strengths increased from 20 K to 40 K, the isosteric adsorption heat at the pressure points of capillary condensation phase transition is 127.47, 117.98, 84.04, 59.16 kJ·mol^-1, respectively.

To explore motion characteristics of indoor respirable particles, Brownian force on particles is considered in motion probability model. Motion characteristics of 0.01 μm, 0.1 μm and 1 μm particles under conditions of up supply with up return and up supply with side return air forms were simulated with lattice Boltzmann method taking Re as 400, 1 000 and 2 000, respectively. It shows that range of particle spacial distribution increases with Re, and smaller particles affected more obvious by air turbulence and diffusion effect. Mean square displacement (MSD) of particles is inversely proportional to Re and diameter of particles. At same Re, MSD of particles is greater under up supply with side return air form, thus lower suspended particles and higher indoor air quality appears in up supply with side return air form.

At late stages of compressible transitional boundary layers over a flat-plate, ejection and sweep flows are investigated by means of direct numerical simulation with high order accuracy. Numerical results with high resolution clearly represent complex phenomena of ejection and sweep flows in the transition process of boundary layers. Close relationships with typical vortex structures and spike structures are revealed. Details of formation mechanisms of ejection and sweep flows are analyzed. Numerical results demonstrate experimental findings. Corresponding mechanisms are analyzed.

An identification method for node parameters and topologies of complex networks is investigated with network heterogeneously splitting in a driving-response model.Modeling of heterogeneously associated networks is presented.Two heterogeneously associated complex networks are disassembled and modeled.With LaSalle invariance theory,Gram matrix and driving-response model,adaptive observers for identifying node dynamical parameters and topological parameters of heterogeneously associated complex networks are obtained.Node parameters and coupling topology parameters can be identified on line using the proposed observers.Numerical simulations indicate that the proposed methods are effective.

In a model for interaction between one atmosphere uniform glow discharge plasma (OAUGDP)and air proposed by Shyy,electric field and electric forces acting on air as body force are obtained.The electric field is calculated with electric potential equation. The influence of OAUGDP on flows around a NACA0015 airfoil at low speed is investigated by solving N-S equation with a body force term.Plasma located upstream of the separation point damps the flow separation around the airfoil efficiently,while plasma located downstream of the separation point has little effect on the flow.It agrees with experimental observations. Effect of plasma on aerodynamics of airfoil and pressure distribution on airfoil surface are investigated.

An in-core fuel management code package HEX-ICFM is developed for Xi'an Pulsed Reactor.The orthogonal design model and the loading pattern optimization program HEX-ORTH are studied.The cell calculation is performed using WIMS-D/4.A two dimensional hexagonal geometry multigroup nodal theory code SIXTUS-2 is used for core diffusion calculation.The core physical parameters are calculated using HEX-ICFM for the first cycle of Xi'an Pulsed Reactor.With the fuel of end of cycle 3 and 30 fresh fuel rods the optimum loading patterns are calculated for the objective function max(K^BOC_eff)using HEX-ORTH.

Through comparing different results of simulating shock reflection flow over a forward step using different schemes including Roe and van Leer schemes,the advanced upwind scheme (AUSM+) is fully studied in terms of its resolution and efficiency for shock waves,expansion vaves,contact discontinuities and their interactions.Its advantages of low numerical dissipation and high resolution for discontinuities are demonstrated.

The critical temperature and the ground state fraction of weakly interacting Bose-Einstein condensation in a harmonic potential trap are calculated with numerical method. The result indicates that both are decreasing with the scattering length increasing. But, there is only a difference about 0.4% between it and that of ideal case.

The combined RK-Rosenbrock(CRKR) methods are presented for a partitioned systems of stiffly large differential equations.Nonstiff and stiff subsystems are integrated in parallel on two processors by an explicit RK method and a Rosenbrock method respectively.Their construction,convergence and numerical stability are studied,and the numerical simulation experiments are conducted on a personal computer and a parallel computer.

With some improvements,the model of Interior Boundary Layer in heating power of Sea wind.Can simulate not only the field of sea wind,but also the flow field of river valley.The numerical computational results are consistent with the values of observation.This paper also brings the discussion on practical application of the preconditional Conjugate Gradient Method (PCG) that is most frequantly used in Multi-level Methods to solve the linear algebra equations,which result from the numerical approximation of Interior Boundary Layer in heating power of sea wind.The Computational time has been shown reduced considenablely.

The theoretic model: there are three interfaces before the tunnel face,the medium between any adjacent interfaces is homogeneous and the wave velocity in one layer of medium may be different from ones in the others.We studied the methods of forecasting,layer by layer,and reversed the arrival time of wave which deduced by the known interfaces and wave velocities to the interfaces. Results are almost concordant with the fact.To the tunnel model and Yuntai Mountain tunnel face, we recorded waveforms respectively by ultrasonic detector and seismic detector. By analysis of severel of waves on record, we obtained the data on the arrival time of the reflection wave, and then reversed the data to the interfaces.The results are in good agreement the true conditions.

The relation between components, pressure, flow velocity and gain coefficient of CW OICL is theoretically analyzed by using a simple model throug digital computation. The operating region concerned is [O₂^*]/[I₂]=10~200,[O₂^*]/[O₂]=2~9,P=1~100 torr,U=5-100 m/s, The necessity of flow conservation equation for a CW OICL is detially discussed, It is shown that the available gain region will by enlarged proportionally to flow velocity while the peak gain value almost remains inveriable when flow velocity is increased.