NUMERICAL METHOD OF MIXED FINITE VOLUME-MODIFIED UPWIND FRACTIONAL STEP DIFFERENCE FOR THREE-DIMENSIONAL SEMICONDUCTOR DEVICE TRANSIENT BEHAVIOR PROBLEMS

Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.

The dimension split element-free Galerkin method for three-dimensional potential problems

This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.

Two kinds of contact problems in three-dimensional icosahedral quasicrystals

Two kinds of contact problems, i.e., the frictional contact problem and the adhesive contact problem, in three-dimensional （3D） icosahedral quasicrystals are dis- cussed by a complex variable function method. For the frictional contact problem, the contact stress exhibits power singularities at the edge of the contact zone. For the adhe- sive contact problem, the contact stress exhibits oscillatory singularities at the edge of the contact zone. The numerical examples show that for the two kinds of contact problems, the contact stress exhibits singularities, and reaches the maximum value at the edge of the contact zone. The phonon-phason coupling constant has a significant effect on the contact stress intensity, while has little impact on the contact stress distribution regu- lation. The results are consistent with those of the classical elastic materials when the phonon-phason coupling constant is 0. For the adhesive contact problem, the indentation force has positive correlation with the contact displacement, but the phonon-phason cou- pling constant impact is barely perceptible. The validity of the conclusions is verified.

Recursive super-convergence computation for multi-dimensional problems via one-dimensional element energy pro jection technique

This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy projection （EEP） technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element （FE） model is reached. This conceptual dimension-by- dimension （D-by-D） discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional （2D） and three-dimensional （3D） problems can be obtained over the domain. Numerical examples of 3D Poisson＇s equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.

A Dimension-Splitting Variational Multiscale Element-Free Galerkin Method for Three-Dimensional Singularly Perturbed Convection-Diffusion Problems

By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.

TWO-DIMENSIONAL RIEMANN PROBLEMS:FROM SCALAR CONSERVATION LAWS TO COMPRESSIBLE EULER EQUATIONS

In this paper we survey the authors＇ and related work on two-dimensional Riemann problems for hyperbolic conservation laws, mainly those related to the compressible Euler equations in gas dynamics. It contains four sections： 1. Historical review. 2. Scalar conservation laws. 3. Euler equations. 4. Simplified models.

Solution Method for Singular Initial Value Problems of One-dimensional Steady Transonic Flow in a Dual-mode Scramjet

Singular initial value problems arise in solving one-dimensional steady transonic flow of dualmode scramjet. The existing solution method has the problems of large initial value errors in principles. This paper puts forward an improved algorithm based on variable transformation, and constructs a nonsingular one-dimensional steady transonic flow equation by defining a new variable. The improved algorithm can eliminate the singularity of the differential equation, and can solve the singular initial value problems of one-dimensional steady transonic flow of dual-mode scramjet.

Boundary Value Problems for Nonlinear Elliptic Equations of Second Order in High Dimensional Domains

This paper mainly concerns oblique derivative problems for nonlinear nondivergent elliptic equations of second order with measurable coefficients in a multiply connected domain. Under certain condition, we derive a priori estimates of solutions. By using these estimates and the fixed-point theorem, we prove the existence of solutions.

A branch-and-bound algorithm for multi-dimensional quadratic 0-1 knapsack problems

In this paper, a branch-and-bound method for solving multi-dimensional quadratic 0-1 knapsack problems was studied. The method was based on the Lagrangian relaxation and the surrogate constraint technique for finding feasible solutions. The Lagrangian relaxations were solved with the maximum-flow algorithm and the Lagrangian bounds was determined with the outer approximation method. Computational results show the efficiency of the proposed method for multi-dimensional quadratic 0-1 knapsack problems.

Eigenfunction expansion method and its application to two-dimensional elasticity problems based on stress formulation

This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.

Exponential Convergence of Finite-dimensional Approximations to Linear Bond-based Peridynamic Boundary Value Problems

In this paper, we demonstrate that the finite-dimensional approximations to the solutions of a linear bond-based peridynamic boundary value problem converge to the exact solution exponentially with the analyticity assumption of the forcing term, therefore greatly improve the convergence rate derived in literature.

Existence Results for m-Point Boundary Value Problem at Resonance with One-dimensional p-Laplacian

By using Mawhin＇s continuation theorem, the existence of a solution for a class of m-point boundary value problem at resonance with one-dimensional p-Laplacian is obtained. An example is given to demonstrate the main result of this paper.

Solution of two-dimensional scattering problem in piezoelectric/piezomagnetic media using a polarization method

