for Journals by Title or ISSN
for Articles by Keywords
  Subjects -> COMPUTER SCIENCE (Total: 1991 journals)
    - ANIMATION AND SIMULATION (27 journals)
    - AUTOMATION (40 journals)
    - COMPUTER ARCHITECTURE (9 journals)
    - COMPUTER ENGINEERING (11 journals)
    - COMPUTER GAMES (16 journals)
    - COMPUTER GRAPHICS (32 journals)
    - COMPUTER PROGRAMMING (19 journals)
    - COMPUTER SCIENCE (1232 journals)
    - COMPUTER SECURITY (44 journals)
    - CYBERNETICS AND ROBOTICS (62 journals)
    - DATA BASE MANAGEMENT (15 journals)
    - DATA MINING (23 journals)
    - INFORMATION SYSTEMS (93 journals)
    - INTERNET (91 journals)
    - SOCIAL WEB (46 journals)
    - SOFTWARE (32 journals)
    - THEORY OF COMPUTING (6 journals)

COMPUTER SCIENCE (1232 journals)                  1 2 3 4 5 6 7 | Last

Showing 1 - 200 of 872 Journals sorted alphabetically
3D Printing and Additive Manufacturing     Full-text available via subscription   (Followers: 10)
Abakós     Open Access   (Followers: 1)
Academy of Information and Management Sciences Journal     Full-text available via subscription   (Followers: 48)
ACM Computing Surveys     Hybrid Journal   (Followers: 19)
ACM Journal on Computing and Cultural Heritage     Hybrid Journal   (Followers: 5)
ACM Journal on Emerging Technologies in Computing Systems     Hybrid Journal   (Followers: 8)
ACM Transactions on Accessible Computing (TACCESS)     Hybrid Journal   (Followers: 2)
ACM Transactions on Algorithms (TALG)     Hybrid Journal   (Followers: 11)
ACM Transactions on Applied Perception (TAP)     Hybrid Journal   (Followers: 4)
ACM Transactions on Architecture and Code Optimization (TACO)     Hybrid Journal   (Followers: 5)
ACM Transactions on Autonomous and Adaptive Systems (TAAS)     Hybrid Journal   (Followers: 5)
ACM Transactions on Computation Theory (TOCT)     Hybrid Journal   (Followers: 10)
ACM Transactions on Computational Logic (TOCL)     Hybrid Journal   (Followers: 2)
ACM Transactions on Computer Systems (TOCS)     Hybrid Journal   (Followers: 14)
ACM Transactions on Computer-Human Interaction     Hybrid Journal   (Followers: 9)
ACM Transactions on Computing Education (TOCE)     Hybrid Journal   (Followers: 2)
ACM Transactions on Design Automation of Electronic Systems (TODAES)     Hybrid Journal   (Followers: 1)
ACM Transactions on Economics and Computation     Hybrid Journal  
ACM Transactions on Embedded Computing Systems (TECS)     Hybrid Journal   (Followers: 3)
ACM Transactions on Information Systems (TOIS)     Hybrid Journal   (Followers: 19)
ACM Transactions on Intelligent Systems and Technology (TIST)     Hybrid Journal   (Followers: 7)
ACM Transactions on Interactive Intelligent Systems (TiiS)     Hybrid Journal   (Followers: 3)
ACM Transactions on Multimedia Computing, Communications, and Applications (TOMCCAP)     Hybrid Journal   (Followers: 6)
ACM Transactions on Reconfigurable Technology and Systems (TRETS)     Hybrid Journal   (Followers: 5)
ACM Transactions on Sensor Networks (TOSN)     Hybrid Journal   (Followers: 6)
ACM Transactions on Speech and Language Processing (TSLP)     Hybrid Journal   (Followers: 8)
ACM Transactions on Storage     Hybrid Journal  
ACS Applied Materials & Interfaces     Full-text available via subscription   (Followers: 16)
Acta Automatica Sinica     Full-text available via subscription   (Followers: 3)
Acta Informatica     Hybrid Journal   (Followers: 4)
Acta Universitatis Cibiniensis. Technical Series     Open Access  
Ad Hoc Networks     Hybrid Journal   (Followers: 10)
Adaptive Behavior     Hybrid Journal   (Followers: 9)
Advanced Engineering Materials     Hybrid Journal   (Followers: 18)
Advanced Science Letters     Full-text available via subscription   (Followers: 4)
Advances in Adaptive Data Analysis     Hybrid Journal   (Followers: 7)
Advances in Artificial Intelligence     Open Access   (Followers: 13)
Advances in Artificial Neural Systems     Open Access   (Followers: 4)
Advances in Calculus of Variations     Hybrid Journal   (Followers: 2)
Advances in Catalysis     Full-text available via subscription   (Followers: 4)
Advances in Computational Mathematics     Hybrid Journal   (Followers: 14)
Advances in Computer Engineering     Open Access   (Followers: 1)
Advances in Computer Science : an International Journal     Open Access   (Followers: 5)
Advances in Computing     Open Access   (Followers: 2)
Advances in Data Analysis and Classification     Hybrid Journal   (Followers: 43)
Advances in Engineering Software     Hybrid Journal   (Followers: 19)
Advances in Geosciences (ADGEO)     Open Access   (Followers: 7)
Advances in Human Factors/Ergonomics     Full-text available via subscription   (Followers: 19)
Advances in Human-Computer Interaction     Open Access   (Followers: 16)
Advances in Materials Sciences     Open Access   (Followers: 10)
Advances in Operations Research     Open Access   (Followers: 9)
Advances in Parallel Computing     Full-text available via subscription   (Followers: 7)
Advances in Porous Media     Full-text available via subscription   (Followers: 3)
Advances in Remote Sensing     Open Access   (Followers: 13)
Advances in Science and Research (ASR)     Open Access   (Followers: 6)
AEU - International Journal of Electronics and Communications     Hybrid Journal   (Followers: 5)
African Journal of Information and Communication     Open Access   (Followers: 7)
African Journal of Mathematics and Computer Science Research     Open Access   (Followers: 2)
AGE     Hybrid Journal   (Followers: 4)
Age and Ageing     Hybrid Journal   (Followers: 54)
AI EDAM     Hybrid Journal  
Air, Soil & Water Research     Open Access   (Followers: 6)
AIS Transactions on Human-Computer Interaction     Open Access   (Followers: 5)
Algebras and Representation Theory     Hybrid Journal  
Algorithmica     Hybrid Journal   (Followers: 5)
Algorithms     Open Access   (Followers: 9)
American Journal of Applied Sciences     Open Access   (Followers: 26)
American Journal of Computational and Applied Mathematics     Open Access   (Followers: 2)
American Journal of Computational Mathematics     Open Access   (Followers: 1)
American Journal of Information Systems     Open Access   (Followers: 4)
American Journal of Sensor Technology     Open Access  
Anais da Academia Brasileira de Ciências     Open Access   (Followers: 2)
Analog Integrated Circuits and Signal Processing     Hybrid Journal   (Followers: 4)
Analysis in Theory and Applications     Hybrid Journal  
Animation Practice, Process & Production     Hybrid Journal   (Followers: 4)
Annals of Combinatorics     Hybrid Journal   (Followers: 3)
Annals of Data Science     Hybrid Journal   (Followers: 4)
Annals of Mathematics and Artificial Intelligence     Hybrid Journal   (Followers: 5)
Annals of Pure and Applied Logic     Open Access   (Followers: 1)
Annals of Software Engineering     Hybrid Journal   (Followers: 8)
Annual Reviews in Control     Hybrid Journal   (Followers: 5)
Applicable Algebra in Engineering, Communication and Computing     Hybrid Journal   (Followers: 1)
Applied and Computational Harmonic Analysis     Full-text available via subscription   (Followers: 2)
Applied Artificial Intelligence: An International Journal     Hybrid Journal   (Followers: 10)
Applied Categorical Structures     Hybrid Journal   (Followers: 1)
Applied Computational Intelligence and Soft Computing     Open Access   (Followers: 9)
Applied Computer Systems     Open Access  
Applied Environmental Education & Communication     Hybrid Journal   (Followers: 13)
Applied Informatics     Open Access  
Applied Intelligence     Hybrid Journal   (Followers: 7)
Applied Mathematics and Computation     Hybrid Journal   (Followers: 30)
Applied Medical Informatics     Open Access   (Followers: 9)
Applied Numerical Analysis & Computational Mathematics     Hybrid Journal   (Followers: 5)
Applied Numerical Mathematics     Hybrid Journal   (Followers: 5)
Applied Soft Computing     Hybrid Journal   (Followers: 11)
Applied Spatial Analysis and Policy     Hybrid Journal   (Followers: 3)
Architectural Theory Review     Hybrid Journal   (Followers: 2)
Archive of Applied Mechanics     Hybrid Journal   (Followers: 4)
Archives and Museum Informatics     Hybrid Journal   (Followers: 97)
Archives of Computational Methods in Engineering     Hybrid Journal   (Followers: 4)
arq: Architectural Research Quarterly     Hybrid Journal   (Followers: 7)
Artifact     Hybrid Journal   (Followers: 2)
Artificial Intelligence and Law     Hybrid Journal   (Followers: 5)
Artificial Intelligence in Medicine     Hybrid Journal   (Followers: 10)
Artificial Intelligence Review     Hybrid Journal   (Followers: 11)
Artificial Life     Hybrid Journal   (Followers: 5)
Asia Pacific Journal on Computational Engineering     Open Access  
Asia-Pacific Journal of Information Technology and Multimedia     Open Access  
Asian Journal of Communication     Hybrid Journal   (Followers: 6)
Asian Journal of Control     Hybrid Journal  
Assembly Automation     Hybrid Journal   (Followers: 1)
at - Automatisierungstechnik     Hybrid Journal   (Followers: 2)
Atlantic Journal of Communication     Hybrid Journal   (Followers: 5)
Australian Educational Computing     Open Access  
Automated Software Engineering     Hybrid Journal   (Followers: 7)
Automatic Control and Computer Sciences     Hybrid Journal   (Followers: 3)
Automatic Documentation and Mathematical Linguistics     Hybrid Journal   (Followers: 4)
Automatica     Hybrid Journal   (Followers: 8)
Automation and Remote Control     Hybrid Journal   (Followers: 2)
Automation in Construction     Hybrid Journal   (Followers: 6)
Autonomous Agents and Multi-Agent Systems     Hybrid Journal   (Followers: 5)
Autonomous Mental Development, IEEE Transactions on     Hybrid Journal   (Followers: 5)
Basin Research     Hybrid Journal   (Followers: 3)
Behaviour & Information Technology     Hybrid Journal   (Followers: 47)
Bioinformatics     Hybrid Journal   (Followers: 165)
Biomedical Engineering     Hybrid Journal   (Followers: 14)
Biomedical Engineering and Computational Biology     Open Access   (Followers: 11)
Biomedical Engineering, IEEE Reviews in     Full-text available via subscription   (Followers: 15)
Biomedical Engineering, IEEE Transactions on     Hybrid Journal   (Followers: 23)
Briefings in Bioinformatics     Hybrid Journal   (Followers: 40)
British Journal of Educational Technology     Hybrid Journal   (Followers: 82)
Broadcasting, IEEE Transactions on     Hybrid Journal   (Followers: 7)
CALCOLO     Hybrid Journal  
Calphad     Hybrid Journal  
Canadian Journal of Electrical and Computer Engineering     Full-text available via subscription   (Followers: 7)
Capturing Intelligence     Full-text available via subscription  
Catalysis in Industry     Hybrid Journal  
CEAS Space Journal     Hybrid Journal  
Cell Communication and Signaling     Open Access   (Followers: 1)
Central European Journal of Computer Science     Hybrid Journal   (Followers: 4)
Cereal Research Communications     Full-text available via subscription   (Followers: 1)
Chaos, Solitons & Fractals     Hybrid Journal   (Followers: 1)
Chemical Communications     Full-text available via subscription   (Followers: 62)
Chemometrics and Intelligent Laboratory Systems     Hybrid Journal   (Followers: 14)
ChemSusChem     Hybrid Journal   (Followers: 6)
China Communications     Full-text available via subscription   (Followers: 5)
Chinese Journal of Catalysis     Full-text available via subscription   (Followers: 1)
CIN Computers Informatics Nursing     Full-text available via subscription   (Followers: 10)
Circuits and Systems     Open Access   (Followers: 9)
Clean Air Journal     Full-text available via subscription   (Followers: 2)
CLEI Electronic Journal     Open Access  
Clin-Alert     Hybrid Journal   (Followers: 1)
Cluster Computing     Hybrid Journal   (Followers: 1)
Cognition and Instruction     Hybrid Journal   (Followers: 9)
Cognitive Computation     Hybrid Journal   (Followers: 4)
Cognitive Processing     Hybrid Journal   (Followers: 4)
Cognitive Science     Hybrid Journal   (Followers: 25)
Cognitive Systems Research     Hybrid Journal   (Followers: 5)
COMBINATORICA     Hybrid Journal  
Combinatorics, Probability and Computing     Hybrid Journal   (Followers: 3)
Combustion Theory and Modelling     Hybrid Journal   (Followers: 10)
Communication Methods and Measures     Hybrid Journal   (Followers: 12)
Communication Quarterly     Hybrid Journal   (Followers: 11)
Communication Reports     Hybrid Journal   (Followers: 7)
Communication Research     Hybrid Journal   (Followers: 16)
Communication Research Reports     Hybrid Journal   (Followers: 6)
Communication Studies     Hybrid Journal   (Followers: 11)
Communication Teacher     Hybrid Journal   (Followers: 4)
Communication Theory     Hybrid Journal   (Followers: 16)
Communications Engineer     Hybrid Journal  
Communications in Algebra     Hybrid Journal   (Followers: 2)
Communications in Computational Physics     Full-text available via subscription   (Followers: 1)
Communications in Partial Differential Equations     Hybrid Journal   (Followers: 1)
Communications of the ACM     Full-text available via subscription   (Followers: 29)
Communications of the Association for Information Systems     Open Access   (Followers: 17)
Communications On Pure & Applied Mathematics     Hybrid Journal   (Followers: 3)
COMPEL: The International Journal for Computation and Mathematics in Electrical and Electronic Engineering     Hybrid Journal   (Followers: 3)
Complex Adaptive Systems Modeling     Open Access  
Complex Analysis and Operator Theory     Hybrid Journal   (Followers: 1)
Complexity     Hybrid Journal   (Followers: 6)
Composite Materials Series     Full-text available via subscription   (Followers: 7)
Computación y Sistemas     Open Access  
Computation     Open Access  
Computational and Applied Mathematics     Hybrid Journal   (Followers: 2)
Computational and Mathematical Methods in Medicine     Open Access   (Followers: 2)
Computational and Mathematical Organization Theory     Hybrid Journal   (Followers: 2)
Computational and Structural Biotechnology Journal     Open Access   (Followers: 1)
Computational and Theoretical Chemistry     Hybrid Journal   (Followers: 10)
Computational Astrophysics and Cosmology     Open Access  
Computational Biology and Chemistry     Hybrid Journal   (Followers: 10)
Computational Biology Journal     Open Access   (Followers: 3)
Computational Chemistry     Open Access   (Followers: 2)
Computational Cognitive Science     Open Access  
Computational Complexity     Hybrid Journal   (Followers: 4)
Computational Condensed Matter     Open Access  
Computational Ecology and Software     Open Access   (Followers: 6)
Computational Economics     Hybrid Journal   (Followers: 9)
Computational Geosciences     Hybrid Journal   (Followers: 12)
Computational Intelligence     Hybrid Journal   (Followers: 5)
Computational Intelligence and Neuroscience     Open Access   (Followers: 8)

        1 2 3 4 5 6 7 | Last

Journal Cover Applied Numerical Mathematics
  [SJR: 1.163]   [H-I: 49]   [5 followers]  Follow
   Hybrid Journal Hybrid journal (It can contain Open Access articles)
   ISSN (Print) 0168-9274 - ISSN (Online) 0168-9274
   Published by Elsevier Homepage  [2970 journals]
  • Editorial Board
    • Abstract: Publication date: August 2016
      Source:Applied Numerical Mathematics, Volume 106

      PubDate: 2016-05-19T15:01:16Z
  • A projection-based variational multiscale method for the optimal control
           problems governed by the stationary Navier–Stokes equations
    • Abstract: Publication date: August 2016
      Source:Applied Numerical Mathematics, Volume 106
      Author(s): Fikriye Yılmaz, Aytekin Çıbık
      In this work, we apply a variational multiscale stabilization (VMS) to the optimal control problems of Navier–Stokes equations. We first obtain the optimality conditions by using Lagrange approach. After stating the continuous optimality system, we formulate the discrete problem by a projection-based stabilized finite element scheme. We give the stability estimates for both state and adjoint state variables. Then, we present the a priori error analysis. The main issue of this paper is to show the efficiency of the variational multiscale stabilization for this optimal control problem. We verify our results by numerical examples.

      PubDate: 2016-05-19T15:01:16Z
  • Fast iterative solvers for large matrix systems arising from
           time-dependent Stokes control problems
    • Abstract: Publication date: Available online 16 May 2016
      Source:Applied Numerical Mathematics
      Author(s): John W. Pearson
      In this manuscript we consider the development of fast iterative solvers for Stokes control problems, an important class of PDE-constrained optimization problems. In particular we wish to develop effective preconditioners for the matrix systems arising from finite element discretizations of time-dependent variants of such problems. To do this we consider a suitable rearrangement of the matrix systems, and exploit the saddle point structure of many of the relevant sub-matrices involved – we may then use this to construct representations of these sub-matrices based on good approximations of their ( 1 , 1 ) -block and Schur complement. We test our recommended iterative methods on a distributed control problem with Dirichlet boundary conditions, and on a time-periodic problem.

      PubDate: 2016-05-19T15:01:16Z
  • High order exponentially fitted methods for Volterra integral equations
           with periodic solution
    • Abstract: Publication date: Available online 17 May 2016
      Source:Applied Numerical Mathematics
      Author(s): A. Cardone, R. D'Ambrosio, B. Paternoster
      The present paper illustrates the construction of direct quadrature methods of arbitrary order for Volterra integral equations with periodic solution. The coefficients of these methods depend on the parameters of the problem, following the exponential fitting theory. The convergence of these methods is analyzed, and some numerical experiments are illustrated to confirm theoretical expectations and for comparison with other existing methods.

      PubDate: 2016-05-19T15:01:16Z
  • A Newton type linearization based two grid method for coupling fluid flow
           with porous media flow
    • Abstract: Publication date: August 2016
      Source:Applied Numerical Mathematics, Volume 106
      Author(s): Peiqi Huang, Mingchao Cai, Feng Wang
      In this paper, we propose a two-grid finite element method for solving the mixed Navier–Stokes/Darcy model with the Beavers–Joseph–Saffman interface condition. After solving a coupled nonlinear problem on a coarse grid, we sequentially solve decoupled and linearized subproblems on a fine grid and then correct the solution on the same grid. Compared with the existing work on the two-grid methods for the coupled model, our two-grid method allows a much higher order scaling between the coarse grid size H and the fine grid size h. Specifically, if a k-th order discretization is applied, by using h = H 2 k + 1 k for k = 1 , 2 and h = H k + 1 k − 1 for k ≥ 3 , the final step two-grid solution errors in the energy norm are still optimal. Numerical experiments are also given to confirm the theoretical analysis.

      PubDate: 2016-05-19T15:01:16Z
  • Block extrapolation methods with applications
    • Abstract: Publication date: August 2016
      Source:Applied Numerical Mathematics, Volume 106
      Author(s): K. Jbilou, A. Messaoudi
      In the present paper we introduce new block extrapolation methods as generalizations of the well known vector extrapolation methods. We give expressions of the obtained approximations via the Schur complement and also propose an efficient implementation of these methods. Applications to linearly generated sequences are given and extensions to nonlinear problems are also given. Applications of the proposed block extrapolation methods to some nonlinear matrix equations are considered and some numerical examples are given.

      PubDate: 2016-05-19T15:01:16Z
  • Analysis of a posteriori error estimates of the discontinuous Galerkin
           method for nonlinear ordinary differential equations
    • Abstract: Publication date: August 2016
      Source:Applied Numerical Mathematics, Volume 106
      Author(s): Mahboub Baccouch
      We develop and analyze a new residual-based a posteriori error estimator for the discontinuous Galerkin (DG) method for nonlinear ordinary differential equations (ODEs). The a posteriori DG error estimator under investigation is computationally simple, efficient, and asymptotically exact. It is obtained by solving a local residual problem with no boundary condition on each element. We first prove that the DG solution exhibits an optimal O ( h p + 1 ) convergence rate in the L 2 -norm when p-degree piecewise polynomials with p ≥ 1 are used. We further prove that the DG solution is O ( h 2 p + 1 ) superconvergent at the downwind points. We use these results to prove that the p-degree DG solution is O ( h p + 2 ) super close to a particular projection of the exact solution. This superconvergence result allows us to show that the true error can be divided into a significant part and a less significant part. The significant part of the discretization error for the DG solution is proportional to the ( p + 1 ) -degree right Radau polynomial and the less significant part converges at O ( h p + 2 ) rate in the L 2 -norm. Numerical experiments demonstrate that the theoretical rates are optimal. Based on the global superconvergent approximations, we construct asymptotically exact a posteriori error estimates and prove that they converge to the true errors in the L 2 -norm under mesh refinement. The order of convergence is proved to be p + 2 . Finally, we prove that the global effectivity index in the L 2 -norm converges to unity at O ( h ) rate. Several numerical examples are provided to illustrate the global superconvergence results and the convergence of the proposed estimator under mesh refinement. A local adaptive procedure that makes use of our local a posteriori error estimate is also presented.

      PubDate: 2016-05-19T15:01:16Z
  • The Prothero and Robinson example: Convergence studies for
           Runge–Kutta and Rosenbrock–Wanner methods
    • Abstract: Publication date: Available online 11 May 2016
      Source:Applied Numerical Mathematics
      Author(s): Joachim Rang
      It is well-known that one-step methods have order reduction if they are applied on stiff ODEs such as the example of Prothero–Robinson. In this paper we analyse the local error of Runge–Kutta and Rosenbrock–Wanner methods. We derive new order conditions and define with them B P R -consistency. We show that for strongly A-stable methods B P R -consistency implies B P R -convergence. Finally we analyse methods from literature, derive new B P R -consistent methods and present numerical examples. The numerical and analytical results show the influence of different properties of the methods and of different order conditions on the numerical error and on the numerical convergence order.

      PubDate: 2016-05-13T14:46:49Z
  • Weak convergence for a stochastic exponential integrator and finite
           element discretization of stochastic partial differential equation with
           multiplicative & additive noise
    • Abstract: Publication date: Available online 9 May 2016
      Source:Applied Numerical Mathematics
      Author(s): Antoine Tambue, Jean Medard T. Ngnotchouye
      We consider a finite element approximation of a general semi-linear stochastic partial differential equation (SPDE) driven by space-time multiplicative and additive noise. We examine the full weak convergence rate of the exponential Euler scheme when the linear operator is self adjoint and also provide the full weak convergence rate for non-self-adjoint linear operator with additive noise. Key part of the proof does not rely on Malliavin calculus. For non-self-adjoint operators, we analyse the optimal strong error for spatially semi discrete approximations for both multiplicative and additive noise with truncated and non-truncated noise. Depending on the regularity of the noise and the initial solution, we found that in some cases the rate of weak convergence is twice the rate of the strong convergence. Our convergence rate is in agreement with some numerical results in two dimensions.

      PubDate: 2016-05-13T14:46:49Z
  • A Maximum-Principle-satisfying finite volume compact-WENO Scheme for
           traffic flow model on networks
    • Abstract: Publication date: Available online 10 May 2016
      Source:Applied Numerical Mathematics
      Author(s): YuFeng Shi, Yan Guo
      In this paper, we apply a maximum-principle-satisfying finite volume compact weighted scheme to numerical modelling traffic flow problems on networks. Road networks can be numerically model as a graph, whose edges are a finite number of roads that join at junctions. The evolution on each road is described by a scalar hyperbolic conservation law, and traffic distribution matrices are used to formulate coupling conditions at the network junctions. In order to achieve maximum-principle of the traffic density on each road, the maximum-principle-satisfying polynomial rescaling limiter is adopted. Numerical results for road networks with rich solution structures are presented in this work and indicate that the finite volume compact weighted scheme produces essentially non-oscillatory, maximum principle preserving and high resolution solutions.

      PubDate: 2016-05-13T14:46:49Z
  • Algorithms for the Haar wavelet based fast evaluation of aggregation
           integrals in population balance equations
    • Abstract: Publication date: Available online 10 May 2016
      Source:Applied Numerical Mathematics
      Author(s): Sabine Le Borne, Lusine Shahmuradyan
      In several production processes, the distribution of particles dispersed in an environmental phase may be mathematically described by the solution of population balance equations. We are concerned with the development of efficient numerical techniques for the aggregation process: It invokes an integral term that is usually numerically expensive to evaluate and often dominates the total simulation cost. We describe an approach on locally refined nested grids to evaluate both the source and the sink terms in almost linear complexity (instead of quadratic complexity resulting from a direct approach). The key is to switch from a nodal to a wavelet basis representation of the density function. We illustrate the numerical performance of this approach, both in comparison to a discretization of piecewise constant functions on a uniform grid as well as to the fixed pivot method on a geometric grid.

      PubDate: 2016-05-13T14:46:49Z
  • Analysis of stretched grids as buffer zones in simulations of wave
    • Abstract: Publication date: Available online 2 May 2016
      Source:Applied Numerical Mathematics
      Author(s): G. Kreiss, B. Krank, G. Efraimsson
      A zone of increasingly stretched grid is a robust and easy-to-use way to avoid unwanted reflections at artificial boundaries in wave propagating simulations. In such a buffer zone there are two main damping mechanisms, dissipation and under-resolution that turns a traveling wave into an evanescent wave. We present analysis in one and two space dimensions showing that evanescent decay through under-resolution is a very efficient way to damp waves. The analysis is supported by numerical computations.

      PubDate: 2016-05-07T14:35:44Z
  • Optical flow with fractional order regularization: Variational model and
           solution method
    • Abstract: Publication date: Available online 4 May 2016
      Source:Applied Numerical Mathematics
      Author(s): Somayeh Gh. Bardeji, Isabel N. Figueiredo, Ercília Sousa
      An optical flow variational model is proposed for a sequence of images defined on a domain in R 2 . We introduce a regularization term given by the L 1 norm of a fractional differential operator. To solve the minimization problem we apply the split Bregman method. Extensive experimental results, with performance evaluation, are presented to demonstrate the effectiveness of the new model and method and to show that our algorithm performs favorably in comparison to another existing method. We also discuss the influence of the order α of the fractional operator in the estimation of the optical flow, for 0 ≤ α ≤ 2 . We observe that the values of α for which the method performs better depend on the geometry and texture complexity of the image. Some extensions of our algorithm are also discussed.

      PubDate: 2016-05-07T14:35:44Z
  • Editorial Board
    • Abstract: Publication date: July 2016
      Source:Applied Numerical Mathematics, Volume 105

      PubDate: 2016-05-02T14:29:12Z
  • An adaptive algorithm based on the shifted inverse iteration for the
           Steklov eigenvalue problem
    • Abstract: Publication date: July 2016
      Source:Applied Numerical Mathematics, Volume 105
      Author(s): Hai Bi, Hao Li, Yidu Yang
      This paper proposes and analyzes an a posteriori error estimator for the finite element multi-scale discretization approximation of the Steklov eigenvalue problem. Based on the a posteriori error estimates, an adaptive algorithm of shifted inverse iteration type is designed. Finally, numerical experiments comparing the performances of three kinds of different adaptive algorithms are provided, which illustrate the efficiency of the adaptive algorithm proposed here.

      PubDate: 2016-05-02T14:29:12Z
  • A path-conservative Osher-type scheme for axially symmetric compressible
           flows in flexible visco-elastic tubes
    • Abstract: Publication date: July 2016
      Source:Applied Numerical Mathematics, Volume 105
      Author(s): Julia Leibinger, Michael Dumbser, Uwe Iben, Isabell Wayand
      Flexible tubes are widely used in modern industrial hydraulic systems as connections between different components like valves, pumps and actuators. For the design and the analysis of the temporal behavior of a hydraulic system, one therefore needs an accurate mathematical model that describes the fluid flow in a compliant duct. Hence, in this paper we want to model the fluid–structure-interaction (FSI) problem given by the axially symmetric flow of a compressible barotropic fluid that flows through flexible tubes made of vulcanized rubber. The material of the tube can be described by using a visco-elastic rheology, which takes into account the strain relaxation of the material. The resulting mathematical model consists in a one-dimensional system of nonlinear hyperbolic partial differential equations (PDE) with non-conservative products and algebraic source terms. To solve this system numerically, we apply the DOT method, which is a generalized path-conservative Osher-type Riemann solver for conservative and non-conservative hyperbolic PDE recently proposed in [23] and [22]. We provide numerical evidence that the proposed DOT Riemann solver is well-balanced for the governing PDE system under consideration. The method is compared to available quasi-exact solutions of the Riemann problem in the case of an elastic wall described by the Laplace law. It is also compared to available experimental data and exact solutions obtained in the frequency domain for a linear visco-elastic wall behavior. In all cases under investigation the proposed path-conservative finite volume scheme based on the DOT Riemann solver is able to produce very accurate results.

      PubDate: 2016-05-02T14:29:12Z
  • Efficient and accurate spectral method using generalized Jacobi functions
           for solving Riesz fractional differential equations
    • Abstract: Publication date: Available online 21 April 2016
      Source:Applied Numerical Mathematics
      Author(s): Zhiping Mao, Sheng Chen, Jie Shen
      We consider numerical approximation of the Riesz Fractional Differential Equations (FDEs), and construct a new set of generalized Jacobi functions, J n − α , − α ( x ) , which are tailored to the Riesz fractional PDEs. We develop optimal approximation results in non-uniformly weighted Sobolev spaces, and construct spectral Petrov–Galerkin algorithms to solve the Riesz FDEs with two kinds of boundary conditions (BCs): (i) homogeneous Dirichlet boundary conditions, and (ii) Integral BCs. We provide rigorous error analysis for our spectral Petrov–Galerkin methods, which show that the errors decay exponentially fast as long as the data (right-hand side function) is smooth, despite that fact that the solution has singularities at the endpoints. We also present some numerical results to validate our error analysis.

      PubDate: 2016-04-23T14:08:13Z
  • Unilateral problem for the Stokes equations: The well-posedness and finite
           element approximation
    • Abstract: Publication date: Available online 1 April 2016
      Source:Applied Numerical Mathematics
      Author(s): Norikazu Saito, Yoshiki Sugitani, Guanyu Zhou
      We consider the stationary Stokes equations under a unilateral boundary condition of Signorini's type, which is one of artificial boundary conditions in flow problems. Well-posedness is discussed through its variational inequality formulation. We also consider the finite element approximation for a regularized penalty problem. The well-posedness, stability and error estimates of optimal order are established. The lack of a coupled Babuška and Brezzi's condition makes analysis difficult. We offer a new method of analysis. Particularly, our device to treat the pressure is novel and of some interest. Numerical examples are presented to validate our theoretical results.

      PubDate: 2016-04-03T13:29:06Z
  • Series expansion method for weakly singular Volterra integral equations
    • Abstract: Publication date: Available online 30 March 2016
      Source:Applied Numerical Mathematics
      Author(s): Zhendong Gu, Xiaojing Guo, Daochun Sun
      We propose series expansion method to solve VIEs (Volterra integral equations) with smooth given functions, including weakly singular VIEs possessing unsmooth solution. The key step in proposed method is to approximate given functions by their own Chebyshev Gauss-Lobatto interpolation polynomials. Lubich's results play an important role in the proposed method for weakly singular VIEs. Convergence analysis is provided for proposed method. Numerical experiments are carried out to confirm theoretical results.

      PubDate: 2016-04-03T13:29:06Z
  • Analysis of numerical solutions to Sommerfeld integral relation of the
           half-space radiator problem
    • Abstract: Publication date: Available online 30 March 2016
      Source:Applied Numerical Mathematics
      Author(s): I. Arun, Murugesan Venkatapathi
      Sommerfeld integrals relate a spherical wave from a point source to a convolution set of plane and cylindrical waves. This relation does not have analytical solutions but it submits to a solution by numerical integration. Among others, it is significant for theoretical studies of many optical and radiation phenomena involving surfaces. This approach is preferred over discretized computational models of the surface because of the many orders of increased computations involved in the latter. One of the most widely used and accurate methods to compute these solutions is the numerical integration of the Sommerfeld integrand over a complex contour. We have analyzed the numerical advantages offered by this method, and have justified the optimality of the preferred contour of integration and the choice of two eigenfunctions used. In addition to this, we have also analyzed four other approximate methods to compute the Sommerfeld integral and have identified their regions of validity, and numerical advantages, if any. These include the high relative permittivity approximation, the short distance approximation, the exact image theory and Fourier expansion of the reflection coefficient. We also finally compare these five methods in terms of their computational cost.

      PubDate: 2016-04-03T13:29:06Z
  • Preface [Special issue: NumAn 2014]
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Georgios D. Akrivis, Vassilios A. Dougalis, Efstratios Gallopoulos, Apostolos Hadjidimos, Ilias S. Kotsireas, Dimitrios Noutsos, Yiannis G. Saridakis, Michael N. Vrahatis

      PubDate: 2016-03-29T13:02:04Z
  • Stochastic evolution equation with Riesz-fractional derivative and white
           noise on the half-line
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Martin P. Arciga Alejandre, Francisco J. Ariza Hernandez, Jorge Sanchez Ortiz
      In this work, we consider an initial boundary-value problem for a stochastic evolution equation with Riesz-fractional spatial derivative and white noise on the half-line, { u t ( x , t ) = D x α u ( x , t ) + N u ( x , t ) + B ˙ ( x , t ) , x > 0 , t ∈ [ 0 , T ] , u ( x , 0 ) = u 0 ( x ) , x > 0 , u x ( 0 , t ) = g 1 ( t ) , t ∈ [ 0 , T ] , where D x α is the Riesz-fractional derivative, α ∈ ( 2 , 3 ) , N is a Lipschitzian operator and B ˙ ( x , t ) is the white noise. To construct the integral representation of solutions we use the Fokas method and Picard scheme to prove existence and uniqueness. Moreover, Monte Carlo methods are implemented to approximate solutions.

      PubDate: 2016-03-29T13:02:04Z
  • A stopping criterion for iterative regularization methods
    • Abstract: Publication date: Available online 25 March 2016
      Source:Applied Numerical Mathematics
      Author(s): G. Landi, E. Loli Piccolomini, I. Tomba
      We present a discrepancy-like stopping criterium for iterative regularization methods for the solution of linear discrete ill-posed problems. The presented criterium terminates the iterations of the iterative method when the residual norm of the computed solution becomes less or equal to the residual norm of a regularized Truncated Singular Value Decomposition (TSVD) solution. We present two algorithms for the automatic computation of the TSVD residual norm using the Discrete Picard Condition. The first algorithm uses the SVD coefficients while the second one uses the Fourier coefficients. In this work, we mainly focus on the Conjugate Gradient Least Squares method, but the proposed criterium can be used for terminating the iterations of any iterative regularization method. Many numerical tests on some selected one dimensional and image deblurring problems are presented and the results are compared with those obtained by state-of-the-art parameter selection rules. The numerical results show the efficiency and robustness of the proposed criterium.

      PubDate: 2016-03-29T13:02:04Z
  • A fitted finite element method for the numerical approximation of void
           electro-stress migration
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Robert Nürnberg, Andrea Sacconi
      Microelectronic circuits usually contain small voids or cracks, and if those defects are large enough to sever the line, they cause an open circuit. A fully practical finite element method for the temporal analysis of the migration of voids in the presence of surface diffusion, electric loading and elastic stress is presented. We simulate a bulk–interface coupled system, with a moving interface governed by a fourth-order geometric evolution equation and a bulk where the electric potential and the displacement field are computed. The method presented here follows a fitted approach, since the interface grid is part of the boundary of the bulk grid. A detailed analysis, in terms of experimental order of convergence (when the exact solution to the free boundary problem is known) and coupling operations (e.g., smoothing/remeshing of the grids, intersection between elements of the two grids), is carried out. A comparison with a previously introduced unfitted approach (where the two grids are totally independent) is also performed, along with several numerical simulations in order to test the accuracy of the methods.

      PubDate: 2016-03-29T13:02:04Z
  • Splitting-methods based on Approximate Matrix Factorization and Radau-IIA
           formulas for the time integration of advection diffusion reaction PDEs
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): S. González-Pinto, D. Hernández-Abreu
      A family of splitting methods for the time integration of evolutionary Advection Diffusion Reaction Partial Differential Equations (PDEs) semi-discretized in space by Finite Differences is obtained. The splitting is performed in the Jacobian matrix by using the Approximate Matrix Factorization (AMF) and by considering up to three inexact Newton Iterations applied to the two-stage Radau IIA method along with a very simple predictor. The overall process allows to reduce the storage and the algebraic costs involved in the numerical solution of the multidimensional linear systems to the level of 1D-dimensional linear systems with small bandwidths. Some specific AMF-Radau methods are constructed after studying the expression for the local error in semi-linear equations, and their linear stability properties are widely studied. The wedge of stability of the methods depends on the number of splittings used for the Jacobian matrix of the spatial semidiscretized ODEs, J h = ∑ j = 1 d J h , j , where h stands for the spatial grid resolution. A-stability is proven for the cases d = 1 , 2 , and A ( 0 ) -stability for any d ≥ 1 . Numerical experiments on a 3D semi-linear advection diffusion reaction test problem and a 2D-combustion model are presented. The experiments show that the methods compare well with standard classical methods in parabolic problems and can also be successfully used for advection dominated problems when some diffusion or stiff reactions are present. In the latter case the stability imposes restrictions on the number of splitting terms (d).

      PubDate: 2016-03-29T13:02:04Z
  • Numerical integration of variational equations for Hamiltonian systems
           with long range interactions
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Helen Christodoulidi, Tassos Bountis, Lambros Drossos
      We study numerically classical 1-dimensional Hamiltonian lattices involving inter-particle long range interactions that decay with distance like 1 / r α , for α ≥ 0 . We demonstrate that although such systems are generally characterized by strong chaos, they exhibit an unexpectedly organized behavior when the exponent α < 1 . This is shown by computing dynamical quantities such as the maximal Lyapunov exponent, which decreases as the number of degrees of freedom increases. We also discuss our numerical methods of symplectic integration implemented for the solution of the equations of motion together with their associated variational equations. The validity of our numerical simulations is estimated by showing that the total energy of the system is conserved within an accuracy of 4 digits (with integration step τ = 0.02 ), even for as many as N = 8000 particles and integration times as long as 106 units.

      PubDate: 2016-03-29T13:02:04Z
  • Obviating the bin width effect of the 1/t algorithm for multidimensional
           numerical integration
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Wanyok Atisattapong, Pasin Maruphanton
      In this work we improve the accuracy and the convergence of the 1 / t algorithm for multidimensional numerical integration. The proposed strategy is to introduce a new approximation method which obviates the bin width effect of the conventional 1 / t algorithm by using the average of y values, which varies as the number of Monte Carlo trials changes, instead of the fixed value of y. The non-convergence of the 1 / t algorithm and the convergence of the new method are proved by theoretical analysis. The potential of the method is illustrated by the evaluation of one-, two- and multi- dimensional integrals up to six dimensions. Our results show that the numerical estimates from our method converge to their exact values without either error saturation or the bin with effect, in contrast with the conventional 1 / t algorithm.

      PubDate: 2016-03-29T13:02:04Z
  • Numerical analysis of Fokas' unified method for linear elliptic PDEs
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): A.C.L. Ashton, K.M. Crooks
      In this paper we examine a numerical implementation of Fokas' unified method for elliptic boundary value problems on convex polygons. Within this setting the unified method provides a reconstruction of the unknown boundary values from the known boundary data for an important class of elliptic boundary value problems. We demonstrate that this approach can be used to study the classical Dirichlet eigenvalue problem for the Laplacian on convex polygons. We show that this problem is equivalent to a nonlinear eigenvalue problem for a semi-Fredholm operator which has holomorphic dependence on the eigenvalue itself. We also study the classical Dirichlet problem for the Helmholtz operator. We provide a new Galerkin scheme for the underlying Dirichlet–Neumann map, and prove that for sufficiently regular Dirichlet data this scheme converges with spectral accuracy.

      PubDate: 2016-03-29T13:02:04Z
  • Coupled symplectic maps as models for subdiffusive processes in disordered
           Hamiltonian lattices
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Chris G. Antonopoulos, Tassos Bountis, Lambros Drossos
      We investigate dynamically and statistically diffusive motion in a chain of linearly coupled 2-dimensional symplectic McMillan maps and find evidence of subdiffusion in weakly and strongly chaotic regimes when all maps of the chain possess a saddle point at the origin and the central map is initially excited. In the case of weak coupling, there is either absence of diffusion or subdiffusion with q > 1 -Gaussian probability distributions, characterizing weak chaos. However, for large enough coupling and already moderate number of maps, the system exhibits strongly chaotic ( q ≈ 1 ) subdiffusive behavior, reminiscent of the subdiffusive energy spreading observed in a disordered Klein–Gordon Hamiltonian. Our results provide evidence that coupled symplectic maps can exhibit physical properties similar to those of disordered Hamiltonian systems, even though the local dynamics in the two cases is significantly different.

      PubDate: 2016-03-29T13:02:04Z
  • Numerical modeling of sediment transport applied to coastal morphodynamics
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): G.V. Kozyrakis, A.I. Delis, G. Alexandrakis, N.A. Kampanis
      A bed-load sediment transport model is used to describe realistic cases of the morphodynamics in coastal areas. The hydrodynamic equations are based on the well-known, two-dimensional depth-averaged non-linear shallow water equations, with bathymetry forces and friction, which are subsequently coupled to the Exner equation to describe the morphological evolution. Different forms of the bed-load transport flux are considered in the Exner equation and certain relations between them are established. The numerical model is expressed in a fully-coupled form where a single system of equations is solved by a high-resolution two-dimensional finite volume scheme of the relaxation type. The relaxation is performed by classical models where neither approximate Riemann solvers nor characteristic decompositions are needed. The overall numerical scheme is validated in benchmark problems, and for a realistic application of a coastal area on the northern side of the island of Crete. The validity of these results is established by comparisons made with the well-known MIKE Software by DHI Group.

      PubDate: 2016-03-29T13:02:04Z
  • Fokas transform method for a brain tumor invasion model with heterogeneous
           diffusion in 1+1 dimensions
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): D. Mantzavinos, M.G. Papadomanolaki, Y.G. Saridakis, A.G. Sifalakis
      Gliomas are among the most aggressive forms of brain tumors. Over the last years mathematical models have been well developed to study gliomas growth. We consider a simple and well established mathematical model focused on proliferation and diffusion. Due to the heterogeneity of the brain tissue (white and grey matter) the diffusion coefficient is considered to be discontinuous. Fokas transform approach for the solution of linear PDE problems, apart from the fact that it avoids solving intermediate ODE problems, yields novel integral representations of the solution in the complex plane that decay exponentially fast and converge uniformly at the boundaries. To take advantage of these properties for the solution of the model problem at hand, we have successfully implemented Fokas transform method in the multi-domain environment induced by the interface discontinuities of our problem's domain. The fact that the integral representation of the solution at any time–space point of our problem's domain is independent on any other points of the domain, except of course on initial data, coupled with a simple composite trapezoidal rule, implemented on appropriately chosen integration contours, yields a fast and efficient analytical–numerical technique capable of producing directly high-order approximations of the solution at any point of the domain requiring no prior knowledge of the solution at any other time instances or space information.

      PubDate: 2016-03-29T13:02:04Z
  • Scientific computations on multi-core systems using different programming
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Panagiotis D. Michailidis, Konstantinos G. Margaritis
      Numerical linear algebra is one of the most important forms of scientific computation. The basic computations in numerical linear algebra are matrix computations and linear systems solution. These computations are used as kernels in many computational problems. This study demonstrates the parallelisation of these scientific computations using multi-core programming frameworks. Specifically, the frameworks examined here are Pthreads, OpenMP, Intel Cilk Plus, Intel TBB, SWARM, and FastFlow. A unified and exploratory performance evaluation and a qualitative study of these frameworks are also presented for parallel scientific computations with several parameters. The OpenMP and SWARM models produce good results running in parallel with compiler optimisation when implementing matrix operations at large and medium scales, whereas the remaining models do not perform as well for some matrix operations. The qualitative results show that the OpenMP, Cilk Plus, TBB, and SWARM frameworks require minimal programming effort, whereas the other models require advanced programming skills and experience. Finally, based on an extended study, general conclusions regarding the programming models and matrix operations for some parameters were obtained.

      PubDate: 2016-03-29T13:02:04Z
  • Multigrid methods for cubic spline solution of two point (and 2D) boundary
           value problems
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Marco Donatelli, Matteo Molteni, Vincenzo Pennati, Stefano Serra-Capizzano
      In this paper we propose a scheme based on cubic splines for the solution of the second order two point boundary value problems. The solution of the algebraic system is computed by using optimized multigrid methods. In particular the transformation of the stiffness matrix essentially in a block Toeplitz matrix and its spectral analysis allow to choose smoothers able to reduce error components related to the various frequencies and to obtain an optimal method. The main advantages of our strategy can be listed as follows: (i) a fourth order of accuracy combined with a quadratic conditioning matrix, (ii) a resulting matrix structure whose eigenvalues can be compactly described by a symbol (this information is the key for designing an optimal multigrid method). Finally, some numerics that confirm the predicted behavior of the method are presented and a discussion on the two dimensional case is given, together with few 2D numerical experiments.

      PubDate: 2016-03-29T13:02:04Z
  • Discontinuous Galerkin methods for fast reactive mass transfer through
           semi-permeable membranes
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Andrea Cangiani, Emmanuil H. Georgoulis, Max Jensen
      A discontinuous Galerkin (dG) method for the numerical solution of initial/boundary value multi-compartment partial differential equation (PDE) models, interconnected with interface conditions, is analysed. The study of interface problems is motivated by models of mass transfer of solutes through semi-permeable membranes. The case of fast reactions is also included. More specifically, a model problem consisting of a system of semilinear parabolic advection–diffusion–reaction partial differential equations in each compartment with only local Lipschitz conditions on the nonlinear reaction terms, equipped with respective initial and boundary conditions, is considered. General nonlinear interface conditions modelling selective permeability, congestion and partial reflection are applied to the compartment interfaces. The interior penalty dG method for this problem, presented recently, is analysed both in the space-discrete and in fully discrete settings for the case of, possibly, fast reactions. The a priori analysis shows that the method yields optimal a priori bounds, provided the exact solution is sufficiently smooth. Numerical experiments indicate agreement with the theoretical bounds.

      PubDate: 2016-03-29T13:02:04Z
  • Editorial Board
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104

      PubDate: 2016-03-29T13:02:04Z
  • Preface [Special issue: NumAn 2012]
    • Abstract: Publication date: June 2016
      Source:Applied Numerical Mathematics, Volume 104
      Author(s): Apostolos Hadjidimos, Ilias S. Kotsireas, Dimitrios Noutsos, Michael N. Vrahatis

      PubDate: 2016-03-29T13:02:04Z
  • Sensitivity analysis for active control of the Helmholtz equation
    • Abstract: Publication date: Available online 21 March 2016
      Source:Applied Numerical Mathematics
      Author(s): M. Hubenthal, D. Onofrei
      In previous works we considered the Helmholtz equation with fixed frequency k outside a discrete set of resonant frequencies, where it is implied that, given a source region D a ⊂ R d ( d = 2 , 3 ‾ ) and u 0 , a solution of the homogeneous scalar Helmholtz equation in a set containing the control region D c ⊂ R d , there exists an infinite class of boundary data on ∂ D a so that the radiating solution to the corresponding exterior scalar Helmholtz problem in R d ∖ D a will closely approximate u 0 in D c . Moreover, it will have vanishingly small values beyond a certain large enough “far-field” radius R. In this paper we study the minimal energy solution of the above problem (e.g. the solution obtained by using Tikhonov regularization with the Morozov discrepancy principle) and perform a detailed sensitivity analysis. In this regard we discuss the stability of the minimal energy solution with respect to measurement errors as well as the feasibility of the active scheme (power budget and accuracy) depending on: the mutual distances between the antenna, control region and far field radius R; value of the regularization parameter; frequency; location of the source.

      PubDate: 2016-03-25T12:42:05Z
  • Inverse coefficient problems for a first order hyperbolic system
    • Abstract: Publication date: Available online 22 March 2016
      Source:Applied Numerical Mathematics
      Author(s): Mansur I. Ismailov, Ibrahim Tekin
      In this paper, the direct and inverse initial boundary value problems for a first order system of two hyperbolic equations are considered. The method of characteristics and the finite difference method are applied to the theoretical and numerical solutions of the direct problem, respectively. Moreover the suitability of the method of characteristics for the inverse problem of finding solely space-dependent coefficients and the finite difference method for solely time-dependent coefficients of the first order hyperbolic system are shown. The stability of the numerical method is supported by the examples.

      PubDate: 2016-03-25T12:42:05Z
  • Compensation of domain modelling errors in the inverse source problem of
           the Poisson equation: Application in electroencephalographic imaging
    • Abstract: Publication date: Available online 21 March 2016
      Source:Applied Numerical Mathematics
      Author(s): Alexandra Koulouri, Ville Rimpiläinen, Mike Brookes, Jari P. Kaipio
      In the inverse source problem of the Poisson equation, measurements on the domain boundaries are used to reconstruct sources inside the domain. The problem is an ill-posed inverse problem and it is sensitive to modelling errors of the domain. These errors can be boundary, structure and material property errors, for example. In this paper, we investigate whether the recently proposed Bayesian approximation error (BAE) approach could be used to alleviate the source estimation errors when an approximate model for the domain is employed. The BAE is based on postulating a probabilistic model for the uncertainties, in this case the geometry and structure of the domain, and to carry out approximate marginalization over these nuisance parameters. We particularly consider electroencephalography (EEG) source imaging as an application. EEG is a diagnostic brain imaging modality, and it can be used to reconstruct neural sources in the brain from electric potential measurements along the scalp. In the feasibility study, we assess to which degree one can recover from the modelling errors that are induced by the use of the three concentric circle head model instead of an anatomically accurate head model. The studied domain modelling errors include errors in the geometry of the exterior boundary and the structure of the interior. We show that, in particular with superficial dipole sources, the BAE yields estimates that can in some cases be considered adequately accurate. This would avoid the need for the extraction of the accurate head features which is conventionally carried out via expensive and time consuming auxiliary imaging modalities such as magnetic resonance imaging.

      PubDate: 2016-03-25T12:42:05Z
  • Generalized quadrature for solving singular integral equations of Abel
           type in application to infrared tomography
    • Abstract: Publication date: Available online 24 March 2016
      Source:Applied Numerical Mathematics
      Author(s): V.S. Sizikov, D.N. Sidorov
      We propose the generalized quadrature methods for numerical solution of singular integral equation of Abel type. We overcome the singularity using the analytic computation of the singular integral. The problem of solution of singular integral equation is reduced to nonsingular system of linear algebraic equations without shift meshes techniques employment. We also propose generalized quadrature method for solution of Abel equation using the singular integral. Relaxed errors bounds are derived. In order to improve the accuracy we use Tikhonov regularization method. We demonstrate the efficiency of proposed techniques on infrared tomography problem. Numerical experiments show that it make sense to apply regularization in case of highly noisy (about 10%) sources only. That is due to the fact that singular integral equations enjoy selfregularization property.

      PubDate: 2016-03-25T12:42:05Z
  • On the robustness of ILU smoothers on triangular grids
    • Abstract: Publication date: Available online 22 March 2016
      Source:Applied Numerical Mathematics
      Author(s): M.A.V. Pinto, C. Rodrigo, F.J. Gaspar, C.W. Oosterlee
      In this work, incomplete factorization techniques are used as smoothers within a geometric multigrid algorithm on triangular grids. A local Fourier analysis is proposed to study the smoothing properties of these methods, as well as the asymptotic convergence of the whole multigrid procedure. With this purpose, two- and three-grid local Fourier analysis are performed. Several two-dimensional diffusion problems, including different kinds of anisotropy are considered to demonstrate the robustness of this type of methods.

      PubDate: 2016-03-25T12:42:05Z
  • The linear sampling method for inhomogeneous medium and buried objects
           from far field measurements
    • Abstract: Publication date: Available online 2 March 2016
      Source:Applied Numerical Mathematics
      Author(s): Haihua Qin, Xiaodong Liu
      We are concerned with the reconstruction of both the penetrable inhomogeneous medium and the buried impenetrable obstacle. Firstly, the classical linear sampling method is used to recover the support of the inhomogeneous medium, and then a modification of the linear sampling method is proposed for objects buried in a known layered medium. The main feature of our method is that it avoids using knowledge of the Green's function for the background media. Finally, some numerical experiments are presented to demonstrate the feasibility and effectiveness of our method.

      PubDate: 2016-03-08T22:40:54Z
  • On a local Fourier analysis for overlapping block smoothers on triangular
    • Abstract: Publication date: Available online 3 March 2016
      Source:Applied Numerical Mathematics
      Author(s): C. Rodrigo, F.J. Gaspar, F.J. Lisbona
      A general local Fourier analysis for overlapping block smoothers on triangular grids is presented. This analysis is explained in a general form for its application to problems with different discretizations. This tool is demonstrated for two different problems: a stabilized linear finite element discretization of Stokes equations and an edge-based discretization of the curl-curl operator by lowest-order Nédélec finite element method. In this latter, special Fourier modes have to be considered in order to perform the analysis. Numerical results comparing two- and three-grid convergence factors predicted by the local Fourier analysis to real asymptotic convergence factors are presented to confirm the predictions of the analysis and show their usefulness.

      PubDate: 2016-03-08T22:40:54Z
  • A hybrid recursive multilevel incomplete factorization preconditioner for
           solving general linear systems
    • Abstract: Publication date: Available online 15 February 2016
      Source:Applied Numerical Mathematics
      Author(s): Yiming Bu, Bruno Carpentieri, Zhaoli Shen, Tingzhu Huang
      In this paper we introduce an algebraic recursive multilevel incomplete factorization preconditioner, based on a distributed Schur complement formulation, for solving general linear systems. The novelty of the proposed method is to combine factorization techniques of both implicit and explicit type, recursive combinatorial algorithms, multilevel mechanisms and overlapping strategies to maximize sparsity in the inverse factors and consequently reduce the factorization costs. Numerical experiments demonstrate the good potential of the proposed solver to precondition effectively general linear systems, also against other state-of-the-art iterative solvers of both implicit and explicit form.

      PubDate: 2016-02-15T21:48:22Z
  • An hp-spectral collocation method for nonlinear Volterra functional
           integro-differential equations with delays
    • Abstract: Publication date: July 2016
      Source:Applied Numerical Mathematics, Volume 105
      Author(s): Chang-tao Sheng, Zhong-qing Wang, Ben-yu Guo
      In this paper, we propose a multistep Legendre–Gauss spectral collocation method for the nonlinear Volterra functional integro-differential equations (VFIDEs) with vanishing delays. This method is easy to implement and possesses the high-order accuracy. We also provide a rigorous convergence analysis of the hp-version of the multistep spectral collocation method under H 1 -norm. Numerical results confirm the theoretical analysis.

      PubDate: 2016-02-12T06:36:05Z
  • A C1-conforming hp finite element method for fourth order singularly
           perturbed boundary value problems
    • Abstract: Publication date: Available online 10 February 2016
      Source:Applied Numerical Mathematics
      Author(s): Pandelitsa Panaseti, Antri Zouvani, Niall Madden, Christos Xenophontos
      We consider a fourth order singularly perturbed boundary value problem (BVP) in one-dimension and the approximation of its solution by the hp version of the Finite Element Method (FEM). The given problem's boundary conditions are not suitable for writing the BVP as a second order system, hence the approximation will be sought from a finite dimensional subspace of the Sobolev space H 2 . We construct suitable C 1 hierarchical basis functions for the approximation and we show that the hp FEM on the Spectral Boundary Layer Mesh yields a robust approximation that converges exponentially in the energy norm, as the number of degrees of freedom is increased. Numerical examples that validate (and extend) the theory are also presented.

      PubDate: 2016-02-12T06:36:05Z
  • Numerical solution of optimization problems for semi-linear elliptic
           equations with discontinuous coefficients and solutions
    • Abstract: Publication date: Available online 11 February 2016
      Source:Applied Numerical Mathematics
      Author(s): A.R. Manapova, F.V. Lubyshev
      In this work we consider optimization problems for processes described by semi-linear partial differential equations of elliptic type with discontinuous coefficients and solutions (with imperfect contact matching conditions), with controls involved in the coefficients. Finite difference approximations of optimization problems are constructed. For the numerical implementation of finite optimization problems differentiability and Lipshitz-continuity of the grid functional of the approximating grid problems are proved. An iterative method for solving boundary value problems of contact for PDEs of elliptic type with discontinuous coefficients and solutions is developed and validated. The convergence of the iterative process is investigated. And the convergence rate of iterations (with calculated constants) is estimated.

      PubDate: 2016-02-12T06:36:05Z
  • Editorial Board
    • Abstract: Publication date: May 2016
      Source:Applied Numerical Mathematics, Volume 103

      PubDate: 2016-02-12T06:36:05Z
  • An NFFT based approach to the efficient computation of dipole–dipole
           interactions under various periodic boundary conditions
    • Abstract: Publication date: July 2016
      Source:Applied Numerical Mathematics, Volume 105
      Author(s): Franziska Nestler
      We present an efficient method to compute the electrostatic fields, torques and forces in dipolar systems, which is based on the fast Fourier transform for nonequispaced data (NFFT). We consider 3d-periodic, 2d-periodic, 1d-periodic as well as 0d-periodic (open) boundary conditions. The method is based on the corresponding Ewald formulas, which immediately lead to an efficient algorithm only in the 3d-periodic case. In the other cases we apply the NFFT based fast summation in order to approximate the contributions of the nonperiodic dimensions in Fourier space. This is done by regularizing or periodizing the involved functions, which depend on the distances of the particles regarding the nonperiodic dimensions. The final algorithm enables a unified treatment of all types of periodic boundary conditions, for which only the precomputation step has to be adjusted.

      PubDate: 2016-02-12T06:36:05Z
  • Numerical Fourier method and second-order Taylor scheme for backward SDEs
           in finance
    • Abstract: Publication date: Available online 15 January 2016
      Source:Applied Numerical Mathematics
      Author(s): M.J. Ruijter, C.W. Oosterlee
      We develop a Fourier method to solve quite general backward stochastic differential equations (BSDEs) with second-order accuracy. The underlying forward stochastic differential equation (FSDE) is approximated by different Taylor schemes, such as the Euler, Milstein, and Order 2.0 weak Taylor schemes, or by exact simulation. A θ-time-discretization of the time-integrands leads to an induction scheme with conditional expectations. The computation of the conditional expectations appearing relies on the availability of the characteristic function for these schemes. We will use the characteristic function of the discrete forward process. The expected values are approximated by Fourier cosine series expansions. Numerical experiments show rapid convergence of our efficient probabilistic numerical method. Second-order accuracy is observed and also proved. We apply the method to, among others, option pricing problems under the Constant Elasticity of Variance and Cox-Ingersoll-Ross processes.

      PubDate: 2016-01-18T11:52:43Z
School of Mathematical and Computer Sciences
Heriot-Watt University
Edinburgh, EH14 4AS, UK
Tel: +00 44 (0)131 4513762
Fax: +00 44 (0)131 4513327
About JournalTOCs
News (blog, publications)
JournalTOCs on Twitter   JournalTOCs on Facebook

JournalTOCs © 2009-2015