Development situations of domestic and overseas computational physics are discussed, such as society organizations, series of intermational conferences, corresponding institutes, university education and new developments in this field.From the above-mentioned developments, it is known that computational physics will have a flourishing development and a nice prospect in the 21st century.

In this paper, we are devoted to a study of the effect of the treatment of Neumann boundary conditions for one-dimensional advection-diffusion equation by the analysis method, and for tow-dimensional Poisson equation by the numerical test method. The reults show that the scheme is first-order accurate in space (i.e. O(h)can be derived first-order accurate difference solution or second-order solution, and the O(h²) or O(h³) scheme is only derived second-order solution.

This paper is in order to solve the numerical computation of the periodic orbits of nonlinear autonomous differential equations. We design an optimal control method of the numerical computation. The solution method is basically a continue shooting method,under the guardianship of nonlinear least squaves that takes advantage of an optimal Central formulation, in which attached to the contral techniques of the penalty function.In this method,we use the Crank-Nicolson discretization schemes and Conjugate Gradient methods.This method has been applied to the solution of nonlinear differential equations modelling chemical reactions,the corresponding numerical results are presented here.

In this paper, we have designed the variational optimal algorithm, and used it to compute simultaneously the limit cycle and period of the Brussela-tor.

We take 15 M_☉ model of Weaver et al. as initial model,the collapse process of stellar core between the onset of collapse and core bounce has been calculated by using the eouation of state of us (EOSⅡ).With the simplified nuclear statistical equilibrium system mentioned above, the treatment of the thermal effects of neutronization, nuclear disintegration, nuclear combination and nuclear excitation has been attemped.The results of us are essential agreed with Van Riper's results.calculated by using the ILPR's equation of state for the same initial model.As the center density of core reaches 1014gcm-3, the maximum velocity of fall and the Ms (mass at sonic point)in model D is greater than that of the typical model of Van Riper by 2.7% and 7% respectively. In model D the neutrino trapping density is less than that of the typical model of Van Riper by 3% (at mass layer 0.15M_☉) but the lepton fraction in model D is greater than that of typical model of Van Riper by 8%(at 0.15M_☉) and by 10% (at 0.5M_☉).

In this paper, the properties of system consisting of various nuclear species and free nucleons are described by a nuclear statistical equilibrium system consisting of a typical heavy nucleus, a particles, free neutrons and protons.starting from a semiempirical expression for the mass of the heavy nucleus and the entropy (energy) of their excited states, an equation of state is developed by a restricted minimization of the total free energy of a fixed volume of matter at given temperature and electron fraction.his scheme gave good agreement with the equation of El Eid & Hillebrandt at density 10⁹-10¹²g cm^-3 and with lamb et al (ILPR's equation of state at density 10¹²-10¹⁴g cm^-3,ye=3,σ/k=1.5k.The imput physics such as the rate of electron capture on matter has been treated also more carefully.

For the computation of saturated-unsaturated water table flow, a kind of difference scheme of irregular convex polygon mesh is developed. Using this method, saturated-unsaturated seepage flow is calculated to AKAI's sand model and KOM-ADA's reservoir model. The results obtained here are compared with results ob -tained from computations in which rectangled difference meshes were used. The applicableness and accuracy of this method are also discussed.

For saturated-unsaturated water table flow, computations of ADEP, ADIP and implicitly alternating direction iterative procedure are made and compared.The i-terations are not convergent and computational results are strongly distortional Based on supplementing a particular term, a kind of completely implicit and altern-ting direction iterative procedure is given The relative results for computational model of different scales are satisfactory. Problems of how to give the particular supplemental term and how to select better parameters are discussed, the computational results are given.

The implicit alternating-direction-iteration method is used to solve question of saturated-unsaturated seepage flow under rapid change of the water surface level. The numerical.results of Akai model are given and compared with experimental results.It is found that the numerical results given in this paper are in good agreement with experimental results and better than results by finite element method.The numerical results of komada model are close to results by finite element method.