A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative ...A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative position,velocity and attitude of two unmanned aerial vehicles (UAVs).The second-order divided difference filter which makes use of multidimensional interpolation formulations to approximate the nonlinear transformations could achieve more accurate estimation and faster convergence from inaccurate initial conditions than standard extended Kalman filter.The filter formulation is based on relative motion equations.The global attitude parameterization is given by quarternion,while a generalized three-dimensional attitude representation is used to define the local attitude error.Simulation results are shown to compare the performance of the second-order divided difference filter with a standard extended Kalman filter approach.展开更多
In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
In this paper, combining the idea of difference method and finite element method, we construct a difference scheme for a self-adjoint problem in conservation form. Its solution uniformly converges to that of the origi...In this paper, combining the idea of difference method and finite element method, we construct a difference scheme for a self-adjoint problem in conservation form. Its solution uniformly converges to that of the original differential equation problem with order h3.展开更多
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 conditi...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.展开更多
In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial deriv...In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.展开更多
For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of calculus of variations, multiplicative commutation r...For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of 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 2 norm are derived to determine the error, in the approximate solution.展开更多
Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploi...Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference 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 change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.展开更多
The research of the miscible oil and water displacement problem with moving boundary values is of great value to the history of oil-gas transport and accumulation in the basin evolution as well as to the rational eval...The research of the miscible oil and water displacement problem with moving boundary values is of great value to the history of oil-gas transport and accumulation in the basin evolution as well as to the rational evaluation in prospecting and exploiting oil-gas resources. The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values. For the twodimensional bounded region, the upwind finite difference schemes are proposed. Some techniques, such as the calculus of variations, the change of variables, and the theory of a priori estimates, are used. The optimal orderl2-norm estimates are derived for the errors in the approximate solutions. The research is important both theoretically and practically for the model analysis in the field, the model numerical method, and the software development.展开更多
In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the ...In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the discretization of p.d.e. via the preconditioned conjugate gradient method. For the second-order partial differential equations with Dirichlel boundary conditions, we prove that the condition number of the preconditioned system is O(1) while the condition number of the original system is O(m 2) Here m is the number of interior gridpoints in each direction. Such condition number produces a linear convergence rale.展开更多
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....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 ate quite interesting and satisfactory.展开更多
In this paper, we first provide a generalized difference method for the 2-dimensional Navier-Stokes equations by combing the ideas of staggered scheme m and generalized upwind scheme in space, and by backward Euler ti...In this paper, we first provide a generalized difference method for the 2-dimensional Navier-Stokes equations by combing the ideas of staggered scheme m and generalized upwind scheme in space, and by backward Euler time-stepping. Then we apply the abstract framework of to prove its long-time convergence. Finally, a numerical example for solving driven cavity flows is given.展开更多
The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechan...The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics （English Edition）, 2007, 28（7）, 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.展开更多
In this short article, the upwind and central compact finite difference schemes for spatial discretization of the first-order derivative are analyzed. Comparison of the schemes is provided and the best discretization ...In this short article, the upwind and central compact finite difference schemes for spatial discretization of the first-order derivative are analyzed. Comparison of the schemes is provided and the best discretization scheme for convection dominated problems is suggested.展开更多
Numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science.The mathematical model is defined by a nonlinear system of initial-boundary problem including four parti...Numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science.The mathematical model is defined by a nonlinear system of initial-boundary problem including four partial differential equations:an elliptic equation for electrostatic potential,two convection-diffusion equations for electron concentration and hole concentration,a heat conduction equation for temperature.The electrostatic potential appears within the concentration equations and heat conduction equation,and the electric field strength controls the concentrations and the temperature.The electric field potential is solved by the conservative block-centered method,and the order of the accuracy is improved by the electric potential.The concentrations and temperature are computed by the upwind blockcentered multistep method,where three different numerical methods are involved.The multistep method is adopted to approximate the time derivative.The blockcentered method is used to discretize the diffusion.The upwind scheme is applied to approximate the convection to avoid numerical dispersion and nonphysical oscillation.The block-centered difference simulates diffusion,concentrations,temperature,and the adjoint vector functions simultaneously.It has the local conservation of mass,which is an important nature in numerical simulation of a semiconductor device.By using the variation,energy estimates,induction hypothesis,embedding theorem and the technique of a priori estimates of differential equations,convergence of the optimal order is obtained.Numerical examples are provided to show the effectiveness and viability.This method provides a powerful tool for solving the challenging benchmark problem.展开更多
The traditional differential quadrature method was improved by using theupwind difference scheme for the convective terms to solve the coupled two-dimensionalincompressible Navier-stokes equations and heat equation. T...The traditional differential quadrature method was improved by using theupwind difference scheme for the convective terms to solve the coupled two-dimensionalincompressible Navier-stokes equations and heat equation. The new method was compared with theconventional differential quadrature method in the aspects of convergence and accuracy. The resultsshow that the new method is more accurate, and has better convergence than the conventionaldifferential quadrature method for numerically computing the steady-state solution.展开更多
In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order converge...In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order convergent, uniformly in the small parameter e , to the solution of problem (1.1). Numerical results are finally provided.展开更多
A high resolution upwind compact streamfunction numerical algorithm for two-dimensional(2D)double-diffusive convection(DDC)is developed.The unsteady Navier-Stokes(N-S)equations in the streamfunction-velocity form and ...A high resolution upwind compact streamfunction numerical algorithm for two-dimensional(2D)double-diffusive convection(DDC)is developed.The unsteady Navier-Stokes(N-S)equations in the streamfunction-velocity form and the scalar temperature and concentration equations are used.An optimized third-order upwind compact(UCD3 opt)scheme with a low dispersion error for the first derivatives is utilized to approximate the third derivatives of the streamfunction in the advection terms of the N-S equations and the first derivatives in the advection terms of the scalar temperature and concentration equations.The remaining first derivatives of the streamfunction(velocity),temperature,and concentration variables used in the governing equations are discretized by the fourth-order compact Pade(SCD4)schemes.With the temperature and concentration variables and their approximate values of the first derivatives obtained by the SCD4 schemes,the explicit fourth-order compact schemes are suggested to approximate the second derivatives of temperature and concentration in the diffusion terms of the energy and concentration equations.The discretization of the temporal term is executed with the second-order Crank-Nicolson(C-N)scheme.To assess the spatial behavior capability of the established numerical algorithm and verify the developed computer code,the DDC flow is numerically solved.The obtained results agree well with the benchmark solutions and some accurate results available in the literature,verifying the accuracy,effectiveness,and robustness of the provided algorithm.Finally,a preliminary application of the proposed method to the DDC is carried out.展开更多
Rotatability is a desirable quality of fitting response surface experimental designs. The propertystates that the variance of the estimated response made from the Taylor’s series expansion areconstant on circles, sph...Rotatability is a desirable quality of fitting response surface experimental designs. The propertystates that the variance of the estimated response made from the Taylor’s series expansion areconstant on circles, spheres and hyper-spheres about the centre of the design. In this article,a measure of rotatability of modified second-order rotatable design is presented. The variancefunction of a second-order response design and an infinite class of supplementary difference setsis used in coming up with the design.展开更多
We extend the SCGS smoothing procedure (Symmetrical Collective Gauss-Seidel relaxation) proposed by S. P. Vanka[4], for multigrid solvers of the steady viscous incompressible Navier-Stokes equations, to corresponding ...We extend the SCGS smoothing procedure (Symmetrical Collective Gauss-Seidel relaxation) proposed by S. P. Vanka[4], for multigrid solvers of the steady viscous incompressible Navier-Stokes equations, to corresponding line-wise versions. The resulting relaxation schemes are integrated into the multigrid solver based on second-order upwind differencing presented in [5]. Numerical comparisons on the efficiency of point-wise and line-wise relaxations are presented展开更多
The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differenti...The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.展开更多
基金Sponsored by the Aerospace Technology Innovation Funding(Grant No. CASC0209)
文摘A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative position,velocity and attitude of two unmanned aerial vehicles (UAVs).The second-order divided difference filter which makes use of multidimensional interpolation formulations to approximate the nonlinear transformations could achieve more accurate estimation and faster convergence from inaccurate initial conditions than standard extended Kalman filter.The filter formulation is based on relative motion equations.The global attitude parameterization is given by quarternion,while a generalized three-dimensional attitude representation is used to define the local attitude error.Simulation results are shown to compare the performance of the second-order divided difference filter with a standard extended Kalman filter approach.
文摘In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
文摘In this paper, combining the idea of difference method and finite element method, we construct a difference scheme for a self-adjoint problem in conservation form. Its solution uniformly converges to that of the original differential equation problem with order h3.
基金supported by National Natural Science Foundation of China(11101244,11271231)National Tackling Key Problems Program(20050200069)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘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.
文摘In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.
基金the Major State Basic Research Program of China(19990328)NNSF of China(19871051,19972039) the Doctorate Foundation of the State Education Commission
文摘For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of 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 2 norm are derived to determine the error, in the approximate solution.
基金supported by the Major State BasicResearch Program of China(19990328)the National Tackling Key Problem Programs(20050200069)+4 种基金the National Natural Science Foundation of China(1077112410372052)the Doctorate Foundation of the Ministryof Education of China(20030422047)Shandong Provance Natural Science Foundation(2R2009AQ12)the Independent Innovation Foundation of Shandong University(2010TS031)
文摘Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference 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 change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.
基金supported by the Major State Basic Research Development Program of China(No.G19990328)the National Key Technologies R&D Program of China (No.20050200069)+1 种基金the National Natural Science Foundation of China (Nos.10771124 and 10372052)the Ph. D. Pro-grams Foundation of Ministry of Education of China (No.20030422047)
文摘The research of the miscible oil and water displacement problem with moving boundary values is of great value to the history of oil-gas transport and accumulation in the basin evolution as well as to the rational evaluation in prospecting and exploiting oil-gas resources. The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values. For the twodimensional bounded region, the upwind finite difference schemes are proposed. Some techniques, such as the calculus of variations, the change of variables, and the theory of a priori estimates, are used. The optimal orderl2-norm estimates are derived for the errors in the approximate solutions. The research is important both theoretically and practically for the model analysis in the field, the model numerical method, and the software development.
文摘In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the discretization of p.d.e. via the preconditioned conjugate gradient method. For the second-order partial differential equations with Dirichlel boundary conditions, we prove that the condition number of the preconditioned system is O(1) while the condition number of the original system is O(m 2) Here m is the number of interior gridpoints in each direction. Such condition number produces a linear convergence rale.
基金supported by the Major State Basic Research Development Program of China(G19990328)National Tackling Key Program(2011ZX05011-004+6 种基金2011ZX0505220050200069)National Natural Science Foundation of China(11101244112712311077112410372052)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘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 ate quite interesting and satisfactory.
基金The project supported by Laboratory of Computational Physics,Institute of Applied Physics & Computational Mathematics,T.O.Box 80 0 9,Beijing 1 0 0 0 88
文摘In this paper, we first provide a generalized difference method for the 2-dimensional Navier-Stokes equations by combing the ideas of staggered scheme m and generalized upwind scheme in space, and by backward Euler time-stepping. Then we apply the abstract framework of to prove its long-time convergence. Finally, a numerical example for solving driven cavity flows is given.
基金Supported by the National Natural Science Foundation of China (Nos.50876114 and 10602043)the Program for New Century Excellent Talents in University,and the Scientific Research Key Project Fund of Ministry of Education (No.106142)
文摘The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics （English Edition）, 2007, 28（7）, 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.
文摘In this short article, the upwind and central compact finite difference schemes for spatial discretization of the first-order derivative are analyzed. Comparison of the schemes is provided and the best discretization scheme for convection dominated problems is suggested.
基金supported by National Natural Science Foundation of China(Grant No.11871312)Natural Science Foundation of Shandong Province(Grant No.ZR2021MA019).
文摘Numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science.The mathematical model is defined by a nonlinear system of initial-boundary problem including four partial differential equations:an elliptic equation for electrostatic potential,two convection-diffusion equations for electron concentration and hole concentration,a heat conduction equation for temperature.The electrostatic potential appears within the concentration equations and heat conduction equation,and the electric field strength controls the concentrations and the temperature.The electric field potential is solved by the conservative block-centered method,and the order of the accuracy is improved by the electric potential.The concentrations and temperature are computed by the upwind blockcentered multistep method,where three different numerical methods are involved.The multistep method is adopted to approximate the time derivative.The blockcentered method is used to discretize the diffusion.The upwind scheme is applied to approximate the convection to avoid numerical dispersion and nonphysical oscillation.The block-centered difference simulates diffusion,concentrations,temperature,and the adjoint vector functions simultaneously.It has the local conservation of mass,which is an important nature in numerical simulation of a semiconductor device.By using the variation,energy estimates,induction hypothesis,embedding theorem and the technique of a priori estimates of differential equations,convergence of the optimal order is obtained.Numerical examples are provided to show the effectiveness and viability.This method provides a powerful tool for solving the challenging benchmark problem.
文摘The traditional differential quadrature method was improved by using theupwind difference scheme for the convective terms to solve the coupled two-dimensionalincompressible Navier-stokes equations and heat equation. The new method was compared with theconventional differential quadrature method in the aspects of convergence and accuracy. The resultsshow that the new method is more accurate, and has better convergence than the conventionaldifferential quadrature method for numerically computing the steady-state solution.
文摘In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order convergent, uniformly in the small parameter e , to the solution of problem (1.1). Numerical results are finally provided.
基金supported by the National Natural Science Foundation of China(Nos.11872151,11372075,and 91330112)。
文摘A high resolution upwind compact streamfunction numerical algorithm for two-dimensional(2D)double-diffusive convection(DDC)is developed.The unsteady Navier-Stokes(N-S)equations in the streamfunction-velocity form and the scalar temperature and concentration equations are used.An optimized third-order upwind compact(UCD3 opt)scheme with a low dispersion error for the first derivatives is utilized to approximate the third derivatives of the streamfunction in the advection terms of the N-S equations and the first derivatives in the advection terms of the scalar temperature and concentration equations.The remaining first derivatives of the streamfunction(velocity),temperature,and concentration variables used in the governing equations are discretized by the fourth-order compact Pade(SCD4)schemes.With the temperature and concentration variables and their approximate values of the first derivatives obtained by the SCD4 schemes,the explicit fourth-order compact schemes are suggested to approximate the second derivatives of temperature and concentration in the diffusion terms of the energy and concentration equations.The discretization of the temporal term is executed with the second-order Crank-Nicolson(C-N)scheme.To assess the spatial behavior capability of the established numerical algorithm and verify the developed computer code,the DDC flow is numerically solved.The obtained results agree well with the benchmark solutions and some accurate results available in the literature,verifying the accuracy,effectiveness,and robustness of the provided algorithm.Finally,a preliminary application of the proposed method to the DDC is carried out.
文摘Rotatability is a desirable quality of fitting response surface experimental designs. The propertystates that the variance of the estimated response made from the Taylor’s series expansion areconstant on circles, spheres and hyper-spheres about the centre of the design. In this article,a measure of rotatability of modified second-order rotatable design is presented. The variancefunction of a second-order response design and an infinite class of supplementary difference setsis used in coming up with the design.
文摘We extend the SCGS smoothing procedure (Symmetrical Collective Gauss-Seidel relaxation) proposed by S. P. Vanka[4], for multigrid solvers of the steady viscous incompressible Navier-Stokes equations, to corresponding line-wise versions. The resulting relaxation schemes are integrated into the multigrid solver based on second-order upwind differencing presented in [5]. Numerical comparisons on the efficiency of point-wise and line-wise relaxations are presented
基金supported the Natural Science Foundation of Shandong Province(ZR2016AM08)Natural Science Foundation of Hunan Province(2018JJ2028)National Natural Science Foundation of China(11871312).
文摘The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.
