An algorithm and implementation of a dual-percolation model based on bond-percolation theory was introduced which combines pore-network with fracture-network. Monte-Carlo numerical simulation was conducted based on statistical analysis of key parameters,i. e. fracture population densities,Fracture angles,fracture lengths,and pore connectivity possibilities. Fractal dimension was presented to quantify fractal connectivity. Numerical results were in agreement with those of Balberg and Berkowitz et al. Finally,appropriate sample scale was discussed considering both calculation precision and computing time.

The finite element method is used to obtain numerical solutions of the Schrödinger equation for the Helium 1s2p-state. The relative errors of approximate energies are 10^-6 for triplet and 10^-4 for singlet which are slightly smaller than J. Shertzer's results for the Helium ground state[9]. The generalized eigenvalue problems obtained by FEM are symmetric for the ground state but unsymmetric for 1s2p-state. It becomes more difficult to solve the problems. Form the graphs of wave functions, it can be seen that it is reasonable to solve the Schrödinger equations in bounded domains.

An approach to the parallelization is proposed on the network of workstation,that simulates the produce of the thin film deposition with a scale of 100×100×100 atomic sites in three-dimension by Monte Carlo method.It uses a efficient parallel algorithm for the thin film deposition simulation using domain decomposition and asynchronous parallel communication,and adopts the reasonable division of domain with the topology geometric of space filling in the film deposition process.The results show that parallel processing generated by these methods can make communication overhead fallen off considerably,and the parallel performance of algorithm improved and computation time for the thin film deposition decreased largely.