A mixed collocation–finite difference method for 3D microscopic heat transport problems

Three-dimensional time-dependent initial-boundary value problems of a novel microscopic heat equation are solved by the mixed collocation–finite difference method in and on the boundaries of a particle when the thickness is much smaller than both the length and width. The collocation method on fixed grid size is used to approximate the space operator, whereas the finite difference scheme is used for time discretization. This new mixed method is applied to a novel heat problem in a particle, in order to compute the temperature distribution in and on the particle's surface. The second derivatives of the basis functions for the spectral approximation are derived. Direct substitution of derivatives in the model transforms the differential equation into a linear system of equations that is solved by the specific preconditioned conjugate gradient method. The high-order accuracy and resolution achieved by the proposed method allows one to obtain engineering-accuracy solution on coarse meshes. The consistency, stability and convergence analysis are provided and numerical results are presented.     

Multirate infinitesimal step methods for atmospheric flow simulation

The numerical solution of the Euler equations requires the treatment of processes in different temporal scales. Sound waves propagate fast compared to advective processes. Based on a spatial discretisation on staggered grids, a multirate time integration procedure is presented here generalising split-explicit Runge-Kutta methods. The advective terms are integrated by a Runge-Kutta method with a macro stepsize restricted by the CFL number. Sound wave terms are treated by small time steps respecting the CFL restriction dictated by the speed of sound.Split-explicit Runge-Kutta methods are generalised by the inclusion of fixed tendencies of previous stages. The stability barrier for the acoustics equation is relaxed by a factor of two.Asymptotic order conditions for the low Mach case are given. The relation to commutator-free exponential integrators is discussed. Stability is analysed for the linear acoustic equation. Numerical tests are executed for the linear acoustics and the nonlinear Euler equations.     

An Implicit Solver for the Time-Dependent Kohn-Sham Equation

The implicit numerical methods have the advantages on preserving the physical properties of the quantum system when solving the time-dependent Kohn-Sham equation. However, the efficiency issue prevents the practical applications of those implicit methods. In this paper, an implicit solver based on a class of Runge-Kutta methods and the finite element method is proposed for the time-dependent Kohn-Sham equation. The efficiency issue is partially resolved by three approaches, i.e., an $h$-adaptive mesh method is proposed to effectively restrain the size of the discretized problem, a complex-valued algebraic multigrid solver is developed for efficiently solving the derived linear system from the implicit discretization, as well as the OpenMP based parallelization of the algorithm. The numerical convergence, the ability on preserving the physical properties, and the efficiency of the proposed numerical method are demonstrated by a number of numerical experiments.     

A mixed and discontinuous Galerkin finite volume element method for incompressible miscible displacement problems in porous media

The incompressible miscible displacement problem in porous media is modeled by a coupled system of two nonlinear partial differential equations, the pressure‐velocity equation and the concentration equation. In this article, we present a mixed finite volume element method for the approximation of pressure‐velocity equation and a discontinuous Galerkin finite volume element method for the concentration equation. A priori error estimates in L^∞(L²) are derived for velocity, pressure, and concentration. Numerical results are presented to substantiate the validity of the theoretical results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

On a Fourier Method of Embedding Domains Using an Optimal Distributed Control

We propose a domain embedding method to solve second order elliptic problems in arbitrary two-dimensional domains. This method can be easily extended to three-dimensional problems. The method is based on formulating the problem as an optimal distributed control problem inside a rectangle in which the arbitrary domain is embedded. A periodic solution of the equation under consideration is constructed easily by making use of Fourier series. Numerical results obtained for Dirichlet problems are presented. The numerical tests show a high accuracy of the proposed algorithm and the computed solutions are in very good agreement with the exact solutions.     

IMPLICIT-EXPLICIT MULTISTEP FINITE ELEMENT-MIXED FINITE ELEMENT METHODS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

The transient behavior of a semiconductor device consists of a Poisson equa-tion for the electric potential and of two nonlinear parabolic equations for the electrondensity and hole density.The electric potential equation is discretized by a mixed finiteelement method. The electron and hole density equations are treated by implicit-explicitmultistep finite element methods. The schemes are very efficient. The optimal order errorestimates both in time and space are derived.     

