可移动半无限垂直圆柱体温度交变振荡时对MHD自由对流的影响

不可压缩粘性流体,流过可移动、温度交变振荡的半无限垂直圆柱体时,对MHD自由对流影响的数值解进行了研究.应用Crank-Nicolson型的隐式有限差分方法,求解无量纲、不稳定、非线性、耦合的偏微分控制方程.在不同参数下研究速度、温度和浓度分布的变化,还分析了局部及平均的表面摩擦力、Nusselt数和Sherwood数,并以图形形式给出.所得结果与其他文献的结果比较,有着很好的一致性.     

FINITE ELEMENT AND DISCONTINUOUS GALERKIN METHOD FOR STOCHASTIC HELMHOLTZ EQUATION IN TWO- AND THREE-DIMENSIONS

In this paper, we consider the finite element method and discontinuous Galerkin method for the stochastic Helmholtz equation in R^d （d = 2, 3）. Convergence analysis and error estimates are presented for the numerical solutions. The effects of the noises on the accuracy of the approximations are illustrated. Numerical experiments are carried out to verify our theoretical results.     

REDUCED BASIS METHOD FOR PARAMETRIZED ELLIPTIC ADVECTION-REACTION PROBLEMS

In this work we consider the Reduced Basis method for the solution of parametrized advection-reaction partial differential equations. For the generation of the basis we adopt a stabilized finite element method and we define the Reduced Basis method in the ＂primal- dual＂ formulation for this stabilized problem. We provide a priori Reduced Basis error estimates and we discuss the effects of the finite element approximation on the Reduced Basis error. We propose an adaptive algorithm, based on the a posteriori Reduced Basis error estimate, for the selection of the sample sets upon which the basis are built; the idea leading this algorithm is the minimization of the computational costs associated with the solution of the Reduced Basis problem. Numerical tests demonstrate the efficiency, in terms of computational costs, of the ＂primal-dual＂ Reduced Basis approach with respect to an ＂only primal＂ one. Parametrized advection-reaction partial differential equations, Reduced Basis method, ＂primal-dual＂ reduced basis approach, Stabilized finite element method, a posteriori error estimation.     

Two-time level ADI finite volume method for a class of second-order hyperbolic problems

In this paper, we develop a two-time level alternating direction implicit (ADI) method for a class of second-order hyperbolic problems on a rectangular domain. The method builds on the finite volume method with biquadratic basis functions for the discretization in space, and a Crank-Nicolson approach for the time stepping. We obtain a second-order error estimation in the H¹ norm. Numerical experiments are performed to demonstrate the theoretical findings.     

Hopscotch method: The numerical solution of the Frank-Kamenetskii partial differential equation

Numerical solutions to the Frank-Kamenetskii partial differential equation modelling a thermal explosion in a cylindrical vessel are obtained using the hopscotch scheme. We observe that a nonlinear source term in the equation leads to numerical difficulty and hence adjust the scheme to accommodate such a term. Numerical solutions obtained via MATLAB, MATHEMATICA and the Crank-Nicolson implicit scheme are employed as a means of comparison. To gain insight into the accuracy of the hopscotch scheme the solution is compared to a power series solution obtained via the Lie group method. The numerical solution is also observed to converge to a well-known steady state solution. A linear stability analysis is performed to validate the stability of the results obtained.     

High order difference schemes for the system of two space second order nonlinear hyperbolic equations with variable coefficients

In this paper, we develop implicit difference schemes of O(k⁴ + k²h² + h⁴), where k > 0, h > 0 are grid sizes in time and space coordinates, respectively, for solving the system of two space dimensional second order nonlinear hyperbolic partial differential equations with variable coefficients having mixed derivatives subject to appropriate initial and boundary conditions. The proposed difference method for the scalar equation is applied for the solution of wave equation in polar coordinates to obtain three level conditionally stable ADI method of O(k⁴ + k²h² + h⁴). Some physical nonlinear problems are provided to demonstrate the accuracy of the implementation.     

ALTERNATING BAND CRANK-NICOLSON METHOD FOR δu／δt=δ^2u／δx^2＋δ^2u／δy^2

the Alternating Segment Crank-Nicolson scheme for one-dimensional diffusion equation has been developed in [1],and the Alternating Block Crank-Nicolson method for two-dimensional problem in [2].The methods have the advantages of parallel computing,stability and good accuracy.In this paper for the two-dimensional diffusion equation,the net region is divided into bands,a special kind of block.This method is called the alternating Band Crank-Nicolson method.     