Using a polarization method, the scattering problem for a two-dimensional inclusion embedded in infinite piezoelectric/piezomagnetic matrices is investigated. To achieve the purpose, the polarization method for a two-dimensional piezoelectric/piezomagnetic ＂comparison body＂ is formulated. For simple harmonic motion, kernel of the polarization method reduces to a 2-D time-harmonic Green＇s function, which is obtained using the Radon transform. The expression is further simplified under conditions of low frequency of the incident wave and small diameter of the inclusion. Some analytical expressions are obtained. The analytical solutions for generalized piezoelectric/piezomagnetic anisotropic composites are given followed by simplified results for piezoelectric composites. Based on the latter results, two numerical results are provided for an elliptical cylindrical inclusion in a PZT-5H-matrix, showing the effect of different factors including size, shape, material properties, and piezoelectricity on the scattering cross-section.

Eigenfunction expansion method of upper triangular operator matrixand application to two-dimensional elasticity problems based onstress formulation

This paper studies the eigenfunction expansion method to solve the two dimensional （2D） elasticity problems based on the stress formulation. The fundamental system of partial differential equations of the 2D problems is rewritten as an upper tri angular differential system based on the known results, and then the associated upper triangular operator matrix matrix is obtained. By further research, the two simpler com plete orthogonal systems of eigenfunctions in some space are obtained, which belong to the two block operators arising in the operator matrix. Then, a more simple and conve nient general solution to the 2D problem is given by the eigenfunction expansion method. Furthermore, the boundary conditions for the 2D problem, which can be solved by this method, are indicated. Finally, the validity of the obtained results is verified by a specific example.

Exact analytical solution to three-dimensional phase change heat transfer problems in biological tissues subject to freezing

Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.

THE RIEMANN PROBLEM FOR A TWO-DIMENSIONAL HYPERBOLIC SYSTEM OF CONSERVATION LAWS WITH NON-CLASSICAL SHOCK WAVES

The Riemann problem for a two-dimensional 2 x 2 nonstrictly hyperbolic system of nonlinear conservation laws has been solved thoroughly for any given initial data which are constant in each quadrant. The non-classical shockwaves, which are labelled as delta-shock waves, appear in some solutions. The solutions have been obtained are not unique. Due to the specific property of the system considered, there are no rarefaction waves in solution. This paper is divided into three parts. The first part constructs Riemann solutions for initial data involving two contact discontinuities while the second considers the case for other initial data. The last part briefly discusses the non-uniqueness of the solutions.

2.5维电阻率成像的新方法

以作者提出的正演方法为基础，给出一种新的２．５维电阻率成像方法，节省了计算时间，使成像迭代过程很稳定，可以对实际上广泛存在的大电性对比度的介质成像，从而提高了电阻率层析成像技术的实用性。

2.5维可控源音频大地电磁法Occam反演理论及应用

Occam反演是一种正则化的光滑模型反演方法,它在寻找最小拟合差的同时追求最光滑模型.因地下介质的电性通常是连续变化的,Occam反演考虑了地下介质横向和纵向的光滑情况,应用Occam反演避免了模型电性参数被随意间断.通过建立层状和存在两个异常体的模型,利用2.5维可控源音频大地电磁有限元法作正演响应计算,对响应数据基于Occam反演理论拟合反演.结果表明,Occam反演对可控源音频大地电磁法响应数据反演稳定收敛,理论模型与反演结果的视电阻率曲线形态基本一致,测点数据的拟合误差小于3%,反演结果能反映出模型电阻率的光滑效果.

2.5维复电阻率反演并行算法设计

近年来2.5维复电阻率反演取得了一些理论和应用方面的成果,但计算过程复杂、计算量巨大,特别是正演计算过程中计算量大导致了反演速度较慢。这里在2.5维复电阻率反演计算基础上,将反演算法中的正演过程的排列和频率循环分割成多个规模较小的程序段,并发执行正演过程中的每个排列和频率循环的求解计算过程,并根据并行计算的节点数目或者计算核心数目,动态调整每个节点或者计算核心的计算量,从而提高了算法的计算效率和适应性。

基于WMTS标准的2.5维电子地图的构建

2.5维数字地图既具有三维数字地图良好的表现效果,又具备数据量小、网上传输速度快的优点,满足了在现有的网络环境下对三维地图的快速浏览访问,在传统二维地图和真三维地图之间起到了很好的承上启下的作用。本文针对阜新市二维矢量数据和属性数据,利用3ds Max软件对研究区域建模并渲染成图,然后,在PhotoShop软件中进行图幅拼接,最后制作成瓦片地图与矢量图叠加发布,以便实现在Web端对2.5维电子地图的浏览、查询和空间分析等功能,并重点研究了基于OGC的WMTS标准和瓦片金字塔技术的标准化切图制作、发布、调用方法。