一种新型压缩预处理CGS算法

CGS算法是求解大型非对称线性方程组的常用算法，然而该算法无极小残差性质，因此它常因出现较大的中间剩余向量而出现典型的不规则收敛行为．本根据IRA方法提出了一种压缩预处理CGS方法，数值实验表明这种算法在一定程度上减小了迭代算法在收敛过程中的剩余问题，从而使得算法具有更好的稳定性，该法构造简单，减少了收敛次数，加快了收敛速度．     

An Efficient and Accurate Spectral Method for Acoustic Scattering in Elliptic Domains

An efficient and accurate method for solving the two-dimensional Helmhokz equation in domains exterior to elongated obstacles is developed in this paper. The method is based on the so called transformed field expansion (TFE) coupled with a spectral-Galerkin solver for elliptical domain using Mathieu functions. Numerical results are presented to show the accuracy and stability of the proposed method.     

A new collocation method for solution of mixed linear integro-differential-difference equations

Numerical solution of mixed linear integro-differential-difference equation is presented using Chebyshev collocation method. The aim of this article is to present an efficient numerical procedure for solving mixed linear integro-differential-difference equations. Our method depends mainly on a Chebyshev expansion approach. This method transforms mixed linear integro-differential-difference equations and the given conditions into matrix equation which corresponds to a system of linear algebraic equation. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments and performed on the computer algebraic system Maple10.     

A weak Galerkin method for second order elliptic problems with polynomial reduction

The second order elliptic equation, which is also know as the diffusion-convection equation, is of great interest in many branches of physics and industry. In this paper, we use the weak Galerkin finite element method to study the general second order elliptic equation. A weak Galerkin finite element method is proposed and analyzed. This scheme features piecewise polynomials of degree $k\geq 1$ on each element and piecewise polynomials of degree $k-1\geq 0$ on each edge or face of the element. Error estimates of optimal order of convergence rate are established in both discrete $H^1$ and standard $L^2$ norm. The paper also presents some numerical experiments to verify the efficiency of the method.     

解Burgers方程的一种显式稳定性方法

拟线性Burgers方程在空间离散后转化成常微分方程,再用指数积分方法求解．数值结果表明指数积分法有显式稳定性,有相应Runge-Kutta方法相同的精度．     

Moving collocation methods for time fractional differential equations and simulation of blowup

A moving collocation method is proposed and implemented to solve time fractional differential equations. The method is derived by writing the fractional differential equation into a form of time difference equation. The method is stable and has a third-order convergence in space and first-order convergence in time for either linear or nonlinear equations. In addition, the method is used to simulate the blowup in the nonlinear equations.     

Lagrange-SQP Techniques for the Control Constrained Optimal Boundary Control for the Burgers Equation

A Lagrange-Newton-SQP method is analyzed for the optimal control of the Burgers equation. Boundary controls are given, which are restricted by pointwise lower and upper bounds. The convergence of the method is proved in appropriate Banach spaces. This proof is based on a second-order sufficient optimality condition and the theory of Newton methods for generalized equations in Banach spaces. For the numerical realization a primal-dual active set strategy is applied. To illustrate the theoretical investigations, numerical examples are included. Moreover, a globalization technique for the SQP method is tested numerically.     

Piecewise parabolic method on local stencil for gasdynamic simulations

A numerical method based on piecewise parabolic difference approximations is proposed for solving hyperbolic systems of equations. The design of its numerical scheme is based on the conservation of Riemann invariants along the characteristic curves of a system of equations, which makes it possible to discard the four-point interpolation procedure used in the standard piecewise parabolic method (PPM) and to use the data from the previous time level in the reconstruction of the solution inside difference cells. As a result, discontinuous solutions can be accurately represented without adding excessive dissipation. A local stencil is also convenient for computations on adaptive meshes. The new method is compared with PPM by solving test problems for the linear advection equation and the inviscid Burgers equation. The efficiency of the methods is compared in terms of errors in various norms. A technique for solving the gas dynamics equations is described and tested for several one-and two-dimensional problems.     