ALTERNATING BAND CRANK-NICOLSON METHOD FOR

The Alternating Segment Crank-Nicolson scheme for one-dimensional diffusion equation has been developed in [ 1 ], and the Alternating Block Crank-Nicolson method for two-dimensional problem in [2]. The methods have the advantages of parallel computing, stability and good accuracy. Tn this paper for the two-dimensional diffusion equation, the net region is divided into bands, a special kind of block. This method is called the alternating Band Crank-Nicolson method.     

An Improved Hybrid Method for Inverse Obstacle Scattering Problems

An improved hybrid method is introduced in this paper as a numerical method to reconstruct the scatterer by far-field pattern for just one incident direction with unknown physical properties of the scatterer. The improved hybrid method inherits the idea of the hybrid method by Kress and Serranho which is a combination of Newton and decomposition method, and it improves the hybrid method by introducing a general boundary condition. The numerical experiments show the feasibility of this method.     

A NUMERICAL STUDY OF UNIFORM SUPERCONVERGENCE OF LDG METHOD FOR SOLVING SINGULARLY PERTURBED PROBLEMS

In this paper, we consider the local discontinuous Galerkin method （LDG） for solving singularly perturbed convection-diffusion problems in one- and two-dimensional settings. The existence and uniqueness of the LDG solutions are verified. Numerical experiments demonstrate that it seems impossible to obtain uniform superconvergence for numerical fluxes under uniform meshes. Thanks to the implementation of two-type different anisotropic meshes, i.e., the Shishkin and an improved grade meshes, the uniform 2p ＋ i-order superconvergence is observed numerically for both one-dimensional and twodimensional cases.     

HOMOTOPY CONTINUATION METHODS FOR STOCHASTIC TWO-POINT BOUNDARY VALUE PROBLEMS DRIVEN BY ADDITIVE NOISES

In this paper the homotopy continuation method for stochastic two-point boundary value problems driven by additive noises is studied. The existence of the solution of the homotopy equation is proved. Numerical schemes are constructed and error estimates are obtained. Numerical experiments demonstrate the effectiveness of the homotopy continu- ation method over other commonly used methods such as the shooting method.     

High-order ADI orthogonal spline collocation method for a new 2D fractional integro-differential problem

We use the generalized L1 approximation for the Caputo fractional derivative, the second-order fractional quadrature rule approximation for the integral term, and a classical Crank-Nicolson alternating direction implicit (ADI) scheme for the time discretization of a new two-dimensional (2D) fractional integro-differential equation, in combination with a space discretization by an arbitrary-order orthogonal spline collocation (OSC) method. The stability of a Crank-Nicolson ADI OSC scheme is rigourously established, and error estimate is also derived. Finally, some numerical tests are given.     

A parameterized Newton method and a quasi-Newton method for nonsmooth equations

This paper presents a parameterized Newton method using generalized Jacobians and a Broyden-like method for solving nonsmooth equations. The former ensures that the method is well-defined even when the generalized Jacobian is singular. The latter is constructed by using an approximation function which can be formed for nonsmooth equations arising from partial differential equations and nonlinear complementarity problems. The approximation function method generalizes the splitting function method for nonsmooth equations. Locally superlinear convergence results are proved for the two methods. Numerical examples are given to compare the two methods with some other methods.This work is supported by the Australian Research Council.     

An exponential high‐order compact ADI method for 3D unsteady convection–diffusion problems

In this article, we develop an exponential high order compact alternating direction implicit (EHOC ADI) method for solving three dimensional (3D) unsteady convection–diffusion equations. The method, which requires only a regular seven‐point 3D stencil similar to that in the standard second‐order methods, is second order accurate in time and fourth‐order accurate in space and unconditionally stable. The resulting EHOC ADI scheme in each alternating direction implicit (ADI) solution step corresponding to a strictly diagonally dominant matrix equation can be solved by the application of the one‐dimensional tridiagonal Thomas algorithm with a considerable saving in computing time. Numerical experiments for three test problems are carried out to demonstrate the performance of the present method and to compare it with the classical Douglas–Gunn ADI method and the Karaa's high‐order compact ADI method. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

The Index of an Infinite Dimensional Implicit System

