With first-principles plane-wave ultrasoft pseudopotential method based on density functional theory (DFT), we studied acceptor level characreristics of anatase titanium dioxide doped and co-doped with 3d metal impurities Cu and non-metallic impurities C, N and F.It shows that Cu-N co-doping system and Cu, N single-doped system are better for visible light.Cu-N co-doping system, compared with Cu, N single-doped system, has a smaller band gap and greater density distribution on shallow acceptor level.It also shows that Cu-N co-doping system has highest absorption and reflectivity.Therefore, this system is the best for visible light.Cu and N acceptor level cooperative action results in the best visible light effect.

To study performance of graphic processing unit (GPU) for computational fluid dynamics, simulation on interactions between shock waves and a flame interface, a typical compressible reactive flow, was carried out on CPU/GPU heterogeneous system. Several optimal strategies are taken to raise GPU code performance. Computational results and acceleration performance of GPU with different grid number were examined. It was found that computational results by parallel GPU are the same as those by traditional CPU based on MPI parallel of 8 threads. Computational times of two parallel methods linearly grow with increase of computational grid numbers. Compuational cost by GPU is less than cost by MPI. As grid number is small(1.6×10⁴), speedup ratio of 8. 6 was achieved on GPU. As grid number grows, speedup ratio decreases. However, a ratio of 5.9 still can be held on GPU when grid number is more (4.2×10⁶). Arithmetic on CPU/GPU heterogeneous system provides a good solution for large-scale computations of compressible reactive flows.

We develop a lattice Boltzmann model for compound Burgers-Korteweg-de Vries (cBKdV) equation. By properly treating dispersive term u_xxx and applying Chapman-Enskog expansion, the governing equation is recovered correctly from lattice Boltzmann equation and local equilibrium distribution functions are obtained. Numerical experiments show that our results agree well with exact solutions and have better numerical accuracy compared with previous numerical results. This hence indicates that the model is satisfactory and efficient.

Since a quantum system is inevitably influenced by some unpredictable perturbations, we thereby conclude that all the experimental realizations of Grover quantum search algorithm reported were, in fact, achieved in a three-dimensional complex subspace. We also prove that in a two-dimensional complex subspace, for any given initial superposition of basis states|γ₀>=cosβ₀|α>+sinβ₀e^iζ|β)(β₀ is a small positive real number, ζ is an arbitrary real number), there exists a set of solutions F_j={(θ_j,θ_j-1,…,θ₁),(φ_j,φ_j-1,…,φ₁)} such that a desired state can be found with certainty for some positive integer j≥2, where the phase rotation angles θ_l andθ_t are real numbers but not equal to 2k'π,1 ≤ 1 ≤ j,k'is an arbitrary integer. If it is only required that a desired state can be found with high success probability, then as the total number of the desired and undesired states in an unsorted database is sufficiently large the above set of solutions F_j can be written in the form  for a relatively small positive integer j.

A high accuracy self-adaptive hp-finite element method(FEM) is used to simulate resistivity logging-while-drilling(LWD) electric field distribution.The algorithm can choose optimal refinement and calculation strategy according to practical formation model and error estimation.Accurate approximations of electrical field can be obtained without an accurate solution in the entire computational domain.Self-adaptive hp-FEM method has an exponential convergence rate.It indicates that voltage phase difference and amplitude ratio curves between receiver coils are well matched with the layer model.It is helpful to interpretate resistivity LWD data.

Ant colony optimization(ACO),a global optimization method,is proposed to analyze ground state energy of quantum mechanical systems.It simulates the way that real ants find a shortest path from nest to food source and back.In order to reduce iterations and improve solutions,ACO algorithm for ground state energy is modified.The proposed method exhibits advantage compared with traditional variation method.Ground state energy of helium atom and hydrogenic donors in GaAs-(Ga,Al)As quantum dot are calculated.The algorithm is demonstrated via comparison with variational method and genetic algorithms (GAs).

A v'²-f model,a k-ε model and an indoor zero equation model are adopted respectively to simulate mixed convections which are typical in ventilation and air conditioning systems.It indicates that there is a separated vortex near the ventilation inlet,and the accuracy for predicting this vortex is a critical factor in the simulation.The separated vortex dimension obtained by experiments is between the values calculated by the v'²-f model and the indoor zero equation model,while the dimension predicted by the k-ε model is a bit smaller.The results predicted by the indoor zero equation fit experimental results best.However,it is not in general from the view of model structure.For calculating temperature field the v'²-f model is better than the indoor zero equation model,and the k-ε model is the worst among the three.

A model and a simulating method for calculating the proton Single Event Upset (SEU) cross section are presented by the analysis of proton reactions in silicon. The energies deposited by protons of various energies are calculated in sensitive volume in memory cell. It is pro posed that high energy proton SEU is greatly attributed to the recoils (heavy ions), produced by proton reactions in silicon, which deposit ener gy and induce charge in sensitive volume. The relations of the SEU cross section vs. proton energy and critical charge have been obtained. The calculation results are in agreement with the experimental data.

The numerical computation is often used to predict or study pool fires.The fluid flow,heat and mass transfer,chemical reaction and their interaction in the burning process are studied by formulating and solving a set of governing equations.The treatment of oil surface is given in details also.Many results of pool fires under various wind speeds are gained,and computational results are obtained.

A new model is established to perform simulation for Single-Event Gate-Rupture of power MOSFETs in use of PSPICE circuit simulation software. The application results have a very good agreement with the corresponding data in published articles.

A simple three-parameter dynamical model in nonlinear quantum theory-the quant ized discrete self trapping equation is used to calculate the highly excited vi brational spectrum of CH stretches of CH₄ molecule with up to total quantum number n=7. The results show the model calculations appear to describe the expe rimental data well.

It investigates the application of solid crystallization process in the optimization algorithm.A special way is also proposed for the parameter "stirring" for the target function to get minimum value of the target function in the whole parameters' region.Meanwhile,the optimization algorithm is independent on the initial value of parameters without the larger calculation time.

A hybrid computational method and the corresponding software for simulating oil well perforator are described.A good agreement has been obtained between computational results and experiments.This research indicates that the software has a higher ratio of effectiveness over cost.

Using Mathematica symbalic manipulation system, the computer program MATHELDCA is developed for calculating one loop Feynman diagrams in renormalizable gauge theories. It is only needed to input the amplitudes of contributed diagrams according to the Feynman rule, so as to calculate the radiative corrections of particle process.

An numerical method was developed for measuring radiant temperature using a set of absorption-foils coupled to soft x ray streak camera, SXRSC, in laser plasma experiments. An numerical code, SCC, was given. An x-ray intensity ratio vs temperature was calculated for Parylene, C₃H₆, Formvar, Mylar and Aluminum. A suitable range of the measuring temperature was discussed to some absorber.

This Paper Presents an approximate method of calculating the ianization degrees, and then an analytic expression of the equation of state for high temperature gases is derived, which can be applied directly to the fluid dynamics calculations. The expression is simple and Its results are fairly accurate.