Exponential splitting time integration for pseudospectral methods on moving meshes

Moving meshes are successfully used in many fields. Here we investigate how a recently proposed approach to combine the Strang splitting method for time integration with pseudospectral spatial discretization by orthogonal polynomials can be extended to include moving meshes. A double representation of a function (by coefficients of polynomial expansion and by values at the mesh nodes associated with a suitable quadrature formula) is an essential part of the numerical integration. Before numerical implementation the original PDE is transformed into a suitable form. The approach is illustrated on the linear heat transfer equation.     

Numerical Treatment of Homogeneous Semi-Markov Processes in Transient Case–a Straightforward Approach

This paper presents the numerical solution of the process evolution equation of a homogeneous semi-Markov process (HSMP) with a general quadrature method. Furthermore, results that justify this approach proving that the numerical solution tends to the evolution equation of the continuous time HSMP are given. The results obtained generalize classical results on integral equation numerical solutions applying them to particular kinds of integral equation systems. A method for obtaining the discrete time HSMP is shown by applying a very particular quadrature formula for the discretization. Following that, the problem of obtaining the continuous time HSMP from the discrete one is considered. In addition, the discrete time HSMP in matrix form is presented and the fact that the solution of the evolution equation of this process always exists is proved. Afterwards, an algorithm for solving the discrete time HSMP is given. Finally, a simple application of the HSMP is given for a real data social security example.     

On the Numerical Solution of Some Eikonal Equations: An Elliptic Solver Approach

The steady Eikonal equation is a prototypical first-order fully nonlinear equation. A numerical method based on elliptic solvers is presented here to solve two different kinds of steady Eikonal equations and compute solutions, which are maximal and minimal in the variational sense. The approach in this paper relies on a variational argument involving penalty, a biharmonic regularization, and an operator-splitting-based time-discretization scheme for the solution of an associated initial-value problem. This approach allows the decoupling of the nonlinearities and differential operators. Numerical experiments are performed to validate this approach and investigate its convergence properties from a numerical viewpoint.     

Exact finite-difference methods based on nonlinear means

The necessary and sufficient condition for an autonomous, ordinary differential equation to be exactly solved by a given Evans-Sanugi, nonlinear, one-step, finite difference method based on ‘classical means’ are derived. A necessary, but not sufficient, condition yields the most general nonlinearity for the differential equation which is independent of the step size. Examples of differential equations for which either nonlinear trapezoidal or nonlinear implicit midpoint methods based on arithmetic, harmonic, contraharmonic, quadratic, geometric, Heronian, centroidal and logarithmic means are exact, are presented. These new exact difference schemes may be useful in future developments of new Denk-Bulirsch, Le-Roux or Kojouhavov-Chen schemes for nonlinear evolution equations with or without blow-up.     

Asymptotic solutions to a discrete airy equation

A nonstandard finite difference scheme for the Airy equation leads to a linear, second—order difference equation for which the theorems of Poincaré and Perron do not apply if asymptotic representations of the solutions are desired. Using the method of dominant balance, a suggested form for the asymptotic solutions is obtained. This relation is then used to construct the required asymptotic representations for the two linearly independent solutions     

Multilevel augmentation methods for solving the sine-Gordon equation

In this paper we develop the multilevel augmentation method for solving nonlinear operator equations of the second kind and apply it to solving the one-dimensional sine-Gordon equation. We first give a general setting of the multilevel augmentation method for solving the second kind nonlinear operator equations and prove that the multilevel augmentation method preserves the optimal convergence order of the projection method while reducing computational cost significantly. Then we describe the semi-discrete scheme and the fully-discrete scheme based on multiscale methods for solving the sine-Gordon equation, and apply the multilevel augmentation method to solving the discrete equation. A complete analysis for convergence order is proposed. Finally numerical experiments are presented to confirm the theoretical results and illustrate the efficiency of the method.