The idea of the index of a differential algebraic equation (DAE) (or implicit differential equation) has played a fundamental role in both the analysis of DAEs and the development of numerical algorithms for DAEs. DAEs frequently arise as partial discretizations of partial differential equations (PDEs). In order to relate properties of the PDE to those of the resulting DAE it is necessary to have a concept of the index of a possibly constrained PDE. Using the finite dimensional theory as motivation, this paper will examine what one appropriate analogue is for infinite dimensional systems. A general definition approach will be given motivated by the desire to consider numerical methods. Specific examples illustrating several kinds of behavior will be considered in some detail. It is seen that our definition differs from purely algebraic definitions. Numerical solutions, and simulation difficulties, can be misinterpreted if this index information is missing.     

WEAK APPROXIMATION OF OBLIQUELY REFLECTED DIFFUSIONS IN TIME-DEPENDENT DOMAINS

In an earlier paper, we proved the existence of solutions to the Skorohod problem with oblique reflection in time-dependent domains and, subsequently, applied this result to the problem of constructing solutions, in time-dependent domains, to stochastic differential equations with oblique reflection. In this paper we use these results to construct weak approximations of solutions to stochastic differential equations with oblique reflection, in time-dependent domains in Rd, by means of a projected Euler scheme. We prove that the constructed method has, as is the case for normal reflection and time-independent domains, an order of convergence equal to 1/2 and we evaluate the method empirically by means of two numerical examples. Furthermore, using a well-known extension of the Feynman-Kac formula, to stochastic differential equations with reflection, our method gives, in addition, a Monte Carlo method for solving second order parabolic partial differential equations with Robin boundary conditions in time-dependent domains.     

A high order finite difference method with Richardsonextrapolation for 3D convection diffusion equation

In this paper, we extend the Sun and Zhang’s [24] work on high order finite difference method, which is based on the Richardson extrapolation technique and an operator interpolation scheme for the one and two dimensional steady convection diffusion equations to the three dimensional case. Firstly, we employ a fourth order compact difference scheme to get the fourth order accurate solution on the fine and the coarse grids. Then, we use the Richardson extrapolation technique by combining the two approximate solutions to get a sixth order accurate solution on coarse grid. Finally, we apply an operator interpolation scheme to achieve the sixth order accurate solution on the fine grid. During this process, we use alternating direction implicit (ADI) method to solve the resulting linear systems. Numerical experiments are conducted to verify the accuracy and effectiveness of the present method.     

CONSTRUCTION OF GEOMETRIC PARTIAL DIFFERENTIAL EQUATIONS FOR LEVEL SETS

Geometric partial differential equations of level-set form are usually constructed by a variational method using either Dirac delta function or co-area formula in the energy functional to be minimized. However, the equations derived by these two approaches are not consistent. In this paper, we present a third approach for constructing the level-set form equations. By representing various differential geometry quantities and differential geometry operators in terms of the implicit surface, we are able to reformulate three classes of parametric geometric partial differential equations （second-order, fourth-order and sixth- order） into the level-set forms. The reformulation of the equations is generic and simple, and the resulting equations are consistent with their parametric form counterparts. We further prove that the equations derived using co-area formula are also consistent with the parametric forms. However, these equations are of much complicated forms than these given by the equations we derived.     

STATIC REGIME IMAGING OF LOCATIONS OF CERTAIN 3D ELECTROMAGNETIC IMPERFECTIONS FROM A BOUNDARY PERTURBATION FORMULA

We are concerned, in a static regime, with a three-dimensional bounded domain of certain an imaging approach of the locations in electromagnetic imperfections. This approach is related to Electrical Impedance Tomography and makes use of a new perturbation formula in the electric fields. We present two localization procedures, from a Current Pro- jection method that deals with the single imperfection context and an inverse Fourier process that is devoted to multiple imperfections configurations. These procedures extend those that were described in our previous work, since operating for a broader class of settings. Namely, the localization is additionally performed for certain purely electric imperfections, as established from numerical simulations.     

An upwind finite volume scheme and its maximum‐principle‐preserving ADI splitting for unsteady‐state advection‐diffusion equations

We develop an upwind finite volume (UFV) scheme for unsteady‐state advection‐diffusion partial differential equations (PDEs) in multiple space dimensions. We apply an alternating direction implicit (ADI) splitting technique to accelerate the solution process of the numerical scheme. We investigate and analyze the reason why the conventional ADI splitting does not satisfy maximum principle in the context of advection‐diffusion PDEs. Based on the analysis, we propose a new ADI splitting of the upwind finite volume scheme, the alternating‐direction implicit, upwind finite volume (ADFV) scheme. We prove that both UFV and ADFV schemes satisfy maximum principle and are unconditionally stable. We also derive their error estimates. Numerical results are presented to observe the performance of these schemes. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 211–226, 2003