The modified method of characteristics with mixed finite element domain decomposition procedures for the transient behavior of a semiconductor device

For the transient behavior of a semiconductor device, the modified method of characteristics with mixed finite element domain decomposition procedures applicable to parallel arithmetic is put forward. The electric potential equation is described by the mixed finite element method, and the electric, hole concentration and heat conduction equations are treated by the modified method of characteristics finite element domain decomposition methods. Some techniques, such as calculus of variations, domain decomposition, characteristic method, energy method, negative norm estimate and prior estimates and techniques are employed. Optimal order estimates in L² norm are derived for the error in the approximation solution. Thus the well‐known theoretical problem has been thoroughly and completely solved.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 353–368 2012     

渗流耦合系统动边值问题特征分数步差分方法

对多层渗流方程耦合系统动边值问题,提出适合并行计算的一类特征分数步差分格式,利用区域变换、变分形式、粗细网格配套、乘积型叁二次插值、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧,得到收敛性的最佳阶l2误差估计．该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中,得到了很好的数值模拟结果．     

THEORY AND APPLICATION OF FRACTIONAL STEP CHARACTERISTIC FINITE DIFFERENCE METHOD IN NUMERICAL SIMULATION OF SECOND ORDER ENHANCED OIL PRODUCTION

A kind of second-order implicit fractional step characteristic finite difference method is presented in this paper for the numerically simulation coupled system of enhanced(chemical) oil production in porous media.Some techniques,such as the calculus of variations,energy analysis method,commutativity of the products of difference operators,decomposition of high-order difference operators and the theory of a priori estimates are introduced and an optimal order error estimates in l~2 norm is derived.This method has been applied successfully to the numerical simulation of enhanced oil production in actual oilfields,and the simulation results are quite interesting and satisfactory.     

多层非线性渗流耦合系统的特征分数步差分方法

对多层非线性渗流耦合系统提出适合并行计算的特征分数步差分格式, 利用变分形式、能量方法、粗细网格配套、分片双二次插值、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶的l₂误差估计. 该方法已成功的应用到多层油资源评估的生产实际中.     

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain are given. These procedures use implicit Galerkin method in the subdomains and simple explicit flux calculation on the interdomain boundaries by integral mean method or extrapolation method to predict the inner‐boundary conditions. Thus, the parallelism can be achieved by these procedures. These procedures are conservative both in the subdomains and across interboundaries. The explicit nature of the flux prediction induces a time‐step limitation that is necessary to preserve stability, but this constraint is less severe than that for a fully explicit method. L²‐norm error estimates are derived for these procedures. Compared with the work of Dawson and Dupont [Math Comp 58 (1992), 21–35], these L²‐norm error estimates avoid the loss of H^?1/2 factor. Experimental results are presented to confirm the theoretical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Nitsche type mortaring for elliptic problems with corner singularities

The paper deals with Nitsche type mortaring as a finite element method (FEM) for treating non‐matching meshes of triangles at the interface of some domain decomposition. The approach is applied to the Poisson equation with Dirichlet conditions for the case that the interface passes re‐entrant corners of the domain and local mesh refinement is applied. Some properties of the finite element scheme and error estimates in a discrete H¹‐like and in the L₂‐norm are proved.     

三维多组分可压缩驱动问题的分数步特征差分方法

本文提出三维多组分驱动问题的分数步长特征差分格式,并应用变分形式,能量方法,粗细网格配套,叁二次插值,高阶差分算子的分解和乘积交换性理论和技巧,得到最佳阶L^2误差估计,该方法已成功应用到油资源评估,强化采油数值模拟和海水入侵预测和防治的数值模拟中。     

三维渗流耦合系统动边值问题迎风差分方法的理论和应用

对多层渗流方程耦合系统动边值问题,提出适合并行计算的两类迎风差分格式,利用区域变换、变分形式、能量方法、隐显格式的相互结合、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧,得到收敛性的l~2误差估计.该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中,得到了很好的数值模拟结果.     

The characteristic finite difference fractional steps methods for compressible two-phase displacement problem

For compressible two-phase displacement problem, a kind of characteristic finite difference fractional steps schemes is put forward and thick and thin grids are used to form a complete set. Some techniques, such as piecewise biquadratic interpolation, of calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates inL ² norm are derived to determine the error in the approximate solution. Project supported by the National Scaling Program, the National Tackling Key Problems Program and the Doctorate Foundation of the State Education Commission of China.     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

FINITE DIFFERENCE FRACTIONAL STEP METHODS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

Characteristic finite difference fractional step schemes are put forward. The electric potential equation is described by a seven-point finite difference scheme, and the electron and hole concentration equations are treated by a kind of characteristic finite difference fractional step methods. The temperature equation is described by a fractional step method. Thick and thin grids are made use of to form a complete set. Piecewise threefold quadratic interpolation, symmetrical extension, calculus of variations, commutativity of operator product, decomposition of high order difference operators and prior estimates are also made use of. Optimal order estimates in l2 norm are derived to determine the error of the approximate solution. The well-known problem is thorongley and completely solred.     

