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.

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.

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.

Guckenheimer structure of solution of Riemann problem with four pieces of constants in two space dimensions for scalar conservation laws

By using the generalized characteristic analysis method, the two-dimensional four-wave Riemann problem for scalar conservation laws, which is nonconvex along the y direction, was studied. Riemann solutions, which involve the Guckenheimer structure, were constructed.

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.

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.

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.

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.

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.

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.

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.

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 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.

基于子循环自适应串行交错时间匹配算法的油浸式变压器绕组瞬态温升计算

针对传统算法存在的步长调整及耦合点选择问题,该文提出基于子循环自适应串行交错(SCAS)的时间匹配算法。首先,在步长调整方面,通过加入自适应(ATS)-启发式(HTS)混合变步长算法,解决步长无法实时调整的问题;然后,针对耦合点选择问题,该文在传统常规串行交错(CSS)时间匹配算法及现有改进方案子循环常规串行交错(SCSS)的基础上,提出SCAS时间匹配算法,采用自适应方案确定计算时间窗及预定耦合点,并引入指数平滑法对其精度进行保证;最后,该文基于110 kV油浸式电力变压器建立二维瞬态单分区分匝绕组模型,并将所提算法与其他传统方法进行对比分析。结果表明:在精度方面,与传统CSS算法相比SCAS算法的流场及温度场最大平均相对误差均不超过0.45%,在可接受误差范围内;在效率方面,SCAS算法总体计算效率较CSS算法提高了近20倍。同时,为了说明所提算法的工程价值,该文搭建基于110 kV变压器绕组模型的温升实验平台,将SCAS时间匹配算法应用于该模型中,实验结果表明:在精度方面,SCAS算法与实验达到稳态时最大绝对误差为2.7 K;在计算效率上,本文所提算法的计算效率较传统算法提升约46倍,计算步数约为传统的1/47,减少了计算冗余。

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.

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.

一种求解高维优化问题的改进灰狼算法

为求解高维优化问题,提出基于反向学习和衰减因子的灰狼优化算法(grey wolf algorithm based on opposition learning and reduction factor,ORGWO).设计一种灰狼反向学习模型,模型考虑问题搜索边界信息和种群历史搜索信息,初始种群阶段增加反向学习,增强种群多样性.根据算法各个阶段不同特征引入衰减因子,平衡全局和局部勘探能力.选取8个高维函数和23个不同特征的优化函数对算法性能进行测试,进一步使用收敛性分析,寻优成功率,CPU时间,Wilcoxon秩和检验来评估改进算法,实验结果表明,ORGWO算法在求解高维问题上具有较好的精度,鲁棒性和更快的收敛速度.

三维二十面体准晶弹性半平面周期接触问题研究

借助平面弹性复变方法,研究了周期刚性压头作用下三维二十面体准晶弹性半平面的周期有限摩擦和周期粘结接触问题.从应力、位移分量的复变函数形式表达式出发,利用半平面Hilbert核积分公式、Plemelj公式,得到两类周期接触问题的解答.针对周期有限摩擦接触问题,得到了三类周期刚性压头(直水平、直倾斜和圆形基底)作用在准晶体半平面上接触应力的封闭解.针对半平面周期粘结接触问题,求得接触边界上作用周期尖劈形位移时接触应力的显式解答.当不考虑相位子场的作用时,文中求得的理论结果可退化到正交各向异性材料平面弹性周期接触问题的对应结果.数值算例说明准晶弹性常数对接触应力分布规律及大小的影响.本文结论可为分析准晶材料压痕实验结果及准晶材料性能提供一定理论依据.

基于改进粒子群算法的木材板材下料方法

木材板材在家具行业应用广泛,以绿色环保、节约能源为目的的木材板材优化下料已经成为研究的热点。木材板材下料优化问题属于二维矩形下料问题,是一种具有高度计算复杂性的问题。本研究主要针对单规格木材板材进行矩形零件下料问题,在木材板材长和宽都大于零件长和宽的情况下,通过建立二维下料的数学模型,采用标准粒子群算法、变邻域搜索算法、粒子群混合变邻域搜索算法分别进行求解,并以某企业的下料实例进行分析计算。首先,利用标准粒子群算法求解单规格板材下料问题;其次,利用变邻域搜索算法求解单规格板材下料问题。在获得局部最优解的基础上改变其邻域结构再进行局部搜索,找到另一个局部最优解,如此不断迭代,直到满足算法的终止条件,获得全局最优解;最后,利用粒子群变邻域搜索混合算法求解单规格板材下料问题。针对粒子群算法局部搜索能力较差、容易过早收敛的问题和具有较好包容性的特点,将变邻域搜索的思想融入粒子群算法中,使结果更加趋向全局最优。结果表明:粒子群变邻域搜索混合算法相比粒子群算法和变邻域算法效率都有显著提升,能显著提高该木材板材的利用率,增加企业经济效益。

基于LSTM的小井眼阵列感应测井响应预测

快速准确地获取测井响应曲线是进行后续测井反演的基础。为了对测井响应进行高效的正演建模,将LSTM神经网络应用于小井眼阵列感应测井响应的预测中。数据集来源于COMSOL软件建模仿真得到的小井眼感应测井响应数据。在离线训练阶段,LSTM模型从测井响应数据库中提取到小井眼感应测井响应的数据特征。给定简单各向同性层状地层和复杂各向同性地层两种地层参数,训练得到的模型可用于预测阵列感应仪器的小井眼测井响应。研究表明,LSTM预测结果与基于测井响应数据库传统多维插值法的结果基本吻合。在保证预测结果精度的前提下,LSTM预测极大地提高了不同地层场景下测井响应的预测时间,且具有良好的泛化能力。