Flow and heat transfer in porous media with settlement distributed porosity is studied. Based on ‘O-shape’ ring theory and in situ measurement, a distribution function of porosity in x direction is established. Considering change of fluid density, Brinkman-Forchheimer-extended Darcy model is introduced. Energy equation is solved with continuity boundary conditions. A model of flow and heat transfer in porous media with settlement distributed porosity is established. The model is discretized with difference method. Gauss-Seidel iteration method is used. Numerical analysis shows that:Fluid flow velocity in porous media with of settlement distributed porosity is greater near wall surface, and it is small at center. Increase of porosity near wall surface reduces low velocity zone and increases high velocity zone. As porosity is small, temperature decreases linearly. As porosity is large, temperature near high and low temperature wall surface reduces rapidly. The decrease becomes slower at center. Heat is transferred through conduction and convection together. Increase of porosity increases average Nusselt number and enhances heat convection.

A PIC/MC model is developed for electrons,atomic ions N⁺ and molecular ions N₂⁺ in a capacitively coupled rf glow discharge.Modeling results are compared with those in a dc glow discharge at similar conditions.It reveals that the mean densities of particles(e,N₂⁺,N⁺) in rf discharge are nearly ten times greater than those in dc discharge,the mean energy near a rf electrode is about three times less than that near a cathode in dc discharge;the ions N⁺ have higher energy near two electrodes,the ions N₂⁺ have higher density,and the density of N₂⁺ is about six times the density of N⁺.The simulated electron energy probability distributions agree with measured results.

In simulation of a charged particle's guiding-center orbit in Tokamak equilibrium configurations, an efficient way to build physical functions such as magnetic field was obtained. Discrete Fourier transform is applied to replace cubic spline interpolation and B spline interpolation methods in calculation of gulding-center's trajectory. It has same accuracy as cubic spline interpolation and facilitates calculations in tearing mode as well. The method holds naturally property of symmetry. Derivatives of 3rd order or above can be obtained easily. Moreover, effect of incorrect data is reduced to the minimal limit compared with cubic spline interpolation.

A numerical investigation is performed to study aerodynamics of a 2D NACA0012 wing and flows in ground effect (IGE). Lift and drag forces are obtained and pressure distribution on airfoil surfaces is recorded at different angles of attack and ground clearances. Viscous flow near ground and air compressibility are taken into account. High lift is obtained as the airfoil is in dose proximity to the ground. The stall angle decreases with reduction of the ground clearance, due to higher adverse pressure gradient. Viscous effect of floor shows less effect on airfoil. Compressibility should be taken into account at lower ground clearances.

Resonance absorption of p-polarized electromagnetic wave by mode conversion in magnetized inhomogeneous plasma is investigated numerically as external magnetic field is perpendicular to density gradient.Different from an unmagnetized inhomogeneous plasma,the mode conversion coefficient in a magnetized inhomogeneous plasma depends on a cooperation of the incident direction,the density gradient direction and the direction of external magnetic fied.

For the interpretation of Azimuthal Resistivity Imager (ARI) in quantity,a modeling code using finite element method (FEM) is developed to simulate ARI responses in three dimensional inhomogeneous formation models.The code has been benchmarked against analytical solutions and 2D finite difference method (FDM) for subsets of the 3D geometry.The differences between the present results and two others are within 5%.The code has been also used to research some typical formation models in 3D environments.The results show it can provide support for the interpretation of ARI in quantity.