Aero-optic effect is investigated in supersonic/hypersonic flow. Flat plate, ramp and cave model in supersonic/hypersonic flow condition are studied with large eddy simulation (LES). Turbulent density pulse is predicted to verify coefficient of density pulse model-mixing length (ML) model. It indicates that LES method simulates turbulent density pulse well in flat plate boundary layer flow, oblique shock wave/boundary layer interaction and strong separation flow. Density pulse data obtained by LES is used to verify ML coefficient and set a quantitative result. The modified ML was used to predict aero-optical effect on missile successfully.

Heat transfer distribution is analyzed on a rough wall in high speed turbulent compressible flow via computational fluid dynamics(CFD) method and analytical correlations, focusing on heat transfer with different roughness number and roughness element shape. It shows that, in all cases, heat flux augmentation predicted with CFD increases with reduction of roughness element density and levels off after roughness shape density is small, which differs with data of three analytical correlations. Predicted heat fluxes are same if same roughness element density and equivalent height are imposed on analytical correlations distinguishing from tendency of CFD results, which changes with roughness element shapes.

Dose coefficients for external photon radiation are widely used for assessment of radiation dose to public and workers due to ground surface contamination and air immersion of radionuclides released in nuclear accidents. Dose coefficients based on Chinese reference voxel phantom were presented. Photons that incident into the cylinder, which is slightly larger than phantom, is simulated using Geant4, including distributions of angle, height and energy of photons; Secondly, photons are sampled on surface of cylinder around Chinese Reference Male/Female Voxel Phantom with MCNPX 2.4.k from the above source, dose coefficients of 20 initial gamma ray energies, from 15keV to 10 MeV, are calculated. Thirdly, dose coefficients of 68 important nuclides are evaluated by combining of decay data from ICRP 107 and photon coefficients using cubic-spline fitting. GB/T 17982-2000 shows an overestimation for ground contamination while an underestimation of high Z nuclides for air immersion compared to our results.

A three-dimensional computational fluid dynamics algorithm is developed to study chemical nonequilibrium flowfield in a standard model, ELECTRE, for which flight data exits. A finite rate chemistry model based on Dunn and Kang's model is implemented. An extensive investigation on catalysis modeling relevance simulations with non-catalytic wall conditions as well as with full catalytic boundary conditions is made. The method developed is used to compute detailed flow features of hypersonic flow around forebody of a sphere-cone vehicle under different attack angles and altitudes. It is shown that reacting gas numerical results are consistent with flight data. It shows that real gas effects are prominent in thin shock wave layers closing to the wall surface. It makes outstanding distance of shock wave short. Heat flux under full catalytic wall boundary condition is higher than that under non-catalytic condition. The greater the attack angle,the more obvious this variance and the less electron number density on the wall. The higher the flight altitude,the lower dissociation of oxygen and nitrogen,and the lower heat flux on stagnation point.

In hypersonic/supersonic flow conditions, numerical simulations, experience relations and experiment analysis are made to study thermal-environment distribution characteristic around local-parts imposed by protuberance interacting body. Comparing flow characteristic structures in different shape, height, breadth ere, we found that interaction-zone size induced by high protuberance depend sensitively on effective breadth and short protuberance interaction-zone size depend on protuberance height. On the other hand, interaction-zone size is imposed by shade of lead edge, whether laminar or turbulence. Square-lead-surface and arc lead surface separation length ratio is 1.414. Sweepback affects separation length. A correlative line is built for size of separation imposed by sweepback.

In productivity evaluation of infill horizontal well interference of oil or water wells are taken into account in a vertical and horizontal well mixed well pattern.With mirror-image method and source-sink theory,a semi-analytical model of productivity evaluation for infill horizontal well in random mixed pattern was built,in which interference of wells and pressure drop caused by wal friction and fluid acceleration were taken into account.It is suitable for infill horizontal well productivity calculation in arbitrary well pattern.

We establish a 2D diffusion mathematic model considering transverse diffusion and feature of channeling reservoir.With Laplace transform an analytical model is solved.A solution with small slug is found.With general-purpose finite element software in a geometrical model numerical solution is solved.2D concentration profile is ploted.It is shown that larger Peclet number leads to higher peak concentration and later break through time.With combination of analytic solution and numerical solution,it is helpful in understanding the essential of mass transfer diffusion.

We simulate water adsorption on PuO₂ surfaces with kinetic Monte Carlo.Desorption activity energies of layers are obtained by fitting experimental data of Statebake and Haschke et al.We predict TDMS of water at different temperature increasing rates,adsorption isothermal and adsorption isotonic curves at different temperatures and water partial pressures.

Collins' formula and its inverse calculation are expressed in a convolution form. Corresponding transfer function is deduced. Sampling theory conditions in fast Fourier transform (FFT) to calculate Collins' formula are discussed. Diffraction field as optic wave passes through an optical system and reconstruction of incident plane are shown.

An algorithm is developed for the flow field in a two-phase cavitation region.The algorithm is implemented in a viscous calculation based on a Reynolds-Averaged Navier-Stokes(RANS) equation solver considering the effect of turbulence.A prior solution of wall detachment point and bubble length is not required.They are simulated by a cavitation model.The liquid/vapor interface is tracked and obtained by an iteration of flow field and interface updating.The algorithm is well validated with a comparison of computational and experimental data available for external flows of a cone/cylinder body and a hemisphere/cylinder body.The liquid/vapor interface shape and bubble length are calculated.The pressure coefficient distributions obtained along the cavitation bubble surface agree well with experimental results.The feasibility and accuracy of the algorithm are illustrated.

To solve the three-dimensional transient problems,the Improved Quasi Static Method(IQS) is adopted to deal with the temporal problem,and an improvement of synthesis method(YN) is also introduced.By factorizing the neutron flux,the time step of this method is enlarged so as to decrease the number of spatial calculations.For spatial calculations,Nodal Green's Function Method(NGFM) is used to determine the distribution of shape function,so the spatial mesh of IQSYN/NGFM can be about twenty times as large as that of finite difference method.