A method of directly solving the heat conduction equation is used to numerically study a controllable thermal cloaking technology with adaptive heat source.We explore heat flow control method and cloaking effect,and deduce general solution of distribution of adaptive heat source in the thermal cloak region with two-dimensional arbitrary symmetrical section under uniform background. Simulation results show that for a cloak with arbitrary symmetrical section,the adaptive heat source makes disturbed temperature field turn to the background temperature field and achieves the purpose of thermal cloaking.

DUGKS is extended to model porous media flow at representative elementary volume scale combined with generalized porous media model. It is verified by several two-dimensional classical problems:Poiseuille flow, Couette flow and cavity flow. Effectiveness of DUGKS for porous media flow is tested and advantage of DUGKS in grid flexibility is demonstrated. A fracture system is modeled by DUGKS for porous media flow.

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.

With semi-classical approximation, thermodynamic properties of Fermi gas trapped in both gravity field and magnetic field are studied. By using theoretical analysis and numerical simulation, influence of gravity field on thermodynamic properties of the system in strong magnetic field is analyzed. It shows that, compared with the case of strong magnetic field only, gravity field makes the energy, chemical potential reduced. With rising temperature, influence of gravity field on chemical potential is gradually enlarged. There is a maximal influence of gravity field on heat capacity. Gravity field makes oscillation of heat capacity almost unchangeable while oscillation center of chemical potential shift down.

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.

We review physical properties and potential applications of two-dimensional boron and carbon based nanostructures such as graphene and BC₃ sheet adsorbed with metal atoms, non-metal atoms or molecules. Pristine graphene is semi-metallic with zero band gap, which is nonmagnetic and has weak spin orbit coupling (SOC). BC₃ sheet is a semiconductor with an indirect band gap. As graphene and BC₃ sheet are adsorbed with metal atoms, non-metal atoms or molecules, graphene systems may have a band gap at Dirac point, strong SOC, magnetic order and superconductivity. Applications in hydrogen storage may be realized in graphene and BC₃ systems. In addition, it has been demonstrated that graphene can be used as gas sensors with significantly high precision.

Effect of conduction delay on dynamical behavior of spiral waves is investigated in Greenberg-Hastings model.It shows that in some cases dynamical behavior of spiral waves is not impacted by conduction delay.Changes in wavelength and wave speed,breathing spiral wave,multi-armed spiral wave,coexistence of spiral wave and anti-spiral wave are observed.Physical mechanisms of these phenomena are briefly analyzed.

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 method is proposed for choosing electron cutoff energy,which is self-adaptive with cell size,background material and electron position.As residuary range of an electron defined by total stopping power is less than distance to the nearest cell boundary,the electron is killed.Compared with the CUT model used in MCNP program,efficiency and precision of self-adaptation method is higher.The most important is that cutoff energy of electron is selected automatically in program in all models,without consideration of energy of particles or dimension of model.

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.

A method with CFD technology is used for accurate prediction of temperature field of supersonic and hypersonic kinetic energy projectiles.It includes equilibrium flow governing equations and difference scheme,equilibrium flow Jacobian matrixs and equilibrium gas state equations.A real gas equilibrium air model is used to calculate temperature in vector flux split in a difference scheme.Thermal environment for typical high-speed kinetic energy projectile is validated.Rationality of the method is considered.Temperature field of a designed hypersonic kinetic energy projectile is simulated.Reliable numerical data for projectile aerodynamic design and aerodynamic heat protection are provided.

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 hybrid optimization algorithm is developed for solving ellipsometric measurement problem.This algorithm is based on the high ability of simulated annealing method in seeking optimum solution and the high computation efficiency of simple shape method.It seeks the global minimum solution of objective function and gets proper parameters of optical model.The computation efficiency is higher than that of simulated annealing algorithm.The selection of inversion controlling parameters is discussed to make this algorithm adaptive.

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.

With the classical Nucleation theory and the combination of Monte Carlo method, the two dimensional model of liquid phases separation for binary system is present. It consists of three parts of nucleation, diffusion and growth The some methods dealing with problem of surface energy in the growth of new phase are raised. The principal procedure is shown for the simulation of phases separation with computer.

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.

The differental equation for the relationship of temperature distribution in solid, interface distance and time during cooling and solidificatian Processing of metal drop has been solved in one dimensional equation of heat conduction in spheriodal coordinate, and the relationships of tempertuer distribution in liquid and time during liquid cooling, and of solidification time at invariant temperature and spherical diameter are inferenced from the solution.

This paper is concerned with the effect of wall-temperature on heat transfer,drag and fore-stagnation-point pressure of a cylinder as it is immersed.into-a thermal plasma cross-flow. Within the wall temperature region from 1000 to 3000K, the 2-D governing equations including actual argon plasma properties are solved numerically by using the SIMPLER, algorithm. Computational results show that the wall temperature has only little effect on the heat flux, the drag and the fore-stagnation-point pressure, although gas properties as well as temperature and velocity gradients at the yclinder sur-face vary significantly.

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.