An Heuristic Finite Difference Scheme on Irregular Plane Regions
|
|
- Maurice Parker
- 5 years ago
- Views:
Transcription
1 Applied Mathematical Sciences, Vol 8, 2014, no 14, HIKARI Ltd, wwwm-hikaricom An Heuristic Finite Difference Scheme on Irregular Plane Regions Francisco J Domínguez-Mota, Pablo M Fernández-Valdez Erika Ruiz-Díaz, Gerardo Tinoco-Guerrero and José G Tinoco-Ruiz Universidad Michoacana de San Nicolás de Hidalgo Facultad de Ciencias Físico Matemáticas Edificio B, Ciudad Universitaria Morelia, México Copyright c 2014 Francisco J Domínguez-Mota et al This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited Abstract In this paper, we present an heuristic finite difference scheme for the second-order linear operator, which is derived from an unconstrained least squares problem defined by the consistency condition on the residuals of order one, two and three in the Taylor expansion of the local truncation error It is based on a non-iterative calculation of the difference coefficients and can be used to solve efficiently Poisson-like equations on non-rectangular domains which are approximated by structured convex grids Mathematics Subject Classification: 65M06, 65M50 Keywords: finite differences, Poisson-like equations, irregular regions, convex structured grids 1 Introduction A very important problem in scientific computing is the numerical solution of Poisson-like problems It is a well studied problem when the domain has a simple geometry; unfortunately, in some specific real-world settings, the
2 672 F Domínguez-Mota et al Figure 1: Example of a typical domain when m = n =4 domain is irregular and very non-symmetric, and in this case the classical finite difference formulas do not produce accurate results As a matter of fact, there are rather few reliable finite difference schemes which can be successfully applied on irregular domains In the following sections we present a sound heuristic scheme based on a least squares problem for this case The domains of interest here are simply connected polygonal domains -mostly irregular- which can not be decomposed into rectangles For such domains, it is possible to generate suitable convex structured grids using the direct optimization method, as discussed in detail in [3, 10, 12, 13, 16] To introduce the required notation for the grids, let m and n be the number of vertical and horizontal numbers of nodes on the sides of a typical domain boundary; the latter is the positively oriented polygonal Jordan curve γ of vertices V = {v 1,,v 2m+n 2) }, and it defines the typical domain Ω IR 2 An example can be seen in figure 1 A set G = {P i,j 1 i m, 1 j n} of points of the plane with the fixed boundary positions given by V is a structured grid with quadrilateral elements for Ω, of order m n; a grid G is convex if and only if each one of the m 1)n 1) quadrilaterals or cells) c i,j of vertices {P i,j,p i+1,j,p i,j+1,p i+1,j+1 },1 i<m,1 j<n, is convex and non-degenerate Fig 2) The functional which was minimized to generate the convex structured grids of the numerical tests for this paper, as implemented in UNAMALLA [4], was
3 An heuristic finite difference scheme 673 G Figure 2: Example of a typical grid G with 4 points per side for the boundary shown in figure 1 the convex linear combination of the area functional S ω G) and length LG) with weight σ =05 See [11] for details) 2 Finite difference schemes for linear operators Let us consider the second order linear operator Lu = Au xx + Bu xy + Cu yy + Du x + Eu y + Fu 1) where A, B, C, D, E and F are smooth functions defined on an open set containing Ω; general finite difference schemes can be obtained by considering a finite set of nodes p 0 =x 0,y 0 ),, p k =x k,y k ), for which it is required to find coefficients Γ 0, Γ 1, Γ 2,, Γ k such that Γ i up i ) [Au xx + Bu xy + Cu yy + Du x + Eu y + Fu] p0 2) i=0 It is important to emphasize that there are few efficient finite difference schemes for numerical solution of partial differential equations on non-rectangular regions using finite differences [1, 2, 5, 6, 7, 9, 14] The limit value in 2) is given by the weak consistency condition [15] Γ i up i ) [Au xx + Bu xy + Cu yy + Du x + Eu y + Fu] p0 0, as p 1, p q p 0, i=0
4 674 F Domínguez-Mota et al whose Taylor expansion up to third order yields [Au xx + Bu xy + Cu yy + Du x + Eu y + Fu] p0 Γ i up i ) = i=0 ) ) F p 0 ) Γ i up 0 )+ Dp 0 ) Γ i Δx i u x p 0 ) + i=0 ) Γ i Δx i ) 2 ) Ep 0 ) Γ i Δy i u y p 0 )+ Ap 0 ) u xx p 0 ) + 2 ) Γ i Δy i ) 2 ) Bp 0 ) Γ i Δx i Δy i u xy p 0 )+ Cp 0 ) u yy p 0 ) + 2 Γ i Δx i ) 3 ) Γ i Δx i ) 2 ) Δy i u xxx p 0 )+ u xxy p 0 ) + 3! 2 Γ i Δx i Δy i ) 2 ) Γ i Δy i ) 3 ) u xyy p 0 )+ Γ i u yyy p 0 ) + 2 3! O max{δx i, Δy i }) 4 where Δx i = x i x 0, Δy i = y i y 0,1 i k In the case of the 3 3 grid defined by p 0,p 1,, p 8, the consistency condition up to third order defines the linear system Δx 1 Δx 8 0 Δy 1 Δy 8 0 Δx 1 ) 2 Δx 8 ) 2 0 Δx 1 Δy 1 Δx 8 Δy 8 0 Δy 1 ) 2 Δy 8 ) 2 0 Δx 1 ) 3 Δx 8 ) 3 0 Δx 1 ) 2 Δy 1 Δx 8 ) 2 Δy 8 0 Δy 1 ) 2 Δx 1 Δy 8 ) 2 Δx 8 0 Δy 1 ) 3 Δy 8 ) 3 Γ 0 Γ 1 Γ 2 Γ 8 = F p 0 ) Dp 0 ) Ep 0 ) 2Ap 0 ) Bp 0 ) 2Cp 0 ) which must be satisfied by the Γ i coefficients One must note that, in general, this system is not well-determined Hence, the main question is how to produce a solution of 3) which provides a consistent scheme In [8], in order to get a second order scheme, we proposed to calculate the coefficients Γ 0, Γ 1,, Γ 8 by solving the constrained optimization problem 3) min z = R R2 7 + R2 8 + R2 9, subject to R i =0,i=0,, 5 4)
5 An heuristic finite difference scheme 675 where R 0 = 8 i=0 Γ i F, R 3 = 8 Γ i Δx i ) 2 2A, R 1 = 8 Γ i Δx i ) D, R 4 = 8 Γ i Δx i )Δy i ) B, R 2 = 8 Γ i Δy i ) E, R 5 = 8 Γ i Δy i ) 2 2C, 5) R 6 = 8 Γ i Δx i ) 3 R 7 = 8 Γ i Δx i ) 2 Δy i ) R 8 = 8 Γ i Δx i )Δy i ) 2 R 9 = 8 Γ i Δy i ) 3 Problem 4) must be solved iteratively for the 3 3 structured subgrid around every inner grid node P i,j defined by G to provide the optimal local approximation to the differential operator 1) Very satisfactory results obtained with this scheme were reported in [8] In this paper, we propose an alternative heuristic scheme based on an unconstrained optimization problem which is closely related to the constrained one First, we separate the first equation of the matrix system 3) 8 Γ i F = 0 7) and then we solve the least squares problem defined by Γ Δx 1 Δx 1 8 Dp Γ 0 ) Δy 1 Δy 2 8 Ep Δx 1 ) 2 Δx 8 ) 2 0 ) = 2Ap 0 ) 8) Δx 1 Δy 1 Δx 8 Δy 8 Bp Δy 1 ) 2 Δy 8 ) 2 0 ) 2Cp Γ 0 ) 8 through the Cholesky factorization of its normal equations 6) M T M Γ=M T β, 9) where M = Δx 1 Δx 8 Δy 1 Δy 8 Δx 1 ) 2 Δx 8 ) 2 Δx 1 Δy 1 Δx 8 Δy 8 Δy 1 ) 2 Δy 8 ) 2 Γ= Γ 1 Γ 2 Γ 8, β = Dp 0 ) Ep 0 ) 2Ap 0 ) Bp 0 ) 2Cp 0 ) Next, Γ 0 is obtained from 7) One must note that, due to the lack of restrictions, this scheme is indeed less
6 676 F Domínguez-Mota et al restrictive than the previous one, and that in some very irregular grids, the residuals of order one and two might not be equal to zero; however, in spite of these facts, a large amount of numerical experiments led us to conclude that the proposed scheme is accurate enough In a similar manner to the rectangular case, the approximation to the differential operator at every inner grid point produces a linear equation whose coefficients are given by the Γ i values, and the set of all these equations can be rearranged into a sparse matrix-vector algebraic system of equations whose solution approximates the solution of 1) 3 Numerical Test For the numerical tests, we have selected 4 polygonal domains which resemble geographical regions: Great Britain ENG), Havana bay HAB), Michoacán MIC) and Ucha UCH) Using these boundaries, we considered two sets of tests: a) Convex grids with 21 points per side were generated in UNAMALLA by minimizing the functional [4] 1/2S ω G)+LG)), these grids were lated uniformly refined in order to generate grids with 41 and 81 points per side The refined grids with 41 points per side for these regions are shown in figure 3 b) Convex grids with 21, 41 and 81 points per side were generated in UNA- MALLA by minimizing the same functional The optimal grids with 41 points per side for these regions are shown in figure 4 The resulting structured grids were used to calculate the Γ i coefficients for the schemes described in the previous section In all cases, the algebraic systems obtained from 1) were assembled and later solved using a sparse Gaussian Elimination routine in MATLAB 72 in a personal computer with Intel Pentium processor B970@230 GHz and 6Gb of RAM memory To test the performance of the new scheme in the solution to the Poisson problem K u) =f, the following values for K and u were selected [5]: 1 First problem Kx, y) = ), u = 2 exp2x + y)
7 An heuristic finite difference scheme 677 Michoacán Havana Bay Ucha Great Britain Figure 3: Refined grids with 41 points per side
8 678 F Domínguez-Mota et al Michoacán Havana Bay Ucha Great Britain Figure 4: Optimal grids with 41 points per side
9 An heuristic finite difference scheme Second problem with and D = P = Kx, y) =P T DP, cos π 8 ) sin π 8 ) sin π 8 ) cos π 8 ) 1+2x 2 + y x 2 +2y 2 u = sinπx) sinπy) Function f inside the domain was chosen in such a way that u was the exact solution in both cases Moreover, for the boundary, the function u was selected as the Dirichlet condition The tests results are summarized in tables 1 and 2 for the tests set a, and in tables 3 and 4 for the test set b In all the tables, the values of the quadratic error norm is calculated as the grid function u U 2 = u i,j U i,j ) 2 A i,j, i,j where U and u are the approximated and the exact solution calculated at the i, j th grid point P i,j =x i,j,y i,j ) respectively, and A i,j is the area of the i, j th - element, calculated as the area of the polygon defined by {P i+1,j,p i,j+1,p i 1,j,P i,j 1 }; e k and ê k are the errors for the unconstrained and the constrained scheme, respectively In an analogous way, t k and ˆt k are the CPU times for the unconstrained and the constrained scheme The empirical order is calculated as log e k 1 log e k ) ), )/ log Nk log N k 1 ), where ek and N k denote the quadratic error and number of points per side, which are calculated for a given grid and the coarser one for the same region The examination of the empirical errors shows that both schemes produce numerical solutions which are very close, despite one of them, the heuristic, was defined by an unconstrained problem However, the algebraic systems assembled for the heuristic scheme are solved faster in all the tests Another important issue is that the presence of strong non convexities in the boundaries causes the generation of elongated grids cells which affect the accuracy of the numerical solution for both schemes, although the unconstrained heuristic scheme turns out to be more robust than the constrained one in most of the numerical tests and always more robust in the case of the refined grids of the test set a) Nevertheless, in a similar way to the case of finite elements on very elongated triangles, convergence is expected to be achieved [17]
10 680 F Domínguez-Mota et al Grid N k e k Order ê k Order t k ˆt k E E ENG E E E E E E HAB E E E E E E MIC E E E E E E UCH E E E E Table 1: Problem 1 data Test set a Grid N k e k Order ê k Order t k ˆt k E E ENG E E E E E E HAB E E E E E E MIC E E E E E E UCH E E E E Table 2: Problem 2 data Test set a
11 An heuristic finite difference scheme 681 Grid N k e k Order ê k Order t k ˆt k E E ENG E E E E E E HAB E E E E E E MIC E E E E E E UCH E E E E Table 3: Problem 1 data Test set b Grid N k e k Order ê k Order t k ˆt k E E ENG E E E E E E HAB E E E E E E MIC E E E E E E UCH E E E E Table 4: Problem 2 data Test set b
12 682 F Domínguez-Mota et al 4 Conclusion As follows from the numerical tests, the proposed heuristic scheme produces satisfactory results and, even though it might be in theory slightly less accurate than the scheme addressed in [8], it was as robust as the scheme defined by a local constrained optimization problem 4 In addition, in all the tests, it was notably faster due to the lack of constraints This leds to conclude that it can be used with ease in order to get reliable approximations to the solution of linear second-order partial equations in irregular regions Besides, its implementation can be made in a very simple direct way Acknowledgements We want to thank CIC-UMSNH 916, and Grant SEP-PROMEP Aplicaciones de la optimización numérica a la solución de diversos problemas de cómputo científico for the financial support for this work References [1] A Chávez-González, A Cortés-Medina and JG Tinoco-Ruiz, A direct finite-difference scheme for solving PDEs over general two-dimensional regions, Applied Numerical Mathematics, ), [2] JE Castillo and RD Grone, A matrix analysis approach to higher-order approximations for divergence and gradients satisfying a global conservation law, SIAM J Matrix Anal Appl, Vol 25, No ), [3] P Barrera-Sánchez, L Castellanos, FJ Domínguez-Mota, GF González- Flores and APérez-Domínguez, Adaptive discrete harmonic grid generation, Math Comput Simulation, ), [4] UNAMALLA, An Automatic Package for Numerical Grid Generation, [5] M Shashkov, Conservative finite difference methods on general grids, Symbolic and Numeric Computation Series, CRC Press, 1996) [6] RB Kellogg, Difference equations for the neutron diffusion equations in hexagonal geometry, Contract AT-11-1-GEN-14 [7] RB Kellogg, Difference equations on a Mesh Arising from a General Triangulation, Mathematics of Computation, Vol 18, No ),
13 An heuristic finite difference scheme 683 [8] FJ Domínguez-Mota, S Mendoza-Armenta and JG Tinoco-Ruiz, Finite Difference Schemes Satisfying an Optimality Condition, MASCOT10 Proceedings, 2011) [9] JE Castillo, JM Hyman, M Shashkov and S Steinberg, Fourth- and sixth-order conservative finite difference approximations of the divergence and gradient, Applied Numerical Mathematics, ), [10] P Knupp and S Steinberg, Fundamentals of Grid Generation, CRC Press, Boca Raton, FL, 1994) [11] P Barrera-Sánchez, FJ Domínguez-Mota, GF González-Flores and JG Tinoco-Ruiz, Generating quality structured convex grids on irregular regions, Electronic Transactions on Numerical Analysis, ), [12] SA Ivanenko, Harmonic mappings, in Handbook of Grid Generation, CRC Press, Boca Raton, FL, 1999), [13] JE Castillo, ed, Mathematical Aspects of Numerical Grid Generation, Frontiers Appl Math 8, SIAM, Philadelphia, 1991) [14] JE Castillo, Mimetic discretization methods, Routledge Chapman & Hal, 2013) [15] M Celia and W Gray, Numerical Methods for Differential Equations, Prentice-Hall, 1992) [16] JG Tinoco-Ruiz and P Barrera-Sánchez, Smooth and convex grid generation over general plane regions, Math Comput Simulation, ), [17] A Hannukainen and S Korotov, The maximum angle condition is not necessary for convergence of the finite element method, Numerische Mathematik, ), [18] S Steinberg and P J Roache, Variational grid generation, Numer Methods Partial Differential Equations, ), Received: December 7, 2013
Multigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
More informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationAntibound State for Klein-Gordon Equation
International Journal of Mathematical Analysis Vol. 8, 2014, no. 59, 2945-2949 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.411374 Antibound State for Klein-Gordon Equation Ana-Magnolia
More informationMorera s Theorem for Functions of a Hyperbolic Variable
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 32, 1595-1600 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2013.212354 Morera s Theorem for Functions of a Hyperbolic Variable Kristin
More informationImprovements in Newton-Rapshon Method for Nonlinear Equations Using Modified Adomian Decomposition Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 39, 1919-1928 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.54124 Improvements in Newton-Rapshon Method for Nonlinear
More informationOn Permutation Polynomials over Local Finite Commutative Rings
International Journal of Algebra, Vol. 12, 2018, no. 7, 285-295 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ija.2018.8935 On Permutation Polynomials over Local Finite Commutative Rings Javier
More informationRemarks on the Maximum Principle for Parabolic-Type PDEs
International Mathematical Forum, Vol. 11, 2016, no. 24, 1185-1190 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/imf.2016.69125 Remarks on the Maximum Principle for Parabolic-Type PDEs Humberto
More informationSymmetry Reduction of Chazy Equation
Applied Mathematical Sciences, Vol 8, 2014, no 70, 3449-3459 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ams201443208 Symmetry Reduction of Chazy Equation Figen AÇIL KİRAZ Department of Mathematics,
More informationA Note of the Strong Convergence of the Mann Iteration for Demicontractive Mappings
Applied Mathematical Sciences, Vol. 10, 2016, no. 6, 255-261 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.511700 A Note of the Strong Convergence of the Mann Iteration for Demicontractive
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationPerformance Comparison of Relaxation Methods with Singular and Nonsingular Preconditioners for Singular Saddle Point Problems
Applied Mathematical Sciences, Vol. 10, 2016, no. 30, 1477-1488 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.6269 Performance Comparison of Relaxation Methods with Singular and Nonsingular
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationA Family of Optimal Multipoint Root-Finding Methods Based on the Interpolating Polynomials
Applied Mathematical Sciences, Vol. 8, 2014, no. 35, 1723-1730 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.4127 A Family of Optimal Multipoint Root-Finding Methods Based on the Interpolating
More informationResearch Article Solution of the Porous Media Equation by a Compact Finite Difference Method
Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2009, Article ID 9254, 3 pages doi:0.55/2009/9254 Research Article Solution of the Porous Media Equation by a Compact Finite Difference
More informationContemporary Engineering Sciences, Vol. 11, 2018, no. 48, HIKARI Ltd,
Contemporary Engineering Sciences, Vol. 11, 2018, no. 48, 2349-2356 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.85243 Radially Symmetric Solutions of a Non-Linear Problem with Neumann
More informationFinite Difference Method of Fractional Parabolic Partial Differential Equations with Variable Coefficients
International Journal of Contemporary Mathematical Sciences Vol. 9, 014, no. 16, 767-776 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.411118 Finite Difference Method of Fractional Parabolic
More informationA Family of Nonnegative Matrices with Prescribed Spectrum and Elementary Divisors 1
International Mathematical Forum, Vol, 06, no 3, 599-63 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/0988/imf0668 A Family of Nonnegative Matrices with Prescribed Spectrum and Elementary Divisors Ricardo
More informationThe Solution of a FEM Equation in Frequency Domain Using a Parallel Computing with CUBLAS
The Solution of a FEM Equation in Frequency Domain Using a Parallel Computing with CUBLAS R. Dominguez 1, A. Medina 1, and A. Ramos-Paz 1 1 Facultad de Ingeniería Eléctrica, División de Estudios de Posgrado,
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationLinearization of Two Dimensional Complex-Linearizable Systems of Second Order Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 9, 2015, no. 58, 2889-2900 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.4121002 Linearization of Two Dimensional Complex-Linearizable Systems of
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationA Note on Product Range of 3-by-3 Normal Matrices
International Mathematical Forum, Vol. 11, 2016, no. 18, 885-891 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2016.6796 A Note on Product Range of 3-by-3 Normal Matrices Peng-Ruei Huang
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationKKM-Type Theorems for Best Proximal Points in Normed Linear Space
International Journal of Mathematical Analysis Vol. 12, 2018, no. 12, 603-609 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2018.81069 KKM-Type Theorems for Best Proximal Points in Normed
More informationThe Shifted Data Problems by Using Transform of Derivatives
Applied Mathematical Sciences, Vol. 8, 2014, no. 151, 7529-7534 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.49784 The Shifted Data Problems by Using Transform of Derivatives Hwajoon
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationFrom Binary Logic Functions to Fuzzy Logic Functions
Applied Mathematical Sciences, Vol. 7, 2013, no. 103, 5129-5138 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2013.36317 From Binary Logic Functions to Fuzzy Logic Functions Omar Salazar,
More informationThe Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive Mappings in Hilbert Spaces
Applied Mathematical Sciences, Vol. 11, 2017, no. 12, 549-560 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.718 The Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive
More informationKey-Words: - fast algorithms, boundary value problems, partial differential equations, Radon transform, MATLAB software
Fast Algorithms and MATLAB Software for Solution of the Dirichlet Boundary Value Problems for Elliptic Partial Differential Equations in Domains with Complicated Geometry ALEXANDRE GREBENNIKOV Faculty
More informationk-pell, k-pell-lucas and Modified k-pell Numbers: Some Identities and Norms of Hankel Matrices
International Journal of Mathematical Analysis Vol. 9, 05, no., 3-37 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ijma.05.4370 k-pell, k-pell-lucas and Modified k-pell Numbers: Some Identities
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationSolving Homogeneous Systems with Sub-matrices
Pure Mathematical Sciences, Vol 7, 218, no 1, 11-18 HIKARI Ltd, wwwm-hikaricom https://doiorg/112988/pms218843 Solving Homogeneous Systems with Sub-matrices Massoud Malek Mathematics, California State
More informationOn Extreme Point Quadratic Fractional. Programming Problem
Applied Mathematical Sciences, Vol. 8, 2014, no. 6, 261-277 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.312700 On Extreme Point Quadratic Fractional Programming Problem Basiya K. Abdulrahim
More informationThe Modified Adomian Decomposition Method for. Solving Nonlinear Coupled Burger s Equations
Nonlinear Analysis and Differential Equations, Vol. 3, 015, no. 3, 111-1 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/nade.015.416 The Modified Adomian Decomposition Method for Solving Nonlinear
More informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
More informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationof a Two-Operator Product 1
Applied Mathematical Sciences, Vol. 7, 2013, no. 130, 6465-6474 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2013.39501 Reverse Order Law for {1, 3}-Inverse of a Two-Operator Product 1 XUE
More informationA New Second Order Finite Difference Conservative Scheme
Divulgaciones Matemáticas Vol. 13 No. 1(005, pp. 107 1 A New Second Order Finite Difference Conservative Scheme Un Nuevo Método Conservativo de Segundo Orden en Diferencias Finitas J. M. Guevara-Jordan
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
More informationResearch Article A Novel Differential Evolution Invasive Weed Optimization Algorithm for Solving Nonlinear Equations Systems
Journal of Applied Mathematics Volume 2013, Article ID 757391, 18 pages http://dx.doi.org/10.1155/2013/757391 Research Article A Novel Differential Evolution Invasive Weed Optimization for Solving Nonlinear
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationA Mimetic Iterative Scheme for Solving Biharmonic Equations
A Mimetic Iterative Scheme for Solving Biharmonic Equations A. Gómez-Polanco a, J.M. Guevara-Jordan b, B. Molina c, a Departamento de Matemáticas, Universidad Central de Venezuela, Caracas 1020, Venezuela,
More informationQuadratic Optimization over a Polyhedral Set
International Mathematical Forum, Vol. 9, 2014, no. 13, 621-629 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.4234 Quadratic Optimization over a Polyhedral Set T. Bayartugs, Ch. Battuvshin
More informationA Numerical Solution of Classical Van der Pol-Duffing Oscillator by He s Parameter-Expansion Method
Int. J. Contemp. Math. Sciences, Vol. 8, 2013, no. 15, 709-71 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2013.355 A Numerical Solution of Classical Van der Pol-Duffing Oscillator by
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationThe Rainbow Connection of Windmill and Corona Graph
Applied Mathematical Sciences, Vol. 8, 2014, no. 128, 6367-6372 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.48632 The Rainbow Connection of Windmill and Corona Graph Yixiao Liu Department
More informationNumerical Solution of Heat Equation by Spectral Method
Applied Mathematical Sciences, Vol 8, 2014, no 8, 397-404 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ams201439502 Numerical Solution of Heat Equation by Spectral Method Narayan Thapa Department
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationA Two-step Iterative Method Free from Derivative for Solving Nonlinear Equations
Applied Mathematical Sciences, Vol. 8, 2014, no. 161, 8021-8027 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.49710 A Two-step Iterative Method Free from Derivative for Solving Nonlinear
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationRecurrence Relations between Symmetric Polynomials of n-th Order
Applied Mathematical Sciences, Vol. 8, 214, no. 15, 5195-522 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.214.47525 Recurrence Relations between Symmetric Polynomials of n-th Order Yuriy
More informationOn the Coercive Functions and Minimizers
Advanced Studies in Theoretical Physics Vol. 11, 17, no. 1, 79-715 HIKARI Ltd, www.m-hikari.com https://doi.org/1.1988/astp.17.71154 On the Coercive Functions and Minimizers Carlos Alberto Abello Muñoz
More informationSPARSE SOLVERS POISSON EQUATION. Margreet Nool. November 9, 2015 FOR THE. CWI, Multiscale Dynamics
SPARSE SOLVERS FOR THE POISSON EQUATION Margreet Nool CWI, Multiscale Dynamics November 9, 2015 OUTLINE OF THIS TALK 1 FISHPACK, LAPACK, PARDISO 2 SYSTEM OVERVIEW OF CARTESIUS 3 POISSON EQUATION 4 SOLVERS
More information+ (u u h. ( h ) 0, + (u u I ( I. + (u I. u h 2( ( I. = CB h. ; h. + ( h. + k, ( u k, ch 2(k 2) 2. = c h. + u h. + f 0,T T T h. F h. ,u h ).
The Open Numerical Methods Journal 009-5 Open Access An Adaptive Least-Squares Mixed Finite Element Method for Fourth- Order Elliptic Equations Gu Haiming * Li Hongwei and Xie Bing Department of Mathematics
More informationAsymptotic for a Riemann-Hilbert Problem Solution 1
Int Journal of Math Analysis, Vol 7, 203, no 34, 667-672 HIKARI Ltd, wwwm-hikaricom htt://dxdoiorg/02988/ijma2033358 Asymtotic for a Riemann-Hilbert Problem Solution M P Árciga-Alejandre, J Sánchez-Ortiz
More informationBasins of Attraction for Optimal Third Order Methods for Multiple Roots
Applied Mathematical Sciences, Vol., 6, no., 58-59 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ams.6.65 Basins of Attraction for Optimal Third Order Methods for Multiple Roots Young Hee Geum Department
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationOn Positive Stable Realization for Continuous Linear Singular Systems
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 8, 395-400 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4246 On Positive Stable Realization for Continuous Linear Singular Systems
More informationDynamical Behavior for Optimal Cubic-Order Multiple Solver
Applied Mathematical Sciences, Vol., 7, no., 5 - HIKARI Ltd, www.m-hikari.com https://doi.org/.988/ams.7.6946 Dynamical Behavior for Optimal Cubic-Order Multiple Solver Young Hee Geum Department of Applied
More informationMULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW
MULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW VAN EMDEN HENSON CENTER FOR APPLIED SCIENTIFIC COMPUTING LAWRENCE LIVERMORE NATIONAL LABORATORY Abstract Since their early application to elliptic
More informationDynamic Model of Space Robot Manipulator
Applied Mathematical Sciences, Vol. 9, 215, no. 94, 465-4659 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.56429 Dynamic Model of Space Robot Manipulator Polina Efimova Saint-Petersburg
More informationResearch Article Residual Iterative Method for Solving Absolute Value Equations
Abstract and Applied Analysis Volume 2012, Article ID 406232, 9 pages doi:10.1155/2012/406232 Research Article Residual Iterative Method for Solving Absolute Value Equations Muhammad Aslam Noor, 1 Javed
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationEfficient Calculation of the Weight. Distributions for Linear Codes over. Large Finite Fields
Contemporary Engineering Sciences, Vol. 9, 0, no. 3, 09-7 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ces.0.8 Efficient Calculation of the Weight Distributions for Linear Codes over Large Finite
More informationThe Expansion of the Confluent Hypergeometric Function on the Positive Real Axis
Applied Mathematical Sciences, Vol. 12, 2018, no. 1, 19-26 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.712351 The Expansion of the Confluent Hypergeometric Function on the Positive Real
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationSymmetric Identities for the Generalized Higher-order q-bernoulli Polynomials
Adv. Studies Theor. Phys., Vol. 8, 204, no. 6, 285-292 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/astp.204.428 Symmetric Identities for the Generalized Higher-order -Bernoulli Polynomials Dae
More informationProx-Diagonal Method: Caracterization of the Limit
International Journal of Mathematical Analysis Vol. 12, 2018, no. 9, 403-412 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2018.8639 Prox-Diagonal Method: Caracterization of the Limit M. Amin
More informationA Note on the Variational Formulation of PDEs and Solution by Finite Elements
Advanced Studies in Theoretical Physics Vol. 12, 2018, no. 4, 173-179 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2018.8412 A Note on the Variational Formulation of PDEs and Solution by
More informationTHE solution of the absolute value equation (AVE) of
The nonlinear HSS-like iterative method for absolute value equations Mu-Zheng Zhu Member, IAENG, and Ya-E Qi arxiv:1403.7013v4 [math.na] 2 Jan 2018 Abstract Salkuyeh proposed the Picard-HSS iteration method
More informationIntroduction into Implementation of Optimization problems with PDEs: Sheet 3
Technische Universität München Center for Mathematical Sciences, M17 Lucas Bonifacius, Korbinian Singhammer www-m17.ma.tum.de/lehrstuhl/lehresose16labcourseoptpde Summer 216 Introduction into Implementation
More informationA Fixed Point Approach to the Stability of a Quadratic-Additive Type Functional Equation in Non-Archimedean Normed Spaces
International Journal of Mathematical Analysis Vol. 9, 015, no. 30, 1477-1487 HIKARI Ltd, www.m-hikari.com http://d.doi.org/10.1988/ijma.015.53100 A Fied Point Approach to the Stability of a Quadratic-Additive
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationThe Solution of the Truncated Harmonic Oscillator Using Lie Groups
Advanced Studies in Theoretical Physics Vol. 11, 2017, no. 7, 327-335 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2017.7521 The Solution of the Truncated Harmonic Oscillator Using Lie Groups
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationCrisp Profile Symmetric Decomposition of Fuzzy Numbers
Applied Mathematical Sciences, Vol. 10, 016, no. 8, 1373-1389 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ams.016.59598 Crisp Profile Symmetric Decomposition of Fuzzy Numbers Maria Letizia Guerra
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationStudy of solar activity by measuring cosmic rays with a water Cherenkov detector
Study of solar activity by measuring cosmic rays with a water Cherenkov detector Angélica Bahena Blas 1 and Luis Villaseñor 2 1 Facultad de ciencias Físico-Matemáticas, Universidad Michoacana de San Nicolás
More informationApplication of Statistical Techniques for Comparing Lie Algebra Algorithms
Int. J. Open Problems Comput. Math., Vol. 5, No. 1, March, 2012 ISSN 2074-2827; Copyright c ICSRS Publication, 2012 www.i-csrs.org Application of Statistical Techniques for Comparing Lie Algebra Algorithms
More informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationRiesz Representation Theorem on Generalized n-inner Product Spaces
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 18, 873-882 HIKARI Ltd, www.m-hikari.com Riesz Representation Theorem on Generalized n-inner Product Spaces Pudji Astuti Faculty of Mathematics and Natural
More informationNumerical Investigation of the Time Invariant Optimal Control of Singular Systems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol. 8, 24, no. 2, 6-68 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.2988/ams.24.4863 Numerical Investigation of the Time Invariant Optimal Control of Singular Systems
More informationA Laplace Type Problems for a Lattice with Cell Composed by Three Quadrilaterals and with Maximum Probability
Applied Mathematical Sciences, Vol. 8, 1, no. 165, 879-886 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/ams.1.11915 A Laplace Type Problems for a Lattice with Cell Composed by Three Quadrilaterals
More informationDiophantine Equations. Elementary Methods
International Mathematical Forum, Vol. 12, 2017, no. 9, 429-438 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/imf.2017.7223 Diophantine Equations. Elementary Methods Rafael Jakimczuk División Matemática,
More informationIterative Methods for Single Variable Equations
International Journal of Mathematical Analysis Vol 0, 06, no 6, 79-90 HII Ltd, wwwm-hikaricom http://dxdoiorg/0988/ijma065307 Iterative Methods for Single Variable Equations Shin Min Kang Department of
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationAN EFFICIENT QUARTER-SWEEP MODIFIED SOR ITERATIVE METHOD FOR SOLVING HELMHOLTZ EQUATION
International Journal of Pure and Applied Mathematics Volume 80 No. 5 2012, 751-764 ISSN: 1311-8080 (printed version) url: http://www.ijpam.eu PA ijpam.eu AN EFFICIENT QUARTER-SWEEP MODIFIED SOR ITERATIVE
More informationOn the Three Dimensional Laplace Problem with Dirichlet Condition
Applied Mathematical Sciences, Vol. 8, 2014, no. 83, 4097-4101 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.45356 On the Three Dimensional Laplace Problem with Dirichlet Condition P.
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More information