结构三角网上抛物方程的有限差分区域分解算法

1引言对于大型科学与工程计算问题,并行计算是必需的．构造高效率的数值并行方法一直是人们关心的问题,并且已有了大量的研究．在三层交替计算方法的研究中出现了许多既具有明显并行性又绝对稳定的差分格式(见[1]-[5])．在只涉及两个时间层的算法研究中,Dawson等人(见[6])首先发展了求解一维热传导方程的区域分解算法,并将其推广到     

Domain decomposition methods for hyperbolic problems

In this paper a method is developed for solving hyperbolic initial boundary value problems in one space dimension using domain decomposition, which can be extended to problems in several space dimensions. We minimize a functional which is the sum of squares of the L ² norms of the residuals and a term which is the sum of the squares of the L ² norms of the jumps in the function across interdomain boundaries. To make the problem well posed the interdomain boundaries are made to move back and forth at alternate time steps with sufficiently high speed. We construct parallel preconditioners and obtain error estimates for the method. The Schwarz waveform relaxation method is often employed to solve hyperbolic problems using domain decomposition but this technique faces difficulties if the system becomes characteristic at the inter-element boundaries. By making the inter-element boundaries move faster than the fastest wave speed associated with the hyperbolic system we are able to overcome this problem.     

非线性渗流耦合系统动边值问题二阶迎风分数步差分方法

对多层非线性渗流方程耦合系统三维动边值问题, 提出适合并行计算的一类二阶迎风分数步差分格式, 利用区域变换、变分形式、能量方法、隐显格式的相互结合、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶l² 误差估计. 该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中, 得到了很好的数值模拟结果.     

多孔介质可混溶流动问题的混合元与区域分裂特征差分格式

1 引 言 用区域分裂方法求微分方程的数值解,是近年来计算数学领域的—个新方法。这种方法通过分裂区域来减少所处理问题的规模,并实现并行计算,因此,特别适用于大范围的工程技术问题和数学物理问题。本文用这种方法处理平面可混溶不可压缩流动问题其中J=[0.T].u为Darcy速度,p为压力,c为浓度,k为渗透率,μ（c)为流体粘性,c为注入井给定浓度,q为外界源汇项,q~+=max{q,0}.φ为孔隙度。D为扩散矩阵,本文D与u无关。即仅考虑分子扩散,边条件可取为第一或第二类边条件。本文考虑第一类初边值问题     

����������������������ŷ���ģ�͵��Ʋ�ǰ���ӯ��ֲ�

本文综述混合效应模型参数估计方面的若干新进展. 平衡混合效应方差分析模型的协方差阵具有一定结构. 对这类模型, 文献[1]提出了参数估计的一种新方法, 称为谱分解法. 新方法的突出特点是, 能同时给出固定效应和方差分量的估计, 前者是线性的, 后者是二次的,且相互独立. 而后, 文献[2--9]证明了谱分解估计的进一步的统计性质, 同时给出了协方差阵对应的估计, 它不仅是正定阵, 而且可获得它的风险函数, 这些文献还研究了谱分解估计与方差分析估计, 极大似然估计, 限制极大似然估计以及最小范数二次无偏估计的关系. 本文综述这一方向的部分研究成果, 并提出一些待进一步研究的问题.     

The upwind finite difference fractional steps method for combinatorial system of dynamics of fluids in porous media and its application

For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L ² norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.     

The unconditional stability and mass‐preserving S‐DDM scheme for solving parabolic equations in three dimensions

In this paper, the unconditional stability and mass‐preserving splitting domain decomposition method (S‐DDM) for solving three‐dimensional parabolic equations is analyzed. At each time step level, three steps (x‐direction, y‐direction, and z‐direction) are proposed to compute the solutions on each sub‐domains. The interface fluxes are first predicted by the semi‐implicit flux schemes. Second, the interior solutions and fluxes are computed by the splitting implicit solution and flux coupled schemes. Last, we recompute the interface fluxes by the explicit schemes. Due to the introduced z‐directional splitting and domain decomposition, the analysis of stability and convergence is scarcely evident and quite difficult. By some mathematical technique and auxiliary lemmas, we prove strictly our scheme meet unconditional stability and give the error estimates in L²‐norm. Numerical experiments are presented to illustrate the theoretical analysis.     

Nonoverlapping domain decomposition characteristic finite differences for three‐dimensional convection‐diffusion equations

One domain decomposition method modified with characteristic differences is presented for non‐periodic three‐dimensional equations by multiply‐type quadratic interpolation and variant time‐step technique. This method consists of reduced‐scale, two‐dimensional computation on subdomain interface boundaries and fully implicit subdomain computation in parallel. A computational algorithm is outlined and an error estimate in discrete l²‐ norm is established by introducing new inner products and norms. Finally, numerical examples are given to illustrate the theoretical results, efficiency and parallelism of this method. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 17‐37, 2012     

The finite difference method for the three-dimensional nonlinear coupled system of dynamics of fluids in porous media

For the three-dimensional coupled system of multilayer dynamics of fluids in porous media, the second-order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in l2 norm are derived to determine the error in the second-order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources.