for Journals by Title or ISSN
for Articles by Keywords
  Subjects -> COMPUTER SCIENCE (Total: 2072 journals)
    - ANIMATION AND SIMULATION (31 journals)
    - ARTIFICIAL INTELLIGENCE (102 journals)
    - AUTOMATION AND ROBOTICS (105 journals)
    - COMPUTER ARCHITECTURE (10 journals)
    - COMPUTER ENGINEERING (11 journals)
    - COMPUTER GAMES (21 journals)
    - COMPUTER PROGRAMMING (26 journals)
    - COMPUTER SCIENCE (1202 journals)
    - COMPUTER SECURITY (46 journals)
    - DATA BASE MANAGEMENT (14 journals)
    - DATA MINING (36 journals)
    - E-BUSINESS (22 journals)
    - E-LEARNING (29 journals)
    - IMAGE AND VIDEO PROCESSING (40 journals)
    - INFORMATION SYSTEMS (107 journals)
    - INTERNET (93 journals)
    - SOCIAL WEB (51 journals)
    - SOFTWARE (33 journals)
    - THEORY OF COMPUTING (8 journals)

COMPUTER SCIENCE (1202 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: 20)
Abakós     Open Access   (Followers: 4)
ACM Computing Surveys     Hybrid Journal   (Followers: 28)
ACM Journal on Computing and Cultural Heritage     Hybrid Journal   (Followers: 8)
ACM Journal on Emerging Technologies in Computing Systems     Hybrid Journal   (Followers: 14)
ACM Transactions on Accessible Computing (TACCESS)     Hybrid Journal   (Followers: 3)
ACM Transactions on Algorithms (TALG)     Hybrid Journal   (Followers: 15)
ACM Transactions on Applied Perception (TAP)     Hybrid Journal   (Followers: 5)
ACM Transactions on Architecture and Code Optimization (TACO)     Hybrid Journal   (Followers: 9)
ACM Transactions on Autonomous and Adaptive Systems (TAAS)     Hybrid Journal   (Followers: 9)
ACM Transactions on Computation Theory (TOCT)     Hybrid Journal   (Followers: 12)
ACM Transactions on Computational Logic (TOCL)     Hybrid Journal   (Followers: 3)
ACM Transactions on Computer Systems (TOCS)     Hybrid Journal   (Followers: 18)
ACM Transactions on Computer-Human Interaction     Hybrid Journal   (Followers: 15)
ACM Transactions on Computing Education (TOCE)     Hybrid Journal   (Followers: 6)
ACM Transactions on Design Automation of Electronic Systems (TODAES)     Hybrid Journal   (Followers: 5)
ACM Transactions on Economics and Computation     Hybrid Journal   (Followers: 1)
ACM Transactions on Embedded Computing Systems (TECS)     Hybrid Journal   (Followers: 4)
ACM Transactions on Information Systems (TOIS)     Hybrid Journal   (Followers: 19)
ACM Transactions on Intelligent Systems and Technology (TIST)     Hybrid Journal   (Followers: 8)
ACM Transactions on Interactive Intelligent Systems (TiiS)     Hybrid Journal   (Followers: 3)
ACM Transactions on Multimedia Computing, Communications, and Applications (TOMCCAP)     Hybrid Journal   (Followers: 9)
ACM Transactions on Reconfigurable Technology and Systems (TRETS)     Hybrid Journal   (Followers: 6)
ACM Transactions on Sensor Networks (TOSN)     Hybrid Journal   (Followers: 8)
ACM Transactions on Speech and Language Processing (TSLP)     Hybrid Journal   (Followers: 9)
ACM Transactions on Storage     Hybrid Journal  
ACS Applied Materials & Interfaces     Hybrid Journal   (Followers: 31)
Acta Automatica Sinica     Full-text available via subscription   (Followers: 2)
Acta Informatica Malaysia     Open Access  
Acta Universitatis Cibiniensis. Technical Series     Open Access  
Ad Hoc Networks     Hybrid Journal   (Followers: 11)
Adaptive Behavior     Hybrid Journal   (Followers: 11)
Advanced Engineering Materials     Hybrid Journal   (Followers: 28)
Advanced Science Letters     Full-text available via subscription   (Followers: 10)
Advances in Adaptive Data Analysis     Hybrid Journal   (Followers: 7)
Advances in Artificial Intelligence     Open Access   (Followers: 15)
Advances in Calculus of Variations     Hybrid Journal   (Followers: 4)
Advances in Catalysis     Full-text available via subscription   (Followers: 5)
Advances in Computational Mathematics     Hybrid Journal   (Followers: 19)
Advances in Computer Engineering     Open Access   (Followers: 4)
Advances in Computer Science : an International Journal     Open Access   (Followers: 14)
Advances in Computing     Open Access   (Followers: 2)
Advances in Data Analysis and Classification     Hybrid Journal   (Followers: 54)
Advances in Engineering Software     Hybrid Journal   (Followers: 28)
Advances in Geosciences (ADGEO)     Open Access   (Followers: 14)
Advances in Human Factors/Ergonomics     Full-text available via subscription   (Followers: 23)
Advances in Human-Computer Interaction     Open Access   (Followers: 20)
Advances in Materials Science     Open Access   (Followers: 14)
Advances in Operations Research     Open Access   (Followers: 12)
Advances in Parallel Computing     Full-text available via subscription   (Followers: 7)
Advances in Porous Media     Full-text available via subscription   (Followers: 5)
Advances in Remote Sensing     Open Access   (Followers: 49)
Advances in Science and Research (ASR)     Open Access   (Followers: 6)
Advances in Technology Innovation     Open Access   (Followers: 6)
AEU - International Journal of Electronics and Communications     Hybrid Journal   (Followers: 8)
African Journal of Information and Communication     Open Access   (Followers: 9)
African Journal of Mathematics and Computer Science Research     Open Access   (Followers: 4)
AI EDAM     Hybrid Journal  
Air, Soil & Water Research     Open Access   (Followers: 12)
AIS Transactions on Human-Computer Interaction     Open Access   (Followers: 6)
Algebras and Representation Theory     Hybrid Journal   (Followers: 1)
Algorithms     Open Access   (Followers: 11)
American Journal of Computational and Applied Mathematics     Open Access   (Followers: 5)
American Journal of Computational Mathematics     Open Access   (Followers: 4)
American Journal of Information Systems     Open Access   (Followers: 5)
American Journal of Sensor Technology     Open Access   (Followers: 4)
Anais da Academia Brasileira de Ciências     Open Access   (Followers: 2)
Analog Integrated Circuits and Signal Processing     Hybrid Journal   (Followers: 7)
Analysis in Theory and Applications     Hybrid Journal   (Followers: 1)
Animation Practice, Process & Production     Hybrid Journal   (Followers: 5)
Annals of Combinatorics     Hybrid Journal   (Followers: 4)
Annals of Data Science     Hybrid Journal   (Followers: 12)
Annals of Mathematics and Artificial Intelligence     Hybrid Journal   (Followers: 12)
Annals of Pure and Applied Logic     Open Access   (Followers: 3)
Annals of Software Engineering     Hybrid Journal   (Followers: 13)
Annual Reviews in Control     Hybrid Journal   (Followers: 8)
Anuario Americanista Europeo     Open Access  
Applicable Algebra in Engineering, Communication and Computing     Hybrid Journal   (Followers: 2)
Applied and Computational Harmonic Analysis     Full-text available via subscription   (Followers: 1)
Applied Artificial Intelligence: An International Journal     Hybrid Journal   (Followers: 12)
Applied Categorical Structures     Hybrid Journal   (Followers: 2)
Applied Computational Intelligence and Soft Computing     Open Access   (Followers: 12)
Applied Computer Systems     Open Access   (Followers: 2)
Applied Informatics     Open Access  
Applied Mathematics and Computation     Hybrid Journal   (Followers: 33)
Applied Medical Informatics     Open Access   (Followers: 10)
Applied Numerical Mathematics     Hybrid Journal   (Followers: 5)
Applied Soft Computing     Hybrid Journal   (Followers: 16)
Applied Spatial Analysis and Policy     Hybrid Journal   (Followers: 5)
Applied System Innovation     Open Access  
Architectural Theory Review     Hybrid Journal   (Followers: 3)
Archive of Applied Mechanics     Hybrid Journal   (Followers: 5)
Archive of Numerical Software     Open Access  
Archives and Museum Informatics     Hybrid Journal   (Followers: 144)
Archives of Computational Methods in Engineering     Hybrid Journal   (Followers: 5)
arq: Architectural Research Quarterly     Hybrid Journal   (Followers: 8)
Artifact     Hybrid Journal   (Followers: 2)
Artificial Life     Hybrid Journal   (Followers: 7)
Asia Pacific Journal on Computational Engineering     Open Access  
Asia-Pacific Journal of Information Technology and Multimedia     Open Access   (Followers: 1)
Asian Journal of Computer Science and Information Technology     Open Access  
Asian Journal of Control     Hybrid Journal   (Followers: 1)
Assembly Automation     Hybrid Journal   (Followers: 2)
at - Automatisierungstechnik     Hybrid Journal   (Followers: 1)
Australian Educational Computing     Open Access   (Followers: 1)
Automatic Control and Computer Sciences     Hybrid Journal   (Followers: 5)
Automatic Documentation and Mathematical Linguistics     Hybrid Journal   (Followers: 5)
Automatica     Hybrid Journal   (Followers: 12)
Automation in Construction     Hybrid Journal   (Followers: 6)
Autonomous Mental Development, IEEE Transactions on     Hybrid Journal   (Followers: 9)
Basin Research     Hybrid Journal   (Followers: 5)
Behaviour & Information Technology     Hybrid Journal   (Followers: 52)
Big Data and Cognitive Computing     Open Access   (Followers: 2)
Biodiversity Information Science and Standards     Open Access  
Bioinformatics     Hybrid Journal   (Followers: 302)
Biomedical Engineering     Hybrid Journal   (Followers: 15)
Biomedical Engineering and Computational Biology     Open Access   (Followers: 13)
Biomedical Engineering, IEEE Reviews in     Full-text available via subscription   (Followers: 21)
Biomedical Engineering, IEEE Transactions on     Hybrid Journal   (Followers: 37)
Briefings in Bioinformatics     Hybrid Journal   (Followers: 48)
British Journal of Educational Technology     Hybrid Journal   (Followers: 142)
Broadcasting, IEEE Transactions on     Hybrid Journal   (Followers: 12)
c't Magazin fuer Computertechnik     Full-text available via subscription   (Followers: 1)
CALCOLO     Hybrid Journal  
Calphad     Hybrid Journal   (Followers: 2)
Canadian Journal of Electrical and Computer Engineering     Full-text available via subscription   (Followers: 15)
Capturing Intelligence     Full-text available via subscription  
Catalysis in Industry     Hybrid Journal   (Followers: 1)
CEAS Space Journal     Hybrid Journal   (Followers: 2)
Cell Communication and Signaling     Open Access   (Followers: 2)
Central European Journal of Computer Science     Hybrid Journal   (Followers: 5)
CERN IdeaSquare Journal of Experimental Innovation     Open Access   (Followers: 3)
Chaos, Solitons & Fractals     Hybrid Journal   (Followers: 3)
Chemometrics and Intelligent Laboratory Systems     Hybrid Journal   (Followers: 14)
ChemSusChem     Hybrid Journal   (Followers: 7)
China Communications     Full-text available via subscription   (Followers: 7)
Chinese Journal of Catalysis     Full-text available via subscription   (Followers: 2)
CIN Computers Informatics Nursing     Full-text available via subscription   (Followers: 11)
Circuits and Systems     Open Access   (Followers: 15)
Clean Air Journal     Full-text available via subscription   (Followers: 1)
CLEI Electronic Journal     Open Access  
Clin-Alert     Hybrid Journal   (Followers: 1)
Cluster Computing     Hybrid Journal   (Followers: 1)
Cognitive Computation     Hybrid Journal   (Followers: 4)
COMBINATORICA     Hybrid Journal  
Combinatorics, Probability and Computing     Hybrid Journal   (Followers: 4)
Combustion Theory and Modelling     Hybrid Journal   (Followers: 14)
Communication Methods and Measures     Hybrid Journal   (Followers: 12)
Communication Theory     Hybrid Journal   (Followers: 22)
Communications Engineer     Hybrid Journal   (Followers: 1)
Communications in Algebra     Hybrid Journal   (Followers: 3)
Communications in Computational Physics     Full-text available via subscription   (Followers: 2)
Communications in Information Science and Management Engineering     Open Access   (Followers: 4)
Communications in Partial Differential Equations     Hybrid Journal   (Followers: 3)
Communications of the ACM     Full-text available via subscription   (Followers: 51)
Communications of the Association for Information Systems     Open Access   (Followers: 16)
COMPEL: The International Journal for Computation and Mathematics in Electrical and Electronic Engineering     Hybrid Journal   (Followers: 3)
Complex & Intelligent Systems     Open Access   (Followers: 1)
Complex Adaptive Systems Modeling     Open Access  
Complex Analysis and Operator Theory     Hybrid Journal   (Followers: 2)
Complexity     Hybrid Journal   (Followers: 6)
Complexus     Full-text available via subscription  
Composite Materials Series     Full-text available via subscription   (Followers: 8)
Computación y Sistemas     Open Access  
Computation     Open Access   (Followers: 1)
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: 2)
Computational and Theoretical Chemistry     Hybrid Journal   (Followers: 9)
Computational Astrophysics and Cosmology     Open Access   (Followers: 1)
Computational Biology and Chemistry     Hybrid Journal   (Followers: 12)
Computational Chemistry     Open Access   (Followers: 2)
Computational Cognitive Science     Open Access   (Followers: 2)
Computational Complexity     Hybrid Journal   (Followers: 4)
Computational Condensed Matter     Open Access  
Computational Ecology and Software     Open Access   (Followers: 9)
Computational Economics     Hybrid Journal   (Followers: 9)
Computational Geosciences     Hybrid Journal   (Followers: 17)
Computational Linguistics     Open Access   (Followers: 23)
Computational Management Science     Hybrid Journal  
Computational Mathematics and Modeling     Hybrid Journal   (Followers: 8)
Computational Mechanics     Hybrid Journal   (Followers: 5)
Computational Methods and Function Theory     Hybrid Journal  
Computational Molecular Bioscience     Open Access   (Followers: 2)
Computational Optimization and Applications     Hybrid Journal   (Followers: 7)
Computational Particle Mechanics     Hybrid Journal   (Followers: 1)
Computational Research     Open Access   (Followers: 1)
Computational Science and Discovery     Full-text available via subscription   (Followers: 2)
Computational Science and Techniques     Open Access  
Computational Statistics     Hybrid Journal   (Followers: 14)
Computational Statistics & Data Analysis     Hybrid Journal   (Followers: 30)
Computer     Full-text available via subscription   (Followers: 98)
Computer Aided Surgery     Open Access   (Followers: 6)
Computer Applications in Engineering Education     Hybrid Journal   (Followers: 8)
Computer Communications     Hybrid Journal   (Followers: 16)
Computer Journal     Hybrid Journal   (Followers: 9)
Computer Methods in Applied Mechanics and Engineering     Hybrid Journal   (Followers: 24)
Computer Methods in Biomechanics and Biomedical Engineering     Hybrid Journal   (Followers: 12)
Computer Methods in the Geosciences     Full-text available via subscription   (Followers: 2)

        1 2 3 4 5 6 7 | Last

Journal Cover
Applied Numerical Mathematics
Journal Prestige (SJR): 0.93
Citation Impact (citeScore): 1
Number of Followers: 5  
  Hybrid Journal Hybrid journal (It can contain Open Access articles)
ISSN (Print) 0168-9274 - ISSN (Online) 0168-9274
Published by Elsevier Homepage  [3161 journals]
  • Stochastic multi-symplectic Runge-Kutta methods for stochastic Hamiltonian
    • Abstract: Publication date: Available online 19 September 2018Source: Applied Numerical MathematicsAuthor(s): Liying Zhang, Lihai JiAbstractIn this paper, we consider stochastic Runge-Kutta methods for stochastic Hamiltonian partial differential equations and present some sufficient conditions for stochastic multi-symplecticity of stochastic Runge-Kutta methods. To present more clearly, we apply these ideas to three dimensional stochastic Maxwell equations driven by multiplicative noise, which play an important role in stochastic electromagnetism and statistical radiophysics areas. Theoretical analysis shows that the methods inherit the energy conservation law of the original system, and preserve the discrete stochastic multi-symplectic conservation law almost surely.
  • Symplectic and symmetric trigonometrically-fitted ARKN methods
    • Abstract: Publication date: Available online 19 September 2018Source: Applied Numerical MathematicsAuthor(s): Jiyong LiAbstractAdapted Runge–Kutta–Nyström (ARKN) methods for solving oscillatory problems q″(t)+w2q(t)=f(q(t),q′(t)) have been investigated by several authors. Recently, Yang et al. [Applied Numerical Mathematics 58 (2008) 1375-1395] proposed trigonometrically-fitted ARKN (TFARKN) methods by introducing frequency depending coefficients into the terms in the internal stages. In applications the function f(q) often does not contain q′ explicitly and satisfies f(q)−w2q=−∇V(q) for some smooth function V(q). Then the problem can be considered as a separable Hamiltonian system. In this paper we investigate the symplecticity and symmetry of TFARKN methods for separable Hamiltonian systems and derive necessary and sufficient conditions for an TFARKN method to be symplectic and symmetric. Based on these conditions, two explicit symplectic and symmetric TFARKN method with order two and four, respectively, are constructed. Some numerical experiments are provided to confirm the theoretical expectations.
  • Convergence and stability of block boundary value methods applied to
           nonlinear fractional differential equations with Caputo derivatives
    • Abstract: Publication date: Available online 18 September 2018Source: Applied Numerical MathematicsAuthor(s): Yongtao Zhou, Chengjian ZhangAbstractIn this paper, by combining the p-order block boundary value methods with the m-th Lagrange interpolation, a class of new numerical methods for solving nonlinear fractional differential equations with the γ-order (0
  • Error estimates of pseudostress-velocity MFEM for optimal control problems
           governed by stokes equations
    • Abstract: Publication date: Available online 18 September 2018Source: Applied Numerical MathematicsAuthor(s): Yanping Chen, Haitao Leng, Wendi YangAbstractIn this paper, we study the error analysis of R-T mixed finite element methods for optimal control problems governed by stationary Stokes equations. To avoid the difficulty induced by the symmetry constraint on the stress tensor, we use the pseudostress proposed by Arnold and Falk [11] to replace it. Then, we obtain the pseudostress-velocity formulation of the optimal control problem. Moreover, we prove a priori and a posteriori error estimates for the new formulation, and an example is provided to confirm our results.
  • A novel kind of efficient symplectic scheme for
           Klein-Gordon-Schrödinger equation
    • Abstract: Publication date: Available online 17 September 2018Source: Applied Numerical MathematicsAuthor(s): Linghua Kong, Meng Chen, Xiuling YinAbstractIn this paper, we construct a family of high order compact symplectic (S-HOC) schemes for the Klein-Gordon-Schrödinger (KGS) equation. The KGS can be cast into a Hamiltonian form. At first, we discretize the Hamiltonian system in space by a high order compact method which has higher convergent rate than general finite difference methods. Then the semi-discretized system is approximated in time by the Euler midpoint scheme which preserves the symplectic structure of the original system. The conserved quantities of the scheme, including symplectic structure conservation law, charge conservation law and energy conservation law, are discussed. The local truncation error and global error of the numerical solvers are investigated. Finally, some numerical verifications are presented to numerically validate the theoretical analysis. The numerical results are persuasive and illustrate the theoretical analysis.
  • Anisotropic bivariate subdivision with applications to multigrid
    • Abstract: Publication date: Available online 13 September 2018Source: Applied Numerical MathematicsAuthor(s): M. Charina, M. Donatelli, L. Romani, V. TuratiAbstractIn this paper, motivated by applications to multigrid, we present families of anisotropic, bivariate, interpolating and approximating subdivision schemes. We study the minimality and polynomial generation/reproduction properties of both families and Hölder regularity of their prominent representatives. From the symbols of the proposed subdivision schemes, we define bivariate grid transfer operators for anisotropic multigrid methods. We link the generation/reproduction properties of subdivision to the convergence and optimality of the corresponding multigrid methods. We illustrate the performance of our subdivision based grid transfer operators on examples of anisotropic Laplacian and biharmonic problems.
  • The Method of Fundamental Solutions for the Helmholtz Equation
    • Abstract: Publication date: Available online 13 September 2018Source: Applied Numerical MathematicsAuthor(s): Zi-Cai Li, Yimin Wei, Yunkun Chen, Hung-Tsai HuangAbstractIn this paper, we study the Helmholtz equation by the method of fundamental solutions (MFS) using Bessel and Neumann functions. The bounds of errors are derived for bounded simply-connected domains, while the bounds of condition number are derived only for disk domains. The MFS using Bessel functions is more efficient than the MFS using Neumann functions. Note that by using Bessel functions, the radius R of the source nodes is not necessarily to be larger than the maximal radius rmax of the solution domain. This is against the well-known rule: rmax
  • A regularizing multilevel approach for nonlinear inverse problems
    • Abstract: Publication date: Available online 12 September 2018Source: Applied Numerical MathematicsAuthor(s): Min Zhong, Wei WangAbstractIn this paper, we propose a multilevel method for solving nonlinear inverse problems F(x)=y in Banach spaces. By minimizing the discretized version of the regularized functionals at different levels, we define a sequence of regularized approximations to the sought solution, which is shown to be stable and globally convergent. The penalty term Θ in regularized functionals is allowed to be non-smooth to include Lp−L1 or Lp−TV (Total Variation) reconstructions, which are significant in reconstructing special features of solutions such as sparsity and discontinuities. Two parameter identification examples are presented to validate the theoretical analysis and to verify the effectiveness of the method.
  • Epsilon penalty method combined with an extension of the Ritz method for
           solving a class of fractional optimal control problems with mixed
           inequality constraints
    • Abstract: Publication date: Available online 11 September 2018Source: Applied Numerical MathematicsAuthor(s): Ali LotfiAbstractIn this paper an approximate scheme based on the combination of epsilon penalty method and an extended form of the well-known Ritz method is developed for solving a general class of fractional optimal control problems. The fractional derivative is in the Caputo sense, and the given fractional optimal control problem contains both state and control inequality constraints. The convergence of the method is discussed extensively and its effectiveness is verified by some illustrative test examples.
  • A mixed virtual element method for a pseudostress-based formulation of
           linear elasticity
    • Abstract: Publication date: Available online 10 September 2018Source: Applied Numerical MathematicsAuthor(s): Ernesto Cáceres, Gabriel N. Gatica, Filánder A. SequeiraAbstractIn this paper we introduce and analyze a mixed virtual element method (mixed-VEM) for a pseudostress-displacement formulation of the linear elasticity problem with non-homogeneous Dirichlet boundary conditions. We follow a previous work by some of the authors, and employ a mixed formulation that does not require symmetric tensor spaces in the finite element discretization. More precisely, the main unknowns here are given by the pseudostress and the displacement, whereas other physical quantities such as the stress, the strain tensor of small deformations, and the rotation, are computed through simple postprocessing formulae in terms of the pseudostress variable. We first recall the corresponding variational formulation, and then summarize the main mixed-VEM ingredients that are required for our discrete analysis. In particular, we utilize a well-known local projector onto a suitable polynomial subspace to define a calculable version of our discrete bilinear form, whose continuous version requires information of the variables on the interior of each element. Next, we show that the global discrete bilinear form satisfies the hypotheses required by the Babusˇka-Brezzi theory. In this way, we conclude the well-posedness of our mixed-VEM scheme and derive the associated a priori error estimates for the virtual solutions as well as for the fully computable projections of them. Furthermore, we also introduce a second element-by-element postprocessing formula for the pseudostress, which yields an optimally convergent approximation of this unknown with respect to the broken H(div)-norm. In addition, this postprocessing formula can also be applied to the postprocessed stress tensor. Finally, several numerical results illustrating the good performance of the method and confirming the theoretical rates of convergence are presented.
  • A reliable, efficient and localized error estimator for a discontinuous
           Galerkin method for the Signorini problem
    • Abstract: Publication date: Available online 5 September 2018Source: Applied Numerical MathematicsAuthor(s): Mirjam WallothAbstractWe present a new residual-type a posteriori error estimator for the discontinuous finite element solution of contact problems. The theoretical results are derived for two and three-dimensional domains and arbitrary gap functions. The estimator yields upper and lower bounds to a suitable error norm which measures the error in the displacements and in a quantity related to the contact stresses and the actual contact zone. In the derivation of the error estimator the local properties of the discontinuous solution are exploited appropriately so that, on the one hand, the error estimator has no contributions related to the non-linearity in the interior of the actual contact zone and, on the other hand, the critical region between the actual and non-actual contact zone can be well refined.
  • Recovering the time-dependent potential function in a multi-term
           time-fractional diffusion equation
    • Abstract: Publication date: Available online 5 September 2018Source: Applied Numerical MathematicsAuthor(s): Liangliang Sun, Yun Zhang, Ting WeiAbstractIn the present paper, we devote our effort to a nonlinear inverse problem for recovering a time-dependent potential term in a multi-term time-fractional diffusion equation from the boundary measured data. First we study the existence, uniqueness and regularity of solution for the direct problem by using the fixed point theorem. Then a stability estimate of inverse coefficient problem is obtained based on the regularity of solution of direct problem and some generalized Gronwall's inequalities. Numerically, we reformulate the inverse potential function into a variational problem, and we use a Levenberg-Marquardt method to find the approximate potential function. Numerical experiments for five examples in one-dimensional and two-dimensional cases are provided to show the effectiveness of the proposed method.
  • Unified smoothing functions for absolute value equation associated with
           second-order cone
    • Abstract: Publication date: Available online 3 September 2018Source: Applied Numerical MathematicsAuthor(s): Chieu Thanh Nguyen, B. Saheya, Yu-Lin Chang, Jein-Shan ChenAbstractIn this paper, we explore a unified way to construct smoothing functions for solving the absolute value equation associated with second-order cone (SOCAVE). Numerical comparisons are presented, which illustrate what kinds of smoothing functions work well along with the smoothing Newton algorithm. In particular, the numerical experiments show that the well known loss function widely used in engineering community is the worst one among the constructed smoothing functions, which indicates that the other proposed smoothing functions can be employed for solving engineering problems.
  • A posteriori error analysis for the scattering by obstacles in a
           homogeneous chiral environment
    • Abstract: Publication date: Available online 3 September 2018Source: Applied Numerical MathematicsAuthor(s): Zhoufeng Wang, Rui GuoAbstractIn this paper we consider the scattering of time-harmonic electromagnetic wave propagation in a homogeneous chiral environment by obstacles. The model is simplified to a two-dimensional scattering problem, and is formulated as a boundary value problem in a bounded domain by introducing the nonlocal boundary conditions associated with Dirichlet-to-Neumann operators. An a posteriori error estimate is established when the truncation of the nonlocal boundary operators takes place. The crucial part of the error analysis is to develop a duality argument and use the Bohren decomposition of the electromagnetic fields. The a posteriori error estimate consists of two parts, finite element approximation error and the truncation error of boundary operators which decays exponentially with respect to the truncation parameter. Numerical experiments are also presented to show the robustness and effectiveness of our numerical algorithm.
  • An efficient two-level preconditioner for multi-frequency wave propagation
    • Abstract: Publication date: Available online 3 September 2018Source: Applied Numerical MathematicsAuthor(s): Manuel Baumann, Martin B. van GijzenAbstractWe consider wave propagation problems that are modeled in the frequency-domain, and that need to be solved simultaneously for multiple frequencies within a fixed range. For this, a single shift-and-invert preconditioner at a so-called seed frequency is applied. The choice of the seed is crucial for the performance of preconditioned multi-shift GMRES and is closely related to the parameter choice for the Complex Shifted Laplace preconditioner. Based on a classical GMRES convergence bound, we present an analytic formula for the optimal seed parameter that purely depends on the original frequency range. The new insight is exploited in a two-level preconditioning strategy: A shifted Neumann preconditioner with minimized spectral radius is additionally applied to multi-shift GMRES. Moreover, we present a reformulation of the multi-shift problem to a matrix equation solved with, for instance, global GMRES. Here, our analysis allows for rotation of the spectrum of the linear operator. Numerical experiments for the time-harmonic visco-elastic wave equation demonstrate the performance of the new preconditioners.
  • Revisiting the Thomas-Fermi Equation: Accelerating Rational Chebyshev
           Series Through Coordinate Transformations
    • Abstract: Publication date: Available online 29 August 2018Source: Applied Numerical MathematicsAuthor(s): Xiaolong Zhang, John P. BoydAbstractWe revisit the spectral solution of the Thomas-Fermi problem for neutral atoms, urr−(1/r)u3/2=0 on r∈[0,∞] with u(0)=1 and u(∞)=0 to illustrate some themes in solving differential equations when there are complications, and also to make improvements in our earlier treatment. By “complications” we mean features of the problem that either destroy the exponential accuracy of a standard Chebyshev series, or render the classic Chebyshev approach inapplicable. The Thomas-Fermi problem has four complications: (i) a semi-infinite domain r∈[0,∞] (ii) a square root singularity in u(r) at the origin (iii) a fractional power nonlinearity and (iv) asymptotic decay as r→∞ that includes negative powers of r with fractional exponents. Our earlier treatment determined the slope at the origin to twenty-five decimal places, but no fewer than 600 basis functions were required to approximate a univariate solution that is everywhere monotonic, and all of the earlier tricks failed to recover an exponential rate of convergence in the truncation of the spectral series N, but only a high order convergence in negative powers of N. Here, using the coordinate z≡r to neutralize the square root singularity as before, we show that accuracy and rate of convergence are significantly improved by solving for the original unknown u(r) instead of the modified unknown v(r)=u(r) used previously. Without a further change of coordinate, a rational Chebyshev basis TLn(z;L) yields twelve decimal place for the slope at the origin, ur(0), with 70 basis functions and twenty-four places with a truncation N=100.True exponential accuracy can be restored by using an appropriate change of coordinates, z=G(Z), where G is some species of exponential. However, the various Chebyshev and Fourier series for the Thomas-Fermi function have “plural asymptotics”, that is, an∼aintermediate(n) for 1≪nͱ...
  • On the optimal CFL number of SSP methods for hyperbolic problems
    • Abstract: Publication date: Available online 28 August 2018Source: Applied Numerical MathematicsAuthor(s): Andrew Giuliani, Lilia KrivodonovaAbstractWe show that the theory for strong stability preserving (SSP) time stepping methods employed with the method of lines-type discretizations of hyperbolic conservation laws may result in overly stringent time step restrictions. We analyze a fully discrete finite volume method with slope reconstruction and a second order SSP Runge-Kutta time integrator to show that the maximum stable time step can be increased over the SSP limit. Numerical examples indicate that this result extends to two-dimensional problems on triangular meshes.
  • Zeros of quadratic quasi-orthogonal order 2 polynomials
    • Abstract: Publication date: Available online 27 August 2018Source: Applied Numerical MathematicsAuthor(s): Claude Brezinski, K.A. Driver, Michela Redivo-ZagliaAbstractCorollary 2 in [1] states that for −32
  • A transversal method of lines for the numerical modeling of vertical
           infiltration into the vadose zone
    • Abstract: Publication date: Available online 25 August 2018Source: Applied Numerical MathematicsAuthor(s): Marco Berardi, Fabio Difonzo, Filippo Notarnicola, Michele VurroAbstractHere some issues are studied, related to the numerical solution of Richards' equation in a one dimensional spatial domain by a technique based on the Transversal Method of Lines (TMoL). The core idea of TMoL approach is to semi-discretize the time derivative of Richards' equation: afterward a system of second order differential equations in the space variable is derived as an initial value problem.The computational framework of this method requires both Dirichlet and Neumann boundary conditions at the top of the column. The practical motivation for choosing such a condition is argued. We will show that, with the choice of the aforementioned initial conditions, our TMoL approach brings to solutions comparable with the ones obtained by the classical Methods of Lines (hereafter referred to as MoL) with corresponding standard boundary conditions: in particular, an appropriate norm is introduced for effectively comparing numerical tests obtained by MoL and TMoL approach and a sensitivity analysis between the two methods is performed by means of a mass balance point of view. A further algorithm is introduced for deducing in a self sustaining way the gradient boundary condition on top in the TMoL context.
  • Numerical analysis for a nonlinear model of elastic strings with moving
    • Abstract: Publication date: Available online 25 August 2018Source: Applied Numerical MathematicsAuthor(s): M.A. Rincon, I.-S. Liu, W.R. Huarcaya, B.A. CarmoAbstractIn this article, the error estimates for semi-discrete and totally discrete problems of a nonlinear model of elastic strings with moving boundary are established. We consider an extension of the Kirchhoff model, that takes into account the change of length of the string during vibration. The existence and uniqueness theorems of the problem, already known in the literature, will be stated for reference. The present numerical analysis is based on finite element method in spatial variable and finite difference method in time with the Newmark's approximation. Since the problem is nonlinear the resulting algebraic system is nonlinear to be solved by Newton's method.Numerical examples are presented for different kinds of moving boundary to verify the efficiency and feasibility of the method and check the coherence with the theoretical analysis. From the numerical results, the rate of convergence are shown to be consistent with the order of convergence expected from the theoretical ones.
  • Iterative methods for the delay Lyapunov equation with T-Sylvester
    • Abstract: Publication date: Available online 23 August 2018Source: Applied Numerical MathematicsAuthor(s): Elias Jarlebring, Federico PoloniAbstractThe delay Lyapunov equation is an important matrix boundary-value problem which arises as an analogue of the Lyapunov equation in the study of time-delay systems x˙(t)=A0x(t)+A1x(t−τ)+B0u(t). We propose a new algorithm for the solution of the delay Lyapunov equation. Our method is based on the fact that the delay Lyapunov equation can be expressed as a linear system of equations, whose unknown is the value U(τ/2)∈Rn×n, i.e., the delay Lyapunov matrix at time τ/2. This linear matrix equation with n2 unknowns is solved by adapting a preconditioned iterative method such as GMRES. The action of the n2×n2 matrix associated to this linear system can be computed by solving a coupled matrix initial-value problem. A preconditioner for the iterative method is proposed based on solving a T-Sylvester equation MX+XTN=C, for which there are methods available in the literature. We prove that the preconditioner is effective under certain assumptions. The efficiency of the approach is illustrated by applying it to a time-delay system stemming from the discretization of a partial differential equation with delay. Approximate solutions to this problem can be obtained for problems of size up to n≈1000, i.e., a linear system with n2≈106 unknowns, a dimension which is outside of the capabilities of the other existing methods for the delay Lyapunov equation.
  • Solving the backward problem in Riesz–Feller fractional diffusion by a
           new nonlocal regularization method ⋆
    • Abstract: Publication date: Available online 21 August 2018Source: Applied Numerical MathematicsAuthor(s): Guang-Hui ZhengAbstractIn this paper, we consider the backward problem introduced in [48] for Riesz–Feller fractional diffusion. To begin with, some basic properties of solution of the corresponding forward problem, such as the Lp estimates, symmetry property and asymptotic estimates, are established by Fourier analysis technique. And then, under various a priori bound assumptions, we give the L2 conditional stability estimates for the solution of backward problem and also its symmetry property. Moreover, in order to overcome the ill-posedness of the backward problem, we propose a new nonlocal regularization method (NLRM) to solve it. That is, the following nonlocal variational functional is introducedJ(φ)=12‖u(φ(x);x,T)−fδ(x)‖2+β2‖[Pα1⁎φ](x)‖2, where β∈(0,1) is a regularization parameter, “⁎” denotes the convolution operation and Pα1(x) is called a convolution kernel with parameter α1, which will be selected properly later. The minimizer of above variational problem is defined as the regularization solution, and the L2 estimates, symmetry property of regularization solution are given. These results actually show the well-posedness of nonlocal variational problem. Our idea is essentially that using this well-posed problem to approximate the backward (ill-posed) problem. Thus, under an a posteriori parameter choice rule, we deduce various convergence rate estimates under different a-priori bound assumptions for the exact solution. Finally, several numerical examples are given to show that the proposed numerical methods are effective and adaptive for different a-priori information.
  • A Posteriori analysis of an IMEX Entropy-Viscosity formulation for
           Hyperbolic conservation laws with dissipation
    • Abstract: Publication date: Available online 20 August 2018Source: Applied Numerical MathematicsAuthor(s): Jehanzeb H. Chaudhry, John N. Shadid, Timothy WildeyAbstractThis study considers adjoint based a posteriori estimation of the error in a quantity of interest computed from numerical solutions based on multistage implicit-explicit (IMEX) time integration schemes and an entropy-viscosity formulation for damped hyperbolic partial differential equations (PDEs). Hyperbolic systems are challenging to solve numerically due to the need to stabilize systems with discontinuous or nearly discontinuous solutions in an attempt to limit non-physical oscillations and provide accurate solutions. The goal of this effort is to provide error estimates based on adjoint operators, variational analysis and computable residuals. The error estimates quantify the total error as well as different contributions to the error arising from the time integration schemes and choices of numerical parameters in the numerical method.
  • Mathematical modeling and numerical simulation of the TGF-β/Smad
           signaling pathway in tumor microenvironments
    • Abstract: Publication date: November 2018Source: Applied Numerical Mathematics, Volume 133Author(s): Adnan Morshed, Prashanta Dutta, Robert H. DillonAbstractMammalian cells respond in a variety of ways to concentrations of activated transforming growth factor (TGF) in the extracellular domain via intracellular Smad signaling pathways. TGF-β/Smad interaction is prevalent in a wide range of tumor environments with both autocrine and paracrine mechanisms driving temporal evolution. TGF-β has been reported to exhibit both pro- and antagonistic roles in tumor progression and survival. It has been argued that tumor cells upregulate TGF-β production by modifying the TGF-β/Smad pathway. Although several numerical studies of the tumor microenvironment have been conducted, most are based on PDEs where the cells are represented as a continuum or on discrete agent-based methods. Here, we develop a hybrid/cells-based model for the tumor microenvironment with cells represented as discrete entities in which diffusion and reaction in the extracellular environment as well as protein/receptor surface interactions are described using immersed interface methods. We develop a model for the cellular TGF-β/Smad pathways with the intracellular processes represented by systems of ODE's. In this study, we investigate the variation in cellular response with different levels of TGF-β in the extracellular environment. Our numerical results show that the model can account for experimental results on tumor cell lines when intracellular production and secretion of TGF-B based on activated Smad concentrations are included in the model system. Additionally, we studied the intracellular and surface levels of TGF-β in two different types of cells and observed variation in the extracellular domain.
  • An hp-version error analysis of the discontinuous Galerkin method
           for linear elasticity
    • Abstract: Publication date: November 2018Source: Applied Numerical Mathematics, Volume 133Author(s): Jianguo Huang, Xuehai HuangAbstractAn hp-version error analysis is developed for the general DG method in mixed formulation for solving the linear elastic problem. First of all, we give the hp-version error estimates of two L2 projection operators. Then incorporated with the techniques in [11], we obtain the hp-version error estimates in energy norm and L2 norm. Some numerical experiments are provided for demonstrating the theoretical results.
  • The upwind hybrid difference methods for a convection diffusion equation
    • Abstract: Publication date: November 2018Source: Applied Numerical Mathematics, Volume 133Author(s): Youngmok Jeon, Mai Lan TranAbstractWe propose the upwind hybrid difference method and its penalized version for the convection dominated diffusion equation. The hybrid difference method is composed of two types of approximations: one is the finite difference approximation of PDEs within cells (cell FD) and the other is the interface finite difference (interface FD) on edges of cells. The interface finite difference is derived from continuity of normal fluxes. The penalty method is obtained by adding small diffusion in the interface FD. The penalty term makes it possible to reduce severe numerical oscillations in the upwind hybrid difference solutions. The penalty parameter is designed to be some power of the grid size. A complete stability is provided. Convergence estimates seems to be conservative according to our numerical experiments. To exposit convergence property and controllability of numerical oscillations several numerical tests are provided.
  • An accurate and asymptotically compatible collocation scheme for nonlocal
           diffusion problems
    • Abstract: Publication date: November 2018Source: Applied Numerical Mathematics, Volume 133Author(s): Xiaoping Zhang, Jiming Wu, Lili JuAbstractIn this paper, we develop and analyze a collocation scheme for solving the linear nonlocal diffusion problem with general kernels. To approximate the nonlocal diffusion operator, we take a classic trapezoidal rule based on the linear interpolation as the starting point, and then carefully derive a new improved quadrature rule, which is not only more accurate but also could avoid the evaluations of singular integrals. We then use this rule to construct a collocation scheme for solving the nonlocal diffusion equations, that produces a symmetric positive definite stiffness matrix with Toeplitz structure. The proposed scheme is rigorously shown to be of second order accurate with respect to the mesh size for the nonlocal problem with fixed horizon, and in particular, it can achieve higher order accuracy for the commonly used kernels in the literature. Furthermore, we also prove that the scheme is asymptotically compatible, i.e., the approximate solution of the nonlocal diffusion problem converges to the exact solution of the corresponding local PDE problem when the horizon and the mesh size both go to zero. Finally, numerical experiments are presented to verify the theoretical results.
  • Mixed generalized Hermite–Fourier spectral method for Fokker–Planck
           equation of periodic field
    • Abstract: Publication date: November 2018Source: Applied Numerical Mathematics, Volume 133Author(s): Guo Chai, Tian-jun WangAbstractIn this paper, we develop a mixed Hermite–Fourier spectral method for the Fokker–Planck equation of periodic field. Radical to the numerical solutions of partial differential equations of the sort, results on mixed generalized Hermite–Fourier orthogonal approximation are established. The convergence of the constructed spectral scheme is proved. Numerical results show the efficiency of this approach and coincide well with theoretical analysis.
  • A stabilized finite volume element method for a coupled Stokes–Darcy
    • Abstract: Publication date: November 2018Source: Applied Numerical Mathematics, Volume 133Author(s): Rui Li, Jian Li, Xiaoming He, Zhangxin ChenAbstractIn this paper, we present a stabilized finite volume element method with the conforming finite element triples P1–P0–P1 and P1–P1–P1 for approximating the velocity, pressure, and hydraulic head of a coupled Stokes–Darcy problem. The proposed method is convenient to implement, computationally efficient, mass conserving, optimally accurate, and able to handle complex geometries; therefore, this method has great potential to be useful for realistic problems involving coupled free flow and porous media flow. To offset the lack of the inf-sup condition of the P1–P0 and P1–P1 elements for the Stokes equation, a parameter free stabilization term is added to the discrete formulation. Stability and optimal error estimates are proved based on a bridge built up between the finite volume element method and the finite element method. An element level implementation of the stabilization term is discussed so that an existing code package can be conveniently modified to handle the stabilization procedures. A series of numerical experiments are provided to illustrate the above features of the proposed method, the theoretical results, and the realistic applications.
  • 20th IMACS World Congress
    • Abstract: Publication date: November 2018Source: Applied Numerical Mathematics, Volume 133Author(s): Xiaoming He, Jie Shen, Chuanju Xu
  • Numerical simulations of the spread of rabies in two-dimensional space
    • Abstract: Publication date: Available online 17 August 2018Source: Applied Numerical MathematicsAuthor(s): Khalaf M. Alanazi, Zdzislaw Jackiewicz, Horst R. ThiemeAbstractWe extend our previous work on the spatial spread of fox rabies from one dimension to two dimensions. We consider the case when the latent period has fixed length. We use the method of lines to replace the spatial derivatives and the integral equations with algebraic approximations, then we apply the explicit continuous Runge-Kutta method of fourth order and discrete Runge-Kutta method of third order with six stages to numerically integrate the resulting systems of ordinary and delay differential equations. We discuss and confirm some of the major results we obtained in earlier work. The asymptotic speeds of spread observed in the two-dimensional simulations and in earlier work are discussed and compared with those found in nature.
  • A class of stochastic one-parameter methods for nonlinear SFDEs with
           piecewise continuous arguments
    • Abstract: Publication date: Available online 15 August 2018Source: Applied Numerical MathematicsAuthor(s): Ying Xie, Chengjian ZhangAbstractThis paper deals with nonlinear stochastic functional differential equations with piecewise continuous arguments (SFDEPCAs). Based on an adaptation to the underlying one-leg θ-methods for ODEs, a class of new one-parameter methods for nonlinear SFDEPCAs are introduced. The mean-square exponential stability criteria of analytical and numerical solutions are derived. Under the suitable conditions, it is proved that the one-parameter methods are convergent with strong order 1/2. Some numerical experiments are given to illustrate the theoretical results and computational advantages of the induced methods.
  • A direct discontinuous Galerkin method for a time-fractional diffusion
           equation with a Robin boundary condition
    • Abstract: Publication date: Available online 15 August 2018Source: Applied Numerical MathematicsAuthor(s): Chaobao Huang, Martin StynesAbstractA time-fractional reaction-diffusion initial-boundary value problem with Robin boundary condition is considered on the domain Ω×[0,T], where Ω=(0,l)⊂R. The coefficient of the zero-order reaction term is not required to be non-negative, which complicates the analysis. In general the unknown solution will have a weak singularity at the initial time t=0. Existence and uniqueness of the solution and pointwise bounds on some of its derivatives are derived. A fully discrete numerical method for computing an approximate solution is investigated; it uses the well-known L1 discretization on a graded mesh in time and a direct discontinuous Galerkin (DDG) finite element method on a uniform mesh in space. Discrete stability of the computed solution is proved. Its error is bounded in the L2(Ω) and H1(Ω) norms at each discrete time level tn by means of a non-trivial projection of the unknown solution into the finite element space. The L2(Ω) bound is optimal for all tn; the H1(Ω) bound is optimal for tn not close to t=0. An optimal grading of the temporal mesh can be deduced from these bounds. Numerical results show that our analysis is sharp.
  • A high-order numerical algorithm for two-dimensional time-space tempered
           fractional diffusion-wave equation
    • Abstract: Publication date: Available online 13 August 2018Source: Applied Numerical MathematicsAuthor(s): Hengfei DingAbstractIn this paper, we focus our attention on the development of the high-order numerical algorithm for the time-space tempered fractional diffusion-wave equation in two spatial dimensions. Based on the fourth-order fractional-compact difference operator, a new difference scheme with convergence order O(τ2+h14+h24) is derived, where τ is the temporal stepsize, h1 and h2 are the spatial stepsizes, respectively. The stability and convergence of the algorithm are investigated by the energy method and numerical experiment is carried out to verify the feasibility of the numerical algorithm.
  • A high order discontinuous Galerkin method with Lagrange multipliers for
           second-order elliptic problems
    • Abstract: Publication date: Available online 11 August 2018Source: Applied Numerical MathematicsAuthor(s): Mi-Young Kim, Dong-wook ShinAbstractA discontinuous Galerkin method with Lagrange multipliers (DGLM) is developed to approximate the solution to the second-order elliptic problems. Lagrange multipliers for the solution and for the flux are considered on the edge/face of each element. The weak gradient and the weak divergence are defined for the elliptic problems. Lagrange multipliers for the solution and for the flux are shown to be the averages of the solutions and the “normal” fluxes at the edge/face, respectively. Unique solvability of the discrete system is proved and an error estimate is derived. The element unknowns are solved in terms of the Lagrange multipliers in element by element fashion. The Schur complement system of the Lagrange multipliers has a block structure, which is kept unchanged while the inside of the blocks gets dense in the higher order approximation. An explanation on algorithmic aspects is given. Some numerical results are presented.
  • Modified Chebyshev collocation method for pantograph-type differential
    • Abstract: Publication date: Available online 9 August 2018Source: Applied Numerical MathematicsAuthor(s): Changqing YangAbstractThe Chebyshev polynomials and a collocation method are applied to the solution of the pantograph equation. A Chebyshev pantograph operational matrix is derived and used to reduce the pantograph equation to a system of algebraic equations. The convergence order of the proposed method is investigated in the L2-norm. Numerical examples are presented to verify the efficiency and accuracy of the proposed method. Results reveal that this method is accurate and easy to implement.
  • Verification of a Low-degree Polynomial Vanishing at Empirical Points
    • Abstract: Publication date: Available online 8 August 2018Source: Applied Numerical MathematicsAuthor(s): Zhe Li, Kai Zheng, Shugong ZhangAbstractGiven a set of distinct empirical points with uniform tolerance, based on the LDP algorithm proposed by Fassino and Torrente, we provide a verification algorithm that computes a polynomial, an admissible perturbed point set with verified error bound, such that the polynomial is guaranteed to vanish at a slightly admissible perturbed point set within computed error bound. The effectiveness of our algorithm is demonstrated in several examples.
  • Comparison results for splitting iterations for solving multi-linear
    • Abstract: Publication date: Available online 24 July 2018Source: Applied Numerical MathematicsAuthor(s): Wen Li, Dongdong Liu, Seak-Weng VongAbstractIt is known that the spectral radius of the iterative tensor can be seen as an approximate convergence rate for solving multi-linear systems by tensor splitting iterative methods. So in this paper, first we give some spectral radius comparisons between two different iterative tensors. Then, we propose the preconditioned tensor splitting method for solving multi-linear systems, which provides an alternative algorithm with the choice of a preconditioner. In particular, also we give some spectral radius comparisons between the preconditioned iterative tensor and the original one. Numerical examples are given to demonstrate the efficiency of the proposed preconditioned methods.
  • Homotopy Perturbation Method for the nonlinear MHD Jeffery–Hamel
           blood flows problem
    • Abstract: Publication date: Available online 24 July 2018Source: Applied Numerical MathematicsAuthor(s): Iftikhar Ahmad, Hira IlyasAbstractIn this paper, Homotopy Perturbation Method is applied to solve the nonlinear MHD Jeffery–Hamel arterial blood flow problem. Primarily, two-dimensional nonlinear Navier–Stokes equations have been converted into third order one-dimensional equation by means of transformation rule. Later the solution of governed equation is obtained by using Homotopy Perturbation Method. The proposed numerical results show a good agreement with reference solution for finite interval and emphasize to understand the human arterial blood flow rate. Further, accuracy and reliability of the proposed method is checked by increasing the iteration process up to third order. Finally, the results showed that product of angle between plates “α” and Reynolds number “Re” is directly proportional to the MHD Jeffery–Hamel flow.
  • An introduction to multi-trace formulations and associated domain
           decomposition solvers
    • Abstract: Publication date: Available online 24 July 2018Source: Applied Numerical MathematicsAuthor(s): X. Claeys, V. Dolean, M.J. GanderAbstractMulti-trace formulations (MTFs) are based on a decomposition of the problem domain into subdomains, and thus domain decomposition solvers are of interest. The fully rigorous mathematical MTF can however be daunting for the non-specialist. The first aim of the present contribution is to provide a gentle introduction to MTFs. We introduce these formulations on a simple model problem using concepts familiar to researchers in domain decomposition. This allows us to get a new understanding of MTFs and a natural block Jacobi iteration, for which we determine optimal relaxation parameters. We then show how iterative multi-trace formulation solvers are related to a well known domain decomposition method called optimal Schwarz method: a method which used Dirichlet to Neumann maps in the transmission condition. We finally show that the insight gained from the simple model problem leads to remarkable identities for Calderón projectors and related operators, and the convergence results and optimal choice of the relaxation parameter we obtained is independent of the geometry, the space dimension of the problem, and the precise form of the spatial elliptic operator, like for optimal Schwarz methods. We illustrate our analysis with numerical experiments.
  • Diamond shaped standing wave patterns of a two-dimensional Boussinesq
    • Abstract: Publication date: Available online 23 July 2018Source: Applied Numerical MathematicsAuthor(s): Shenghao Li, Min ChenAbstractWe prove the existence of a large family of two-dimensional diamond shaped standing waves for a Boussinesq system which describes two-way propagation of water waves in a channel. Our proof uses the Lyapunov-Schmidt method to find the bifurcation standing waves.
  • Localized kernel-based approximation for pricing financial options under
           regime switching jump diffusion model
    • Abstract: Publication date: Available online 18 July 2018Source: Applied Numerical MathematicsAuthor(s): Reza Mollapourasl, Majid Haghi, Ruihua LiuAbstractIn this paper, we consider European and American option pricing problems under regime switching jump diffusion models which are formulated as a system of partial integro-differential equations (PIDEs) with fixed and free boundaries. For free boundary problem arising in pricing American option, we use operator splitting method to deal with early exercise feature of American option. For developing a numerical technique we employ localized radial basis function generated finite difference (RBF-FD) approximation to overcome the ill-conditioning and high density issues of discretized matrices. The proposed method leads to linear systems with tridiagonal and diagonal dominant matrices. Also, in this paper the convergence and consistency of the proposed method are discussed. Numerical examples presented in the last section illustrate the robustness and practical performance of the proposed algorithm for pricing European and American options.
  • A second-order, uniquely solvable, energy stable BDF numerical scheme for
           the phase field crystal model
    • Abstract: Publication date: Available online 18 July 2018Source: Applied Numerical MathematicsAuthor(s): Qi Li, Liquan Mei, Bo YouAbstractIn this paper, we propose a second-order time accurate convex splitting scheme for the phase field crystal model. The temporal discretization is based on the second-order backward differentiation formula (BDF) and a convex splitting of the energy functional. The mass conservation, unconditionally unique solvability, unconditionally energy stability and convergence of the numerical scheme are proved rigorously. Mixed finite element method is employed to obtain the fully discrete scheme due to a sixth-order spatial derivative. Numerical experiments are presented to demonstrate the accuracy, mass conservation, energy stability and effectiveness of the proposed scheme.
  • A fast numerical method for two-dimensional Riesz space fractional
           diffusion equations on a convex bounded region
    • Abstract: Publication date: Available online 17 July 2018Source: Applied Numerical MathematicsAuthor(s): S. Chen, F. Liu, I. Turner, V. AnhAbstractFractional differential equations have attracted considerable attention because of their many applications in physics, geology, biology, chemistry, and finance. In this paper, a two-dimensional Riesz space fractional diffusion equation on a convex bounded region (2D-RSFDE-CBR) is considered. These regions are more general than rectangle or circular domains. A novel alternating direction implicit method for the 2D-RSFDE-CBR with homogeneous Dirichlet boundary conditions is proposed. The stability and convergence of the method are discussed. The resulting linear systems are Toeplitz-like and are solved by the preconditioned conjugate gradient method with a suitable circulant preconditioner. By the fast Fourier transform, the method only requires a computational cost of O(nlog⁡n) per time step. These numerical techniques are used for simulating a two-dimensional Riesz space fractional FitzHugh-Nagumo model. The numerical results demonstrate the effectiveness of the method. These techniques can be extended to three spatial dimensions, which will be the topic of our future research.
  • A novel finite difference discrete scheme for the time fractional
           diffusion-wave equation
    • Abstract: Publication date: Available online 11 July 2018Source: Applied Numerical MathematicsAuthor(s): Zhengguang Liu, Aijie Cheng, Xiaoli LiAbstractIn this article, we consider initial and boundary value problems for the diffusion-wave equation involving a Caputo fractional derivative(of order α, with 1
  • Numerical analysis of a collocation method for functional integral
    • Abstract: Publication date: Available online 10 July 2018Source: Applied Numerical MathematicsAuthor(s): Adson M. Rocha, Juarez S. Azevedo, Saulo P. Oliveira, Maicon R. CorreaAbstractWe study the numerical approximation of functional integral equations, a class of nonlinear Fredholm-type integral equations of the second kind, by the collocation method with piecewise continuous basis functions. The resulting nonlinear algebraic system is solved with the Picard iteration method. Starting from the analysis of the continuous problem in L∞([a,b]), we prove the convergence of numerical solution and, under an additional regularity assumption, provide an a priori error estimate. Numerical examples illustrate the predicted theoretical results.
  • Estimates of the discrete van Cittert deconvolution error in approximate
           deconvolution models of turbulence in bounded domains
    • Abstract: Publication date: Available online 4 July 2018Source: Applied Numerical MathematicsAuthor(s): Argus A. DuncaAbstractRecent turbulence models such as the Approximate Deconvolution Model (ADM) or the Leray-deconvolution model are derived from the Navier–Stokes equations using the van Cittert approximate deconvolution method. As a consequence, the numerical error in the approximation of the Navier–Stokes weak solution with discrete solutions of the above models is influenced also by the discrete deconvolution error u−DNu‾h caused by the approximate deconvolution method. Here u is the flow field, u‾ is its average and DN is the N-th order van Cittert deconvolution operator.It is therefore important to analyze the deconvolution error u−DNu‾h in terms of the mesh size h, the filter radius α and the order N of the deconvolution operators used in the computation.This problem is investigated herein in the case of bounded domains and zero-Dirichlet boundary conditions. It is proved that on a sequence of quasiuniform meshes the L2 norm of the discrete deconvolution error convergences to 0 in the order of hk+1+KNh provided that the filter radius is in the order of the mesh size and the flow field has enough regularity. Here K
  • Fourth order wave equation with nonlinear strain and logarithmic
    • Abstract: Publication date: Available online 27 June 2018Source: Applied Numerical MathematicsAuthor(s): Runzhang Xu, Wei Lian, Xiangkun Kong, Yanbing YangAbstractThe main goal of this work is to investigate the initial boundary value problem of fourth order wave equation with nonlinear strain and logarithmic nonlinearity at three different initial energy levels, i.e., subcritical energy E(0)d. First, we prove the local existence of weak solution by Galerkin method. In the framework of potential well, we obtain the global existence and infinite time blow up of the solution with sub-critical initial energy. Moreover by the scaling technique, we obtain global existence and infinite time blow up of the solution with critical initial energy. Also, a high energy infinite time blow up result is established.
  • Global existence and blow up of the solution for nonlinear Klein–Gordon
           equation with general power-type nonlinearities at three initial energy
    • Abstract: Publication date: Available online 6 June 2018Source: Applied Numerical MathematicsAuthor(s): Yongbing Luo, Yanbing Yang, Md Salik Ahmed, Tao Yu, Mingyou Zhang, Ligang Wang, Huichao XuAbstractThis paper investigates the local existence, global existence and finite time blow up of the solution to the Cauchy problem for a class of nonlinear Klein–Gordon equation with general power-type nonlinearities. We give some sufficient conditions on the initial data such that the solution exists globally or blows up in finite time with low initial energy and critical energy. Further a finite time blow up result of the solution with high initial energy is proved.
  • Manakov model with gain/loss terms and N-soliton interactions: Effects of
           periodic potentials
    • Abstract: Publication date: Available online 30 May 2018Source: Applied Numerical MathematicsAuthor(s): V.S. Gerdjikov, M.D. TodorovAbstractWe analyze the dynamical behavior of the N-soliton train in the adiabatic approximation of the perturbed nonlinear Schrödinger equation (NLSE) and the Manakov model. The perturbations include the simultaneous by a periodic external potential, and linear and nonlinear gain/loss terms. We derive the corresponding perturbed complex Toda chain (PCTC) models for both NLSE and Manakov model. We show that the soliton interactions dynamics for the PCTC models compares favorably to full numerical results of the original perturbed NLSE and Manakov model.
  • Stability band structure for periodic states in periodic potentials
    • Abstract: Publication date: Available online 18 May 2018Source: Applied Numerical MathematicsAuthor(s): Bin Liu, Lu Li, Boris A. MalomedAbstractA class of periodic solutions of the nonlinear Schrödinger equation with non-Hermitian potentials are considered. The system may be implemented in planar nonlinear optical waveguides carrying an appropriate distribution of local gain and loss, in a combination with a photonic-crystal structure. The complex potential is built as a solution of the inverse problem, which predicts the potential supporting required periodic solutions. The main subject of the analysis is the spectral structure of the linear (in)stability for the stationary spatially periodic states in the periodic potentials. The stability and instability bands are calculated by means of the plane-wave-expansion method, and verified in direct simulations of the perturbed evolution. The results show that the periodic solutions may be stable against perturbations in specific Floquet–Bloch bands, even if they are unstable against small random perturbations.
  • Spectral stability of travelling wave solutions in a Keller–Segel
    • Abstract: Publication date: Available online 17 May 2018Source: Applied Numerical MathematicsAuthor(s): P.N. Davis, P. van Heijster, R. MarangellAbstractWe investigate the point spectrum associated with travelling wave solutions in a Keller–Segel model for bacterial chemotaxis with small diffusivity of the chemoattractant, a logarithmic chemosensitivity function and a constant, sublinear or linear consumption rate. We show that, for constant or sublinear consumption, there is an eigenvalue at the origin of order two. This is associated with the translation invariance of the model and the existence of a continuous family of solutions with varying wave speed. These point spectrum results, in conjunction with previous results in the literature, imply that in these cases the travelling wave solutions are absolute unstable if the chemotactic coefficient is above a certain critical value, while they are transiently unstable otherwise.
  • Long wave approximation using conformal mapping for large-amplitude
           internal waves in a two-fluid system
    • Abstract: Publication date: Available online 11 April 2018Source: Applied Numerical MathematicsAuthor(s): Sunao MurashigeAbstractThis paper describes a new type of long wave model for periodic internal waves propagating in permanent form at the interface between two immiscible inviscid fluids. This model for irrotational plane motion of these waves is derived in the complex velocity potential planes where the flow domains are conformally mapped. Since no smallness assumption of wave amplitude is made and the wave elevation at the interface is represented by a single-valued function of the velocity potential, this model is applicable to large-amplitude motions of which wave profile may overhang. Numerical examples demonstrate that the proposed model can produce overhanging solutions, and variations of solutions with wavelength or wave amplitude are qualitatively similar to those of the full Euler system. It is also pointed out that the kinematic condition at the interface is exactly satisfied in the proposed model for all wave amplitudes, but not in an existing long wave model derived in the physical plane.
  • Error analysis of finite element approximations of the optimal control
           problem for stochastic Stokes equations with additive white noise
    • Abstract: Publication date: Available online 5 March 2018Source: Applied Numerical MathematicsAuthor(s): Youngmi Choi, Hyung-Chun LeeAbstractFinite element approximation solutions of the optimal control problems for stochastic Stokes equations with the forcing term perturbed by white noise are considered. To obtain the most efficient deterministic optimal control, we set up the cost functional as we proposed in [20]. Error estimates are established for the fully coupled optimality system using Green's functions and Brezzi–Rappaz–Raviart theory. Numerical examples are also presented to examine our theoretical results.
  • An adaptive nonlinear elimination preconditioned inexact Newton algorithm
           for highly local nonlinear multicomponent PDE systems
    • Abstract: Publication date: Available online 12 February 2018Source: Applied Numerical MathematicsAuthor(s): Haijian Yang, Feng-Nan HwangAbstractThis work aims to develop an adaptive nonlinear elimination preconditioned inexact Newton method as the numerical solution of large sparse multi-component partial differential equation systems with highly local nonlinearity. A nonlinear elimination algorithm used as a nonlinear preconditioner has been shown to be a practical technique for enhancing the robustness and improving the efficiency of an inexact Newton method for some challenging problems, such as the transonic full potential problems. The basic idea of our method is to remove some components causing troubles in order to decrease the impact of local nonlinearity on the global system. The two key elements of the method are the valid identification of the to-be-eliminated components and the choice of subspace correction systems, respectively. In the method, we employ the point-wise residual component of nonlinear systems as an indicator for selecting these to-be-eliminated components adaptively and build a subspace nonlinear system consisting of the components corresponding to the bad region and an auxiliary linearized subsystem to reduce the interfacial jump pollution. The numerical results demonstrate that the new approach significantly improves performance for incompressible fluid flow and heat transfer problems with highly local nonlinearity when compared to the classical inexact Newton method.
  • On a parallel, 3-dimensional, finite element solver for viscous,
           resistive, stationary magnetohydrodynamics equations: Velocity–current
    • Abstract: Publication date: Available online 3 February 2018Source: Applied Numerical MathematicsAuthor(s): K.D. Brauss, A.J. MeirAbstractWe describe a parallel implementation for the numerical approximation of solutions to the three-dimensional viscous, resistive magnetohydrodynamics (MHD) equations using a velocity–current formulation. In comparison to other formulations, the velocity–current formulation presented in this paper is an integro-differential system of equations that incorporates nonideal boundaries and nonlinearities due to induction. The solution to the equations is approximated using a Picard iteration, discretized with the finite element method, and solved iteratively with the Krylov subspace method GMRES. Effective preconditioning strategies are required to numerically solve the resulting equations with Krylov solvers [12]. For GMRES convergence, the system matrix resulting from the discretization of the velocity–current formulation is preconditioned using a simple, block-diagonal Schur-complement preconditioner based on [14]. The MHD solver is implemented using freely available, well-documented, open-source, libraries deal.II, p4est, Trilinos, and PETSc, capable of scaling to tens of thousands of processors on state-of-the-art HPC architectures.
  • Optimal error estimates of both coupled and two-grid decoupled methods for
           a mixed Stokes–Stokes model
    • Abstract: Publication date: Available online 31 January 2018Source: Applied Numerical MathematicsAuthor(s): Yuhong Zhang, Haibiao Zheng, Yanren Hou, Li ShanAbstractIn this paper, we provide a coupled algorithm and a two-grid decoupled algorithm for a mixed Stokes–Stokes model, which is coupled by a nonlinear interface transmission condition. The coupled algorithm is to discretize the mixed model directly by standard finite element method. For the two-grid decoupled algorithm, we first solve the mixed model on a coarse grid, and update the solution on a fine grid by two separated Stokes problems. Under a hypothesis about the regularity of analytical solutions, optimal error estimates for two algorithms are achieved. Several numerical tests are given to verify our theoretical results.
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
Home (Search)
Subjects A-Z
Publishers A-Z
Your IP address:
About JournalTOCs
News (blog, publications)
JournalTOCs on Twitter   JournalTOCs on Facebook

JournalTOCs © 2009-